GoogleEarthEngine(GEE)实例代码学习六——⼟地利⽤分类(Landcov。。
。
标题
利⽤MODSI全球⼟地覆盖产品显⽰⼟地利⽤分类
⾸先⼤致介绍⼀下MODIS⼟地覆盖产品MCD12Q1
MCD12Q1 V6产品按六种不同的分类⽅案按年间隔(2001-2016年)提供全球⼟地覆盖类型。 它是使⽤MODIS Terra和Aqua反射数据的监督分类得出的。 然后,对监督分类结果进⾏额外的后处理,这些后处理结合了先验知识和辅助信息以进⼀步完善分类。
分为1-17表⽰不同的类别
LC_Type1表⽰IGBP分类系统
LC_Type表⽰不同的分类系统,可通过查看该产品⾃⾏查阅
实例代码如下:
var SCALE=500;//在投影时强制⽐例为500⽶/像素
var image1 = ee.Image('MODIS/051/MCD12Q1/2001_01_01');// 加载MODIS⼟地覆盖产品MCD12Q1
var image2 = image1.select(['Land_Cover_Type_1']);// 选择IGBP分类系统
var image3 = project('EPSG:4326',null,SCALE);// 重投影为WGS84坐标系
var image4 = image3.focal_mode();//平滑图像
//采⽤不同窗⼝平滑图像⽬的是为了减⼩分类为很⼩的⼟地类型
var image5 = image4.focal_max(3).focal_min(5).focal_max(3);
// 将上⾯的影像再次投影
var image6 = project('EPSG:4326',null,SCALE);
// 设定⼀个⼟地分类调⾊板
var PALETTE=[
'aec3d4',// water
'152106','225129','369b47','30eb5b','387242',// forest
'6a2325','c3aa69','b76031','d9903d','91af40',// shrub, grass, savannah
'111149',// wetlands
'cdb33b',// croplands
'cc0013',// urban
'33280d',// crop mosaic
'd7cdcc',// snow and ice
'f7e084',// barren
'6f6f6f'// tundra
].join(',');
var vis_params ={min:0, max:17, palette:PALETTE};//类别有17类
// 将每步计算结果都进⾏显⽰
landMap.setCenter(116,39.9,6);
Map.addLayer(image2, vis_params,'IGBP classification');
Map.addLayer(image3, vis_params,'Reprojected');
Map.addLayer(image4, vis_params,'Mode');
Map.addLayer(image5, vis_params,'Smooth');
Map.addLayer(image6, vis_params,'Smooth');
代码运⾏结果如下,可以看到全球⼟地类别
这个是进⾏平滑后的图像
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论