4. 遥感蚀变异常提取方法

蚀变异常提取是遥感找矿的核心环节。其基本依据是：热液蚀变过程中，含铁矿物（褐铁矿、赤铁矿、黄钾铁矾等）和含羟基矿物（高岭石、绢云母、绿泥石、明矾石等）会在可见光—短波红外波段产生特征性吸收与反射。铁染矿物在0.45～0.55 μm附近有吸收、在0.85～1.0 μm附近反射增强；羟基矿物则在2.2 μm附近形成强吸收谷。提取方法要做的，就是把这些微弱的光谱特征从背景岩性、植被、阴影中分离出来，形成可叠加分析的铁染、羟基异常图层。

一、主成分分析（PCA）法：用统计正交性分离蚀变信息

PCA的原理是把多波段数据投影到新的正交主成分空间，使方差最大的方向成为PC1，依次递减。对蚀变信息而言，关键不是“做不做PCA”，而是判断哪个主成分承载了目标矿物的光谱特征。以Landsat 8 OLI为例，铁染常用波段组合为B2（蓝）、B4（红）、B5（近红外）、B6（短波红外1）；羟基常用B5、B6、B7（短波红外2）加B2或B4。

判断依据是主成分载荷矩阵：铁染异常通常表现为B2与B4载荷符号相反、B5与B6载荷符号相反，且B6载荷较大；羟基异常则表现为B6与B7载荷符号相反，B7载荷显著。选定主成分后，还要判断异常方向——若目标矿物在该主成分上为低值，需对结果取反，才能得到“高值=强蚀变”的图层。

实际操作流程为：辐射定标与大气校正（建议FLAASH或QUAC）→ 掩膜水体、云、雪、浓密植被 → 按组合波段做PCA → 读取载荷矩阵确定目标主成分 → 按载荷符号决定是否取反 → 用均值±k倍标准差（常用k=2或2.5）做阈值分割 → 得到分级异常图层。示例参数：某景OLI数据经PCA后，PC4中B6载荷为+0.71、B7为−0.63，判定PC4为羟基主成分，取反后以均值+2σ为一级异常、+2.5σ为二级异常。

二、波段比值法：用比值压制地形与照度差异

比值法通过波段相除突出吸收—反射反差，同时压制地形阴影和整体照度差异。常用组合包括：铁染用B4/B2（红/蓝）或B5/B4；羟基用B6/B7（SWIR1/SWIR2）或B7/B6。比值结果接近1表示无明显吸收，偏离越大异常越强。

比值法的短板是容易受植被和云影干扰，因此常配合“比值+阈值”或“比值+主成分”使用。实操中建议先做比值，再用直方图观察分布：若直方图呈明显右偏长尾，说明存在高值异常群；以均值+2σ或累积频率97.5%分位作为切割阈值，可得到较稳定的异常范围。示例：B6/B7比值图中，高岭石化区比值普遍大于1.35，取1.35为下限可圈出羟基异常带。

三、光谱角填图（SAM）：用光谱形状匹配识别蚀变矿物

SAM把每个像元光谱视为n维向量，计算其与参考光谱之间的夹角，夹角越小越相似。它只比较光谱形状，不受整体亮度影响，因此对反照率差异不敏感，适合识别特定蚀变矿物组合。使用前提是数据已完成大气校正，并准备好端元参考光谱——可来自USGS光谱库、野外实测光谱或影像中纯净像元。

操作要点：将参考光谱重采样到影像波段中心波长；设定夹角阈值（常用0.05～0.12弧度，约3°～7°）；输出每个像元的SAM角度图，角度越小异常越强。示例：以明矾石USGS光谱为参考，对OLI B2～B7做SAM，阈值取0.08 rad，得到明矾石蚀变异常图，再与羟基PCA结果叠加，可区分高级泥化带与中级泥化带。

四、三种方法的适用条件与配合策略

PCA适合多波段数据、对蚀变组合整体敏感，但主成分解释依赖载荷分析，不同影像需重新判定；比值法简单稳健、可解释性强，但单比值对矿物种类区分能力有限；SAM对特定矿物识别精度高，但强依赖参考光谱质量和大气校正精度。工程上常用“比值初筛 + PCA分级 + SAM精判”的组合：先用比值快速圈定异常带，再用PCA生成分级异常图层，最后用SAM在重点区区分矿物组合，为后续靶区排序提供依据。

五、生成铁染与羟基异常图层的关键控制点

• 波段对应关系必须核对：不同传感器波段编号不同，OLI的B6/B7对应SWIR1/SWIR2，ASTER则需用B4/B5/B6/B8/B9，套用组合前先查波段中心波长。
• 取反方向要统一：PCA与比值法都可能出现“低值=强异常”，图层符号必须统一为“高值=强蚀变”，否则后续叠加会相互抵消。
• 分级而非二值：建议输出一级、二级、三级异常，保留强度信息，便于靶区排序时区分优先级。
• 统一空间基准：所有异常图层重采样到同一栅格大小与投影，叠加前做掩膜一致性检查。

完成上述步骤后，铁染与羟基异常图层即可作为后续GIS叠加分析的输入层。需要强调的是，异常图层反映的是“光谱异常”而非“矿化确认”，其可靠性取决于大气校正质量、掩膜完整性和阈值合理性，必须与地质、构造、化探信息联合解释，才能进入靶区排序环节。