Tengo dos capas. Una capa de forma de polígono con muchos mosaicos y una capa ráster que contiene Cobertura de suelo CORINE 2006 con muchas categorías en un mapa de colores. Quiero obtener para cada polígono en la capa de formas una suma de cada categoría de cobertura terrestre de la capa ráster.
Por ejemplo, hay un polígono con id '2' y quiero atributos como este para este polígono (en porcentaje o en metros cuadrados):
- Tierra arable: 15%
- Bosque: 11%
- Calles: 2% (... y así uno)
Intenté hacerlo en grass, qgis (sin función), saga (solo resume todo a un valor total) r (suma total), pero todavía no encontré solución. La mayoría de los complementos (estadísticas zonales en qgis) solo admiten 0-1 capas ráster. v.rast.stats tampoco ayudó. Estoy abierto a cualquier solución buena e inteligente !. Tal vez incluso utilicé un enfoque equivocado o cometí errores.
En Arcgis, esta tarea es bastante fácil, si me acuerdo bien, pero todavía me falta una buena solución para su usuario de Linux diario.
Estoy ejecutando un sistema Debian Linux y por eso solo puedo usar programas para este sistema operativo.
EDITAR: Porque esta pregunta todavía tiene tantas visitas y visitantes: Escribí un complemento de QGIS, que también es capaz de calcular la cobertura terrestre de la capa ráster. No he codificado una superposición de polígonos todavía, pero definitivamente planeado. Encuentre el complemento aquí e instale la biblioteca Scipy primero.