极轨卫星遥感监测第9部分地表温度.doc

ICS 07.060A 47DB 21地 方 标 准DB 21/ T2015—2012极轨卫星遥感监测 第 9 部分：地表温度The remote sensing monitoring of polar-orbiting satelliteNinth part:Land Surface Temperature（报批稿）（本稿完成日期：2012.1.302012 – 08 – 13 发布 2012 – 09 – 13 实施辽 宁 省 质 量 技 术 监 督 局       发 布DB21/ XXXXX—XXXXI目  次前言 ..................................................................................II引言 .................................................................................III1　 范围 ................................................................................12　 术语和定义 ..........................................................................13　 监测方法 ............................................................................23.1　 分裂窗算法 ......................................................................23.1.1　 AVHRR 分裂窗算法 ............................................................23.1.2　 MODIS 分裂窗算法 ............................................................33.2　 MERSI 单通道算法 ................................................................63.2.1　 中间参数计算 ................................................................6附录 A（资料性附录）　 NOAA/AVHRR 地表温度监测通道相关参数表 .............................7附录 B（资料性附录）　 EOS/MODIS 地表温度监测通道相关参数表 ..............................8附录 C（资料性附录）　 FY3A/MERSI 地表温度监测通道相关参数表 .............................9表 1　 MODIS 第 31 和 32 波段的大气透过率估计方程 ..........................................4表 2　大气透过率的温度校 正函数 .........................................................5DB21/ XXXXX—XXXXII前  言本部分依据GB/T 1.1-2009的编写规则起草本部分的附录A、附录B、附录C为资料性附录。

本部分由辽宁省气象局归口本部分起草单位：中国气象局沈阳大气环境研究所本部分主要起草人：武晋雯、张玉书、冯锐、纪瑞鹏、张淑杰、于文颖、姜鹏DB21/ XXXXX—XXXXIII引  言地表温度是研究地表与大气之间能量与物质交换的重要参数，在气候、水文、生态、环境等研究中都有重要应用利用卫星遥感技术反演的大面积地表温度信息已被广泛的应用在农作物估产、作物长势和农业旱情监测等方面地表温度的反演精度直接影响到遥感技术在生态环境上应用的效果地表温度遥感反演有单通道法、单通道多角度法、多通道多角度法、分裂窗法，而目前比较成熟的分裂窗法算法，又演变成大约17种之多因此需要制定极轨卫星遥感监测地方标准，以便提高辽宁地表温度遥感监测的精度和效率DB21/ XXXXX—XXXX1极轨卫星遥感监测 第 9 部分：地表温度1　范围本部分是极轨卫星遥感监测第9部分，其它8个部分是术语、干旱灾害、洪涝灾害、森林火灾、作物长势、大雾、沙尘、海冰均已发布，本部分规定地表温度遥感监测算法本部分适用于辽宁地区NOAA/AVHRR、EOS/MODIS、 FY3A/MERSI卫星星载光学仪器的地表温度遥感监测2　术语和定义2.1　地表温度 Land Surface Temperature地表温度就是指地球表面的温度,包括陆地表面和海洋表面，是地表各种地物的综合温度。

2.2　单通道算法 Single-channel Algorithm单通道算法是针对卫星传感器上单独一个热红外波段而提出来的单波段温度反演算法该算法以辐射传导方程为基础，结合大气和比辐射率的修正，从而得出地表温度2.3　分裂窗算法 Split-Window Algorithm分裂窗方法是两个不同的通道在同一个大气窗口上吸收的差异来消除大气的影响，由两通道的亮温线性运算而得2.4　归一化植被指数 Normalized Difference Vegetation Index归一化植被指数是反映土地覆盖植被状况的一种遥感指标，定义为近红外通道与可见光通道反射率之差与之和的商2.5　植被覆盖度 Vegetation cover植被覆盖度是指单位面积内植被的垂直投影面积所占百分比2.6　大气透过率 Atmospheric TransmissibilityDB21/ XXXXX—XXXX2大气透过率是遥感传感器接收到的热辐射能与地表真实辐射能的比值2.7　地表比辐射率 Surface Emissivity是物体与黑体之间辐射能量对比的一种度量，其定义为物体辐射与同温度下黑体辐射的比值，与物体的温度、波长、水分含量，表面粗糙度，植被的结构以及观测角度等有关。

