xarray Rainbow5 backend

Contents

xarray Rainbow5 backend#

In this example, we read Rainbow5 data files using the xradar rainbow xarray backend.

[1]:

import glob
import gzip
import io
import wradlib as wrl
import warnings

warnings.filterwarnings("ignore")
import matplotlib.pyplot as plt
import numpy as np
import xradar as xd
import xarray as xr

try:
    get_ipython().run_line_magic("matplotlib inline")
except:
    plt.ion()

Load Rainbow5 Volume Data#

[2]:

fpath = "rainbow/2013051000000600dBZ.vol"
f = wrl.util.get_wradlib_data_file(fpath)
vol = xd.io.open_rainbow_datatree(f, reindex_angle=False)

Downloading file 'rainbow/2013051000000600dBZ.vol' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/rainbow/2013051000000600dBZ.vol' to '/home/runner/work/wradlib-notebooks/wradlib-notebooks/wradlib-data'.

Georeferencing#

[6]:

swp = vol["sweep_0"].ds.copy()
swp = swp.assign_coords(sweep_mode=swp.sweep_mode)
swp = swp.wrl.georef.georeference()

Create simple plot#

Using xarray features a simple plot can be created like this. Note the sortby('time') method, which sorts the radials by time.

For more details on plotting radar data see under Visualization.

[8]:

swp.DBZH.sortby("time").plot(x="range", y="time", add_labels=False)

[8]:

<matplotlib.collections.QuadMesh at 0x7f3e1b30f7d0>

../../../_images/notebooks_fileio_backends_rainbow_backend_16_1.png

[9]:

fig = plt.figure(figsize=(5, 5))
pm = swp.DBZH.wrl.vis.plot(crs={"latmin": 3e3}, fig=fig)

../../../_images/notebooks_fileio_backends_rainbow_backend_17_0.png