全球高分辨率数字高程模型研究进展与展望

点击图片上方蓝色字体“慧天地”即可订阅

李振洪1,2  李鹏2,3  丁咚2,3  王厚杰2, 3 

1.英国纽卡斯尔大学工程学院, 英国 纽卡斯尔, NE1 7RU

2.中国海洋大学海洋地球科学学院海底科学与探测技术教育部重点实验室, 山东 青岛, 266100

3.青岛海洋科学与技术国家实验室海洋地质过程与环境功能实验室, 山东 青岛, 266061

李振洪

博士，教授，主要研究方向为影像大地测量及其在地质灾害、精准农业与人工建筑物稳定性监测方面的应用。   

摘要

简要介绍了数字高程模型（digital elevation model，DEM）的起源与定义，根据4种不同的观测平台分类介绍了DEM 数据获取方法，给出目前国际上发布的高分辨率全球DEM 的主要性质和特点。

重点介绍了9大类全球DEM，分析了DEM 质量评估相关的评定方法和精度指标。

论述了DEM 在地质灾害监测、海岸带脆弱性分析方面的应用，以美国地质勘探局和德国航空太空中心正在开展的DEM 项目为例，讨论了高精度、高分辨率全球同质DEM 和地形测深高程模型的最新需求，最后总结展望全球高分辨率DEM 的发展趋势。

引用 

李振洪, 李鹏, 丁咚, 王厚杰. 全球高分辨率数字高程模型研究进展与展望[J]. 武汉大学学报 ● 信息科学版, 2018, 43(12): 1927-1942

重要表格

全文 

地球表面的地球物理过程通常受到地形的控制或强烈影响，大尺度自然灾害与气候环境变化相关的研究（如地震、火山、滑坡、地面沉降、洪水淹没等）均需要区域或全球尺度的精密地形信息。

构建具有长期稳定大地测量参考框架（包含海岸带浅水测深数据）的高质量、高分辨全球数字高程模型（digital elevation model, DEM）是全球地球观测系统（Global Earth Observation System of Systems, GEOSS）过去十年的主要任务之一，对其正在发挥其作用的环境、气候相关的社会福利领域（societal benefit areas，SBAs）具有极其重要的价值。

本文主要介绍全球高分辨率DEM 的主要特征、数据获取与质量评估方法，分析其在地质灾害与海岸带等方面的应用，讨论其最新需求与挑战及未来发展趋势。

DEM的起源和定义

地形（topography）是指地球表面的高低起伏形态，高程是描述地表起伏形状的最基本的几何量。美国麻省理工学院摄影测量学实验室的Miller和Laflamme通过摄影测量技术采样地形数据，用于解决计算机辅助道路设计相关的工程问题，最早提出了数字地形/地面模型（digital terrain model，DTM）的概念，即利用一个任意坐标场中选择大量已知X、Y、Z的坐标点对连续地面进行简单统计表示，得到描述地球表面形态多种信息空间分布特征的有序数值阵列。基于高程数据得到的派生产品，如坡度、坡向、曲率等地貌因子，均可以作为DTM的第三维分量用于描述地形特征。

仅用于描述地面高程信息，采用有限的离散点高程实现对地形曲面的数字化模拟，以数字形式按照一定结构组织表示实际地形特征空间分布的DTM 模型称为DEM。DEM 原本是DTM 的一个子集。

由于地球表面被各种自然或人工地物覆盖，难以直接获取精确的地形形态。特别是基于摄影测量与遥感技术观测得到的DEM，通常包含了森林植被、建筑物等地物分布形态，此类DEM 又称为数字地表模型（digital surface model，DSM）。

给定平面坐标（Xi，Yi）、相应离散点高程Zi（规则或不规则），DEM 可以采用如下函数表示地面起伏形态的数字集合：

Vi =（Xi，Yi，Zi） i =1,2…n（1）

若平面位置为规则格网形式，式（2）可以简化为一维向量序列：

｛Zi ,i =1,2…n｝ （2）

20世纪90年代以来，DEM 被列为中国国家空间数据基础设施（National Spatial Data Infrastructure, NSDI）和数字地球空间数据框架（Digital Geospatial Data Framework, DGDF）的标准产品，同时也是地理国情监测的重要产品之一。作为地学分析与建模的基础数据（如海岛礁及周边复杂环境），DEM 逐步替代等高线成为地形描述与分析的重要数据源；同时，数字地球和智慧地球的建设和发展促进了DEM 与遥感学科等的一体化进程。

