DB52∕T 1373-2018 极轨卫星遥感监测地表温度(贵州省) 18页

ICS07.040M50DB52贵州省地方标准DB52/T1373—2018极轨卫星遥感监测地表温度Theremotesensingmonitoringofpolar-orbitingsatellitelandsurfacetemperature2018-11-28发布2019-05-28实施贵州省市场监督管理局发布\n\nDB52/T1373—2018目次前言..............................................................................II1范围..............................................................................12术语和定义........................................................................13监测方法..........................................................................3附录A（资料性附录）气象卫星/多光谱通道的扫描辐射仪(NOAA/AVHRR)地表温度监测通道相关参数表...................................................................................8附录B（资料性附录）地球观测系统/中分辨率成像光谱仪（EOS/MODIS）地表温度监测通道相关参数表...................................................................................9附录C（资料性附录）风云三号气象卫星/中分辨率成像光谱仪（FY3A/MERSI）地表温度监测通道相关参数表............................................................................10I\n\nDB52/T1373—2018前言本标准按照GB/T1.1-2009《标准化工作导则第1部分：标准的结构和编写》给出的规则起草。请注意：本标准的某些内容可能涉及专利内容，本文件的发布机构不承担识别这些专利的责任。本标准由贵州省机械电子产品质量监督检验院提出。本标准由贵州省质量技术监督局归口。本标准起草单位：贵州省机械电子产品质量监督检验院、遵义市产品质量检验检测院。本标准主要起草人：王海、叶鑫、郜愿愿、潘积文、彭菲菲、王玉琴、禹忠、付悦、戚婉君、谭克刚、高元淑、郑殿民。II\n\nDB52/T1373—2018极轨卫星遥感监测地表温度1范围本标准规定地表温度遥感监测的术语和定义、监测方法等。本标准适用于贵州省气象卫星/多光谱通道的扫描辐射仪(NOAA/AVHRR)、地球观测系统/中分辨率成像光谱仪（EOS/MODIS）、风云三号气象卫星/中分辨率成像光谱仪（FY3A/MERSI）卫星星载光学仪器的地表温度遥感监测。2术语和定义下列术语和定义适用于本文件。2.1极轨卫星离地面约840公里的轨道上运行，运行轨道通过地球的南北极，与太阳同步的卫星。2.2遥感监测通过航空或卫星等收集环境的电磁波信息对远离的环境目标进行监测识别环境质量状况的技术。2.3气象卫星（NOAA）是近极地、与太阳同步的卫星，高度为833km～870km，轨道倾角98.7°，成像周期12小时。2.4多光谱通道的扫描辐射仪（AVHRR）一种多光谱通道的扫描辐射仪，两条轨道可以覆盖我国大部分国土，三条轨道可完全覆盖我国全部国土。2.5地球观测系统（EOS）由多颗卫星组成和为实行多学科（大气、海洋、陆面、生物、化学等）综合研究而建立的全球卫星观测体系。1\nDB52/T1373—20182.6中等分辨率成像分光仪（MODIS）一台使用可见光和红外线对地球（陆地、海洋和大气）进行一系列综合性观测的成像分光仪。2.7FY-3A卫星我国第二代极轨气象卫星“风云三号气象卫星”的第一颗卫星。2.8中分辨率成像光谱仪（MERSI）探测来自大气系统的电磁辐射，通过成像，可以实现植被、地表覆盖分类以及积雪覆盖等陆表特性全球遥感监测。2.9地表温度landsurfacetemperature地球表面的温度,包括陆地表面和海洋表面，是地表各种地物的综合温度。2.10单通道算法single-channelalgorithm针对卫星传感器上单独一个热红外波段而提出来的单波段温度反演算法。该算法以辐射传导方程为基础，结合大气和比辐射率的修正，从而得出地表温度。2.11分裂窗算法split-windowalgorithm两个不同的通道在同一个大气窗口上吸收的差异来消除大气的影响，由两通道的亮温线性运算而得。2.12归一化植被指数normalizeddifferencevegetationIndex反映土地覆盖植被状况的一种遥感指标，定义为近红外通道与可见光通道反射率之差与之和的商。2.13植被覆盖度vegetationcover单位面积内植被的垂直投影面积所占百分比。2.14大气透过率atmospherictransmissibility遥感传感器接收到的热辐射能与地表真实辐射能的比值。2\nDB52/T1373—20182.15地表比辐射率surfaceemissivity物体与黑体之间辐射能量对比的一种度量，其定义为物体辐射与同温度下黑体辐射的比值，与物体的温度、波长、水分含量，表面粗糙度，植被的结构以及观测角度等有关。3监测方法3.1分裂窗算法和单通道算法3.1.1气象卫星/多光谱通道的扫描辐射仪(NOAA/AVHRR)、地球观测系统/中分辨率成像光谱（EOS/MODIS）卫星监测地表温度采用分裂窗算法。3.1.2风云三号气象卫星/中分辨率成像光谱仪（FY3A/MERSI）卫星监测地表温度采用单通道算法。3.2基本要求提供卫星监测的数据源。3.3分裂窗算法3.3.1分裂窗地表温度算法按（1）式计算。TSA0A1T热红外1-A2T热红外2...............................(1)式中：TS——地表温度；T热红外1、T热红外2——热红外通道的亮度温度（AVHRR的4、5通道数据源；MODIS的31、32通道数据源）；A0、A1、A2——算法的参数变量（分别由AVHRR、MODIS提供的通道数据源计算确定）。3.3.2AVHRRA0、A1、A2的参数确定基础计算数据来源：1、2为卫星可见光、近红外通道反射率，AVHRR的1、2通道。计算过程：21NDVI.....................................(2)21式中：NDVI——归一化植被指数。2NDVI2.0Pv...................................(3).009式中：Pv——像元的植被构成比例。地表比辐射率计算如下：3\nDB52/T1373—2018.0981-.00421，NDVI2.0.0971.0018Pv，0.2NDVI0.5.........................(4)0.989，NDVI0.5.0003.00291，NDVI2.0.00061Pv，2.0NDVI5.0.........................(5)0，NDVI5.0式中：——地表比辐射率；——通道4、通道5的地表比辐射率差。中间参数计算如下：2P1.0156161.0482..........................(6)2M.626.398138.33..........................(7)A0、A1、A2的参数:A0.1274T4P1..................................(8)A11MP2...................................(9)A2MP2....................................(10)式中：T4——通道4的温度。3.3.3MODISA0、A1、A2的参数确定基础计算数据来源：.0940、.0865为地表反射率，MODIS卫星第19和2波段。计算过程：2-ln0.9400.865...............................(11)式中：-2——大气水汽含量（g·cm）；、——常量，.002、.0651。地表比辐射率的计算如下：iPvRviv1PvRsisd.............................(12)NDVINDVIminPv................................(13)NDVImax-NDVImin4\nDB52/T1373—20180，Pv0或Pv1.0003796Pv，0Pv5.0d..........................(14).0001898，Pv5.0.00037961Pv，5.0Pv1式中：i——波段i的地表比辐射率；计算得到的iiv时，则取iiv，iis时，则取iis；对于水体像元，iiw31.099683、32.0992324；iv——植被在波段i的地表比辐射率，31v=0.98672、32v=0.98990；is——裸土在波段i的地表比辐射率，31s=0.96767、32s=0.97790；Pv——像元的植被覆盖度；NDVI——归一化植被指数,通常贵州省NDVImin=0.79，NDVImax=0.87；Rv、Rs——植被、裸土的辐射比率，Rv.092762.007033Pv，Rs.099782.008362Pv；d——热辐射校正。注：wi是指若是水体像元时iiw（常数代号）（31.099683，32.0992324），此时i不由公式计算，直接取常数31或32。中间参数计算如下：Ciii......................................(15)Di1i11ii..............................(16)式中：i——卫星的热红外波段（MODIS的31、32通道值）；i——视角为的大气透过率，见表1；i——波段i的地表比辐射率。A0、A1、A2的参数:A[D1(CD）/(DCDC)]a[D1(CD)/(DCDC)]a032313132313132313132323231313232...(17)A[D1(CD）/(DCDC)]a[D1(CD)/(DCDC)]a.....(18)132313132313132313132323231313232AD/(DCDC)[D1(CD)/(DCDC)]b.............(19)231323131323132323231313232注1：a31=-64.60363、b31=0.440817、a32=-68.72575、b32=0.473453。注2：大气透过率估算方程与校正：大气透过率与大气水汽含量之间的关系，可以通过大气辐射传输模型模拟来确定，见表1，根据式(16)得到大气水汽含量后，按照表1所示的关系方程即可得到大气透过率。大气透过率还受到遥感器视角和大气剖面温度的影响，还需要进行校正，见表2。5\nDB52/T1373—2018表1MODIS第31和32波段的大气透过率估计方程22大气水汽含量(g/cm)大气透过率方程SEERF31.101636-.0103460.003440.991361605.70.4～2.032.102144-.0139270.003730.994392480.131.111795-.0155360.001560.9998398912.92.0～4.032.109361-0.179800.002810.9995840526.831.077313-0.074040.020050.96637316.14.0～5.432.065166-0.073540.020890.96313287.431.1101089-0.096560.002890.992781100.7冬季0.4～1.432.097022-0.080570.012490.8407473.9表2大气透过率的温度校正函数波段大气透过率方程温度区间(K)31T.008T31318MODIS3131T.005.000325T31278278T3131831T.005T3127832T.0095T32318MODIS3232T.0065.0004T32278278T3231832T.0065T322783.3MERSI单通道算法基础计算数据来源：.0940、.0865——地表反射率，MERSI卫星第18和第16波段；Lsensor——传感器接收到的热辐射亮度；Tsensor——传感器反映的亮度温度。计算过程：-2——大气水汽含量（g·cm）的计算和MODIS相同，见公式(11)。1、2、3——大气水汽含量的函数计算：321.001652-.000331.013380.099251....................(20)322-.010768-.033642-.191671.023652...................(21)323-0.05494.039078.099000-.008928...................(22)中间参数计算如下：14c2Lsensor12Lsensor...........................(23)Tsensorc1-LsensorTsensor.................................(24)6\nDB52/T1373—2018式中：——有效波长,=11.25m；84214c1、c2——常数，c1.11910410Wmmsr，c2.14387710mK；Ts——地表温度计算：1Ts1Lsensor23.........................(25)7\nDB52/T1373—2018AA附录A（资料性附录）气象卫星/多光谱通道的扫描辐射仪(NOAA/AVHRR)地表温度监测通道相关参数表表A.1气象卫星/多光谱通道的扫描辐射仪(NOAA/AVHRR)地表温度监测通道相关参数表通道波长（μm）波段星下点分辨率（km）10.58～0.68可见光1.120.70～1.10近红外1.1410.30～11.30远红外1.1511.50～12.50远红外1.18\nDB52/T1373—2018BB附录B（资料性附录）地球观测系统/中分辨率成像光谱仪（EOS/MODIS）地表温度监测通道相关参数表表B.1地球观测系统/中分辨率成像光谱仪（EOS/MODIS）地表温度监测通道相关参数表通道波长（μm）星下点分辨率（km）10.62～0.670.2520.84～0.870.25190.92～0.9713110.78～11.2813211.77～12.2719\nDB52/T1373—2018CC附录C（资料性附录）风云三号气象卫星/中分辨率成像光谱仪（FY3A/MERSI）地表温度监测通道相关参数表表C.1风云三号气象卫星/中分辨率成像光谱仪（FY3A/MERSI）地表温度监测通道相关参数表通道波长（μm）星下点分辨率（km）30.60～0.700.2540.82～0.920.2558.75～13.750.25160.85～0.891180.92～0.961_________________________________10\n\nDB52/T1373-2018