从Landsat8影像反演地表温度的劈窗算法研究.docx

第 39 卷第 4 期2014 年 8 月测绘地理信息Journal of GeomaticsVol.39 No.4 Aug. 2014文章编号：2095-6045（2014）04-0073-05文献标志码：B从 Landsat8 影像反演地表温度的劈窗算法研究杨 槐 1（1 厦门精图信息技术股份有限公司，厦门市吕岭路 1819 号，361008）Research of Split_Window Algorithm for Retrieval of Land Surface Temperature from Landsat8YANG Huai1（1 Xiamen Kingtop Information Technology Co.，Ltd.，1819 Lüling Road，Xiamen 361008，China）摘要：以 厦门市为研究区域 ，以 2013 年的 Landsat8 为 研究数据 ，研 究亮度温度反演以及地表温度反演的劈窗 算法 ，推 导出该算法的演算过程 ，根 据 Landsat8 的特性， 确定该算法的参数取值。通过 ENVI/IDL 软件编程实现 该算法 ，采 用厦门市遥感数据反演地表温度 ，并 且利用 厦门市气象局实测温度验证地表温度反演结果的有效 性。关键词：陆地卫星；亮度温度；地表温度；劈窗算法 中图法分类号：P237Abstract：Taking Xiamen as the study area and Landsat8 in 2013 as the research data，this paper studies the split_window algorithm for retrieval of land surface temperature，specifies the derivation process，determines the parameter selection based on characteristics of Landsat8. It also programs and implements with the ENVI/IDL software，verifies the effectiveness of the surface inversion results by the measured temperature from Xia? men Bureau of Meteorology.Key words：landsat；bright temperature；land surface tempera? ture；split_window algorithm地表温度是监测地球资源环境动态变化的重要 指标之一，对水文、生态、环境和生物地球化学等研 究有重要意义[1]。 因此 ，获 取区域地表温度空间差 异，并进而分析其对区域资源环境变化的影响，是区 域资源环境动态监测的重要内容。目前，针对不同 的热红外波段提出了多种地表温度反演方法，如劈 窗算法[2_3] 、单窗算法[4] 、单通道算法[5] 等。其中，劈窗 算 法 主 要 是 针 对 NOAA/AVHRR 的 第 4、5 热 红 外 波段提出来的 ，并 被首先运用到海面温度反演 ，有 许多不同改进算法应用于其他多个热红外波段进行 地表温度的反演 ，例 如 MODIS 影像数据等。 对于Landsat8 数据，由于波谱传感器不同，两个热红外波 段存在温度差异大等特点，上述算法直接应用则不 适合。 为了解决上述问题 ，笔 者针对 Landsat8 数据 的特点 ，把 用于 AVHRR 的劈窗算法改进成适合于 Landsat8 数据的劈窗算法， 并且采用 ENVI/IDL 实 现了该算法。1 反演地表温度的劈窗算法笔者提出的算法是在现有劈窗算法的基础上， 针对 Landsat8 的波段特征 ，先 对 Planck 方程进行线 性简化，接着给出了地表比辐射率的估计，并且通过 建立大气水汽含量与热波段透过率的关系，由此可 以反演得到地表温度。1.1 算法流程算法流程如图 1 所示。包括了从数据获取到专 题图制作等工作。几何校正辐射校正专题图制作获取卫星影像Landsat8（2013.08.04）影像校正亮度温度反演植被指数计算水体指数计算矢量数据统计分析地表温度反演建立掩膜图 1 算法流程图1.2 关键内容1.2.1 亮度温度反演根据普朗克定律将辐射亮度值转化为亮度温 度，计算公式如下[4]：T =K2 ln( K1 + 1 )（1）Lλ式中，T 为亮度温度，单位为 K ；Lλ 为辐射亮度值； K1 和 K2 参数可从影像头文件中获取。通常 T 可以 转换为摄氏度温度 T 来表示：T = T - 273.15（2）201