1

My research area needs multiple landsat images(About 15 images), so I want to mosaic images where they are in the same path first (they have the same acquisition date). Howerver, I only know the mosaic operation for the specified date

ee.Image(collection.filterDate('2018-01-01','2018-01-02').mosaic())

So how do I get the date of each image in the image collection? Or is there any better iterative way to mosaic the same date images?

1 Answer 1

0

I typically iterate over the dates you need your mosaic for:

var startDate = '2022-05-01'
var endDate = '2022-06-01'

var collection = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2')
  
var numberOfDays = ee.Date(endDate).difference(ee.Date(startDate), 'days')
var dayOffsets = ee.List.sequence(0, numberOfDays.subtract(1))
var dailyMosaics = ee.ImageCollection(dayOffsets
  .map(function (dayOffset) {
    var date = ee.Date(startDate).advance(ee.Number(dayOffset), 'days')
    var images = collection.filterDate(date.getRange('day'))
    return ee.ImageCollection(images)
      .map(preProcess)
      .mosaic()
      .set('system:time_start', date.millis())
  })
)

Map.addLayer(dailyMosaics.first(), {bands: 'SR_B4,SR_B3,SR_B2', min: 0, max: 3000, gamma: 1.5})


function preProcess(image) {
  var qa = image.select('QA_PIXEL')
  var cloudShadow = bitwiseExtract(qa, 4)
  var snow = bitwiseExtract(qa, 5)
  var cloud = bitwiseExtract(qa, 6).not()
  return image
    .select('SR_B.*')
    .multiply(2.75e-05)
    .add(-0.2)
    .multiply(10000)
    .int16()
    .updateMask(cloudShadow.not().and(snow.not()).and(cloud.not()))
}

function bitwiseExtract(value, fromBit, toBit) {
  if (toBit === undefined)
    toBit = fromBit
  var maskSize = ee.Number(1).add(toBit).subtract(fromBit)
  var mask = ee.Number(1).leftShift(maskSize).subtract(1)
  return value.rightShift(fromBit).bitwiseAnd(mask)
}

https://code.earthengine.google.com/e08db00830e500decf7b55ded7642864

Alternatively, you could use a join:

var dates = ee.FeatureCollection(dayOffsets
  .map(function (dayOffset) {
      var date = ee.Date(startDate).advance(ee.Number(dayOffset), 'days')
        .advance(12, 'hours')
      return ee.Feature(null, {'system:time_start': date.millis()})
  })
)
var twelveHoursInMillis = 12 * 60 * 60 * 1000
var dailyMosaics = ee.ImageCollection(
  ee.Join.saveAll('images').apply({
    primary: dates,
    secondary: collection,
    condition: ee.Filter.maxDifference({
      difference: twelveHoursInMillis,
      leftField: 'system:time_start',
      rightField: 'system:time_start'
    })
  })
).map(function (feature) {
  return ee.ImageCollection(ee.List(feature.get('images')))
    .map(preProcess)
    .mosaic()
    .set('system:time_start', feature.get('system:time_start'))
})

https://code.earthengine.google.com/411502413e90e52f36e356461e738ffc

1
  • Thanks for your help, this is exactly what I want
    – yiboliu
    Jun 9 at 7:47

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.