Google Earth Engine（GEE）提取不透水面的方法

最新推荐文章于 2024-02-27 09:00:00 发布

生态遥感监测笔记

最新推荐文章于 2024-02-27 09:00:00 发布

阅读量3k

点赞数 2

分类专栏： GEE 文章标签：云计算 javascript 大数据 chrome

本文链接：https://blog.csdn.net/guaicaicexu/article/details/125871736

版权

GEE 专栏收录该内容

22 篇文章 81 订阅

订阅专栏

1.目标：

提取山西省不透水面，使用的数据为清华大学宫鹏老师的不透水面数据集
“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.

2.数据集介绍：

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

3.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: '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);
});