Las operaciones de reasignación o reclasificación, consisten en la modificación total o parcial de los valores de un ráster. La reclasificación permite agrupar en grupos los valores actuales de un ráster.
Estas operaciones pueden ser de utilidad en análisis de los diferentes tipos de ráster, Modelos Digitales de Elevación, ráster discretos, imágenes de satélite.
Ejemplos en ráster de variables continúas:
Denominamos raster continuos aquellos que representan una variable continua en el espacio como elevación, precipitación, ruido, entre otros, también comprende imágenes de satélite o producto generados.
1.- Modelo Digital de Elevación (MDE):
Podemos utilizar la reclasificación en un MDE para generar un mapa de hipsometría o en un raster de pendiente para agrupar por rangos.
Reclasificar utilizando el PlugIn GRASS. Hipsometría:
Para conocer las estadísticas básicas de un Raster se debe acceder a Propiedades de la capa > Histograma o Propiedades de la capa >Información, en el ejemplo los valores de altura del MDE oscilan 58 a 1635 msnm.
QGIS 3 no cuenta con una herramienta específica para reclasificar para ello podemos utilizar las herramientas de GRASS.
En el menú Procesos active la Caja de Herramientas, luego ubique GRASS > r.reclass, en la caja de dialogo se selecciona el MDE, en la caja de texto Reclass Rules colocamos la expresión que define los intervalos:
- 0 thru 100 = 1
- 100 thru 200 = 2
- 200 thru 300 = 3
- 300 thru 400 = 4
- 400 thru 600 = 5
- 600 thru 1000 = 6
- 1000 thru 1600 = 7
- 1600 thru 2000 = 8
Se ejecuta la herramienta y genera el ráster de salida, que por defecto tiene el nombre Reclassified. Este procedimiento facilitara la cuantificación de la superficie de cada clase y la representación cartográfica.
Finalmente se puede alterar el texto en las Etiquetas y guardar el estilo.
Reclasificar utilizando la Calculadora Raster. Pendiente:
El ráster de pendiente en porcentaje generado del MDE anterior presenta valores de 0 a 238, se aplica una reclasificación en rangos que permitan una mejor interpretación del relieve:
- 0 – 25% Terrenos planos a inclinados
- 25 – 50% Terrenos ondulados
- 50 – 75% Terrenos escarpados
- 75 – 100% Terrenos muy escarpados
Para ello se introduce en la Calculadora Raster la expresión:
("Pendiente@1"<= 25) * 25 + (("Pendiente@1" > 25) AND ("Pendiente@1" <= 50)) * 50 + (("Pendiente@1"> 50) AND ("Pendiente@1" <= 75)) * 75 + ("Pendiente@1" > 75) * 100
Explicación: el operador de * asigna un valor a los pixels que cumplen la condición
Como el ráster generado cuenta con solo 4 categorías se puede aplicar un estilo por Único Valor.
Generando estadísticas: para cuantificar la superficie ocupada por cada rango de pendiente seleccione:
Caja de Herramientas de Procesos > Análisis Raster > Informe de valores únicos capa raster
2.- Ejemplo en imagines de satélite. Reclasificando ráster punto flotante:
La reclasificación de un NDVI puede ayudar al proceso de interpretación de una imagen de satélite, también puede ser útil para evaluar el resultado de una clasificación.
Los valores del NDVI son relacionados con la presencia o ausencia de vegetación, valores bajos generalmente coinciden con suelos desnudos o cuerpos de agua.
Comparando los valores del NDVI con la imagen de satélite se definen los rangos
- < 0 agua
- 0 a 0,3 suelos desnudos zonas urbanas
- 0,3 a 0,65 herbazales cultivos
- 0,65 a 1 cultivos densos bosques
Los valores del NDVI son de tipo punto flotante (Float) con 14 decimales, la Calculadora Raster tiene la limitación que no realiza transformaciones a enteros.
Antes de realizar la reclasificación se multiplica el NDVI por una constante, en este caso se utilizó 1000:
Expresión: ndvi@1*1000 llamamos al Raster de salida ndvi1000
Luego seleccione el menú Raster > Conversión > Traducir (Convertir Formato), para convertir el ráster generado a entero.
En la caja de dialogo en Capa de entrada seleccionamos ndvi1000, configuramos el SRC, en Parámetros Avanzados > Tipo de Datos de Salida seleccionamos Int32, ejecutamos.
Se genera un ráster denominado Convertido cuyos valores oscilan de 0 a 874.
Realizamos la reclasificación introduciendo la siguiente expresión en la Calculadora Raster:
("Convertido@1"<= 0) * 1 + (("Convertido@1" > 0) AND ("Convertido@1" <= 300)) * 2 + (("Convertido@1"> 300) AND ("Convertido@1" <= 650)) * 3 + ("Convertido@1" > 650) * 4
Finalmente configuramos el estilo y calculamos las estadísticas
Como hacerlo con un script, tengo q reclasificar 1200 imagenes, luego transformarlas a polígonos para saber el area.
Saludos Cristhian, de la forma como muestro en el articulo no es viable en tu caso
Por otro lado, tengo entendido que la versión 3,2 Bonn cuenta con la herramienta reclasificar, también puedes conseguir el modulo reclasificar en SAGA y GRASS, como muestro en el articulo
En estos módulos tienes la opción de trabajar en Batch, tengo entendido que esto se refiere a realizar de una vez múltiples corridas puedes probar
Hola, muy buena guia
Por error en la expresión introducida en la calculadora raster para proceder a la reclasificación del ndvi en cuatro clases se repite la clase “2” dos veces, la tres es la que debiera ser la 4 y la 4 no aparece
debería ser
(“Convertido@1” 0) AND (“Convertido@1” 300) AND (“Convertido@1” 650) * 4
saludos
Saludos, muchas gracias por tus comentarios Alvaro, es agradable saber que la información es útil y reproducida, tienes toda la razón, la corrección ya fue realizada
Muy buenos dias, muy útil el ejemplo.
¿Cómo puedo restar el valor mínimo? En esta versión min(raster) no funciona como en las anteriores
Gracias por la ayuda