GEE学习总结（8）——依据土地利用分类结果进行数据掩膜提取

根据土地利用分类结果对GEE中的产品数据进行掩膜提取

• 获取土地利用分类结果数据的方法 ☞ GEE获取土地利用分类数据

以下通过MODIS土地利用分类数据的草地类型，提取MODIS的NDVI产品数据

var geometry = ee.FeatureCollection("users/13363961879/shanmeng");
var dataset = ee.ImageCollection('MODIS/006/MOD13Q1')
                  .filterDate('2000-06-01', '2000-09-01')
                  .select('NDVI')
                  .max().divide(10000)
                  .clip(geometry);

//调用土地利用数据对影像进行分类掩膜提取
var modisLandcover =  ee.ImageCollection('MODIS/006/MCD12Q1')
    .filterDate('2020-01-01', '2020-12-31')
    .first()
    .select('LC_Type1')
    .subtract(1)       //将波段排序从0开始排列
    .clip(geometry);   //切割出研究区范围
    
var mask=function(image,my_type){
  var mask=image.eq(my_type);
  var masked_image=image.updateMask(mask);
  return masked_image;
}
//调用mask函数，筛选出对应波段数据
//植被像元
var modismask = mask(modisLandcover,9)  
var grass = dataset.updateMask(modismask);
var modismask2 = mask(modisLandcover,11)  
var crop = dataset.updateMask(modismask2);
var modismask3 = mask(modisLandcover,3)  
var forest1 = dataset.updateMask(modismask3);
var modismask4 = mask(modisLandcover,4)  
var forest2 = dataset.updateMask(modismask4);
var modismask5 = mask(modisLandcover,8)  
var savannas = dataset.updateMask(modismask5);
var imageCollection=ee.ImageCollection([grass,crop,forest1,forest2,savannas])
var roi=imageCollection.mosaic()
//裸土像元
var modismask6 = mask(modisLandcover,15)  
var soil = dataset.updateMask(modismask6);

var igbpLandCoverVis = {
  min: 1.0,
  max: 17.0,
  palette: [
    '05450a', '086a10', '54a708', '78d203', '009900', 'c6b044', 'dcd159',
    'dade48', 'fbff13', 'b6ff05', '27ff87', 'c24f44', 'a5a5a5', 'ff6d4c',
    '69fff8', 'f9ffa4', '1c0dff'
  ],
};

var ndviVis = {
  min: 0.0,
  max: 1,
  palette: [
    'FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718', '74A901',
    '66A000', '529400', '3E8601', '207401', '056201', '004C00', '023B01',
    '012E01', '011D01', '011301'
  ],
};

//
Map.centerObject(geometry,7);
Map.addLayer(grass, ndviVis, 'NDVI'); 
Map.addLayer(modismask, igbpLandCoverVis, 'IGBP Land Cover')

• 草地覆盖范围
  
• 依据草地范围提取的NDVI数据