3　监测方法应用NOAA/AVHRR、EOS/MODIS卫星监测地表温度采用分裂窗算法,应用FY3A/MERSI卫星监测地表温度采用单通道算法3.1　分裂窗算法...............................(1)外210ST-AT式中：—地表温度；ST、 —热红外通道的亮度温度（MODIS的31、32通道值，AVHRR的4、5通道；）外12外、 、 —分裂窗算法的参数变量，AVHRR见式(2)～(4)、MODIS见式(11)～(13)0A3.1.1　AVHRR 分裂窗算法...................................(2)1PT1.274A0....................................(3)/M......................................(4)2式中：、 、 —式(1)中的参数变量；0A12—通道4的亮度温度；T、 —中间参数，见式(5)～(10)MP3.1.1.1　中间参数计算...........................(5)2/0.48/ε10.56DB21/ XXXXX—XXXX3...........................(6)2/38./ε13.96.2M式中：—地表比辐射率，见式(7)；—通道4、通道5的地表比辐射率差，见式(8)。

3.1.1.2　地表比辐射率计算..........................(7)0.5NDVI0.980.5I外2.1P7.,4ρ.εv，.........................(8).I.I.6,ρ3Δεv1，...................................(9)0.92NDVPv‘....................................(10)12I式中：—像元的植被构成比例；’vP—归一化植被指数；NDVI、 —为卫星可见光、近红外通道反射率，AVHRR的1、2通道，MERSI的3、4通道123.1.2　MODIS 分裂窗算法...........(11)3213232313213231320 )]/()([)]/(([ aCDCDaCDA ）...........(12)11 b............(13)321323231231322 )]/()([)/(式中：、 、 —式(1)中的参数变量；0A12、 、 、 —常量， =-64.60363、 =0.440817、 =-3ab331a31b32a68.72575、 =0.473453；2DB21/ XXXXX—XXXX4、 —中间参数，见式(14)、式(15)。

iCiD3.1.2.1　中间参数计算......................................(14)iiC..............................(15)θτε1θτDiii 式中：—卫星的热红外波段（MODIS 的 31、32 通道值）；i—视角为 的大气透过率；i—波段 的地表比辐射率ii3.1.2.2　大气透过率的计算...................................(16)2865.094ln-式中：—大气水汽含量（g·cm-2）；、 —常量， =0.02、 =0.651；、 —MODIS卫星第19和2波段的地表反射率,MERSI卫星第18和16波段的地表反射率940.865.3.1.2.2.1　大气透过率估算方程与校正大气透过率与大气水汽含量之间的关系，可以通过大气辐射传输模型模拟来确定，见表1，根据式(16)得到大气水汽含量后，按照表1所示的关系方程即可得到大气透过率大气透过率还受到遥感器视角和大气剖面温度的影响，还需要进行校正，见表2表 1　MODIS 第 31 和 32 波段的大气透过率估计方程大气水汽含量(g/cm 2)大气透过率方程 SEE R2 F0.4～2.0 10346.-6.1392720.003440.003730.991360.994391605.72480.12.0～4.0 5.7.18-32 0.001560.002810.999830.9995898912.940526.8表 1（续）DB21/ XXXXX—XXXX54.0～5.4 .074-31.0315652 0.020050.020890.966370.96313316.1287.4冬季0.4～1.4.968.-3 0.002890.012490.992780.840741100.773.9表 2　大气透过率的温度校正函数波段 大气透过率方程 温度区间(K)MODIS 31 08.31T275.3131 . 318T2731MODIS 32 09322784.65.332TT832T73.1.2.3　地表比辐射率的计算..............................(17)dRPisvivi 1................................(18)minaxNDVIIv...........................(19)1P,0.510.37968,.。