Wie kann man zwei DEMs vergleichen und ihren Unterschied zeigen?

10

Ich möchte meine DEMs miteinander vergleichen (ASTER GDEM2 und SRTM4.1). Ich weiß, dass das GDEM einige "Löcher" mit falschen Höhen hat (ich habe 3 davon in meinem Arbeitsbereich gefunden). Jetzt muss ich wissen, ob es mehr davon gibt, um sie schließlich mit dem SRTM zu korrigieren. Ich benutze SAGA, GRASS und QGIS.

KARTE
quelle
Vergleiche sie mit was? miteinander? einer ist deine Grundlinie und der andere ist fehlerhaft?
Nickves
Recht! Ich möchte nur den Unterschied zwischen ihnen ziehen.
Karte

Antworten:

5

Nehmen wir an, Raster1 ist das Richtige, und Raster2 ist das, von dem Sie vermuten, dass es Fehler aufweist.

Sie sind auch beide in den gleichen Einheiten (da Sie sagen, dass beide Raster DEMs sind, sollte der z-Wert in Metern sein - ich gehe davon aus, dass wir korrekt sind und wir vergleichen die gleichen Dinge)

Sie können Rast Calc aus beiden Paketen (QGIS, Gras, Saga) verwenden, um sie zu vergleichen, z.

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

Jede Zelle aus dem resultierenden Raster zeigt den Prozentsatz an, um den Raster2 von Raster1 abweicht.

Danach können Sie Abweichung_rast basierend auf Ihren Anforderungen neu klassifizieren: z

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

um problematische Bereiche zu identifizieren - und um zu entscheiden, ob Sie eingreifen möchten oder nicht

Nickves
quelle
8
Das Verhältnis ist keine gute Idee für DEMs mit Höhen nahe dem Meeresspiegel: Sie erhalten überall dort, wo raster2 == 0 ist, viele falsche Ergebnisse und völlig undefinierte Ergebnisse. Der richtige Ansatz für alle DEMs besteht darin, nur ein Raster vom anderen zu subtrahieren und zu screenen für Ausreißer (was auf viele, viele Arten möglich ist: konsultieren Sie stats.stackexchange.com für einige Methoden).
whuber
11

Sie könnten damit beginnen, die Differenz der beiden DEMs zu ermitteln. QGIS verfügt über ein Rasterrechner-Tool, das sich als nützlich erweisen sollte. Holen Sie sich einfach Kacheln beider DEMs, die denselben Bereich abdecken, und subtrahieren Sie die Werte eines DEM vom anderen. Dadurch erhalten Sie eine schöne Rasterebene, die die Höhenunterschiede zwischen den beiden DEMs anzeigt.

RK
quelle
5
Genau, subtrahieren Sie einfach ein DEM vom anderen. Ich würde es dann in QGIS darstellen oder sogar neu klassifizieren, wobei ich eine Toleranz um Null verwenden würde. Etwas wie Wenn -5m <Ergebnis <5m, dann bedenken Sie, dass DEMs in dieser bestimmten Zelle den gleichen Wert haben.
Alexandre Neto
6
(+1) Das ist ein guter Anfang - und auch einfach. Fehler in DEMs sind etwas kompliziert: Sie resultieren aus Höhenfehlern und Positionsfehlern. Wenn kein Positionsfehler vorliegt, ist ein Bildschirm, der auf einer konstanten Toleranz basiert, normalerweise gut - er sollte jedoch nicht unbedingt um Null liegen, da ein DEM möglicherweise eine (kleine) konstante vertikale Verschiebung relativ zu einem anderen aufweist. Verwenden Sie stattdessen eine Toleranz um den Medianunterschied. In hügeligen Gebieten führt der Positionsfehler zu einer größeren Variation, wenn die Steigungen steiler sind. Daher muss die Toleranz ungefähr proportional zur Tangente der Steigung erhöht werden.
whuber
Ich habe ein ähnliches Problem, aber meine Raster überlappen sich nur teilweise (Raster B deckt nur einen Teil von Raster A ab) und das Werkzeug "Rasterunterschied" funktioniert nicht. Irgendwelche Ideen, was ich tun kann?
H. Wiener