基于单窗算法反演地表温度的ENVI操作教程(word文档良心出品).docx

单窗算法反演地表温度教程

1.1 算法原理

1.1.1 单窗算法

单窗算法（ MW 算法）是覃志豪于2001 年提出的针对TM 数据只有一个热红外波段的地面温度反演算法。经过众多学者验证，单窗算法具有很高的反演精度，且同样适用于 ETM+ 和 landsat 8数据。公

式如下：

T s(a6 (1 C 6 D 6 ) (b6 (1 C6D6 ) C 6 D 6 )T sensor D 6T a ) / C6

式中，LST 为地表温度（K），T sensor是传感器上的亮度温度（ K），T a是大气平均温度（ K）；a、b 为参考系数，当地表温度为 0-70℃时，a = -

67.355351，b = 0.458606；C 、D 为中间变量，计算公式为：

式中，为地表比辐射率，为地面到传感器的大气总透射率。

因此单窗算法反演地表温度的关键是计算得到亮度温度 T senso、地表比辐射率、大气透射率和大气平均作用温度 T a。

1.1.2 参数计算

1.1.

2.1 辐射亮温计算

利用 Planck 公式将图像像元对应传感器辐射强度值转换为对应的亮度温度值。公式如下

式中，T senso为亮度温度值；L影像预处理后得到的光谱辐射值，

单位为w /(m

2

sr m)

，K1 、K2 为常量，可由数据头文件获取。

计算图像辐射亮温之前，需采用辐射定标参数将像元灰度值 DN 转换为热辐射强度值，公式如下：

式中， M L为增益参数， A L为偏移参数，该参数可直接在影像通文件数据中获取，且ENVI 软件中已经集成，不需要自己在查找。

1.1.

2.2 地表比辐射率计算

根据覃志豪针对TM 影像提出的混合像元分解法来确定区域地表

福辐射率。对于城市区域，我们简单的将其分为水体、自然表面和建

筑表面三种，因此针对混合像元尺度上的地表比辐射率通过下式来估算：

式中，为混合像元的地表比辐射率；P V为植被覆盖率； R V为植被的温度比率； R M为建筑表面的温度比率；V 表示植被法地表比辐射率，m 表示建筑表面的地表比辐射率； d 表示辐射校正项。

根据覃志豪经验公式，V =0.986；m=0.972。

根据 Aster 提供的常用地物比辐射率光谱库 Landsat 8 数据特点，对于 B10 波段，V =0.98672；m=0.96767。对于 B11 波段，V =0.98990；m=0.977515。

在地表相对较平整的区域,一般可取dε =0;在地表起伏较大的区域,d ε可根据植被的构成比例估计。研究表明热辐射相互作用在植被

与裸土分别占一半时达到最大 ,所以可以根据如下经验公式来估计 dε：当Pv≤0.5时,d ε=0.003 8Pv

当Pv>0.5 时,d ε=0.003 8 (1-Pv)

当Pv=0.5 时, d ε=0.0019

RV 和 Rm 的确定采用覃志豪等提出的经验公式计算:

RV=0.9332+0.0585*PV

Rm=0.9886+0.1287*PV

PV 为植被覆盖率可用下式计算:

PV=(NDVI-NDVI S)/(NDVI V -NDVI S)

其中 NDVI 是归一化植被指数 ,NDVIV 、NDVIS 分别是植被和

裸土的 NDVI 值,由于没有详细的区域植被和土壤光谱资料 ,采用覃志

豪经验值以 NDVIV=0.70 和 NDVIS=0.05 来进行植被覆盖度的近似估计。

1.1.

2.3 大气平均作用温度的计算

对于大气平均温度，覃志豪等也提供了一个经验公式，该公式是根据 Modtran 所提供的标准大气推导出来的，如下表所示，其中 T0为近地表温度 (K) 。

一般研究区在重庆的，选择中纬度夏季平均大气。

1.1.

2.4 大气透射率的计算

由于单窗算法对输入的参数较敏感，因此需要较精确的大气透射率参数来反演地表温度。影响大气透射率的因素较多,如气压、气温、气溶胶含量、大气水分含量、O3、CO2、CO、NH4 等对热辐射传导均有不同程度的作用,从而使地表的热辐射在大气中的传导产生衰减。因此 ,准确的大气透射率求算比较复杂,需要较详细的大气剖面数据，但在实际应用中很难获取详细的大气剖面数据。研究表明大气水分含量对大气透射率的变化影响较大,而其它因素由于很难在短时间内发生较大变化而对大气透射率的影响较小。因此,水分含量就成为大气

透射率估计的主要考虑因素。

覃志豪等根据大气中的水汽含量动态变化对大气透射率的影响

程度进行了研究，并对大气透射率进行了估计，如表 3.1。其中 :w 是大气水汽含量。以重庆为研究区的研究一般选用是气温较高时的估算

公式。

大气透射率估计方程

水分含量大气透射率估计方相关系数标准误大气剖面

程平方 (R2)差(SEE) w/(g · cm-2)

0.4-1.6τ 6=0.974290-

0.996110.002368 0.08007w

高气温

τ 6=1.031412-

1.6-3.00.998270.002539

0.11536w

0.4-1.6τ 6=0.982007-

0.994630.003340 0.09611w

低气温

τ 6=1.053710-

1.6-3.00.998990.002375

0.14142w

接下来只需要计算出大气水分含量便可计算得到大气透射率。学者杨景梅根据地面、高空气象要输资料，提出利用地面湿度计算整层大气可降水量，由此计算得出整层大气水汽含量的方法（杨景梅邱金恒，用地面湿度参量计算我国整层大气可降水量及有效水汽含量的方法的研究，大气科学， 2002.01）公式如下：

式中，ω为大气水汽含量（g/cm2），W 为整层大气可降水量（cm）；C0 C1 为经验系数，以重庆为研究对象，一般可选取成都的经验值来