DEM数据获取方法

根据式（1）可知，建立DEM 实现地表重建的前提是对连续地球表面进行有限离散点的精确采样和三维坐标测量。目前，DEM 数据获取方法主要包括野外地面测量、地形图数字化、航空航天摄影测量与遥感（激光雷达扫描（light detection and ranging, LiDAR）、合成孔径雷达干涉测量（interferometric synthetic aperture radar, InSAR）、立体光学影像）等，本文根据观测平台的差异来分类介绍DEM 数据来源。

2.1地基模式

地基模式主要指野外地面测量方法。

1）传统地面测量方法，如全球导航卫星系统静态或动态测量、全站仪测量等。这类方法采用单点采集方式耗时费力，但是观测精度较高，主要适用于小范围、大比例尺测图任务，服务于工程建设项目等。

2）地面三维激光扫描仪（如加拿大Optech、奥地利Rigel、瑞士Leica等）采用无合作目标激光测距仪与角度测量系统组成自动化快速测量系统，可以快速获得海量点云数据，有效数据采集速率可达几十万点/s，可以构建高精度（mm 级）、高分辨率DTM。点云数据采集包括场地踏勘、控制网布设、靶标布设、扫描作业等步骤，点云数据处理包括配准、分析、优化压缩、分层、等高线拟合、三维建模、纹理分析等方面。另外，集成惯性测量单元（inertial measurement unit, IMU）与全球导航卫星系统等组合单元的车载移动测图系统可以支持激光扫描仪进行动态数据采集。

3）基于数字化仪通过交互式跟踪数字化或半自动扫描数字化和栅格矢量化，可以实现已有地形图要素的数字化，是目前常见的DEM 数据采集方法之一。由于原始地形图、数字化、栅格化、插值等步骤均存在误差，该方法效率和精度相对较低，但对历史记录或存档信息提取而言，可能是唯一选项。

2.2船载模式

广泛的海底地质过程导致了海底深度和粗糙度的变化，进而影响了海洋环流、生物多样性分布、海洋油气资源及地震、海啸等。由于高分辨率海底测图成本高、效率低、难度大，导致全球海底地形仍然存在许多空白区域。

船基海底地形测量是目前最常用的高分辨率海底地形测量技术。传统的船载海底地形测量主要是借助船载单波束、多波束回声测深仪开展水深测量，同步开展潮位、定位和声速测量。目前，多波束全覆盖测深技术已经成为海洋测量的常规技术，传统的点线式海底地形测量模式已发展成面状的全覆盖作业模式。

无人船、水下载人潜航器、水下自治机器人以及水下遥控机器人等平台可以集成单波束、多波束测深仪等水深测量系统，结合水下声学定位技术能够开展高效、高精度、高分辨率的水下地形测量。

2.3机载模式

机载模式主要包括航空摄影测量、机载激光LiDAR、机载激光测深系统、机载InSAR 以及无人机等方法。

1）航空摄影测量方法是利用大幅面数码航空相机拍摄的立体像对重建地表三维立体模型、量测三维空间坐标的方法，经历了模拟法、解析法、数字化3个发展阶段。目前，基于数字摄影测量工作站（如Pixel Factory Neo、VirtuoZo等）的全数字自动化或者交互式方法相对比较成熟，可用于获取高精度、高保真DEM。

2）机载LiDAR 扫描系统是集激光测距仪、全球导航卫星系统、惯性导航系统、CCD（charge-coupled device）相机、计算机数据处理等于一体的新型测量技术，可以获得地面高密度三维位置点云数据，通过多次回波滤波、分类与建模可以同时得到DSM 和DTM。机载LiDAR 技术在活断层识别、滑坡监测、海岸带脆弱性分析等领域得到了广泛应用。

