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

是否有任何方法可以使用R编程“将IMG格式图像从WGS84转换为NAD83”

输入图像：

名称：

LST\u 2011-03-30\u WGS.img

格式：

IMG（ERDAS）

投影：

UTM，12区

球体：

WGS 84

基准：

WGS 84

---

输出图像：

名称：

LST\u 2011-03-30\u NAD.img

格式：

IMG（ERDAS）

投影：

UTM，12区

球体：

GRS 1980

数据：

nad83

您可以使用包

gdalUtils

，该包假定用户的系统上有一个工作的GDAL。如果设置了“gdalUtils\u gdalPath”选项（通常由gdal\u setInstallation设置），则将使用在该路径中找到的gdal。如果未找到任何内容，将执行gdal_setInstallation以尝试查找具有使用“of”（输出格式）参数指定的正确驱动程序的工作gdal

更多关于

这里的函数

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光栅和其他空间参考系统对其进行了测试，它成功了。

您可以做到这一点

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也许会有帮助。你真的需要指定所有的参数，是吗？令人惊讶的是，没有比crs（map，name=“NAD_1983_StatePlane_Pennsylvania_South_FIPS_3702_Feet”）更优雅的解决方案了，它不需要我知道所有这些。知道
“NAD_1983_StatePlane_Pennsylvania_South_FIPS_3702_Feet”怎么容易或更容易？不管怎样，这个网站都很容易发现：它确实很长。但至少对我来说，命名的CRS对我的特定应用程序总是有很多锚——当我在密歇根州制作地图时，它是
NAD\u 1983\u StatePlane\u Michigan…
。（对我来说）这比上面提到的列表更直观。无论如何，你提供的链接是一个很好的参考；谢谢