4. 几何标定基础

几何标定的目标，是让影像上的每一个像元都能被准确映射到地面位置。要做到这一点，必须先理解影像几何误差从哪里来，再用一套参数模型把这些误差描述出来，最后用地面控制点求解参数并评价精度。本节围绕这条主线展开。

一、几何畸变的来源

遥感影像的几何畸变可分为两大类。第一类是传感器内部引起的畸变：镜头径向畸变使像点沿径向偏移，切向畸变由镜头与焦平面不平行引起；CCD/CMOS 探元的排列间距误差、线阵推扫时的扫描非线性、多波段配准误差，都会让像点偏离理想位置。第二类是外部几何条件引起的畸变：平台姿态（俯仰、横滚、偏航）变化、轨道位置漂移、地球自转与曲率、地形起伏引起的投影差，以及侧视成像带来的透视变形。前者与传感器自身结构有关，后者与成像瞬间的位姿和地面条件有关。

区分这两类畸变很重要：内畸变在实验室可标定，外畸变需结合轨道与姿态数据或地面控制点解算。实际标定中，两者往往耦合在一起，需要统一建模。

二、内方位元素与外方位元素

内方位元素描述传感器内部的成像几何，包括主距 f（镜头中心到像平面的距离）、像主点坐标 (x₀, y₀)（主光轴与像平面的交点），以及镜头畸变系数（径向 k₁, k₂, k₃，切向 p₁, p₂）。它们决定了像空间坐标系中光线如何汇聚，是恢复成像光束形状的基础。

外方位元素描述成像瞬间传感器在物方空间中的位置和姿态，共 6 个：三个线元素 (Xₛ, Yₛ, Zₛ) 表示投影中心在地面坐标系中的坐标，三个角元素 (φ, ω, κ) 表示像平面相对于地面坐标系的旋转。线元素决定“从哪里拍”，角元素决定“朝哪个方向拍”。

内、外方位元素共同构成共线条件方程的基础。若内方位元素已知且外方位元素可解，就能建立像点与地面点的严格对应关系。

三、几何标定模型

最经典的模型是共线方程模型：地面点、投影中心、像点三点共线。其形式为

x - x₀ = -f · [a₁(X-Xₛ)+b₁(Y-Yₛ)+c₁(Z-Zₛ)] / [a₃(X-Xₛ)+b₃(Y-Yₛ)+c₃(Z-Zₛ)]
y - y₀ = -f · [a₂(X-Xₛ)+b₂(Y-Yₛ)+c₂(Z-Zₛ)] / [a₃(X-Xₛ)+b₃(Y-Yₛ)+c₃(Z-Zₛ)]

其中 aᵢ, bᵢ, cᵢ 是由三个角元素构成的旋转矩阵元素。对于推扫式线阵传感器，常采用有理多项式模型（RFM），用有理多项式拟合严格几何关系，形式统一、便于分发，但系数本身不具物理意义。此外还有仿射变换模型，适用于视场角小、地形平坦的场景，参数少、求解稳定。

模型选择取决于传感器类型、地面起伏和精度要求。严格模型精度最高但需要完整轨道姿态数据；RFM 通用性强但需足够控制点；仿射模型简单但适用范围有限。

四、地面控制点与几何精度评价

地面控制点（GCP）是已知地面坐标且在影像上可精确定位的点，用于求解标定参数。选取 GCP 应满足：分布均匀、覆盖影像四角与中心、地形代表性好、点位清晰易识别（如道路交叉口、建筑角点）。数量上，严格模型一般需要至少 6 个点，实际常用 15～30 个，并留出 20%～30% 作为检查点。

求解过程通常采用最小二乘平差：先由 GCP 建立误差方程，迭代解算内外方位元素或 RFM 系数，再用检查点评估精度。精度指标包括平面中误差（RMSE）和高程中误差，以及最大残差。若检查点 RMSE 明显大于 GCP 拟合残差，说明模型过拟合或控制点分布不合理。

几何精度评价还应关注误差的空间分布：若残差在影像边缘明显增大，可能是镜头畸变未充分建模；若沿扫描方向呈系统性偏移，可能是扫描非线性或时间同步问题。通过残差图可以直观判断误差来源，指导模型改进。

综上，几何标定是一个“理解误差来源—建立参数模型—利用控制点解算—评价并反馈优化”的闭环过程。内、外方位元素是模型的物理基础，GCP 是解算的输入，精度评价则是判断标定是否可用的依据。只有把这几环都做扎实，后续的影像几何校正和产品生产才有可靠保障。