Smoke from forest fires and saharan dust seen with VIIRS on Suomi NPP

On October 17, 2017, an air mass with a heavy load of particles and gases stemming from extensive and devastating forest/grass fires on the Iberian Peninsula reached Northern Europe. Probably the same air mass also caried a significant amount of saharan dust from outbreaks from Western Africa a few days earlier.

The dirty air mass is clearly seen with the VIIRS sensor after correcting for the atmospheric background signal. Pyspectral enables correction of Rayleigh scattering by the background (climatological) atmospheric atoms and absorption by aerosols.

In [1]:
from glob import glob
from satpy.scene import Scene
scn = Scene(
        filenames=glob("/data/lang/proj/safworks/adam/oktober17smoke/npp/lvl1/npp_20171017_1216_30945/*"),
        reader='viirs_sdr')

Load the SDR data

Now load the SDR data needed to generate a pan sharpened (using the I1 band) True Color RGB.
Here pyspectral is being used to do the atmospheric correction.

In [2]:
scn.load(['true_color'])
No handlers could be found for logger "satpy.scene"

Resample the data onto a map-projected area

Here we resample the data onto a cartographic map-projection. We use a predefined area called eurol. For this pyresample and pykdtree are used.

In [3]:
newscn = scn.resample('europe2km')

Show the image

So, now it is time to show the image using trollimage. We also add coastlines and political borders with pycoast.

In [4]:
newscn.show('true_color', overlay={'coast_dir': '/home/a000680/data/shapes/', 'color': (255, 140, 100)})
/home/a000680/usr/src/satpy/satpy/enhancements/__init__.py:51: RuntimeWarning: invalid value encountered in log10
  np.log10(chn.data, out=chn.data)