Support for multibanded rasters

[RyanConway91]

Hello,

I have ~300 single band rasters with the same geometry (extent, cell size) and a shapefile with ~100000 zones.

Based on current testing, running zonal_stats on all the rasters is going to take several days.

If I could combine the rasters into a single multibanded raster, would the run time be reduced? Does zonal_stats work with multibanded rasters?

[debboutr]

Here's an example that might help you out. You can use the multiprocessing package to speed things up!