3）机载激光测深系统（Laser Airborne Depth Sounder, LADS 或Airborne LiDAR Bathymetry, ALB）分别向海面发射两种波段的激光，其中红外波段被海面反射，蓝绿波段透射海水后反射，利用激光传播时间差得到相对瞬时海面的相对深度，根据潮汐资料和潮汐模型获得瞬时海面高程和海底绝对高程。该方法适合用于近岸、浅海、岛礁等近海浅水区域的水深测量，最小探测深度可达0.15m，最大测深在50m 左右，测点密度可达dm 级。目前，可以采用单一蓝绿波段激光（如美国EAARL 系统、奥地利RIEGLVQ 系统等）开展沿岸陆地与浅水的一体化无缝地形测量。

4）机载InSAR 技术适用于测绘困难地区（如多云、雾地区）的高精度、高分辨率、实时地表信息获取。美国航空航天局（National Aeronautics and Space Administration, NASA）喷气动力实验室（Jet Propulsion Laboratory, JPL）最早研发了一款机载干涉合成孔径雷达（synthetic aperture radar, SAR）测图系统（AIRSAR），又称为Topographic SAR（TOPSAR），1988-2004年搭载于DC-8喷气式飞机之上，具备P/L/C 波段全极化能力，幅宽为10km，距离向分辨率分别为7.5m/3.75 m/1.875 m。最具标志性意义的是，2000 年美国航天飞机雷达地形测绘任务（ShuttleRadarTopography Mission,SRTM）利用机载InSAR 技术，仅11d便获取了全球80%陆地1″分辨率的三维地形信息。

5）无人机（unmanned aerial vehicles, UAV）或微型无人机（mini UAV, mUAV）也可用于搭载数码相机、SAR、LiDAR 等设备，特别适用于小范围、低成本项目（如局部侵蚀、沉积、城市规划等）的DEM 获取。近年来发展起来的应用于构造地貌、林业、农业、考古等领域的运动重建技术（structure from motion, SfM）仅需要一台搭载于传统航空摄影测量平台或UAV、风筝、气球、竿子等平台的相机对感兴趣区域的影像进行重叠拍摄，基于点云后处理软件即可生成高精度、高分辨率DEM。

2.4星载模式

星载模式主要包括卫星遥感立体影像、星载InSAR、雷达测高、激光测高LiDAR 等，有利于获取全球范围质量一致的DEM。

1）高分辨率光学立体影像近年来已成为生成DEM 产品（特别是全球DEM）的首要数据源，例如法国SPOT 系列、美国Terra ASTER（Advanced Spaceborne Thermal Emissionand Reflection Radiometer）、日本ALOS（Advanced Land Observing Satellite）PRISM（PanchromaticRemoteSensingInstrumentforStereoMapping）、中国资源三号（ZY-3）等。目前，超高分辨率（very high resolution, VHR）光学遥感卫星已经越来越普遍，例如WorldView-3/4（0.31m）、WorldView-1/2（0.46 m）、GeoEye-1（0.46m）、Pleiades-1A/B（0.5 m）、KOMPSAT-3/3A（0.7/0.55 m）、Quickbird（0.65 m）、Gaofen-2（0.8 m）、TripleSat（0.8 m）、IKONOS（0.82m）、SkySat-1/2（0.8m）等。

2）星载InSAR 技术一般为重复轨道干涉测量模式，主要面向地表形变监测应用。高精度DEM 的生成容易受时间、空间去相关及大气延迟和轨道误差等影响。双星编队分布式ＳＡＲ 系统（例如德国航空太空中心发射的TerraSAR-X 与TanDEM-X 卫星）能够有效克服上述问题，已经得到了覆盖全球、空间分辨率12 m 的World DEM ，在精度方面（相对精度2 m、绝对精度4m）超过了其他全球DEM。

3）卫星LiDAR 是利用卫星平台搭载激光测高仪，向地面发射激光脉冲并接收回波信号，根据激光往返时间、卫星轨道等信息获得足印点（footprint）高程的方法，是获取月球及小行星三维地形地貌的主要方法。NASA 于2003年发射的ICESat卫星（Ice, Cloud, and land Elevation Satellite）搭载有唯一有效载荷GLAS（Geoscience Laser Altimeter System），是首个全球连续观测的星载LiDAR 传感器，可以提供全球地表高程及反演森林植被特征。GLAS每隔170m左右生成约70m直径的地面激光光斑，其地表高程观测精度为15cm。其后继卫星ICESat2于2018 年9 月发射，搭载唯一载荷ATLAS（Advanced Topographic Laser Altimeter System）以每隔70cm 进行高程观测。另外，2010年4月欧空局发射的新一代测高卫星CryoSat-2搭载的合成孔径干涉雷达高度计（Synthetic Aperture Interferometric Radar Altimeter, SIRAL）则特别适合极地冰盖与海冰的高精度高程监测。

