Points solar radiation in an hourly basis for a whole year.

JimenoFonseca

Hi there,

I have been working lately on setting the points solar radiation routine to calculate the 8760 values of solar radiation in the year for a set of points in a single run.

The routine in ArcGIS 10.1 does it already for as much as an hourly scale but on a daily basis and only delivering apparently those hours with solar irradiance (not exactly the 24 hours)

To solve this, it is necessary to run the algorithm for every day of the year + identifying the hours of daylight of the day and do the correction to obtain 8760. Is there a way to edit the source code or do it so in a simpler manner?

Thanks,

Jimeno

ShaunWalbridge

Hello Jimeno,

I'm not sure exactly what you're looking for, but you should be able to use the solar functions, perhaps in conjunction with other packages, to answer your question. I see that as you mention, when you run Points Solar Radiation on a single day, you get back the time steps but not the actual time that it corresponds to -- one way to get back the corresponding actual hours would be to use a script such as this one calculating sunrise and sunset or use a package like PyEphem to get more precise calculations.

From there, you should be able to combine the solar observations from ArcGIS with 'actual time', so you can correspond the output T0 with the sunrise time and the final Tn with the sunset time. So my psuedo-code for this on a single day would be something like:

    import arcpy
    import ephem

    julian_day = 165
    
    # run solar calculation for all sites on this day; defaults used here
    arcpy.PointsSolarRadiation("elevation", "my_sites", 
        "output_global_radiation", "", 35, 200, 
        arcpy.sa.TimeWithinDay(julian_day, 0, 24))

    # figure out sunrise and sunset. if the points are close together, do this 
    # once for all sites, if they're geographically spread over a large area, you'll need
    # to make these calculations for each site.
    
    obs = ephem.Observer()
    obs.lat = ' 34' # obtained from my_sites
    obs.long = '-120.0' # obtained from my_sites
    obs.elev = 15 # extract from elevation raster at site location

    sunrise = obs.previous_rising(ephem.Sun()) # sunrise
    noon = obs.next_transit(ephem.Sun(), start=sunrise) # solar noon
    sunset = obs.next_setting(ephem.Sun()) # sunset
   
    # then, map those times to the T0...Tn time steps, which are by default in half-hour intervals.
    # output this result to a CSV or the like.

Does that help get you started? I think something like this is easier than trying to dig up the internals to figure out the corresponding times.

cheers,
Shaun

jbromp_1-1715221057029.png