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.00341802149指定输出文件名和位置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_fidfid ; 温度转换 envi_doit, math_doit, fidfid, $ exp(b1*0.00341802149)-273.15, $ out_namefile_LST_C endforeach end5.3 与其他数据源对比将Landsat LST数据与其他来源如MODIS进行对比时应注意空间分辨率的差异过境时间的差异算法原理的不同典型对比结果数据源平均差异主要差异原因MODIS2-3°C分辨率差异ASTER±1°C算法不同地面测量±2°C尺度不匹配6. 常见问题解答QUSGS Landsat Collection 2地表温度产品的精度如何A根据USGS官方验证产品绝对精度约为1-2°C相对精度同一区域不同时间可达0.5°C以内。Q为什么我的研究区域有大片缺失数据A这通常是由于ASTER GED地表发射率数据缺失造成的。可以尝试使用相邻像元插值或选择其他日期的数据。Q能否直接使用这些数据进行气候变化研究A可以但建议进行严格的质量控制和一致性检查特别是对于长时间序列分析。Q除了ENVI还能用什么软件处理这些数据AQGIS、ArcGIS、PythonGDAL/rasterio等常见地理空间软件都能处理这些数据。
NASA大气参数计算器停服后，手把手教你用USGS Landsat Collection 2地表温度成品数据（附ENVI实操）
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.00341802149指定输出文件名和位置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_fidfid ; 温度转换 envi_doit, math_doit, fidfid, $ exp(b1*0.00341802149)-273.15, $ out_namefile_LST_C endforeach end5.3 与其他数据源对比将Landsat LST数据与其他来源如MODIS进行对比时应注意空间分辨率的差异过境时间的差异算法原理的不同典型对比结果数据源平均差异主要差异原因MODIS2-3°C分辨率差异ASTER±1°C算法不同地面测量±2°C尺度不匹配6. 常见问题解答QUSGS Landsat Collection 2地表温度产品的精度如何A根据USGS官方验证产品绝对精度约为1-2°C相对精度同一区域不同时间可达0.5°C以内。Q为什么我的研究区域有大片缺失数据A这通常是由于ASTER GED地表发射率数据缺失造成的。可以尝试使用相邻像元插值或选择其他日期的数据。Q能否直接使用这些数据进行气候变化研究A可以但建议进行严格的质量控制和一致性检查特别是对于长时间序列分析。Q除了ENVI还能用什么软件处理这些数据AQGIS、ArcGIS、PythonGDAL/rasterio等常见地理空间软件都能处理这些数据。
