Projection of the plot: satellite viewΒΆ
[1]:
%matplotlib inline
# for figures in notebook
# import & initialize epygram
import epygram
epygram.init_env()
import os
INPUTS_DIR = os.path.join('..', 'inputs')
[2]:
r = epygram.formats.resource(os.path.join(INPUTS_DIR, 'grid.arome-forecast.guyane0025+0012:00.grib'), 'r')
[3]:
f = r.readfield('shortName:2r')
[4]:
from cartopy.crs import NearsidePerspective
proj = NearsidePerspective(central_longitude=-60.,
central_latitude=10.,
satellite_height=1e6)
fig, ax = f.cartoplot(projection=proj, extent='global', colormap='viridis_r')
[ ]: