“…When multiple sources are at play, surrounded by extended reflecting surfaces, the resulting observed radiance can be computed by adding the individual emitter patch contributions, resorting to the whole conceptual and calculation toolbox provided by linear systems theory, as described in detail elsewhere [2,[27][28][29]37]. The same can be said of the computation of integral indicators like, e.g., the average sky radiance or the horizontal irradiance [30,38]. The use of a single-scattering approach is justified for our problem under most practically interesting atmospheric conditions, since the short source-detector distances here considered translate into a low probability of double scattering events.…”