3. Landsat 9数据产品解读
在USGS EarthExplorer或Earthdata上下载Landsat 9影像时，最先遇到的岔路口就是产品级别选择：Collection 2的Level-1（L1）还是Level-2（L2）。很多初学者凭直觉选“更高级”的L2，但在矿区场景中，这个选择直接决定了后续能否做辐射一致性分析、能否跨时相拼接、以及分类阈值是否可迁移。本节把两类产品的物理含义、文件构成、QA波段编码和元数据字段逐一拆开，最后给出一张可执行的矿区选型判断表。
一、L1与L2的本质差异：DN值与地表反射率
L1产品提供的是经过几何精校正和辐射定标、但未做大气校正的量化DN值（Digital Number），存储为16位无符号整型。它保留了传感器端的原始辐射信息，包括大气散射、气溶胶影响和地形阴影。L1的价值在于“可回溯”：任何大气校正算法都可以从L1重新出发，且不同算法之间的差异可以被显式控制。对于矿区这种地表反射率动态范围极大（裸露矸石、煤堆、水体、植被混杂）的场景，L1是自建大气校正流程、做辐射一致性对比实验的起点。
L2产品则是在L1基础上，经LaSRC（Landsat 8/9）算法完成大气校正后生成的地表反射率（Surface Reflectance, SR）产品，像元值为0～10000的整型缩放反射率，真实反射率 = DN / 10000。同时L2还附带地表温度（ST）产品（仅热红外波段，缩放因子0.00341802，偏移量149.0，单位开尔文）。L2省去了大气校正环节，适合快速出图、大范围时序监测和跨传感器比较（Landsat 8/9 L2 SR具有一致的算法基线）。
二、QA波段结构：像素级的质量说明书
L1和L2都附带QA_PIXEL（质量评估像素波段）和QA_RADSAT（辐射饱和波段），L2额外提供QA_AEROSOL（气溶胶质量波段）。QA_PIXEL是一个16位整型波段，每一位（bit）代表一种质量标志，必须用位运算解码，不能当作普通数值读取。以Collection 2的QA_PIXEL为例，常用位定义如下：
| 位（bit） | 含义 | 取值说明 |
|---|---|---|
| 0 | 填充（Fill） | 1=该像元无数据 |
| 1 | 膨胀云（Dilated Cloud） | 1=云区膨胀范围 |
| 2 | 卷云（Cirrus） | 1=检测到卷云 |
| 3 | 云（Cloud） | 1=云 |
| 4 | 云阴影（Cloud Shadow） | 1=云阴影 |
| 5 | 雪/冰（Snow） | 1=雪或冰 |
| 6 | 晴空（Clear） | 1=晴空（无云） |
| 7 | 水体（Water） | 1=水体 |
| 8～9 | 云置信度 | 00=未检测，01=低，10=中，11=高 |
| 10～11 | 云阴影置信度 | 同上编码 |
| 12～13 | 雪/冰置信度 | 同上编码 |
| 14～15 | 卷云置信度 | 同上编码 |
实际解码时，例如要提取“晴空且非云”的像元，可用 (QA & (1<<6)) != 0 判断晴空位；要排除云和云阴影，则用 (QA & ((1<<3)|(1<<4))) == 0。在Python中配合rasterio读取后，用numpy位运算即可生成掩膜。矿区影像中，矸石山和煤堆的亮温特征有时会被云检测算法误判为云，因此不能只看bit 3，还要结合bit 8～9的云置信度，把“高置信度云”作为硬掩膜，“中置信度”作为人工复核区。
三、元数据关键字段：别只看文件名
Landsat 9的元数据以MTL（Metadata）文本文件形式随产品提供，L1为 *_MTL.txt，L2为 *_MTL.xml（也可转txt）。矿区预处理最需要关注的字段集中在辐射定标、太阳几何和产品标识三类。
辐射定标字段：L1的 RADIANCE_MULT_BAND_x 与 RADIANCE_ADD_BAND_x 用于DN转辐亮度，REFLECTANCE_MULT_BAND_x 与 REFLECTANCE_ADD_BAND_x 用于DN转大气顶层反射率（TOA）。以Band 4为例，若 REFLECTANCE_MULT_BAND_4=2.0000E-05、REFLECTANCE_ADD_BAND_4=-0.100000，则 TOA反射率 = DN × 2.0e-5 − 0.1，再除以 SUN_ELEVATION 的正弦值完成太阳高度角校正。L2的SR波段则直接给出 REFLECTANCE_MULT_BAND_x=0.0000275、REFLECTANCE_ADD_BAND_x=-0.2，即 SR = DN × 0.0000275 − 0.2，无需再做太阳角校正（已内置）。
太阳几何字段：SUN_ELEVATION（太阳高度角）和 SUN_AZIMUTH（太阳方位角）是地形校正和阴影分析的基础。矿区地形起伏大，若做地形校正（如C校正、SCS校正），必须从MTL读取这两个值，而不是凭经验估计。同一景影像内太阳角是常数，但跨季节时序分析时，太阳高度角差异会显著影响阴影长度和反射率，需要在预处理阶段记录并纳入一致性评估。
产品标识字段：LANDSAT_PRODUCT_ID 编码了传感器、路径行号、采集日期和处理级别，例如 LC09_L2SP_123032_20230615_20230621_02_T1 中，LC09 表示Landsat 9 OLI/TIRS，L2SP 表示Level-2地表反射率产品，123032 为Path/Row，02 为Collection 2，T1 为Tier 1（几何精度最高）。矿区若做变化检测，应优先选用Tier 1产品，避免Tier 2的几何偏差引入伪变化。
四、矿区场景下的产品决策与实操建议
把上述知识落到矿区实际任务，可以按“任务目标—产品选择—关键字段—QA策略”四步走。以“矿区植被覆盖度年际变化监测”为例：目标是提取每年夏季植被指数（NDVI）并比较，此时应选L2 SR产品，因为NDVI对大气校正敏感，L2的SR已消除大气影响，且跨年基线一致；读取 REFLECTANCE_MULT/ADD 完成缩放，用QA_PIXEL的bit 3、bit 4、bit 6生成云和云影掩膜，只保留晴空像元；同时记录每景的 SUN_ELEVATION，若两年间差异超过10°，需在分析中说明其对NDVI的潜在影响。
再以“矸石山自燃热异常识别”为例：目标是检测地表温度异常，应选L2 ST产品（直接可用）或L1热红外波段自反演。若选L2 ST，注意其缩放因子0.00341802和偏移149.0，反演后单位为开尔文；同时用QA_PIXEL排除云像元，避免云顶低温掩盖真实热异常。若矿区存在浓烟，L2的QA_AEROSOL可能标记为高气溶胶，此时SR质量下降，热红外相对受影响较小，但仍需在报告中注明。
最后给出一份可直接套用的矿区产品选型检查表，供下载后逐项勾选：
| 检查项 | 判断依据 | 推荐选择 |
|---|---|---|
| 是否需要自建大气校正 | 是 | L1 |
| 是否跨L8/L9长时序 | 是 | L2 SR |
| 是否存在煤尘/自燃烟雾 | 是 | L1优先，L2核查QA_AEROSOL |
| 是否做热异常反演 | 是 | L2 ST 或 L1 热红外 |
| 几何精度要求 | 变化检测 | Tier 1 |
| 云掩膜方式 | 所有任务 | QA_PIXEL位运算 |
掌握L1与L2的物理差异、QA位域解码和MTL关键字段后，后续的辐射定标、大气校正和云掩膜步骤才能建立在正确的数据起点上。矿区遥感预处理的第一原则是“先定产品、再谈算法”，产品选错，后面所有处理都可能在错误的基线上运行。