Höhensprünge ausschließen

Gespeichert von Christine am
im Forum ArcGIS
Hallo Forum,

ich habe folgende Daten vorliegen

Länge Breite Höhe

Meine Datenpunkte sind 172m voneinander entfernt und sollen mir eigentlich schöne Höhendaten liefern. Nun werden die Höheninformationen aber durch Wolken gestört und verursachen Höhensprüngen von unterschiedlichen Ausmaßen. Kurz...innerhalb von 172 m kann sich die Topographie nicht um 500m ändern (zumindest in meinem Untersuchungsgebiet nicht).

Hier nun meine Frage, wie kann ich solche Sprünge in der Höhe ausschließen? Also ungefähr so....wenn die Höhendifferenz von einem Punkt zum Nächsten > 100m ist, schließe die nachfolgenden Daten bis zur nächsten Höhendifferenz von 100m aus.

Ein Beispiel:

Länge Breite Höhe
-13,887327 -71,263701 8,523
-13,888487 -71,265205 8,517
-13,889652 -71,266709 8,979
-13,890822 -71,268213 8,482
-13,891997 -71,269716 208,455
-13,893175 -71,271218 208,666
-13,894356 -71,272718 208,729
-13,895539 -71,274215 208,784
-13,896723 -71,275711 208,430
-13,897905 -71,277208 208,554
-13,899084 -71,278707 208,548
-13,900258 -71,280208 208,605
-13,901428 -71,281711 208,737
-13,902592 -71,283215 8,544
-13,903749 -71,284720 8,641
-13,904898 -71,286227 8,685
-13,906040 -71,287734 8,802


Ich möchte dann alle Punkte mit den Höhen 208,xxx von den nachfolgenden Untersuchungen ausschließen.

Kann man sowas in ArcGIS 9.1 machen? Wenn ja, wie?

Danke schonmal....

Christine

Gespeichert von Obiwan am Mi., 30.05.2007 - 08:54

Permalink
Bin nicht sicher, ob ich das jetzt richtig verstanden habe....
Du kannst doch die Datenpunkte in ArcGIS laden. Du hast dann ein Punktshape. Die Punkte, die eine "zu große" Höhendifferenz haben kannst Du löschen!
Aus deinen Punkten kannst Du natürlich auch ein Raster erstellen und dann die nicht gewünschten Rasterzellen auf "NoData" setzen!

Habe ich das jetzt verstanden oder etwa nicht?!?!?

Gespeichert von Christine am Mi., 30.05.2007 - 09:18

Permalink
hallo obiwan,

doch doch...du hast es richtig verstanden, nur meine datei ist seeehhhhrrrr groß und ich würde diesen vorgang sehr gern automatisieren. vielleicht mit einem queri oder so. leider weiß ich nicht, wie man sagt..."suche eine höhendifferenz von über 50m(oder was auch immer) und lösche die nachfolgenden punkte bis zum nächsten sprung von 50m". geht sowas?

grüße

Gespeichert von Obiwan am Mi., 30.05.2007 - 09:45

Permalink
Nur so eine Idee, bin mir nicht sicher ob es praktikabel ist.
Aus einer Rasterkarte der Höhenangeben kannst Du Dir die Hangneigung (Spatial Analyst Slope) rechnen lassen. Wenn es hier dann große Sprünge gibt, kannst du diese mit einer konditionalen Abfrage (Spatial Analyst - Con) auf NoData setzen lassen. Allerdings bleibt das Problem, dass Du so nur die Sprünge ganz einfach und automatisch identifizieren kannst. Um die Zellen zwischen zwei Sprüngen zu identifizieren muss ich noch etwas grübeln.... :o

Gespeichert von Maxilla am Mi., 30.05.2007 - 10:18

Permalink
Mal dazwischen gefragt: Ist es nicht so, dass Werte über - sagen wir mal - 150m in deinem Untersuchungsgebiet sowieso immer als unrealistisch einzustufen sind oder kommen solche Werte auch "richtigerweise" vor? Falls nicht, könntest du diesen Schwellenwert ja einfach abfragen und die gefundenen Punkte löschen. Falls doch ... lassen wir Obi erst mal weitergrübeln ;o)

Gespeichert von Christine am Mi., 30.05.2007 - 10:33

Permalink
solche sprünge sollten nicht vorkommen, doof ist nur, wenn ich die höhen nachfolgender punkte vergleiche, dann erkenne ich zwar einen höhensprung von 8m auf 208m, jedoch keinen mehr von 208m auf 208m...aber alle daten (keiner weiß genau wie viele es sind) mit 208m sollen gelöscht werden. genau das ist nämlich mein problem....(ich hoffe ich habe es verständlich formuliert)

