天宫二号宽波段成像仪青藏高原专题数据集

李雷娟，任海根，李盛阳，覃帮勇*

1.中国科学院空间应用工程与技术中心，中国科学院太空应用重点实验室，北京100094

2.中国载人航天工程办公室，北京 100071

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

数据库（集）名称 天宫二号宽波段成像仪青藏高原专题数据集数据作者 李盛阳、张万峰、刘志文、覃帮勇、于海军、王博、邵雨阳、刘康、李雷娟、黑保琴、李轩、任海根数据通信作者 覃帮勇（qinby@csu.ac.cn）数据时间范围 2016年9月至2018年11月地理区域青藏高原地理位置：北纬 26°00′-39°47′，东经 73°19′-104°47′。主要由昆仑山脉、唐古拉山脉、喜马拉雅山脉等山脉，青海湖、色林错、纳木错等湖泊，长江、黄河、澜沧江、怒江等河流组成。空间分辨率 （1）可见近红外谱段（100 m）；（2）短波红外谱段（200 m）；（3）热红外谱段（400 m）。数据量 2.0 TB数据格式 影像数据为Geotiff格式（*.tif），影像缩略图为png格式，影像拇指图为jpg格式，辅助文件为xml格式数据服务系统网址 http://www.msadc.cn/t2/t2web/data/data_list.jsp http://www.sciencedb.cn/dataSet/handle/701

基金项目 国家科技部国家科技基础条件平台项目“国家基础科学数据共享服务平台”（DKA2018-12-02-23）。数据库（集）组成数据集由 4部分数据组成：（1）可见近红外谱段数据产品；（2）短波红外谱段数据产品；（3）热红外谱段数据产品；（4）专题产品：地表反射率产品、地表温度产品、地表亮温产品、归一化植被指数产品、归一化积雪指数产品、归一化水体指数产品。

引 言

青藏高原是世界屋脊、亚洲水塔和地球第三极，位于北纬 26°00′-39°47′，东经 73°19′-104°47′之间，东西长约2800公里，南北宽约300-1500公里，总面积约250万平方公里。境内高山大川密布，地势险峻多变，其平均海拔4000 m以上，是中国最大、世界海拔最高的高原。高原腹地年平均温度在0℃以下，最暖月平均温度也不足10℃，从南往北横跨亚热带湿润气候到寒带干旱高原气候等10个气候区，是北半球气候变化的启张器和调节器。其独特的地质构造、地形地貌和高原气候造就了特有的自然资源与生态环境系统，成为我国重要的战略资源储备基地和生态安全屏障[1-2]。

天宫二号空间实验室宽波段成像仪（以下简称宽波段成像仪）是国内首次在单台仪器上，实现多光谱大视场推扫成像功能的宽刈幅（约300公里）光学遥感器，包含3个谱段：可见近红外谱段（0.4 μm-1.0 μm，14 个通道）、短波红外谱段（1.0 μm-1.7 μm，2 个通道）和热红外谱段（8.0 μm-10.0 μm，2个通道）。3个谱段的空间分辨率分别为100 m、200 m和400 m[3-5]。自发射以来，宽波段成像仪通过对青藏高原的持续观测，获得了大量青藏高原数据产品，并生成了相应的专题产品，包括地表反射率产品、地表温度产品、地表亮温产品、归一化植被指数产品、归一化积雪指数产品、归一化水体指数产品，数据量总计超过2.0 TB。数据覆盖范围如图1所示。

图1 青藏高原专题数据集覆盖图

青藏高原专题数据集包含数据产品与专题产品，其中数据产品有可见近红外谱段数据产品、短波红外谱段数据产品和热红外谱段数据产品，专题产品有地表反射率产品、地表温度产品、归一化植被指数产品及归一化积雪指数产品等。本数据集可作为青藏高原科学研究与综合考察的参考数据，应用于气候变化、生物多样性、土地覆盖类型及湖泊水文监测变化等领域。

1 数据采集和处理方法

1.1 宽波段成像仪数据产品

宽波段成像仪采用推扫成像模式，可见近红外、短波红外与热红外3个谱段单独成像，可见近红外有3个视场，短波红外和热红外各有2个视场。数据下行后，首先进行解格式处理，得到各谱段格式化原始数据；再经过视场拼接和均一性校正，将多个视场组合成一幅宽刈幅影像；接着按照成像时间进行分景处理，将一轨图像分割成多个连续的标准景；标准景图像数据经过相对及绝对辐射校正后，生成入瞳处的辐亮度图像，再经过传感器校正和几何投影生成带有地理坐标的图像数据产品[3,5]。处理流程如图2。

图2 宽波段成像仪数据产品处理流程图

1.2 专题产品

宽波段成像仪专题产品包括地表反射率产品、地表温度产品、地表亮温产品、归一化植被指数产品、归一化积雪指数产品、归一化水体指数产品。在生成的标准数据产品的基础上，利用辐射传输模型或暗像元法进行大气校正，进而得到地表反射率和地表亮温产品，在此基础上进一步生成地表温度、归一化植被指数、归一化积雪指数、归一化水体指数专题产品。植被指数公式为：

