如何将longlat格式的DEM数据传输到utm
How to transfer longlat format DEM data to utm
我下载了一个数字高程数据集。该文件为“.tif”格式。我已在 R 中加载。但是当我尝试将坐标参考系统从 long/latitude 转换为本地 UTM 时，出现了问题。代码和错误如下。
dem=raster('srtm_58_04.tif')
proj4string(dem)="+proj=longlat +ellps=GRS80 +no_defs"
myproj="+proj=utm +zone=48 +north +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"
library(rgdal)
datutm=spTransform(dem,CRS(myproj))
Error in spTransform(dem, CRS(myproj)) :
load package rgdal for spTransform methods
我该如何解决这个问题？
使用raster::projectRaster() or (not in R) gdalwarp转换光栅。
我下载了一个数字高程数据集。该文件为“.tif”格式。我已在 R 中加载。但是当我尝试将坐标参考系统从 long/latitude 转换为本地 UTM 时，出现了问题。代码和错误如下。
dem=raster('srtm_58_04.tif')
proj4string(dem)="+proj=longlat +ellps=GRS80 +no_defs"
myproj="+proj=utm +zone=48 +north +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"
library(rgdal)
datutm=spTransform(dem,CRS(myproj))
Error in spTransform(dem, CRS(myproj)) :
load package rgdal for spTransform methods
我该如何解决这个问题？
使用raster::projectRaster() or (not in R) gdalwarp转换光栅。