Fast Fourier-transform calculation of artificial night sky brightness maps
Abstract
Light pollution poses a growing threat to optical astronomy, in addition to its detrimental impacts on the natural environment, the intangible heritage of humankind related to the contemplation of the starry sky and, potentially, on human health. The computation of maps showing the spatial distribution of several light pollution related functions (e.g. the anthropogenic zenithal night sky brightness, or the average brightness of the celestial hemisphere) is a key tool for light pollution monitoring and control, providing the scientific rationale for the adoption of informed decisions on public lighting and astronomical site preservation. The calculation of such maps from satellite radiance data for wide regions of the planet with sub-kilometric spatial resolution often implies a huge amount of basic pixel operations, requiring in many cases extremely large computation times. In this paper we show that, using adequate geographical projections, a wide set of light pollution map calculations can be reframed in terms of two-dimensional convolutions that can be easily evaluated using conventional fast Fourier-transform (FFT) algorithms, with typical computation times smaller than 10-6 s per output pixel.
- Publication:
-
Journal of Quantitative Spectroscopy and Radiative Transfer
- Pub Date:
- January 2020
- DOI:
- 10.1016/j.jqsrt.2019.106658
- arXiv:
- arXiv:1907.02891
- Bibcode:
- 2020JQSRT.24006658B
- Keywords:
-
- Light pollution;
- Atmospheric optics;
- Photometry;
- Radiometry;
- Fourier transforms;
- Astrophysics - Instrumentation and Methods for Astrophysics;
- Electrical Engineering and Systems Science - Image and Video Processing
- E-Print:
- Author formatted text of the accepted version of the paper published in Journal of Quantitative Spectroscopy &