3. ENVI数据准备与预处理

蚀变信息提取本质上是在反射率空间中寻找与羟基、铁染等蚀变矿物相关的微弱光谱特征。如果输入数据仍是原始DN值，传感器响应差异、太阳高度角、地形阴影和大气散射都会叠加在像元亮度上，导致蚀变异常无法与背景区分。因此，本节的任务是把一景原始遥感影像加工成“物理量明确、空间位置准确、范围聚焦”的反射率数据，为后续波段比值、主成分分析和异常阈值分割提供可靠输入。整个流程可以概括为：辐射定标把DN值转换为辐亮度或表观反射率，大气校正把表观反射率还原为地表反射率，几何精校正把像元位置配准到统一地理坐标系，裁剪则把研究区之外的冗余数据剔除。

一、辐射定标：从DN值到辐亮度

辐射定标是后续所有定量处理的基础。以Landsat 8 OLI为例，ENVI中通常使用Radiometric Calibration工具，输入多光谱波段文件，选择Radiance（辐亮度）输出类型，并勾选Apply FLAASH Settings，使输出量纲为µW/(cm²·sr·nm)，这正是FLAASH大气校正所要求的输入单位。定标的核心公式为：

Lλ = Gain × DN + Offset
ρ_TOA = (π × Lλ × d²) / (ESUNλ × cosθs)

其中Gain和Offset来自影像元数据中的RADIANCE_MULT_BAND和RADIANCE_ADD_BAND字段，d为日地距离，θs为太阳天顶角。实际操作中不必手工计算，但需要理解两点：一是定标后的辐亮度仍包含大气影响，不能直接用于蚀变提取；二是定标输出应为浮点型，若误选整型输出会造成反射率量化损失。对于Sentinel-2 L2A产品，若已包含地表反射率，则可跳过大气校正，但仍需检查是否已做辐射定标和几何配准。

二、大气校正：FLAASH与QUAC的选择

大气校正是预处理中最影响蚀变提取效果的一步。ENVI提供FLAASH和QUAC两种常用工具。FLAASH基于MODTRAN辐射传输模型，需要输入传感器类型、飞行高度、像元大小、成像时间、大气模型、气溶胶模型和能见度等参数，输出地表反射率，精度较高但参数较多；QUAC（Quick Atmospheric Correction）是一种经验性方法，仅需输入辐亮度数据，自动估计大气参数，速度快但精度略低，适合缺乏详细气象数据或快速筛查的场景。

以FLAASH为例，关键参数设置如下：传感器类型选择Landsat-8 OLI；地面高程根据研究区DEM平均值填写，例如高原地区可填2000 m；大气模型按成像月份和纬度选择，如中纬度夏季选择Mid-Latitude Summer；气溶胶模型在植被覆盖区选Rural，城市周边选Urban；初始能见度可填40 km，若影像雾霾较重可降至20 km。FLAASH运行后会输出反射率影像和一个日志文件，日志中会给出各波段的平均反射率估计值，可用于检查校正是否合理。需要强调的是，FLAASH对短波红外波段（如Landsat 8的B6、B7）的校正效果直接影响羟基蚀变信息提取，因为这些波段是识别Al-OH吸收特征的关键窗口。

三、几何精校正与裁剪

几何精校正的目的是消除影像的几何形变，使蚀变异常位置与地质图、矿权边界准确套合。ENVI中可使用Image-to-Map Registration工具，以已有正射影像或地形图为基准，选取不少于15个均匀分布的控制点，采用二次多项式或Delaunay三角网校正模型，重采样方法选择Cubic Convolution，输出像元大小与基准影像一致。控制点残差（RMS）应控制在1个像元以内，若超过则需剔除误差较大的点重新选点。对于Landsat 8 L1T产品，几何精度通常已较好，但在山区或地形起伏大的区域仍需做精校正，否则蚀变异常会沿地形方向偏移。

裁剪环节使用Subset Data via ROIs或Resize Data工具，以研究区矢量边界或矩形范围裁剪校正后的反射率影像。裁剪时应保留足够的边缘缓冲区，避免后续滤波和主成分分析在边界处产生伪异常。输出文件建议采用ENVI标准格式（.dat），并记录波段顺序、波长信息和反射率缩放系数（如FLAASH输出通常乘以10000存储为整型），以便后续蚀变提取时正确还原反射率量纲。

四、结果检查与常见问题

预处理完成后不能直接进入蚀变提取，应先做三项检查。第一，光谱合理性检查：在影像中选取水体、浓密植被和裸土三类典型地物，查看其光谱曲线。水体在近红外波段反射率应接近0，植被在0.8 µm附近应有明显反射峰，裸土曲线应平滑无异常尖峰。第二，波段间相关性检查：若短波红外波段与可见光波段相关性异常高，可能说明大气校正不充分或定标单位错误。第三，空间套合检查：将反射率影像与已知地质界线或矿权边界叠加，观察是否出现系统性偏移。

常见问题中，最典型的是FLAASH报错“内存不足”或“波长文件缺失”。前者可通过分块处理或降低空间分辨率解决，后者需确认定标输出时已正确写入波长信息。另一个常见问题是裁剪后影像出现黑边，这通常是因为裁剪范围超出了校正后影像的有效区域，应在裁剪前先用统计工具确认影像有效范围。此外，若后续蚀变提取使用波段比值法，需确保参与比值的波段均已正确校正，否则比值结果会被大气程辐射放大。

完成上述步骤后，得到的反射率数据应满足：空间位置与基准影像套合误差小于1个像元，典型地物光谱曲线形态合理，短波红外波段保留矿物吸收特征，研究区范围完整且无系统性黑边。这份数据即可作为后续波段比值、主成分分析和蚀变异常圈定的标准输入。