ENVI：RPC 正射校正：参考影像和 DEM 的选择，自带 DEM 为什么不推荐

一、RPC 模型与正射校正的数学基础

1.1 有理多项式模型的基本形式

RPC（Rational Polynomial Coefficients）模型用两个有理多项式分别拟合影像行、列坐标与地面经纬度、高程之间的映射关系。其标准形式为：

Line = P1(Lat, Lon, Height) / P2(Lat, Lon, Height)
Samp = P3(Lat, Lon, Height) / P4(Lat, Lon, Height)

其中 P 为三次多项式，含 20 个系数，
分母多项式通常归一化，使 P2、P4 的常数项为 1。

RPC 本质是对严格几何模型（共线方程）的高精度拟合，其优势在于用户无需掌握卫星轨道、姿态等敏感参数即可完成几何处理。本文评述：RPC 的“黑箱”特性既是便利也是陷阱——当输入高程存在系统性偏差时，误差会以近似线性的方式传递到平面坐标，而用户往往难以察觉。

1.2 正射校正的误差传播链

正射校正的核心是逐像元求解地面点高程，再依据 RPC 反算其对应的影像坐标并重采样。误差来源可分解为四类：

其中 DEM 垂直误差的平面影响与地形坡度正相关。在 30° 坡面上，10 m 高程误差可造成约 17 m 的平面位移（tan30°≈0.577，位移=高程误差/tan坡度）。这意味着在山区，DEM 选型的重要性甚至超过控制点采集精度。

1.3 有理多项式与严格模型的精度对比

根据 Grodecki 与 Dial（2003）的经典研究，在无地面控制点条件下，RPC 模型对 IKONOS 影像的拟合精度可达 0.1 像元以内。国内学者张力等人（2022）对资源三号、高分七号的测试表明，RPC 拟合残差普遍在 0.3 像元以下。本文评述：RPC 本身的拟合精度已足够支撑 1:10000 测图，工程中出现的精度问题绝大多数来自 DEM 与控制数据，而非 RPC 模型本身。

二、ENVI 两条正射校正技术路线的本质差异

2.1 RPC Orthorectification Workflow

该工作流仅依赖 RPC 文件与 DEM，通过自动匹配或手动采集地面控制点（GCP）优化 RPC 参数。其核心步骤包括：

• 1) 加载影像与 RPC 文件，自动读取元数据中的高程基准信息；
• 2) 选择 DEM 文件，设置重采样方法（推荐双线性或三次卷积）；
• 3) 可选：采集 GCP 进行 RPC 精化，通常 6–12 个均匀分布点即可；
• 4) 设置输出投影、像元大小与重采样方式，执行校正。

该路线的精度完全取决于 DEM 质量与 GCP 精度，适合有实测控制资料的生产项目。

2.2 RPC Orthorectification Using Reference Image

该路线引入参考影像（通常为已校正的正射影像）作为几何基准，通过自动匹配连接点反算 RPC 修正参数。其优势在于无需实测 GCP，适合无控区或境外区域。但风险同样明显：参考影像自身的几何误差会被“继承”到输出产品中。

本文评述：参考影像路线的本质是“相对几何一致性”而非“绝对几何精度”。若项目要求绝对定位精度优于 5 m，必须引入外部控制或高精度参考底图，否则只是把误差从一个坐标系搬到另一个坐标系。

2.3 两条路线的适用场景对比

三、自带 DEM 的精度真相：从 SRTM 到 ASTER GDEM

3.1 常见全球 DEM 的精度指标

ENVI 自带或可自动下载的 DEM 主要包括 SRTM、ASTER GDEM、ALOS World 3D 等。其公开精度指标如下：

上述数据来源为各机构官方精度报告，经笔者整理。需要强调的是，这些指标为全球统计值，在陡峭山区、植被密集区或城市高层建筑区，实际误差可能显著放大。

3.2 自带 DEM 不推荐的四个技术原因

原因一：垂直精度不足以支撑大比例尺测图。以 SRTM 为例，其 6–16 m 的垂直误差在 1:10000 测图中对应平面误差可达 10–28 m（按 30° 坡面估算），远超规范要求的 5 m 中误差。本文评述：很多用户抱怨“校正后影像偏移”，根源往往就在此处。

原因二：高程基准不统一。SRTM 采用 EGM96 大地水准面，而国产卫星 RPC 通常基于 WGS84 椭球高。两者在中国区域差异可达 10–30 m。若未进行基准转换，相当于给所有像元叠加了一个系统性高程偏差。

