Skip to article frontmatterSkip to article content

GEE Map!

DSE EcoTech Series (2024)

The Eric and Wendy Schmidt Center for Data Science & Environment
University of California Berkeley

license: CC-BY-4.0


IMPORTS

import geemap
import ee
ee.Initialize()

GEE IMPORTS

FIRES = ee.FeatureCollection("users/brookie/CALFIRES/Boundaries")
S2 = ee.ImageCollection("COPERNICUS/S2_SR_HARMONIZED")
Loading...

CONFIG

START_DATE = '2021-01-01'
END_DATE = '2021-04-01'
REGION = ee.Geometry.Polygon(
    [[[-124.20390625, 41.05342459748544],
      [-124.643359375, 40.25328623978284],
      [-123.588671875, 39.375667620490226],
      [-122.006640625, 36.393606666819615],
      [-119.809375, 34.460099309324185],
      [-118.0515625, 33.73228855323828],
      [-117.17265625, 32.6289197967587],
      [-114.184375, 32.62891979675872],
      [-114.184375, 34.60490865998548],
      [-119.897265625, 38.83003685201642],
      [-119.98515625, 42.10531374703366],
      [-124.291796875, 42.04007296056778]]])
S2_RESCALE = 1e-4
Loading...
ZOOM = 6
S2_TITLE = '(medain) Sentinel-2'
CALFIRES_TITLE = 'Historical Fire Boundaries [1898 - 2021]'
S2_VIS = dict(
  min=0.0,
  max=0.3,
  bands=['B4', 'B3', 'B2'])
CALFIRES_VIS = dict(color='f720b9')
CALFIRES_ALPHA = 0.75
Loading...
S2_TITLE += ' ['+START_DATE+' to '+END_DATE+']'
DATA_FILTER = ee.Filter.And(
    ee.Filter.bounds(REGION),
    ee.Filter.date(START_DATE, END_DATE))
Loading...

PROCESS SENTINEL-2 DATA AND VISUALIZE

  1. filter by date and region
  2. create median composite
  3. rescale
  4. add to map
S2 = S2.filter(DATA_FILTER).median().multiply(S2_RESCALE)
Map = geemap.Map(center=[40,-100], zoom=4)
Map.centerObject(REGION, ZOOM);
Map.addLayer(S2, S2_VIS, S2_TITLE);
Map.addLayer(FIRES, CALFIRES_VIS, CALFIRES_TITLE, True, CALFIRES_ALPHA)
Map
Loading...