1. 矿区遥感预处理任务全景

在矿区开展遥感监测，第一步不是急着做分类或提取，而是把原始影像处理成“可信、可比、可量测”的基础数据。Landsat 9 于 2021 年 9 月发射，搭载 OLI-2 与 TIRS-2 两台传感器，延续了 Landsat 8 的轨道与波段设计，这为矿区长期时序分析提供了难得的连续性。理解它的数据特性，是设计预处理流程的起点。

Landsat 9 的数据特性

Landsat 9 运行在约 705 km 的太阳同步轨道，重访周期 16 天，幅宽约 185 km，空间分辨率方面：全色波段 15 m，多光谱波段 30 m，热红外波段 100 m（重采样后常以 30 m 产品分发）。OLI-2 提供 9 个波段，包括海岸/气溶胶波段（B1）、可见光与近红外（B2–B5）、短波红外（B6、B7）、全色（B8）和卷云波段（B9）；TIRS-2 提供 B10、B11 两个热红外波段。数据以 16 位量化，动态范围远大于 8 位产品，这意味着辐射信息丰富，但也要求我们在拉伸与显示时格外小心。

更关键的是产品级别。Landsat Collection 2 提供 Level-1（L1TP，几何校正但未做大气校正）和 Level-2（L2SP，表面反射率与表面温度）两类产品。对矿区而言，L2 的表面反射率产品能显著减少大气影响，但并非所有历史时相都齐全，且云掩膜质量参差。因此预处理流程必须能同时兼容 L1 与 L2 输入，并在元数据中记录来源级别。

矿区场景带来的特殊挑战

矿区不是均匀的自然地表。露天采场、排土场、尾矿库、沉陷区与工业广场在空间上高度异质，光谱上表现为裸岩、煤尘、水体、植被与人工建筑的混合像元。这带来三类典型问题：其一，地形起伏与台阶状采场造成强烈的阴影与坡向效应，同一岩性在阳坡与阴坡反射率差异明显；其二，煤尘与铁氧化物改变地表光谱，使常规植被指数或水体指数阈值失效；其三，采动沉陷导致地表湿度与水体分布快速变化，时序上出现“同物异谱、异物同谱”。此外，矿区常伴随扬尘与局地雾霾，气溶胶空间分布不均，单一的大气校正参数难以覆盖全幅。

预处理流程的整体结构

面向矿区的预处理可以组织为一条“输入—校正—增强—成果”的链路，但各环节并非简单串行，而是存在反馈与分支。下图给出本节建议的任务全景，重点表达数据级别分流、几何与辐射校正的先后关系，以及成果输出前的质量门控。

图中把“级别判定”放在几何校正之前，是因为 L2 产品已包含大气校正与地形校正，若再对 L2 做一次大气校正会造成重复处理；而 L1 产品必须先完成大气校正才能进入反射率分析。几何精校正与辐射校正的先后顺序在实践中可以互换，但建议先几何后辐射：几何重采样会轻微改变像元光谱统计，若先做辐射再重采样，统计参数与最终像元不完全对应。

成果要求与判断依据

矿区预处理成果至少要满足三条：可比、可量测、可追溯。可比指同一矿区不同时相影像在同一投影、同一分辨率、同一辐射基准下可直接比较；可量测指反射率与温度值具有物理意义，能支撑指数计算与阈值判定；可追溯指每一步参数、软件版本、输入级别都记录在案。具体到质量门控，建议设定：几何配准残差不超过 0.5 个像元（以 30 m 多光谱计约 15 m），残云率低于 5%，反射率有效值落在 0–1 区间且异常值有标记，热红外温度与当地气象站观测的偏差在合理范围内。

以一个露天煤矿为例：选取 2023 年 6 月与 2024 年 6 月两期 L2SP 产品，先检查 QA_PIXEL 云掩膜，发现 2024 年期西部排土场有薄云残留，人工勾画补充掩膜；随后用同一 DEM 做正射复核，配准残差 0.3 像元；再计算归一化植被指数与归一化水体指数，发现采场内部因煤尘导致指数偏低，于是在成果中附加“煤尘影响区”标记图层，供后续分析时区分。这个例子说明，预处理不是机械套流程，而是围绕矿区实际地表条件做判断与补充。