13_geemap学习笔记 | 按组分区统计

原文：zonal statistics by group

import ee
import geemap
geemap.set_proxy(10809)
geemap.show_youtube('cORcGGH03gg')

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

#add NLCD Data
dataset = ee.Image('USGS/NLCD/NLCD2016')
landcover = ee.Image(dataset.select('landcover'))
Map.addLayer(landcover, {}, 'landcover')

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

Map.add_legend(builtin_legend='NLCD')

Calculate land cover composite of each US states

geemap.zonal_statistics_by_group(landcover, states, 'G:/learnpy/shp/NLCD_stats1.csv', statistics_type='PERCENTAGE', denominator=1000000, decimal_palce = 2)

geemap.create_download_link('G:/learnpy/shp/NLCD_stats1.csv')

MCD12Q1.006 MODIS Land Cover Type Yearly Global 500m

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

para_vis = {
    'min': 1,
    'max': 17,
    'palette': ['006633', 'E5FFCC', '662A00', 'D8D8D8', 'F5F5F5']}

landcover_glb = ee.Image('MODIS/006/MCD12Q1/2019_01_01').select('LC_Type1')
Colormap='terrain'
Map.addLayer(landcover_glb, para_vis, 'landcover_glb')
Map.add_legend(builtin_legend = 'MODIS/006/MCD12Q1')

China = "G:/learnpy/shp/China.kml"
China = geemap.kml_to_ee(China)
Map.addLayer(China, {}, 'China')

geemap.zonal_statistics_by_group(landcover_glb, China, 'G:/learnpy/shp/China_lc_stats.csv',statistics_type='PERCENTAGE', denominator=1000000, decimal_places = 2)