Google Earth Engine(GEE)实例代码学习六——土地利用分类（Landcover Cleanup）

标题

利用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 = image2.reproject('EPSG:4326', null, SCALE);// 重投影为WGS84坐标系
var image4 = image3.focal_mode();//平滑图像
//采用不同窗口平滑图像目的是为了减小分类为很小的土地类型
var image5 = image4.focal_max(3).focal_min(5).focal_max(3);
// 将上面的影像再次投影
var image6 = image5.reproject('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类
// 将每步计算结果都进行显示
Map.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');

代码运行结果如下，可以看到全球土地类别

这个是进行平滑后的图像