数字高程模型和表面模型（全球高分辨率数字高程模型研究进展与展望）

数字高程模型和表面模型（全球高分辨率数字高程模型研究进展与展望）20世纪90年代以来，DEM被列为中国国家空间数据基础设施(National Spatial Data Infrastructure NSDI)和数字地球空间数据框架(Digi- tal Geospatial Data Framework DGDF)的标准产品，同时也是地理国情监测的重要产品之一[4-7]。作为地学分析与建模的基础数据(如海岛礁及周边复杂环境[8]，DEM逐步替代等高线成为地形描述与分析的重要数据源;同时，数字地球和智慧地球的建设和发展促进了DEM与遥感学科等的一体化进程[9-10]。若平面位置为规则格网形式，式(1)可以简化为一维向量序列：摘 要：简要介绍了数字高程模型(digital elevation model，DEM)的起源与定义，根据4种不同的观测平台分类介绍了DEM数据获取方法，给出目前国际上发布的高分辨率全球DEM的主要性质和特点。重点介绍了9大类全球D

来源：《武汉大学学报·信息科学版》2018年12期

作者：李振洪， 李鹏，丁咚， 王厚杰

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

提示：其他平台转载请务必注明原始出处及作者

摘 要：简要介绍了数字高程模型(digital elevation model，DEM)的起源与定义，根据4种不同的观测平台分类介绍了DEM数据获取方法，给出目前国际上发布的高分辨率全球DEM的主要性质和特点。重点介绍了9大类全球DEM，分析了DEM质量评估相关的评定方法和精度指标。论述了DEM在地质灾害监测、海岸带脆弱性分析方面的应用，以美国地质勘探局和德国航空太空中心正在开展的DEM项目为例，讨论了高精度、高分辨率全球同质DEM和地形测深高程模型的最新需求，最后总结展望全球高分辨率DEM的发展趋势。

关键词：数字高程模型 合成雷达干涉测量 光学立体像对 地质灾害 质量评估

地球表面的地球物理过程通常受到地形的控制或强烈影响，大尺度自然灾害与气候环境变化相关的研究(如地震、火山、滑坡、地面沉降、洪水淹没等)均需要区域或全球尺度的精密地形信息。构建具有长期稳定大地测量参考框架(包含海岸带浅水测深数据)的高质量、高分辨全球数字高程模型(digital elevation model DEM)是全球地球观测系统(Global Earth Observation System of Systems GEOSS)过去十年的主要任务之一，对其正在发挥其作用的环境、气候相关的社会福利领域(societal benefit areas SBAs)具有极其重要的价值。

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

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

2 DEM数据获取方法

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

2.1 地基模式

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

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

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

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

2.2 船载模式

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

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

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

2.3 机载模式

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

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

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

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

4) 机载InSAR技术适用于测绘困难地区(如多云、雾地区)的高精度、高分辨率、实时地表信息获取[33-35]。美国航空航天局(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波段全极化能力，幅宽为10 km，距离向分辨率分别为7.5 m/3.75 m/1.875 m[36-37]。最具标志性意义的是，2000年美国航天飞机雷达地形测绘任务(Shuttle Radar Topography Mission SRTM)利用机载InSAR技术 仅11 d便获取了全球80%陆地1″分辨率的三维地形信息。

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

2.4 星载模式

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

1) 高分辨率光学立体影像近年来已成为生成DEM产品(特别是全球DEM)的首要数据源[42-43]，例如法国SPOT系列、美国Terra ASTER(Advanced Spaceborne Thermal Emission and Reflection Radiometer)、日本ALOS(Advanced Land Observing Satellite) PRISM(Panchromatic Remote Sensing Instrument for Stereo Mapping)、中国资源三号(ZY-3)等[44-46]。目前，超高分辨率(very high resolution VHR)光学遥感卫星已经越来越普遍，例如WorldView-3/4(0.31 m)、WorldView-1/2(0.46 m)、GeoEye-1(0.46 m)、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.82 m)、SkySat-1/2(0.8 m)等。

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

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

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

3 全球DEM

早期大范围DEM(全球、大陆或区域尺度)主要基于已有地图资料的重新编译和数据融合，空间分辨率较为粗糙(一般大于1 km)，例如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)、WorldDEMTM、NEXTMAP等。表 1给出了主要全球DEM数据的分辨率、高程基准、存储格式等特征 所有产品水平基准均为WGS84。

表 1 主要全球DEM产品汇总

3.1 ETOPO系列

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

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

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

3.2 GTOPO30

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

3.3 GMTED2010

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

3.4 SRTM DEM

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

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

3) SRTMGL3 (Global 3″ V003)。2013年11月，NASA发布了3″分辨率的SRTM v3(SRTM Plus)，使用了ASTER GDEM v2、USGS GMTED2010或美国国家高程数据集(National Elevation Dataset NED)填补数据空洞[61]。

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

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

3.5 ASTER GDEM

2009年6月，日本经济、贸易与工业部(Ministry of Economy Trade and Industry METI)与NASA合作发布了覆盖全球83°N~83°S陆地、空间格网为1″的ASTER GDEM v1。该数据利用了约150万景ASTER光学影像，通过立体像对相关(stereo-correlation)技术得到[63-64]。

2011年10月，METI与NASA发布了ASTER GDEM v2，新版本重新处理了2000-2010年约150万景ASTER影像(同时增加了2008年9月以后的26万景新影像)，在受云持续覆盖的区域采用SRTM v3等数据源替代缺失数据，通过减小归一化相关配准窗口、校正高程系统偏差(-5 m)、提升水体探测能力等措施改进其空间分辨率和高程精度[65]。

3.6 AW3D系列

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

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

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

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

3.7 WorldDEMTM系列

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

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

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

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

3.8 NEXTMAP系列

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

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

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

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

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

3.9 AIRBUS Elevation系列

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

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

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

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

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

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

4 DEM质量评估

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

4.1 评定方法

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

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

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

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

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

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

4.2 精度指标

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

1) 误差均值

2) 标准差

3) RMSE

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

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

1) 中位数

2) 绝对偏差中位数(median absolute deviation MAD)

3) 归一化绝对偏差中位数(normalized median NMAD)

4) 90%分位数绝对偏差(LE90)

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

5 DEM应用

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

5.1 地质灾害应用

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

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

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

5.2 海岸带应用

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

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

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

6 DEM最新需求与发展趋势

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

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

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

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

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

6.2 地形测深高程模型

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

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

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

6.3 高分辨率DEM发展趋势

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

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

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

7 结语

自然界地球物理与地质过程导致了地球表面地形起伏变化，而地形变化是研究地球科学、社会经济等的关键影响因子。本文以观测平台为例，详细阐述了DEM数据获取方法、全球DEM数据性质和特征、质量评估方法与精度指标，介绍了高精度DEM数据在地质灾害监测与海岸带全球变化响应方面的应用。通过分析USGS与DLR正在开展的DEM数据项目，讨论了DEM数据的最新需求与发展趋势。观测技术与观测平台的发展加快了高精度、高分辨率、全球同质、海陆无缝衔接DEM数据的构建，对国家海洋安全、科学研究、工程建设等具有重要的意义。

致谢: 感谢DLR TerraSAR-X(No. IDEM_CALVAL0059)与TanDEM-X(No. DEM_CALVAL1423)项目对本文前期研究的支持。

参考文献：略

更多阅读：CGCS2000 国土三调 HEC-RAS 注册测绘师