4. 数字地形模型（DTM/DEM）

数字地形模型（DTM）与数字高程模型（DEM）是地形匹配导航系统的“基准底图”。地形匹配导航的核心思想是将飞行器实时测得的地形剖面与机载存储的参考地形进行相关匹配，从而确定飞行器位置。因此，参考地形数据的质量直接决定了匹配算法的收敛概率、定位精度和鲁棒性。本章围绕地形数据的来源、分辨率与精度、坐标系与格式标准、数据预处理方法四个维度展开，帮助学员建立从原始数据获取到导航可用地形数据库的完整工程链路。

4.1 地形数据来源：从卫星到地面的多平台获取

地形匹配导航所需的地形数据来源可分为三大类：卫星遥感、航空测绘和地面/机载激光雷达。不同来源的数据在覆盖范围、分辨率、精度和成本上差异显著，工程实践中需要根据任务剖面和匹配算法需求进行选择。

4.1.1 卫星遥感数据

卫星遥感是目前获取全球尺度地形数据的主要手段。最具代表性的是SRTM（航天飞机雷达地形测绘任务）和ASTER GDEM。SRTM采用C波段和X波段干涉合成孔径雷达（InSAR）技术，在2000年完成了全球60°N至56°S范围的地形测绘，公开版本分辨率分别为30m（SRTM1）和90m（SRTM3）。ASTER GDEM由日本METI和美国NASA联合发布，基于光学立体像对提取，分辨率30m，但垂直精度略低于SRTM。

对于地形匹配导航而言，卫星数据的优势在于全球覆盖、成本低、获取便捷；劣势在于垂直精度有限（LE90通常为6~16m），且在高山、峡谷等剧烈地形区域容易出现数据空洞或异常值。因此，卫星数据更适合作为战略级导航数据库的基础底图，或用于任务规划阶段的匹配区域初选。

4.1.2 航空测绘与机载LiDAR

当任务区域范围有限（如数百平方公里）且对定位精度要求较高时，航空测绘和机载LiDAR成为首选。航空摄影测量通过立体像对重建地形，分辨率可达0.5~2m，垂直精度LE90约1~3m。机载LiDAR（激光雷达）通过发射激光脉冲并测量回波时间，直接获取地面点云数据，在植被稀疏区域垂直精度可达0.1~0.3m，是当前高精度地形数据获取的标杆技术。

LiDAR数据的核心优势在于能够穿透植被冠层获取真实地表高程（DTM），而非植被表面高程（DSM）。这一特性对于地形匹配导航至关重要——匹配算法需要的是物理地表形态，而非植被覆盖层。但LiDAR数据获取成本高、处理复杂，且需要地面控制点进行精度校验，通常仅用于重点任务区域的高精度数据库建设。

4.1.3 数据源选择策略

4.2 分辨率与精度：匹配性能的决定性因素

地形匹配导航对DEM数据的分辨率和精度有明确要求。分辨率决定了地形特征的空间表达粒度，精度则决定了地形高程值的可信程度。两者共同影响匹配算法的捕获概率（C_p）和定位精度（CEP）。

4.2.1 分辨率对匹配性能的影响

DEM分辨率（网格间距）直接影响地形剖面的采样密度。设飞行器雷达高度表测量间隔为Δx，则DEM网格间距应满足Δg ≤ Δx/2（奈奎斯特采样定理），否则地形剖面中的高频特征将被丢失，导致匹配模糊。例如，若高度表每50m输出一个测点，则DEM网格间距不应大于25m。对于TERCOM算法，通常要求匹配窗口内至少包含20×20个有效网格点，因此分辨率过低会直接缩小有效匹配窗口，降低匹配概率。

但分辨率并非越高越好。过高的分辨率（如1m）会带来三个问题：数据存储量急剧增大（1°×1°区域在1m分辨率下约需3.6GB）、匹配计算量呈平方级增长、以及高频噪声被放大。工程实践中，应根据飞行器高度、速度、高度表测量间隔和机载计算机算力综合确定最优分辨率。

4.2.2 精度指标与误差预算

DEM精度通常用LE90（90%线性误差）或RMSE（均方根误差）表示。垂直精度对匹配性能的影响体现在两个方面：一是高程误差会降低地形剖面的信噪比，导致相关峰模糊；二是系统性的高程偏差（如基准面不一致）会引入匹配位置偏移。水平精度则直接影响匹配结果的定位精度。

在误差预算设计中，需要将DEM误差、高度表测量误差、惯导漂移误差三者统筹考虑。经验公式为：系统总误差σ_total = √(σ_DEM² + σ_alt² + σ_INS²)。当σ_DEM > 2σ_alt时，地形数据误差成为主导因素，此时提升高度表精度对匹配性能改善有限，应优先改进DEM质量。

4.3 坐标系与格式标准：数据互操作的基础

地形匹配导航系统需要将DEM数据、飞行器位置、高度表测量值统一到同一坐标框架下。坐标系不统一是数据预处理中最常见且最危险的错误源。

4.3.1 坐标基准体系

