Radiación solar de puntos (Spatial Analyst)

Disponible con una licencia de Spatial Analyst.

Resumen

Calcula la radiación solar entrante para ubicaciones específicas en una tabla de ubicación o clase de entidad de puntos.

Más información acerca del cálculo de la radiación solar

Uso

  • Las ubicaciones de entrada pueden ser una clase de entidad de puntos o una tabla de coordenadas de puntos. La tabla puede ser una tabla de geodatabase, un archivo .dbf, una tabla INFO o un archivo de tabla de texto. Los valores pueden ser de tipo entero largo, flotante o doble.

  • Cuando se introducen ubicaciones por tabla, se debe especificar una lista de ubicaciones con una coordenada x,y. Si usa un archivo de coordenadas ASCII, cada línea debe contener un par x,y separado por una coma, un espacio o un tabulador. A continuación se muestra un ejemplo delimitado por espacios:

    X Y
    325541.218750 4314768.5
    325169.250000 4313907.0
    325874.031250 4313134.0
    325825.093750 4314181.5

    Como alternativa, puede especificar la pendiente (grados) y la orientación en la tabla de ubicación. Junto con la coordenada x,y, el archivo debe contener el valor de pendiente y de orientación para cada ubicación en cualquier orden. A continuación se muestra un ejemplo delimitado por comas:

    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
  • Para configuraciones de tiempo de varios días, el rango máximo de días es un total de un año (365 días o 366 días para los años bisiestos). Si el día de inicio es mayor que el día de fin, los cálculos de tiempo continuarán en el año siguiente.

    Por ejemplo, [día de inicio, día de fin] = [365, 31] representa del 31 de diciembre al 31 de enero del año siguiente. Para un ejemplo [1, 2], el tiempo se incluye para el primer día desde las 0:00 horas (1 de enero) a las 0:00 (2 de enero). El día de inicio y el día de fin no pueden ser iguales.

  • El valor del año para la configuración de la hora se utiliza para determinar un año bisiesto. No tiene ninguna otra influencia en el análisis de radiación solar ya que los cálculos son una función del período de tiempo determinado por días julianos.

  • Para las configuraciones de tiempo dentro de un día, el rango máximo de tiempo es un día (24 horas). Los cálculos no se realizarán entre días (por ejemplo, desde las 12:00 p.m. a las 12:00 p.m. del día siguiente). La hora de inicio debe ser menor a la hora de finalización.

  • El uso de un factor z es esencial para corregir los cálculos cuando las unidades de la superficie z se expresan en unidades diferentes de las unidades x,y de terreno. Para obtener resultados exactos, las unidades z deben ser las mismas que las unidades x,y de terreno. Si las unidades no son las mismas, utilice un factor z para convertir las unidades z en unidades x,y. Por ejemplo, si las unidades x,y son metros y las unidades z son pies, podría especificar un factor z de 0,3048 para convertir los pies a metros.

  • Se recomienda tener los datos en un sistema de coordenadas proyectadas con unidades de metros. Si elige ejecutar el análisis con un sistema de coordenadas esféricas, necesitará especificar un factor z apropiado para esa latitud. A continuación se muestra una lista de algunos de los factores z apropiados para usar si las unidades x,y son grados decimales y las unidades z son metros:

        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
  • El desplazamiento de altura sólo se debe especificar en metros.

  • La latitud para el área del sitio (unidades: grados decimales, positiva para el hemisferio norte y negativa para el hemisferio sur) se utiliza en cálculos como la declinación solar y la posición solar. Debido a que el análisis solar se ha diseñado para escalas de paisaje y escalas locales, es aceptable utilizar un valor de latitud para todo el DEM. Para regiones geográficas más amplias, es necesario dividir el área de estudio en zonas con diferentes latitudes.

  • Para rásteres de superficie de entrada que contengan una referencia espacial, el valor medio de la latitud se calcula de manera automática; en caso contrario, la latitud estará a 45 grados de forma predeterminada. Al utilizar una capa de entrada, se utiliza la referencia espacial del marco de datos.

  • El tamaño del cielo es la resolución de los rásteres de la cuenca visual, el mapa celeste y el mapa solar que se utilizan en los cálculos de la radiación (unidades: celdas por lado). Estas son representaciones del cielo del ráster hemisféricas, de visión ascendente, y no tienen un sistema de coordenadas geográficas. Estos rásteres son cuadrados (mismo número de filas que de columnas).

    Al aumentar el tamaño del cielo aumenta la exactitud del cálculo, pero también aumenta considerablemente el tiempo de cálculo.

  • Cuando la configuración del intervalo de día es pequeña (por ejemplo, < 14 días), se debería utilizar un tamaño de cielo más grande. Durante el análisis, se utiliza el mapa del sol (determinado por el tamaño del cielo) para representar las posiciones del sol (recorridos) para períodos de tiempo particulares para calcular la radiación directa. Con intervalos de día más pequeños, si la resolución del tamaño del cielo no es suficientemente grande, los recorridos del sol podrían superponerse, dando como resultado valores de radiación cero o inferiores para ese recorrido. Aumentar la resolución proporciona un resultado más exacto.

  • El valor máximo de tamaño del cielo es 10.000. Un valor de 200 es el predeterminado y es suficiente para DEM completos con intervalos de día grandes (por ejemplo, > 14 días). Un valor de tamaño del cielo de 512 es suficiente para los cálculos en las ubicaciones de puntos donde el tiempo de cálculo es menos problemático. Con intervalos de día menores (por ejemplo, < 14 días), se recomienda utilizar valores más grandes. Por ejemplo, para calcular la insolación para una ubicación en el ecuador con un intervalo de día = 1, se recomienda utilizar un tamaño de cielo de 2.800 o superior.

  • Se recomiendan los intervalos de días mayores que 3, dado que los recorridos del sol dentro de tres días habitualmente se superponen, dependiendo del tamaño del cielo y de la época del año. Para cálculos de todo el año con intervalos mensuales, el intervalo de día está deshabilitado y el programa utiliza intervalos mensuales de calendario a nivel interno. El valor predeterminado es 14.

  • Debido a que el cálculo de la cuenca visual puede ser altamente intenso, los ángulos del horizonte solo se trazan para el número de direcciones de cálculo especificadas. Los valores válidos deben ser múltiplos de 8 (8, 16, 24, 32, etc.) Generalmente, un valor de 8 o 16 es adecuado para las áreas con una topografía suave, mientras que un valor de 32 es el adecuado para una topografía compleja. El valor predeterminado es 32.

  • El número de direcciones de cálculo necesario está relacionado con la resolución del DEM de entrada. El terreno natural con una resolución de 30 metros es generalmente bastante suave, de forma que es suficiente un menor número de direcciones para la mayoría de las situaciones (16 o 32). Con DEM más finos y, en particular, con estructuras creadas por el hombre incorporadas en los DEM, es necesario aumentar el número de direcciones. Incrementar el número de direcciones aumentará la exactitud pero también aumentará el tiempo de cálculo.

  • La casilla de verificación Crear salida para cada intervalo proporciona la flexibilidad para calcular la insolación integrada durante un período de tiempo especificado o la insolación para cada intervalo en una serie de tiempo. Por ejemplo, para el período de tiempo dentro del día con un intervalo de hora de uno, si se marca esta casilla se crearán valores de insolación cada hora; en caso contrario, se calculará la insolación integrada para todo el día.

  • El parámetro Crear salida para cada intervalo afecta al número de atributos para las entidades de salida. Cuando se marca para el análisis de radiación de puntos, la clase de entidad de salida incluye atributos adicionales (t0, t1, t2 y así sucesivamente) que indican los valores de radiación o duración para cada intervalo de tiempo (intervalo de hora cuando la configuración de tiempo es menor a un día o intervalo de día cuando son varios días).

  • La cantidad de radiación solar que recibe la superficie es solo una parte de lo que se podría recibir desde fuera de la atmósfera. La transmisividad es una propiedad de la atmósfera que se expresa como la relación de la energía (promedio de todas las longitudes de onda) que llega a la superficie terrestre con respecto a la que se recibe en el límite superior de la atmósfera (extraterrestre). Los valores varían de 0 (sin transmisión) a 1 (transmisión completa). Los valores que se observan habitualmente son 0,6 o 0,7 para condiciones de cielo muy claro y 0,5 solo para un cielo generalmente claro.

    El valor de la energía recibida en la superficie terrestre se encuentra en la ruta más corta a través de la atmósfera (es decir, el sol se encuentra en el cénit, o directamente encima) y para el nivel del mar. Para las áreas más allá del Trópico de Capricornio y el Trópico de Cáncer, el sol nunca se puede encontrar exactamente en el cénit, incluso tampoco a mediodía. Sin embargo, este valor se sigue refiriendo al momento en el que el sol se encuentra en el cénit. Debido a las correcciones de los algoritmos a efectos de la elevación, siempre se debería dar la transmisividad a nivel del mar.

    La transmisividad tiene una relación inversa con el parámetro de proporción difusa.

  • Consulte Entornos de análisis y Spatial Analyst para obtener detalles adicionales sobre los entornos de geoprocesamiento que se aplican a esta herramienta.

