GEE：掩膜掉不符合条件的像素

掩膜掉不符合条件的像素主要用两行代码

// 生成 DEM 掩膜
var DEM_mask = elevation.select('elevation').gt(4700)

//使用DEM掩膜
var img = img_max.updateMask(DEM_mask)

---

例一：按掩膜提取，去除NDVI数据中DEM小于4700m的像素。相当于用DEM大于4700米的区域做研究区掩膜。

源代码链接：https://code.earthengine.google.com/c38f4a4e50a6de445ad1d36460ff3485?noload=true

// 添加试验区矢量
var geo=ee.Geometry.Rectangle(90, 35, 95, 30)
print("geo",geo)
Map.setCenter(93, 32,6)
Map.addLayer(geo, {}, 'geo');

//DEM 数据下载
var DEM = ee.Image("NASA/ASTER_GED/AG100_003")
var elevation = DEM.select('elevation')
 .clip(geo)

// 生成 DEM 掩膜
var DEM_mask = elevation.select('elevation').gt(4700)

// 加载 来自MODIS的NDVI 数据
var collection = ee.ImageCollection("MODIS/006/MOD13Q1")
 .filterDate('2000-01-01', '2018-01-01')
 .select('NDVI');

//NDVI按最大值合成
var img_max=collection.select("NDVI").max();

//掩膜掉DEM小于4700m的像素
var img = img_max.updateMask(DEM_mask)

// 可视化参数
var args = {
  crs: 'EPSG:3857',  // Maps Mercator
  dimensions: '300',
  region: geo,
  min: -2000,
  max: 10000,
  palette: 'black, blanchedalmond, green, green',
  framesPerSecond: 12,
};

// 打印影像属性
print(img,"img");

// 加载影像
Map.addLayer(img, args, 'img');

---

例二：去除NDVI数据中NDVI小于1000的像素。相当于保留NDVI大于1000的区域做研究区掩膜。

源代码链接：https://code.earthengine.google.com/be4cc632fab0fef352a34b2b80e68e3e?noload=true

// 添加试验区矢量
var geo=ee.Geometry.Rectangle(90, 35, 95, 30)
print("geo",geo)
Map.setCenter(93, 32,6)
Map.addLayer(geo, {}, 'geo');

// 加载 来自MODIS的NDVI 数据
var collection = ee.ImageCollection("MODIS/006/MOD13Q1")
 .filterDate('2000-01-01', '2018-01-01')
 .select('NDVI');

//NDVI按最大值合成
var img_max=collection.select("NDVI").max();

// 生成 NDVI掩膜
var NDVI_mask = img_max.select('NDVI').gt(3000)

//掩膜掉DEM小于3000m的像素
var img = img_max.updateMask(NDVI_mask).clip(geo)

// 可视化参数
var args = {
  crs: 'EPSG:3857',  // Maps Mercator
  dimensions: '300',
  region: geo,
  min: -2000,
  max: 10000,
  palette: 'black, blanchedalmond, green, green',
  framesPerSecond: 12,
};

print(img,"img");
Map.addLayer(img, args, 'img');

例三：合并例一和例二的掩膜，共同作为掩膜提取NDVI

源代码链接：https://code.earthengine.google.com/f75f938e5a2d064e863fb7c37c802fb4?noload=true

// 添加试验区矢量
var geo=ee.Geometry.Rectangle(90, 35, 95, 30)
print("geo",geo)
Map.setCenter(93, 32,6)
Map.addLayer(geo, {}, 'geo');

//DEM 数据下载
var DEM = ee.Image("NASA/ASTER_GED/AG100_003")
var elevation = DEM.select('elevation')
 .clip(geo)

// 生成 DEM 掩膜
var DEM_mask = elevation.select('elevation').gt(4700)

// 加载 来自MODIS的NDVI 数据
var collection = ee.ImageCollection("MODIS/006/MOD13Q1")
 .filterDate('2000-01-01', '2018-01-01')
 .select('NDVI');

//NDVI按最大值合成
var img_max=collection.select("NDVI").max();

// 生成 NDVI掩膜
var NDVI_mask = img_max.select('NDVI').gt(3000)

//掩膜掉DEM小于4700m的像素
var img = img_max.updateMask(DEM_mask)

img = img.updateMask(NDVI_mask)

// 可视化参数
var args = {
  crs: 'EPSG:3857',  // Maps Mercator
  dimensions: '300',
  region: geo,
  min: -2000,
  max: 10000,
  palette: 'black, blanchedalmond, green, green',
  framesPerSecond: 12,
};

print(img,"img");
Map.addLayer(img, args, 'img');