4）基于卫星测高重力异常与船测水深数据，或者利用多光谱或SAR 卫星遥感影像构建反演模型，也可以实现大面积海底地形的反演，但在反演精度、空间分辨率等方面仍有待提高。

全球DEM

早期大范围DEM（全球、大陆或区域尺度）主要基于已有地图资料的重新编译和数据融合，空间分辨率较为粗糙（一般大于1km），例如ETOPO5（earth topography five minute grid）、TerrainBase、GLOBE（Global Land One-km Base Elevation）、GTOPO30（Global 30-arc-second Digital Elevation Models）等，但在众多领域仍然发挥了重要作用。目前常见的高分辨率全球DEM 的生成主要得益于全球卫星遥感、卫星测高、船载测深等资料的获取，例如SRTM、ASTER GDEM、AW3D30 （ALOS World 3D-30 m ）、World DEM、NEXTMAP等｡表1给出了主要全球DEM 数据的分辨率、高程基准、存储格式等特征，所有产品水平基准均为WGS84。

3.1ETOPO系列

1）ETOPO5。1988年5月，美国国家海洋与大气管理局（National Oceanic and Atmospheric Administration, NOAA）国家地球物理数据中心（National Geophysical Data Center, NGDC）发布了第一个全球陆地和海底高程数据集ETOPO5，其经纬度格网为5′（约10km），在缺少海底观测数据的亚洲、南美洲、非洲等区域，其真实空间分辨率在1°左右。

2）ETOPO2。2001年9月，NGDC 发布了经纬度格网2′（约4km）的ETOPO2，该版本融合了船载声学测量与卫星测高结果以及来自GLOBE 项目的陆地新高程数据。

3）ETOPO1。2008 年8 月，NGDC 发布了ETOPO2 v2的改进版本ETOPO1，用于支持海啸预警与建模、海洋环流建模与地球可视化，其经纬度格网为1′（约1.8km），可以获取两个不同格网版本：南极和格陵兰冰原顶部冰表面（ice surface）和冰盖底部基岩（bedrock）。

3.2GTOPO30

为了满足NASA 地球观测系统（Earth ObservingSystem, EOS）及其他全球变化研究项目对高质量全球DEM 产品的需求，1996年，美国地质调查局（United States Geological Survey, USGS）发布了全球范围空间分辨率约1km（经纬度格网为30″）的GTOPO30。由于当时没有全球覆盖的一致数据源，GTOPO30是由8种不同栅格、矢量地形数据源编译而成，其主要数据源为美国国家地理空间情报局（National Geospatial Intelligence Agency, NGA）生产的数字地面高程数据（digital terrain elevation data, DTED）和世界数字地图（digital chart of the world, DCW）。该数据发布后，成为GLOBE DEM 和ETOPO2等其他全球高程产品和大尺度应用的首要数据源。

3.3GMTED2010

在GTOPO30发布之后，全球多数地区有了新高程数据（如DTED､SRTM､ICESat等），USGS与NGA 利用11种栅格高程数据源（主要数据源为1″SRTM DTED2）融合得到了全球多分辨率地面高程数据（global multi resolution terrain elevation data2010, GMTED2010），是GTOPO30、GLOBE 及其他30″分辨率全球DEM替代版本。该产品具有30″（约1km）、15″（约500m）和7.5″（约250m）等３个不同分辨率，基于聚合方法（aggregation method）得到7类高程子产品（最小值、最大值、均值、中值、标准差、系统二次抽样、断线强调），范围主要覆盖84°N~56°S陆地部分。

3.4SRTM DEM

1）SRTMv1&v2.1。2000年2月，NASA 联合NGA､DLR､意大利航天局开展了SRTM，利用C 波段InSAR 技术收集了60°N~56°S陆地地形信息。SRTM DEM 在美国本土分辨率为1″（约30m），在全球其他地区则为3″（约90 m）。该数据是InSAR 形变监测､重力场建模等应用的关键数据源。