Sintaxis

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})
ParámetroExplicaciónTipo de datos
in_surface_raster

Ráster de superficie de elevación de entrada.

Raster Layer
in_points_feature_or_table

La clase o la tabla de entidad de puntos de entrada especifican las ubicaciones para analizar la radiación solar.

Feature Layer; Table View
out_global_radiation_features

La clase de entidad de salida que representa la radiación global o la cantidad de insolación solar entrante (directa + difusa) que se calcula para cada ubicación.

La salida tiene unidades de vatios hora por metro cuadrado (WH/m2).

Feature Class
height_offset
(Opcional)

La altura (en metros) sobre la superficie del DEM para la que se llevan a cabo los cálculos.

El desplazamiento de altura se aplicará a todas las ubicaciones de entrada.

Double
latitude
(Opcional)

La latitud para el área del sitio. Las unidades son grados decimales, con valores positivos para el hemisferio norte y negativos para el sur.

Para rásteres de superficie de entrada que contengan una referencia espacial, el valor medio de la latitud se calcula de manera automática; en caso contrario, la latitud estará a 45 grados de forma predeterminada.

Double
sky_size
(Opcional)

La resolución o tamaño del cielo para los rásteres de cuenca visual, mapa del cielo y mapa del sol. Las unidades son celdas.

