SEVIRI RGB with the derived 3.9 micron reflectance

>>> from mpop.satellites import GeostationaryFactory
>>> from datetime import datetime
>>> tslot = datetime(2015, 4, 20, 10, 0)
>>> glbd = GeostationaryFactory.create_scene("Meteosat-10", "", "seviri", tslot)

Load only the bands needed in order to create the snow RGB:

>>> glbd.load(glbd.image.snow.prerequisites)
>>> lcd = glbd.project('euro4')

And now we derive the snow rgb according to the MSG Interpretation Guide. It requires pyspectral:

>>> img = lcd.image.snow()
>>> img.add_overlay()

Necessary configuration settings

For this tutorial template config files (see Installation and configuration) can be used. These are located in the etc dir of the mpop source. Copy mpop.cfg.template, areas.def.template and Meteosat-9.cfg.template to another dir and remove the .template extension. In the config file Meteosat-9.cfg locate the section severi-level1 and modify the defined dir to point to the dir of your uncompressed HRIT data.

Set PPP_CONFIG_DIR to the directory containing your modified mpop config files.