CPL_transform中的错误(x，crs，aoi，管道，反向)：OGRCreateCoordinateTransformation()返回NULL: PROJ可用？

Question

我使用的是ubuntu18.04，下面的代码正在生成一个错误

library(sf)
library(tmap)
library(dplyr)
library(raster)
#sudo apt install libproj-dev
#devtools::install_github("robinlovelace/geocompr")
library(spDataLarge)
if(!file.exists("e.tif"))
  download.file("https://github.com/geocompr/geocompkg/releases/download/0.1/e.tif",
                "e.tif")
elev = raster("e.tif")
urban = spData::urban_agglomerations %>% 
  filter(year == 2030) %>% 
  dplyr::select(population_millions) 
summary(urban)

tm_shape(elev) +
  tm_raster(breaks = c(-10000, 0, 10, 50, 100, 10000)) +
  tm_shape(urban) +
  tm_dots(size = "population_millions", scale = 0.5)

我得到以下错误

Error in CPL_transform(x, crs, aoi, pipeline, reverse) : 
  OGRCreateCoordinateTransformation() returned NULL: PROJ available?
In addition: Warning message:
In CPL_transform(x, crs, aoi, pipeline, reverse) :
  GDAL Error 1: No PROJ.4 translation for source SRS, coordinate transformation initialization has failed.

如果我使用

sudo apt-get install proj-bin

上面说我有最新版本

proj-bin is already the newest version (5.2.0-1~bionic0).

任何帮助都将不胜感激。

Answer 1

多亏了sf的开发人员，我现在知道了解决方案(以及问题)：

参见此处的github问题。

概述：

• 如果您用更新的GDAL版本保存save，然后在具有较旧版本的GDAL的系统上尝试st_transform，则投影信息无法正确读取(至少在Ubuntu18上发布的版本中是这样的，这也是我遇到这个问题的地方)。
• 解决方案是重新设置投影：

st_crs(data) <- 4326  # or whatever projection your data is in

如果您有多个几何列，则可能需要分别设置它：

st_crs(data$areacolumn) <- 4326