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编号。

诊断性吸收特征(示意) 波长 λ (μm) 反射率 ρ 0.5 0.9 1.4 1.9 2.2 2.35 Fe³⁺ 吸收 0.85–0.95 μm Al-OH 吸收 2.17–2.21 μm H₂O 1.9 比值响应与异常方向 铁染异常:Band4 / Band2(TM3/TM1) 分子取蓝端吸收肩,分母取 Fe³⁺ 吸收谷 蚀变区比值 ↑ 高值 → 暖色(红/橙) 典型对象:赤铁矿化、褐铁矿化、黄钾铁矾 羟基异常:Band6 / Band7(TM5/TM7) 分子取 2.1 μm 肩部,分母取 Al-OH 吸收谷 蚀变区比值 ↑ 高值 → 青/蓝或按方案配色 典型对象:绢云母化、高岭石化、绿泥石化 判读前提 比值高 ≠ 一定有蚀变:需排除植被、水体、 云影、雪、盐碱壳与地形阴影的伪异常。

三、从比值到异常图:关键处理环节

比值本身只是中间产品,要得到可用的蚀变异常图,通常要经过“辐射定标—大气校正—比值计算—干扰掩膜—阈值分割—分级配色”这条链。其中两个环节最影响地质意义:

其一,大气校正不能省。比值虽能部分抵消乘性光照差异,但大气程辐射是加性的,会压缩吸收对比度。以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波段宽度不同,必须按中心波长重配);忽略比值的方向约定(同一组合分子分母颠倒,异常从高值变低值,配色和阈值全部要改);把比值高值直接当矿(比值只指示“存在该吸收特征的地表物质”,不等于矿化,必须结合地质构造、已知矿点和野外查证)。理解到这一层,比值才真正成为有地质意义的增强工具,而不是一张好看的彩色图。