基于Landsat数据的三亚市陆表温度反演数据集

谷艳春，胡蝶，张颖，张琳琳

1.福州大学福建省空间信息工程研究中心，福州 350000

2.海南省地球观测重点实验室，海南三亚 572029

3.中国科学院遥感与数字地球研究所，北京 100101

4.中国科学院大学资源与环境学院，北京 100049

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

引 言

陆表温度即陆地表面温度（Land Surface Temperature，LST）是环境监测、气候、水文、农业等陆表研究分析和应用模型创建中的重要输入参数，也是其他很多研究和应用的基础，如土壤水分监测、地热位置判别等均离不开陆表温度的测定[1]。反演陆表温度一直是遥感定量化应用的重要内容，在具体研究中选择最佳的陆表温度反演算法，是提高反演精度的首要因素。

目前，针对陆表温度反演已提出了诸多不同算法，如基于影像的反演算法[2]、单通道算法[3-4]、单窗算法[5-7]、劈窗算法[8-13]和多通道算法[14-16]。其中最常用的是由Becker等在1990年提出的劈窗算法，此算法经诸多改进，已广泛应用于NOAA系列卫星热红外遥感数据。但由于陆地表面组成物质的不均一性、复杂性以及大气影响，遥感影像反演陆表温度始终是一个重要而复杂的研究课题[17]。RSTAR辐射传输模型充分考虑到这些问题，弥补了长序列陆表温度反演有较多不足之处，如大气参数未考虑；大气传输物理意义不明确；只有在晴空无云且精确估算地表发射率的情况下精度高；对NOAA系列卫星热红外遥感数据效果最佳；一般的辐射传输方程方法参数难确定、温度反演过程繁杂。随着遥感卫星数量的逐年增加，遥感影像也随之愈加丰富，因此充分发挥遥感能够长时间序列的温度反演，可实现区域温度监测的优势。

本文利用Landsat 5和Landsat 8影像数据，基于RSTAR辐射传输模型计算大气廓线，采用覃志豪[18]等提出的简单易行且精度高的地表发射率估算方法，计算出三亚市2000年、2004年、2008年、2012年、2016年和2018年夏季(6-8月)共6景（无云或少云）陆表温度，以便从整体监测三亚市陆表温度变化趋势。

本文的研究区三亚市（18°09′34″-18°37′27″N，108°56′30″-109°48′28″E）位于海南岛的最南端，东邻陵水县，西接乐东县，北毗保亭县，南临南海。气候属于热带海洋季风气候，年均气温25.7°C。

1 数据采集和处理方法

本数据集构建过程包括数据下载、计算大气上下行辐射和大气透过率参数、温度反演代码实现、裁剪和输出几个步骤，如图1。

其中，数据处理过程中使用到的IDL是一种数据分析和图像化应用程序及编程语言，有灵活的数据输入输出、快速实现可视化、支持外部语言接口等优势，是遥感常用编程方式处理影像的最佳选择，本数据集用到的是IDL 8.5版本。

1.1 数据来源与数据概括

Landsat系列影像数据以相对较高的空间分辨率优势（如Landsat 8热红外波段的空间分辨率为100 m），有助于提高陆表温度反演结果的精度。1984年3月1日，USGS成功发射的Landsat 5卫星有一个热红外波段Band 6（10.40 -12.5μm）；2013年2月11日，NASA成功发射的Landsat 8卫星有两个热红外波段Band 10（10.6-11.2μm）和Band 11（11.2-12.5μm），因Band 11运行出现问题，故本文利用Band 10单波段反演地表温度。本文所用Landsat系列影像数据均可在地理空间数据云网站(http://www.gscloud.cn/search)获取。

图1 数据集构建及处理流程

1.2 数据处理

1.2.1 大气剖面参数计算

大气剖面参数计算主要是依据大气辐射传输模式计算大气上行辐射、下行辐射和大气透过率参数，大气传输模式为RSTAR辐射传输模式。本文根据具体影像的过境时间、中心经纬度、所属季度（夏季或冬季）和传感器类型，直接在大气校正参数计算器网站（https://atmcorr.gsfc.nasa.gov/）计算得到每幅影像的大气上行辐射、下行辐射和大气透过率参数。

1.2.2 可见光波段处理

根据Chander等研究[19]，比较来自不同传感器图像时，应使用大气表观反射率（Top of Atmosphere Reflectance，TOA reflectance），具体做法如下：首先，消除不同太阳天顶角余弦值产生的误差；其次，利用TOA reflectance补偿光谱带差引起的太阳辐照度的不同值；最后，用TOA reflectance校正由不同数据采集时间导致的日地距离的变化。该过程的关键代码（图2）是在IDL 8.5中通过代码完成的。

图2 可见光波段处理关键代码

1.2.3 地表发射率计算

基于覃志豪[18]等提出的地表发射率反演算法，对水体掩膜后，将地表分为自然（植被+裸土）和城镇（植被+建筑表面）两类。算法主要有三步：首先，计算植被覆盖度；其次，城市表面可以看作是由不同比例的植被和建筑物组成，自然表面可以看作是由不同比例的植被和裸露土壤组成；最后，生成整个研究区影像的地表发射率图。该过程的关键代码如图3所示。

图3 计算地表发射率关键代码

1.2.4 陆表温度反演

本小节主要分为三步：首先，将大气剖面参数计算得到的大气上行辐射、下行辐射和大气透过率参数输入图4所示的代码里；其次，根据可见光波段处理、地表发射率计算的数据计算得到反演的温度图；最后选择温度反演结果的存储位置。该过程的关键代码如图4所示。

图4 陆表温度反演关键代码

2 数据样本描述

在多种陆表温度反演算法中选择辐射传输方程计算方法，基于 RSTAR辐射传输模型计算大气廓线，采用覃志豪[18]提出的地表发射率估算方法，计算出2000年、2004年、2008年、2012年、2016年和2018年共6景（无云或少云）三亚市陆表温度,部分结果如图5所示。

图5 2018年三亚市陆表温度反演结果

3 数据质量控制

本数据集通过以下几方面进行质量控制：

（1）数据源质量控制。在选取Landsat卫星数据时，以有数据、无云或者少云、无条带为选取标准。因为高质量的数据源是得到高精度结果的必要条件。

（2）温度反演结果质量控制。本文为验证陆表温度产品精度，依据以往三亚市温度数据统计分析反演得到的温度与公布的同地同时地表温度的相关性，其相关系数R2为0.91，并以此作为验证陆表温度产品数据的可靠依据。

4 数据价值

本研究提出了适合于城市热环境特征的地表温度反演技术，为遥感服务于城市建设提供技术支持。数据成果具有很强的实用性、普适性，如反演产品可为城市生态环境研究提供数据产品服务，无论区域大小均可适用。同时数据成果可为环保部门和城建部门开展城市热环境变化监测提供检测方式，具有很强的应用价值和实用性。

通过遥感手段反演城市热空间信息参量，充分发挥了遥感监测的周期性和时效性，便于在短时间快速分析判断城市生态环境时空演化特征。