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});

代码运行结果
在这里插入图片描述
发现检测出有云的位置像元值为空

  • 6
    点赞
  • 74
    收藏
    觉得还不错? 一键收藏
  • 10
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 10
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值