12_geemap学习笔记 | 分区统计

12 zonal statistics

import ee
import geemap
import os
geemap.set_proxy(10809)
geemap.show_youtube('ou-Xm3CLitM')

Map = geemap.Map(center=(38,105))
Map

# 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('LE7_TOA_5YEAR/1999_2003')

landsat_vis = {
    'band':['B4', 'B3', 'B2'],
    'gamma':1.4
}

Map.addLayer(landsat,landsat_vis,'LE7_TOA_5YEAR/1999_2003')

provences = "G:/learnpy/shp/province.shp"
provences_ee = geemap.shp_to_ee(provences)
Map.addLayer(provences_ee,{},'Provences')

out_dem_stats = 'G:/learnpy/data/dem_stats.csv'
geemap.zonal_statistics(dem, provences_ee, out_dem_stats, statistics_type='MEAN', scale = 1000)

# Computing statistics ...
# Generating URL ...
# Downloading data from https://earthengine.googleapis.com/v1alpha/projects/earthengine-legacy/tables/9025e99d374e6a8d142bbb2ff18be2fa-aaf1bcaac1aac1e6c156c8982ff63ce8:getFeatures
# Please wait ...
# Data downloaded to G:\learnpy\data\dem_stats.csv

out_landsat_stats = 'G:/learnpy/data/landsat_stats.csv'
geemap.zonal_statistics(landsat, provences_ee, out_landsat_stats,statistics_type='SUM', scale = 1000)

# out_dem_stats = 'G:/learnpy/data/dem_stats.csv'
# geemap.zonal_statistics(dem, provences_ee, out_dem_stats, statistics_type='MEAN', scale = 1000)

# # Computing statistics ...
# # Generating URL ...
# # Downloading data from https://earthengine.googleapis.com/v1alpha/projects/earthengine-legacy/tables/9025e99d374e6a8d142bbb2ff18be2fa-aaf1bcaac1aac1e6c156c8982ff63ce8:getFeatures
# # Please wait ...
# # Data downloaded to G:\learnpy\data\dem_stats.csv

geemap.create_download_link(out_dem_stats)

geemap.create_download_link(out_landsat_stats)