Question

有没有办法“使用R编程将IMG格式图像从WGS84转换为NAD83” ？

输入图片：

名称：LST_2011-03-30_WGS.img

格式：IMG (ERDAS)

投影：UTM, Zone 12

Spheroid：WGS 84

基准：WGS 84

输出图片：

名称：LST_2011-03-30_NAD.img

格式：IMG (ERDAS)

投影：UTM, Zone 12

Spheroid：GRS 1980

基准：NAD 83

Answer 1

您可以使用假设用户在其系统上具有正常工作GDAL的包gdalUtils。如果＆＃34; gdalUtils_gdalPath＆＃34;选项已设置（通常通过gdal_setInstallation），将使用在该路径中找到的GDAL。如果找不到任何内容，将执行gdal_setInstallation以尝试查找具有正确驱动程序的工作GDAL，该驱动程序使用＆＃34;＆＃34;指定。（输出格式）参数。

更多关于inside-r.org。

这里函数gdalwarp()似乎完成了这项工作：

gdalwarp(srcfile="/your/path/LST_2011-03-30_WGS.img", #source file
         dstfile="/your/path/LST_2011-03-30_NAD.img", #destination file
         s_srs="+proj=utm +zone=12 +datum=WGS84 +no_defs +ellps=WGS84", #input spatial reference
         t_srs="+proj=utm +zone=12 +datum=NAD83 +no_defs +ellps=GRS80") #output spatial reference

N.b。：由于我没有你的图片，我无法确保它真的适用于上面提供的参数。然而，我用另一个.tif栅格和其他空间参考系统对它进行了测试，并且它有效。

Answer 2

你可以这样做

library(raster)
r <- raster('LST_2011-03-30_WGS.img')

## crs is normally defined, see
r
## but if it is not, you can set it
## crs(r) <- "+proj=utm +zone=12 +datum=WGS84 +no_defs +ellps=WGS84"

# set up an output RasterLayer
x <- raster(r)
crs(x) <- "+proj=utm +zone=12 +datum=NAD83 +no_defs +ellps=GRS80"

# compute    
x <- projectRaster(r, x)

如何使用R将WGS84转换为NAD83？

2 个答案: