Paradise

Import libraries¶

In [ ]:

import ee
import geemap

import ee import geemap

Create an interactive map¶

In [ ]:

Map = geemap.Map(center=[40, -100], zoom=4)

Map = geemap.Map(center=[40, -100], zoom=4)

Add Earth Engine Python script¶

In [ ]:

# Add Earth Engine dataset
image = ee.Image("USGS/SRTMGL1_003")

#*** Start of imports. If edited, may not auto-convert in the playground. ***#
imageVisParam = {"opacity":1,"bands":["swir2","nir","red"],"min":7749,"max":22215,"gamma":1},
    point =

    # shown: False #
    ee.Geometry.Point([-121.55188822809693, 39.794570152006976])
#**** End of imports. If edited, may not auto-convert in the playground. ****#
landsat8 = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2') \
    .select(
      ['SR_B2', 'SR_B3', 'SR_B4', 'SR_B5', 'SR_B6', 'SR_B7'],
      ['blue', 'green', 'red', 'nir', 'swir1', 'swir2'])

preImage = landsat8 \
    .filterBounds(point) \
    .filterDate('2018-10-01', '2018-10-30') \
    .sort('CLOUD_COVER', True) \
    .first()
postImage = landsat8 \
    .filterBounds(point) \
    .filterDate('2018-12-01', '2019-04-30') \
    .sort('CLOUD_COVER', True) \
    .first()

Map.centerObject(point, 10)
Map.addLayer(preImage, imageVisParam, 'pre')
Map.addLayer(postImage, imageVisParam, 'post')

# Add Earth Engine dataset image = ee.Image("USGS/SRTMGL1_003") #*** Start of imports. If edited, may not auto-convert in the playground. ***# imageVisParam = {"opacity":1,"bands":["swir2","nir","red"],"min":7749,"max":22215,"gamma":1}, point = # shown: False # ee.Geometry.Point([-121.55188822809693, 39.794570152006976]) #**** End of imports. If edited, may not auto-convert in the playground. ****# landsat8 = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2') \ .select( ['SR_B2', 'SR_B3', 'SR_B4', 'SR_B5', 'SR_B6', 'SR_B7'], ['blue', 'green', 'red', 'nir', 'swir1', 'swir2']) preImage = landsat8 \ .filterBounds(point) \ .filterDate('2018-10-01', '2018-10-30') \ .sort('CLOUD_COVER', True) \ .first() postImage = landsat8 \ .filterBounds(point) \ .filterDate('2018-12-01', '2019-04-30') \ .sort('CLOUD_COVER', True) \ .first() Map.centerObject(point, 10) Map.addLayer(preImage, imageVisParam, 'pre') Map.addLayer(postImage, imageVisParam, 'post')

Display the interactive map¶

In [ ]:

Map

Map