保姆级教程：用ENVI 5.6和Landsat8数据，一步步反演城市地表温度（附单窗算法公式详解）

从零开始ENVI 5.6与Landsat8数据的地表温度反演实战指南当城市热岛效应日益成为环境研究的焦点地表温度LST反演技术正从专业领域走向大众视野。对于环境监测、城市规划甚至农业研究而言掌握这项技能意味着能够直接看见热量在地表的分布规律。本文将带您用ENVI 5.6和Landsat8数据像搭积木一样逐步构建完整的地表温度反演流程——无需任何遥感基础只需跟随操作您就能将卫星数据转化为直观的温度分布图。1. 环境准备与数据获取工欲善其事必先利其器。在开始温度反演前我们需要准备好软件环境和数据资源。ENVI 5.6作为遥感处理的主流工具其界面友好性和算法可靠性已得到业界广泛验证。建议安装时勾选所有扩展模块特别是Thermal Analysis工具包这对后续热红外波段处理至关重要。数据获取方面Landsat8的Level 1级数据是最佳选择它包含11个波段其中波段1010.6-11.19μm和波段1111.5-12.51μm是专门用于地表温度反演的热红外波段可见光与近红外波段1-7波段则用于辅助计算植被指数推荐通过USGS EarthExplorer平台下载数据其优势在于提供完整元数据包括成像时间、太阳高度角等关键参数数据经过系统辐射校正支持自定义区域裁剪减少下载体积注意下载时务必选择Landsat 8-9 OLI/TIRS C2 L1产品级别并确认数据质量标签为9最高质量2. 数据预处理从原始数据到辐射亮度2.1 辐射定标将DN值转化为物理量打开ENVI后通过File Open加载下载的Landsat8数据通常为.tar.gz压缩包。右键点击热红外波段B10或B11选择Radiometric Calibration进入定标界面。关键参数设置如下参数项推荐值物理意义Calibration TypeRadiance转换为辐射亮度值Output InterleaveBIL优化大文件处理效率Scale Factor0.1保持数值精度Output Data TypeFloat避免计算精度损失定标后生成的新图层其数值范围通常在5-11 W/(m²·sr·μm)之间。此时可通过Display Color Mapping进行伪彩色渲染初步观察辐射强度分布。2.2 辐射亮温计算理解K1/K2常数辐射亮度温度Tb是地表温度反演的基础计算公式为Tb K2 / ln(K1/Lλ 1)其中Lλ辐射定标后的亮度值单位W/(m²·sr·μm)K1/K2Landsat8特有的定标常数波段10K1774.89K21321.08波段11K1480.89K21201.14在ENVI中实现该计算打开Band Math工具输入公式1321.08/alog(774.89/b1 1)指定b1为辐射定标后的波段输出结果单位即为开尔文K提示不同传感器版本的K1/K2值可能微调建议通过NASA官网查询最新参数3. 地表比辐射率计算破解发射率难题地表比辐射率ε反映不同地物发射热辐射的能力其计算需要结合植被覆盖信息。完整流程如下3.1 NDVI计算植被健康度指标对可见光波段B4红波段B5近红外波段进行辐射定标使用Transform NDVI工具计算植被指数NDVI (NIR - Red) / (NIR Red)通过Statistics Compute Statistics验证数值范围理想值应在-0.2到0.8之间3.2 植被覆盖度FVC估算基于NDVI结果采用混合像元分解模型计算植被覆盖比例FVC (NDVI - NDVI_soil) / (NDVI_veg - NDVI_soil)经验阈值建议NDVI_soil取累计分布5%处的值如-0.2NDVI_veg取累计分布95%处的值如0.6在ENVI中可通过Band Math实现(b1 lt -0.2)*0 (b1 gt 0.6)*1 (b1 ge -0.2 and b1 le 0.6)*((b10.2)/(0.60.2))3.3 比辐射率最终计算结合FVC结果采用面积加权法计算每个像元的比辐射率ε 0.004*FVC 0.986该公式基于以下假设植被比辐射率0.99裸土比辐射率0.986建筑表面比辐射率0.9704. 单窗算法实现温度反演核心步骤4.1 大气参数获取单窗算法需要三个关键大气参数大气透射率τ通过NASA的Atmcorr工具输入成像时间、经纬度获取大气平均作用温度Ta可通过近地面气温估算Ta ≈ 0.8*T0 16近地面气温T0从气象站数据或再分析资料获取实战技巧若无实时大气数据可采用历史同期平均值但会引入约1-2℃误差4.2 中间参数计算构建两个辅助参数C和DC ε * τ D (1-τ) * [1 (1-ε)*τ]在ENVI中使用Band Math分别计算C计算式b1 * 0.8假设τ0.8D计算式(1-0.8)*(1 (1-b1)*0.8)4.3 温度反演最终公式将各参数代入单窗算法核心公式LST [a*(1-C-D) (b*(1-C-D)CD)*Tb D*Ta] / C其中经验系数当温度在0-70℃范围时a -62.735657b 0.434036ENVI实现步骤新建Band Math窗口输入完整公式(-62.735657*(1-b1-b2) (0.434036*(1-b1-b2)b1b2)*b3 b2*300)/b1假设Ta300K指定变量b1 Cb2 Db3 辐射亮温Tb4.4 温度单位转换与结果验证最后一步将开尔文温度转为摄氏度LST_Celsius LST_Kelvin - 273.15验证反演结果合理性的三个指标水体温度应在15-30℃之间夏季城市建筑区温度通常比郊区高3-8℃植被覆盖区呈现明显的低温特征5. 结果分析与常见问题排查5.1 典型错误与修正方案问题现象可能原因解决方案温度值异常高100℃辐射定标错误检查K1/K2值是否匹配波段图像出现条带噪声数据下载不完整重新下载并校验MD5值水体温度偏高比辐射率设置不当调整水体区域ε至0.96-0.98城市区域温度偏低大气参数不准确使用实时大气廓线数据5.2 进阶优化方向劈窗算法同时利用B10和B11波段可减少大气影响LST A0 A1*T10 A2*T11系数A0-A2需通过大气模拟确定时空融合技术结合多时相数据提升空间分辨率机器学习方法利用随机森林等算法建立NDVI-温度统计关系在完成首次反演后建议保存处理流程为ENVI的Model文件后续只需更换输入数据即可批量处理。对于城市热岛研究可进一步将温度数据导出至GIS软件进行空间格局分析和热点检测。