Sonneneinstrahlung (Punkte) (Spatial Analyst)

Mit der Spatial Analyst-Lizenz verfügbar.

Zusammenfassung

Leitet die Sonneneinstrahlung für bestimmte Positionen in einer Point-Feature-Class oder einer Positionstabelle ab.

Weitere Informationen zur Berechnung der Sonneneinstrahlung

Verwendung

  • Die Eingabepositionen können eine Point-Feature-Class oder eine Tabelle mit Punktkoordinaten sein. Die Tabelle kann eine Geodatabase-Tabelle, eine .dbf-Datei, eine INFO-Tabelle oder eine Text-Tabellendatei sein. Die Werte können den Typ "Short Integer", "Float" oder "Double" aufweisen.

  • Bei der Eingabe von Positionen nach Tabelle muss eine Liste der Positionen mit einer XY-Koordinate angegeben werden. Unter Verwendung einer ASCII-Koordinatendatei sollte jede Linie ein durch ein Komma, ein Leerzeichen oder eine Registerkarte getrenntes XY-Paar enthalten. Nachfolgend finden Sie ein Beispiel für eine Trennung durch Leerzeichen:

    X Y
    325541.218750 4314768.5
    325169.250000 4313907.0
    325874.031250 4313134.0
    325825.093750 4314181.5

    Alternativ können Sie Neigung (Grad) und Ausrichtung in der Positionstabelle angeben. Neben den XY-Koordinaten sollte die Datei für jeden Standort den Neigungs- und Ausrichtungswert in beliebiger Reihenfolge enthalten. Nachfolgend finden Sie ein Beispiel für eine Trennung durch Kommas:

    x, y, slope, aspect
    325541.218750, 4314768.5, 15.84516716, 310.2363586
    325169.250000, 4313907.0, 39.39801788,   2.03503442
    325874.031250, 4313134.0, 16.10847282, 223.8308563
    325825.093750, 4314181.5,  8.89850712, 205.2011261
  • 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 Schaltjahres 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.

  • 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.

  • Die Daten sollten in einem projizierten Koordinatensystem mit den Einheiten Meter vorliegen. Wenn Sie die Analyse mit einem sphäroidischen Koordinatensystem ausführen, müssen Sie einen geeigneten Z-Faktor für diesen Breitengrad angeben. Es folgt eine Liste einiger geeigneter Z-Faktoren, die verwendet werden können, wenn die XY-Einheiten als Dezimalgrade und die Z-Einheiten als Meter angegeben sind:

        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 Höhenversatz sollte nur in Meter angegeben werden.

  • Der Breitengrad der Standortfläche (Einheiten: Dezimalgrad, für die Nordhalbkugel positiv und für die Südhalbkugel negativ) wird u. a. zur Berechnung der Sonnenneigung 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 Eingabe-Oberflä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 (Einheiten: 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 Raster 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 0 oder negativ sind. Mit einer höheren Auflösung wird ein genaueres Ergebnis erzielt.

  • Der maximale Himmelsgrößenwert beträgt 10.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 Sonneneinstrahlung 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 Metern 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 Ausgaben für jedes Intervall erstellen bietet die Flexibilität, um die integrierte Sonneneinstrahlung über einen 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 Sonneneinstrahlungswerte erzeugt; andernfalls wird die integrierte Sonneneinstrahlung für den ganzen Tag berechnet.

  • Das Kontrollkästchen Ausgaben für jedes Intervall erstellen wirkt sich auf die Anzahl der Attribute für Ausgabe-Features aus. Wenn die Option für die Punktstrahlungsanalyse aktiviert ist, enthält die Ausgabe-Feature-Class zusätzliche Attribute (t0, t1, t2 usw.), die für jedes Zeitintervall (Stundenintervall bei einer Zeitkonfiguration von weniger als einem Tag, Tagesintervall bei mehreren Tagen) Strahlungs- bzw. Dauerwerte angeben.

  • Von der Oberfläche wird nur ein geringer Teil der Sonneneinstrahlung empfangen, die außerhalb der Atmosphäre auf die Erde einwirkt. Die Durchlässigkeit ist eine Eigenschaft der Atmosphäre. Sie wird definiert als die Energie (Durchschnitt aller Wellenlängen), die die Erdoberfläche erreicht, im Verhältnis zu der Energie, die an der oberen Grenze der Atmosphäre (extraterrestrial) ankommt. Die Werte liegen zwischen 0 (keine Übertragung) und 1 (vollständige Übertragung). Typische Werte liegen bei 0,6 oder 0,7 (sehr klare Bedingungen) bzw. 0,5 (klarer Himmel allgemein).

    Der Wert der Energie, die die Erdoberfläche erreicht, ist der kürzeste Pfad durch die Atmosphäre (d. h. mit der Sonne im Zenit oder direkt über der Oberfläche) auf Meereshöhe. In Gebieten außerhalb des südlichen und nördlichen Wendekreises kann die Sonne nie direkt im Zenit stehen, auch nicht zur Mittagszeit. Dieser Wert verweist dennoch auf den Moment, in dem die Sonne im Zenit steht. Da der Algorithmus Höheneffekte korrigiert, sollte die Durchlässigkeit immer auf Meereshöhe angegeben werden.

    Die Durchlässigkeit ist umgekehrt proportional zum diffusen Verhältnisparameter.

  • Weitere Informationen zur Geoverarbeitung von Umgebungen mit diesem Werkzeug finden Sie unter Analyseumgebungen und Spatial Analyst.

Syntax

PointsSolarRadiation(in_surface_raster, in_points_feature_or_table, out_global_radiation_features, {height_offset}, {latitude}, {sky_size}, {time_configuration}, {day_interval}, {hour_interval}, {each_interval}, {z_factor}, {slope_aspect_input_type}, {calculation_directions}, {zenith_divisions}, {azimuth_divisions}, {diffuse_model_type}, {diffuse_proportion}, {transmittivity}, {out_direct_radiation_features}, {out_diffuse_radiation_features}, {out_direct_duration_features})
ParameterErklärungDatentyp
in_surface_raster

Eingabe-Höhenoberflächen-Raster.

Raster Layer
in_points_feature_or_table

Die Eingabe-Point-Feature-Class oder die Tabelle mit den Positionen zur Analyse der Sonneneinstrahlung.

Feature Layer; Table View
out_global_radiation_features

Die Ausgabe-Feature-Class, die die globale Strahlung bzw. die Menge der Sonneneinstrahlung (direkt + diffus) darstellt, die für jede Position berechnet wird.

Ausgabeeinheiten: Wattstunden pro Quadratmeter (WH/m2).

Feature Class
height_offset
(optional)

Die Höhe (in Meter) über der DEM-Oberfläche, für die Berechnungen durchgeführt werden sollen.

Der Höhenversatz wird auf alle Eingabepositionen angewendet.

Double
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 Eingabe-Oberflä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-Karten-Raster. Einheiten: 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, TimeMultipleDays, TimeSpecialDays und TimeWholeYear.

Formate:

  • TimeWithinDay({day},{startTime},{endTime})
  • TimeMultipleDays({year},{startDay},{endDay})
  • TimeSpecialDays()
  • TimeWholeYear({year})

Die standardmäßige Zeitkonfiguration lautet TimeMultipleDays mit startDay 5 und endDay 160 (aktuelles julianisches Jahr).

Time configuration
day_interval
(optional)

Das Zeitintervall für das ganze Jahr (Einheiten: 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 (Einheiten: 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-Sonneneinstrahlungswert für alle Standorte oder mehrere Werte berechnet werden sollen.

  • NOINTERVALEin einzelner Gesamtstrahlungswert wird für die ganze Zeitkonfiguration berechnet. Dies ist die Standardeinstellung.
  • INTERVALFür jedes Zeitintervall der ganzen Zeitkonfiguration werden mehrere Strahlungswerte berechnet. Die Anzahl der Ausgaben hängt vom Stunden- oder Tagesintervall ab. Für ein ganzes Jahr mit monatlichen Intervallen umfasst das Ergebnis beispielsweise 12 Ausgabestrahlungswerte für jeden Standort.
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 Eingabe-Oberfläche unterscheiden. Die Z-Werte der Eingabe-Oberfläche werden bei der Berechnung der endgültigen Ausgabe-Oberflä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. Dies ist 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.

  • FROM_DEM Die Neigungs- und Ausrichtungs-Raster werden anhand des Eingabe-Oberflächen-Rasters berechnet. Dies ist die Standardeinstellung.
  • FLAT_SURFACE Für Neigung und Ausrichtung wird konstant der Wert 0 verwendet.
  • FROM_POINTS_TABLE Werte für Neigung und Ausrichtung können zusammen mit den XY-Koordinaten in der Positionen-Datei angegeben werden.
String
calculation_directions
(optional)

Die Anzahl der azimutalen Richtungen, die beim Berechnen des Sichtfeldes 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 0 und kleiner als die Hälfte des Himmelsgrößenwertes 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 0 und kleiner als 160 sein.

Long
diffuse_model_type
(optional)

Typ des diffusen Strahlungsmodells.

  • UNIFORM_SKY Einheitliches diffuses Modell. Die eingehende diffuse Strahlung ist aus allen Himmelsrichtungen gleich. Dies ist die Standardeinstellung.
  • STANDARD_OVERCAST_SKY Diffuses Standardstrahlungsmodell bei Bewölkung. Der eingehende diffuse Strahlungsfluss ändert sich mit dem Zenitwinkel.
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 Bedingungen.

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_features
(optional)

Die Ausgabe-Feature-Class, die die direkte Sonneneinstrahlung für jede Position darstellt.

Ausgabeeinheiten: Wattstunden pro Quadratmeter (WH/m2).

Feature Class
out_diffuse_radiation_features
(optional)

Die Ausgabe-Feature-Class, die die Sonneneinstrahlung für jede Position darstellt, die diffus ist.

Ausgabeeinheiten: Wattstunden pro Quadratmeter (WH/m2).

Feature Class
out_direct_duration_features
(optional)

Die Ausgabe-Feature-Class, die die Dauer der direkten Sonneneinstrahlung darstellt.

Ausgabeeinheiten: Stunden.

Feature Class

Codebeispiel

PointsSolarRadiation – Beispiel 1 (Python-Fenster)

Das folgende Skript veranschaulicht, wie dieses Werkzeug im Python-Fenster verwendet wird.

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
PointsSolarRadiation("elevation", "observers.shp", 
                     "c:/sapyexamples/output/outglobalrad1.shp", "", 35, 200, 
                     TimeMultipleDays(2009, 91, 212), 14, 0.5,"NOINTERVAL", 
                     1, "FROM_DEM", 32, 8, 8,"STANDARD_OVERCAST_SKY", 0.3, 0.5, 
                     "c:/sapyexamples/output/outdirectrad1.shp", 
                     "c:/sapyexamples/output/outdiffuserad1.shp", 
                     "c:/sapyexamples/output/outduration1.shp")
PointsSolarRadiation – Beispiel 2 (eigenständiges Skript)

Berechnen Sie die Menge der Sonneneinstrahlung für bestimmte Punktpositionen.

# PointsSolarRadiation_Example02.py
# Description: For all point locations, calculates total global, direct,
#    diffuse and direct duration solar radiation for a whole year.
# Requirements: Spatial Analyst Extension

# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "C:/sapyexamples/data"

# Set local variables
inRaster = "elevation"
inPntFC = "observers.shp"
outFeatures = "c:/sapyexamples/output/outglobal1.shp"
latitude = 35.75
skySize = 200
timeConfig = TimeMultipleDays(2009, 91, 212)
dayInterval = 14
hourInterval = 0.5
zFactor = 0.3048
calcDirections = 32
zenithDivisions = 8
azimuthDivisions = 8
diffuseProp = 0.3
transmittivity = 0.5
outDirectRad = "C:/sapyexamples/output/outdirectrad1.shp"
outDiffuseRad = "C:/sapyexamples/output/outdiffuserad1.shp"
outDirectDur = "C:/sapyexamples/output/outduration1.shp"

# Execute PointsSolarRadiation...
PointsSolarRadiation(inRaster, inPntFC, outFeatures, "", latitude, skySize, 
                     timeConfig, dayInterval, hourInterval, "INTERVAL", 
                     zFactor, "FROM_DEM", calcDirections, zenithDivisions, 
                     azimuthDivisions,"STANDARD_OVERCAST_SKY", diffuseProp, 
                     transmittivity, outDirectRad, outDiffuseRad, outDirectDur)

Lizenzinformationen

  • Basic: Erfordert Spatial Analyst
  • Standard: Erfordert Spatial Analyst
  • Advanced: Erfordert Spatial Analyst

Verwandte Themen