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

单窗算法反演地表温度教程 算法原理 单窗算法 单窗算法（MW算法）是覃志豪于2001年提出的针对TM数据只有一个热红外波段的地面温度反演算法。经过众多学者验证，单窗算法具有很高的反演精度，且同样适用于ETM+和landsat 8数据。公式如下： 式中，LST为地表温度（K），Tsensor是传感器上的亮度温度（K），Ta是大气平均温度（K）；a、b为参考系数，当地表温度为0-70℃时，a = -67.355351，b = 0.458606；C 、D为中间变量，计算公式为： 式中，为地表比辐射率，为地面到传感器的大气总透射率。 因此单窗算法反演地表温度的关键是计算得到亮度温度Tsenso、地表比辐射率、大气透射率和大气平均作用温度Ta。 参数计算 辐射亮温计算 利用Planck公式将图像像元对应传感器辐射强度值转换为对应的亮度温度值。公式如下 式中，Tsenso为亮度温度值；影像预处理后得到的光谱辐射值，单位为，K1 、K2为常量，可由数据头文件获取。 计算图像辐射亮温之前，需采用辐射定标参数将像元灰度值DN转换为热辐射强度值，公式如下： 式中，ML为增益参数，AL为偏移参数，该参数可直接在影像通文件数据中获取，且ENVI软件中已经集成，不需要自己在查找。 地表比辐射率计算 根据覃志豪针对TM影像提出的混合像元分解法来确定区域地表福辐射率。对于城市区域，我们简单的将其分为水体、自然表面和建筑表面三种，因此针对混合像元尺度上的地表比辐射率通过下式来估算： 式中，为混合像元的地表比辐射率；PV为植被覆盖率；RV为植被的温度比率；RM为建筑表面的温度比率；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 当Pv0.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-NDVIS)/(NDVIV-NDVIS) 其中NDVI是归一化植被指数,NDVIV、NDVIS分别是植被和裸土的NDVI值,由于没有详细的区域植被和土壤光谱资料,采用覃志豪经验值以NDVIV=0.70和NDVIS=0.05来进行植被覆盖度的近似估计。 大气平均作用温度的计算 对于大气平均温度，覃志豪等也提供了一个经验公式，该公式是根据Modtran所提供的标准大气推导出来的，如下表所示，其中T0为近地表温度(K)。 一般研究区在重庆的，选择中纬度夏季平均大气。 大气透射率的计算 由于单窗算法对输入的参数较敏感，因此需要较精确的大气透射率参数来反演地表温度。影响大气透射率的因素较多,如气压、气温、气溶胶含量、大气水分含量、O3、CO2、CO、NH4等对热辐射传导均有不同程度的作用,从而使地表的热辐射在大气中的传导产生衰减。因此,准确的大气透射率求算比较复杂,需要较详细的大气剖面数据，但在实际应用中很难获取详细的大气剖面数据。研究表明大气水分含量对大气透射率的变化影响较大,而其它因素由于很难在短时间内发生较大变化而对大气透射率的影响较小。因此,水分含量就成为大气透射率估计的主要考虑因素。 覃志豪等根据大气中的水汽含量动态变化对大气透射率的影响程度进行了研究，并对大气透射率进行了估计，如表3.1。其中:w是大气水汽含量。以重庆为研究区的研究一般选用是气温较高时的估算公式。 大气透射率估计方程 大气剖面 水分含量w/(g·cm-2) 大气透射率估计方程 相关系数平方(R2) 标准误差(SEE) 高气温 0.4-1.6 τ6=0.974290-0.08007w 0.99611 0.002368 1.6-3.0 τ6=1.031412-0.11536w 0.99827 0.002539 低气温 0.4-1.6 τ6=0.982007-0.09611w 0.99463 0.003340 1.6-3.0 τ6=1.053710-0.14142w 0.99899 0.002375 接下来只需要计算出大气水分含量便可计算得到大气透射率。学者杨景梅根据地