Mehr

Richtiges Zusammenführen von DEM-Rastern und Eliminieren von Kanteneffekten

Richtiges Zusammenführen von DEM-Rastern und Eliminieren von Kanteneffekten


Ich verwende 10-m-DEM-Raster vom NRCS Data Gateway für einen großen Teil von Kalifornien. Das Problem ist eine seltsame Kachelung, die nach dem Zusammenführen auftritt, aber in den Originaldaten (vor dem Zusammenführen) nicht vorhanden ist. Die Gitter erzeugen Tiefpunkte, so dass Flussansammlungen usw. in sie hineingezogen werden.

Ich habe versucht, Mosaik zu einem neuen Raster zu erstellen, einen Raster-Katalog zu erstellen und in ArcGIS 10.2 zu filtern, vrt in Qgis2.4 zusammenzuführen und zu erstellen, aber das Problem besteht weiterhin.

Die Gitter stellen flache Tiefen im Krümmungs-Raster an den zusammengeführten Kanten der DEM-Raster dar. Beachten Sie, dass diese Linien die Kanten von BEIDEN des mosaikierten DEM sind, wodurch 2 Linien erzeugt werden.

Die Linien sind flache Tiefs, was ein Problem darstellt, wenn man versucht, jede Art von Hangstabilität oder hydrologischer Analyse durchzuführen ... was der springende Punkt dieser Übung ist.

Unmosaiced DEM-Raster:

Krümmungs-Raster, das aus Mosaiced-Raster mit bilinearer Interpolation generiert wurde:

BEARBEITEN: Die Daten- und Mosaikmethode waren die Probleme: Die Daten selbst haben Kanteneffekte eingebaut und ich habe den nächsten Nachbarn anstelle der bilinearen Interpolation verwendet


Es scheint, dass hier zwei Dinge passieren. Erstens sind Ihre Kacheln nicht nahtlos, da in den Überlappungsbereichen an den Kanten der Kacheln die Erhebungen nicht identisch sind. Ich kann dies bestätigen, da ich mehrere Punkte entlang des Überlappungsbereichs digitalisiert und die Höhen in beiden Raster-DEMs extrahiert und Folgendes gefunden habe:

Im Fall der beiden von Ihnen hochgeladenen Kacheln scheinen die DEMs im unteren Bereich sehr ähnlich zu sein, aber wenn Sie sich zum oberen Rand bewegen, nahmen die Höhenabweichungen auf fast einen halben Meter zu, was ziemlich signifikant ist. Sie können nicht wirklich erwarten, dass ein mosaikiertes DEM keine Kantenartefakte aufweist, wenn die zusammengesetzten Kacheln im Überlappungsbereich nicht genau gleich sind. Das heißt, es sei denn, Sie führen eine Mosaikierung mit irgendeiner Art von Überblendung durch, und selbst dann ist es wahrscheinlich, dass einige Kanteneffekte aufgrund der großen Abweichungen erhalten bleiben. Welchen der Mosaik-Typen haben Sie beim Mosaikieren der Fliesen gewählt? ArcGIS bietet einen 'Misch'-Modus.

Das zweite Problem scheint damit zu tun zu haben, wie ArcGIS die Daten in der Nähe der Kanten der Verbundkacheln während des Mosaikierens tatsächlich verarbeitet. Ich fürchte, mehr kann ich Ihnen nicht sagen, ohne den Quellcode von Arc untersuchen zu können, was nie passieren wird. Sehen Sie sich jedoch die folgenden Bilder an, das resultierende Plankrümmungsbild, das aus dem Mosaik Ihrer Kacheln in einem anderen Open-Source-GIS-Softwarepaket abgeleitet wurde:

Die Kanten wurden mit einem überlagerten Vektor hervorgehoben, um die Kantenpositionen anzuzeigen (einige der Punkte, die ich zuvor verwendet habe, sind auch darin enthalten). Jetzt ohne diese Kantenbegrenzungsvektoren überlagert:

Beachten Sie, dass die Kantenartefakte, die Sie in Ihrer Frage hervorheben, in diesem Mosaik, das auch keine Mischtechniken anwendet, nicht so offensichtlich sind. Es handelt sich um dieselben DEM-Quelldaten, der einzige Unterschied besteht im spezifischen Mosaikalgorithmus.

Zusammenfassend würde ich sagen, 1) Ihre Kacheln sind nicht nahtlos, also sollten Sie ohne einen ausgefeilteren Mosaikierungsalgorithmus mit Auslaufen einige Kanteneffekte erwarten, und 2) zu einem gewissen Grad ist das Ausmaß der Artefakte, die Sie beobachten, das Ergebnis wie der ArcGIS-Mosaikierungsalgorithmus die Kanten der zusammengesetzten Kacheln behandelt. Ich hoffe das hilft.


Schlagen Sie vor, dass Sie versuchen, ein virtuelles Raster aus den DEM-TIFFs zu erstellen: Raster -> Verschiedenes -> Virtuelles Raster erstellen (Katalog). Es scheint mit Ihren Beispieldaten zu funktionieren. Entschuldigung, ich musste das Ergebnis nur ein wenig schattieren. N.