Google Earth Engine合并Landsat集报错：B5波段匹配失败求助

问题根源分析 #

你遇到的Image.select: Pattern 'B5' did not match any bands错误，核心原因有三个：

1. 重复定义函数：两次命名maskL8sr函数，后定义的会覆盖前一个，导致Landsat 5的掩膜逻辑完全失效
2. 波段编号不兼容：Landsat 5和8的波段编号对应关系不同（L5 B5是近红外，L8 B5是短波红外），直接合并后用统一的波段选择会出现匹配错误
3. 空影像未过滤：L5与L8的过渡阶段（如2012-2013年）可能没有覆盖研究区的有效影像，生成的空影像在计算EVI时找不到指定波段

修复步骤 #

1. 分离掩膜函数，分别处理L5和L8 #

给L5和L8的掩膜函数起不同名称，确保各自的逻辑生效：

// Landsat 5 掩膜函数
function maskL5sr(image) {
  var cloudShadowBitMask = 1 << 3;
  var cloudsBitMask = 1 << 5;
  var qa = image.select('QA_PIXEL');
  var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0)
      .and(qa.bitwiseAnd(cloudsBitMask).eq(0));
  return image.updateMask(mask).divide(10000)
      .copyProperties(image, ["system:time_start"]);
}

// Landsat 8 掩膜函数
function maskL8sr(image) {
  var cloudShadowBitMask = 1 << 3;
  var cloudsBitMask = 1 << 5;
  var qa = image.select('QA_PIXEL');
  var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0)
      .and(qa.bitwiseAnd(cloudsBitMask).eq(0));
  return image.updateMask(mask).divide(10000)
      .copyProperties(image, ["system:time_start"]);
}

2. 统一波段命名，对齐光学波段 #

将L5和L8的波段重命名为通用名称（nir/red/blue等），避免编号混淆：

// 处理L5，重命名波段
var collectionL5 = ee.ImageCollection('LANDSAT/LT05/C02/T1_TOA')
  .filterBounds(studyArea)
  .map(maskL5sr)
  .select(
    ['B1', 'B2', 'B3', 'B4', 'B5', 'B7'],
    ['blue', 'green', 'red', 'nir', 'swir1', 'swir2']
  );

// 处理L8，重命名波段
var collectionL8 = ee.ImageCollection('LANDSAT/LC08/C02/T1_TOA')
  .filterBounds(studyArea)
  .map(maskL8sr)
  .select(
    ['B2', 'B3', 'B4', 'B5', 'B6', 'B7'],
    ['blue', 'green', 'red', 'nir', 'swir1', 'swir2']
  );

3. 过滤空影像，避免无效合成 #

在生成年度合成前，先判断该年份是否有有效影像，没有则返回null或跳过：

var filterCollection = stepList.map(function(year){
  var startDate = ee.Date.fromYMD(year,5,1);
  var endDate = ee.Date.fromYMD(year,9,15);
  var yearCollection = imageCollection.filterDate(startDate, endDate);
  // 判断是否有影像，无则返回null
  return ee.Algorithms.If(
    yearCollection.size().gt(0),
    yearCollection.median().set('system:time_start', startDate),
    null
  );
});

// 过滤掉null值，得到有效合成影像集
var yearlyComposites = ee.ImageCollection(filterCollection).filter(ee.Filter.notNull());

4. 优化EVI计算，适配统一波段 #

基于重命名后的通用波段计算EVI，避免依赖原始波段编号：

function evi(img){
  var eviImg = img.expression(
    '2.5 * ((NIR - RED) / (NIR + 6 * RED - 7.5 * BLUE + 1))', {
      'NIR': img.select('nir'),
      'RED': img.select('red'),
      'BLUE': img.select('blue')
    }).rename('EVI');
  return img.addBands(eviImg);
}

完整修复代码 #

// 假设studyArea已定义，替换为你的研究区范围
var studyArea = ee.Geometry.Rectangle([xxx, xxx, xxx, xxx]);

// Landsat 5 掩膜函数
function maskL5sr(image) {
  var cloudShadowBitMask = 1 << 3;
  var cloudsBitMask = 1 << 5;
  var qa = image.select('QA_PIXEL');
  var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0)
      .and(qa.bitwiseAnd(cloudsBitMask).eq(0));
  return image.updateMask(mask).divide(10000)
      .copyProperties(image, ["system:time_start"]);
}

// Landsat 8 掩膜函数
function maskL8sr(image) {
  var cloudShadowBitMask = 1 << 3;
  var cloudsBitMask = 1 << 5;
  var qa = image.select('QA_PIXEL');
  var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0)
      .and(qa.bitwiseAnd(cloudsBitMask).eq(0));
  return image.updateMask(mask).divide(10000)
      .copyProperties(image, ["system:time_start"]);
}

// 处理L5并统一波段命名
var collectionL5 = ee.ImageCollection('LANDSAT/LT05/C02/T1_TOA')
  .filterBounds(studyArea)
  .map(maskL5sr)
  .select(
    ['B1', 'B2', 'B3', 'B4', 'B5', 'B7'],
    ['blue', 'green', 'red', 'nir', 'swir1', 'swir2']
  );

// 处理L8并统一波段命名
var collectionL8 = ee.ImageCollection('LANDSAT/LC08/C02/T1_TOA')
  .filterBounds(studyArea)
  .map(maskL8sr)
  .select(
    ['B2', 'B3', 'B4', 'B5', 'B6', 'B7'],
    ['blue', 'green', 'red', 'nir', 'swir1', 'swir2']
  );

// 合并影像集
var imageCollection = collectionL5.merge(collectionL8);

// 生成年份序列
var stepList = ee.List.sequence(1984,2021);

// 生成年度合成，过滤空影像
var filterCollection = stepList.map(function(year){
  var startDate = ee.Date.fromYMD(year,5,1);
  var endDate = ee.Date.fromYMD(year,9,15);
  var yearCollection = imageCollection.filterDate(startDate, endDate);
  return ee.Algorithms.If(
    yearCollection.size().gt(0),
    yearCollection.median().set('system:time_start', startDate),
    null
  );
});

var yearlyComposites = ee.ImageCollection(filterCollection).filter(ee.Filter.notNull());
print(yearlyComposites, 'Masked and Filtered Composites');

// 计算EVI
function evi(img){
  var eviImg = img.expression(
    '2.5 * ((NIR - RED) / (NIR + 6 * RED - 7.5 * BLUE + 1))', {
      'NIR': img.select('nir'),
      'RED': img.select('red'),
      'BLUE': img.select('blue')
    }).rename('EVI');
  return img.addBands(eviImg);
}

yearlyComposites = yearlyComposites.map(evi);
print(yearlyComposites, 'With EVI as Band');

// 提取EVI波段集合
var eviCollection = yearlyComposites.select('EVI');

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