2）SRTM v4.1。2008年9月，国际农业研究磋商组织下属的空间信息协会（Consultative Group on International Agricultural Research, Consortium for Spatial Information, CGIARCSI）发布了3″分辨率的全球无缝高程数据集SRTM v4.1，该版本采用高级插值与补洞方法，与之前版本相比，高程质量有显著改进。

3）SRTMGL3（Global3″ V003）。2013年11月，NASA 发布了3″分辨率的SRTM v3（SRTMPlus），使用了ASTER GDEM v2､USGS GMTED2010或美国国家高程数据集（National Elevation Dataset, NED）填补数据空洞。

4）SRTMGL1（Global1″ V003）。2014 年9月，NGA 宣布1″分辨率的SRTM 数据逐步向全球用户免费开放。自2015年7月起，中国大陆用户可以免费下载。

5）X-SAR SRTM。2010年12月，DLR 发布了X 波段获取的1″分辨率SRTM DEM，覆盖全球约40%陆地表面，高程精度约为6 m。由于采用高分窄幅模式观测，X-SAR 影像条带间存在大量数据空隙。

3.5ASTER GDEM

2009年6 月，日本经济､贸易与工业部（Ministry of Economy, Trade,and Industry, METI）与NASA 合作发布了覆盖全球83°N~83°S陆地､空间格网为1″的ASTER GDEM v1。该数据利用了约150 万景ASTER 光学影像，通过立体像对相关（stereo-correlation）技术得到。2011年10 月，METI与NASA 发布了ASTER GDEM v2，新版本重新处理了2000-2010年约150万景ASTER 影像（同时增加了2008年9月以后的26 万景新影像），在受云持续覆盖的区域采用SRTM v3 等数据源替代缺失数据，通过减小归一化相关配准窗口､校正高程系统偏差（-5m）､提升水体探测能力等措施改进其空间分辨率和高程精度。

3.6AW3D系列

日本宇航局（Japan Aerospace Exploration Agency，JAXA）与日本遥感技术中心（the Remote Sensing Technology Center of Japan, RESTEC）､日本NTT DATA公司自2014年起合作开发了AW3D系列产品。

1）AW3D Standard。基于ALOS卫星搭载的PRISM 进行三维立体像对地形观测，是首个覆盖全球陆地的5m分辨率DSM/DTM，高程精度RMSE（root mean square error）为5m；

2）AW3D Enhanced。基于数字地球（Digital Globe）公司的WorldView立体像对生成3种空间分辨率（0.5 m/1 m/2 m）的DSM/DTM，无地面控制点高程精度RMSE为2m，另外提供高精度正射影像､城市､建筑物､电信､机场等栅格或矢量产品；

3）AW3D30。2016 年5 月发布了基于5 m分辨率的AW3D 标准产品得到的覆盖全球82°N~82°S 的30 m 分辨率DSM 数据集AW3D30 Version1，2018年4月发布了基于ICESat高程校正及条带误差改正后的新版本Version2.1。

3.7WorldDEM系列

WorldDEMTM 系列产品是由DLR 与法国空中客车防务与空间公司（Airbus Defense and Space, ADS）采用公私伙伴关系（Public Private Partnership, PPP）联合发布，基于TanDEM-X与TerraSAR-X雷达卫星组成双站InSAR观测模式，在2010年12月至2014年6月对全球所有陆地进行了至少两次完整覆盖。DLR负责卫星在轨运行、数据获取与干涉处理，ADS 负责World DEM数据的商业化。该系列主要包括如下产品。

1）WorldDEMcore。该产品是未经任何后处理精化的干涉输出结果DSM，包含了地表建筑物、基础设施与植被，也包含了InSAR相关的误差源、空洞等；

2）WorldDEMTM 。该产品是对地形（如异常值移除、空洞填充、噪声滤波、海岸带负值抬升等）、水文（如海洋、湖泊、河流、岛屿、海岸线、桥、堤坝、船只、海洋平台等）和机场进行编辑、精化后的DSM；

3）WorldDEM DTM。该产品是经过剔除建筑物和植被之后不含有地物特征的裸地地形。

3.8NEXTMAP系列

美国Intermap公司利用多传感器数据源（机载InSAR、星载InSAR、地面测量等）推出了全球多尺度无缝商业DSM 与DTM 三维高程数据服务NEXTMAP｡该系列主要包括如下产品｡

