Zonal statistics using geopandas GeoDataFrame and/or shapely objects?

[vidurmithal]

Most use-cases of zonal statistics would have the geometries stored in a geopandas GeoDataFrame, rather than as an xarray DataArray. Converting a GeoDataFrame to an DataArray with "zones" is relatively simple using the regionmask package.