Sun position calculator (SPC) for Landsat imagery with geodetic latitudes

Landsat imagery comes with sun position information such as azimuth and sun elevation, but they are available only at the center of a scene. To aid in the use of Landsat imagery for various solar radiation applications such as topographic correction, solar power, urban heat island, agriculture, climate and vegetation, it is necessary to calculate the sun position information at every pixel. This research developed a PC application that creates sun position data layers in ArcGIS at every pixel in a Landsat scene. The SPC program is composed of two major routines - converting universal transverse Mercator (UTM) projection coordinates to geographic longitudes and latitudes, and calculating sun position information based on the Meeus' routine. For the latter, an innovative method was also implemented to account for the Earth's flattening on an ellipsoid. The Meeus routine implemented in this research showed about 0.2′ of mean absolute difference from the National Renewable Energy Laboratory (NREL) Solar Position Algorithm (SPA) routine when solar zenith and azimuth angles were tested with every 30 min data at four city locations (Fairbanks, Atlanta, Sydney and Rio Grande) on June 30, 2014. The Meeus routine was about ten times faster than the SPA routine. Professionals who need the Sun's position information for Landsat imagery will benefit from the SPC application.