查看原文
其他

【干货】DEM知识及获取方式(赠全球DEM数据、技术规范))

达北云 多普云 2022-12-04

本文主要介绍了DEM成图方法及生产工艺流程、相关术语与定义、精度要求及讲解、DEM如何获取、DEM如何使用。希望可以帮助到您。

1.1 成图方法及生产工艺流程

    大于1:5万(包括1:5万)比例尺的数字高程模型应从地形图、航空象片或实测数据生成,小于1:5万比例尺的数字高程模型应从大于该比例尺的数字高程模型综合。

1.2术语和定义

  • 属性数据  attribute data

描述确定地理实体的类别、级别等质量特征和数量特征的数字信息。 

  • 拓扑数据  topological data

描述地图上点、线、面状要素之间关联、邻接、包含等空间关系的数据。

  • 空间数据  spatial data

用来表示地理实体的位置、形状、大小和分布特征的信息,以及实体的空间关系(如拓扑关系)信息的数据。

  • 元数据  metadata

关于数据的数据,是数据和信息资源的描述性信息。如有关数据源、数据分层、产品归属、空间参考系、数据质量、数据更新以及图幅接边等的数据。

  • 数字高程模型  digital elevation model (DEM)

定义在X、Y域离散点(矩形或三角形)上以绝对高程表示地面起伏形态的数据集合。

  • 原始地形数据  original terrain data

利用摄影测量方法直接获取,符合统一分类编码要求,并经拼接截幅处理的图幅地物、地貌空间数据。

  • 数字地形图入库数据  the input-oriented data for database of digital topographic maps

在原始地形数据的基础上,经规格化处理、拓扑关系处理、属性赋值、地名指针挂接等,且满足矢量地形图数据库结构要求的矢量数据。

  • 数字地形图出版数据  the print-oriented data for database of digital topographic maps

对数字地形图入库数据进行符号化、分版、图廓整饰等处理并符合军用数字地形图出版要求的矢量数据。

  • 地图投影、大地坐标系统、高程系统和深度基准

1:10 000~1:100 000比例尺军用数字地形图采用克拉索夫斯基椭球,高斯—克吕格投影,1:10 000军用数字地形图按经差3º分带,1:25 000~1:100 000其他比例尺军用数字地形图按经差6º分带。大地坐标系统采用1954年北京坐标系。高程系统采用1985国家高程基准。以理论最低潮面作为海洋深度基准面。

  • 图幅分幅与编号

1:10 000~1:100 000比例尺军用数字地形图,按国际分幅法分幅,特殊情况也可合幅成图。图幅分幅与编号执行GJB 4447-2002。

  • 地形类别与基本等高距

    A.地形类别

地形类别根据图幅内的地面坡度和高差分为平地、丘陵地、山地、高山地。

a)平地:地势平坦,图幅内大部分坡度小于2º,高差在40m以下;

b)丘陵地:地势起伏较缓,图幅内大部分坡度在2º~6º,高差在40 m~300m;

c)山地:图幅内大部分坡度在6º~25º,高差在300 m~1000m;

d)高山地:图幅内大部分坡度在25º以上,高差在1000m以上;

e)当高差与坡度发生矛盾时,以坡度划分为准。

     B.基本等高距

    a)基本比例尺地形图的基本等高距根据地形类别确定,具体内容见表1。

    表1基本等高距与成图比例尺的                  单位为m

成图比例尺

1:10   000

1:25   000

1:50   000

1:100   000

平  地

2.5

5

10

20

丘陵地

2.5

5

10

20

山  地

5

10

20

40

高山地

10

20

40

80

b)一幅图内只采用一种基本等高距,基本等高距不能充分显示地貌特征时,可加测间曲线和助曲线。

  • 高程注记点密度

高程注记点应均匀分布,其密度在图上每100cm2内平地、丘陵地为10个~20个;山地、高山地为8个~15个。

1.3 精度要求

  • 中误差

    数字高程模型对邻近高程控制点的高程中误差不得大于下表的规定

  • 最大误差

    数字高程模型格网点高程误差大于中误差的点数不得超过总点数的10%,其最大高程误差不得大于中误差的2倍

  • 接边精度

    各数字高程模型接边时,不应出现漏洞、重叠,起始结尾数据应当吻合。

1.4 DEM数据来源库

地形DEM数据
国内:中科院-地理空间数据云
美国:USGS - National Map
全球:USGS - Earth Explorer
USGS - Global DataExplorer
SRTM - 90m DEM Digital Elevation Database
海深DEM数据
全球:GEBCO - Gridded Bathymetric Data

1.5  获取方式

  1. 使用地理空间数据云平台

1)地址:http://www.gscloud.cn/search#dlv=Wzg4LFswLDEwLDEsMF0sW1siZGF0YWRhdGUiLDBdXSxbXSw5OV0%3D

2)获取方式:

第一步 : 如何选择数据类型?

打开网页,点击右上角完成注册之后,点击左上角“数据集”,选择DEM数字高程数据,30M分辨率数字高程数据(精度比90M更高),确定。当然该网站还提供其他一系列卫星遥感数据,比如LANDSAT8,之后我们会推出文章解释什么是LANDSAT8,如何利用LANDSAT8来做地表覆被类型分析。

第二步 : 选择下载范围

你可以选择行政区确定下载边界,比如北京市,海淀区,网站会自动加载一个mask,来依据该边界下载mask内部的数据。

 

或者,你也可以划定多边形范围来下载数据。选择“地图选择”,可以用矩形自定义框选北京市区范围。然后点击绿色按钮“搜索”。


第三步 : 下载数据