这里NDVI为归一化植被指数，NIR为位于近红外波段的遥感通道所得到的反射率，RED为位于可见光波段的红光通道得到的反射率[6]，NIR与RED分别对应宽波段成像仪可见近红外谱段的通道5和通道9。

利用积雪在绿光波段反射率较高，在短波红外波段反射率较低的独特变化特征，构建的归一化积雪指数能有效识别积雪[7]。公式为：

这里NDSI为归一化积雪指数，Green为位于可见光波段的绿光通道得到的反射率，SWI为位于短波红外波段的遥感通道得到的反射率，Green与SWI分别对应宽波段成像仪可见近红外谱段通道11和短波红外谱段通道2。根据经验，一般NDSI≥0.4区域可以定义为冰雪覆盖区。

归一化水体指数是基于绿波段与近红外波段反射率的归一化比值指数[8]，其公式为：

这里NDWI为归一化水体指数，Green与NIR分别选取宽波段可见近红外谱段的通道11和通道5进行归一化差值处理，突出显示影像水体信息。

2 数据样本描述

宽波段成像仪数据产品以辐亮度值作为灰度值存储地面空间信息与光谱特征。本数据集采用对地观测领域标准数据格式，并配套有完整的元数据信息，方便数据读取和使用。

数据产品文件名以“飞行器_载荷名称_产品标识_数据类型_数据采集开始时间_数据采集结束时间_产品级别_分段号_分景号_产品生成时间_版本号.文件后缀”的格式命名[3]。

按照前后字段顺序，具体描述如表1。

表1 数据产品文件名字段描述表

可见近红外谱段数据产品是使用8、10、12三个通道合成的近似真彩色影像（图3），图片名称T2_MWI_VNI_IMG_20180121122946_20180121124748_L2_000_8_20181111143510_V211.tif。

图3 可见近红外谱段影像及xml文件

短波红外谱段数据产品通道1如图4所示，图片名称T2_MWI_SWI_IMG_20180121122946_201 80121124748_L2_000_8_20181111121833_V211.tif。

图4 短波红外谱段影像及xml文件

热红外谱段数据产品通道1如图5所示，图片名称T2_MWI_INF_IMG_20180121122946_20180 121124748_L2_000_8_20181111122009_V211.tif。

图5 热红外谱段影像及xml文件

地表反射率产品如图6、图7所示，图片名称分别为T2_MWI_VNI_LSR-RTM_20180121122946_20180121124748_L5_000_8_20181127091444_V211.tif、T2_MWI_SWI_LSR-RTM_20180121122946_20180121124748_L5_000_8_20181127092919_V211.tif。

图6 可见近红外谱段地表反射率产品（通道1）及xml文件

图7 短波红外谱段地表反射率产品（通道1）及xml文件

地表温度产品如图8所示，图片名称为T2_MWI_INF_LST_20180121122946_20180121124748_L5_000_8_20181127123217_V211.tif。

图8 地表温度产品（通道1）及xml文件

3 数据质量控制和评估

在数据处理的各个环节，对数据产品进行质量控制，以确保数据产品的质量和精度。针对传感器在轨辐射性能的变化，2018年8月在敦煌和青海湖开展了宽波段成像仪外场辐射定标试验，通过外场同步测量的方式，计算得到传感器各通道的绝对辐射定标系数，从而有效改善和提高了数据产品的辐射精度。宽波段成像仪的可见近红外和短波红外谱段的绝对辐射定标精度优于10%，热红外谱段的绝对辐射定标精度优于2K。

在数据产品质量评估方面，针对数据格式[9]、辐射质量、几何精度、云量等进行了评价，并根据评价结果进行了筛选，将满足质量要求的数据产品发布共享。

数据产品格式检查包括数据格式正确性、数据完整性及数据有效性检查。

（1）数据格式检查主要对数据的命名、数据内容及数据结构等的正确性进行检查；

（2）数据完整性检查主要对数据文件的完备性、数据内容的完整性等进行检查；

（3）数据有效性主要对数据及参数取值是否有效进行检查。

数据产品辐射质量要求地表反射率及相关的各类指数产品与 MODIS的地表反射率及专题产品的交叉验证平均误差在10%以内，地表亮温与MODIS地表亮温产品交叉验证平均误差在2K以内，地表温度产品的交叉验证平均误差在2.5K以内。数据产品的平均几何定位精度在10个像元以内。此外，为保障数据产品的可用性，在数据发布前对图像中的云量进行检测。对于云量较多的图像数据，由于可用性较低，在发布的数据集中进行了剔除。

4 数据价值

宽波段成像仪获取的中等分辨率多谱段数据产品，在青藏高原区域中大尺度地物目标监测、气候变化、生物多样性变化、土地覆盖类型及湖泊水文监测变化等方面具有广泛的应用前景。可见近红外谱段数据产品及地表反射率产品等专题产品可应用于气候变化监测[10-12]、土地覆盖类型调查[2]、冰川及积雪分布演变[13-16]、高原湖泊提取[17]、生态环境评价及突发性灾害监测等应用研究；短波红外谱段数据产品可用于土壤湿度[18-19]、农作物长势监测、云监测等应用研究；热红外谱段数据产品、地表温度产品可用于地表温度反演与分析、火灾监测、城市热岛效应等应用研究。