xagg
A package to aggregate gridded data in xarray
to polygons in geopandas
using area-weighting from the relative area overlaps between pixels and polygons.
The simplest code run, involving raster data in an xarray
Dataset
ds
and polygons in a geopandas
GeoDataFrame
gpd
, is:
import xagg as xa
# Get overlap between pixels and polygons
weightmap = xa.pixel_overlaps(ds,gdf)
# Aggregate data in [ds] onto polygons
aggregated = xa.aggregate(ds,weightmap)
aggregated
can then be turned into an xarray
Dataset
, a geopandas
GeoDataFrame
, or directly exported to a CSV (for use in, e.g., STATA), NetCDF, or Shapefile.
Contents
- Intro
- Installation
- Base run
- Detailed code run
- Tips
- xagg