Extraiga valores ráster dentro de shapefile con pygeoprocessing o gdal
Me gustaría saber cómo obtener todos los valores ráster dentro de un polígono usando gdal o pygeoprocessing, sin leer toda la cuadrícula como una matriz. pygeoprocessing y gdal pueden hacer estadísticas zonales, pero solo la función min, max, mean, stdev o count están disponibles desde dicha...