原因三：数据时效性差。SRTM 数据采集于 2000 年，ASTER GDEM 主要基于 2000–2013 年数据。对于近十年快速城市化的区域，地表高程已发生显著变化，用旧 DEM 校正新影像会引入“伪位移”。

原因四：分辨率与影像不匹配。30 m DEM 对 0.5–2 m 分辨率的国产高分影像而言过于粗糙，地形细节丢失严重。在山区，DEM 网格间的线性插值会平滑掉真实地形起伏，导致校正后影像出现“阶梯状”畸变。

3.3 什么情况下自带 DEM 勉强可用

并非所有项目都需要高精度 DEM。以下场景可考虑使用自带 DEM：

• 1) 平原地区（坡度 < 5°），高程误差对平面精度影响极小；
• 2) 1:50000 及更小比例尺的概查、普查任务；
• 3) 仅需相对几何一致性、不要求绝对定位精度的变化检测；
• 4) 作为快速预览或算法验证的临时数据。

但即便如此，仍建议进行高程基准转换，并在精度报告中明确说明 DEM 来源与已知局限。

四、参考影像的选择准则：几何、辐射与时效

4.1 几何精度：参考影像的“天花板”效应

参考影像的几何精度决定了输出产品的精度上限。若参考影像本身存在 5 m 中误差，无论后续匹配多精确，输出产品的绝对精度都不可能优于 5 m。本文评述：这是参考影像路线最容易被忽视的“天花板效应”。

选择参考影像时，应优先考虑：

• 1) 已有实测控制点校正的正射影像产品；
• 2) 公开的全国 DOM 成果（如天地图、地理国情监测成果）；
• 3) 同源或同轨道的高精度正射影像。

4.2 辐射一致性：匹配成功率的隐形门槛

ENVI 的自动匹配算法对辐射差异敏感。若参考影像与待校正影像的成像季节、太阳高度角、传感器波段差异过大，匹配点数量会急剧下降。根据笔者经验，以下情况匹配成功率通常低于 30%：

4.3 时效性：变化区域的“伪匹配”陷阱

在快速变化区域（如新建城区、填海造陆区），参考影像与待校正影像的地物可能完全不同。此时自动匹配算法可能找到错误的同名点，导致局部几何畸变。本文评述：时效性差异带来的不是随机误差，而是系统性粗差，必须通过人工检查剔除。

五、DEM 选型：分辨率、垂直精度与基准匹配

5.1 分辨率匹配原则

DEM 分辨率与影像分辨率的匹配关系直接影响校正精度。经验准则如下：

• 1) 1:10000 测图：DEM 分辨率应优于 5 m，垂直精度优于 2 m；
• 2) 1:50000 测图：DEM 分辨率应优于 10 m，垂直精度优于 5 m；
• 3) 1:250000 及更小：30 m DEM 可接受。

上述准则参考了国家测绘行业标准 CH/T 3007-2011 与笔者工程实践。需要说明的是，在陡峭山区，即使满足上述指标，仍建议使用更高精度 DEM。

5.2 垂直精度与坡度的影响

DEM 垂直误差对平面精度的影响可用下式近似估算：

平面位移 ≈ 高程误差 / tan(坡度)

示例：高程误差 10 m，坡度 30°
平面位移 ≈ 10 / 0.577 ≈ 17.3 m

该公式为简化估算，实际误差还受坡向、卫星侧视角等因素影响。但足以说明：在山区，DEM 垂直精度每提升 1 m，平面精度可改善 1.7 m 以上。

5.3 高程基准转换：最容易被忽略的一步

国产卫星 RPC 通常基于 WGS84 椭球高，而多数 DEM 采用 EGM96 或 EGM2008 大地水准面高。两者差异在中国区域可达 10–30 m。ENVI 中可通过以下步骤进行转换：

• 1) 使用 ENVI 的 Geoid 工具或外部 EGM2008 模型计算大地水准面差距；
• 2) 将 DEM 高程加上大地水准面差距，转换为椭球高；
• 3) 或在 RPC 精化时通过 GCP 自动吸收该系统性偏差。

本文评述：基准转换是正射校正中最容易被跳过、也最容易导致系统性误差的步骤。建议在项目开始时即建立基准转换清单，明确各数据源的基准信息。

六、工程操作路径：从控制点到精度报告

6.1 数据准备清单

6.2 控制点采集与 RPC 精化

控制点采集是正射校正中最耗时的环节，也是精度提升最显著的环节。推荐流程：

• 1) 在影像四角与中心区域均匀选取 9–15 个控制点；
• 2) 优先选择道路交叉口、建筑物角点等特征明显、不易变化的地物；
• 3) 控制点坐标应来自实测 GPS 或高精度正射影像；
• 4) 在 ENVI 中执行 RPC Refinement，检查残差是否小于 1 像元；
• 5) 剔除残差过大的点，迭代优化。

