GDAL

Polygonize Raster and Compute Zonal-Statistics in Python

Computing zonal statistics over raster using GDAL and Python