GDAL 更新和 rgdal 错误答案

【问题标题】：GDAL Update and rgdal ErrorsGDAL 更新和 rgdal 错误
【发布时间】：2018-09-24 04:47:57
【问题描述】：

我最近更新了我的 GDAL 版本，但遇到了一些 rgdal 的问题，请参阅 here。我以为我修好了，但现在我收到了这个奇怪的警告信息。当我使用 proj4string +proj=omerc +lat_0=35 +lonc=-75 +alpha=40 +k=0.9996 +x_0=0 +y_0=0 +gamma=40 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs 将栅格加载到 R 中时，我收到以下警告：

Warning message:
In .newCRS(value) :
  +proj=omerc +lat_0=0 +lonc=0 +alpha=0 +k=0 +x_0=0 +y_0=0 +gamma=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs is not a valid PROJ.4 CRS string

由于某种原因，k 的值被更改为 0。没什么大不了的，我可以手动更改它。令人沮丧的是，当我尝试使用 writeRaster 时，我收到相同的错误消息，并且 CRS 没有写入输出 GeoTiff 文件。因此，当我将 GeoTiff 加载到 QGIS 或 R 中时，它不知道将它放在哪里 b/c 它不再是地理参考。

proj4string(input) = CRS('+proj=omerc +lat_0=35 +lonc=-75 +alpha=40 +k=0.9996 +x_0=0 +y_0=0 +gamma=40 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs')
writeRaster(input, '/home/usr/Desktop/output.tif', 'GTiff', overwrite = TRUE)
Warning message:
In .newCRS(value) :
  +proj=omerc +lat_0=0 +lonc=0 +alpha=0 +k=0 +x_0=0 +y_0=0 +gamma=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs is not a valid PROJ.4 CRS string

我还尝试在命令行上使用gdalwarp，它运行时没有警告或错误，但也不会将 CRS 添加到输出 GeoTiff 文件中。奇怪……

【问题讨论】：

其实好像proj字符串中的所有参数都设置为0，而不仅仅是k。你检查你的GDAL安装了吗？也许它缺少数据目录（包含所有 .csv 文件以及基准定义等的目录）。
@S.Monteleone 我尝试遵循post 中的建议，但没有奏效。使用here 的说明重新安装所有内容后，仍然存在问题。甚至我的 QGIS 发行版也无法找到适用于 Windows 的地理参考 tiff 文件的 CRS。
恐怕这确实是一个安装问题，因此您的问题可能与 SO 无关。您可以尝试通过 gis.stackexchange.com 获得一些反馈，但如果问题是特定于您的系统设置，则不太可能。

标签： r raster gdal rgdal

【解决方案1】：

检查工作目录您可能已将工作目录设置为您的评估者文件所在的位置。但是 rgdal 可能安装在其他地方。

【讨论】：