Google Earth Engine（GEE）批量下载不透水面数据

目标：

提取山西省不透水面数据并批量下载，使用的数据为清华大学宫鹏老师的不透水面数据集“Tsinghua/FROM-GLC/GAIA/v10”

参考文献：Gong, P., Li, X., Wang, J., Bai, Y., Chen, B., Hu, T.,& Zhou, Y. (2020). Annual maps of global artificial impervious area (GAIA) between 1985 and 2018. Remote Sensing of Environment, 236, 111510.

数据集介绍：

数据集利用谷歌地球引擎平台上30米分辨率的landsat数据，绘制了1985年至2018年的年度GAIA。通过辅助数据集，包括夜间灯光数据和Sentinel-1合成孔径雷达数据，评估了1985年、1990年、1995年、2000年、2005年、2010年和2015年的GAIA数据，平均总体精度高于90%。

GAIA的数据集可以从

http://data.ess.tsinghua.edu.cn 免费下载

GEE实现代码：

var roi = ee.FeatureCollection("users/lilei655123/shanxi");
Map.centerObject(roi,7) 
var styling = {color:"red",fillColor:"00000000"};
Map.addLayer(roi.style(styling),{},"geometry")
//建立时间序列
var years = ee.List.sequence(1, 34);
var data_list=[]
years.evaluate(function(y) 
{
 data_list.push( y.map(function(y) {
    var year = y+1984
var image = ee.Image('Tsinghua/FROM-GLC/GAIA/v10')
.select("change_year_index").eq(y).clip(roi);
var visParam = {
 min: 0,
 max: 34,
 palette: '#56613B'
};
Map.addLayer(image, visParam,year+"year_image");
//每一年的像素个数 
 var imagecount = image.selfMask().reduceRegion({
  reducer: ee.Reducer.count(),               // 像素统计
  geometry: roi,
  scale: 30,
  maxPixels:1e13
  });
//逐年的面积计算
var unit_area = ee.Number(0.0009)
var year_area = imagecount.getNumber("change_year_index").multiply(unit_area);
// 研究区像素统计
 var studyarea = image.reduceRegion({
  reducer: ee.Reducer.count(),               
  geometry: roi,
  scale: 30,
  maxPixels:1e13
  });
 var study_area = studyarea.getNumber("change_year_index").multiply(unit_area);
//逐年的不透水面积所占比例
var areapercent = year_area.divide(study_area);
  var feature=ee.Feature(null,null);
  feature = feature.set('year', year);
  feature = feature.set('areapercent', areapercent);
  feature = feature.set('year_area', year_area);
//影像下载
  Export.image.toDrive({
    image: image,
    description: year+"year影像下载",
    folder: '影像',
    region:roi,
    scale: 30,
    maxPixels:1e13
  });
  return feature;
  }) )
 // 定义图表
 var featureCollection = ee.FeatureCollection(data_list[0]);
   Export.table.toDrive({
    collection: featureCollection,
    description: 'year_area',
    fileFormat: 'CSV'
  });
 var color = {
  high: 'ff0000',
  low: '0000ff'
};
var areaChart =
  ui.Chart.feature.byFeature(featureCollection, 'year', ['year_area'])
    .setChartType('ColumnChart')
    .setOptions({
      title: 'year-area',
      vAxis: {
        title: 'year_area'
      },
          lineWidth: 5,
          colors: ['de2d26'],
          curveType: 'function'
    });
print(areaChart);
var percent =
  ui.Chart.feature.byFeature(featureCollection, 'year', ['areapercent'])
    .setChartType('LineChart')
    .setOptions({
      title: 'year-areapercent',
      vAxis: {
        title: 'areapercent'
      },
          lineWidth: 5,
          colors: ['67000d'],
          curveType: 'function'
    });
print(percent);
});

结果展示：

统计结果：

1985年至2018年不透水面面积

1985年至2018年逐年的不透水面积所占比例

感谢关注，欢迎转发！

声明：仅供学习使用！

GEE demo：

https://code.earthengine.google.com/3b6eb772bea496fd65677c3122076c49

如果对你有帮助的话记得给小编点个赞！