学习日记01.1 GEE下载Sentinel-2A数据(优化版)

其实本想在原文章上直接修改,但考虑到可能有不同的需求,还是决定另起一篇做单独介绍。

基本内容和前提条件不再赘述,可在之前的文章中进行查看。

本次优化主要体现于数据导入和出图部分

话不多说,直接上代码

Map.centerObject(table,7);
//这里可调整加载时影像的显示大小
Map.addLayer(table)
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).divide(10000);
}
var dataset = ee.ImageCollection('COPERNICUS/S2_SR')
                  .filterBounds(table)
                  .filterDate('2022-07-05', '2022-07-07')
                  // Pre-filter to get less cloudy granules.
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 10))
                  .map(maskS2clouds);
                  
var rgbVis = {
  min: 0.0,
  max: 0.3,
  bands: ['B4', 'B3', 'B2'],
};


Map.addLayer(dataset.median(), rgbVis, 'RGB');                 
var exportdataset =  ee.ImageCollection('COPERNICUS/S2_SR')
                  .filterBounds(table)
                  .filterDate('2022-07-05', '2022-07-07')
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 10))
                  .map(maskS2clouds)
                  .select(['B4', 'B3', 'B2']);
//镶嵌裁剪函数,减少之后的工作量
var mosaic = exportdataset.mosaic().clip(table);
  
  Export.image.toDrive({
      image:mosaic,
      description:'2022',
      scale:10,
      maxPixels: 1e13,
      region:table,
      fileFormat: 'GeoTIFF',
      formatOptions: {
        cloudOptimized: true
      }
    });

运行代码时仅需将shp导入,然后修改你所需的时间和波段,点击运行即可,最后下载后的影像是已经过shp裁剪过的影像,无需后续再使用Arcmap等软件裁剪。

 

 

  • 3
    点赞
  • 24
    收藏
    觉得还不错? 一键收藏
  • 2
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值