Google Earth Engine(GEE)实例代码学习三十四——利用Landsat质量评估波段(QA)去云(Landsat457 Surface Reflectance))

影像去云处理

本文分享Landsat4，5，7利用其从CFMASK算法生成的质量评估波段QA进行去云处理。
质量评估波段（QA）是从CFMASK算法生成的关于像元质量的波段单位为二进制。
主要操作是按位与运算bitwiseAnd（）筛选像元
其不同位数表示含义如下
Bit 0: Fill
Bit 1: Clear
Bit 2: Water
Bit 3: Cloud Shadow
Bit 4: Snow
Bit 5: Cloud
Bits 6-7: Cloud Confidence
0: None
1: Low
2: Medium
3: High

// 创建一个函数用来掩模云层以及云阴影
var cloudMaskL457 = function(image) {
  var qa = image.select('pixel_qa');
//云层表示为第五位，云层置信度为6-7位，云阴影为第三位
//选择出有云并且云层置信度为中等，以及有云阴影覆盖的像元。
  var cloud = qa.bitwiseAnd(1 << 5)
          .and(qa.bitwiseAnd(1 << 7))
          .or(qa.bitwiseAnd(1 << 3))
  // 移除边界像元
  var mask2 = image.mask().reduce(ee.Reducer.min());
  //将检测有关云像元置为0，掩模保留位置不为0的数据。
  return image.updateMask(cloud.not()).updateMask(mask2);
};

// M选取Landsat5地表反射率数据
var collection = ee.ImageCollection('LANDSAT/LT05/C01/T1_SR')
    .filterDate('2010-04-01', '2010-07-30')
//对影像数据集每张影像进行云掩膜处理，并求出中值
var composite = collection
    .map(cloudMaskL457)
    .median();

// 影像显示
Map.setCenter(-6.2622, 53.3473, 12);
Map.addLayer(composite, {bands: ['B3', 'B2', 'B1'], min: 0, max: 3000});

代码运行结果

发现检测出有云的位置像元值为空