Ich habe ein Höhenraster und möchte einige grundlegende Geländevariablen wie Neigung, Robustheit, Aspekt oder Krümmung erhalten.
Das Problem ist, dass ich anscheinend nicht die Möglichkeit habe, die Krümmung mit dem Terrain-Analysten zu berechnen:
Kann mir jemand sagen, warum dieser Indikator nicht in der Liste ist? Wie kann ich es anders berechnen?
Antworten:
Die Krümmung ist eine komplexe Geländeableitung, die zu berechnen ist. Die von Ihnen verwendete Gleichung hängt von der Auflösung Ihrer Eingabedaten ab, da Sie sicherstellen müssen, dass die von Ihnen berechneten Krümmungsergebnisse von Rauschen in den Daten unterschieden werden können.
In letzter Zeit wurden viele Untersuchungen zu Krümmungsberechnungen an hochauflösenden LiDAR-Daten durchgeführt, die zeigten, dass eine Skalierungsunterbrechung bei einer Auflösung von etwa 2 bis 3 Metern vorliegt und ab diesem Punkt mehr verschiedene Algorithmen (mit denen ich nicht so vertraut bin) verwendet werden müssen . Die besten Informationen zur Berechnung der topografischen Krümmung stammen wahrscheinlich von Hurst et al. 2012 und den darin enthaltenen Referenzen.
Das Grundprinzip der Krümmungsberechnung besteht, wie bei Steigung und Aspekt, darin, ein sich bewegendes Fenster über die Höhenoberfläche zu führen und die Höhenwerte an eine 6-Term-Polynomfunktion anzupassen, deren Koeffizienten die Steigung, den Aspekt und die Krümmung des Zentrums ergeben Zelle des beweglichen Fensters.
ArcGIS verwendet ein 3x3-Suchfenster, das nur in Gebieten ohne Vegetation gute Ergebnisse liefert. Dies macht das Tool ziemlich unbrauchbar, es sei denn, die Benutzer kennen diese Einschränkung. Dies könnte darauf hindeuten, dass es in QGIS nicht vorhanden ist.
Die Mathematik stammt ursprünglich (glaube ich) aus Evans (1980) und wurde auf einigen Seiten in Prinzipien geographischer Informationssysteme (Amazon Link) vereinfacht, die ich als gute Anleitung für diese Art der Geländeanalyse auf einer grundlegenden Ebene empfehlen kann.
Eine Möglichkeit, die Krümmung eines DEM zu berechnen, besteht darin, das DEM in ein ASCII-Raster zu konvertieren, es in ein Numpy-Array zu lesen und dann die Polynomanpassung an einem sich bewegenden Fenster durchzuführen, das durch die Daten verläuft. Dies ist recht einfach, aber sehr langsam auszuführen und bedarf einer gewissen Optimierung (diese Art von Operationen werden häufig nach c ++ portiert, um sie zu beschleunigen).
Um die Operation in QGIS auszuführen, können Sie das GRASS-Plugin r.slope.aspect verwenden, das auch durch das feste 3x3-Fenster begrenzt ist.
Mir ist klar, dass dies nicht die einfache Antwort ist, auf die Sie zweifellos gehofft haben, aber ich hoffe, dass Sie verstehen, dass es komplex ist, die Krümmung auf sinnvolle Weise abzuleiten. Alles Gute.
Evans, I. S. (1980), An integrated system of terrain analysis and slope mapping, Z. Geomorphol., 36, 274–295.
Die ESRI-Version von Raster Analysis zur Berechnung der Krümmung kann hilfreich sein, um ein Plugin für QGIS zu entwickeln.
Für jede Zelle ein Polynom vierter Ordnung der Form:
Z = Ax²y² + Bx²y + Cxy² + Dx² + Ey² + Fxy + Gx + Hy + I
Wird an eine Oberfläche angepasst, die aus einem 3x3-Fenster besteht. Die Koeffizienten a, b, c usw. werden aus dieser Oberfläche berechnet.
Die Beziehungen zwischen den Koeffizienten und den neun Höhenwerten für jede Zelle, die wie im Diagramm gezeigt nummeriert ist, lauten wie folgt: Diagramm der Krümmungswerte Diagramm der Krümmungswerte
A = [(Z1 + Z3 + Z7 + Z9) / 4 - (Z2 + Z4 + Z6 + Z8) / 2 + Z5] / L4
B = [(Z1 + Z3 - Z7 - Z9) / 4 - (Z2 - Z8) / 2] / L3
C = [(-Z1 + Z3 - Z7 + Z9) / 4 + (Z4 - Z6)] / 2] / L3
D = [(Z4 + Z6) / 2 - Z5] / L2
E = [(Z2 + Z8) / 2 - Z5] / L2
F = (-Z1 + Z3 + Z7 - Z9) / 4L2
G = (-Z & sub4; + Z & sub6;) / 2L
H = (Z2 - Z8) / 2L
I = Z5
Die Ausgabe des Werkzeugs "Krümmung" ist die zweite Ableitung der Oberfläche, z. B. die Steigung der Steigung.
Krümmung = -2 (D + E) * 100
Vollständige Informationen und Quelle:
http://help.arcgis.com/de/arcgisdesktop/10.0/help/index.html#//00q90000000t000000
-2(D+E)*100als (100/L2)*(3*Z5 - [Z2+Z4+Z6+Z8+Z5])zeigt, wie dieser Wert als dreifacher Wert Z5abzüglich einer Nachbarschaftssumme Z2+Z4+Z6+Z8+Z5(unter Verwendung eines Kreises mit Radius 1) erhalten wird, wobei alle Werte um skaliert werden 100/L2. Das sind nur drei einfache Netzoperationen.
Die Krümmung kann mit dem SAGA-Modul 'Geländeanalyse - Morphometrie ---> Steigung, Aspekt, Krümmung' berechnet werden.
Die Berechnung könnte auf der Grundlage eines der folgenden Algorithmen erfolgen:
LandSerf kann dies tun. Sie können die Fenstergröße definieren (3 * 3, 5 * 5, 7 * 7, 11 * 11, ...), aber es muss eine ungerade Zahl sein. Multiskalenanalysen. Sie können die Skalierungsabhängigkeit von Neigung, Aspekt und Krümmung berücksichtigen. http://www.landserf.org/ LandSerf passt die quadratische Polynomgleichung an bestimmte vordefinierte Fenster an. Wenn Sie jedoch eine große Fenstergröße wie 50 * 50 definieren, dauert dies lange. Dies hängt von der von Ihnen definierten Raster- und Fenstergröße ab. Jo Wood schrieb Landserf für seine Doktorarbeit. Es ist in Java geschrieben.
Kann in LandSerf in verschiedenen Maßstäben berechnet werden
Sie können auch kostenlos SAGA GIS (http://sourceforge.net/apps/trac/saga-gis/wiki) oder TAS (http://www.uoguelph.ca/~hydrogeo/TAS/index.html) ausprobieren.