GEE中Sentinel-2影像中位数合成与导出

给学生的授课教程，一点点记录和发布。

目的：筛选2022年10月至2023年5月清澜港红树林保护区Sentinel-2无云像，进行中位数合成，并导出。便于后续在eCognition或ENVI中分析、建模用。

下图为：shp导入 和 限定导出矩形边界

/**
 * Function to mask clouds using the Sentinel-2 QA band
 * @param {ee.Image} image Sentinel-2 image
 * @return {ee.Image} cloud masked Sentinel-2 image
 */
function maskS2clouds(image) {
  var qa = image.select('QA60');

  // Bits 10 and 11 are clouds and cirrus, respectively.
  var cloudBitMask = 1 << 10;
  var cirrusBitMask = 1 << 11;

  // Both flags should be set to zero, indicating clear conditions.
  var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
      .and(qa.bitwiseAnd(cirrusBitMask).eq(0));

  return image.updateMask(mask).divide(10000);
}

var dataset = ee.ImageCollection('COPERNICUS/S2_SR_HARMONIZED')
                  .filterBounds(HN) //利用海南岛边界限定筛选范围
                  .filterDate('2022-10-01', '2023-05-30') 
                  // Pre-filter to get less cloudy granules.
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE',20))
                  .map(maskS2clouds);

//计算影像集中位数，合成为image
var dataset_median = dataset.median();
print(dataset_median);

//提取所需波段
var dataset_mul = dataset_median.select('B2', 'B3', 'B4', 'B5', 'B6', 'B7','B8', 'B8A','B11', 'B12');  
print(dataset_mul);

var visualization = {
  min: 0.0,
  max: 0.3,
  bands: ['B4', 'B3', 'B2'],
};

Map.setCenter(110.7886, 19.62, 12);
Map.addLayer(dataset.median(), visualization, 'RGB');

Export.image.toDrive({
  image: dataset_mul,
  description: "2022_DZGS2_MUL_Median",
  fileNamePrefix: "2022_DZGS2_MUL_Median",
  region: QLG, //限定导出范围，需要自己设置
  scale: 10,   //scale:10意思是分辨率10米
  crs: "EPSG:4326", //Geodetic CRS: WGS 84, EPSG:32649 is WGS84 UTM Zone 49N
  folder: "GEE",
  maxPixels: 1e13
});