DEM数据介绍与获取方法总结

数字高程模型（Digital Elevation Model， DEM）是GIS常用的数据，对于表现一个地区的地形地貌有着很重要的作用。在这里对常见DEM数据做一个总结。

DEM介绍

介绍几种常用的DEM：

点击蓝色标题查看具体推文介绍

• 30米分辨率，全球陆地覆盖
• 下载地址：https://search.earthdata.nasa.gov/search?q=ASTER&ok=ASTER
• GEBCO海床数据,包含全球15秒分辨率的陆地和海深数据
• 下载网址：https://download.gebco.net/

GEBCO海深数据

上面的推文中附有数据获取方法，不再赘述。

此外，还有ENVI软件中附带的全球DEM数据，可以用于大范围的高程研究。

<<< 左右滑动见更多 >>>

DEM获取方法

前面的推文中给出了全球陆地和全球陆地海洋数据获取的方法，下面再介绍几个常用的：

国内的DEM下载资源

• 地理空间数据云http://www.gscloud.cn/home#page1/3
• 主要是国内的DEM
• 优势是可以根据行政区在线切割

地理空间数据云可以在线切割DEM，获取所需行政区范围内的

• 国家地球系统科学数据中心http://www.geodata.cn/data/index.html?word=DEM
• 各种DEM，数据产品，拿来即用

直接搜索，拿来即用

• 国家地球系统科学数据中心详细介绍：

如何找到靠谱的科研数据？免费的国家数据中心为你服务！

QGIS下载DEM

QGIS中可以使用SRTM-Downloader插件获取DEM

• 需要注册earthdata.nasa账号，登录
• 下载完之后仍需拼接

SRTM-Downloader插件

rgee获取DEM

这个是我目前最推荐的方法，超级灵活，范围可以定制，直接运行完代码就能出拼接好的数据，方便！

library(rgee)
library(sf)

ee_Initialize(drive = TRUE)

# Define a region of interest with sf
ee_roi <-  read_sf("./SHP/JJJ.shp")%>%
  sf_as_ee()

#选取DEM数据
srtdem <- ee$Image('CGIAR/SRTM90_V4')$   
  select('elevation')

#DEM拼接和研究区范围裁剪
mosaic = srtdem$clip(ee_roi)

imageVisParam <- list(min = 0, max = 2000, palette = c(
  '2c7bb6', 'abd9e9', 'ffffbf', 'fdae61', 'd7191c'))
Map$addLayer(mosaic, imageVisParam, 'mosaic')+
  Map$addLayer(ee_roi, visParams = list(palette = "red"), "roi")
Map$centerObject(ee_roi)

s2r <- ee_as_raster(
  image = mosaic,
  region = ee_roi$geometry(),
  scale = 90,
  dsn = './DEM/Hebei.tif',
  via = "drive"
)

京津冀DEM获取

点击阅读原文参加培训班课程学习