4. 数字地形模型

在无人机地形匹配辅助导航中,地形数据是“基准地图”,没有它就无法完成地形轮廓与实时测高的相关匹配。工程上最常用的地形表达方式是数字地形模型(Digital Terrain Model,DTM),它用规则或不规则采样点的高程集合来离散描述地表形态。理解 DTM、DEM、DSM 三者的区别,是后续做匹配区选择、基准图生成和精度评估的前提。

一、DTM、DEM、DSM 的概念边界

DTM 是广义概念,指对地表形态的数字化描述,除高程外还可包含坡度、坡向、曲率等派生属性。DEM(Digital Elevation Model)是 DTM 中最核心的高程子集,通常以规则格网存储地面高程。DSM(Digital Surface Model)则记录地表最上层表面的高程,包括植被冠层、建筑物屋顶、桥梁等人工与自然附着物。

对无人机地形匹配而言,雷达高度表或激光测距得到的是相对地表(含地物)的斜距,而气压高度、GNSS 提供的是绝对高度。若基准图使用 DEM 而实测包含树冠、屋顶,就会产生系统性偏差。因此,匹配区应优先选择地物稀少、地表裸露的 DEM 区域,或明确采用 DSM 作为基准并保持一致性。

核心结论:DEM 描述“裸地”,DSM 描述“地表最上层”,DTM 是包含高程与地形属性的总称。地形匹配导航的基准图必须与实时测高所对应的表面类型一致。

二、数据来源与采样方式

地形数据来源主要分为实测采集与遥感反演两类。实测方式包括野外测量、机载激光雷达(LiDAR)和合成孔径雷达干涉测量(InSAR);遥感方式包括立体像对摄影测量、卫星雷达和光学影像匹配。不同来源在覆盖范围、成本、精度和更新周期上差异明显。

  • 机载 LiDAR:点云密度高,可穿透植被缝隙获取地面点,适合生成高精度 DEM,但成本高、覆盖有限。
  • 摄影测量:利用无人机或卫星立体像对生成 DSM,效率高、成本低,但植被区地面高程精度受限。
  • InSAR:大范围快速获取高程,适合全球或区域尺度,但受植被和地形坡度影响较大。
  • 公开数据库:如 SRTM、ASTER GDEM、ALOS World 3D 等,覆盖广、易获取,但分辨率和精度有限。

采样方式决定数据组织形态。规则格网采样将区域划分为等间距网格,每个格网存储一个高程值,便于快速索引和插值,是导航基准图最常用的形式。不规则三角网(TIN)用三角形面片拟合地形,能更好适应地形起伏变化,但存储和检索复杂。等高线采样则以矢量线表达高程,适合制图但不便于直接匹配计算。

DTM / DEM / DSM 概念与数据来源关系 DTM 数字地形模型 高程 + 坡度 + 坡向 + 曲率 广义地形描述 导航基准图的上层概念 DEM 数字高程模型 裸地地面高程 规则格网 / TIN 地形匹配首选基准 DSM 数字表面模型 含植被、建筑、桥梁 地表最上层高程 与实测表面需一致 数据来源与采样方式 机载 LiDAR 点云密度高 可穿透植被 高精度 DEM 摄影测量 立体像对匹配 效率高成本低 多生成 DSM InSAR 干涉 大范围快速 受植被坡度影响 区域尺度 DEM 公开数据库 SRTM / ASTER ALOS World 3D 覆盖广易获取 采样方式与分辨率、精度关系 规则格网:等间距、索引快、导航基准图常用 TIN:三角形面片、适应起伏、存储检索复杂 等高线:矢量表达、适合制图、不便于直接匹配 分辨率:格网间距越小越精细 精度:RMSE、LE90、相对/绝对精度 匹配需分辨率与精度匹配

三、分辨率与精度指标

分辨率通常指格网间距,如 30 m、10 m、5 m、1 m。间距越小,地形细节越丰富,但数据量和计算量也越大。对无人机地形匹配而言,基准图分辨率应与飞行高度、测高设备波束宽度和匹配算法窗口相匹配。若基准图过粗,匹配相关峰不明显;若过细,则实时图与基准图之间的几何畸变会放大误差。

精度指标常用 RMSE(均方根误差)和 LE90(90% 线性误差)表示。RMSE 反映整体偏差水平,LE90 表示 90% 的检查点误差不超过该值。工程上还需区分绝对精度与相对精度:绝对精度决定地形匹配的绝对定位能力,相对精度影响相邻区域匹配的一致性。

指标 含义 典型要求 对匹配的影响
格网间距 相邻采样点水平距离 5~30 m 决定地形细节与相关峰锐度
高程 RMSE 高程误差均方根 ≤ 5 m(区域级) 影响匹配定位偏差
LE90 90% 点误差上限 ≤ 10 m 评估极端误差风险
绝对精度 相对大地基准的偏差 ≤ 10 m 决定绝对定位可信度
相对精度 局部区域内部一致性 ≤ 3 m 影响匹配区切换稳定性

四、常见地形数据库与工程选用

公开地形数据库中,SRTM 覆盖全球 60°N~56°S,分辨率 1 弧秒(约 30 m),适合区域级导航基准;ASTER GDEM 分辨率约 30 m,覆盖更广但噪声较大;ALOS World 3D 提供约 30 m 和 5 m 版本,精度较好;Copernicus DEM 提供 30 m 和 10 m 全球数据,更新较新。国内常用还有资源三号、高分七号等测绘卫星生成的 DEM 产品。

选用原则:先看覆盖范围,再看分辨率是否满足匹配窗口,最后核对精度指标与更新年份。植被茂密区优先选 LiDAR 生成的 DEM,城市区若实测含屋顶则需用 DSM 或做地物剔除。

五、常见误区与复盘

误区一:把 DSM 当 DEM 用。在森林或城区,DSM 与 DEM 可相差十几米,直接用于地形匹配会引入系统性高度偏差。误区二:只看分辨率不看精度。1 m 格网但高程 RMSE 达 10 m 的数据,匹配效果可能不如 10 m 格网、RMSE 2 m 的数据。误区三:忽略基准与实测的表面一致性。若实时测高包含植被,而基准为裸地 DEM,应在匹配前做地物高度补偿或改用 DSM。

地形数据检查清单:
1. 覆盖范围是否包含任务航线全部区域
2. 格网间距与飞行高度、测高波束是否匹配
3. 高程基准是否统一(EGM96 / EGM2008 / WGS84)
4. 绝对精度与相对精度是否满足指标
5. 数据年份是否与当前地物状态一致
6. 表面类型(DEM/DSM)是否与实测一致

工程实践中,建议在任务规划阶段建立地形数据质量档案,记录来源、分辨率、精度、基准和表面类型。匹配前用已知控制点做高程残差检验,若残差均值明显非零,应优先排查基准不一致问题,而不是盲目调整匹配算法参数。