Funktionsweise von "Schummerung"

Das Werkzeug Schummerung berechnet die hypothetische Beleuchtung einer Oberfläche durch die Bestimmung von Beleuchtungswerten für jede Zelle in einem Raster. Dies geschieht durch Festlegen der Position einer hypothetischen Lichtquelle und Berechnen der Beleuchtungswerte für jede Zelle in Bezug zu den benachbarten Zellen. Mithilfe dieser Funktion kann die Visualisierung einer Oberfläche für die Analyse oder die grafische Anzeige stark verbessert werden, insbesondere bei Verwendung von Transparenz.

Licht und Schatten sind standardmäßig Graustufen, die mit ganzen Zahlen von 0 bis 255 verknüpft sind (zunehmend von schwarz zu weiß).

Parameter "Schummerung"

Der primäre Faktor beim Erstellen einer Schummerungskarte für einen bestimmten Standort ist die Position der Sonne am Himmel.

Azimut

Der Azimut ist die Winkelrichtung der Sonne, gemessen von Norden im Uhrzeigersinn von 0 bis 360 Grad. Ein Azimut von 90 Grad steht für Osten. Der Standardazimut ist 315 Grad (NW).

Sonnenazimut von 315 Grad
Der Standardsonnenazimut (Richtung) für Schummerung ist 315 Grad.

Höhe

Die Höhe ist der Winkel der Beleuchtungsquelle über dem Horizont. Die Einheiten werden in Grad von 0 (am Horizont) bis 90 Grad (direkt über dem Objekt) gemessen. Die Standardeinstellung ist 45 Grad.

Sonnenazimut von 45 Grad
Die Standardsonnenhöhe für Schummerung beträgt 45 Grad.

Beispiel für eine Schummerung

Im nachfolgenden Beispiel hat die Schummerung einen Azimut von 315 Grad und eine Höhe von 45 Grad:

Beispiel für eine Schummerungsausgabe
Beispiel für eine Schummerungsausgabe

Verwendung von Schummerungen für die Anzeige

Durch die Platzierung eines Höhen-Rasters auf einem Schummerungs-Raster und die Anpassung der Transparenz des Höhen-Rasters können Sie einfach eine optisch ansprechende Reliefkarte einer Landschaft erstellen. Weitere Informationen zum Anpassen der Anzeige und des Aussehens finden Sie unter Darstellung von Bilddaten.

Verwenden von Transparenz zur Kombination von Höhen-Raster mit Schummerung
Erstellen eines geschummerten Relief-Rasters durch Kombination von Höhen- und Schummerungs-Rastern mit Transparenz.

Sie können weitere Layer wie Landnutzungsarten, Vegetation, Straßen oder Wasserläufe hinzufügen, um den Informationsgehalt der Anzeige zu erhöhen.

Verwenden von Schummerungen in Analysen

Durch Modellieren des Schattens (Standardeinstellung) können Sie die lokale Beleuchtung berechnen und prüfen, ob die Zelle im Schatten liegt oder nicht.

Ferner lässt sich durch die Modellierung des Schattens jede Zelle ermitteln, die zu einer bestimmten Tageszeit im Schatten einer anderen Zelle liegt. Zellen im Schatten einer anderen Zelle werden mit 0 kodiert. Alle anderen Zellen werden mit ganzen Zahlen von 1 bis 255 kodiert. Sie können alle Werte größer 1 als 1 reklassifizieren, wodurch sich ein binäres Ausgabe-Raster ergibt. Im folgenden Beispiel befinden sich die schwarzen Bereiche im Schatten. Der Azimut ist in allen Abbildungen identisch, der Sonnenwinkel (die Höhe) wurde jedoch geändert.

Schatten mit niedrigem Sonnenwinkel
Ein niedriger Sonnenwinkel führt zu mehr Schattenflächen.
Schatten mit hohem Sonnenwinkel
Ein hoher Sonnenwinkel führt zu weniger Schattenflächen.

Berechnung der Schummerung

Um den Schummerungswert zu berechnen, werden zunächst Höhe und Azimut der Beleuchtungsquelle benötigt. Diese Werte werden mit den Berechnungen für Neigung und Ausrichtung verarbeitet, um den endgültigen Schummerungswert für jede Zelle im Ausgabe-Raster zu bestimmen.

