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

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()

Point of interest.

[4]:

point = ee.Geometry.Point([-75.92, 2.86])

Get and filter the Sentinel-3 Radiance image collection and filter it by region and time.

[5]:

S3 = (ee.ImageCollection('COPERNICUS/S3/OLCI')
      .filterBounds(point)
      .filterDate('2020-01-01','2021-01-01'))

Clouds Masking

The clouds in Sentinel-3 are masked according to the Bright Pixels Bitmask in the quality_flags band using the maskClouds() method.

[6]:

S3masked = (S3.maskClouds()
            .scaleAndOffset() # NEW! Note that we also can scale this product!
            .first())

Visualization

Let’s define the RGB visualization.

[7]:

rgbVis = {
    'min':0,
    'max':6,
    'bands':['Oa08_radiance', 'Oa06_radiance', 'Oa04_radiance'],
    'gamma':1.5
}

Use geemap to display results:

[9]:

Map.centerObject(point,10)
Map.addLayer(S3.scale().first(),rgbVis,'Sentinel-3 Original')
Map.addLayer(S3masked,rgbVis,'Sentinel-3 With Cloud Mask')
Map

Cloudless MOD09Q1 Median Composite

Masking Clouds and Shadows in VIIRS Products