mtStHelen

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. ***#
point =

    # shown: False #
    ee.Geometry.Point([-122.19174741047361, 46.20040199038881])
#**** End of imports. If edited, may not auto-convert in the playground. ****#
landsat2 = ee.ImageCollection('LANDSAT/LM02/C02/T2') \
    .select(['B4', 'B5', 'B6', 'B7'],
            ['green', 'red', 'nir1', 'nir2'])

preImage = landsat2 \
    .filterBounds(point) \
    .filterDate('1979-08-01', '1979-10-30') \
    .sort('CLOUD_COVER', True) \
    .first()
postImage = landsat2 \
    .filterBounds(point) \
    .filterDate('1980-04-01', '1981-10-30') \
    .sort('CLOUD_COVER', True) \
    .first()

visParam = {
  'opacity': 1,
  'bands': ['nir1', 'red', 'green'],
  'min': 0,
  'max': 128,
  'gamma':1
}

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

# Add Earth Engine dataset image = ee.Image("USGS/SRTMGL1_003") #*** Start of imports. If edited, may not auto-convert in the playground. ***# point = # shown: False # ee.Geometry.Point([-122.19174741047361, 46.20040199038881]) #**** End of imports. If edited, may not auto-convert in the playground. ****# landsat2 = ee.ImageCollection('LANDSAT/LM02/C02/T2') \ .select(['B4', 'B5', 'B6', 'B7'], ['green', 'red', 'nir1', 'nir2']) preImage = landsat2 \ .filterBounds(point) \ .filterDate('1979-08-01', '1979-10-30') \ .sort('CLOUD_COVER', True) \ .first() postImage = landsat2 \ .filterBounds(point) \ .filterDate('1980-04-01', '1981-10-30') \ .sort('CLOUD_COVER', True) \ .first() visParam = { 'opacity': 1, 'bands': ['nir1', 'red', 'green'], 'min': 0, 'max': 128, 'gamma':1 } Map.centerObject(point, 10) Map.addLayer(preImage, visParam, 'pre') Map.addLayer(postImage, visParam, 'post')

Display the interactive map¶

In [ ]:

Map

Map