Incheon

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 = ee.Geometry.Point([126.46560976562499, 37.520194763339695]),
    preIncheon = {"opacity":1,"bands":["red","nir1","green"],"min":3.9,"max":143.1,"gamma":0.711},
    postIncheon = {"opacity":1,"bands":["red","nir","green"],"min":7196.82,"max":16980.18,"gamma":0.814}
#**** 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'])
landsat3 = ee.ImageCollection('LANDSAT/LM03/C01/T2') \
    .select(['B4', 'B5', 'B6', 'B7'], ['green', 'red', 'nir1', 'nir2'])

preImage = landsat3 \
    .filterBounds(point) \
    .filterDate('1981-01-01', '1981-12-30') \
    .sort('CLOUD_COVER', True) \
    .first()
postImage = landsat8 \
    .filterBounds(point) \
    .filterDate('2020-01-01', '2020-12-30') \
    .sort('CLOUD_COVER', True) \
    .first()

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

# Add Earth Engine dataset image = ee.Image("USGS/SRTMGL1_003") #*** Start of imports. If edited, may not auto-convert in the playground. ***# point = ee.Geometry.Point([126.46560976562499, 37.520194763339695]), preIncheon = {"opacity":1,"bands":["red","nir1","green"],"min":3.9,"max":143.1,"gamma":0.711}, postIncheon = {"opacity":1,"bands":["red","nir","green"],"min":7196.82,"max":16980.18,"gamma":0.814} #**** 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']) landsat3 = ee.ImageCollection('LANDSAT/LM03/C01/T2') \ .select(['B4', 'B5', 'B6', 'B7'], ['green', 'red', 'nir1', 'nir2']) preImage = landsat3 \ .filterBounds(point) \ .filterDate('1981-01-01', '1981-12-30') \ .sort('CLOUD_COVER', True) \ .first() postImage = landsat8 \ .filterBounds(point) \ .filterDate('2020-01-01', '2020-12-30') \ .sort('CLOUD_COVER', True) \ .first() Map.centerObject(point, 10) Map.addLayer(preImage, preIncheon, 'pre') Map.addLayer(postImage, postIncheon, 'post')

Display the interactive map¶

In [ ]:

Map

Map