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

        今天连续发了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
      }
    });

二、结果展示

        我平时用GEE比较少,大量时间都是使用Python去处理数据,但是用的少也还是会一点的。令人失望的是网上有很多教程都是VIP文章或者时间周期太长,导致代码不能使用,对于没接触过代码的新手来说很不友好。我希望能在力所能及的范围内尽可能多地去分享一些GEE的基本操作,如果大家感兴趣也可以一起留言交流。

  • 22
    点赞
  • 32
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 8
    评论
要在Google Earth Engine (GEE) 中下载Sentinel-1数据,可以使用 `ee.ImageCollection` 和 `.filter()` 函数来筛选和获取所需的影像。以下是一个示例代码,展示如何从Sentinel-1 GRD产品中获取一幅影像的方法: ``` var collection = ee.ImageCollection("COPERNICUS/S1_GRD") .filter(ee.Filter.date('YYYY-MM-DD', 'YYYY-MM-DD')) // 根据日期筛选影像 .filterBounds(geometry) // 根据地理范围筛选影像 .select(['VV', 'VH']); // 选择VV和VH波段 var image = collection.first(); // 获取第一幅影像 // 可以使用Export函数将影像下载到Google Drive Export.image.toDrive({ image: image, description: 'Sentinel-1_Image', folder: 'Sentinel-1_Images', fileNamePrefix: 'Sentinel-1_Image', region: image.geometry(), scale: 10, // 设置分辨率 maxPixels: 1e12 // 设置最大像素数 }); ``` 上述代码中,使用了 `ee.ImageCollection` 函数获取了Sentinel-1 GRD产品的图像集合,然后根据日期和地理范围进行筛选,最后选择了VV和VH波段。通过 `Export.image.toDrive` 函数将所需影像下载到Google Drive。需要注意的是,代码中的`'YYYY-MM-DD'`需要替换为具体的日期,并且`geometry`需要替换为所需的地理范围。 请注意,为了使用GEE下载Sentinel-1数据,您需要在GEE环境中具有相关权限。具体操作细节可以参考官方文档和教程。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* [GEE Sentinel-1影像数据下载及处理](https://blog.csdn.net/XinemaChen/article/details/123024261)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *3* [ddl8s1ts:使用Landsat 8和Sentinel-1时间序列检测并绘制森林干扰图](https://download.csdn.net/download/weixin_42098892/15206452)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 8
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

RS迷途小书童

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值