4. 遥感数据预处理实战

铁矿区多位于基岩裸露、地形起伏较大的山地丘陵地带，遥感影像上既有大面积铁染蚀变异常，也有植被、水体、云影和地形阴影的干扰。要把原始 DN 值影像变成可用于提取铁染异常、羟基异常和线性构造的定量底图，预处理必须按“辐射定标—大气校正—几何精校正—影像裁剪”的顺序逐级推进。这一顺序不是形式上的排列：辐射定标把无量纲 DN 值转换为辐亮度或表观反射率，是大气校正的输入前提；大气校正得到地表反射率，才能支撑后续的铁染指数（如 B4/B3、B2/B1 比值）和主成分分析；几何精校正保证异常位置能与已知矿点、地质图套合；裁剪则把数据量压缩到矿区范围，减少无效计算和边缘畸变影响。下面以一幅 Landsat 8 OLI 铁矿区影像（示例假设为 2023 年 10 月、行列号 128/039、云量约 8%）为例，说明各环节的操作要点与判据。

辐射定标：从 DN 到辐亮度与表观反射率

辐射定标的实质是用定标系数把传感器记录的整型 DN 值还原为物理量。对 Landsat 8 OLI，定标公式为：

Lλ = ML × Qcal + AL
ρλ' = (Mρ × Qcal + Aρ) / sin(θSE)

其中 Lλ 为辐亮度（W·m⁻²·sr⁻¹·μm⁻¹），Qcal 为像元 DN 值，ML、AL 为辐亮度乘性、加性定标系数，Mρ、Aρ 为反射率定标系数，θSE 为太阳高度角。实际操作中，ML、AL、Mρ、Aρ 和太阳高度角均从影像元数据 MTL.txt 中读取，例如示例中某波段 ML=2.0×10⁻²、AL=-100，太阳高度角约 45°，则 sin(θSE)≈0.707。定标后应检查：辐亮度值是否落在传感器动态范围内、有无负值或异常饱和像元；表观反射率是否大致在 0～1 之间。若出现大面积负值，通常是定标系数读取错误或对 DN 值做了不当的偏移处理。需要强调，辐射定标只消除传感器响应差异，并未去除大气影响，因此不能直接用于铁染指数计算。

大气校正：获取地表反射率

大气校正的目标是去除瑞利散射、气溶胶散射和水汽吸收，把表观反射率转换为地表反射率。铁矿区找矿信息提取对短波红外和可见光波段比值敏感，大气程辐射会显著压低铁染指数对比度，因此这一步不可省略。工程上常用 FLAASH 或 6S 模型，关键输入包括：影像中心经纬度、成像时间、传感器类型、大气模式（示例取中纬度秋季）、气溶胶模式（矿区周边有裸露地表，可取大陆型）、能见度（示例假设 20 km）以及气溶胶反演波段。校正后应做两项检查：一是选取影像中的深清水体或浓密植被作为暗像元，其地表反射率应接近 0～0.05；二是对比校正前后同一铁染蚀变点的波段曲线，校正后应更接近典型铁氧化物光谱形态（可见光强、近红外相对低）。若校正后影像出现整体偏蓝或反射率普遍大于 1，多为气溶胶模式或能见度设置不当，应回退调整。

几何精校正：让异常落到正确的地面位置

Landsat 8 的 L1 级产品已做过系统几何校正，但铁矿区地形起伏大，山脊和沟谷处仍存在明显的地形位移，直接叠加地质图会出现异常“漂移”。几何精校正以已有正射影像或 1:5 万地形图为基准，采用二次多项式或有理函数模型，控制点应优先选在道路交叉口、水库坝体、矿山平台边缘等稳定且易识别的地物上，避免选在植被边界和阴影边缘。示例中在矿区范围内布设 25 个控制点，均匀分布于四角和中部，要求控制点残差（RMS）小于 0.5 个像元，总误差控制在 1 个像元以内。重采样方式对找矿信息提取有直接影响：计算比值指数前宜用双线性或三次卷积，保持光谱连续性；若后续要做分类或变化检测，则用最近邻以保留原始光谱值。校正完成后，用未参与校正的检查点复核，并叠加已知矿点，确认矿点落在异常晕圈内部而非边缘。

影像裁剪：限定矿区分析范围

裁剪不是简单“切图”，而要兼顾分析范围与边缘质量。建议按矿区地质单元外扩 2～3 km 建立感兴趣区，既覆盖蚀变晕和控矿构造，又避免把大量非矿区地物带入统计。裁剪方式上，按掩膜裁剪可保留不规则边界，按矩形范围裁剪便于后续分块处理；无论哪种，都应确保裁剪边界不穿过主要异常区，否则边缘像元会因重采样和混合像元影响指数计算。裁剪后应核对：影像行列数、投影信息、像元大小是否与原始一致，NoData 值是否设置正确，避免后续统计把填充值当作有效反射率。

操作闭环与常见问题

把四个环节串起来看，预处理的质量控制应形成闭环：每一步的输出既是下一步的输入，也是上一步的检验依据。例如大气校正后若发现铁染蚀变点反射率曲线异常，应回查辐射定标是否用了正确的太阳高度角；几何精校正后若矿点系统性偏移，应检查控制点是否集中在影像一侧导致外推误差。常见问题中，最容易被忽视的是重采样方式与指数计算不匹配——先做比值再重采样，与先重采样再算比值，结果在异常边界处会有差异，工程上应统一为“先重采样到目标分辨率，再计算指数”。另一个高频问题是裁剪时未统一 NoData 值，导致后续统计把填充值计入均值，压低异常对比度。建议在预处理结束后输出一份简短的质检记录，列明各环节参数、控制点残差、暗像元检验结果和裁剪范围，作为后续异常提取的可追溯依据。