Mining

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([-70.54473946976444, -13.046675490210111]),
    mininPostVis = {"opacity":1,"bands":["swir2","nir","red"],"min":7749,"max":22215,"gamma":1},
    miningPreVis = {"opacity":1,"bands":["swir2","nir","red"],"min":569.5106854708367,"max":3002.543077970024,"gamma":0.649}
#**** 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'])
landsat5 = ee.ImageCollection('LANDSAT/LT05/C01/T2_SR') \
    .select(
      ['B1', 'B2', 'B3', 'B4', 'B5', 'B7'],
      ['blue', 'green', 'red', 'nir', 'swir1', 'swir2'])

preImage = landsat5 \
    .filterBounds(point) \
    .filterDate('1985-01-01', '2002-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, miningPreVis, 'pre')
Map.addLayer(postImage, mininPostVis, '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([-70.54473946976444, -13.046675490210111]), mininPostVis = {"opacity":1,"bands":["swir2","nir","red"],"min":7749,"max":22215,"gamma":1}, miningPreVis = {"opacity":1,"bands":["swir2","nir","red"],"min":569.5106854708367,"max":3002.543077970024,"gamma":0.649} #**** 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']) landsat5 = ee.ImageCollection('LANDSAT/LT05/C01/T2_SR') \ .select( ['B1', 'B2', 'B3', 'B4', 'B5', 'B7'], ['blue', 'green', 'red', 'nir', 'swir1', 'swir2']) preImage = landsat5 \ .filterBounds(point) \ .filterDate('1985-01-01', '2002-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, miningPreVis, 'pre') Map.addLayer(postImage, mininPostVis, 'post')

Display the interactive map¶

In [ ]:

Map

Map