4. 波段组合分析原理
ASTER 在可见光—短波红外—热红外区间共设 14 个波段，其中 SWIR 的 6 个波段（B4–B9）是识别绢云母、高岭石、明矾石、绿泥石、方解石等热液蚀变矿物的核心窗口。但原始单波段影像上，蚀变矿物微弱的光谱吸收特征往往被地形阴影、反照率差异和大气残余所掩盖。波段组合分析的本质，就是通过代数运算或统计变换，把“某一吸收特征的有无与强弱”从多波段亮度中分离出来，形成一幅对比度高、可判读的蚀变异常图。本节围绕四类方法展开其数学逻辑与增强机制。
一、波段比值：用除法消去乘性干扰
波段比值的基本形式为 R = Bi / Bj。其增强机制在于：地形坡度、照度差异、颗粒度变化在相邻波段上近似为乘性因子 k，即观测值 Li ≈ k·ρi，相除后 k 被约去，只留下反射率的相对形状。因此比值图能显著压制阴影、突出吸收特征。
对蚀变识别而言，关键是让分子落在吸收谷、分母落在吸收肩部或相邻无吸收波段。例如绢云母在 2.20 μm（B6）有强 Al–OH 吸收，而 B5（2.16 μm）与 B7（2.26 μm）位于肩部，故 B5/B6 或 B6/B7 的比值越高，指示 2.20 μm 吸收越深。明矾石、高岭石的吸收位置略有偏移，需相应调整波段对。比值法的局限是：当分母接近零或噪声较大时结果不稳定，且比值丢失了绝对反射率信息，难以区分“真吸收”与“整体暗化”。
二、波段运算：以加减构造吸收深度
更稳健的做法是用“肩部均值减谷值”构造吸收深度指数，例如 (B5 + B7)/2 − B6，或归一化形式 (B5 + B7 − 2B6)/(B5 + B7 + 2B6)。前者是绝对吸收深度，后者是相对吸收深度，可抑制反照率差异。这类运算的物理含义直接对应光谱学中的吸收深度（band depth），比单纯比值更接近矿物丰度的单调关系。实际应用中常配合阈值分割或假彩色合成（如 R= B6/B8、G= B5/B6、B= B4/B2）快速圈定蚀变带。
三、主成分分析：让蚀变信息落到某一分量
PCA 对多波段协方差矩阵做正交分解，把相关性强的波段能量集中到前几个主成分。其增强机制不是“消去干扰”，而是“重新分配方差”：若某矿物吸收使 B6 相对 B5、B7 变暗，则这种“B6 与其余波段反号”的结构会集中出现在某个主成分的载荷向量上。经典做法是选取 SWIR 四波段（B5–B8）做 PCA，检查各分量载荷：当某分量在 B6 上为负、在 B5/B7 上为正时，该分量即“含 Al–OH 蚀变信息的主成分”，常取 PC3 或 PC4。对含铁矿物则用 VNIR 的 B1–B4 做 PCA，选取 B2 与 B4 反号的分量。PCA 的优点是无需预设波段对，能同时利用多波段协同信息；缺点是结果依赖输入波段子集与场景统计，不同影像的主成分编号不可直接比较，需逐景检查载荷。
四、光谱指数构建：把物理量归一化
光谱指数是把吸收深度写成归一化形式，使其对不同光照、不同传感器增益具有可比性。以绢云母指数为例：
SI_sericite = (B5 + B7 - 2*B6) / (B5 + B7 + 2*B6)
# 值越大，2.20 μm 吸收越深，绢云母/伊利石类越富集
# 示例阈值（需按场景统计调整）：SI > 0.03 判为异常候选类似地可构造碳酸盐指数（利用 B8 附近 2.35 μm 吸收）、绿泥石指数（利用 B7/B8 与 Fe 吸收）。指数法的关键是：分子分母使用同一组波段，保证量纲一致；对负值或近零分母做掩膜处理；阈值必须基于本景直方图或已知矿点标定，而非套用固定数值。
五、方法选择与验证
四类方法并非互斥：实践中常先用比值或指数快速扫描，再用 PCA 交叉验证异常的空间一致性，最后用已知矿点或野外查证点检验命中率。判断增强是否成功，可观察三点：异常图是否与地层、断裂、岩体接触带等控矿要素空间吻合；异常是否呈条带状或环带状而非随机噪点；同一蚀变类型在不同景的指数值分布是否可比。若比值图噪声大，改用吸收深度运算；若 PCA 分量难以解释，回到载荷符号逐波段核对。只有把数学运算与矿物光谱机理对应起来，波段组合才真正成为蚀变信息增强的工具，而不是盲目的图像变换。