如何將longlat格式的DEM數據轉換爲utm

我下載了一個數字高程數據集。該文件是'.tif'格式。我已經加載了R.但是當我嘗試將座標參考系統從長/緯度轉換爲本地UTM時，會出現問題。代碼和錯誤如下所示。如何將longlat格式的DEM數據轉換爲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

我該如何解決這個問題？

來源

2015-02-07 Xiaoguang Ouyang

請爲您的問題添加一些格式。 – Jost 2015-02-07 11:49:04

使用['raster :: projectRaster（）']（http://www.inside-r.org/packages/cran/raster/docs/projectRaster）或（不在R）[gdalwarp]（http：// www .gdal.org/gdalwarp.html）來轉換柵格。 – EDi 2015-02-07 20:55:12

使用raster::projectRaster()或（不在R）gdalwarp來轉換光柵。

來源

2015-02-07 20:58:10 EDi

如何將longlat格式的DEM數據轉換爲utm

回答

相關問題