全球30米分辨率火烧迹地产品

全球30米分辨率火烧迹地产品
作者：张兆明　唐朝　何国金　龙腾飞　魏明月
2020年4月30日

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

张兆明, 唐朝，何国金, 等. 全球30米分辨率火烧迹地产品[J/OL]. 中国科学数据, 2020, 4(4). (2020-04-27). DOI: 10.11922/csdata.2020.0019.zh.

摘要&关键词

摘要：火烧迹地是全球变化和碳循环等研究领域的一个重要参量。在全球高精度样本库基础上，基于Landsat 8时序卫星数据和火烧迹地敏感光谱参量，利用机器学习算法生产并发布全球30米分辨率火烧迹地产品。全球30米分辨率火烧迹地产品能够有效反映面积较小的燃烧斑块，同时在火烧迹地位置确定和面积量算上具有优势，可应用于全球火灾监测和评估、碳排放计算、生态环境保护等领域。

关键词：Landsat 8；火烧迹地；全球；机器学习

Abstract & Keywords

Abstract: Burned area (BA) is an important parameter in the field of global change and carbon cycle. Firstly, high quality global sample database was constructed, then Landsat 8 time series satellite data and several sensitive spectral parameters of burned area were inputted to the machine learning algorithm, and 30 meter resolution global burned area product was produced and released. The 30 meter resolution global BA product can effectively detect small burned patches, and has advantages in determining fire location and area of the burned patches. It can be applied to global fire monitoring and disaster assessment, carbon emission calculation, ecological environment protection and other fields.

Keywords: Landsat 8; burned area; global; machine learning algorithm

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

数据库（集）名称	全球30米分辨率火烧迹地产品
数据作者	张兆明、唐朝、何国金、龙腾飞、魏明月
数据通信作者	何国金（hegj@aircas.ac.cn）
数据时间范围	2015年1月至2015年12月
地理区域	60°S–80°N，180°W–180°E
空间分辨率	30 m
数据量	5.39 GB
数据格式	*.TIF
数据服务系统网址	ftp://bigrs-info.com/GABAM/burned area/2015/ http://www.sciencedb.cn/dataSet/handle/976
基金项目	国家重点研发计划课题（2016YFA0600302）
数据库（集）组成	数据集包含经纬度10度乘以10度分块的数据文件，文件格式为TIF格式，共504个文件。

Dataset Profile

Title	Generation of 30 meter resolution global burned area products based on Landsat 8 images
Data corresponding author	Guojin He (hegj@aircas.ac.cn)
Data authors	Zhaoming Zhang, Chao Tang, Guojin He, Tengfei Long, Mingyue Wei
Time range	From January 2015 to December 2015
Geographical scope	60°S–80°N，180°W–180°E
Spatial resolution	30 m
Data volume	5.39 GB
Data format	*.TIF
Data service system	<ftp://bigrs-info.com/GABAM/burned area/2015/> <http://www.sciencedb.cn/dataSet/handle/976>
Source of funding	National Key Research and Development Program of China（2016YFA0600302）
Dataset composition	The data set consists of 10*10 degree files, 504 files in total. The file format is TIF.

引 言

火烧迹地是全球及区域碳循环和气候变化等研究所需的重要参量之一。美国国家航空航天局（NASA）、欧洲空间局（ESA）和国际地圈生物圈计划（IGBP）等都曾发布过各种基于不同卫星遥感数据获取的全球尺度火烧迹地信息产品，空间分辨率从250米至0.25度不等。美国陆地卫星Landsat系列数据是目前存档时间最长、使用最广泛的中高空间分辨率卫星数据，为全球火烧迹地信息高精度提取提供了理想的数据源，利用Landsat系列卫星数据研发高空间分辨率和精度的大尺度火烧迹地信息产品具有迫切的应用需求，也是国际上的一个热点研究课题。美国地质调查局于2017年发布了覆盖美国本土的Landsat 30米分辨率火烧迹地产品^[1]，也制定出研发全球30米分辨率火烧迹地产品的科学计划。然而，相比MODIS等中低分辨率卫星数据，Landsat系列卫星数据具有较少的光谱波段和更长的重返周期，现有的中低分辨率火烧迹地信息提取方法无法直接应用。因此，基于Landsat数据的全球尺度火烧迹地信息的高精度自动提取是一个具有挑战性的研究方向。

近年来，机器学习方法在遥感大数据分析和智能信息提取领域获得了成功的应用。通过利用全球火烧迹地样本库和基于生态系统分区的机器学习方法，实现了基于Landsat数据的全球火烧迹地信息高精度自动化提取^[1]，生产了全球30米空间分辨率火烧迹地遥感产品。

1 数据采集和处理方法

1.1 数据采集方法

主要利用Landsat 8卫星OLI传感器数据。按照该卫星的运行计划，可以覆盖60°S–80°N的所有陆地区域。制作年度全球火烧迹地产品，需要利用当年获取的全部Landsat 8 OLI数据，以尽可能减少漏分现象，同时需要利用前一年获取的全部Landsat 8 OLI数据，去除前一年发生且未恢复的火烧迹地，制作一期全球火烧迹地产品需要利用约40万景Landsat 8 OLI影像，数据量超过300TB，在产品生产时利用了Google Earth Engine（简称GEE）云平台上的Landsat 8 OLI地表反射率或者星上反射率^[2]。

