#!/usr/bin/env python
# coding: utf-8

# <h1>Read VIIRS SDR data, resample, and make a true color RGB image</h1>
# 
# Created with:
# 
# Package | Version
# --- | ---
# SatPy | 0.9.0
# PyResample | 1.10.1
# Trollimage | 1.5.3
# PyKdtree | 1.3.1
# 

# In[1]:


from satpy.scene import Scene
from glob import glob


# In[2]:


scn = Scene(
        filenames=glob("/home/a000680/data/polar_in/direct_readout/npp/lvl1/npp_20150420_1202_18019/*"),
        reader='viirs_sdr')


# In[3]:


from satpy.resample import get_area_def


# In[4]:


areaid = 'eurol'


# In[5]:


areadef = get_area_def(areaid)


# In[6]:


composite = 'true_color_lowres'


# In[7]:


scn.load([composite])


# In[8]:


newscn = scn.resample(areadef)


# In[9]:


newscn.show(composite)


# <img src='true_color_rayleigh_only_201504201203_eurol_thumb.png'/>

# In[10]:


prfx = newscn.start_time.strftime('%Y%m%d%H%M')


# In[11]:


newscn.save_dataset(
        composite, './true_color_rayleigh_only_{0}_{1}.png'.format(prfx, areaid))