完整word版MODIS数据反演地表温度.doc

(完整word版)MODIS数据反演地表温度 (完整word版)MODIS数据反演地表温度 PAGE / NUMPAGES (完整word版)MODIS数据反演地表温度 表 1 MODIS 部分波段及其参数 [14] 波段 光谱范围 信噪比 主要用途 分辨率 1 620～ 670nm 128 陆地、云界限 250m 2 841 ～876nm 201 陆地、云界限 250m 19 915～ 965nm 250 大气水汽 1000m 31 ～μm 地球表面和 1000m 32 ～μm 云顶温度 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数据反演地表温度的劈窗算法使用的公式以下 : Ts A0 AT A T （ 1） 1 31 2 32 此中： Ts 是地表温度， T31 、 T32 分别是 MODIS 第 31、32 通道的亮温。 A0 、 A1 、 A2 是系数，分别 定义以下： A0 64.60363 E1 68.72575 E2 A1 1 A 0.440817 E1 A2 A 0.473453 E2 A D31 /( D32C31 D31C32 ) E1 D 32 (1 C31 D31) /( D32C31 D31C32 ) E2 D31 (1 C32 D32 ) /( D32C31 D31C32 ) D31 (1 31 )[1 (1 31 ) 31] D32 (1 32 )[1 (1 32 ) 32 ] C31 31 31 C32 32 32 此中： A 、 E1 、 E2 、 D31 、 D32 、 C31 、 C32 为中间变量，可迭代除去； 31 、 32 分别为 31、 32 波 段的地表比辐射率； 31 、 32 为 31、 32 波段的大气透过率。 参数计算 亮度温度计算 亮温是指辐射出与观察物体相等辐射能量的黑体温度， 能够依据普朗克 （ Planck ）公式计算获得， MODIS 第 31 和 32 波段的亮度温度 T 由下式计算： T31 、 32 Ti C 2 C1 ) （ 2） ln(1 i i 5 Ri 此中 : Ti 是 MODIS第 i(i=31,32) 波段的亮度温度 , i 是波段 i 的中心波长 , 针对 MODIS的第 31 和 32 波段 , 其值可分别取 31 m 和 32 m； C1 和 C2 是常量 , 分别取 C1=1× 10-16W· m2 和 C2 =1.4387685 × 104 m· K。因为需要注意 C1、 Ri 、 i 的单位转变问题，为了便于计算，可将该式进行 简化，设 Ki ,2 C2 / i ， K i ,1 C / 5 。则上式可转变为： 1 i Ti K i ,2 （2b） K i ,1 ) ln(1 Ri 此中 K i ,1， Ki ,2 为常量，关于 31 波段： K 31,1 ， K 31,2 ； 2 关于 32 波段： K32,1 474.684780 ， K 32,2 1196.978785 。 Ri 是 MODIS第 i(i=31,32) 波段的热辐射强度 ,