影像去云处理
本文分享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});
代码运行结果
发现检测出有云的位置像元值为空