El valor predeterminado crea un ráster de 200 por 200 celdas.

Long
time_configuration
(Opcional)

Especifica la configuración de la hora (período) utilizada para calcular la radiación solar.

Los objetos de la clase Time se utilizan para especificar la configuración de la hora.

Los diferentes tipos de configuraciones de hora disponibles son TimeWithinDay, TimeMultipleDays, TimeSpecialDays y TimeWholeYear.

Las formas son las siguientes:

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

La configuración de tiempo predeterminada es TimeMultipleDays con startDay de 5 y endDay de 160 para el año juliano actual.

Time configuration
day_interval
(Opcional)

El intervalo de tiempo a lo largo de los años (unidades: días) utilizado para el cálculo de los sectores de cielo para el mapa del sol.

El valor predeterminado es 14 (bisemanal).

Long
hour_interval
(Opcional)

El intervalo de tiempo a lo largo del día (unidades: horas) utilizado para el cálculo de los sectores de cielo para mapas del sol.

El valor predeterminado es 0,5.

Double
each_interval
(Opcional)

Especifica cuándo calcular un valor de insolación total único para todas las ubicaciones o diversos valores para el intervalo de día y hora especificado.

  • NOINTERVALSe calculará un valor de radiación total único para toda la configuración de tiempo. Esta es la opción predeterminada.
  • INTERVALSe calcularán diversos valores de radiación para cada intervalo de tiempo a lo largo de toda la configuración de tiempo. El número de resultados dependerá del intervalo de hora o de día. Por ejemplo, para un año completo con intervalos mensuales, el resultado contendrá 12 valores de radiación de salida para cada ubicación.
Boolean
z_factor
(Opcional)

El número de unidades x, y de suelo en una superficie de unidades z.

El factor z ajusta las unidades de medida para las unidades z cuando son diferentes de las unidades x, y de la superficie de entrada. Los valores z de la superficie de entrada se multiplican por el factor z al calcular la superficie de salida final.

