#!/usr/bin/python
# coding: utf8
from mpop.satellites import GeostationaryFactory
from mpop.projector import get_area_def
import datetime
time_slot = datetime.datetime(2014, 05, 01, 15, 30)
global_data = GeostationaryFactory.create_scene("meteosat", "10", "seviri", time_slot)
europe = get_area_def("EuropeCanary")
global_data.load(global_data.image.natural.prerequisites, area_extent=europe.area_extent)
print global_data["VIS006"].data
global_data.image.channel_image(0.6).save("./channelVIS006.png")
print global_data["VIS008"].data
global_data.image.channel_image(0.8).save("./channelVIS008.png")
print global_data["IR_016"].data
global_data.image.channel_image(1.6).save("./channelIR_016.png")

nums = [0.6,
        0.8,
        1.6]
 
names = ['VIS006',
         'VIS008',
         'IR_016']
 
for num, nam in zip(nums, names):
    filename = './channel{}.png'.format(nam)
    global_data.image.channel_image(num).save(filename)