如何在现有S2代码中添加Landsat 8获取2012年以来湖泊表面积

完全可以通过加入Landsat系列传感器来将时间序列回溯至2012年，甚至更早。以下是具体实现思路和关键步骤：

1. 选择合适的Landsat数据集 #

• Landsat 8（2013年至今）：数据集ID为LANDSAT/LC08/C02/T1_L2（地表反射率产品，与Sentinel-2的产品类型一致）
• Landsat 7（1999年至今，2003年后存在扫描线误差）：数据集ID为LANDSAT/LE07/C02/T1_L2
• Landsat 5（1984年-2012年）：数据集ID为LANDSAT/LT05/C02/T1_L2，刚好覆盖你需要的2012年时段

2. 统一波段命名与水体指数计算 #

不同传感器的波段编号不同，需要先统一命名，确保后续分析逻辑一致：

• 编写通用波段重命名函数，将各传感器的绿、近红外、短波红外波段统一命名为green、nir、swir（匹配原有Sentinel-2代码中的水体指数计算逻辑，比如NDWI或MNDWI）
• 示例函数（以Landsat 8为例）：

function renameL8Bands(image) {
  return image.select(
    ['SR_B3', 'SR_B5', 'SR_B6'], // Landsat8的绿、近红外、短波红外波段
    ['green', 'nir', 'swir']
  ).copyProperties(image, ['system:time_start']);
}

• 对Landsat7、Landsat5编写类似的重命名函数，对应正确的波段编号即可

3. 统一云掩膜处理 #

各传感器的云检测依赖不同的QA波段，需编写对应的掩膜函数：

• Landsat 8/7/5的C02产品均包含QA_PIXEL波段，可通过位运算提取云、云阴影、雪等干扰像素：

function maskLandsatClouds(image) {
  var qa = image.select('QA_PIXEL');
  // 云（bit 3）、云阴影（bit 4）、雪（bit 5）掩码
  var mask = qa.bitwiseAnd(1 << 3).eq(0)
              .and(qa.bitwiseAnd(1 << 4).eq(0))
              .and(qa.bitwiseAnd(1 << 5).eq(0));
  return image.updateMask(mask);
}

4. 整合多传感器影像集合 #

将处理后的Sentinel-2和各Landsat影像集合合并，按时间排序：

// 处理Landsat8
var l8Collection = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2')
  .filterDate('2013-01-01', '2024-01-01')
  .filterBounds(yourPointCollection)
  .map(maskLandsatClouds)
  .map(renameL8Bands);

// 处理Landsat5（覆盖2012年）
var l5Collection = ee.ImageCollection('LANDSAT/LT05/C02/T1_L2')
  .filterDate('2012-01-01', '2013-01-01')
  .filterBounds(yourPointCollection)
  .map(maskLandsatClouds)
  .map(renameL5Bands); // 需提前定义Landsat5的波段重命名函数

// 合并原有Sentinel-2集合与Landsat集合
var combinedCollection = s2Collection.merge(l8Collection).merge(l5Collection)
  .sort('system:time_start');

5. 适配原有分析流程 #

由于已经统一了波段命名和基础处理，你原有代码中针对Sentinel-2的湖泊面积计算逻辑（如缓冲区裁剪、水体指数阈值提取、面积统计）可以直接复用在合并后的集合上，无需大幅修改。

注意事项：

• Landsat 7在2003年5月之后存在扫描线缺失问题，若需覆盖2003-2012年，可结合Landsat5和Landsat7的数据
• 不同传感器的空间分辨率不同（Sentinel-2为10/20m，Landsat为30m），统计面积时需注意分辨率差异带来的误差，可统一重采样至相同分辨率再分析

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