Extrahieren Sie Rasterwerte innerhalb des Shapefiles mit Pygeoprocessing oder GDAL

Ich möchte wissen, wie man alle Rasterwerte innerhalb eines Polygons mit gdal oder pygeoprocessing erhält, ohne das gesamte Raster als Array zu lesen. pygeoprocessing und gdal können zonale Statistiken erstellen, aber von einer solchen Funktion sind nur die Werte min, max, mean, stdev oder count...