4. 遥感影像预处理与增强

用地变化判读的本质，是在两个时相（或多时相）影像上识别同一地物是否发生了类型转换。判读结论是否可靠，很大程度上不取决于判读员的经验，而取决于进入判读环节的影像是否"可比"。原始遥感影像受传感器响应差异、大气散射与吸收、地形起伏、成像几何、云量以及分辨率限制等因素影响，直接拿两期原始影像做对比，很容易把辐射差异误判为地表变化。因此，预处理与增强的目标只有一个：让影像在辐射、几何、空间分辨率三个维度上具备可比性，并让目标地物在视觉上更易辨识。

一、辐射定标：把DN值还原为物理量

传感器记录的原始像元值（DN值）是无量纲的整数，不同传感器、不同增益设置下同一地物的DN值并不相同。辐射定标的作用是建立DN值与入瞳辐亮度或表观反射率之间的定量关系。以常见的Landsat 8 OLI为例，其定标公式为：

Lλ = ML × Qcal + AL
ρλ' = (π × Lλ × d²) / (ESUNλ × cosθs)

其中 ML、AL 为乘性、加性定标系数（来自元数据MTL文件），Qcal 为量化DN值，d 为日地距离（天文单位），ESUNλ 为大气顶太阳辐照度，θs 为太阳天顶角。定标后得到的是大气顶表观反射率（TOA），它已经消除了传感器增益差异，但尚未消除大气影响。这一步是后续大气校正的输入前提，不能跳过。

二、大气校正：消除"同物异谱"的辐射干扰

大气校正把TOA反射率转换为地表反射率。大气散射会使暗地物（如水体、阴影）偏亮，气溶胶和水汽吸收会削弱特定波段信号，两期影像成像时大气条件不同，就会造成同一地物反射率出现系统性偏移。工程上常用FLAASH、6S、Sen2Cor或暗像元法（DOS）完成校正。以DOS法为例，其核心假设是影像中存在反射率接近0的暗地物（如清洁深水体），将该暗像元在各波段的最小值视为大气程辐射，逐波段减去即可：

ρ_surface(i,j) = ρ_TOA(i,j) − ρ_path(i)
ρ_path(i) ≈ min[ρ_TOA(i)] − 0.01

DOS法简单快速，适合快速排查，但对气溶胶空间变化剧烈或缺少暗地物的场景精度有限；涉及定量反演或精细变化检测时，应优先采用基于辐射传输模型的FLAASH或6S。判断校正是否合理，可检查植被在近红外波段的反射率是否明显高于红光波段、水体在近红外是否接近0，若出现异常高值则需回查参数。

三、几何校正与影像配准：让像元"对得上"

几何校正解决两件事：一是消除传感器姿态、地形、地球曲率造成的几何畸变（正射校正）；二是让两期影像的同一地物落在同一像元位置上（影像间配准）。用地变化判读对配准精度要求较高，一般要求均方根误差（RMSE）控制在0.5个像元以内，否则变化检测会在边界处产生大量伪变化。实际操作中，以一期影像为基准，在另一期上选取不少于20～30个分布均匀、特征清晰的控制点（道路交叉口、堤坝拐角、独立建筑物角点），采用二次多项式或有理函数模型进行校正，重采样优先选用双线性或三次卷积，保持光谱信息连续性。

四、影像融合与裁剪：兼顾空间细节与光谱信息

当全色影像空间分辨率高于多光谱影像时（如全色15 m、多光谱30 m），可通过融合（Pan-sharpening）获得兼具高空间细节与多光谱信息的影像，常用方法有Gram-Schmidt、Brovey、NNDiffuse等。融合后应检查光谱保真度，避免出现明显色偏。裁剪则按评价范围或行政区边界进行，建议在评价范围外扩一定缓冲距离（如500 m～1 km）再裁剪，避免边缘效应影响判读。裁剪后统一投影坐标系与像元大小，确保两期影像严格对齐。

五、影像增强：让变化"看得见"

增强处理不改变地物的物理属性，而是改善目视判读效果。常用手段包括：线性拉伸与直方图均衡化提升整体对比度；波段组合（如标准假彩色 5-4-3 突出植被、4-5-3 突出水体与建设用地）强化目标地物；主成分分析（PCA）与缨帽变换（TCT）压缩冗余信息、突出亮度—绿度—湿度差异；对两期影像做差值或比值运算生成变化增强图，使新增建设用地、植被减少等变化在灰度上直接显现。增强参数应记录在案，保证判读过程可复现。

六、结果判断与常见问题

预处理是否到位，可以用几个可观察的指标来检验：两期影像同名地物（如硬化路面、深水体）的地表反射率差异应控制在合理范围内；配准残差在目视检查中不应出现明显的"重影"或边界错位；融合影像在植被与水体交界处不应出现明显的光谱失真。若判读阶段发现大量沿道路、地块边界分布的细碎"变化图斑"，通常提示配准精度不足；若整幅影像出现系统性亮度偏移，则多与大气校正缺失或定标系数误用有关，应回到对应环节复核，而不是在判读阶段强行解释。

综上，辐射定标与大气校正保证"同物同谱"，几何校正与配准保证"同地同元"，融合与裁剪保证"同尺度同范围"，增强处理则让变化在视觉上凸显。四类处理环环相扣，任何一环缺失或参数不当，都会以伪变化的形式传导到最终的用地变化判读结论中。