NASA大气参数计算器停服后的应急方案：USGS Landsat地表温度数据全流程解析
当NASA大气参数计算器在2024年1月停止服务后，许多依赖该工具进行地表温度(LST)研究的用户陷入了困境。本文将详细介绍如何利用USGS Landsat Collection 2地表温度成品数据作为替代方案，从数据获取到ENVI处理的完整流程。
1. USGS Landsat Collection 2地表温度数据概述
USGS Landsat Collection 2地表温度产品是NASA大气参数计算器停服后的理想替代方案。该数据集已经完成了复杂的大气校正和温度反演过程，用户可以直接使用成品数据，无需自行计算大气参数。
主要优势：
- 省去了复杂的大气参数计算过程
- 数据质量经过USGS严格验证
- 覆盖全球范围，时间序列完整
- 提供多种辅助数据产品
USGS Landsat Collection 2地表温度产品包含多个文件，其中最关键的是*_ST_B10.tif文件，该文件存储了反演完成的地表温度数据，单位为开尔文(K)。
2. 数据获取与下载指南
2.1 访问USGS EarthExplorer平台
- 打开浏览器，访问 USGS EarthExplorer
- 注册或登录USGS账号
- 在搜索条件中设置感兴趣的区域和时间范围
2.2 筛选Landsat Collection 2数据
在数据集选项中，选择：
Landsat > Landsat Collection 2 Level-2 > Landsat 8-9 OLI/TIRS C2 L2重要提示：确保选择的是"Level-2"产品，而不是"Level-1"数据。
2.3 下载地表温度产品
在搜索结果页面，点击感兴趣的场景，然后选择下载以下文件：
*_ST_B10.TIF（地表温度数据）*_MTL.txt（元数据文件）
3. ENVI数据处理流程
3.1 数据导入与初步检查
- 打开ENVI软件
- 使用
File > Open导入下载的*_ST_B10.TIF文件 - 右键点击图层，选择
Quick Stats检查数据范围
常见问题：
- 数据值异常大（这是正常现象，因为USGS使用了缩放存储）
- 某些区域显示为NoData（可能是云覆盖或数据质量问题）
3.2 温度值转换与校正
USGS存储的温度数据使用了缩放系数和偏移量，需要进行转换才能得到真实的开尔文温度值。根据USGS文档，转换公式为：
真实温度(K) = DN值 × 0.00341802 + 149在ENVI中执行此转换：
- 打开
Band Math工具 - 输入公式：
b1*0.00341802+149 - 指定输出文件名和位置
3.3 开尔文转摄氏度
将温度从开尔文转换为摄氏度：
- 再次打开
Band Math - 输入公式：
b1 - 273.15 - 使用上一步的结果作为输入
示例代码：
; 第一步：缩放校正 corrected_k = b1 * 0.00341802 + 149 ; 第二步：单位转换 final_c = corrected_k - 273.154. 数据质量评估与应用建议
4.1 质量检查方法
USGS提供了*_ST_QA.TIF文件用于质量评估：
- 值为0：高质量数据
- 值为1：中等质量
- 值为2：低质量/可能有问题
在ENVI中可以通过以下步骤进行质量筛选：
- 导入QA文件
- 使用
Band Math创建掩膜：b1 eq 0（仅保留高质量数据） - 将掩膜应用到温度数据
4.2 与历史数据对比
根据用户反馈，USGS Collection 2地表温度产品与使用NASA计算器得到的结果具有良好的一致性：
- 平均差异：0.5-1.5°C
- 空间分布模式高度相似
- 极端值出现频率相近
4.3 应用注意事项
- 空值处理：某些区域可能由于云覆盖或地表发射率数据缺失而出现空值
- 时间一致性：不同季节的数据可能因大气条件差异而存在系统性偏差
- 空间分辨率：热红外波段的空间分辨率约为100米，不适合小尺度精细分析
5. 高级技巧与替代方案
5.1 处理缺失数据问题
当研究区域存在数据缺失时，可考虑以下解决方案：
- 使用相邻像元插值：
; 使用3x3窗口均值填充空值 filled = focal_stat(b1, /mean, size=[3,3], /ignore_nan)使用同一区域其他日期的数据替代
手动计算地表发射率（适用于高级用户）
5.2 批量处理多个场景
对于需要处理大量数据的用户，可以编写ENVI批处理脚本：
pro process_landsat_lst, filelist foreach file, filelist do begin ; 数据导入 envi_open_file, file, r_fid=fid ; 温度转换 envi_doit, 'math_doit', fid=fid, $ exp='(b1*0.00341802+149)-273.15', $ out_name=file+'_LST_C' endforeach end5.3 与其他数据源对比
将Landsat LST数据与其他来源（如MODIS）进行对比时，应注意：
- 空间分辨率的差异
- 过境时间的差异
- 算法原理的不同
典型对比结果：
| 数据源 | 平均差异 | 主要差异原因 |
|---|---|---|
| MODIS | +2-3°C | 分辨率差异 |
| ASTER | ±1°C | 算法不同 |
| 地面测量 | ±2°C | 尺度不匹配 |
6. 常见问题解答
Q：USGS Landsat Collection 2地表温度产品的精度如何？
A：根据USGS官方验证，产品绝对精度约为1-2°C，相对精度（同一区域不同时间）可达0.5°C以内。
Q：为什么我的研究区域有大片缺失数据？
A：这通常是由于ASTER GED地表发射率数据缺失造成的。可以尝试使用相邻像元插值或选择其他日期的数据。
Q：能否直接使用这些数据进行气候变化研究？
A：可以，但建议进行严格的质量控制和一致性检查，特别是对于长时间序列分析。
Q：除了ENVI，还能用什么软件处理这些数据？
A：QGIS、ArcGIS、Python（GDAL/rasterio）等常见地理空间软件都能处理这些数据。
