A21a 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")

#  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#  Chapter:      A2.1 Groundwater Monitoring with GRACE
#  Checkpoint:   A21a
#  Authors:      A.J. Purdy, J.S. Famiglietti
#  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Import Basins.
basins = ee.FeatureCollection("USGS/WBD/2017/HUC04")

# Extract the 3 HUC 04 basins for the Central Valley.
codes = ["1802", "1803", "1804"]
basin = basins.filter(ee.Filter.inList("huc4", codes))

# Add the basin to the map to show the extent of our analysis.
Map.centerObject(basin, 6)
Map.addLayer(basin, {"color": "green"}, "Central Valley Basins", True, 0.5)

landcover = (
    ee.ImageCollection("USDA/NASS/CDL")
    .filter(ee.Filter.date("2019-01-01", "2019-12-31"))
    .select("cultivated")
)

Map.addLayer(landcover.first().clip(basin), {}, "Cropland", True, 0.5)

# This table was generated using the index from the CDEC website
res = ee.FeatureCollection("projects/gee-book/assets/A2-1/ca_reservoirs_index")
# Filter reservoir locations by the Central Valley geometry
res_cv = res.filterBounds(basin)
Map.addLayer(res_cv, {"color": "blue"}, "Reservoirs")

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

# Add Earth Engine dataset image = ee.Image("USGS/SRTMGL1_003") # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Chapter: A2.1 Groundwater Monitoring with GRACE # Checkpoint: A21a # Authors: A.J. Purdy, J.S. Famiglietti # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Import Basins. basins = ee.FeatureCollection("USGS/WBD/2017/HUC04") # Extract the 3 HUC 04 basins for the Central Valley. codes = ["1802", "1803", "1804"] basin = basins.filter(ee.Filter.inList("huc4", codes)) # Add the basin to the map to show the extent of our analysis. Map.centerObject(basin, 6) Map.addLayer(basin, {"color": "green"}, "Central Valley Basins", True, 0.5) landcover = ( ee.ImageCollection("USDA/NASS/CDL") .filter(ee.Filter.date("2019-01-01", "2019-12-31")) .select("cultivated") ) Map.addLayer(landcover.first().clip(basin), {}, "Cropland", True, 0.5) # This table was generated using the index from the CDEC website res = ee.FeatureCollection("projects/gee-book/assets/A2-1/ca_reservoirs_index") # Filter reservoir locations by the Central Valley geometry res_cv = res.filterBounds(basin) Map.addLayer(res_cv, {"color": "blue"}, "Reservoirs") # ----------------------------------------------------------------------- # CHECKPOINT # -----------------------------------------------------------------------

Display the interactive map¶

In [ ]:

Map

Map