GEE:哨兵数据下载,将不同景按同一天日期合成

03
六月
2021

输入:

在这里插入图片描述

输出:

在这里插入图片描述
源代码:https://code.earthengine.google.com/dc1c75bdbeb215d4671db79862acbffb?noload=true

如图所示,输入的影像集合里面有四幅影像,其中三幅是S2卫星同一天(5月7号)拍摄的,将同一天拍摄的s2数据按中值合成为一副。7月21号的合成为一副。最终输出两幅影像。


//#    哨兵数据下载,将不同景按同一天日期合成 
// date: 2021-06-03
// author: Ruan xiaohao  

var xcqbuf = table

function NDVI(img) {

var nir = img.select("B8");

var red = img.select("B4");

var ndvi = nir.subtract(red).divide(nir.add(red));

return ndvi;
}

var s2img = sen2.filterBounds(xcqbuf)

.filterDate('2020-12-01','2020-12-31')

.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 10))

.map(NDVI).sort('system:index')

print(s2img)  
            
Map.addLayer(s2img.first())

// 影像集合导出方法
function exportImageCollection(imgCol) {

  var indexList = imgCol.reduceColumns(ee.Reducer.toList(), ["system:index"])
                        .get("list");
                        
  // print(indexList)
  indexList.evaluate(function(indexs) {
  
    for (var i=0; i<indexs.length; i++) {
      
      var image = s2img.filter(ee.Filter.stringStartsWith('system:index', indexs[i].slice(0,31))).median()
      
      if(i>0)
      
      {
      var a = indexs[i].slice(0,31)
      
      var b = indexs[i-1].slice(0,31)}
      
      print(a)
      
      print(b)
      
      if (i>0&a==b)
      
        {continue;}

      Export.image.toDrive({
      
        image: image.clip(xcqbuf),
        
        description: indexs[i],
        
        fileNamePrefix: indexs[i],
        
        folder: 'shl',
        
        region: xcqbuf,
        
        scale: 10,
        
        // crs: "EPSG:32652",
        
        maxPixels: 1e13
      });
    }
  });
}

exportImageCollection(s2img);

Map.addLayer(xcqbuf,{},'roi_Boundary');

Map.centerObject(xcqbuf, 7);
TAG

网友评论

共有访客发表了评论
请登录后再发布评论,和谐社会,请文明发言,谢谢合作! 立即登录 注册会员