1）NEXTMAP ONETM ｡基于2016 至今的星载InSAR 资料，可提供全球所有国家和地区1m空间分辨率的DSM 与DTM，绝对高程精度RMSE<1m；

2）NEXTMAP5TM ｡利用Intermap 公司机载InSAR 平台，获取了全球超过40 个国家5 m空间分辨率的DSM 和DTM，绝对高程精度RMSE 为1m；

3）NEXTMAP World10TM ｡利用多源数据获取方式（机载、星载、地基等），可提供全球范围1/3弧秒（约10m 空间分辨率）的DSM 和DTM，绝对高程精度RMSE<10m；

4）NEXTMAP World30TM ｡通过校正、配准、滤波、融合多源数据集，得到覆盖全球范围1弧秒的DSM（约30m 空间分辨率），基于ASTERGDEM、SRTM 等数据源填补数据空洞和改正内部误差，利用上百万个独立控制点（如ICESat、地面控制点、机载数据等）对全部数据（包括水体与海岸线）进行三维建模校正，绝对高程精度RMSE<10m｡

3.9AIRBUS Elevation系列

ADS公司利用不同分辨率、不同数据来源的影像产品得到了多尺度DEM 产品。

1）Elevation1。基于Pleiades立体像对和三线阵光学立体影像资料生成的具有1 m格网、1.5m高程精度的DEM。

2）Elevation4。基于Pleiades立体像对和三线阵光学立体影像资料生成的具有4 m格网、2m 高程精度的DEM。

3）Elevation8。基于SPOT6立体像对和三线阵光学立体影像资料生成的具有8 m格网、3m高程精度的DEM。

4）Elevation10。基于TerraSAR-X雷达卫星数据生成的具有10 m格网、5 m高程精度的DEM。

5）Elevation30。主要基于SPOT5光学卫星数据及TerraSAR-X雷达卫星资料进行漏洞填补得到的具有30m格网、8m高程精度的DEM。

表1  主要全球DEM产品汇总

DEM质量评估

已有许多学者研究了DEM 质量相关的评估方法及高精度DEM 建模方法，系统分析了DEM 误差的影响因素（如采样、插值等）与精度建模等，本文不再赘述。本文简要介绍全球DEM质量评估通常采用的评定方法与精度指标。

4.1评定方法

DEM 精度是DEM 面上的点高程与其真值的接近程度,可以通过描述系统偏差评估其系统误差，通过估计高差的离散情况来评估其随机误差。

全球DEM 通常给出了该产品官方声称的全局精度指标，但由于DEM 数据质量在区域尺度上存在空间变化，因此需要对其定量评估并谨慎使用。本文给出DEM 质量评估的一般思路。

1）利用目视检查和剖面图，初步识别DEM的系统误差和随机误差，了解DEM 的可见细节程度。

2）利用连续运行参考站（Continuously Operating Reference Station, CORS），或者高精度（cm~dm 级）实时动态测量（real time kinematic, RTK），或者ICESat地面参考点位数据来估计ＤＥＭ 绝对高程误差。参考点的高程精度需要至少优于DEM 高程精度3 倍，通过双线性内插方法对参考点进行格网插值。

3）利用参考DEM（如LiDAR、机载InSAR或者其他已有高精度DEM 资料），以像素所在的面状区域为基础，评估待测DEM 的高程精度。

4）由于绝对高程误差与局部地貌坡度通常存在强相关，土地覆盖是DEM 精度的关键影响因子，因此需要将DEM 与参考数据按照坡度或者坡向、地物覆盖类型进行划分，用于确定高程误差空间分布与地形、地物的相关性。

4.2精度指标

假定DEM 高程误差服从正态分布，高程精度可用如下指标进行衡量。

在非开阔区域，由于滤波与插值误差导致激光扫描或数字摄影测量技术得到的DEM 误差较少服从正态分布，因此需要采用高程误差直方图和分位数Q-Q图诊断其是否服从正态分布。

考虑到系统偏差μ 和标准差σ 等精度估计会受到异常值和误差非正态分布的影响，因此DEM精度评估还需要考虑如下稳健统计方法。

NMAD 通常用于估计重尾分布标准差σ，对于正态分布（68.3% 置信区间）而言，NMAD=σ，对于90%置信区间的线性误差，LE90=1.65σ。

