中巴经济走廊2013-2018年地表温度高分辨率反演数据集

赵国辉，张耀南，康建芳

1. 中国科学院西北生态环境资源研究院，寒旱区科学大数据中心，兰州 730000

2. 中国科学院大学，北京 100049

3. 中国科学院超级计算兰州分中心，兰州 730000

数据库（集）基本信息简介

数据库（集）组成国家科技基础条件平台“特殊环境特殊功能观测研究台站共享服务平台”项目（Y719H71001）；中国科学院信息化专项“寒旱区环境演变研究‘科技领域云’的建设与应用”项目（XXH13506）。数据集按照年份共分为6个目录，在每个年份目录下按照月份划分为12个子目录，数据命名规则为 CPEC_LST_Xm_PPPRRR _YYYYMMDD.tif，其中 X为分辨率，PPPRRR为 WRS2行列号，YYYYMMDD 为年月日标识，例如CPEC_LST_30m_152042_20170603.tif。

引 言

地表温度（Land Surface Temperature，LST）作为陆地表面热辐射状态的度量，综合了地表与大气相互作用以及大气和陆地之间能量交换的全部结果，是表征地表水热平衡和气候变化的关键参数、地表循环过程研究的重要组成部分。准确地反演区域地表温度时空变化不仅对气候变化、生态保护以及资源评价等相关领域的研究至关重要，而且对环境监测、灾害防治、工程建设等也具有重要应用价值[1]。

传统地表温度监测技术很难获取其时空分布特征数据，在区域研究和应用中受极大的限制，而热红外遥感探测技术能够探测到地表直接发射能量，成为获取区域地温的重要途径[2]。基于遥感手段的陆表参数反演为地表温度时空数据制备提供了直接的技术手段。目前常用的地表温度遥感数据多数来自MODIS卫星，其虽具有较好的区域性和实效性，对于地表变化连续观测具有明显优势，但是无法满足局部区域环境变化监测的空间分辨率。Landsat系列卫星数据凭借其长期连续、全球覆盖、高分辨率、科学的数据存档与分发策略等优势，成为地表特征和地球系统科学研究中最有效的遥感数据之一[3]。自2013年Landsat 8卫星成功发射以来，提供免费的中分辨率数据以及16天的重放周期，成为区域热红外遥感的重要数据源，目前已拥有5年多的热红外波段存档数据，为地表温度反演提供了重要的数据支撑。

目前针对不同的应用领域以及区域，国内外学者发展了不同的地表温度反演模型。为了准确地反演地表温度的时空变化，兼顾模型的适用性，构建了地表温度反演模型集，其中包括单通道反演模型、双通道劈窗反演模型以及数据融合模型。依据遴选标准，自动生成制备了中巴经济走廊区域2013-2018年Landsat 8地表温度遥感反演数据集（图1），以期为中巴经济走廊自然环境变时空化监测研究、工程建设以及社会服务提供基础数据资料。

1 数据采集和处理方法

1.1 基础数据

Landsat 8卫星携带的2个主要传感器陆地成像仪（OLI）和热红外传感器（TIRS），在扫描方式、波段设置、辐射分辨性能等方面较之前的 Landsat系列卫星均有很大的改进，尤其是热红外遥感器TIRS，可探测10.6-11.2 μm和11.5-12.5 μm两个热红外波段的辐射，瞬时视场角15°，地面幅宽185 km，可以更好地对热红外数据进行大气校正，进而能够对地表温度进行准确地反演。为此，本数据集基于Landsat 8的两个热红外通道（第10与11波段），反演制备地表温度产品。基础数据来源于美国USGS EarthExplorer（https://earthexplorer.usgs.gov/）官方网站，获取了自2013年到目前所有的Landsat 8 OLI/TIRS C1 Level-1数据产品（表1），具体行列号见图1。

图1 中巴经济走廊区域2013-2018年地表温度遥感反演数据集覆盖范围

表1 Landsat 8 OLI/TIRS C1 Level-1数据产品描述

1.2 处理方法

本数据集的数据处理流程如图2，具体步骤如下：

图2 中巴经济走廊地表温度数据集制备流程

步骤 1首先对Landsat 8的基础数据按照数据质量进行筛选，把云量大于20%的数据剔除，然后分别对OLI和TIRS传感器数据进行处理。

步骤 2按照Landsat 8数据手册[4]，分别计算TIRS第10、11波段的辐射亮度值Lλ，并基于辐射亮度值计算星上亮度温度BT。依据裂窗协方差比算法[5,10]反演大气水汽含量ω，以修正大气和地表比辐射率的影响，算法如下：

式中τi，τj分别为第10波段和第11波段的大气透过率，εi与εj为其对应的比辐射率，为其移动窗口内第10波段、第11波段的亮温均值，k为窗口内像元的索引，N为窗口大小，c0、c1与c2为基于大气辐射传输模型 MODTRAN和大气廓线数据库 TIGR（Thermodynamic Initial Guess Retrieval database）模拟产生的回归系数。

步骤 3地表发射率是衡量表面散热效率的指标，不同的下垫面有不同的发射率值（0-1）。本数据集采用植被覆盖度加权法[6-7]，基于Landsat 8 OLI第4、5波段获取地表发射率ε：

式中：植被组分发射率εv与地表背景发射率εs数据来源于光谱数据库；表示像元内组分间多次散射而形成的腔体效应参数。FVC为植被覆盖度，NDVIs和NDVIv分别为裸土和浓密植被的归一化植被指数NDVI值，ρ4与ρ5为经过大气辐射校正模型COST（Cosine of solar zenith angle）修正过的Landsat 8红波段Band 4与近红外波段Band 5的地表反射率。

