【GEE】基于GEE可视化和下载Sentinel2 L2A数据（去云、镶嵌、裁剪、筛选波段）

RS迷途小书童

已于 2024-03-30 12:08:24 修改

阅读量3.6k

点赞数 32

分类专栏： Google Earth Engine 文章标签： GEE 云计算 javascript

于 2024-02-01 12:16:14 首次发布

本文链接：https://blog.csdn.net/m0_56729804/article/details/135971019

版权

Google Earth Engine 专栏收录该内容

9 篇文章

订阅专栏

本文作者RS迷途小书童分享了如何在GoogleEarthEngine(GEE)中使用代码下载并处理Sentinel-2L2A数据，包括时间筛选、云量过滤、图像镶嵌和波段选择，以便高效下载所需数据。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

今天连续发了4篇关于Landsat8数据的下载代码，本来都不想再发GEE这个专栏的文章了，但是又想了想都快过年了，赶紧把手里的代码余货都分享出去吧，省的心里有压力。本篇文章主要分享了GEE可视化和下载Sentinel2 L2A数据。

原创作者：RS迷途小书童

博客地址：https://blog.csdn.net/m0_56729804?type=blog

代码主要包含了时间、云量筛选数据集，对数据进行中值合成、镶嵌、裁剪等操作，最后加入筛选波段的代码，仅导出我们需要的波段，以免数据过大。

PS：即使镶嵌后，下载时也会分4-5幅，这是因为GEE的限制问题，哨兵数据太大了一个市将近3个G，Landsat数据就没事。

一、代码部分

//作者：RS迷途小书童
//博客：https://blog.csdn.net/m0_56729804?type=blog

Map.centerObject(table,10);
//这里可调整加载时影像的显示大小
var style_set = {color:"red",fillColor:"00000000"};
Map.addLayer(table.style(style_set),{}, 'ROI')

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')//S2==L1C
                  .filterBounds(table)
                  .filterDate('2020-04-01', '2020-10-31')
                  // Pre-filter to get less cloudy granules.
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE',10))
                  .map(maskS2clouds)
                  .median()
                  .clip(table);
                  //.select(['B4', 'B3', 'B2']);
print(dataset);

var mergedImage = dataset.select('B1','B2','B3','B4','B5','B6','B7','B8','B8A','B9','B11','B12');
print(mergedImage);

var rgbVis = {
  min: 0.0,
  max: 0.3,
  bands: ['B4', 'B3', 'B2'],
};
 
Map.addLayer(mergedImage, rgbVis, 'RGB');                 
// //镶嵌裁剪函数，减少之后的工作量
// var mosaic = exportdataset.mosaic().clip(table);
  
  Export.image.toDrive({
      image:mergedImage,
      description:'2023',//Tasks名称
      folder: 'Sentinel2-L2A',//云盘文件夹
      scale:10,
      maxPixels: 1e13,
      region:table,
      fileFormat: 'GeoTIFF',
      crs: "EPSG:4326",//投影坐标系
      formatOptions: {
        cloudOptimized: true
      }
    });