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

# # Hurricane Trami,  September 28th 2018

# In[1]:


from satpy import Scene, find_files_and_readers
from satpy.resample import get_area_def
from datetime import datetime


# In[2]:


files = find_files_and_readers(start_time=datetime(2018, 9, 28, 6, 0),
                               end_time=datetime(2018, 9, 28, 6, 10),
                               base_dir="/data/h8",
                               reader='ahi_hrit')


# In[3]:


scn = Scene(filenames=files)


# In[4]:


rgbname = 'true_color'
scn.load([rgbname])


# In[5]:


new_scn = scn.resample(scn.min_area(), resampler='native')
new_scn.show(rgbname)


# In[6]:


from pyresample.utils import get_area_def
area_id = 'east_asia'
x_size = 2407
y_size = 1655
area_extent = (-3603877.326193321, -2254860.2505908236, 1203877.3261933187, 1050373.6358697552)
projection = '+proj=laea +lat_0=37.5 +lon_0=137.5 +ellps=WGS84'
description = "East Asia"
proj_id = 'laea_137.5_37.5'
areadef = get_area_def(area_id, description, proj_id, projection,x_size, y_size, area_extent)


# In[7]:


local_scene = scn.resample(areadef)


# In[8]:


local_scene.show(rgbname)


# In[ ]:





# In[ ]:





# In[ ]: