中国区域大气–地形–BRDF全要素地表反射率产品

中国区域大气–地形–BRDF全要素地表反射率产品
作者：张兆明　何国金　彭燕　龙腾飞　冷宛春
2020年4月30日

本作品收录于《中国科学数据》

张兆明, 何国金, 彭燕, 等. 中国区域大气–地形–BRDF全要素地表反射率产品[J/OL]. 中国科学数据, 2020. (2020-04-27). DOI: 10.11922/csdata.2020.0020.zh.

摘要&关键词

摘要：在Landsat 5高精度正射影像的基础上，挑选高质量、无云的影像数据，经过大气–地形–BRDF耦合校正处理，消除大气散射和吸收、地形起伏以及地表二向反射（Bidirectional Reflectance Distribution Function，BRDF）对地表反射率反演的影响，生成了2009年中国区域的全要素地表反射率产品。数据空间分辨率为30 m，包含完整的数据文件、质量文件以及元数据文件。作为即得即用（RTU）产品的一种，大气–地形–BRDF全要素地表反射率产品可应用于地表覆盖变化调查、遥感地表参数定量反演、全球变化研究等领域。

关键词：地表反射率；定量反演；陆地卫星；大气–地形–BRDF耦合校正

Abstract & Keywords

Abstract: In order to eliminate the influences of atmospheric scattering and absorption, terrain fluctuation, and solar-viewing geometry, the coupled atmospheric-topographic-BRDF correction was performed to Landsat 5 ortho-rectified images covering the whole area of China in 2009. In this way, a new land surface reflectance product was produced and released. The spatial resolution of this product is 30 meter. The file folder contains data file, quality attribute file, and metadata file. This new product is valuable for land cover change monitoring, high precision land surface parameter retrieval, global change studies, etc.

Keywords: land surface reflectance; quantitative retrieval; Landsat satellite; atmospheric-terrain-BRDF coupled correction

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

数据库（集）名称	中国区域大气–地形–BRDF全要素地表反射率产品
数据作者	张兆明、何国金、彭燕、龙腾飞、冷宛春
数据通信作者	何国金（hegj@aircas.ac.cn）
数据时间范围	2009年
地理区域	中国
空间分辨率	30 m
数据量	349 GB
数据格式	*.TIF
数据服务系统网址	ftp://bigrs-info.com/public/China-LSR-2009/ http://www.sciencedb.cn/dataSet/handle/975
基金项目	中国科学院A类先导专项（XDA19090300）
数据库（集）组成	数据集共包括540个文件夹，每个文件夹包含数据文件（分波段存储，TIF格式）、质量文件（QA）以及元数据。

Dataset Profile

Title	An atmospheric-terrain-BRDF coupled corrected surface reflectance product over China
Data corresponding author	He Guojin (hegj@aircas.ac.cn)
Data authors	Zhang Zhaoming , He Guojin, Peng Yan, Long Tengfei, Leng Wanchun
Time range	2009
Geographical scope	China
Spatial resolution	30 m
Data volume	349 GB
Data format	*.TIF
Data service system	< ftp://bigrs-info.com/public/China-LSR-2009/> <http://www.sciencedb.cn/dataSet/handle/975>
Source of funding	Strategic Priority Research Program of the Chinese Academy of Sciences (XDA19090300).
Dataset composition	The dataset consists of 540 file folders in total. Each file folder contains data file in TIFF format, Quality Attribute file, and metadata file in XML format.

引 言

对地观测领域已经步入大数据时代，每天接收的卫星遥感数据呈爆炸式增长。通常，卫星遥感数据以原始影像数据产品（Digital Number值，即DN值）的形式提供给用户。DN值没有明确的物理意义，不同时间、不同地点、不同传感器的DN值不具有可比性，遥感数据用户难以直接使用。

为充分挖掘海量卫星数据的应用价值，迫切需要在原始影像数据的基础上进行辐射定量化处理，研发用户能直接使用的“即得即用”（Ready To Use，RTU)产品^[1]。地表反射率（Land Surface Reflectance，LSR）是最基本的遥感RTU产品之一，是开展遥感土地利用分类、长时序地表覆盖变化检测等应用的基础。

为了从原始卫星影像数据得到精确的地表反射率，需要综合考虑传感器辐射定标、大气散射和吸收、地形起伏以及地表二向反射（Bidirectional Reflectance Distribution Function，BRDF）的影响。目前已有的地表反射率产品通常只考虑了其中两项或者三项因素的影响^[2]，为了得到更精确的地表反射率RTU产品，需要研发大气–地形–BRDF耦合校正技术，从而生成大气–地形–BRDF全要素地表反射率产品。

1 数据采集和处理方法

1.1 数据采集方法

选择中国遥感卫星地面站接收的Landsat 5卫星数据，覆盖中国全境，数据获取时间为2009年，尽可能选择无云、高质量的影像数据，所有数据均做过高精度的正射校正处理。

1.2 数据处理方法

在对卫星原始影像各波段数据做高精度辐射定标的基础上，针对卫星传感器参数特点，利用6S辐射传输模型逐像元计算的方式模拟大气散射和吸收影响。6S模型需要的输入参数利用NCEP（National Centers for Environmental Prediction，美国国家环境预报中心）再分析资料或者与卫星过境时间同步的MODIS（Moderate Resolution Imaging Spectroradiometer，中分辨率成像光谱仪）大气参数产品（大气水蒸汽含量、气溶胶光学厚度、臭氧含量等），具体处理过程参见文献[3]。地形校正利用改进C地形校正算法，该算法具有较好的地形影响去除效果^[3]。BRDF校正利用C因子法^[4]，将地表反射率调整到同一“太阳-地表-传感器”观测几何下，从而减小太阳角度及卫星观测角度变化对地表反射率反演的影响。通过上述处理，最终实现传感器辐射定标、大气、地形和BRDF全要素耦合校正，得到地表反射率。

