zonal_stats function returns 'infinity' value

[shubham-scisar]

zonal_stats function returns infinity when calculating mean for some geometries.

[ircwaves]

zonal_stats function returns infinity when calculating mean for some geometries.

If there's an inf within the portion of the raster delineated by the polygon, that seems mathematically correct.

[shubham-scisar]

zonal_stats function returns infinity when calculating mean for some geometries.

If there's an inf within the portion of the raster delineated by the polygon, that seems mathematically correct.

There is no inf value in the raster.Also, I think inside the main script of rasterstats where mean is being calculated as float(masked.mean()), if it is changed to np.mean(masked.astype('float64')) (as the data type of my raster is float32, the problem vanishes). So, I think this limitation of numpy should be taken careof in the versions to come.

[perrygeo]

@shubhamsharma1609 can you provide a minimal example with code and data to reproduce this?