DEM应用

高分辨率、高精度DEM 数据在科学、工程、军事、社会经济等领域发挥着基础重要作用，特别是在地质灾害、全球变化与应对等方面，例如滑坡监测、冰川变化、海岸带淹没等。

5.1地质灾害应用

地质灾害通常定义为对生命、财产和环境具有潜在风险的地质状态和过程，包括自然灾害（如地震、滑坡、火山、海啸和洪水等）与人类活动相关的灾害（如地下水抽取、矿产开采等导致的地面沉降、水污染、大气污染等）。

滑坡灾害是一种重大的地质灾害，每年造成大量人员伤亡和经济财产损失，近年来其发生频率随着气候变化和城市化的扩张而显著增加。高分辨率、高精度DEM 是滑坡识别、滑坡监测、滑坡灾害分析与预测等阶段的关键数据集，可以用于茂密植被地区的滑坡现象与滑坡机制目视判读、滑坡目录制图、滑坡形变量求解、环境因子与触发因子分析等研究。

目前，高分辨率卫星遥感影像或LiDAR 技术被广泛用于高地质风险区域的宏观地物覆盖、裸地DEM 和地表形变获取。为了更好地融合不同观测平台的数据与模型，需要进一步结合地面实测数据（如地面三维激光扫描、地基InSAR、近景摄影测量等）和局部传感器网络开展DEM 数据的验证与补充。

5.2海岸带应用

根据联合国政府间气候变化专业委员会2014年发布的第5次气候变化评估报告，过去110a间全球海平面上升了19cm，在全球变暖大背景下，海平面上升仍呈加速趋势，预计到2100年全球海平面最大升幅为0.98m。

海平面上升会淹没滨海低地，加剧风暴潮、洪涝、海岸侵蚀、海水入侵与土壤盐渍化等灾害，威胁沿海基础设施安全。沿海城市在城市规划和重大基础设施的设计参数制定中，需要充分考虑海平面上升等气候变化增量因素。过去二十多年，欧美国家及中国已经针对沿海地区开展了海平面上升影响评估和脆弱区划研究，根据经济社会发展程度采取防护、后退和顺应等策略应对海平面上升。

高程数据是确定海岸带地区自然灾害脆弱性的主要变量，特别是与洪水和海平面上升相关的淹没灾害，但是对于DEM 高程不确定性对海平面上升评估结果的影响研究仍然较少。因此，需要了解底层高程数据质量，才能正确建模潜在影响，进而得到可靠的灾害评估结果。

DEM最新需求与发展趋势

本文以USGS 和DLR 正在开展的两类DEM 项目为例，分析全球DEM 领域的最新需求与发展趋势。

6.1高精度、高分辨率全球同质DEM

对于全球部分国家或地区而言，仅有来自多种不同数据源、不同数据处理方法得到的粗糙、不一致或者不完整的DEM。过去十年间，DLR 致力于开发全球空间分辨率12m×12m、相对高程精度2m的TanDEM-X DEM（即WorldDEM），可以为众多科学与商业应用提供不可或缺的全球同质DEM。

TanDEM-X卫星采用波长3.1cm的X波段干涉SAR，仅能穿透少部分森林树冠（大约树冠以下10%~20%左右，该比例取决于森林类型）。因此，该DEM 可代表陆地上层表面模型（DSM），包含如密集森林植被上层、城市顶层、建筑物屋顶等高程，而DTM 描述的是植被下部的地形变化。尽管TanDEM-X 采用不同基线多个相同季节观测相同区域两次以上，对于存在厚植被的农作物地区仍然无法有效描述。利用卫星遥感或激光测高（如ICESat）等方法编译得到的森林高度图，移除SRTM、ASTER GDEM或TanDEM-X DSM的树高偏移量，有望得到全球裸地DEM。

Tandem-L卫星将作为TanDEM-X 卫星的扩展任务，采用波长为24cm 的L波段SAR信号，能够穿透植被获取DSM，利用极化干涉SAR（PolInSAR）技术获取植被下部地形（DTM），基于相同波长与观测几何的双基SAR模式获取与TanDEM-X相同空间分辨率和精度的全球DEM 产品。该卫星预计2022年发射，2023年起可获得全球DSM和DTM，并且对其保持季度和年度更新。