1.2 数据处理方法

在全球高精度样本库基础上，基于Landsat 8时序卫星数据和火烧迹地敏感光谱参量〔波段反射率、归一化燃烧指数（Normalized Burned Ratio，NBR），归一化差值植被指数（Normalized Difference Vegetation Index，NDVI），全球环境监测指数（Global Environmental Monitoring Index，GEMI），中红外燃烧指数（Mid-Infrared Burn Index，MIRBI），火烧迹地指数（Burned Area Index，BAI），修正的土壤调节植被指数（Soil-Adjusted Vegetation Index，SAVI）和归一化差值湿度指数（Normalized Difference Moisture Index，NDMI）等〕，利用机器学习算法（随机森林模型）进行样本训练和学习，得到火烧迹地识别规则和疑似火烧迹地种子点。对疑似火烧迹地种子点进行一系列过滤和优化，得到确定的火烧迹地种子点。在种子点周围进行区域生长，生成最终的火烧迹地空间分布结果。总体技术流程如图1所示。

图片

图1 火烧迹地提取技术流程

2 数据样本描述

火烧迹地数据产品利用经纬度投影，分辨率为0.00025°（约30米），以10°×10°分块，GeoTiff格式存储，以分块图像左上角像素的经纬度作为文件名，像素值中以1代表火烧迹地像元，以0代表非火烧迹地像元。例如文件名为S20W100.TIF的分块文件代表左上角像素坐标为南纬20度、西经100度的火烧迹地数据文件。

3 数据质量控制和评估

利用随机分层抽样的方式对全球火烧迹地产品进行精度验证和评估。为确保精度验证的全面性和代表性，精度验证样区的选择兼顾不同地表覆盖类型和火行为特征。依据MODIS地表覆盖类型产品将全球概括为7个地表类型（Broadleaved Evergreen（常绿阔叶林），Broadleaved Deciduous（落叶阔叶林），Coniferous（针叶林），Mixed Forest（混交林），Shrub（灌丛），Rangeland（草地），和Agriculture（农田），），同时根据GFED 4（Global Fire Emission Database 4，第四版本的全球火烧排放数据库）2015年的火烧迹地密度数据将全球均匀划分为5个密度级，最终在全球范围内选择了80个验证样区，这些样区覆盖了全部的地表类型和火行为分区。

精度验证利用的数据源主要包括Landsat 8、GF 1、CBERS 4和美国MTBS（Monitoring Trends in Burn Severity，火烧强度趋势监测）火烧迹地产品。其中中国境内的验证样区使用GF 1数据，南美区域使用CBERS 4，美国区域使用MTBS和Landsat 8，全球其他区域使用Landsat8。对于Landsat8数据，验证样区的大小为185 km×185 km；对于CBERS4数据，验证样区的大小为120 km×120 km；对于GF1数据，验证样区的大小为100 km×100 km。

通过收集验证样区2015年整年的时序卫星数据，进行正射校正等处理。然后目视寻找2015年内新增的火烧迹地，确定火烧前后的卫星影像对。基于卫星影像对手动选择火烧迹地和非火烧迹地样本点，利用支持向量机分类器进行分类，最终得到2015年火烧迹地验证的参考数据。

在每个验证样区分别进行验证，然后得到全球火烧迹地精度验证结果，选用误分率（Commission Error）、漏分率（Omission Error）和整体精度（Overall Accuracy）3个指标来定量表征火烧迹地信息提取的精度，最终得到2015年全球火烧迹地产品的误分率、漏分率和整体精度分别为13.17%、30.13%和93.92%^[3]。

4 数据价值

火烧迹地是资源环境监测、全球变化和碳循环研究等领域的一个重要参数，全球30米分辨率火烧迹地产品能够有效反映面积较小的燃烧斑块，同时在火烧迹地位置确定和面积量算上具有优势，可应用于全球火灾监测和评估、碳排放计算、生态环境保护等领域。

参考文献

1. ^ ^{1.0 1.1} Hawbaker T J, Vanderhoof M K, Beal Y J, et al. Mapping burned areas using dense time-series of Landsat data[J]. Remote Sensing of Environment, 2017, 198: 504-522.
2. LONG T F, ZHANG Z M, HE G J, et al. 30 m Resolution Global Annual Burned Area Mapping Based on Landsat Images and Google Earth Engine[J]. Remote Sensing, 2019, 11: 489-519.
3. GORELICK N, HANCHER M, DIXON M, et al. Google earth engine: Planetary-scale geospatial analysis for everyone[J]. Remote Sensing of Environment, 2017, 202: 18-27.

数据引用格式

张兆明, 唐朝, 何国金, 等. 全球30米分辨率火烧迹地产品[DB/OL]. Science Data Bank, 2020. (2020-04-15). DOI: 10.11922/sciencedb.976.