目标:
提取山西省不透水面数据并批量下载,使用的数据为清华大学宫鹏老师的不透水面数据集“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
如果对你有帮助的话记得给小编点个赞!