Extraire des valeurs raster dans un fichier de formes avec pygeoprocessing ou gdal

Je voudrais savoir comment obtenir toutes les valeurs raster dans un polygone en utilisant gdal ou pygeoprocessing, sans lire la grille entière comme un tableau. pygeoprocessing et gdal peuvent faire des statistiques zonales mais seuls les min, max, moyenne, stdev ou count sont disponibles à partir...