¿Cálculo de trama continua que representa la proporción de tierra cubierta de raster categórica utilizando QGIS?

Question

Tengo un ráster (5 x 5 m de resolución), que representa la cobertura de la tierra. Cada celda contiene un valor entero (de 1 a 10), que representa un total de diez clases de cobertura de la tierra (por ejemplo, 1-forestal; 2-zonas agrícolas; 3-las áreas urbanas;(...)10-matorral).

Basado en esta categórica de la tierra cubierta de la trama, me gustaría producir diez rásteres con una resolución de 10 m x 10 m que representa la proporción de cada clase de cobertura terrestre por celda. En otra palabras, me gustaría producir un mapa de bits que representa la proporción de bosque, otro ráster que representa la proporción de áreas agrícolas, etc.

Esta cifra se explica exactamente lo que yo quiero hacer (para fines ilustrativos, yo sólo representan el 3 clases en la figura).

He estado buscando preguntas similares en internet, pero por lo general la gente tenía ligeramente diferentes preguntas (por ejemplo, calcular la proporción de cobertura de la tierra dentro de un polígono).

¿Alguien sabe cómo hacer esto en QGIS?

Answer 1

No tengo acceso a QGIS en este portátil, pero aquí es más o menos cómo iba a acercarse a ella:

1. Tome su 5m de uso de la tierra de trama y de uso Raster calculator a la producción de un nuevo booleano trama que tiene un valor de 1 para una clase dada, y 0 para todos los otros valores. Por ejemplo, para la clase de bosque: 5m_landuse = 1
2. El uso de la r.resamp.stats herramienta de la GRASS GIS/Raster cuadro de herramientas para agregar los nuevos booleano raster a 10 millones de células. Seleccione la resolución de salida de hasta 10m y el método de agregación a average. Por cada 10m de píxeles, este tomará el promedio de los cuatro contribuyendo 1s (bosque) y 0 (no forestales) en su 5m píxeles, lo cual termina siendo igual a una proporción (por ejemplo, promedio de 1, 1, 1 y 0 = 0.75).

Usted tendrá que repita esto para cada uno de sus diez clases de uso de la tierra, pero esto probablemente puede ser automatizado mediante el graphical modeller herramienta. Buena suerte!

Answer 2

Estoy de acuerdo con el método propuesto en la primera respuesta (primero convierta a un ráster binario por clase, tome el promedio). Sin embargo, creo que sería un poco más rápido con las herramientas gdal (de raster> conversion si no te gustan las líneas de comando).

 gdal_translate in_binary_5m.tif out_proportion_10m.tif -r average -outsize 50% 50%
 

Observación: tamaño 50% 50% establece la salida en una resolución del 50% (equivalente a 1/50% (= doble) el tamaño de los píxeles)