3. 矿山扰动遥感监测的波段选择与指数构建

矿山开采扰动在地表留下的痕迹，本质上是一组光谱特征的组合：植被被剥离后裸露的岩土、堆积的废石场、沉降形成的积水洼地、压占形成的硬化地面，它们在可见光—近红外—短波红外区间表现出彼此不同又相互重叠的反射行为。波段选择的任务，就是找出能把这几类地物从背景中"分出来"的光谱位置；指数构建的任务，则是把这些波段组合成对扰动敏感、对光照和地形相对稳健的数值。二者共同决定了后续分类与变化检测的信息上限。

一、波段特性：从可见光到短波红外

Sentinel-2 MSI 在矿山监测中常用的波段可归为四组。可见光（B2 蓝、B3 绿、B4 红）主要反映铁氧化物、有机质和部分矿物染色，其中 B4 与 B3 的差异对植被叶绿素吸收敏感；红边波段（B5–B7，约 705–783 nm）位于植被反射陡坡区，对植被胁迫和稀疏植被覆盖的响应比红波段更细腻，是识别"扰动边缘半退化植被"的重要补充；近红外 B8（842 nm）是植被细胞结构反射的高值区，也是水体强吸收区，因此同时承担植被增强与水体抑制两种角色；短波红外 B11（1610 nm）与 B12（2190 nm）对含水量、粘土矿物和铁铝氧化物敏感，裸岩、废石堆与压实地面在此区间反射率明显抬升，而水体几乎全吸收。

需要强调的是空间分辨率的差异：B2–B4、B8 为 10 m，B5–B7、B11、B12 为 20 m。当用 20 m 波段与 10 m 波段组合成指数时，应先将所有波段重采样到统一网格（通常重采样到 10 m，采用最近邻或双线性），否则像元错位会在扰动边界处产生虚假高值。此外，B11/B12 在 20 m 尺度上对小型采坑、窄道路的响应会被混合像元稀释，这是指数在细碎扰动上灵敏度下降的客观原因。

二、四类核心指数的构造与适用性

归一化差值形式 NDX = (A − B)/(A + B) 是这些指数的共同骨架，其优点是比值化后部分抵消了光照与地形差异，值域被压缩到 [−1, 1]，便于跨时相、跨景比较。四个指数只是把不同的波段对代入同一骨架：

NDVI = (B8  - B4 ) / (B8  + B4 )     植被增强，扰动区取低值
NDWI = (B3  - B8 ) / (B3  + B8 )     水体增强（McFeeters 形式）
BSI  = ((B11 + B4) - (B8 + B2)) / ((B11 + B4) + (B8 + B2))   裸土增强
NDBI = (B11 - B8 ) / (B11 + B8 )     建成/硬化地面增强

NDVI 是最基础的扰动指示量：开采剥离、废石压占、道路修建都会使 NDVI 从背景的 0.5–0.8 骤降到 0.1–0.3，因此"NDVI 显著下降"常作为扰动初筛条件。但 NDVI 在浓密植被区易饱和，在裸土与硬化地面之间几乎无区分力，单靠它无法判断扰动类型。

NDWI 用于识别扰动区内的积水：采坑积水、尾矿库水面、沉降洼地。需要注意 Sentinel-2 的 B8 与 B11 都可构造水体指数，B3/B8 形式对浅水更敏感，而 (B8−B11)/(B8+B11) 形式对浑浊矿坑水抑制更好。实际使用时应明确采用哪一种，并在全文保持一致，避免同一项目内混用导致阈值不可比。

BSI 通过"短波红外+红"与"近红外+蓝"的对比突出裸土：岩土裸露区 BSI 升高，植被区 BSI 为负。它对废石场、排土场、裸露边坡响应良好，但对已压实硬化的场地与裸岩区分不足。

NDBI 利用 B11 高于 B8 的特性突出硬化地面，适合识别矿区建筑、选矿厂、压实运输道路。其局限在于：部分干燥裸土的 B11 同样高于 B8，会造成 NDBI 与 BSI 同时升高的"混淆区"，这正是需要多指数联合判读的根本原因。

三、指数组合的判读逻辑

单一指数只能回答"是不是"，多指数联合才能回答"是什么"。下图给出四类典型矿山扰动在 NDVI、NDWI、BSI、NDBI 四个维度上的响应组合与判读路径，可作为阈值设定与分类规则设计的依据。

从图中可以看出，判读并非简单的阈值串联，而是"先降维、再分流"：NDVI 负责把扰动候选区从大面积植被背景中剥离出来，NDWI、BSI、NDBI 负责在候选区内区分积水、裸土与硬化地面，红边波段则用于处理边界处的半退化植被。四类响应并非互斥，实际像元常落在两类之间，此时应结合形状特征（采坑多为不规则洼地、道路呈线状、排土场呈台阶状）与邻接关系辅助判定。

四、阈值标定的实际操作

指数值本身没有绝对意义，阈值必须落到具体研究区与具体时相。一种可复现的做法是：选取扰动发生前的一期质量较好的影像（云量低于 10%），在稳定植被区、稳定裸地区各布设若干样本多边形，统计 NDVI 的均值 μ 与标准差 σ，取 μ − 2σ 作为"疑似扰动"的下限；再在已知扰动类型（如历史采坑、排土场）上统计各指数的分布区间，形成分类规则的上下界。示例性地，若某研究区稳定植被 NDVI 为 0.72 ± 0.08，则初筛阈值约为 0.56；而实际扰动区 NDVI 常降至 0.15–0.35，两者之间留有充分间隔，说明该阈值在该区是可行的。这些数值仅为说明标定流程的假设示例，实际项目必须用本地样本重新计算。

还需处理两类系统性干扰。其一是地形阴影：阴影区可见光与近红外同时被压低，NDVI 可能异常偏低而被误判为扰动，应结合坡度坡向或太阳高度角生成阴影掩膜。其二是物候差异：若对比的两期影像分属不同季节，植被本身的 NDVI 就会变化，此时应优先选择同季节影像，或改用变化向量、时间序列残差等方法而非固定阈值。指数构建只是提供了敏感的光谱变量，真正决定监测精度的，是波段组合、阈值标定与干扰掩膜三者是否与本地条件匹配。