4. 正射校正原理与农业地形适配
卫星或无人机在成像瞬间,传感器姿态、地形起伏和地球曲率都会让像点偏离其真实的地面正射位置。正射校正的任务,就是把这些"投影误差"逐像素地消除,生成具有统一比例尺、可直接量测的正射影像。农业应用对这一步尤其敏感:地块边界、垄沟走向、梯田坎线都需要在影像上保持几何真实,否则后续的面积统计、农机路径规划与长势分区都会出现系统性偏差。
一、共线方程:几何校正的物理基础
共线方程描述的是"投影中心—像点—地面点"三点共线这一成像几何关系。设地面点坐标为 (X, Y, Z),投影中心坐标为 (XS, YS, ZS),像点坐标为 (x, y),焦距为 f,则严格成像模型可写为:
x - x0 = -f · [a1(X-Xs) + b1(Y-Ys) + c1(Z-Zs)] / [a3(X-Xs) + b3(Y-Ys) + c3(Z-Zs)]
y - y0 = -f · [a2(X-Xs) + b2(Y-Ys) + c2(Z-Zs)] / [a3(X-Xs) + b3(Y-Ys) + c3(Z-Zs)]
其中 (x0, y0) 为主点偏移,ai、bi、ci 是由三个外方位角元素(俯仰 φ、横滚 ω、偏航 κ)构成的旋转矩阵元素。这组方程把传感器内方位(f、主点)与外方位(位置、姿态)同地面坐标联系起来,是摄影测量严格正射校正的核心。实际处理时,先由 POS 数据或空三解算给出初值,再用控制点做光束法平差,把外方位元素精化到足够精度,然后逐像素迭代求解其对应的地面坐标,完成重采样。
共线方程的优点是物理意义明确、可解释性强,缺点是必须已知内外方位元素,且对 POS 精度依赖较大。当缺少高质量 POS 或传感器参数不公开时,工程上常改用有理多项式模型。
二、有理多项式模型(RPC)
有理多项式模型用两个多项式之比来逼近严格成像几何,形式为:
Sample = P1(X, Y, Z) / P2(X, Y, Z)
Line = P3(X, Y, Z) / P4(X, Y, Z)
式中 P1~P4 为含 20 个系数的三元三次多项式,X、Y、Z 与 Sample、Line 均经过归一化处理。RPC 文件通常随高分卫星影像一同提供,用户无需知道传感器内部结构即可完成正射纠正。其精度在平坦地区可达亚像素级,但在高差剧烈区域,由于 RPC 本身是对严格模型的近似,需要引入 DEM 参与迭代计算地面点高程,否则山脊与沟谷处会出现明显的投影差残留。因此,无论采用共线方程还是 RPC,DEM 都是地形相关误差的决定性输入。
三、DEM 在平原、丘陵、梯田中的差异化应用
平原区高差通常在数米以内,投影差 Δr ≈ Δh·r/H 中 Δh 很小,即使使用 30 m 分辨率的 SRTM 或 ASTER GDEM 也能满足 1:10000 成图要求。此时 DEM 的主要作用是消除地球曲率与传感器姿态引起的低频误差,而非地形投影差。控制点可优先布设在影像四角与中心,形成包围式构架。
丘陵区坡向变化剧烈,同一景影像内迎坡面与背坡面的投影差方向相反。若 DEM 分辨率不足,山脊两侧会出现"错位"或"拉花"。建议采用 10~30 m 的 DEM,并在坡面中上部、坡脚各布设控制点,使平差能捕捉到高程方向的系统性偏差。以某丘陵区无人机作业为例(示例参数):航高 300 m、焦距 35 mm、地面分辨率 5 cm,当 Δh = 60 m 时,边缘像点投影差可达约 1.0 m,相当于 20 个像素,若不纠正将直接导致地块边界错位。
梯田区是农业遥感中最棘手的地形。田坎高差常在 1~3 m 内突变,而坎线宽度仅数十厘米。若 DEM 网格大于坎宽,坎顶与坎脚被平滑到同一高程,正射影像上会出现"坎线模糊"或"田面倾斜"。此时应使用 5 m 以内的高精度 DEM(如无人机倾斜摄影生成的 DSM 或激光雷达点云),并在坎顶、坎脚分别布设控制点,必要时对坎线附近做局部加密。需要提醒的是,DSM 包含植被冠层高度,用于正射纠正时可能把树冠"拉"到地面位置,农业地块内若作物较高,宜优先使用滤波后的 DTM。
四、控制点布设策略与精度要求
控制点(GCP)的质量直接决定正射成果的绝对定位精度。布设时应遵循以下原则:
- 数量:平原区每景不少于 6~9 个,丘陵区 12~15 个,梯田区建议 15 个以上并保证每个地形单元均有覆盖。
- 分布:采用"四角+中心+地形特征点"的包围式布局,避免全部集中在影像一侧;控制点应跨越不同高程带。
- 点位:选择田块拐角、道路交叉口、沟渠交汇等影像上清晰可辨、不随时间变化的地物,避免选在作物冠层顶部或阴影区。
- 精度:平面中误差应优于 0.5 个地面分辨率像素,高程中误差优于 1/3 等高距;检查点(CKP)应独立于控制点,数量不少于控制点的 20%。
平差完成后,用检查点计算 RMSE。以 5 cm 分辨率的无人机影像为例(示例判据):平面 RMSE 应控制在 2.5 cm 以内,最大残差不超过 5 cm;若某检查点残差明显偏大,应回查该点是否位于地形突变处或 DEM 误差较大的区域,必要时补充局部控制点重新平差。只有残差分布均匀、无明显方向性偏移,才能认为正射成果可用于地块量测与面积统计。
综上,正射校正不是简单的"套公式",而是共线方程/RPC、DEM 与控制点三者协同的过程。理解地形如何通过 DEM 影响投影差,再据此调整控制点布设密度与 DEM 分辨率,才能让正射影像在平原、丘陵、梯田等不同农业地形上都保持可靠的几何精度。