Generate a snow RGB using the 3.7 reflectance derived with PySpectral

In [8]:
from satpy.scene import Scene
from glob import glob
FILES = glob("/home/a000680/data/hrit/20150420/*201504201000*__")
In [9]:
glbscene = Scene(filenames=FILES, reader='hrit_msg')
In [3]:
glbscene.load(['snow'])
/home/a000680/usr/src/satpy/satpy/readers/hrit_msg.py:1000: RuntimeWarning: divide by zero encountered in reciprocal
  data.data[:] **= -1
No handlers could be found for logger "pyspectral.near_infrared_reflectance"
In [4]:
lscn = glbscene.resample('germ')
lscn.show('snow')

In [ ]: