如何使用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，该GDAL具有使用"of"(输出格式)参数指定的正确驱动程序。

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

的被害者。:因为我没有你的图像，我不能确保它真的与上面提供的参数工作。然而，我用另一个。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)