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

Question

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

输入图像：

姓名：LST_2011-03-30_WGS.img

格式：IMG (ERDAS)

投影：UTM, Zone 12

球体：WGS 84

基准：WGS 84

---

输出图像：

姓名：LST_2011-03-30_NAD.img

格式：IMG (ERDAS)

投影：UTM, Zone 12

球体：GRS 1980

基准：NAD 83

Answer 1

您可以使用gdalUtils假设用户在他们的系统上有一个工作的 GDAL 的包。如果已设置“gdalUtils_gdalPath”选项（通常由 gdal_setInstallation 设置），则将使用在该路径中找到的 GDAL。如果没有找到，gdal_setInstallation 将被执行以尝试找到一个工作 GDAL，它具有由“of”（输出格式）参数指定的正确驱动程序。

更多关于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

Nb：由于我没有您的图像，我无法确定它是否真的适用于上述提供的参数。但是，我使用另一个 .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)