Google Earth Engine（GEE）利用森林变化数据集分析

Google Earth Engine（GEE）利用森林变化数据集分析

今天来分享一下如何在GEE中利用森林变化数据集进行分析，所使用的数据集为UMD/hansen/global_forest_change_2020_v1_8 GEE调用代码如下：

//以山西省为研究区var roi = ee.FeatureCollection("users/lilei655123/shanxi");var GFC = ee.Image("UMD/hansen/global_forest_change_2020_v1_8")Map.centerObject(roi,5)var treeCover = GFC.select(['treecover2000']).clip(roi);var loss = GFC.select(['loss']).clip(roi);var gain = GFC.select(['gain']).clip(roi);// 绿色显示森林覆盖范围Map.addLayer(treeCover,{palette: ['000000', '00FF00'], max: 100}, 'Forest Cover')// 红色代表损失Map.addLayer(loss.updateMask(loss),{palette: ['FF0000']}, 'Loss')// 蓝色代表增益Map.addLayer(gain.updateMask(gain),{palette: ['0000FF']}, 'Gain')// 统计每年损失的面积var lossAreaImage = loss.multiply(ee.Image.pixelArea())var lossYear = GFC.select(['lossyear'])var lossByYear = lossAreaImage.addBands(lossYear).reduceRegion({reducer: ee.Reducer.sum().group({groupField: 1}),geometry: roi,scale: 30,maxPixels: 1e13})print(lossByYear)var statsFormatted = ee.List(lossByYear.get('groups')).map(function(year) {var d = ee.Dictionary(year);return [ee.Number(d.get('group')).format("20%02d"), d.get('sum')]})var statsDictionary = ee.Dictionary(statsFormatted.flatten())print(statsDictionary)var chart = ui.Chart.array.values({array: statsDictionary.values(),axis: 0,xLabels: statsDictionary.keys()}).setChartType('ColumnChart').setOptions({title: 'Yearly Forest Loss',hAxis: {title: 'Year', format: '####'},vAxis: {title: 'Area (square meters)'},legend: { position: "none" },lineWidth: 1,pointSize: 3})print(chart)

分析结果

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