6.2地形测深高程模型

在未来几十年，海岸线将对普遍预测的海平面上升、风暴潮和灾难性事件等导致的沿海洪水泛滥做出响应。沿海环境物理过程受控于地上地形地貌与水下测深地貌，许多地理空间数据应用需要详细的近岸地形和水深测绘信息。

USGS在2013-2023年核心科学系统策略报告中，将关键带（critical zone）的描述、理解与测图作为未来十年3大目标之一。通过获取更高精度和空间分辨率的三维陆地高程与海岸带水下地形数据，来更好地描述和理解人类活动影响与自然地表过程的关系。

自2003 年以来，USGS 非常重视利用LiDAR（机载、地面等平台）技术获取高分辨率、高精度的海岸带地形地貌信息，用于海平面上升相关的海岸带脆弱性分析与建模研究。USGS 通过开展海岸带国家高程数据库项目（Coastal National Elevation Database, CoNED），过去十年致力于开发全国无缝衔接的地形测深与高程模型（Topobathymetric Elevation Models, TBDEMs），定期发布其境内海岸带陆地、部分内陆湖泊、河流区域的潜在淹没风险制图。

6.3高分辨率DEM发展趋势

1）观测技术与观测平台趋于多样化。机载InSAR、LiDAR 等技术在工程建设领域逐渐取代传统野外测图方式，UAV、无人船、无人车、无人潜航器等观测平台的使用越来越普遍，基于多源遥感影像数据（如LiDAR、SAR、光学遥感等）融合已成为全球无缝、多尺度DEM 生成的重要方法，SfM 技术使得网络相片等免费资源利用成为可能，人工智能与机器学习方法正在加速这种趋势，将极大地降低DEMs的更新换代成本。

2）高精度、高分辨率DEM 更新换代周期逐渐缩短。得益于观测技术与观测平台的发展，全球尺度DEM 与局部DEM 的精度与空间分辨率差距逐渐缩小。基于地物变化探测应用驱动的高分辨率DEM 数据的更新换代周期逐渐缩短，从以年为周期逐渐发展到天为周期。

3）高精度、高分辨率、空间基准一致的全球无缝地形数据获取逐渐成为可能。国内外正在加快构建基于潮汐模型与海洋基准面模型的海洋垂直基准转换关系，实现沿海海域陆地与海洋高程与深度基准的转换与统一。

结语

自然界地球物理与地质过程导致了地球表面地形起伏变化，而地形变化是研究地球科学、社会经济等的关键影响因子。

本文以观测平台为例，详细阐述了DEM 数据获取方法、全球DEM 数据性质和特征、质量评估方法与精度指标，介绍了高精度DEM 数据在地质灾害监测与海岸带全球变化响应方面的应用。

通过分析USGS与DLR 正在开展的DEM 数据项目，讨论了DEM 数据的最新需求与发展趋势。观测技术与观测平台的发展加快了高精度、高分辨率、全球同质、海陆无缝衔接DEM 数据的构建，对国家海洋安全、科学研究、工程建设等具有重要的意义。

点击文末“阅读原文”即可查看原文章。

内容转载、商务活动、投稿等合作请联系

微信号：huitiandi321

邮箱：[e-mail]

欢迎关注慧天地同名新浪微博：

ID：慧天地_geomaticser

往期精彩推荐

《慧天地》敬告

《慧天地》公众号聚焦国内外时空信息科技前沿、行业发展动态、跨界融合趋势，探索企业核心竞争力，传播测绘地理信息文化，为测绘、地信、遥感等相关专业的同学提供日常学习、考研就业一站式服务，旨在打造政产学研用精准对接的平台。《慧天地》高度重视版权，对于原创、委托发布的稿件，会烦请作者、委托方亲自审核通过后才正式推发；对于来自网站、期刊、书籍、微博、微信公众号等媒介的稿件，会在作者栏或者文章开头显著标明出处，以表达对作者和推文引用平台版权的充分尊重和感谢；对于来源于网络作者不明的优质作品，转载时如出现侵权，请后台留言，我们会及时删除。感谢大家一直以来对《慧天地》的关注和支持！——《慧天地》运营团队

编辑：张自轩 审核：富裕 姚超

指导：万剑华教授