水平坐标基准方面，WGS-84是全球导航卫星系统（GNSS）和绝大多数DEM数据的默认基准。但部分历史数据可能基于北京54、西安80或CGCS2000等局部基准，使用时必须进行基准转换。垂直基准方面，DEM高程通常基于EGM96大地水准面或WGS-84椭球面。EGM96与椭球面之间的差距在全球范围内可达-107m至+85m，若混用将产生系统性高程偏差。

投影方式方面，地形匹配导航通常采用UTM投影或等距圆柱投影。UTM投影在6°带内变形小，适合中纬度区域；等距圆柱投影计算简单，适合全球统一存储。但需注意，UTM投影在相邻带边缘存在变形不连续问题，跨带飞行时需进行投影转换。

4.3.2 数据格式标准

军事领域最常用的地形数据格式是DTED（Digital Terrain Elevation Data），由NGA制定。DTED Level 0、1、2分别对应约900m、90m、30m分辨率，每个文件覆盖1°×1°区域，采用二进制存储，包含高程数据块和元数据头。民用领域则广泛使用GeoTIFF、ASCII Grid、NetCDF等格式。GeoTIFF在TIFF图像基础上嵌入地理参考信息，支持浮点高程值，是目前最灵活的格式。

在系统设计中，建议将原始数据统一转换为内部标准格式（如NetCDF或自定义二进制格式），并建立元数据管理机制，记录数据源、采集时间、分辨率、精度、坐标基准、投影参数等关键信息。元数据缺失的数据应视为不可用数据。

// DTED Level 2 文件头解析示例（C语言伪代码）
typedef struct {
    char  sentinel[8];      // "DTED02  "
    char  origin_lat[8];    // 原点纬度（毫秒）
    char  origin_lon[8];    // 原点经度（毫秒）
    char  lat_interval[4];  // 纬度间隔（毫秒）
    char  lon_interval[4];  // 经度间隔（毫秒）
    char  num_lat_lines[4]; // 纬度方向行数
    char  num_lon_lines[4]; // 经度方向列数
    // ... 其他字段
} DTED_Header;

// 高程数据读取（单位：米，有符号16位整数）
int16_t elevation = (int16_t)((buffer[0] & 0xFF) << 8 | (buffer[1] & 0xFF));

4.4 数据预处理方法：从原始数据到导航数据库

原始DEM数据通常包含粗差、空洞、噪声和坐标系不一致等问题，必须经过系统化预处理才能用于地形匹配导航。预处理流程包括：坐标统一、粗差剔除、空洞填补、滤波平滑、精度验证五个核心步骤。

4.4.1 坐标统一与重采样

首先将所有数据统一到WGS-84坐标系和EGM96垂直基准。若原始数据使用其他基准，需通过七参数转换或大地水准面模型进行修正。随后根据导航需求确定目标分辨率和投影方式，进行重采样。重采样方法包括最邻近插值（保留原始值，适合分类数据）、双线性插值（平滑过渡，适合连续地形）和三次卷积插值（保留边缘特征，计算量大）。对于地形匹配导航，推荐使用双线性插值，在平滑噪声和保留地形特征之间取得平衡。

4.4.2 粗差剔除与空洞填补

粗差（异常值）通常表现为孤立的高程突变点，可通过中值滤波或坡度阈值检测识别。具体做法是：计算每个网格点与周围8邻域的中值之差，若差值超过设定阈值（如3倍局部标准差），则判定为粗差并剔除。空洞（数据缺失区域）则需通过插值填补。常用的插值方法包括：反距离加权（IDW）、克里金插值（Kriging）和样条插值。对于地形匹配导航，克里金插值能更好地保持地形自相关性，但计算复杂度高；IDW简单快速，适合实时处理。

4.4.3 滤波平滑与特征增强

滤波的目的是去除高频噪声，同时保留对匹配有用的地形特征。常用的滤波方法包括高斯滤波、中值滤波和各向异性扩散滤波。高斯滤波平滑效果好但会模糊边缘；中值滤波能有效去除脉冲噪声但可能损失细节；各向异性扩散滤波在平滑噪声的同时保留边缘特征，是较优选择。滤波窗口大小需根据DEM分辨率和地形特征尺度确定，一般取3×3或5×5。

4.4.4 精度验证与质量评估

预处理完成后，必须进行独立的精度验证。验证方法包括：① 检查点法——使用GNSS实测点或ICESat-2激光测高点作为独立参考，计算DEM高程与参考值之差；② 交叉验证法——从原始数据中随机抽取10%的点作为验证集，其余90%用于插值，比较插值结果与验证集差异；③ 地形特征一致性检验——对比DEM提取的山脊线、山谷线与高分辨率影像或实测数据的一致性。精度验证结果应形成质量报告，作为导航数据库的组成部分随产品交付。

综上所述，数字地形模型是地形匹配导航系统的基石。数据源选择需权衡覆盖范围、分辨率、精度和成本；分辨率与精度需根据匹配算法需求和系统误差预算确定；坐标系与格式标准是数据互操作的基础；预处理流程则是将原始数据转化为可用导航数据库的关键环节。工程实践中，应建立“数据源评估→预处理→精度验证→数据库发布”的完整质量管理体系，确保地形数据满足导航任务要求。