xarray GAMIC backend

Contents

xarray GAMIC backend#

In this example, we read GAMIC (HDF5) data files using the xradar gamic backend.

import glob
import wradlib as wrl
import wradlib_data
import warnings
from IPython.display import display

import matplotlib.pyplot as plt
import numpy as np
import xradar as xd
import xarray as xr

warnings.filterwarnings("ignore")

Load GAMIC HDF5 Volume Data#

fpath = "hdf5/DWD-Vol-2_99999_20180601054047_00.h5"
f = wradlib_data.DATASETS.fetch(fpath)
vol = xd.io.open_gamic_datatree(f)

Downloading file 'hdf5/DWD-Vol-2_99999_20180601054047_00.h5' from 'https://github.com/wradlib/wradlib-data/raw/main/data/hdf5/DWD-Vol-2_99999_20180601054047_00.h5' to '/home/docs/.cache/wradlib-data'.

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.

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

<matplotlib.collections.QuadMesh at 0x7638d1ec67b0>

../../../_images/f784367a46264dcf3f15f09c5ec90528a0f76304e538ddfced9fc11af01c74e8.png

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

../../../_images/59985cb1e0b6c634615fe2b1d8b2807744861a99f7b38a3cc306d11625a223cb.png