Das Werkzeug SchummerungSchummerung 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).
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.
Beispiel für eine Schummerung
Im nachfolgenden Beispiel hat die Schummerung einen Azimut von 315 Grad und eine Höhe von 45 Grad:
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.
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.
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.
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.