Schummerungsalgorithmus

Der Algorithmus zur Berechnung des Schummerungswertes lautet wie folgt:

(1) Hillshade = 255.0 * ((cos(Zenith_rad) * cos(Slope_rad)) + (sin(Zenith_rad) * sin(Slope_rad) * cos(Azimuth_rad - Aspect_rad)))

Beachten Sie: Wenn aus der Berechnung des Schummerungswertes ein Ergebnis kleiner als 0 hervorgeht, ergibt sich daraus, dass der Ausgabezellwert gleich 0 ist.

Berechnen des Beleuchtungswinkels

Die Höhe der Beleuchtungsquelle wird in Graden über der Horizontalen angegeben. Für die Berechnungsformel des Schummerungswertes ist jedoch eine Darstellung des Winkels in Radianten erforderlich. Außerdem muss er der Ablenkung von der Vertikalen entsprechen. Die Richtung von der Oberfläche direkt nach oben (direkt darüber) wird als "Zenit" bezeichnet. Der Zenitwinkel wird vom Zenitpunkt in Richtung der Beleuchtungsquelle gemessen und entspricht dem 90-Grad-Komplement der Höhe. Um den Beleuchtungswinkel zu berechnen, konvertieren Sie zuerst den Höhenwinkel in den Zenitwinkel. Als Nächstes konvertieren Sie den Winkel in Radianten.

Höhenwinkel in Zenitwinkel konvertieren:

(2) Zenith_deg = 90.0 - Altitude

In Radianten konvertieren:

(3) Zenith_rad = Zenith_deg * pi / 180.0

Berechnen der Beleuchtungsrichtung

Die Richtung der Beleuchtungsquelle, Azimut, wird in Graden angegeben. Für die Schummerungsformel muss dieser Winkel in Radianteneinheiten vorliegen. Ändern Sie zunächst den Azimutwinkel von der geographischen Einheit (Kompassrichtung) in eine mathematische Einheit (rechter Winkel). Als Nächstes konvertieren Sie den Azimutwinkel in Radianten.

Azimutwinkel anpassen:

(4) Azimuth_math = 360.0 - Azimuth + 90.0

Beachten Sie, dass, wenn Azimuth_math >= 360.0 Folgendes zutrifft:

(5) Azimuth_math = Azimuth_math - 360.0

In Radianten konvertieren:

(6) Azimuth_rad = Azimuth_math * pi / 180.0

Berechnen der Neigung und Ausrichtung

Ein sich bewegendes 3x3-Fenster durchläuft jede Zelle im Eingabe-Raster. Für jede Zelle in der Mitte des Fensters werden ein Ausrichtungs- und ein Neigungswert anhand eines Algorithmus berechnet, der die Werte der acht angrenzenden Zellen mit berücksichtigt. Die Zellen werden als Buchstaben a bis i identifiziert, wobei e die Zelle darstellt, für die die Ausrichtung berechnet wird.

Die Änderungsrate in der x-Richtung für Zelle e wird mit dem folgenden Algorithmus berechnet:

(7) [dz/dx] = ((c + 2f + i) - (a + 2d + g)) / (8 * cellsize)

Die Änderungsrate in der y-Richtung für die Zelle e wird mit dem folgenden Algorithmus berechnet:

(8) [dz/dy] = ((g + 2h + i) - (a + 2b + c)) / (8 * cellsize)

Das steilste Gefälle für eine Zelle in der Oberfläche ist die Neigung. Der Algorithmus zur Berechnung der Neigung in Radianten, einschließlich Z-Faktor, lautet wie folgt:

(9) Slope_rad = ATAN (z_factor * √ ([dz/dx]2 + [dz/dy]2))

Die Richtung der steilsten Neigungsrichtung ist die Ausrichtung. Die Ausrichtung in Radianten wird im Bereich von 0 bis 2 Pi definiert, wobei 0 Osten ist. Die Ausrichtung wird anhand der Regeln mit folgendem Algorithmus ermittelt:

