Using intake to access specMACS products
Posted on Mi 29 Juli 2020 in code examples • by Tobias Kölling
In [1]:
# upgrade some packages aggressively
#import sys
#!{sys.executable} -m pip install --user --upgrade intake netcdf4 pydap
#!{sys.executable} -m pip install --user --upgrade git+https://github.com/leifdenby/intake-xarray#egg=intake-xarray
In [2]:
%pylab inline
In [3]:
from intake import open_catalog
In [4]:
cat = open_catalog("https://raw.githubusercontent.com/eurec4a/eurec4a-intake/master/catalog.yml")
ds = cat.specmacs.cloudmask_20200205.to_dask()
In [5]:
import urllib.request
import yaml
flight = yaml.load(urllib.request.urlopen("https://raw.githubusercontent.com/eurec4a/halo-flight-phase-separation/master/flight_phase_files/EUREC4A_HALO_Flight-Segments_20200205_v1.0.yaml"), Loader=yaml.SafeLoader)
In [6]:
seg = [s for s in flight["segments"] if s["segment_id"] == "HALO-0205_sl1"][0]
In [7]:
selection = ds.sel(time=slice(seg["start"], seg["end"]))
In [8]:
fig, ax = plt.subplots(figsize=(16,4))
selection.swir_radiance.T.plot(cmap="gray")
Out[8]:
In [ ]: