Google Earth Engine——基于新的Landsat SR数据集去云处理

根据GEE官方公告，明年原来的Landsat/LT05/C01/T1_SR和Landsat/LC08/C01/T1_SR数据集将停止更新，并提供了新的地表反射率数据，就是LANDSAT/LT05/C02/T1_L2和LANDSAT/LC08/C02/T1_L2。新的数据集直接提供LST产品，方便获取地表温度，同时也给出官方代码去云处理。

参考：USGS Landsat 5 Level 2, Collection 2, Tier 1 (google.com)

去云代码如下：

/****************************************L5 去云****************************************/
/***************************************************************************************/
// 使用Landsat5 Collection 2,Level 2 QA_PIXEL波段（CFMask）去云
function maskL5sr(image) {
// Bit 0 - fill
// Bit 1 - Dilated Cloud
// Bit 2 - Unused
// Bit 3 - Cloud
// Bit 4 - Cloud Shadow
var qaMask = image.select('QA_PIXEL').bitwiseAnd(parseInt('11111', 2)).eq(0);
var saturationMask = image.select('QA_RADSAT').eq(0);

// 将缩放系数应用于适当的波段
var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);
var thermalBand = image.select('ST_B6').multiply(0.00341802).add(149.0);

// 用缩放后的波段替换，并应用云掩膜
return image.addBands(opticalBands, null, true)
.addBands(thermalBand, null, true)
.updateMask(qaMask)
.updateMask(saturationMask);
}

//处理一年的数据
var collection = ee.ImageCollection('LANDSAT/LT05/C02/T1_L2')
.filterDate('2010-01-01', '2011-01-01')
.map(maskL5sr);
//中值合成
var composite = collection.median();

// 地图展示
Map.addLayer(composite, {bands: ['SR_B3', 'SR_B2', 'SR_B1'], min: 0, max: 0.3});
/***************************************************************************************/

/****************************************L8 去云****************************************/
/***************************************************************************************/
// 使用Landsat8 Collection 2,Level 2 QA_PIXEL波段（CFMask）去云 

function maskL8sr(image) {
  // Bit 0 - Fill
  // Bit 1 - Dilated Cloud
  // Bit 2 - Cirrus
  // Bit 3 - Cloud
  // Bit 4 - Cloud Shadow
  var qaMask = image.select('QA_PIXEL').bitwiseAnd(parseInt('11111', 2)).eq(0);
  var saturationMask = image.select('QA_RADSAT').eq(0);

  // 用缩放后的波段替换，并应用云掩膜
  var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);
  var thermalBands = image.select('ST_B.*').multiply(0.00341802).add(149.0);

  // 用缩放后的波段替换，并应用云掩膜
  return image.addBands(opticalBands, null, true)
      .addBands(thermalBands, null, true)
      .updateMask(qaMask)
      .updateMask(saturationMask);
}

// 处理一年的数据
var collection = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2')
                     .filterDate('2020-01-01', '2021-01-01')
                     .map(maskL8sr);
//中值合成
var composite = collection.median();

// 地图展示
Map.addLayer(composite, {bands: ['SR_B4', 'SR_B3', 'SR_B2'], min: 0, max: 0.3});
/***************************************************************************************/