Landsat C2L2数据预处理与遥感分析指南

1. Landsat C2L2数据基础解析

Landsat Collection 2 Level 2（简称C2L2）是美国地质调查局（USGS）提供的经过大气校正的二级产品数据集，包含地表反射率（Surface Reflectance）和地表温度（Land Surface Temperature）两大核心衍生数据产品。与原始Level 1数据相比，C2L2产品已经过辐射定标、大气校正和几何校正处理，可直接用于定量遥感分析。

这个数据集主要来自Landsat 8和Landsat 9卫星的OLI（ Operational Land Imager）和TIRS（Thermal Infrared Sensor）传感器。其中OLI提供可见光到短波红外的多光谱数据（波段1-7），TIRS提供热红外数据（波段10-11）。每个标准场景覆盖185×180公里的区域，空间分辨率多为30米（热红外波段为100米，但已重采样到30米）。

关键提示：Collection 2相比前代Collection 1进行了多项改进，包括更好的几何校正、更精确的地形校正以及改进的大气补偿算法。如果是长期变化监测研究，建议不要混用不同Collection版本的数据。

2. 数据获取与初步检查

2.1 官方数据下载途径

USGS提供多种数据获取渠道，最常用的是EarthExplorer平台（https://earthexplorer.usgs.gov/）。下载前需要注册免费账号，搜索时需注意以下关键筛选条件：

• 数据集选择"Landsat Collection 2 Level-2"
• 产品类型可选"L2SP"（含反射率和温度）或"L2SR"（仅含反射率）
• 根据需求选择Tier级别："T1"（最高质量）或"T2"（次级质量）

下载的数据包通常包含多个GeoTIFF文件，每个文件对应一个波段或辅助数据。命名规则示例：
LC08_L2SP_123045_20220520_20220527_02_T1_SR_B4.tif
其中各部分含义：

• LC08：Landsat 8
• L2SP：处理级别
• 123045：路径/行号
• 20220520：获取日期
• SR_B4：波段4地表反射率

2.2 数据质量评估

在预处理前必须检查数据质量，主要关注元数据中的以下指标：

1. IMAGE_QUALITY_OLI：OLI影像质量（1-9，越高越好）
2. CLOUD_COVER：云量百分比
3. GEOMETRIC_RMSE_MODEL：几何校正误差（应<30米）
4. SUN_ELEVATION：太阳高度角（影响辐射计算）

特别要注意QA波段（QA_PIXEL），它采用位掩码标记各种质量状况。例如，要提取云掩膜可使用以下位运算：

python复制# 提取云掩膜（位3为1表示云）
cloud_mask = (qa_band & (1 << 3)) != 0  

3. 像元值到物理量的转换

3.1 地表反射率计算

Landsat C