地表温度反演算法.doc

地表温度作为地球环境分析的重要指标, 而遥感技术作为现代重要的对地观测手段, 使得基于遥感图像的地表温度反演的研究越来越多。主要的地表温度反演方法有:大气校正法,单窗算法,单通道法等等。本文介绍用辐射传输方程法对地表温度进行反演。技术流程: 例子数据为 2002 年9月2 日的襄樊市 Landsat ETM+ 数据。根据数据的特点以及地表温度反演研究的技术要求,采用的技术路线为:先对 Landsat ETM+ 数据进行预处理:数据读取、辐射定标、大气校正、襄樊区域裁剪, 利用大气校正,即: 辐射传输方程法对其影像热红外波段数据进行操作反演, 实现襄樊市地区的地表真实温度的反演研究。具体的处理流程如下: 具体的实现步骤如下: 第一步:准备数据热红外数据使用的是 Landsat 的第六波段,已经做了传感器定标、几何校正、工程区裁剪,详细流程参考上面的流程图。文件为 TM6-rad-subset-jz- 。由 TM 影像( 已经过大气校正) 生成的 NDVI 数据, 已经利用主菜单->Basic Tools->Resize Data ( SFatial/SFectral )重采样为 60 米分辨率,与 TMi6 数据保持一致,文件名为: TM-NDVI- 。第二步:地表比辐射率计算物体的比辐射率是物体向外辐射电磁波的能力表征。它不仅依赖于地表物体的组成,而且与物体的表面状态( 表面粗糙度等) 及物理性质( 介电常数、含水量等) 有关,并随着所测定的波长和观测角度等因素有关。在大尺度上对比辐射率精确测量的难度很大, 目前只是基于某些假设获得比辐射率的相对值,本文主要根据可见光和近红外光谱信息来估计比辐射率。(一) 植被覆盖度计算计算植被覆盖度 Fv 采用的是混合像元分解法, 将整景影像的地类大致分为水体、植被和建筑,具体的计算公式如下: F V= (NDVI- NDVI S)/ (NDVI V- NDVI S)(2) 其中, NDVI 为归一化差异植被指数,取 NDVI V= 和 NDVI S= , 且有, 当某个像元的 NDVI 大于 时, F V 取值为 1 ;当 NDVI 小于 ,F V 取值为 0。利用 ENVI 主菜单->Basic Tools->Band Math ,在公式输入栏中输入: ( b1 gt ) *1+ ( b1 lt 0.) *0+ ( b1 ge0 and b1 le )*( (b1-)/(-) ) b1: 选择 NDVI 图像得到植被盖度图像。(二) 地表比辐射率计算根据前人的研究, 将遥感影像分为水体、城镇和自然表面 3 种类型。本专题采取以下方法计算研究区地表比辐射率:水体像元的比辐射率赋值为 ,自然表面和城镇像元的比辐射率估算则分别根据下式( 3 )( 4 )进行计算: ε surface = + V- V 2(3) ε building = + V- V 2(4) 式中, ε surface 和ε building 分别代表自然表面像元和城镇像元的比辐射率。利用 ENVI 主菜单->Basic Tools->B