4. 蚀变信息提取的波谱机理与方法选择
蚀变矿物在可见光—短波红外（VNIR-SWIR）与热红外波段具有诊断性吸收特征，这些特征由Fe³⁺、Fe²⁺、Al-OH、Mg-OH、CO₃²⁻等电子跃迁与分子振动过程产生。例如，赤铁矿在0.85～0.95 μm附近出现宽缓吸收，明矾石、高岭石在2.17 μm与2.20 μm附近形成Al-OH吸收双谷，绿泥石、绿帘石则在2.33～2.35 μm出现Mg-OH吸收。蚀变信息提取的本质，是把这些“位置—深度—宽度”可量化的吸收特征从影像光谱中稳定地分离出来。不同方法对光谱形状、波段连续性、背景干扰和先验知识的依赖程度不同，因此方法选择必须与研究区矿物组合、传感器波段设置和地质背景匹配，而不能机械套用某一种算法。
四类方法的机理与适用条件
比值法利用蚀变矿物在特征吸收波段与邻近反射峰波段的反射率比值增强微弱异常。以ASTER为例，常用B4/B2增强铁染（赤铁矿、褐铁矿），B5/B6或B7/B6增强Al-OH类蚀变（高岭石、明矾石、绢云母），B7/B5则对Mg-OH类（绿泥石、绿帘石）敏感。比值法计算简单、物理意义直观、对地形阴影和照度差异有一定抑制能力，适合大面积快速扫面和初筛。但其局限也很明显：比值只利用两个波段，难以区分吸收位置相近的矿物；当研究区存在植被、水体、云影或岩性背景干扰时，比值异常容易“假阳性”；此外，比值结果对波段配准和定标精度较敏感。
主成分分析（PCA）通过正交变换把多波段信息压缩到少数主成分上，使蚀变信息在某一主成分中与背景分离。经典做法是选取包含诊断吸收波段的波段子集（如ASTER的B1、B3、B4、B6或B5、B6、B7、B8），对子集做PCA，再根据特征向量符号判断哪一主成分代表蚀变增强方向。PCA能同时利用多个波段、抑制相关性强的背景信息，适合波段较多、蚀变类型较单一且背景相对均一的区域。其不足在于：主成分的地质含义需人工判读，特征向量符号易受场景统计分布影响，不同子区结果可能不一致；当研究区矿物组合复杂、多种蚀变并存时，单一主成分往往难以同时突出所有目标。
光谱角制图（SAM）把像元光谱与参考光谱视为高维空间向量，以两者夹角衡量相似性，夹角越小越相似。SAM对整体反射率高低不敏感，能较好适应照度与地形差异，适合有可靠端元或实测光谱、蚀变类型明确的研究区。其关键在于参考光谱的代表性与波段匹配：若参考光谱来自USGS库而传感器波段响应不同，需先做波段重采样；阈值设置过严会漏提，过松则混入非蚀变地物。SAM通常作为“已知目标识别”手段，而非无先验的异常发现手段。
匹配滤波（MF）属于部分解混方法，在已知目标光谱的前提下，抑制背景端元影响、突出目标丰度。它比SAM更强调“目标在混合像元中的相对含量”，适合蚀变矿物与围岩混合严重、端元光谱可获取的场景，如斑岩铜矿外围的绢英岩化带识别。MF对端元质量和背景统计估计较敏感，端元选取偏差会直接导致丰度图偏移；同时它输出的是相对丰度，需结合地质资料解释，不能简单当作绝对含量。
按矿物组合确定技术路线
方法选择的核心依据是“目标矿物的诊断吸收位置是否落在传感器可用波段内，以及背景干扰能否被有效压制”。以斑岩铜矿外围为例，典型蚀变分带由内向外为钾化带、绢英岩化带、青磐岩化带。绢英岩化以绢云母、高岭石为主，Al-OH吸收位于2.17～2.22 μm，ASTER B5、B6、B7可覆盖，适合先用B5/B6比值圈定绢英岩化带，再用PCA在B5、B6、B7、B8子集上分离绢云母与高岭石的相对强弱；青磐岩化以绿泥石、绿帘石为主，Mg-OH吸收位于2.33～2.35 μm，B7/B5比值更敏感，但B7与B5波段宽度较大，需结合SAM用实测绿泥石光谱确认。若研究区植被覆盖较高，比值法易受植被红边干扰，应优先采用PCA或MF压制背景，并在预处理阶段做植被掩膜。
对于高光谱数据（如HyMap、CASI/SASI），波段连续、光谱分辨率高，SAM与MF的优势更明显：SAM可直接用实测端元做逐像元夹角匹配，MF可输出绢云母、高岭石、绿泥石等端元的相对丰度图，适合精细分带。对于多光谱数据（ASTER、Landsat 8/9），波段少且宽，比值法与PCA更稳妥，SAM和MF需依赖波段重采样后的端元，精度受限于波段设置。若研究区无实测光谱，可先用比值法或PCA发现异常，再结合USGS光谱库选取端元做SAM验证，避免“无先验硬套SAM”。
实操中的关键控制点
无论采用哪种方法，预处理一致性是前提：辐射定标、大气校正（如FLAASH或QUAC）、坏波段剔除、云与水体掩膜必须统一，否则比值与PCA结果会在不同景之间漂移。以ASTER为例，B4饱和、B9受水汽影响，通常不参与蚀变提取；B1与B2用于铁染，B5～B7用于Al-OH与Mg-OH。PCA子集选择要围绕目标吸收位置，而非随意取波段；特征向量判读时，应检查目标波段载荷符号是否与吸收方向一致，必要时对主成分取反。SAM阈值建议先在已知蚀变露头区做试算，取夹角分布的低分位作为初值，再结合野外验证调整；MF则需保证背景统计样本纯净，避免把蚀变像元混入背景协方差估计。
综上，方法选择不是“哪个先进用哪个”，而是由矿物组合、波段覆盖、背景干扰和先验光谱共同决定。比值法适合快速定位，PCA适合压制背景，SAM适合目标确认，MF适合丰度量化；将四者按“扫面—净化—确认—量化”的顺序组合，才能形成与研究区地质条件匹配、可验证、可复现的蚀变提取技术路线。