GEE合并不同MODIS影像集计算归一化差值波段识别错误如何解决

问题根本原因 #

你错误使用了ee.ImageCollection.merge()方法，这个方法的作用是将两个影像集合的影像条目堆叠拼接：比如b01集合有100张仅含sur_refl_b01的影像，b11集合有100张仅含sur_refl_b11的影像，merge之后得到的是200张单波段影像的集合，前100张只有b01、后100张只有b11，自然会出现找不到对应波段、集合波段不一致的报错。你需要的是按规则匹配两个集合的对应影像，将两个波段合并到同一张影像中。

修复方案 #

可以直接用GEE内置的combine方法合并两个等长、顺序一致的集合的波段，也可以用时间匹配的join方法保证匹配准确性，修复后完整代码如下：

// 定义研究区
var roi = ee.Geometry.Polygon(
        [[[-122.0, 45.0],
          [-122.1, 32.0],
          [-108.9, 32.0],
          [-108.9, 45.0]]]);

// 加载MOD09GA的b01波段集合
var modis_b01_500m = ee.ImageCollection('MODIS/006/MOD09GA')
                   .select(['sur_refl_b01'])
                   .filterDate('2001-02-28','2001-12-31')
                   .filterBounds(roi)
                   .map(function(image){return image.clip(roi)});

// 加载MODOCGA的b11波段集合，重采样为500m
var modis_b11_1km = ee.ImageCollection('MODIS/006/MODOCGA')
                  .select(['sur_refl_b11'])
                  .filterDate('2001-02-28','2001-12-31')
                  .filterBounds(roi)
                  .map(function(image){return image.clip(roi)});

var modis_b11_resampled = modis_b11_1km.map(function(image) {
  return image.reproject(image.projection(), null, 500)
  .copyProperties(image, ['system:time_start']); // 保留时间属性方便后续校验
})

// 直接合并两个集合的对应波段，要求两个集合影像数量一致、顺序完全对应
var modis_combined = modis_b01_500m.combine(modis_b11_resampled);

// 计算归一化差值
var cci = modis_combined.map(function(image){
  return image.normalizedDifference(['sur_refl_b11','sur_refl_b01']).rename('cci');
})

// 结果校验与可视化
print('合并后集合', modis_combined);
print('CCI结果集合', cci);
Map.addLayer(modis_combined, {bands:['sur_refl_b01'], min:0, max:10000}, 'b01波段');
Map.addLayer(modis_combined, {bands:['sur_refl_b11'], min:0, max:10000}, 'b11波段');
Map.addLayer(cci, {min:-1, max:1, palette:['blue','white','red']}, 'CCI结果');

注意事项 #

• 如果两个集合的影像数量、顺序不能保证完全对应，可以改用ee.Join.saveFirst按system:time_start属性匹配对应影像再合并波段，避免出现波段错位
• 重采样操作注意保留原影像的关键属性，方便后续计算与校验

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