Ekstrak nilai raster di dalam shapefile dengan pygeoprocessing atau gdal

Saya ingin tahu cara mendapatkan semua nilai raster dalam poligon menggunakan gdal atau pygeoprocessing, tanpa membaca seluruh grid sebagai array. pygeoprocessing dan gdal dapat melakukan statistik zona tetapi hanya min, maks, rata-rata, stdev atau hitungan tersedia dari fungsi tersebut. Karena...