I have a shapefile that I have loaded as a GeoPandas GeoDataFrame, and it has a geometry column containing polygons and multipolygons.
My end goal is to do some zonal computations on a GeoTiff raster, in particular, I want to compute the mean value within each polygon, and also count the number of pixels that contributed to that polygon's zonal mean. Here's a link to another GIS-SE question of mine that outlines my attempts to do that so far using GDAL and rasterstats.
I am just wondering if GeoPandas is up to zonal calculations yet. If I open a raster in GDAL/numpy/rasterio and load the shapefile in a GeoPandas GeoDataFrame, is there a way to compute the number of raster pixels inside each polygon?