本文评述：控制点数量并非越多越好，均匀分布比数量更重要。笔者经验表明，9 个均匀分布的控制点通常优于 30 个集中在某一区域的控制点。

6.3 DEM 重采样与裁剪

ENVI 支持多种 DEM 重采样方法，推荐使用双线性插值或三次卷积。对于大区域，建议先裁剪 DEM 至影像覆盖范围，以减少处理时间。操作步骤：

File → Open → DEM
Basic Tools → Resize Data → Spatial Subset
选择影像范围 → 输出裁剪后 DEM
Basic Tools → Layer Stacking → 可选：与影像叠加检查

6.4 执行正射校正

在 ENVI 中执行 RPC Orthorectification Workflow 的关键参数设置：

6.5 精度验证与报告

校正完成后，必须进行精度验证。推荐使用检查点法：

• 1) 独立采集 20–30 个检查点，不参与 RPC 精化；
• 2) 计算平面中误差 RMS = sqrt(Σ(Δx²+Δy²)/n)；
• 3) 绘制误差分布图，检查是否存在系统性偏移；
• 4) 输出精度报告，包含 DEM 来源、参考影像来源、控制点数量与残差。

七、精度验证与常见问题排查

7.1 常见问题与解决方案

7.2 精度验证的量化指标

根据国家测绘行业标准 CH/T 3007-2011，1:10000 正射影像的平面位置中误差应优于 5 m（平地）或 7.5 m（山地）。本文评述：实际工程中，建议将内控标准提高 30%，即平地 3.5 m、山地 5 m，以留出足够的质量余量。

八、前沿进展与趋势预判

8.1 深度学习辅助的 DEM 超分辨率

近三年，基于深度学习的 DEM 超分辨率重建成为研究热点。例如，Zhang 等人（2023）提出的 SRGAN-DEM 方法可将 30 m DEM 提升至 10 m 等效精度。本文评述：这类方法在纹理丰富的山区效果显著，但在平原地区可能引入“伪地形”，需谨慎使用。

8.2 无控正射校正的自动化趋势

随着参考影像库的完善与匹配算法的进步，无控正射校正的精度已从早期的 10–20 m 提升至 3–5 m。ESA 的 Copernicus DEM 与 Sentinel-2 正射影像的组合，为全球无控校正提供了新的可能。

8.3 多源数据融合的几何一致性

未来正射校正将不再局限于单一影像，而是多源数据的几何一致性处理。本文评述：这要求工程师不仅掌握 ENVI 操作，还需理解不同传感器、不同基准之间的转换关系。

九、结论与选型速查表

9.1 核心结论

• 1) 自带 DEM（SRTM/ASTER GDEM）仅适用于平原、小比例尺或预览场景；
• 2) 1:10000 及以上测图必须使用优于 5 m 分辨率、2 m 垂直精度的 DEM；
• 3) 参考影像的几何精度决定了输出产品的精度上限；
• 4) 高程基准转换是必做步骤，不可跳过；
• 5) 控制点均匀分布比数量更重要，9–15 个即可。

9.2 选型速查表

9.3 拓展学习资源

• ENVI 官方教程：RPC Orthorectification Workflow
• USGS DEM 下载：EarthExplorer
• Copernicus DEM：ESA Copernicus DEM
• 视频教程：YouTube ENVI RPC 正射校正

主要参考文献

1. Grodecki J, Dial G. Block adjustment of high-resolution satellite images described by rational polynomials. Photogrammetric Engineering & Remote Sensing, 2003, 69(1): 59–68.
2. 张力, 张继贤, 等. 高分七号卫星立体影像几何精度评估. 测绘学报, 2022, 51(6): 889–899.
3. Zhang Y, et al. DEM super-resolution using generative adversarial networks for orthorectification. ISPRS Journal of Photogrammetry and Remote Sensing, 2023, 195: 123–138.
4. NASA JPL. SRTM Version 3.0 Global DEM Validation Report. 2013.
5. NASA/METI. ASTER GDEM Version 3 Validation Report. 2019.
6. ESA. Copernicus DEM Product Handbook. 2022.
7. 国家测绘局. CH/T 3007-2011 数字航空摄影测量 测图规范. 2011.
8. Reinartz P, et al. Orthorectification of satellite imagery using reference images. Remote Sensing, 2021, 13(4): 678.
9. 李德仁, 王密, 等. 高分辨率光学卫星遥感影像高精度几何处理. 武汉大学学报·信息科学版, 2023, 48(1): 1–12.