Extraia valores raster dentro do shapefile com pygeoprocessing ou gdal

Gostaria de saber como obter todos os valores de varredura em um polígono usando gdal ou pygeoprocessing, sem ler toda a grade como uma matriz. pygeoprocessing e gdal podem fazer estatísticas zonais, mas apenas o min, max, mean, stdev ou count estão disponíveis nessa função. Como as estatísticas...