Sentinel数据S2_SR_HARMONIZED连续云掩膜+中位数合成

在GEE中实现时,发现简单的QA60是无法去云的,最近S2地表反射率数据集又进行了更新,原有的属性集也进行了变化,现在的SR数据集名称是"S2_SR_HARMONIZED"。那么:

要想得到研究区无云的图像,可以参考执行以下几个步骤,

  1. 遥感图像中筛选云占比不超过10%的图像

    // 加载Sentinel-2 L2A图像

    var s2 = ee.ImageCollection('COPERNICUS/S2_SR_HARMONIZED')

    .filterDate('2021-09-01', '2021-10-31')

    .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE',10))

    .filterBounds(basin_shp);

  2. QA60云掩膜

    // 第一次云掩膜:使用QA60

    function maskS2clouds(image) {

    var qa = image.select('QA60');

    var cloudBitMask = 1 << 10;

    var cirrusBitMask = 1 << 11;

    var mask = qa.bitwiseAnd(cloudBitMask).eq(0)

    .and(qa.bitwiseAnd(cirrusBitMask).eq(0));

    return image.updateMask(mask);

    }

    var s2_masked = s2.map(maskS2clouds);

  3. SCL在最新的S2_SR_HARMONIZED产品中是存在的,判断云和云阴影

    //SCL Cloud mask

    function maskS2cloudsSCL(image) {

    var scl = image.select('SCL');

    // 去除云(8)和云阴影(9)

    var mask = scl.neq(8).and(scl.neq(9));

    return image.updateMask(mask);

    }

    var s2_masked = s2_masked.map(maskS2cloudsSCL);

  4. 继续用像元云概率进行判断,要求每个像元是云的概率小于10%

    // MSK_CLDPRB Mask(Sen2Cor 云概率),可直接用阈值掩膜,云概率低于 10% 保留:

    function maskCloudProb(image) {

    var prob = image.select('MSK_CLDPRB');

    return image.updateMask(prob.lt(10)); // 云概率低于 20% 保留

    }

    var s2_masked = s2_masked.map(maskCloudProb);

    最后中位数合成,避免反射率高值影响图像像元整体性。

    var s2_max = s2_masked.median();

相关推荐
非凡ghost2 天前
Hard Disk Sentinel(固态硬盘监控) 多语便携版
windows·sentinel·软件需求
柳贯一(逆流河版)4 天前
Sentinel 深度解析:限流与熔断降级的微服务稳定性保障实践
微服务·架构·sentinel
耳东哇5 天前
sentinel docker gateway k8s 集群 主从
docker·gateway·sentinel
遥感之家6 天前
AWS下载sentinel-2原始影像
云计算·sentinel·aws
月夕·花晨6 天前
Gateway-断言
java·开发语言·分布式·spring cloud·微服务·nacos·sentinel
齐 飞8 天前
Spring Cloud Alibaba快速入门-Sentinel熔断规则
spring boot·spring cloud·sentinel
齐 飞9 天前
Spring Cloud Alibaba快速入门-Sentinel流量控制(FlowRule)
spring cloud·微服务·sentinel
Takumilove10 天前
Spring Boot 接入 Redis Sentinel:自动主从切换与读写分离实战(修复单机多实例与 Sentinel 配置坑)
spring boot·redis·sentinel
波波烤鸭19 天前
Redis 高可用实战源码解析(Sentinel + Cluster 整合应用)
数据库·redis·sentinel
hzzzzzo021 天前
微服务保护全攻略:从雪崩到 Sentinel 实战
数据库·微服务·sentinel