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

10

Quiero comparar mis DEM entre sí (ASTER GDEM2 y SRTM4.1). Sé que el GDEM tiene algunos "agujeros" con elevaciones incorrectas (encontré 3 de ellos en mi área de trabajo). Ahora necesito saber si hay más de ellos para corregirlos eventualmente con el SRTM. Yo uso SAGA, GRASS y QGIS.

— MAPA
fuente

compararlos con qué? ¿juntos? uno es tu línea de base y el otro es defectuoso?

— nickves

¡derecho! Solo quiero dibujar la diferencia entre ellos.

— MAPA

5

Digamos que Raster1 es el correcto, y Raster2 es el que sospecha que tiene fallas.

Ambos también están en las mismas unidades (dado que usted dice que ambos rásteres son DEM, el valor z debe estar en metros; supongo que estamos en lo correcto y comparamos las mismas cosas)

Puede usar rast calc de cualquier paquete (qgis, grass, saga) para compararlos, por ejemplo:

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

cada celda del ráster resultante mostrará el porcentaje que ráster2 se desvía del ráster1.

Después de eso, 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 áreas problemáticas y elegir si desea intervenir o no

— nickves
fuente

8

La relación no es una buena idea para los DEM con elevaciones cerca del nivel del mar: obtendrá muchos resultados falsos y resultados completamente indefinidos donde sea raster2 == 0. El enfoque correcto, para todos los DEM, es restar un ráster del otro y la pantalla para valores atípicos (que se pueden hacer de muchas, muchas formas: consulte stats.stackexchange.com para conocer algunos métodos).

— whuber

11

Puede comenzar obteniendo la diferencia de los dos DEM. QGIS tiene una herramienta de calculadora ráster que debería ser útil. Simplemente obtenga mosaicos de ambos DEM que cubran la misma área y reste los valores de un DEM del otro. Eso debería proporcionarle una buena capa ráster que muestra las diferencias de elevación entre los dos DEM.

— RK
fuente

55

Exactamente, solo resta un DEM del otro. Luego lo representaría en QGIS o incluso lo reclasificaría, usando cierta tolerancia alrededor de Cero. Algo así como Si -5m <Resultado <5m, entonces considere que los DEM tienen el mismo valor en esa celda en particular.

— Alexandre Neto

66

(+1) Ese es un buen comienzo, y también simple. Los errores en los DEM son un poco complicados: son el resultado de errores de elevación y de errores posicionales. Si no hay un error posicional, normalmente una pantalla basada en una tolerancia constante es buena, pero no necesariamente debe estar alrededor de cero, porque un DEM podría tener un desplazamiento vertical constante (pequeño) con respecto a otro. Use una tolerancia alrededor de la diferencia mediana en su lugar. En áreas montañosas, el error posicional agrega más variación donde las pendientes son más pronunciadas, por lo que la tolerancia debe incrementarse aproximadamente en proporción a la tangente de la pendiente.

— whuber

Tengo un problema similar, pero mis rásteres se superponen solo parcialmente (el ráster B cubre solo una parte del ráster A) y la herramienta "diferencia de cuadrícula" no funcionará. ¿Alguna idea de lo que puedo hacer?

— H.Wiener