(10) If [dz/dx] is non-zero: Aspect_rad = atan2 ([dz/dy], -[dz/dx]) if Aspect_rad < 0 then Aspect_rad = 2 * pi + Aspect_rad If [dz/dx] is zero: if [dz/dy] > 0 then Aspect_rad = pi / 2 else if [dz/dy] < 0 then Aspect_rad = 2 * pi - pi / 2 else Aspect_rad = Aspect_rad

Beispiel für die Berechnung der Schummerung

In diesem Beispiel wird der Schummerungswert der mittleren Zelle des beweglichen Fensters berechnet.

Eingabe-Höhenraster

Die Zellengröße beträgt 5 Einheiten. Für die Höhe wird der Standardwert 45 Grad und für den Azimut 315 Grad verwendet.

  • Beleuchtungswinkel

    Die Berechnung des Zenitwinkels aus Gleichung 2 lautet wie folgt:

    (2) Zenith_deg = 90.0 - Altitude = 90.0 - 45.0 = 45.0

    Und konvertiert wie folgt in Radianten aus Gleichung 3:

    (3) Zenith_rad = Zenith_deg * pi / 180.0 = 45.0 * 3.1415926536 / 180.0 = 0.7853981634

  • Beleuchtungsrichtung

    Die Berechnung für die Umwandlung des Azimutwinkels von einem geographischen in einen mathematischen Winkel mit Gleichung 4 lautet wie folgt:

    (4) Azimuth_math = 360.0 - Azimuth + 90.0 = 360.0 - 315.0 + 90.0 = 135.0

    Die Umwandlung von Azimutwinkel in Radianten mit Gleichung 6 lautet wie folgt:

    (6) Azimuth_rad = Azimuth_math * pi / 180.0 = 135.0 * 3.1415926536 / 180 = 2.3561944902

  • Neigung und Ausrichtung

    Die Berechnung für die Änderungsrate in der x-Richtung für die mittlere Zelle e lautet wie folgt:

    (7) [dz/dx] = ((c + 2f + i) - (a + 2d + g)) / (8 * cellsize) = ((2483 + 4966 + 2477) - (2450 + 4904 + 2447)) / (8 * 5) = (9926 - 9801) / 40 = 3.125

    Die Berechnung für die Änderungsrate in der y-Richtung für die mittlere Zelle e lautet wie folgt:

    (8) [dz/dy] = ((g + 2h + i) - (a + 2b + c)) / (8 * cellsize) = (2447 + 4910 + 2477) - (2450 + 4922 + 2483) / (8 * 5) = (9834 - 9855) / 40 = -0.525

    Der Neigungswinkel wird wie folgt berechnet:

    (9) Slope_rad = ATAN ( z_factor * √ ([dz/dx]2 + [dz/dy]2)) = atan(1 * sqrt((3.125 * 3.125) + (-0.525 * -0.525))) = atan(1 * sqrt(10.04125 + 0.275625)) = atan(1 * 3.1687931457) = 1.2651101670

    Die Berechnung für den Winkel Aspect_rad von Regel 10 lautet wie folgt (da dz/dx in diesem Beispiel ungleich Null ist):

    Aspect_rad = atan2 ([dz/dy], -[dz/dx]) = atan2(-0.525, -3.125) = -2.9751469600

    Da der Wert niedriger als 0 ist, wird dieser Teil der Regel angewendet:

    Aspect_rad = 2 * pi + Aspect_rad = 2 * 3.1415926536 + -2.9751469600 = 3.3080383471

  • Schummerung

    Die endgültige Berechnung der Schummerung lautet wie folgt:

    Hillshade = 255.0 * ((cos(Zenith_rad) * cos(Slope_rad)) + (sin(Zenith_rad) * sin(Slope_rad) * cos(Azimuth_rad - Aspect_rad))) = 255.0 * ((cos(0.7857142857) * cos(1.26511)) + (sin(0.7857142857) * sin(1.26511) * cos(2.3571428571 - 3.310567))) = 153.82

Wenn das Ausgabe-Raster vom Typ "Integer" ist, lautet der Schattenwert für die mittlere Zelle e = 154.

Referenzen

Burrough, P. A., und R. A. McDonell. 1998. Principles of Geographical Information Systems. New York: Oxford University Press.