Cloudless MOD09Q1 Median Composite

Tutorial created by **David Montero Loaiza**: GitHub | Twitter

GitHub Repo: https://github.com/davemlz/eemont
PyPI link: https://pypi.org/project/eemont/
Conda-forge: https://anaconda.org/conda-forge/eemont
Documentation: https://eemont.readthedocs.io/
More tutorials: https://github.com/davemlz/eemont/tree/master/docs/tutorials

Let’s start!

If required, please uncomment:

[1]:

#!pip install eemont
#!pip install geemap

Import the required packages.

[2]:

import ee, eemont, geemap

Authenticate and Initialize Earth Engine and geemap.

[3]:

Map = geemap.Map()

Get, filter and mask clouds and shadows in the MOD09Q1 MODIS Product. Also, compute the median composite.

[4]:

MOD09Q1 = (ee.ImageCollection('MODIS/006/MOD09Q1')
          .filterDate('2020-01-01','2021-01-01')
          .maskClouds()
          .scaleAndOffset()
          .median())

Visualization

Let’s define the visualization.

[5]:

visParams = {
    'min':0,
    'max':1,
    'bands':['sur_refl_b02', 'sur_refl_b02', 'sur_refl_b01']
}

Use geemap to display results:

[6]:

Map.addLayer(MOD09Q1,visParams,'MOD09Q1 Cloudless Median Composite')
Map

Clouds/Shadows Masking on Sentinel-2 Surface Reflectance Product

Clouds Masking on Sentinel-3 Radiance (Bright Pixels Masking)