We present a method to map the artificial sky brightness across large territories in astronomical photometric bands with a resolution of approximately 1 km. This is of use in quantifying the situation regarding night sky pollution, recognizing potential astronomical sites and allowing future monitoring of trends. The artificial sky brightness present in the chosen direction at a given position on the surface of the Earth is obtained by the integration of the contributions produced by every surface area in the surroundings. Each contribution is computed via detailed models for the propagation in the atmosphere of the upward light flux emitted by the area. The light flux is measured with top-of-atmosphere radiometric observations made by the Defense Meteorological Satellite Program (DMSP) Operational Linescan System.

We have applied the described method to Europe, obtaining maps of artificial sky brightness in the V and B bands.