¿Cómo comparar dos MDE y mostrar su diferencia?

Question

Quiero comparar mis MDEs entre sí (ASTER GDEM2 y SRTM4.1). Sé que el GDEM tiene algunos "agujeros" con elevaciones erróneas (he encontrado 3 de ellos en mi área de trabajo). Ahora necesito saber, si hay más de ellos para eventualmente corregirlos con el SRTM. Utilizo SAGA, GRASS y QGIS.

Answer 1

Podrías empezar por obtener la diferencia de los dos MDE. QGIS tiene una herramienta de cálculo de rásteres que debería ser útil. Sólo tienes que obtener mosaicos de ambos MDE que cubran la misma área y restar los valores de un MDE del otro. Así obtendrás una bonita capa ráster que muestra las diferencias de elevación entre los dos MDE.

Answer 2

Digamos que el Raster1 es el correcto, y el Raster2 es el que sospechas que tiene fallos.

Además, ambos están en las mismas unidades (como dices que ambos rásteres son MDE, el valor de z debería estar en metros, supongo que estamos en lo cierto y comparamos lo mismo)

Puede utilizar el rast calc de cualquiera de los paquetes (qgis, grass, saga) para compararlos, por ejemplo:

deviation_rast  = abs( (raster1 - raster2) / raster1 * 100 ) 

cada celda del raster resultante mostrará el porcentaje que el raster2 se desvía del raster1.

Después puede reclasificar deviation_rast en función de sus necesidades: por ejemplo

red -> deviation_rast_value > 75
orange -> deviation_rast_value > 50
yellow ->  deviation_rast_value > 25
green -> deviation_rast_value  > 0

para identificar las áreas problemáticas - y para elegir si quiere intervenir o no