基于单窗算法反演地表温度的ENVI操作教程.docxVIP

单窗算法反演地表温度教程 1.1算法原理 1.1.1单窗算法 单窗算法(MW算法)是覃志豪于2001年提出的针对TM数据只 有一个热红外波段的地面温度反演算法。经过众多学者验证,单窗算 法具有很高的反演精度，且同样适用于ETM+和land s at8数据。公 式如下： T* =(?6(1-C6 -￡)6) + (Z?6(l-C6 -Z)6) + C6 +D6)Tww,r-D67^)/C6 式中上S T为地表温度(K),7\.g?是传感器上的亮度温度(K),几是 大气平均温度(K) ;a、b为参考系数，当地表温度为0-7 0°C时,a =?6 7.3 5 5 351, b = 0.4 5 860 6 ； C 、D为中间变量,计算公式为： C = Z) = (l-r)[l + (l-f)r] 式中，为地表比辐射率,r为地而到传感器的大气总透射率。 因此单窗算法反演地表温度的关键是计算得到亮度温度Ts^s。、地表 比辐射率￡、大气透射率r和大气平均作用温度TaQ 1.1.2参数计算 1.1.2.1辐射亮温计算 利用Pla n ck公式将图像像元对应传感器辐射强度值转换为对应的 亮度温度值。公式如下 ln（l+才） 式中，Ts沁为亮度温度值;乙影像预处理后得到的光谱辐射值，单 位为“，/（〃，』??““），ki 、K2为常量，可由数据头文件获取。 计算图像辐射亮温之前，需采用辐射定标参数将像元灰度值DN 转换为热辐射强度值，公式如下： La =虫 乂DN+A] 式中,M为增益参数3厶为偏移参数，该参数可直接在影像通文件 数据中获取，且ENVI软件中己经集成,不需要自己在查找。 1.1.2.2地表比辐射率计算 根据覃志豪针对TM影像提出的混合像元分解法来确定区域地 表福辐射率。对于城市区域，我们简单的将其分为水体、自然表面和 建筑表面三种，因此针对混合像元尺度上的地表比辐射率通过下式来 估算： S = PvRvSv + （1 - E ）豔 +￡ 式中Q为混合像元的地表比辐射率；Pv为植被覆盖率；Rv为植 被的温度比率;Rm为建筑表而的温度比率;表示植被法地表比辐射 率,乩表示建筑表而的地表比辐射率;ds表示辐射校正项。 根据覃志豪经验公式，￡\ = 0 .986：m=0.97 2。 根据Aster提供的常用地物比辐射率光谱库Lan d s at 8 数据特点，对于B 1 0波段,￡v= 0 .9 8 6 7 2 ;尸0?96767。对于Bl 1 波段,￡=0.9899 0； m=O. 977 5 15o 在地表相对较平整的区域，一般可取d沪0;在地表起伏较大的区 域,d￡可根据植被的构成比例估计。研究表明热辐射相互作用在植被 与裸土分别占一半时达到最大,所以可以根据如下经验公式来估计d 当 P v 0. 5 时,d￡=0.0 0 3 8Pv 当 P v0. 5 时，血=0.00 3 8 (1-Pv) 当 P\=0.5 时，d 8=0.0 0 1 9 RV和Rin的确定采用覃志豪等提岀的经验公式计算： RV=0.93 3 2 +0.058 5 *PV Rm=0. 9 886+0. 1 2 87*PV PV为植被覆盖率可用下式计算： P V= (ND\7- N DVI s )/(NDVIv-N D Vis) 其中NDVI是归一化植被指数,NDV I V、NDV I S分别是植 被和裸土的ND VI值，由于没有详细的区域植被和土壤光谱资料，采 用覃志豪经验值以NDVW=0.7 0和NDVI S=0. 05来进行植被覆盖 度的近似估计。 1.1.2.3大气平均作用温度的计算 对于大气平均温度，覃志豪等也提供了一个经验公式，该公式是 根据Mo d trail所提供的标准大气推导出来的，如下表所示，其中To 为近地表温度(K)。 表3?4大气平均作用温度与地面附近气温的关系固 人气的剖面类型 7；与珀的关系 美国1976年平均大气 =25.9396+0.88045 热带平均人气 Za=17.9769+0.21715 7^ 中纬度夏季半均大气 7； =16.0110+0.92621 70 中纬度冬季平均人气 7； =19.2704+0.91118 一般研究区在重庆的，选择中纬度夏季平均大气。 1.1.2.4大气透射率的计算 由于单窗算法对输入的参数较敏感，因此需要较精确的大气透射 率参数来反演地表温度。影响大气透射率的因素较多，如气压、气温、 气溶胶含量、大气水分含量、03、CO2、CO、NH4等对热辐射传导 均有不同程度的作用，从而使地表的热辐射在大气中的传导产生衰减。 因此,准确的大气透射率求算比较复杂，需要较详细的大气剖面数据， 但在实际应用中很难获取详细的大气剖面数据。研究表明大气水分含 量对大气透射率的变化影响较大,而其它因素由于很难在短时