MODIS数据反演地表温度.docVIP

表1 MODIS部分波段及其参数[14] 波段 光谱范围 信噪比 主要用途 分辨率 1 620～670nm 128 陆地、云边界 250m 2 841～876nm 201 陆地、云边界 250m 19 915～965nm 250 大气水汽 1000m 31 10.780～11.280μm 0.05 地球表面和 云顶温度 1000m 32 11.770～12.270μm 0.05 1000m 劈窗算法介绍 McMillin（1975年）最早提出了劈窗算法，最先是用于海面温度的反演，这种方法是利用2个相邻的热红外窗口大气水汽吸收特性的差异，把海面温度表达成2个热红外窗口亮度温度的线性组合。Price（1984年）最先把劈窗算法推广到陆面温度的反演，通过引入比辐射率改正项来减小因陆地表面比辐射率变化而引起的误差。Becker从理论上证明了用分裂窗技术反演地表温度的可行性，并且第一次从理论上给出了使用分裂窗技术时大气和比辐射率对地表温度反演的影响。Becker和Li根据热辐射传导的地方性特征，提出了著名的局地劈窗算法，已得到了较广泛的应用。Wan和Dozier在Becker和Li的研究基础上，于1996年提出了一种广义的地表温度反演劈窗算法。Sobrino和Becker用Lowtran 7对不同的大气、观测角度以及地表参数进行模拟，得出了各参数的表达式。在这些表达式里，大气和比辐射率的作用是耦合在一起的。而Sobrino等则通过某些近似把这2种作用分开了，通过对大气向下热辐射的近似解和对Planck辐射函数的线性化。覃志豪等推导了劈窗算法，该算法仅需要2个因素来进行地表温度的演算，即大气透过率和地表比辐射率 [15][ 16]。在众多的劈窗算法中，覃志豪等提出的算法由于需要参数少、计算简单且精度较高，被认为是较好的算法之一。本文主要针对这一算法进行介绍。 覃志豪[15]等提出的针对MODIS数据反演地表温度的劈窗算法使用的公式如下: （1） 其中：是地表温度，、分别是MODIS第31、32通道的亮温。、、是系数，分别定义如下： 其中：、、、、、、为中间变量，可迭代消除；、分别为31、32波段的地表比辐射率；、为31、32波段的大气透过率。 参数计算 亮度温度计算 亮温是指辐射出与观测物体相等辐射能量的黑体温度，可以根据普朗克（Planck）公式计算得到，MODIS第31和32波段的亮度温度、由下式计算： （2） 其中: 是MODIS第i(i=31,32)波段的亮度温度, 是波段i的中心波长,针对MODIS的第31和32波段,其值可分别取=11.28m和=12.02m；C1和C2是常量,分别取C1=110-16W·和C2=1.4387685×104 m·K。由于需要注意C1、、的单位转化问题，为了便于计算，可将该式进行简化，设，。则上式可转化为： （2b） 其中，为常量，对于31波段：，； 对于32波段：，。 是MODIS第i(i=31,32)波段的热辐射强度,可根据下式得到： （3） 式中radiance_offset为截距，radiance_scale为斜率，均可以从HDF格式的MODIS图像的头文件中直接查出，DN为遥感图像第31和32波段的实际保存数值。 大气透过率的计算 大气透过率是地表辐射、反射透过大气到达传感器的能量与地表辐射能、反射能的比值,它与大气状况、高度等因素有关。对于热红外波段,最重要的大气变化是大气温