2. 遥感岩性识别物理基础
遥感岩性识别的本质,是把传感器记录的地表反射或发射能量,还原为矿物组成与岩石类型的信息。要做到这一点,必须先理解一条主线:太阳辐射(或地表热辐射)进入地表后,被矿物中的电子跃迁、分子振动等微观过程选择性吸收,形成具有诊断性的光谱特征;传感器在不同波段窗口记录这些特征,我们再用波段位置、吸收深度、吸收宽度等参数反演矿物。本节围绕可见光—近红外(VNIR)、短波红外(SWIR)、热红外(TIR)三个关键谱段,讲清“哪个波段看什么矿物、光谱特征由什么物理机制产生、实际识别时怎么用”。
一、光谱特征形成的物理机制
矿物光谱的选择性吸收主要来自两类过程。第一类是电子过程,发生在可见光—近红外区。过渡金属离子(Fe³⁺、Fe²⁺、Mn²⁺等)的 d 轨道电子在晶体场中被分裂,吸收特定能量光子发生跃迁,产生宽而强的吸收带。例如赤铁矿(Fe₂O₃)在 0.53 μm 附近有强吸收,使岩石呈红色;含 Fe²⁺ 的辉石、橄榄石在 1.0 μm 附近出现宽吸收。第二类是分子振动过程,主要发生在短波红外区。OH、H₂O、CO₃²⁻ 等基团的伸缩与弯曲振动产生倍频和合频吸收,位置稳定、形态尖锐,是识别含羟基硅酸盐、碳酸盐、硫酸盐的关键。热红外区则以硅氧四面体骨架振动(Reststrahlen 效应)为主,吸收位置随 SiO₂ 含量升高向短波方向移动,是区分酸性、中性、基性岩的重要依据。
需要强调的是,矿物光谱不是孤立的一条曲线,而是诊断吸收位置 + 吸收深度 + 吸收宽度 + 整体反射率高低的组合。同一吸收位置在不同矿物中可能重叠,因此实际识别往往依赖多个特征波段的组合判断,而不是单一波段。
二、三个谱段各自“看什么”
VNIR(0.4–1.0 μm)主要看铁质矿物与岩石色调。这一区间由电子跃迁主导,吸收带通常宽而缓。赤铁矿、针铁矿、褐铁矿等铁氧化物在 0.48–0.55 μm 和 0.85–0.95 μm 附近形成特征吸收,是识别铁染蚀变、氧化带和红层的关键。含 Fe²⁺ 的镁铁质矿物(辉石、角闪石、橄榄石)在 1.0 μm 附近有宽吸收,可用于辅助判别基性—超基性岩。实际应用中,常用波段比值(如红/蓝比值)增强铁染信息,但要注意该谱段受大气散射和植被影响较大。
SWIR(1.0–2.5 μm)是岩性识别信息量最大的窗口。OH、H₂O、CO₃²⁻ 的合频与倍频吸收集中于此,位置稳定、形态尖锐,是识别黏土矿物、碳酸盐、硫酸盐、云母类矿物的“主战场”。典型诊断位置包括:高岭石 2.17/2.21 μm 双吸收、绢云母(白云母类)2.20 μm 尖锐吸收、明矾石 2.17 μm 强吸收、方解石与白云石 2.34/2.32 μm 宽吸收、石膏 1.75 μm 与 2.22 μm 吸收。由于这些特征与热液蚀变分带高度对应,SWIR 成为矿产勘查中蚀变矿物填图的核心谱段。
TIR(8–14 μm)主要看硅酸盐骨架与岩类划分。硅氧四面体骨架振动使岩石在热红外区产生 Reststrahlen 吸收,吸收中心随 SiO₂ 含量增加向短波方向移动:酸性岩(花岗岩、流纹岩)约在 8–9 μm,中性岩(安山岩)约 9–10 μm,基性岩(玄武岩、辉长岩)约 10–11 μm。因此 TIR 特别适合在植被覆盖较少、岩石裸露区做岩类快速划分和 SiO₂ 含量半定量反演,且受光照条件影响小,可昼夜成像。
三、实际识别中的关键操作与判断
以一条实测(或示例)光谱曲线为例,识别流程可以落到具体参数上。首先做连续统去除:在吸收带两侧选取包络点,用直线或多项式拟合连续统,将原反射率除以连续统值,得到归一化光谱。这样做的目的是消除整体反射率高低和地形阴影的影响,突出吸收形态。随后计算吸收深度:以吸收中心波长处的连续统值为基准,吸收深度 D = 1 − R_center / R_continuum,D 越大通常表示该矿物丰度越高或吸收越强。再计算吸收宽度(半高宽)和吸收面积,用于辅助判断结晶度与混合程度。
判断时要注意几个易错点。其一,吸收位置重叠:2.20 μm 附近绢云母、伊利石、蒙脱石都可能出现吸收,需结合 2.35 μm 附近碳酸盐特征和 1.4/1.9 μm 水汽吸收综合判断。其二,混合像元效应:一个像元内多种矿物叠加会使吸收变浅、位置偏移,需用光谱解混或端元提取处理。其三,大气与植被干扰:1.4 μm 和 1.9 μm 附近大气水汽强吸收,应避开或做大气校正;植被覆盖会引入叶绿素与水分吸收,削弱矿物特征。其四,示例参数需以实际数据为准:上述吸收位置为矿物学典型值,实际影像受分辨率、信噪比和定标精度影响,需用已知样本点验证后再推广。
把上述物理机制与波段对应关系掌握牢固,后续选择传感器、设计波段组合、解释分类结果时就有了判断依据:知道为什么 2.2 μm 对绢云母敏感、为什么 8–14 μm 能区分酸性岩与基性岩,才能在训练样本制作和岩性填图中做出合理的光谱特征选择,而不是机械套用波段比值。