Area Solar Radiation (Spatial Analyst)
Zusammenfassung
Leitet die Sonneneinstrahlung von einer Raster-Oberfläche ab.
Verwendung
-
Da Sonneneinstrahlungsberechnungen zeitaufwendig sein können, ist es wichtig, sicherzustellen, dass alle Parameter korrekt sind. Die Berechnung für ein großes digitales Höhenmodell (DEM) kann Stunden dauern, bei einem sehr großen DEM sogar Tage.
-
Die Ausgabestrahlungs-Raster weisen immer den Typ "Gleitkomma" auf und verwenden als Einheit Wattstunden pro Quadratmeter (WH/m2). Die Raster-Ausgabe für die direkte Dauer ist ganzzahlig, mit Stunden als Einheit.
-
Der Breitengrad der Standortfläche (Einheit: Dezimalgrad, für die Nordhalbkugel positiv und für die Südhalbkugel negativ) wird u. a. zur Berechnung der Sonnendeklination und -position verwendet.
Die Analyse wurde nur für lokale Querformatmaßstäbe entwickelt, weshalb die Verwendung eines einzigen Breitengradwertes für das ganze DEM im Allgemeinen akzeptabel ist. Bei größeren Datasets (d. h. Staaten, Ländern oder Kontinenten) unterscheiden sich die Einstrahlungsergebnisse je nach Breitengrad (Unterschied von mehr als 1 Grad) erheblich. Um umfangreichere geographische Regionen zu analysieren, muss das Untersuchungsgebiet in Zonen mit verschiedenen Breitengraden unterteilt werden.
-
Für mehrtägige Zeitkonfigurationen beträgt der maximale Zeitraum ein Jahr (365 Tage bzw. 366 Tage bei Schaltjahren). Wenn der Starttag größer als der Endtag ist, werden die Zeitberechnungen im folgenden Jahr fortgesetzt.
[start day, end day] = [365, 31] steht beispielsweise für den Zeitraum vom 31. Dezember bis zum 31. Januar des Folgejahres. Im Beispiel [1, 2] beginnt der Zeitraum am ersten Tag um 0:00 Uhr (1. Januar) und endet am 2. Januar um 0:00 Uhr (inklusiv). Starttag und Endtag können nicht gleich sein.
Der Jahreswert für die Zeitkonfiguration wird zur Bestimmung eines Schaltjahrs verwendet. Er hat keine weiteren Auswirkungen auf die Analyse der Sonneneinstrahlung, da die Berechnung eine Funktion des Zeitraums ist, der von den julianischen Tagen bestimmt wird.
-
Bei Zeitkonfigurationen innerhalb eines Tages beträgt der maximale Zeitraum einen Tag (24 Stunden). Berechnungen werden nicht über den Tageswechsel durchgeführt (z. B. von 12:00 Uhr mittags des einen bis 12:00 mittags des nächsten Tages). Die Startzeit muss kleiner als die Endzeit sein.
-
Bei Zeitkonfigurationen innerhalb eines Tages werden die Start- und Endzeiten als Sonnenzeit (Einheit: Dezimalstunden) angezeigt. Verwenden Sie das Dialogfeld zur Zeitkonvertierung, um die lokale Standardzeit und lokale Sonnenzeit (HMS) zu konvertieren. Beim Konvertieren der lokalen Standardzeit in Sonnenzeit berücksichtigt das Programm die Zeitgleichung.
-
Der Z-Faktor ist für die Korrektur von Berechnungen wichtig, wenn die Z-Oberflächeneinheiten andere Maßeinheiten als die XY-Geländeeinheiten aufweisen. Um präzise Ergebnisse zu erzielen, müssen die Z-Einheiten den XY-Geländeeinheiten entsprechen. Sind die Einheiten nicht identisch, rechnen Sie die Z-Einheiten mit einem Z-Faktor in XY-Einheiten um. Beispiel: Wenn die Z-Einheiten in Fuß und die XY-Einheiten in Metern angegeben sind, wählen Sie den Z-Faktor "0,3048", um die Z-Einheiten von Fuß in Meter umzurechnen.
-
Es wird empfohlen, die Daten in einem Projektionskoordinatensystem (Einheit: Meter) vorliegen zu haben. Wenn Sie die Analyse allerdings mit einem sphäroidischen Koordinatensystem ausführen, müssen Sie einen geeigneten Z-Faktor für diesen Breitengrad angeben. Wenn die XY-Einheiten in Dezimalgrad und die Z-Einheiten in Meter angegeben werden, können Sie zum Beispiel die folgenden Z-Faktoren für bestimmte Breitengrade verwenden:
Latitude Z-factor 0 0.00000898 10 0.00000912 20 0.00000956 30 0.00001036 40 0.00001171 50 0.00001395 60 0.00001792 70 0.00002619 80 0.00005156 -
Der Breitengrad der Standortfläche (Einheit: Dezimalgrad, für die Nordhalbkugel positiv und für die Südhalbkugel negativ) wird u. a. zur Berechnung der Sonnendeklination und -position verwendet. Da die Sonnenanalyse für Querformatmaßstäbe und lokale Maßstäbe entwickelt wurde, ist die Verwendung eines einzigen Breitengradwertes für das ganze DEM akzeptabel. Für umfangreichere geographische Regionen muss das Untersuchungsgebiet in Zonen mit verschiedenen Breitengraden unterteilt werden.
-
Für Eingabeoberflächen-Raster, die einen Raumbezug enthalten, wird der mittlere Breitengrad automatisch berechnet; andernfalls wird der Breitengrad standardmäßig auf 45 Grad festgelegt. Bei Verwendung eines Eingabe-Layers wird der Raumbezug des Datenrahmens verwendet.
-
Die Himmelsgröße ist die Auflösung der Sichtfeld-, Himmelskarten- und Sonnenkarten-Raster, die in den Strahlungsberechnungen (Einheit: Zellen pro Seite) verwendet werden. Dabei handelt es sich um nach oben hin offene, halbkugelförmige Raster-Repräsentationen des Himmels ohne geographisches Koordinatensystem. Diese Gitter sind quadratisch (gleiche Anzahl von Zeilen und Spalten).
Durch Vergrößern der Himmelsgröße vergrößert sich die Berechnungsgenauigkeit, aber auch die Berechnungszeit nimmt beachtlich zu.
-
Bei einer kleinen "Tagesintervall"-Einstellung (z. B. < 14 Tage) sollte eine größere Himmelsgröße verwendet werden. Während der Analyse wird die Sonnenkarte (entsprechend der Himmelsgröße) verwendet, um Sonnenpositionen (Spuren) für bestimmte Zeiträume darzustellen und die direkte Strahlung zu berechnen. Bei geringen Tagesintervallen und einer zu kleinen Himmelsgrößenauflösung überschneiden sich die Sonnenspuren möglicherweise, was für die betreffende Spur zu Strahlungswerten führt, die null oder negativ sind. Mit einer höheren Auflösung wird ein genaueres Ergebnis erzielt.
-
Der maximale Himmelsgrößenwert beträgt 4.000. Der Standardwert 200 reicht für vollständige DEMs mit großen Tagesintervallen (z. B. > 14 Tage) aus. Ein Himmelsgrößenwert von 512 ist ausreichend für Berechnungen an Punktpositionen, bei denen die Berechnungszeit weniger problematisch ist. Bei kleineren Tagesintervallen (z. B. < 14 Tage) wird empfohlen, höhere Werte zu verwenden. Um beispielsweise die Insolation für eine Position am Äquator mit Tagesintervall = 1 zu berechnen, wird empfohlen, eine Himmelsgröße von 2.800 oder mehr zu verwenden.
-
Tagesintervalle über 3 werden empfohlen, da sich die Sonnenspuren innerhalb von drei Tagen in der Regel je nach Himmelsgröße und Jahreszeit überschneiden. Für Berechnungen über das ganze Jahr mit monatlichem Intervall wird das Tagesintervall deaktiviert; das Programm verwendet dann intern Kalendermonatsintervalle. Der Standardwert ist 14.
-
Da die Sichtfeldberechnung sehr ressourcenintensiv sein kann, werden Horizontwinkel nur für die angegebene Anzahl von Berechnungsrichtungen aufgezeichnet. Gültige Werte müssen ein Vielfaches von 8 (8, 16, 24, 32 usw.) sein. In der Regel ist ein Wert von 8 oder 16 für Flächen mit sanfter Topografie geeignet, wohingegen der Wert 32 für komplexe Topografie angemessen ist. Der Standardwert ist 32.
-
Die Anzahl der benötigten Berechnungsrichtungen ist mit der Auflösung des Eingabe-DEMs verknüpft. Natürliches Terrain mit einer Auflösung von 30 m ist in der Regel relativ glatt, sodass in den meisten Situationen weniger Richtungen (16 oder 32) ausreichend sind. Bei feineren DEMs und insbesondere bei künstlichen Strukturen, die in die DEMs integriert wurden, muss die Anzahl der Richtungen erhöht werden. Mit zunehmender Zahl der Richtungen steigt die Genauigkeit, aber auch die Berechnungszeit.
-
Das Kontrollkästchen Create outputs for each interval bietet die Flexibilität, um die integrierte Insolation über einem angegebenen Zeitraum oder "für jedes Intervall" in einer Zeitserie zu berechnen. Wenn beispielsweise für den Zeitraum innerhalb eines Tages mit einem Stundenintervall von 1 dieses Kontrollkästchen aktiviert wird, werden stündliche Insolationswerte erzeugt; andernfalls wird die integrierte Insolation für den ganzen Tag berechnet.
-
Der Kontrollkästchenparameter For each interval wirkt sich auf das Format und die Anzahl von Ausgabestrahlungsdateien aus. Überprüfen Sie bei der Gebietsanalyse immer, dass genügend Speicherplatz ist verfügbar ist, bevor Sie Berechnungen starten, da dabei mehrere Ausgaben erstellt werden. Wenn dieses Kontrollkästchen aktiviert ist, ist das Standardausgabeformat ein Esri GRID Stack, der mehrere Bänder enthält, die den Strahlungs- oder Dauerwerten für die einzelnen Zeitintervalle entsprechen (Stundenintervall bei einer Zeitkonfiguration von weniger als einem Tag, Tagesintervall bei mehreren Tagen).
-
Der diffuse Anteil ist der Bruchteil des globalen normalen Strahlungsflusses, der diffus ist. Die Werte liegen zwischen 0 und 1. Dieser Wert sollte entsprechend den atmosphärischen Bedingungen festgelegt werden. Typische Werte liegen bei 0,2 bei sehr klaren Himmelsbedingungen bzw. 0,3 bei allgemein klaren Himmelsbedingungen.
-
Von der Oberfläche wird nur ein geringer Teil der Sonneneinstrahlung empfangen, die außerhalb der Atmosphäre auf die Erde einwirkt. Die Transmittivity ist eine Eigenschaft der Atmosphäre. Sie ist definiert als der Durchschnitt der Energie aller Wellenlängen, die an der Oberkante der Atmosphäre auftrifft, im Verhältnis zu der Energie, die die Erdoberfläche auf dem kürzesten Weg (in Zenitrichtung) erreicht. Die Werte liegen zwischen 0 (keine Übertragung) und 1 (vollständige Übertragung). Typische Werte liegen bei 0,6 oder 0,7 (sehr klare Himmelsbedingungen) bzw. 0,5 (klarer Himmel allgemein).
Da die Algorithmen Höheneffekte korrigieren, sollte die Transmittivity immer auf Meereshöhe angegeben werden. Die Transmittivity ist umgekehrt proportional zum diffusen Verhältnisparameter.
Syntax
| Parameter | Erläuterung | Datentyp |
in_surface_raster |
Eingabe-Höhenoberflächen-Raster. | Raster Layer |
latitude (optional) |
Der Breitengrad der Standortfläche. Die Angabe erfolgt in Dezimalgrad (für die Nordhalbkugel positiv und für die Südhalbkugel negativ). Für Eingabeoberflächen-Raster, die einen Raumbezug enthalten, wird der mittlere Breitengrad automatisch berechnet; andernfalls wird der Breitengrad standardmäßig auf 45 Grad festgelegt. | Double |
sky_size (optional) |
Die Auflösung oder Himmelsgröße für das Sichtfeld, die Himmelskarte und Sonnen-Kartengitternetze. Einheit: Zellen. Standardmäßig wird ein Raster von 200 x 200 Zellen erzeugt. | Long |
time_configuration (optional) |
Gibt die Zeitkonfiguration (den Zeitraum) an, die zum Berechnen der Sonneneinstrahlung verwendet wird. Die Time-Class-Objekte dienen zur Angabe der Zeitkonfiguration. Die folgenden Zeitkonfigurationsarten sind verfügbar: TimeWithinDay, TimeMultiDays, TimeSpecialDays und TimeWholeYear. Formate:
Die standardmäßige Zeitkonfiguration lautet TimeMultiDays mit start_day 5 und end_day 160 (aktuelles julianisches Jahr). | Time configuration |
day_interval (optional) |
Das Zeitintervall für das ganze Jahr (Einheit: Tage), das zur Berechnung von Himmelssektoren für die Sonnenkarte verwendet wird. Der Standardwert ist 14 (zweiwöchentlich). | Long |
hour_interval (optional) |
Das Zeitintervall für den ganzen Tag (Einheit: Stunden), das zur Berechnung von Himmelssektoren für Sonnenkarten verwendet wird. Der Standardwert ist 0.5. | Double |
each_interval (optional) |
Gibt an, ob für das angegebene Stunden- und Tagesintervall ein einziger Gesamt-Insolationswert für alle Standorte oder mehrere Werte berechnet werden sollen.
| Boolean |
z_factor (optional) |
Die Anzahl der XY-Geländeeinheiten in einer Z-Oberflächeneinheit. Durch den Z-Faktor werden die Maßeinheiten der Z-Einheiten angepasst, falls sie sich von den XY-Einheiten der Eingabeoberfläche unterscheiden. Die Z-Werte der Eingabeoberfläche werden bei der Berechnung der endgültigen Ausgabeoberfläche mit dem Z-Faktor multipliziert. Falls die XY-Einheiten und die Z-Einheiten in denselben Maßeinheiten ausgedrückt sind, lautet der Z-Faktor 1. Hierbei handelt es sich um die Standardeinstellung. Wenn die XY- und Z-Einheiten unterschiedliche Maßeinheiten aufweisen, muss der Z-Faktor entsprechend festgelegt werden, da andernfalls falsche Ergebnisse erzielt werden. Beispiel: Wenn die Z-Einheiten in Fuß und die XY-Einheiten in Metern angegeben sind, müssen Sie den Z-Faktor "0,3048" wählen, um die Z-Einheiten von Fuß in Meter umzurechnen (1 Fuß = 0,3048 Meter). | Double |
slope_aspect_input_type (optional) |
Gibt an, wie Neigungs- und Ausrichtungsinformationen für die Analyse abgeleitet werden.
| String |
calculation_directions (optional) |
Die Anzahl der azimutalen Richtungen, die beim Berechnen des Sichtfelds verwendet werden. Gültige Werte müssen ein Vielfaches von 8 (8, 16, 24, 32 usw.) sein. Der Standardwert liegt bei 32 Richtungen, was für komplexe Topografie angemessen ist. | Long |
zenith_divisions (optional) |
Die Anzahl der Abschnitte, die zum Erstellen von Himmelssektoren in der Himmelskarte verwendet werden. Der Standard beträgt 8 Abschnitte (relativ zum Zenit). Die Werte müssen größer als null und kleiner als die Hälfte des Himmelsgrößenwerts sein. | Long |
azimuth_divisions (optional) |
Die Anzahl der Abschnitte, die zum Erstellen von Himmelssektoren in der Himmelskarte verwendet werden. Der Standard beträgt 8 Abschnitte (relativ zur nördlichen Richtung). Gültige Werte müssen ein Vielfaches von 8 sein. Außerdem müssen sie größer als null und kleiner als 160 sein. | Long |
diffuse_model_type (optional) |
Typ des diffusen Strahlungsmodells.
| String |
diffuse_proportion (optional) |
Der diffuse Anteil des globalen normalen Strahlungsflusses. Die Werte liegen zwischen 0 und 1. Dieser Wert sollte entsprechend den atmosphärischen Bedingungen festgelegt werden. Der Standardwert beträgt 0,3 bei allgemein klaren Himmelsbedingungen. | Double |
transmittivity (optional) |
Der Strahlungsanteil, der durch die Atmosphäre dringt (Durchschnitt aller Wellenlängen). Die Werte liegen zwischen 0 (keine Übertragung) und 1 (vollständige Übertragung). Der Standardwert beträgt 0,5 bei allgemein klarem Himmel. | Double |
out_direct_radiation_raster (optional) |
Das Ausgabe-Raster, das die direkte Sonneneinstrahlung für jede Position darstellt. Ausgabeeinheit: Wattstunden pro Quadratmeter (WH/m2). | Raster Dataset |
out_diffuse_radiation_raster (optional) |
Das Ausgabe-Raster, das die diffuse Sonneneinstrahlung für jede Position darstellt. Ausgabeeinheit: Wattstunden pro Quadratmeter (WH/m2). | Raster Dataset |
out_direct_duration_raster (optional) |
Das Ausgabe-Raster, das die Dauer der direkten Sonneneinstrahlung darstellt. Ausgabeeinheit: Stunden. | Raster Dataset |
Rückgabewert
| Name | Erläuterung | Datentyp |
| out_global_radiation_raster |
Das Ausgabe-Raster, das die globale Strahlung oder die Gesamtmenge der Sonneneinstrahlung (direkt + diffus) darstellt, die für jede Position der Eingabeoberfläche berechnet wird. Ausgabeeinheit: Wattstunden pro Quadratmeter (WH/m2). | Raster |
Codebeispiel
Das folgende Skript im Python-Fenster veranschaulicht die Verwendung des Werkzeugs AreaSolarRadiation.
import arcpy
from arcpy.sa import *
from arcpy import env
env.workspace = "C:/sapyexamples/data"
outGlobalRadiation = AreaSolarRadiation("dem30", "", "400", TimeMultipleDays(2008,91,152))
outGlobalRadiation.save("C:/sapyexamples/output/glob_rad")
Berechnen Sie die Menge an der Sonneneinstrahlung in einem geographischen Bereich.
# Name: AreaSolarRadiation_example02.py
# Description: Derives incoming solar radiation from a raster surface.
# Outputs a global radiation raster and optional direct, diffuse and direct duration rasters
# for a specified time period. (April to July).
#
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/output"
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Set local variables
inRaster = "C:/sapyexamples/data/solar_dem"
latitude = 35.75
skySize = 400
timeConfig = TimeMultipleDays(2008, 91, 212)
dayInterval = 14
hourInterval = 0.5
zFactor = 0.3048
calcDirections = 32
zenithDivisions = 16
azimuthDivisions = 16
diffuseProp = 0.7
transmittivity = 0.4
outDirectRad = ""
outDiffuseRad = ""
outDirectDur = Raster("C:/sapyexamples/output/dir_dur")
# Execute AreaSolarRadiation
outGlobalRad = AreaSolarRadiation(inRaster, latitude, skySize, timeConfig,
dayInterval, hourInterval, "NOINTERVAL", zFactor, "FLAT_SURFACE",
calcDirections, zenithDivisions, azimuthDivisions, "UNIFORM_SKY",
diffuseProp, transmittivity, outDirectRad, outDiffuseRad, outDirectDur)
# Save the output
outGlobalRad.save("C:/sapyexamples/output/glob_rad")