Si las unidades z y las unidades x,y están en las mismas unidades de medida, el factor z es 1. Esta es la opción predeterminada.

Si las unidades z y las unidades x,y están en diferentes unidades de medida, el factor z se debe establecer en el factor adecuado o los resultados serán incorrectos.

Por ejemplo, si las unidades z son pies y las unidades x, y son metros, debe utilizar un factor z de 0,3048 para convertir las unidades z de pies a metros (1 pie = 0,3048 metros).

Double
slope_aspect_input_type
(Opcional)

Cómo se deriva la información de la pendiente y de la orientación para el análisis.

  • FROM_DEM Los rásteres de pendiente y orientación se calculan a partir del ráster de superficie de entrada. Esta es la opción predeterminada.
  • FLAT_SURFACE Se utilizan valores constantes de cero para la pendiente y la orientación.
  • FROM_POINTS_TABLE Se pueden especificar los valores para la pendiente y la orientación a lo largo de las coordenadas x, y en el archivo de ubicaciones.
String
calculation_directions
(Opcional)

El número de direcciones acimutales utilizadas al calcular la cuenca visual.

Los valores válidos deben ser múltiplos de 8 (8, 16, 24, 32, etc.) El valor predeterminado es 32 direcciones, lo que es adecuado para una topografía compleja.

Long
zenith_divisions
(Opcional)

El número de divisiones utilizadas para crear sectores de cielo en el mapa de cielo.

El valor predeterminado es ocho divisiones (en relación con el cénit). Los valores deben ser mayores que cero y menores que la mitad del valor del tamaño del cielo.

Long
azimuth_divisions
(Opcional)

El número de divisiones utilizadas para crear sectores de cielo en el mapa de cielo.

El valor predeterminado es ocho divisiones (en relación con el norte). Los valores válidos deben ser múltiplos de 8. Los valores deben ser mayores que cero y menores que 160.

Long
diffuse_model_type
(Opcional)

Tipo de modelo de radiación difusa.

  • UNIFORM_SKY Modelo de difusión uniforme. La radiación difusa entrante es la misma desde todas las direcciones del cielo. Esta es la opción predeterminada.
  • STANDARD_OVERCAST_SKY Modelo de difusión nublada estándar. El flujo de radiación difusa entrante varía con el ángulo cénit.
String
diffuse_proportion
(Opcional)

La proporción del flujo de radiación normal global que es difusa. Los valores varían de 0 a 1.

Este valor se debe establecer de acuerdo con las condiciones atmosféricas. El valor predeterminado es 0,3 para condiciones de cielo generalmente claro.

Double
transmittivity
(Opcional)

La fracción de la radiación que pasa a través de la atmósfera (promediada para todas las longitudes de onda). Los valores varían de 0 (sin transmisión) a 1 (transmisión completa).

El valor predeterminado es 0,5 para un cielo generalmente claro.

Double
out_direct_radiation_features
(Opcional)

La clase de entidad de salida que representa la radiación solar entrante directa para cada ubicación.

La salida tiene unidades de vatios hora por metro cuadrado (WH/m2).

Feature Class
out_diffuse_radiation_features
(Opcional)

La clase de entidad de salida que representa la radiación solar entrante directa para cada ubicación que es difusa.

La salida tiene unidades de vatios hora por metro cuadrado (WH/m2).

Feature Class
out_direct_duration_features
(Opcional)

La clase de entidad de salida que representa la duración de radiación solar entrante directa.

La salida tiene unidades de horas.

Feature Class

Muestra de código

Ejemplo 1 de PointsSolarRadiation (ventana de Python)

El siguiente script de la ventana de Python muestra cómo utilizar esta herramienta.

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")
Ejemplo 2 de PointsSolarRadiation (secuencia de comandos independiente)

Calcula la cantidad de radiación solar entrante para ubicaciones de punto específicas.

# 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)

Información de licenciamiento

  • Basic: Requiere Spatial Analyst
  • Standard: Requiere Spatial Analyst
  • Advanced: Requiere Spatial Analyst

Temas relacionados