Как передать данные DEM формата longlat в utm

Я скачал цифровой набор данных о высоте. Файл имеет формат .tif. Я загрузил в R. Но когда я пытаюсь преобразовать систему отсчета координат из долготы/широты в локальную 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

Как я могу решить эту проблему?


r geospatial rgdal
person Xiaoguang Ouyang    schedule 07.02.2015    source источник
comment
Пожалуйста, добавьте форматирование к вашему вопросу.   -  person Jost    schedule 07.02.2015
comment
Используйте raster::projectRaster() или (не в R) gdalwarp для преобразования растра.   -  person EDi    schedule 07.02.2015

Ответы (1)


arrow_upward
0
arrow_downward

Используйте raster::projectRaster() или (не в R) gdalwarp для преобразования растра.

person EDi    schedule 07.02.2015