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