Tengo un archivo raster sobre el uso de la tierra en el norte de Italia (cultivos / áreas urbanas, etc.)
He definido algunas zonas de amortiguación alrededor de algunos puntos dados de una capa vectorial y me gustaría calcular el% de esas áreas ocupadas por cada valor del archivo ráster para hacer una estimación del uso promedio de las áreas seleccionadas (como 70% de cultivos, 30% urbanos, etc.).
Sé que en ArcGIS hay una función de análisis espacial que se llama " Área de tabulación ", pero no puedo encontrar dicha función en QGIS.
Si alguien me puede ayudar, estaría muy agradecido.
Respuestas:
intente utilizar "Clasificación cruzada y tabulación" de SAGA GIS del complemento Sextante.
Espero que esto ayude
fuente
Codifiqué un complemento para QGIS llamado LecoS (mira en Plugin Downloader, complementos experimentales, ¡instala Scipy primero!). Puede lograr lo que desea (porcentajes de cobertura del suelo para un determinado ráster), sin embargo, hasta ahora no he encontrado el tiempo para codificar una superposición de polígonos en Python y, por lo tanto, debe recortar su ráster individualmente por búfer y función vectorial de antemano. Vea un caso de uso aquí .
fuente
Para este caso, puede usar la Clasificación cruzada y la Tabulación como se mencionó anteriormente. Para usarlo, primero debe tener una capa ráster con valores int. Dependiendo de su capa de uso, necesitaría realizar una reclasificación o recodificación para tener listo su ráster.
La otra premisa es que sus polígonos (sus zonas de amortiguamiento) deben tener algún valor int. Si no, debe crearlo (con la fórmula de la calculadora $ id + 1).
Rasterice su capa vectorial usando el código int como valor ráster. Aquí debe poner algunas condiciones: use el mismo valor de píxel y la misma extensión que el ráster en el que desea extraer la información.
Ahora use la clasificación cruzada y la tabulación de las herramientas SAGA. Debe tener en cuenta que esto creará una tabla xy con el recuento de píxeles de cada categoría de uso en cada polígono. Aquí debe insertar el "número máximo de clases", lo que significa que, por ejemplo, si tiene un ráster con 10 categorías y 25 clases de polígonos, su máximo será 25.
Aquí está el resultado final: una tabla donde tienes el recuento de píxeles de cada categoría en cada polígono. Si necesita convertirlo en área, simplemente multiplique el número de celdas por la resolución de su ráster.
fuente
La herramienta de áreas tabuladas se encuentra en la caja de herramientas de Spatial Analyst en la sección de herramientas zonales para ArcGIS 10. No estoy seguro de qué versión está ejecutando. Tenga en cuenta que debe tener una licencia separada para la extensión Spatial Analyst para usar esas herramientas. Parece que esa herramienta hará lo que buscas.
Espero que esto ayude
fuente