GEE：基于MCD64A1的GlobFire火灾斑块检测数据集

作者：CSDN @ _养乐多_

本文将介绍 Google Earth Engine（GEE）上存档的 JRC/GWIS/GlobFire/v2/FinalPerimeters 数据集。

---

---

一、背景介绍

1.1 MCD64A1数据集

MCD64A1是MODIS（Moderate Resolution Imaging Spectroradiometer）卫星遥感数据的一部分，专门用于监测全球火灾热点和火灾烧迹。该数据集提供了每天、每16天和每月的火灾热点信息，以及全球范围内的火灾烧迹。

1.2 GlobFire

GlobFire是一个项目，旨在使用卫星数据和遥感技术来监测全球火灾活动。该项目的目标之一是提供全球火灾事件的及时和高质量信息，以帮助研究人员、政府和公众更好地理解火灾的发生、影响和趋势。

1.3 参考文献

Artés, T., Oom, D., De Rigo, D., Durrant, T. H., Maianti, P., Libertà, G., & San-Miguel-Ayanz, J. (2019). A global wildfire dataset for the analysis of fire regimes and fire behaviour. Scientific data, 6(1), 1-11. doi:10.1038/s41597-019-0312-2

1.4 GEE上存档

ee.FeatureCollection("JRC/GWIS/GlobFire/v2/FinalPerimeters")

基于MODIS数据集MCD64A1的火灾边界。这些数据是通过一种算法计算得出的，该算法依赖于将烧毁区域的空时关系编码为图结构。

每个火灾事件都有一个唯一的编号来标识。

名称	类型	描述
area	Double	火灾面积，平方米
FinalDate	Int	最终火灾日期，自1970-01-01以来的毫秒数
Id	Int	火灾的数值标识符
InitialDate	Int	初始火灾日期，自1970-01-01以来的毫秒数

二、案例代码

2.1 查看所有数据ID

// 一系列表格的文件夹名称。
var folder = 'JRC/GWIS/GlobFire/v2/DailyPerimeters';

// 使用异步回调的ee.data.listAssets列出可用的表格。
function printAssetList(listAssetsOutput) {
  print('文件列表:', listAssetsOutput['assets']);
}
ee.data.listAssets(folder, {}, printAssetList);

结果如下图所示，可以看出来数据集中有2000年——2021年的数据。

2.2 可视化研究区火灾斑块的面积

代码链接：

https://code.earthengine.google.com/695f5479fcfacf6585a794cb40cafc01?noload=true

var roi = geometry;// 画一个多边形作为研究区边界

// 定义从可用表格列表中识别出的表格名称（表格 ID）。
var tableName = 'JRC/GWIS/GlobFire/v2/DailyPerimeters/2020';

var computeArea = function(f) {
  return f.set({'area': f.area()});
};

// 将选定的表格作为 FeatureCollection 导入。
var features = ee.FeatureCollection(tableName).filterBounds(roi).map(computeArea);

// 用于线性火灾面积渐变的可视化参数。
var visParams = {
  palette: ['f5ff64', 'b5ffb4', 'beeaff', 'ffc0e8', '8e8dff', 'adadad'],
  min: 0,
  max: 600000000,
  opacity: 0.8,
};

// 使用计算得到的火灾面积作为值属性，将火灾边界绘制到图像上。
var image = ee.Image().float().paint(features, 'area');

// 将图像显示到地图上（包括用于使用 Inspector 探索的要素）。
Map.addLayer(image, visParams, 'GlobFire 2020');
Map.addLayer(features, null, '用于检查的要素', false);
Map.setCenter(126.116, 44.452, 7);

筛选出来研究区的火灾斑块，如图，

2.3 可视化全球火灾斑块的面积

该代码将可视化每个火灾区域的面积。通过不同的颜色展示火灾斑块的面积。

代码如下，

// 定义从可用表格列表中识别出的表格名称（表格 ID）。
var tableName = 'JRC/GWIS/GlobFire/v2/DailyPerimeters/2020';

var computeArea = function(f) {
  return f.set({'area': f.area()});
};

// 将选定的表格作为 FeatureCollection 导入。
var features = ee.FeatureCollection(tableName).map(computeArea);

// 用于线性火灾面积渐变的可视化参数。
var visParams = {
  palette: ['f5ff64', 'b5ffb4', 'beeaff', 'ffc0e8', '8e8dff', 'adadad'],
  min: 0,
  max: 600000000,
  opacity: 0.8,
};

// 使用计算得到的火灾面积作为值属性，将火灾边界绘制到图像上。
var image = ee.Image().float().paint(features, 'area');

// 将图像显示到地图上（包括用于使用 Inspector 探索的要素）。
Map.addLayer(image, visParams, 'GlobFire 2020');
Map.addLayer(features, null, '用于检查的要素', false);
Map.setCenter(126.116, 44.452, 7);

结果如下图所示，