ICS 07.040 M 50 DB52 贵 州 省 地 方 标 准 DB52/T 1373—2018 极轨卫星遥感监测地表温度 The remote sensing monitoring of polar-orbiting satellite land surface temperature 2018 - 11 - 28 发布 贵州省市场监督管理局 2019 - 05 - 28 实施 发 布 DB52/T 1373—2018 目 次 前 言 .............................................................................. II 1 范围 .............................................................................. 1 2 术语和定义 ........................................................................ 1 3 监测方法 .......................................................................... 3 附录 A（资料性附录） 气象卫星/多光谱通道的扫描辐射仪(NOAA/AVHRR)地表温度监测通道相关参数 表 ................................................................................... 8 附录 B（资料性附录） 地球观测系统/中分辨率成像光谱仪（EOS/MODIS）地表温度监测通道相关参数 表 ................................................................................... 9 附录 C（资料性附录） 风云三号气象卫星/中分辨率成像光谱仪（FY3A/MERSI）地表温度监测通道相 关参数表 ............................................................................ 10 I DB52/T 1373—2018 前 言 本标准按照GB/T 1.1-2009《标准化工作导则 第1部分：标准的结构和编写》给出的规则起草。 请注意：本标准的某些内容可能涉及专利内容，本文件的发布机构不承担识别这些专利的责任。 本标准由贵州省机械电子产品质量监督检验院提出。 本标准由贵州省质量技术监督局归口。 本标准起草单位：贵州省机械电子产品质量监督检验院、遵义市产品质量检验检测院。 本标准主要起草人：王海、叶鑫、郜愿愿、潘积文、彭菲菲、王玉琴、禹忠、付悦、戚婉君、谭克 刚、高元淑、郑殿民。 II DB52/T 1373—2018 极轨卫星遥感监测地表温度 1 范围 本标准规定地表温度遥感监测的术语和定义、监测方法等。 本标准适用于贵州省气象卫星/多光谱通道的扫描辐射仪(NOAA/AVHRR)、地球观测系统/中分辨率成 像光谱仪（EOS/MODIS）、风云三号气象卫星/中分辨率成像光谱仪（FY3A/MERSI）卫星星载光学仪器的 地表温度遥感监测。 2 术语和定义 下列术语和定义适用于本文件。 2.1 极轨卫星 离地面约840公里的轨道上运行，运行轨道通过地球的南北极，与太阳同步的卫星。 2.2 遥感监测 通过航空或卫星等收集环境的电磁波信息对远离的环境目标进行监测识别环境质量状况的技术。 2.3 气象卫星（NOAA） 是近极地、与太阳同步的卫星，高度为833 km～870 km，轨道倾角98.7°，成像周期12小时。 2.4 多光谱通道的扫描辐射仪（AVHRR） 一种多光谱通道的扫描辐射仪，两条轨道可以覆盖我国大部分国土，三条轨道可完全覆盖我国全部 国土。 2.5 地球观测系统（EOS） 由多颗卫星组成和为实行多学科（大气、海洋、陆面、生物、化学等）综合研究而建立的全球卫星 观测体系。 1 DB52/T 1373—2018 2.6 中等分辨率成像分光仪（MODIS） 一台使用可见光和红外线对地球（陆地、海洋和大气）进行一系列综合性观测的成像分光仪。 2.7 FY-3A 卫星 我国第二代极轨气象卫星“风云三号气象卫星”的第一颗卫星。 2.8 中分辨率成像光谱仪（MERSI） 探测来自大气系统的电磁辐射，通过成像，可以实现植被、地表覆盖分类以及积雪覆盖等陆表特性 全球遥感监测。 2.9 地表温度 land surface temperature 地球表面的温度,包括陆地表面和海洋表面，是地表各种地物的综合温度。 2.10 单通道算法 single-channel algorithm 针对卫星传感器上单独一个热红外波段而提出来的单波段温度反演算法。该算法以辐射传导方程为 基础，结合大气和比辐射率的修正，从而得出地表温度。 2.11 分裂窗算法 split-window algorithm 两个不同的通道在同一个大气窗口上吸收的差异来消除大气的影响，由两通道的亮温线性运算而 得。 2.12 归一化植被指数 normalized difference vegetation Index 反映土地覆盖植被状况的一种遥感指标，定义为近红外通道与可见光通道反射率之差与之和的商。 2.13 植被覆盖度 vegetation cover 单位面积内植被的垂直投影面积所占百分比。 2.14 大气透过率 atmospheric transmissibility 遥感传感器接收到的热辐射能与地表真实辐射能的比值。 2 DB52/T 1373—2018 2.15 地表比辐射率 surface emissivity 物体与黑体之间辐射能量对比的一种度量，其定义为物体辐射与同温度下黑体辐射的比值，与物体 的温度、波长、水分含量，表面粗糙度，植被的结构以及观测角度等有关。 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.2 AVHRR A0、A1、A2 的参数确定 基础计算数据来源：  1 、  2 为卫星可见光、近红外通道反射率，AVHRR的1、2通道。 计算过程： NDVI   2  1 ..................................... (2)  2  1 式中： NDVI ——归一化植被指数。 Pv  NDVI  0.22 0.09 ................................... (3) 式中： Pv ——像元的植被构成比例。 地表比辐射率计算如下： 3 DB52/T 1373—2018 0.981 - 0.042  1，NDVI  0.2  0.2  NDVI  0.5   0.971  0.018Pv， 0.989，NDVI  0.5     ......................... (4)    0.003  0.029  1，NDVI  0.2  0.2  NDVI  0.5    0.0061  Pv ， 0，NDVI  0.5     ......................... (5)   式中：  ——地表比辐射率；  ——通道4、通道5的地表比辐射率差。 中间参数计算如下： P  1  0.15616 1      0.482   2 .......................... (6) M  6.26  3.98 1      38.33   2 .......................... (7) A0、A1、A2的参数: A0  1.274  T 4P  1 .................................. (8) A1  1  M  P  2 ................................... (9) A2  M  P  2 .................................... (10) 式中： T 4 ——通道4的温度。 3.3.3 MODIS A0、A1、A2 的参数确定 基础计算数据来源：  0.940 、  0.865 为地表反射率，MODIS卫星第19和2波段。 计算过程：  - ln 0.940  0.865    2 ............................... (11) 式中：  ——大气水汽含量（g·cm-2）；  、  ——常量，   0.02 、   0.651 。 地表比辐射率的计算如下： i  PvRviv  1  Pv Rsis  d ............................. (12) Pv  4 NDVI  NDVImin ................................ (13) NDVImax - NDVImin DB52/T 1373—2018 0，Pv  0或Pv  1 0.003796 Pv， 0  Pv  0.5  d   0.001898，Pv  0.5 0.0037961  Pv ， 0.5  Pv  1     ..............