12 zonal statistics

Uncomment the following line to install geemap if needed.

In [ ]:

# !pip install geemap

# !pip install geemap

In [ ]:

import ee
import geemap
import os

import ee import geemap import os

In [ ]:

geemap.show_youtube('ou-Xm3CLitM')

geemap.show_youtube('ou-Xm3CLitM')

In [ ]:

Map = geemap.Map()
Map

Map = geemap.Map() Map

In [ ]:

# Add Earth Engine dataset
dem = ee.Image('USGS/SRTMGL1_003')

# Set visualization parameters.
dem_vis = {
    'min': 0,
    'max': 4000,
    'palette': ['006633', 'E5FFCC', '662A00', 'D8D8D8', 'F5F5F5'],
}

# Add Earth Engine DEM to map
Map.addLayer(dem, dem_vis, 'SRTM DEM')

# Add Landsat data to map
landsat = ee.Image('LANDSAT/LE7_TOA_5YEAR/1999_2003')

landsat_vis = {'bands': ['B4', 'B3', 'B2'], 'gamma': 1.4}
Map.addLayer(landsat, landsat_vis, "LE7_TOA_5YEAR/1999_2003")

states = ee.FeatureCollection("TIGER/2018/States")
Map.addLayer(states, {}, 'US States')

# Add Earth Engine dataset dem = ee.Image('USGS/SRTMGL1_003') # Set visualization parameters. dem_vis = { 'min': 0, 'max': 4000, 'palette': ['006633', 'E5FFCC', '662A00', 'D8D8D8', 'F5F5F5'], } # Add Earth Engine DEM to map Map.addLayer(dem, dem_vis, 'SRTM DEM') # Add Landsat data to map landsat = ee.Image('LANDSAT/LE7_TOA_5YEAR/1999_2003') landsat_vis = {'bands': ['B4', 'B3', 'B2'], 'gamma': 1.4} Map.addLayer(landsat, landsat_vis, "LE7_TOA_5YEAR/1999_2003") states = ee.FeatureCollection("TIGER/2018/States") Map.addLayer(states, {}, 'US States')

In [ ]:

out_dir = os.path.join(os.path.expanduser('~'), 'Downloads')
out_dem_stats = os.path.join(out_dir, 'dem_stats.csv')

if not os.path.exists(out_dir):
    os.makedirs(out_dir)

# Allowed output formats: csv, shp, json, kml, kmz
# Allowed statistics type: MEAN, MAXIMUM, MINIMUM, MEDIAN, STD, MIN_MAX, VARIANCE, SUM
geemap.zonal_statistics(dem, states, out_dem_stats, statistics_type='MEAN', scale=1000)

out_dir = os.path.join(os.path.expanduser('~'), 'Downloads') out_dem_stats = os.path.join(out_dir, 'dem_stats.csv') if not os.path.exists(out_dir): os.makedirs(out_dir) # Allowed output formats: csv, shp, json, kml, kmz # Allowed statistics type: MEAN, MAXIMUM, MINIMUM, MEDIAN, STD, MIN_MAX, VARIANCE, SUM geemap.zonal_statistics(dem, states, out_dem_stats, statistics_type='MEAN', scale=1000)

In [ ]:

out_landsat_stats = os.path.join(out_dir, 'landsat_stats.csv')
geemap.zonal_statistics(
    landsat, states, out_landsat_stats, statistics_type='SUM', scale=1000
)

out_landsat_stats = os.path.join(out_dir, 'landsat_stats.csv') geemap.zonal_statistics( landsat, states, out_landsat_stats, statistics_type='SUM', scale=1000 )

In [ ]:

geemap.create_download_link(out_dem_stats)

geemap.create_download_link(out_dem_stats)

In [ ]:

geemap.create_download_link(out_landsat_stats)

geemap.create_download_link(out_landsat_stats)

---

Last update: 2023-04-06
Created: 2020-03-31