步骤 4地表热辐射传输方程是遥感反演地表温度的基础，依据地表热辐射理论，建立地表温度LST的反演计算模型集，包括辐射传输方程单通道反演模型（SC1）[8]、单通道发射率校正模型（SC2）[9]、双通道劈窗算法模型（SW）[10]以及数据融合算法（DF），其具体推导过程及参数不再一一赘述，用户可具体参考每个算法对应的文献资料，主要公式如下：

1) SC1模型

式中：ψ1、ψ2与ψ3是大气函数，利用大气水汽含量ω近似获取，Lλ为星上辐射亮度值，bi,0、bi,1与bi,2为回归参数，K为常数，其他符号同上。

2) SC2模型

3) SW模型

式中：BT11为Landsat 8第11波段亮度温度，bn（n= 0,1,...7）为模型回归系数，ε和∆ε分别为TIRS两个通道的发射率均值与差值。

4) DF模型

依据空间相似性，每个MOD11A1数据产品像元对应不同地表类型Landsat 8像元的地表温度表达式为：

式中：PM为MOD11A1像元值，N为一个MOD11A1像元对应Landsat像元的个数，m为N个像元中所对应的高分辨率地表覆盖类型总数，PL,j为m个地物类别中j类别的像元值，fj标识不同地表类型之间高低分辨率之间的对应关系，ai,j(i= 0,1,...4) 为MODIS与Landsat混合像元法求解出的回归参数。

步骤 5遴选地表温度产品。以上4种反演模型产品以及第10与第11波段反演得到的亮温均值产品共计5种地表温度产品输出。为了与MODIS MOD11A1地温数据产品空间网格单元一一对应，每种模型反演产品分别进行坐标转换与重采样处理，然后逐像元进行统计比较，依据评价指标，选择最优反演结果纳入数据集。

2 数据样本描述

地表温度产品数据集命名格式为CPEC_LST_Xm_PPPRRR_YYYYMMDD.tif，其中X表示地表温度产品的分辨率，Landsat反演的地温产品为30 m，YYYY表示年份，MM表示月份，DD表示天，PPP是 Landsat WRS2全球参考系统的条带号，RRR 为其行号，例CPEC_LST_30m_152042_20170603.tif，中巴经济走廊区域中的数据如图3所示。目前数据起止时间为2013-2018年，并实时更新，目前数据共17 018景，数据容量3 822 GB，数据文件按照时间进行分类，数据格式为Geotiff文件，坐标系统为UTM坐标系统，数据类型为Single，温度单位为℃，可使用通过空间分析软件直接读取显示、融合以及分析（图4）。

图3 中巴经济走廊地表温度产品单幅数据示例（CPEC_LST_30m_149038_20170918.tif）

3 数据质量控制和评估

3.1 质量控制

数据集生产过程中对模型反演结果进行了逐景比较，每种模型反演的数据产品进行空间融合后与MODIS MOD11A1 V6地表温度产品逐像元比较，按照评价指标从中选择最优反演产品。本数据集的数据遴选指标包括纳什效率系数（NSE）、一致性指数（d）、Kling-Gupta系数（KGE）[11]、相关系数（R）与百分比偏差（Pbias）5种指标对数据产品进行质量控制。

3.2 数据评估

数据集质量评估采用整体评估和局部评估相结合的方法。整体评估采用随机抽样的方法，从全部数据产品中随机选取6组数据进行对比实验，数据样本分布如图5所示，每组数据产品的对比评价结果见表2。

图4 中巴经济走廊全区域地表温度产品数据拼接融合示例（CPEC_LST_1km_201807.tif）

图5 中巴经济走廊地表温度数据集随机抽样数据空间分布图（红色表示）

表2 中巴经济走廊地表温度数据集抽样数据整体统计特征

其次，为确保数据统计的随机性，开展了局部抽样评估，即再次从每组数据产品中随机选取子区域，进行对比分析，验证其空间相关性（图6），结果表明R2均在0.66以上，与整体趋势吻合。

图6 中巴经济走廊地表温度案例数据子区域随机采样统计分析

4 数据价值

本数据集在整体上与MODIS地表温度产品保持较高的空间相关性，由于Landsat 8较高的空间分辨率，可以反演出更加清晰的地表纹理特征，例如水体、居住区、植被等地物，表达了更加详尽的地表温度细节与特征（图7），可直接用于表征中巴经济走廊区域陆表温度时空变化，开展探索其时空分异特征与规律研究，为研究冻融灾害、气候变化、灾害预测、热岛效应、工程安全等科学问题提供基础数据资料，服务于中巴两国的线路建设、生态保护、资源调查、工程建设以及防灾减灾等社会经济发展需求。

图7 数据产品局部对照图（a为区域遥感影像，b为MODIS地表温度产品，c为本数据集产品）

5 数据使用方法和建议

本数据集可通过国家特殊环境、特殊功能观测研究台站共享服务平台（http://www.crensed.ac.cn）申请下载。由于数据量比较大也可联系数据作者，访问中国科学院超级计算兰州分中心的计算与数据平台，进行在线数据下载、共享、分析、处理及应用。数据文件均为GeoTIFF格式，可通过GIS相关的软件如ENVI、GRASS、ArcGIS等，直接进行查看与应用。

由于地表温度的遥感反演受地表发射率、大气透过率以及下垫面状况的影响，目前针对地表温度的遥感反演产品适宜于基于相对温度的科学应用，地表温度的真实值需要依据每个像元的发射率、地物分类以及大气条状况对本数据产品进行修正。

致 谢

感谢国家特殊环境、特殊功能观测研究台站共享服务平台给与的项目支持。感谢 USGS EarthExplorer网站提供免费的Landsat 8与MODIS数据产品。