如何在Google Earth Engine中仅加载MODIS特定土地覆盖类型

嘿，我来帮你搞定这个筛选土地覆盖类别的问题！你之前尝试的.select('LC_Type1', [1,5])思路不对——select()方法是用来选择影像的波段，而不是筛选单个波段里的像素值。咱们用正确的方式来保留1-5类（各类森林）的像素：

核心思路：利用掩码筛选像素 #

我们可以通过创建布尔掩码，标记出属于1-5类的像素，然后用这个掩码过滤原始影像，把不符合条件的像素设为NoData。

完整实现代码 #

// 假设你已经定义了感兴趣区域roi（替换成你的实际区域）
var roi = ee.Geometry.Rectangle([100, 20, 110, 30]); 

// 加载MODIS土地覆盖数据集
var landcover = ee.Image("MODIS/061/MCD12Q1/2002_01_01")
                  .select('LC_Type1')
                  .clip(roi);

// 定义需要保留的土地覆盖类别（1-5类森林）
var desiredClasses = ee.List([1, 2, 3, 4, 5]);

// 创建掩码：筛选出属于目标类别的像素
var forestMask = landcover.inList(desiredClasses);

// 应用掩码，得到仅包含目标类别的土地覆盖影像
var landcoverFiltered = landcover.updateMask(forestMask);

// 自定义可视化参数（只保留1-5类对应的颜色）
var LCVisFiltered = {
  min: 1.0,
  max: 5.0,
  palette: [
    '05450a', '086a10', '54a708', '78d203', '009900' // 分别对应1-5类的颜色
  ],
};

// 添加到地图展示
Map.centerObject(roi, 8);
Map.addLayer(landcoverFiltered, LCVisFiltered, 'Filtered Forest Land-Cover');

代码解释 #

• landcover.inList(desiredClasses)：生成一个布尔影像，像素值为true表示该像素属于1-5类，false则不属于。
• updateMask(forestMask)：将布尔影像中false对应的像素设为NoData，最终影像就只显示你需要的森林类别了。

后续用于NPP提取的拓展 #

你提到要提取这些区域的NPP值，直接用这个掩码过滤NPP影像即可：

// 示例：加载MODIS NPP数据集（替换成你实际使用的NPP数据集）
var npp = ee.Image("MODIS/061/MOD17A3HGF/2002_01_01").select('Npp');

// 应用森林掩码，仅保留森林区域的NPP值
var nppForest = npp.updateMask(forestMask);

这样处理后，你就可以基于nppForest来绘制年度NPP时序图表啦～

内容的提问来源于stack exchange，提问作者frog11235