4.1 成像几何模型：像素为什么歪了
4.1 成像几何模型：像素为什么歪了
本节摘要：影像几何畸变有三宗来源：传感器的姿态与成像方式、地球的自转与曲率、地形起伏的投影位移。本节拆解每宗畸变的机理，建立中心投影与推扫成像的几何模型，解释为什么严格成图必须做正射纠正、为什么这项工序离不开数字高程模型。
像素为什么歪了
拿到一幅零级影像，它的几何通常"歪"在三个层面，诊断思路也按这三层展开。第一层是传感器内部：镜头畸变让直线拍成弯曲，探测器阵列的像元尺寸与排列存在制造误差，这些参数在实验室定得极精、随时间缓慢漂移，属于"知根知底"的系统性畸变。第二层是平台与地球的运动：卫星姿态的微小抖动直接平移与旋转成像视线；推扫式成像每行时刻不同，几秒的成像窗口里地球自转让地面目标向东偏移，轨道前进方向与探测器排列方向的夹角还会造成影像的斜切——这些因素叠加，让原始影像的边缘常常比中心偏出几十个像元。第三层是地形起伏：中心投影下，山头顶点的成像位置由"投影中心—山顶"连线决定，山越高、偏离像主点越远，位置位移越大，一座八百米的山在影像边缘可偏出十余个像元；河谷与山脊的相对位置被系统性地扭曲，这就是"影像上量出来的山谷偏了"的原因。
三层畸变的处理策略不同。内部畸变靠实验室定标参数修正，属于"出厂配置"；平台与地球运动畸变靠精密轨道与姿态数据加地面控制点联合建模；地形畸变必须引入数字高程模型做正射纠正。三步都做完，影像才达到"每个像素都能对上地图"的正射级精度——这是与矢量数据叠加、多期影像对比、面积量算的共同前提。
图 4-1：地形起伏引起的投影位移
侧视雷达的几何是另一套账
第 3 章留下的 SAR 畸变伏笔在这里接上：SAR 成像是斜距投影，记录的是"回波先后"而不是"角度方向"，于是山地地形出现三种光学影像没有的畸变。透视收缩让朝向雷达的坡在影像上被压缩——山坡的真实长度在斜距上被"折叠"；坡度陡过雷达视角时，山顶回波比山脚先到，山体在影像上整体"倒向"传感器一侧，形成叠掩，影像上叠掩区的亮度因回波叠加而异常高；背向雷达的坡完全收不到回波，形成雷达阴影。处理思路上，SAR 正射纠正同样依赖数字高程模型，但模型要按斜距几何重算，模拟每个地形点的回波时刻再反推位置；业务上还要评估阴影区占比，阴影区没有信息，制图时必须标记出来而不是假装看见过。理解这套几何，还有一个现实收益：挑 SAR 数据时看入射角——陡峭地形选大入射角减少叠掩，平缓区域选小入射角换取更强的回波。
几何模型怎么用：从坐标变换到重采样
几何校正的数值流程可以概括为"建模 + 重采样"。建模阶段建立"影像坐标到地图坐标"的函数关系：有轨道与姿态参数的影像用严格成像模型（共线方程或多项式轨道模型），参数缺失的老影像退而用多项式或薄板样条，靠散布全幅的地面控制点拟合。控制点的质量是精度天花板——应选取道路交叉口、桥梁、田埂交角这类角点清晰、全年稳定的地物，均匀分布在整幅影像上，数量按模型自由度的三倍以上配置；控制点打在易变地物（树冠、水边、阴影边缘）上是新手最常犯的错。重采样阶段按模型把原图像素搬到目标格网，最近邻法保原值但锯齿明显，双线性与三次卷积更平滑但会轻微修改 DN 值——做分类与反演的数据建议最近邻或双线性，并在日志里记下所用方法，让下游知道像素值被动过哪一种手术。
# 几何校正的账面演算：控制点数量与残差的关系（最小二乘思想演示） import math # 一次仿射变换有 6 个参数，至少 3 个控制点；评估残差用多余的点 gcp_total, params_used = 12, 6 dof = gcp_total * 2 - params_used # 每个控制点提供 x、y 两个方程 print("自由度:", dof) rms_budget = {1.0: "可用：常规制图", 0.5: "良好：多期对比", 0.2: "严格：工程量算"} for r in rms_budget: print(f"目标 RMS {r} 像元（10m 数据 = {r*10:.0f} m） -> {rms_budget[r]}") # 变化检测的经验线：配准 RMS 应小于目标变化像元的 1/5， # 否则两期影像的"错位"会伪装成"变化"
常见问答
问：正射纠正与几何校正到底差在哪？ 几何校正把影像整体"摆正"到地图坐标，适合平坦地区；正射纠正额外引入数字高程模型逐点消除地形投影差，山区与高分辨率数据不可省。判断标准很简单：影像范围内的地形起伏导致的位移超过了你的精度预算，就必须正射。
问：控制点从哪里来？ 三个来源按优先级排：正射影像底图上自动匹配（效率最高）、外业实测的像控点（精度最高、成本最高）、旧版地形图与已校正影像（免费但注意时效）。混合使用时以高精度源为基准，避免"多套基准互相打架"。
问：重采样方法怎么选才不后悔？ 记一句口诀：分类反演用最近邻，目视制图用卷积。最近邻不改动原始 DN 值，保住定量合法性；三次卷积视觉最平滑，适合出图。混用的下场是定量结果带着重采样污染不自知——这类事故在质量事故档案里常年占据前排。
排错手记：几何问题的三种现场
现场一：镶嵌后道路错位成"台阶"。 两期影像配准残差超标，接缝两侧各按各的坐标走。处置：以高精度一景为基准重配准，复查控制点分布是否覆盖接缝带。
现场二：山区正射后影像边缘拉花。 数字高程模型分辨率远粗于影像，高程内插失真。处置：换更高分辨率高程数据，或把成图范围裁到高程数据可靠覆盖区内。
现场三：变化检测在所有高对比边界上报变化。 不是地物变了，是两期影像错位（第 5 章的伪变化清单头名）。处置：先做配准诊断——把两期影像做成假彩色叠加，边界泛红即错位证据；再按五分之一预算重配。几何问题的共同教训：查变化，先查配准。
要点回顾：几何畸变分传感器内部、平台与地球运动、地形起伏三层，处理手段分别是定标参数、严格模型加控制点、正射纠正；推扫影像是几十秒的接力拍摄，自转偏斜与姿态抖动是它的招牌畸变；SAR 斜距几何带来透视收缩、叠掩与阴影，正射处理与阴影标记是雷达数据的标准工序；控制点选角点、忌选易变地物，配准误差应控制在目标变化尺度的五分之一以内。下一章进入辐射账本——从 DN 值到地表反射率的完整换算。
附：几何工序速记卡
工序速记四条。一、三级精度对三级工序：粗纠（多项式）够目视，精纠（严格模型）够对比，正射（加高程）够量算——先定用途再定工序。二、控制点是天花板：再好的模型也救不了烂控制点，数量按自由度三倍起、分布按全覆盖、选点按"角点加稳定"。三、重采样定生死：定量流程最近邻，制图流程卷积，日志必须记。四、雷达是另一套账：斜距几何、地形模拟、阴影标记，套用光学流程必翻车。把四条贴在处理流水线的门口，每天进门看一眼——几何事故没有一起是知识不足造成的，全是纪律松懈造成的。
延伸：从畸变表读懂一颗卫星
每个数据产品的说明文档都会附几何精度指标（绝对定位精度、相对精度、波段配准精度），教你把它们读成"畸变表"：绝对定位精度反映轨道与控制点质量，决定影像能不能直接叠地图；相对精度（影像内部）反映姿态稳定度，决定内部量算可信度；波段间配准精度反映探测器对齐质量，差一像元以上时指数计算会在边缘发花。三个数字各自对应不同工序的输入要求——几何校正的精度预算不能超过绝对定位精度，指数计算要求波段配准在半像元内。读产品说明时把这三个数字圈出来，与任务预算对表，五分钟就能判断"这份数据配不配我的任务"。这份阅读能力不依赖任何工具，却能让新人在评审会上第一次发言就说到点上。