3. 遥感影像预处理与波段配准

波段比值（如铁染指数 B4/B2、羟基指数 B5/B7，或 ASTER 的 B4/B6、B5/B8）本质上是一个逐像元的除法运算。除法对分子分母的误差极其敏感：若两个波段之间存在辐射量级差异、几何错位或大气影响不一致，比值结果就会出现虚假的条带、边缘晕圈和"伪异常"。因此，在进入比值计算之前，必须完成辐射定标、大气校正、几何精校正与波段间配准四项工作，使各波段在辐射尺度、空间位置、成像时刻状态三个维度上保持一致。

一、辐射定标：把DN值转换为有物理意义的辐亮度或反射率

原始影像记录的是无量纲的DN值，受传感器增益、积分时间、偏移量影响，不同波段、不同景之间的DN值不可直接比较。辐射定标的任务是将其转换为辐亮度（W·m⁻²·sr⁻¹·μm⁻¹）或表观反射率。以Landsat 8 OLI为例，常用定标公式为：

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

ML、AL：定标系数（元数据RADIANCE_MULT_BAND_x / RADIANCE_ADD_BAND_x）
Qcal：像元DN值；d：日地距离（天文单位）；θs：太阳天顶角
ESUNλ：波段大气外太阳辐照度

实操中应逐波段读取元数据中的增益与偏移，避免用统一系数处理所有波段。定标后需检查：各波段辐亮度是否落在合理量级（如OLI可见光波段典型值在10～100 W·m⁻²·sr⁻¹·μm⁻¹），有无负值或饱和像元。若使用ENVI的Radiometric Calibration工具，应选择"Radiance"或"Reflectance"输出，并确认输出数据类型为浮点型，防止后续比值运算被整型截断。

二、大气校正：消除波段间大气影响的不一致性

大气散射在短波波段（蓝、绿）贡献大，在短波红外贡献小。若不做大气校正，B2/B4这类跨可见光—近红外的比值会把气溶胶差异误判为地物异常。常用方法是FLAASH或6S模型，输入参数包括成像时间、太阳几何、气溶胶模型、能见度、地面高程。以FLAASH为例，关键设置：气溶胶模型按区域选择（如大陆型/城市型），初始能见度可设为40 km，水汽反演波段选940 nm附近。校正后应输出地表反射率，并检查植被、水体、裸地三类典型地物的光谱曲线是否符合常识——植被在近红外应明显高于红光，水体在近红外应接近0。

三、几何精校正与波段间配准：让同一像元对应同一地面位置

几何精校正解决整景影像与地理坐标系的对应关系，波段间配准则解决同一景内各波段之间的相对错位。二者不能互相替代：即使整景已正射校正，若波段间存在0.5～1个像元的残余错位，比值图在边缘处仍会出现明显"描边"效应。配准精度要求通常为优于0.5个像元，对高分辨率数据建议优于0.3像元。

操作要点：以全色波段或信噪比最高的波段作为参考波段，其余波段用二次多项式或局部三角网纠正；控制点应选在道路交叉、建筑角点、水体边界等稳定特征上，每景不少于20～30个，且均匀分布。配准后必须做残差检查：查看控制点RMS误差，并用"闪烁对比"或差值图目视检查边缘是否重影。下表给出一个可执行的配准质量检查记录示例（示例数据，非实测）：

四、一致性核验：进入比值计算前的最后一道关

四项处理完成后，还需做一次整体一致性核验：各波段是否同为反射率产品、空间分辨率与像元对齐是否一致、有效值范围是否统一（如统一为0～1浮点）、NoData值是否一致。建议在比值计算前先做一次波段间散点图，观察两波段是否存在合理的线性相关；若散点呈明显分层或截断，往往说明定标或校正环节存在波段间不一致。

上图把四项处理与一致性核验串成一条链，并标出了不合格时的回退路径。需要强调的是，回退不是失败：配准残差超标时回到控制点选取环节重新布点，往往比在比值图上做后处理更有效。实际项目中，建议把上表的检查项固化为处理记录，随成果一并归档，便于后续异常图的可追溯性。

完成上述步骤后，各波段已处于统一的辐射尺度与空间框架下，波段比值计算才具备可比性。此时再进入比值运算与异常分级，得到的异常图才真正反映地物光谱差异，而非预处理残留误差。