Comparar dos extensiones raster en QGIS

Question

Tengo dos cuadrículas, la cuadrícula A representa las extensiones observadas de inundaciones y la cuadrícula B representa las extensiones simuladas de un modelo hidráulico. Los rásteres tienen el mismo alcance y tamaño de cuadrícula.

Quiero poder comparar cada píxel y crear un nuevo ráster con 4 valores únicos (A/B/C/D) equivalente a la tabla a continuación para mostrar la diferencia entre los dos rásteres:

Si esto fuera un problema de una calculadora de campo de shapefile, usaría una lógica de CASO como la siguiente, pero estoy luchando por identificar un equivalente basado en ráster de QGIS. ¿Alguien puede sugerir algo que pueda intentar?

CASO
    CUANDO "Modelado" > 0 Y "Observado" > 0 LUEGO D
    CUANDO "Modelado" = 0 Y "Observado" > 0 LUEGO C
    CUANDO "Modelado" > 0 Y "Observado" = 0 LUEGO B
    CUANDO "Modelado" = 0 Y "Observado" = 0 LUEGO A
FIN

Answer 1

Puedes usar la calculadora de ráster de QGIS para esto, utilizando máscaras booleanas. Un ejemplo:

(("grid1@1">0) AND ("grid2@1">0)) * 1
+ (("grid1@1"=0) AND ("grid2@1">0)) * 2
+ (("grid1@1">0) AND ("grid2@1"=0)) * 3
+ (("grid1@1"=0) AND ("grid2@1"=0)) * 4

donde puedes cambiar 1, 2, 3 y 4 a los valores que desees.

La forma en que funciona esto es que la primera parte de cada término es una máscara booleana. ("grid1@1">0) AND ("grid2@1">0)) es Verdadero cuando se cumple tu primera condición, y puede ser multiplicado por tu valor categórico D (yo usé 1). Los siguientes tres términos coinciden con tus siguientes tres condiciones.

Nota

Steven Kay señala en los comentarios que la misma expresión se puede realizar sin usar AND reemplazándolo con un *. En este caso, el primer término se vería así

(("grid1@1">0) * ("grid2@1">0)) * 1