我最近更新了我的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文件中。奇怪...
答案 0 :(得分:0)
检查工作目录 您可能已将工作目录设置为评估者文件所在的位置。但是rgdal可能安装在其他地方。