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 datatree as xt
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/wradlib/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 0x7f37444037d0>

../../../_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