搜索之后,可以看到我们选择的数据范围横跨了两个卫星数据采集单元。依次点击下载按钮获取这两个数据。下载后可以看到格式为tiff,后文中我们会讲解怎么在GIS中生成等高线。

  B. 使用Global Mapper 

  C.使用LocaSpaceViewer

      上百种地图随意下载,包括谷歌地球影像、谷歌地图、谷歌地形、谷歌历史影像、百度、高德、腾讯、天地图、天地图分省、OSM、地质图、海图等等,无水印、无偏移、无限制下载。

D.使用http://earthexplorer.usgs.gov/
使用EarthExplore下载DEM数据,类似的,我们需要选择我们感兴趣的数据下载范围。进入页面之后在地图上直接点击就会出现红色的坐标点。点击三个以上的坐标点即可生成mask。这边我们以马萨诸塞州的波士顿为例。



或者,你也可以以行政边界选择所需区域。在Predefined Area标签中,选择马萨诸塞州作为mask边界。



在Data Set标签中,选择ASTER GLOBAL DEM,该数据集精度为30m。点击result就可以下载啦。


 E.使用https://search.asf.alaska.edu

还有一些其他数据下载来源,就不一一赘述,在这里放上链接。

全球地形DEM数据
https://scihub.copernicus.eu/
https://www.asf.alaska.edu/sar-data/palsar/

USGS - Global Data Explore

https://gdex.cr.usgs.gov/gdex/

SRTM - 90m DEM Digital Elevation Database 

http://srtm.csi.cgiar.org/srtmdata/


全球海深DEM数据

GEBCO - Gridded Bathymetric Data

https://www.gebco.net/data_and_products/gridded_bathymetry_data/

1.6  如何使用DEM

由于DEM是带有地理坐标信息的栅格数据,我们可以在ArcGIS中将其进行矢量计算, 得出等高线。但是由于各个国家、地区在不同历史时期数据采集方式的不同,我们在利用dem数据生成等高线的时候需要选择适合该场地的地理坐标系统。

打开ArcGIS,会弹出一个对话框,我们在my templates标签中选择blank map。

在左边的layer右击, 选择add data,找到你下载下来的tif格式的dem数据,点击add。

然后我们可以看到Boston harbour的文件。

那么现在就离生成等高线只有一步之遥了,但是在这之前我们还要设置地理坐标系统。现在让我们理解一下什么是地理坐标系统(coordinate system)。

右键左上角layer,在Coordinate System标签中,我们会看到有“Geographic” 和“Projected”两种坐标系统。他们有什么区别呢?



从这张图示可以看出, Geographic Coordinate System是三维的坐标系统,但是如果我们要把三维的坐标数据变成二维平面,就需要用投影的方法。这就是Projected Coordinate System。

 

而投影又有许多不同的方式,他们的变形程度是不一样的。比如你可以把地球塞进一个横过来的圆柱体里,然后在地球中心放一盏灯,想象一下将地球上的所有的坐标投到圆柱体表面。然后将圆柱体展开,你就可以得到我们最常用的UTM 投影坐标系统(Universal Transverse Mercatorcoordinate system)



还有很多很神奇的投影方式,比如把地球塞到一个圆锥里面然后展开。。。



再比如把地球像剥柚子一样剥开。。。




那么在实际操作中,如果你画的是平面,需要后续在rhino中生成地形,我们会强调一定要使用Projected coordinate system,因为rhino的是不能很好的识别Geographic Coordinate System产生的数据的。那怎么选择Projected coordinate system呢?通常我们会根据场地所处的UTM带号来选择。同学们可以在这张map中找到你场地的UTM带号,以Boston,MA为例,我们可以找到是在北半球19号投影带中。



在coordinate system的选择路径中,选择projected coordinate system - UTM - NAD (1983) -NAD 1983 UTMZone 19N这里的NAD是North American Datum的缩写,是针对北美各个带号为投影中心的投影方法,使得北美洲的投影变形最小。



应用之后,我们的图像是不是发生了变化?这是因为图像四角根据NAD 1983 UTM Zone 19N重新投影了。



接下来就可以生成等高线了。点击菜单栏上的Arctoolbox,在弹出的侧边栏中选择spatial analysttools-surface-contour。



在弹出的对话框中选择tif,将contour interval 设置为你需要的等高距,然后重点!在environment中选择same as display,意思是你导出的等高线投影方式和显示是一样的(否则默认same as layer的话,刚才修改投影方式就白做了)点击ok就大功告成了!



这里再介绍一个小技巧。如果你只需要局部的等高线怎么办呢?step1,将图放大到你场地局部。Step2,在environment中,将processing extent选择为same as display。其他同上。等高距可以根据需要设置的小一点,这里设置为5m,点击ok。



为了导出带高程的CAD,在生成的contour 图层上右键,点击open attribute table,有一列Contour就是高程数值。但是由于CAD只认Elevation字段,所以我们在table下add field,添加一列elevation数据。




在生成的elevation标题上右击选择Field calculator,选择contour,来赋值给elevation。



可以看到赋值完成了。那么我们导出cad。



右键你的contour图层,选择export - export to cad,在output file中修改保存路径。



导入rhino,成功!



有了等高线我们可以生成地形,研究雨洪径流模拟!



可以结合街道网格研究积水情况!


诸多应用情况可以按自己的需要使用。

如若有其它的下载方式或需求,可以留言告诉小编哦!

关注公众号,回复【资料9】,即可获取DEM数据SRTMv3_1 30米精度全国范围和DEM行业规范指南,更多干货,敬请期待

点击阅读原文,了解更多精彩内容

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存