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

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的数据集可以从

免费-

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+1984var 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

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

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。