Reducing image collection

Reducing an ImageCollection

To composite images in an ImageCollection, use imageCollection.reduce(). This will composite all the images in the collection to a single image representing, for example, the min, max, mean or standard deviation of the images. (See the Reducers section for more information about reducers). For example, to create a median value image from a collection:

Create an interactive map¶

In [ ]:

import ee
import geemap

In [ ]:

Map = geemap.Map()
Map

Compute a median image¶

In [ ]:

# Load a Landsat 8 collection for a single path-row.
collection = ee.ImageCollection('LANDSAT/LC08/C01/T1_TOA') \
    .filter(ee.Filter.eq('WRS_PATH', 44)) \
    .filter(ee.Filter.eq('WRS_ROW', 34)) \
    .filterDate('2014-01-01', '2015-01-01')

# Compute a median image and display.
median = collection.median()
Map.setCenter(-122.3578, 37.7726, 12)
Map.addLayer(median, {'bands': ['B4',  'B3',  'B2'], 'max': 0.3}, 'median')

In [ ]:

collection.size().getInfo()

In [ ]:

collection.aggregate_array("system:id").getInfo()

Use median reducer¶

At each location in the output image, in each band, the pixel value is the median of all unmasked pixels in the input imagery (the images in the collection). In the previous example, median() is a convenience method for the following call:

In [ ]:

# Reduce the collection with a median reducer.
median = collection.reduce(ee.Reducer.median())

# Display the median image.
Map.addLayer(median,
             {'bands': ['B4_median',  'B3_median',  'B2_median'], 'max': 0.3},
             'also median')
Map

Create an image composite¶

In [ ]:

states = ee.FeatureCollection('TIGER/2018/States')
Map.addLayer(states, {}, "US States")

In [ ]:

ca = states.filter(ee.Filter.eq("NAME", "California"))
Map.addLayer(ca, {}, "California")

In [ ]:

collection = ee.ImageCollection('LANDSAT/LC08/C01/T1_TOA') \
    .filterBounds(ca) \
    .filterDate('2020-01-01', '2021-01-01')

In [ ]:

collection.size().getInfo()

In [ ]:

image = collection.median().clip(ca)
Map.addLayer(image, {'bands': ['B4',  'B3',  'B2'], 'max': 0.3}, 'Landsat 2020')

---

Last update: 2021-03-17