F60h Checkpoint

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. ***#
geometry = ee.Geometry.MultiPoint()
#**** End of imports. If edited, may not auto-convert in the playground. ****#
#  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#  Chapter:      F6.0 Advanced Raster Visualization
#  Checkpoint:   F60h
#  Authors:      Gennadii Donchyts, Fedor Baart
#  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Include packages.
palettes = require('users/gena/packages:palettes')
text = require('users/gena/packages:text')

point =  ee.Geometry.Point([-
    106.15944300895228, -74.58262940096245
])

rect =
    ee.Geometry.Polygon(
        [
            [
                [-106.19789515738981, -74.56509549360152],
                [-106.19789515738981, -74.78071448733921],
                [-104.98115931754606, -74.78071448733921],
                [-104.98115931754606, -74.56509549360152]
            ]
        ], None, False)

# Lookup the ice palette.
palette = palettes.cmocean.Ice[7]

# Show it in the console.
palettes.showPalette('Ice', palette)

# Center map on geometry.
Map.centerObject(point, 9)

# Select S1 images for the Thwaites glacier.
images = ee.ImageCollection('COPERNICUS/S1_GRD') \
    .filterBounds(rect) \
    .filterDate('2021-01-01', '2021-03-01') \
    .select('HH') \
    .filter(ee.Filter.isContained({
        'leftValue': rect,
        'rightField': '.geo'
    })) \
    .sort('system:time_start')

# Print number of images.
print(images.size())

#  -----------------------------------------------------------------------
#  CHECKPOINT
#  -----------------------------------------------------------------------

# Render images.
vis = {
    'palette': palette,
    'min': -15,
    'max': 1
}

scale = Map.getScale()
textProperties = {
    'outlineColor': '000000',
    'outlineWidth': 3,
    'outlineOpacity': 0.6
}


def func_kbg(i):
    # Use the date as the label.
    label = i.date().format('YYYY-MM-dd')
    labelImage = text.draw(label, point, scale,
        textProperties)

    return i.visualize(vis) \
        .blend(labelImage) \
        .set({
            'label': label
        }); # Keep the text property.

imagesRgb = images.map(func_kbg)












Map.addLayer(imagesRgb.first())
Map.addLayer(rect, {'color':'blue'}, 'rect', 1, 0.5)

#  -----------------------------------------------------------------------
#  CHECKPOINT
#  -----------------------------------------------------------------------

# Add Earth Engine dataset image = ee.Image("USGS/SRTMGL1_003") #*** Start of imports. If edited, may not auto-convert in the playground. ***# geometry = ee.Geometry.MultiPoint() #**** End of imports. If edited, may not auto-convert in the playground. ****# # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Chapter: F6.0 Advanced Raster Visualization # Checkpoint: F60h # Authors: Gennadii Donchyts, Fedor Baart # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Include packages. palettes = require('users/gena/packages:palettes') text = require('users/gena/packages:text') point = ee.Geometry.Point([- 106.15944300895228, -74.58262940096245 ]) rect = ee.Geometry.Polygon( [ [ [-106.19789515738981, -74.56509549360152], [-106.19789515738981, -74.78071448733921], [-104.98115931754606, -74.78071448733921], [-104.98115931754606, -74.56509549360152] ] ], None, False) # Lookup the ice palette. palette = palettes.cmocean.Ice[7] # Show it in the console. palettes.showPalette('Ice', palette) # Center map on geometry. Map.centerObject(point, 9) # Select S1 images for the Thwaites glacier. images = ee.ImageCollection('COPERNICUS/S1_GRD') \ .filterBounds(rect) \ .filterDate('2021-01-01', '2021-03-01') \ .select('HH') \ .filter(ee.Filter.isContained({ 'leftValue': rect, 'rightField': '.geo' })) \ .sort('system:time_start') # Print number of images. print(images.size()) # ----------------------------------------------------------------------- # CHECKPOINT # ----------------------------------------------------------------------- # Render images. vis = { 'palette': palette, 'min': -15, 'max': 1 } scale = Map.getScale() textProperties = { 'outlineColor': '000000', 'outlineWidth': 3, 'outlineOpacity': 0.6 } def func_kbg(i): # Use the date as the label. label = i.date().format('YYYY-MM-dd') labelImage = text.draw(label, point, scale, textProperties) return i.visualize(vis) \ .blend(labelImage) \ .set({ 'label': label }); # Keep the text property. imagesRgb = images.map(func_kbg) Map.addLayer(imagesRgb.first()) Map.addLayer(rect, {'color':'blue'}, 'rect', 1, 0.5) # ----------------------------------------------------------------------- # CHECKPOINT # -----------------------------------------------------------------------

Display the interactive map¶

In [ ]:

Map

Map