假设太阳天顶角为z，太阳方位角为φ，倾斜坡面的坡度和坡向分别为S和A，坡面太阳入射角为i，改进C地形校正模型计算过程如下：

T\(=\frac{\mathit{cos}z+c{h}_{0}^{-1}}{\mathit{cos}i+c{h}_{0}^{-1}h}\) （1）

图片 （2）

图片 （3）

图片 （4）

其中，T为地形校正系数，c=b/m, b和m可以通过在图像上自动选取样点，对地形校正前地表反射率和cosi之间的关系进行回归分析得到。

在实际计算时，利用NDVI（Normalized Difference Vegetation Index，归一化差值植被指数）区分植被和非植被地表，以NDVI=0.4为阈值，把影像分为植被区和非植被区，同时将坡度分级，分别统计c系数。

BRDF校正模型表达式如下：

图片 （5）

ρ 为对应于波长λ、观测角度Ω（观测天顶角和方位角）和太阳照射角度Ωˊ（太阳天顶角和方位角）的光谱反射率。Kvol(Ω, Ωˊ)和Kgeo(Ω, Ωˊ)分别为体散射（Ross Thick kernel）和几何光学散射核（Li Sparse Reciprocal Kernel），这两项只与太阳-传感器观测几何(Ω, Ωˊ)有关，是太阳天顶角、观测天顶角及相对方位角的三角函数。fiso 为各向同性散射系数，等于太阳天顶照射、传感器天顶观测时的地表反射率；fvol 和fgeo 为权重系数，分别表示体散射和几何光学散射所占的比例，后三个参数和波长相关。具体校正过程参见文献[5]。

2 数据样本描述

数据产品按照影像行列号（path/row）和获取时间存放在不同的文件夹中，文件夹中共包含9个文件，其中数据文件为TIFF格式，分波段存储，共6个文件，文件名末尾为B1-B5以及B7。数据质量文件为TIFF格式，以QA（Quality Attribute）标识。元数据文件为txt格式。数据文件的空间分辨率为30 m，为了降低存储空间，将地表反射率反演结果（浮点型）乘以10000变成16位整型，背景填充值为-9999。例如文件夹名为L5-TM-119-036-20090113-LSR-BRDF-TC，表示轨道号为119-036，获取时间为2009年1月13日的Landsat5全要素地表反射率产品，L5-TM-119-036-20090113-LSR-BRDF-TC-B1.TIF表示第一波段的数据文件，L5-TM-119-036-20090113-LSR-CLOUD-QA.TIF表示云质量文件，L5-TM-119-036-20090113-PIXEL-QA.TIF表示像元质量文件，L5-TM-119-036-20090113-MTL.txt表示元数据文件。

3 数据质量控制和评估

本节以具体例子通过对比分析来展示大气-地形-BRDF全要素地表反射率校正的效果。

通过图1对比可以看出，经过地形校正，地形起伏对地表反射率的影响得到很大程度的抑制，山体阴阳坡同种地物的地表反射率近似相等，地表反射率图像变得更加“平坦”和真实。

图片（a）

图片（b）

图片（c）

图片（d）

图片（e）

图片（f）

图1 地形校正效果对比a, c和e为校正前，b, d和f为校正后，a-b为整景影像对比，c-f为局部放大对比，Landsat 5卫星影像，RGB：543, path/row:120/42, 获取日期：2010-12-09。

为了检验BRDF校正的效果，选取地表反射率不随季节变化的地表（大片均质的裸地等，简称不变地表），选择冬季（1月）和夏季（7月）的卫星影像，由于冬季和夏季太阳高度角和卫星观测角的变化，这些不变地表的地表反射率随季节出现一定的波动（图2a），经过BRDF校正，减小了太阳高度角及卫星观测角变化对地表反射率反演的影响，冬季和夏季影像上这些不变地表的地表反射率变得更为接近。

图片

图2 BRDF校正效果对比a为校正前，b为校正后，path/row:123/32, 横轴为1999年1月14日影像，纵轴为1999年7月25日影像，地表反射率缩放系数为0.0001。

4 数据价值

中国全境大气–地形–BRDF全要素地表反射率产品可为地表覆盖变化调查、遥感地表参数反演、全球变化研究、生态环境监测等遥感应用提供高精度的基础数据产品。

参考文献

1. He G J, Zhang Z M, Jiao W L, et al. Generation of ready to use (RTU) products over China based on Landsat series data[J]. Big Earth Data, 2018, 2 (1): 56-64. DOI: 10.1080/20964471.2018.1433370.
2. Zhang Z M, He G J, Zhang X M, et al. A coupled atmospheric and topographic correction algorithm for remotely sensed satellite imagery over mountainous terrain[J]. GIScience & Remote Sensing, 2017, 54(5): 1-17.
3. FRANTZ D, RODER A, STELLMES M, et al. An Operational Radiometric Landsat Preprocessing Framework for Large-Area Time Series Applications[J]. IEEE Transactions on Geoscience & Remote Sensing, 2016, 54(7): 3928-3943.
4. ROY D P, ZHANG H K, JU J, et al. A general method to normalize Landsat reflectance data to nadir BRDF adjusted reflectance[J]. Remote Sensing of Environment, 2016, 176: 255-271.

数据引用格式

张兆明, 何国金, 彭燕, 等. 中国区域大气–地形–BRDF全要素地表反射率产品[DB/OL]. Science Data Bank, 2020. (2020-04-14). DOI: 10.11922/sciencedb.975.