Google Earth Engine（GEE）对指定地点Sentinel-2 Level1C数据的读取及云量处理

Google Earth Engine (GEE) 是由谷歌公司开发的众多应用之一。借助谷歌公司超强的服务器运算能力以及与NASA的合作关系，GEE平台将Landsat/Sentinel等可以公开获取的遥感图像数据存储在谷歌的磁盘阵列中，使得GEE用户可以方便的提取、调用和分析海量的遥感大数据资源。
 

本文基于GEE平台，实现对指定地点Sentinel-2 Level1C数据的读取及云量的处理，示例地点为湖南省长沙市。代码如下：

（JS）

/**
 * 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
 */
var cs = China.filterBounds(point).geometry();
print(cs)
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);
}

// Map the function over one year of data and take the median.
// Load Sentinel-2 TOA reflectance data.
var dataset = ee.ImageCollection('COPERNICUS/S2')
                  .filterDate('2018-01-01', '2019-12-31')
                  // Pre-filter to get less cloudy granules.
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
                  .map(maskS2clouds);

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

Map.addLayer(dataset.median().clip(cs), rgbVis, 'RGB');

注：①使用前，需要import我国的行政区划文件，其次需要选取例如湖南省长沙市境内的一点（point），最后import导入S2图                 像。

       ②关于我国行政区划文件，分享百度网盘：链接：https://pan.baidu.com/s/1vrT5AmRTX8KmMMbpGtTAcg 
                                                                           提取码：o5vt       
效果：

参考文章：

1.https://developers.google.com/earth-engine/datasets/catalog/COPERNICUS_S2