Hallo! Bin neu hier und muss mich gleich mit einer Frage an euch wenden:
Ich muss verschiedene Polygone mit Rasterthemen (DHM Derivate) verschneiden und mir dann für jedes Polygon einen Mittelwert berechnen lassen. Das funktioniert mit Zonal Statistics vom Spatial Analyst auch hervorragend (z.b. für Seehöhen oder Hangneigungen).
Jetzt habe ich aber das Problem, dass ich auch die Exposition für meine Polygone brauche. Hier tuts eine einfacher Mittelwert (wie ich ihn mit den Zonal Statistic bekomme) nicht mehr, da Expositionen in einem Polygon ja z.B. die Werte 2, 356, 3 haben können. Das Polygon wäre damit logischerweise nach Norden ausgerichtet. Mein Mittelwert aus diesen drei Werte wäre aber 120,3. Das ist leider Blödsinn.
Hat dieses Problem schon jemand gehabt bzw. gelöst?
Lg, roland.
Ich muss verschiedene Polygone mit Rasterthemen (DHM Derivate) verschneiden und mir dann für jedes Polygon einen Mittelwert berechnen lassen. Das funktioniert mit Zonal Statistics vom Spatial Analyst auch hervorragend (z.b. für Seehöhen oder Hangneigungen).
Jetzt habe ich aber das Problem, dass ich auch die Exposition für meine Polygone brauche. Hier tuts eine einfacher Mittelwert (wie ich ihn mit den Zonal Statistic bekomme) nicht mehr, da Expositionen in einem Polygon ja z.B. die Werte 2, 356, 3 haben können. Das Polygon wäre damit logischerweise nach Norden ausgerichtet. Mein Mittelwert aus diesen drei Werte wäre aber 120,3. Das ist leider Blödsinn.
Hat dieses Problem schon jemand gehabt bzw. gelöst?
Lg, roland.
- Anmelden oder Registieren, um Kommentare verfassen zu können
Gespeichert von Maxilla am Di., 18.05.2010 - 18:23
PermalinkAtan2 (Summe(Cosinus(Exposition)), Summe(Sinus(Exposition)))
Diese Formel muss man jetzt nur noch ins GIS übersetzen - was ich selber allerdings auch noch nicht ausprobiert habe.
Erstmal die Sinus- und Cosinusgrids errechnen. Das dürfte im Raster Calculator einfach sein:
Sin([Exposition])
Cos([Exposition])
Die beiden Grids summiert man nun über Summarize und bindet die Felder an die Polygontabelle. Jetzt müsste man nur noch die Atan2-Funktion auf die beiden Felder loslassen, aber ich wüsste nicht, wie man diese Funktion im Field Calculator umschreibt (würde mich aber sehr interessieren :)
Der Spatial-Analyst unterstützt die Atan2-Funktion allerdings. Falls deine Polygone also einem Raster entsprechen, könntest du die summierten Sin und Cos-Werte in Grids umwandeln und dann mit dem ATan2-Tool aus der ArcToolbox die mittlere Exposition berechnen.
Gespeichert von RolandG am Mi., 19.05.2010 - 09:14
PermalinkCosinus- und Sinusgrids berechnen. Zonal Statistics - Zum Polygonlayer joinen. Hier mit dem Fieldcalculator die Werte auslesen und in neuem Feld berechnen.
Hab die ATAN2 Funktion im Field Calculator so umgesetzt:
Muss aber die Ergebnisse noch überprüfen.
'========================'
x= [COS_Stat.SUM]
y = [SIN_Stat.SUM]
'========================'
If y > 0 Then
If x >= y Then
Atan2 = Atn(y / x)
ElseIf x <= -y Then
Atan2 = Atn(y / x) + Pi
Else
Atan2 = Pi / 2 - Atn(x / y)
End If
Else
If x >= -y Then
Atan2 = Atn(y / x)
ElseIf x <= y Then
Atan2 = Atn(y / x) - Pi
Else
Atan2 = -Atn(x / y) - Pi / 2
End If
End If
Gespeichert von Maxilla am Mi., 19.05.2010 - 10:07
PermalinkGespeichert von RolandG am Do., 27.05.2010 - 14:44
Permalink