Quickstart¶

Loading data¶

Changed in version 2.0.0-alpha.1: New syntax

To work with weather satellite data, one has to create an instance of the Scene class. In order for mpop to get access to the data, either the current wording directory has to be set to the directory containing the data files, or the base_dir keyword argument has to be provided on scene creation:

>>> import os
>>> os.chdir("/home/a001673/data/satellite/Meteosat-10/seviri/lvl1.5/2015/04/20/HRIT")
>>> from mpop import Scene
>>> from datetime import datetime
>>> time_slot = datetime(2015, 4, 20, 10, 0)
>>> global_scene = Scene(platform_name="Meteosat-10", sensor="seviri", start_time=datetime(2015, 4, 20, 10, 0))

or:

>>> from mpop.scene import Scene
>>> from datetime import datetime
>>> time_slot = datetime(2015, 4, 20, 10, 0)
>>> global_scene = Scene(platform_name="Meteosat-10", sensor="seviri", start_time=datetime(2015, 4, 20, 10, 0), base_dir="/home/a001673/data/satellite/Meteosat-10/seviri/lvl1.5/2015/04/20/HRIT") 
>>>

For some platforms, it might be necessary to also specify an end_time:

>>> Scene(platform_name="SNPP", sensor="viirs", start_time=datetime(2015, 3, 11, 11, 20), end_time=datetime(2015, 3, 11, 11, 26)) 

Loading weather satellite data with mpop is as simple as calling the Scene.load() method:

>>> global_scene.load([0.6, 0.8, 10.8])
>>> print global_scene

seviri/IR_108:
        area: On-the-fly area
        start_time: 2015-04-20 10:00:00
        units: K
        wavelength_range: (9.8, 10.8, 11.8) μm
        shape: (3712, 3712)
seviri/VIS006:
        area: On-the-fly area
        start_time: 2015-04-20 10:00:00
        units: %
        wavelength_range: (0.56, 0.635, 0.71) μm
        shape: (3712, 3712)
seviri/VIS008:
        area: On-the-fly area
        start_time: 2015-04-20 10:00:00
        units: %
        wavelength_range: (0.74, 0.81, 0.88) μm
        shape: (3712, 3712)

As you can see, this loads the visible and IR channels provided as argument to the load() method as a list of wavelengths in micrometers. Another way to load the channels is to provide the names instead:

>>> global_scene.load(["VIS006", "VIS008", "IR_108"])
>>> print global_scene

To have a look at the available bands you should be able to load with your Scene object, you can call the available_datasets() method:

>>> global_scene.available_datasets()

[u'HRV',
 u'IR_108',
 u'IR_120',
 u'VIS006',
 u'WV_062',
 u'IR_039',
 u'IR_134',
 u'IR_097',
 u'IR_087',
 u'VIS008',
 u'IR_016',
 u'WV_073']

To access the loaded data:

>>> print global_scene[0.6]

or:

>>> print global_scene["VIS006"]

To visualize it:

>>> global_scene.show(0.6)

To combine them:

>>> global_scene["ndvi"] = (global_scene[0.8] - global_scene[0.6]) / (global_scene[0.8] + global_scene[0.6])
>>> global_scene.show("ndvi")

Generating composites¶

The easiest way to generate composites is to load them:

>>> global_scene.load(['overview'])
>>> global_scene.show('overview')

To get a list of all available composites for the current scene:

>>> global_scene.available_composites()

[u'overview_sun',
 u'airmass',
 u'natural',
 u'night_fog',
 u'overview',
 u'green_snow',
 u'dust',
 u'fog',
 u'natural_sun',
 u'cloudtop',
 u'convection',
 u'ash']

To save a composite to disk:

>>> global_scene.save_dataset('overview', 'my_nice_overview.png')

One can also specify which writer to use for filenames with non-standard extensions

>>> global_scene.save_dataset('overview', 'my_nice_overview.stupidextension', writer='geotiff')

Resampling¶

Until now, we have used the channels directly as provided by the satellite, that is in satellite projection. Generating composites thus produces views in satellite projection, i.e. as viewed by the satellite.

Most often however, we will want to resample the data onto a specific area so that only the area of interest is depicted in the RGB composites.

Here is how we do that:

>>> local_scene = global_scene.resample("eurol")
>>>

Now we have resampled channel data and composites onto the “eurol” area in the local_scene variable and we can operate as before to display and save RGB composites:

>>> local_scene.show('overview')
>>> local_scene.save_dataset('overview', './local_overview.tif')

The image is automatically saved here in GeoTiff format.

Making custom composites¶

Building custom composites makes use of the RGBCompositor class. For example, building an overview composite can be done manually with:

>>> from mpop.composites import RGBCompositor
>>> compositor = RGBCompositor("myoverview", "bla", "")
>>> composite = compositor([local_scene[0.6],
...                         local_scene[0.8],
...                         local_scene[10.8]])
>>> from mpop.writers import to_image
>>> img = to_image(composite)
>>> img.invert([False, False, True])
>>> img.stretch("linear")
>>> img.gamma(1.7)
>>> img.show()

One important thing to notice is that there is an internal difference between a composite and an image. A composite is defined as a special dataset which may have several bands (like R, G, B bands). However, the data isn’t stretched, or clipped or gamma filtered until an image is generated.