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);
});

4.结果展示：

显示结果

统计结果：

1985年至2018年不透水面面积
1985年至2018年逐年的不透水面积所占比例
可以一键下载

感谢关注，欢迎转发！

声明：仅供学习使用！如果对你有帮助的话记得给小编点个赞

**更多内容请关注微信公众号“生态遥感监测笔记”

Google Earth Engine（GEE）提取不透水面的方法相关推荐

使用google earth engine(GEE)提取2000年到2019年长江下游水体（河流、湖泊）数据
我最近想通过GEE直接统计长时间序列下长江下游流域的水体面积变化情况,如果通过传统做法很复杂,于是想到了使用GEE.下面是提取水体的效果图,时间是2000年的,你也可以设置显示2000-2019年中任 ...
基于google earth engine(GEE)下载研究区域影像
基于google earth engine(GEE)下载研究区域影像当研究需要Landsat数据时,我们可以通过USGS官网或者地理空间数据云平台下载.由于地理空间数据云目前无法下载到较新的数据,可 ...
Google Earth Engine(GEE)——User memory limit exceeded（2）
上一次我们已经知道如何去进行避免这种错误的发生,有关详细内容,如果单单只是解决这个问题我们用到的是limit 和 first,上一次的博客在这里: (207条消息) Google Earth Engi ...
Google Earth Engine(GEE) 01-中输入提示快捷键Ctrl+space无法使用的问题
Google Earth Engine(GEE) 01-中输入提示快捷键Ctrl+space无法使用的问题 GEE中 Ctrl+space组合键用于代码输入快捷提示,能够提高编码的准确度和速度,但是, ...
Google Earth Engine(GEE)批量下载代码（以 NDVI数据为例）
下载数据先准备工作(具体细节都能查到):科学上网谷歌邮箱,谷歌邮箱注册GEE账号. 一:导入需要下载边界shp文件. 标题找到自己的shp文件,导入除了sbx文件的所有文件. 导入成功命名ass ...
使用Google Earth Engine (GEE)实现MODIS数据批量下载
使用Google Earth Engine GEE实现MODIS数据批量下载前言下载数据代码批量执行run任务关注公众号,分享GIS知识.ArcGIS教程.SCI论文与科研日常等前言上图是 ...
Google Earth Engine(GEE)——可视化动态图
代码: var geometry = /* color: #d63000 *//* shown: false *//* displayProperties: [{"type": & ...
Google Earth Engine (GEE) ——卫星影像的监督分类（svm）
问题 GEE 提供哪些机器学习技术? 如何对卫星图像进行监督分类? 如何评估分类器的准确性? 如何手动创建自己的几何图形? 目标练习查找无云图像和使用手绘几何导入学习训练和应用分类算法所需的基本功 ...
关于google earth engine(GEE)的一些想法与大胆预测
我接触GEE有两年了,GEE留给我的印象是:无所不能. 不管是从庞大的数据量,还是包含遥感的各类算法:随机森林.SVM.CNN,都让人惊讶. 从GEE的云端操作来看,传统遥感需要几个月做出来的全国ND ...
google earth engine GEE批量 run下载插件
在Google earth engine中,批量导出文件经常遇到多个文件下载的情况,只能一个个点击RUN.这里介绍一款插件,Open Earth Engine extension .该插件的作者是Ma ...

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