I am trying to apply the rasterio.features.rasterize function in order to obtain a zonal statistics from a vector file containing multiple vectors.
Despite my efforts, I have been unsuccessful in this task. All I can do is rasterize one single shapefile-feature-file per time.
Would anyone know how to properly rasterize multiple features using rasterio?
My ultimate goal would be an algorithm so that:
For each iteration, each vector feature would produce a single matriz mask in memory. This memory-matriz would then be applied over my raster data with same dimensions as that one. Zonal statistics as (mean, maximum, minimum, std, var) would then be inserted as new properties of each rasterized feature of my original vector (shapefile) file.