Gespeichert von Maxilla am Mi., 30.05.2007 - 10:58

Permalink
Ob Sprung oder nicht, wäre hier ja irrelevant. Frage wäre in diesem Falle lediglich: Kommen "realistische" Höhen über 200m in deinem Gebiet vor, oder sind das grundsätzlich Fehldaten?

Gespeichert von Christine am Mi., 30.05.2007 - 11:02

Permalink
ach so...tschuldigung...sie kommen vor...hab alles zwischen 0 und 1800m

Gespeichert von Maxilla am Mi., 30.05.2007 - 11:27

Permalink
Schade, aber das wäre ja auch zu einfach gewesen. Okay, dann was anderes. In Avenue wüsste ich, wie ich die Tabelle dahingehend abfragen würde, aber mein gebrochenes VBA reicht dafür nicht aus. Daher - solange Obi noch grübelt - noch ne andere Sache, die mir eingefallen ist, als ich mir deine Daten so vorgestellt habe. Keine Ahnung, in wie weit das praktikable Ergebnisse liefert:

- Aus deinen Punkten ein Raster erzeugen
- Das Raster "auf den Kopf stellen" (mit -1 multiplizieren)
- Deine Fehldaten werden dadurch zu "Löchern"
- Flow Direction berechnen
- Sink berechnen
- Dies sollte deine Fehldaten als abflusslose Senken ausgeben
- Das Sink-Grid als Maske zur Zuweisung von NoDatas an dein Original verwenden

Aber wie gesagt: reines Experiment

Gespeichert von Obiwan am Mi., 30.05.2007 - 11:45

Permalink
Hier läuft ja ein ganz schönes Geschnatter während ich mir (leider unnötig und ohne dolles Ergebnis) den kopf zerbreche ;)!

Hört sich (natürlich!) gut an, Mx!

Gespeichert von Christine am Mi., 30.05.2007 - 14:20

Permalink
hallo ihr beiden...

so ein geschnatter kann aber sehr helfend sein...also...es ist zwar ziemlich krass wie es sich anhört, aber es scheint zu funktionieren...zumindest habe ich abflusslose senken bekommen...jetzt muss ich noch die nodata generieren.
melde mich dann wieder...aber bis jetzt siehts gut aus....

grüße

Gespeichert von Christine am Mi., 30.05.2007 - 16:48

Permalink
so...ich habe dann jetzt ein sink-raster erstellt...nur an meinem vorhaben die nodatas zu zuweisen bin ich gescheitert. mir möchte nicht recht in den sinn kommen, wie ich das machen soll...

Gespeichert von Maxilla am Mi., 30.05.2007 - 17:49

Permalink
Im Raster Calculator:

Con(IsNull([Sink]), [hoehe])

wobei [Sink] dein Sink-Grid und [hoehe] dein Original Höhengrid sein sollte. Das Statement sagt: "Wo Sink gleich Null (NoData) ist, da weise das Höhengrid zu, ansonsten weise NoData zu."

Anschließend kannst du noch versuchen, die NoData-Bereiche - wenn sie nicht zu groß sind - über ein FocalMean aus den interpolierten Werten der Nachbarschaft zu schließen:

Con(IsNull([hoehe]), FocalMean([hoehe], rectangle, 5, 5), [hoehe])

wobei du mit der Nachbarschaft ("rectangle, 5, 5") noch experimentieren kannst.

Gespeichert von Maxilla am Do., 31.05.2007 - 08:11

Permalink
Hmm, zu dumm. Die ganze Sache hat natürlich noch ein Manko, das ich gestern gar nicht bedacht habe: Denn natürlich kann ein solches, auf den Kopf gestelltes DEM auch "abflusslose Senken" haben, die nicht auf Fehldaten beruhen (also z.B. Bergkuppen u.ä.). Diese Bereiche müsste man aus dem Sink-Grid entfernen oder von der NoData-Zuweisung ausschließen, um eine unnötige Datenverschlechterung zu vermeiden. Eventuell wäre hier die von Obi vorgeschlagene Slope-Berechnung ein Weg...

Gespeichert von Christine am Do., 31.05.2007 - 10:38

Permalink
ich werde alles mal ausprobieren....aber ich danke euch für die riesen große hilfe...
wie das mit dem bergkuppen aussieht, muss ich schauen...aber soviele hab ich da nicht ;o)....

besten dank noch mal...

christine