4. 波段比值的地质意义与物理基础
波段比值是遥感蚀变信息提取中最经典、也最容易被误用的增强手段。它的本质不是“把两个波段相除”这么简单,而是利用蚀变矿物在可见光—近红外—短波红外区间形成的诊断性吸收特征,通过两个波段反射率的相对关系,把微弱的吸收差异放大成可判读的异常。要真正用好比值,必须先理解“矿物为什么在某个波长吸收”“比值为什么能把这个吸收变成高值或低值”“哪些干扰会让比值失真”。
一、物理基础:吸收特征来自电子跃迁与分子振动
在0.4~2.5 μm区间,地质遥感关注的吸收机制主要有两类。第一类是铁离子的电子跃迁:Fe³⁺在0.85~0.95 μm附近形成宽缓吸收,Fe²⁺在1.0~1.1 μm附近产生吸收,二者叠加使含铁矿物(赤铁矿、针铁矿、黄钾铁矾、绿泥石等)在可见光—近红外段呈现特征性的反射率下降。第二类是羟基和水分子团的振动吸收:OH⁻、H₂O在1.4 μm、1.9 μm附近产生强吸收,Al-OH在2.17~2.21 μm、Mg-OH在2.30~2.35 μm附近产生诊断性窄吸收,这是黏土矿物、绢云母、绿泥石、滑石等羟基蚀变矿物的“指纹”。
关键点在于:吸收是相对的。单一波段的反射率高低同时受矿物含量、粒度、表面粗糙度、光照和地形影响,绝对亮度并不可靠;但吸收谷两侧的反射率之比,能在相当程度上抵消整体亮度差异,突出“谷”的存在。这就是比值增强的物理出发点。
二、比值响应:把吸收谷变成可判读的异常
设某矿物在波长λ₁处为吸收谷、λ₂处为吸收肩(反射相对高),则比值 R = ρ(λ₂)/ρ(λ₁) 会显著大于1;若把分子分母颠倒,则得到小于1的低值异常。工程上通常约定“肩/谷”方向,使蚀变区呈高值,便于统一配色和阈值分割。以铁染为例,常用 TM3/TM1(或OLI Band4/Band2)增强Fe³⁺在蓝端的强吸收;以羟基为例,常用 TM5/TM7(OLI Band6/Band7)增强2.2 μm附近Al-OH吸收谷相对于2.1 μm肩部的差异。注意:不同传感器波段位置不同,比值组合必须按实际波段中心波长重新对应,不能照搬TM编号。
三、从比值到异常图:关键处理环节
比值本身只是中间产品,要得到可用的蚀变异常图,通常要经过“辐射定标—大气校正—比值计算—干扰掩膜—阈值分割—分级配色”这条链。其中两个环节最影响地质意义:
其一,大气校正不能省。比值虽能部分抵消乘性光照差异,但大气程辐射是加性的,会压缩吸收对比度。以OLI为例,若不做FLAASH或QUAC校正,Band2受瑞利散射抬升明显,TM3/TM1类比值会被系统性压低,铁染异常边界模糊。实操中可用暗像元法粗校,但正式制图建议用辐射传输模型校正,并检查校正后水体、深阴影的反射率是否接近0。
其二,干扰掩膜要按“比值敏感对象”设计。植被在0.67 μm有叶绿素强吸收、在0.8 μm有陡坡,会同时污染铁染和羟基比值;水体在近红外极低,任何含近红外波段的比值都会出现极端值;云影造成整体压暗但比值可能不变,需用亮度阈值单独剔除。一个可参考的示例掩膜规则(示例参数,需按影像统计调整)为:NDVI > 0.35 判为植被;Band5(0.86 μm)反射率 < 0.08 判为水体;全波段均值 < 0.05 判为深阴影;Band6 > 0.45 且 Band7/Band6 接近1 判为雪或盐碱壳。掩膜后再统计比值直方图,用均值加k倍标准差(如k=2.0~2.5)或累积频率2%~5%作为异常下限,比固定阈值更稳健。
四、比值组合的对应关系与常见误区
| 目标蚀变 | 诊断吸收 | 常用比值(示例) | 异常方向 |
|---|---|---|---|
| 铁染(Fe³⁺) | 0.85~0.95 μm 宽吸收 | Band4/Band2 | 高值 |
| 羟基(Al-OH) | 2.17~2.21 μm 窄吸收 | Band6/Band7 | 高值 |
| 镁羟基(Mg-OH) | 2.30~2.35 μm 吸收 | Band6/Band7 辅助 Band7/Band6 对比 | 需结合波段位置判断 |
| 铁氧化物+羟基复合 | 两处吸收并存 | 铁染比值与羟基比值做彩色合成 | 双高值区优先 |
最常见的三个误区是:照搬波段编号(ASTER、OLI、GF-5波段宽度不同,必须按中心波长重配);忽略比值的方向约定(同一组合分子分母颠倒,异常从高值变低值,配色和阈值全部要改);把比值高值直接当矿(比值只指示“存在该吸收特征的地表物质”,不等于矿化,必须结合地质构造、已知矿点和野外查证)。理解到这一层,比值才真正成为有地质意义的增强工具,而不是一张好看的彩色图。