GEE Map!
DSE EcoTech Series (2024)
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¶
- filter by date and region
- create median composite
- rescale
- 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...