RADOLAN Product Showcase#
In this notebook an overview over the existing RADOLAN products is given.
Setup Environment#
Import modules, filter warnings to avoid cluttering output with DeprecationWarnings and use matplotlib inline or interactive mode if running in ipython or python respectively.
[1]:
import wradlib as wrl
import matplotlib.pyplot as pl
import matplotlib as mpl
import warnings
import io
import os
import glob
import tarfile
warnings.filterwarnings("ignore")
try:
get_ipython().run_line_magic("matplotlib inline")
except:
pl.ion()
import numpy as np
/home/runner/micromamba-root/envs/wradlib-tests/lib/python3.11/site-packages/tqdm/auto.py:22: TqdmWarning: IProgress not found. Please update jupyter and ipywidgets. See https://ipywidgets.readthedocs.io/en/stable/user_install.html
from .autonotebook import tqdm as notebook_tqdm
Get RADOLAN Grid and RADOLAN Extended Grid.
[2]:
# Get coordinates
radolan_grid_xy = wrl.georef.get_radolan_grid(900, 900)
radolan_egrid_xy = wrl.georef.get_radolan_grid(1500, 1400)
radolan_wgrid_xy = wrl.georef.get_radolan_grid(1100, 900)
x = radolan_grid_xy[:, :, 0]
y = radolan_grid_xy[:, :, 1]
xe = radolan_egrid_xy[:, :, 0]
ye = radolan_egrid_xy[:, :, 1]
xw = radolan_wgrid_xy[:, :, 0]
yw = radolan_wgrid_xy[:, :, 1]
Define data reading function and plotting function.
[3]:
def read_radolan(radfile):
radfile = wrl.util.get_wradlib_data_file("radolan/showcase/" + radfile)
return wrl.io.read_radolan_composite(radfile)
[4]:
def plot_radolan(data, attrs, grid, clabel=None):
fig = pl.figure(figsize=(10, 8))
ax = fig.add_subplot(111, aspect="equal")
x = grid[:, :, 0]
y = grid[:, :, 1]
pm = ax.pcolormesh(x, y, data, cmap="viridis")
cb = fig.colorbar(pm, shrink=0.75)
cb.set_label(clabel)
pl.xlabel("x [km]")
pl.ylabel("y [km]")
pl.title(
"{0} Product\n{1}".format(attrs["producttype"], attrs["datetime"].isoformat())
)
pl.xlim((x[0, 0], x[-1, -1]))
pl.ylim((y[0, 0], y[-1, -1]))
pl.grid(color="r")
RADOLAN Composite#
A few products including RW and SF are available free of charge at this DWD FTP Server. A full list of RADOLAN products can be found in the DWD RADOLAN Produktübersicht. Specific details on the RADOLAN Format can be retrieved from the DWD RADOLAN/RADVOR-OP Kompositformat.
Currently, most of the RADOLAN composites have a spatial resolution of 1km x 1km, with the National Composites (R-, S- and W-series) being 900 x 900 km grids, and the European Composites 1500 x 1400 km grids. The polar-stereographic projection is described in the chapter RADOLAN Grid.
One difference is the extended National Composite (only WX) with a 1100 x 900 km grid.
Also the PG/PC-Product with 460 x 460 km grid and runlength-coding is shortly described.
National Composites#
ID |
INT |
avail |
Description |
---|---|---|---|
RX/WX |
5 min |
5 min |
original radardata in qualitative RVP6-units (1 byte coded) |
RZ |
5 min |
5 min |
radardata after correction of PBB converted to rainrate with improved Z-R-relation |
RY |
5 min |
5 min |
radardata after correction with Quality-composit (QY) |
RH |
1 h |
5 min |
1 h summation of RZ-composit |
RB |
1 h |
hh:50 |
1 h summation with preadjustment |
RW |
1 h |
hh:50 |
1 h summation with standard adjustment “best of two” |
RL |
1 h |
hh:50 |
1 h summation with adjustment by Merging |
RU |
1 h |
hh:50 |
1 h summation with standard and merging adjustment “best of three” |
SQ |
6 h |
hh:50 |
6 h summation of RW |
SH |
12 h |
hh:50 |
12 h summation of RW |
SF |
24 h |
hh:50 |
24 h summation of RW |
W1 |
7 d |
05:50 |
7 d summation of RW |
W2 |
14 d |
05:50 |
14 d summation of RW |
W3 |
21 d |
05:50 |
21 d summation of RW |
W4 |
30 d |
05:50 |
30 d summation of RW |
RADOLAN RX Product#
Load data from data source.
[5]:
data, attrs = read_radolan("raa01-rx_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-rx_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-rx_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data and apply scale and offset
[6]:
data = np.ma.masked_equal(data, -9999) / 2 - 32.5
[7]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="dBZ")

RADOLAN RZ Product#
Load data from data source
[8]:
data, attrs = read_radolan("raa01-rz_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-rz_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-rz_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[9]:
data = np.ma.masked_equal(data, -9999)
[10]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * 5min-1")

RADOLAN RY Product#
Load data from data source
[11]:
data, attrs = read_radolan("raa01-ry_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-ry_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-ry_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[12]:
data = np.ma.masked_equal(data, -9999)
[13]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * 5min-1")

RADOLAN RH Product#
Load data from data source
[14]:
data, attrs = read_radolan("raa01-rh_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-rh_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-rh_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[15]:
data = np.ma.masked_equal(data, -9999)
[16]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * h-1")

RADOLAN RB Product#
Load data from data source
[17]:
data, attrs = read_radolan("raa01-rb_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-rb_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-rb_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[18]:
data = np.ma.masked_equal(data, -9999)
[19]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * h-1")

RADOLAN RL Product#
Load data from data source
[20]:
data, attrs = read_radolan("raa01-rl_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-rl_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-rl_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[21]:
data = np.ma.masked_equal(data, -9999)
[22]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * h-1")

RADOLAN RW Product#
Load data from data source
[23]:
data, attrs = read_radolan("raa01-rw_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-rw_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-rw_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[24]:
data = np.ma.masked_equal(data, -9999)
[25]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * h-1")

RADOLAN RU Product#
Load data from data source
[26]:
data, attrs = read_radolan("raa01-ru_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-ru_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-ru_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[27]:
data = np.ma.masked_equal(data, -9999)
[28]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * h-1")

RADOLAN SQ Product#
Load data from data source
[29]:
data, attrs = read_radolan("raa01-sq_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-sq_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-sq_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[30]:
data = np.ma.masked_equal(data, -9999)
[31]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * 6h-1")

RADOLAN SH Product#
Load data from data source
[32]:
data, attrs = read_radolan("raa01-sh_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-sh_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-sh_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[33]:
data = np.ma.masked_equal(data, -9999)
[34]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * 12h-1")

RADOLAN SF Product#
Load data from data source
[35]:
data, attrs = read_radolan("raa01-sf_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-sf_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-sf_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[36]:
data = np.ma.masked_equal(data, -9999)
[37]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * 24h-1")

RADOLAN W1 Product#
Load data from data source
[38]:
data, attrs = read_radolan("raa01-w1_10000-1408110550-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-w1_10000-1408110550-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-w1_10000-1408110550-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[39]:
data = np.ma.masked_equal(data, -9999)
[40]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * 7d-1")

RADOLAN W2 Product#
Load data from data source
[41]:
data, attrs = read_radolan("raa01-w2_10000-1408110550-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-w2_10000-1408110550-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-w2_10000-1408110550-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[42]:
data = np.ma.masked_equal(data, -9999)
[43]:
plot_radolan(data, attrs, radolan_grid_xy, clabel="mm * 14d-1")

RADOLAN WX Product#
[44]:
data, attrs = read_radolan("raa01-wx_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-wx_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-wx_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[45]:
data = np.ma.masked_equal(data, -9999) / 2 - 32.5
[46]:
plot_radolan(data, attrs, radolan_wgrid_xy, clabel="dBZ")

Extended RADOLAN Composites#
The common central european products with a range of 1500 km by 1400 km are presented in the following table:
ID |
INT |
avail |
Description |
---|---|---|---|
EX |
5 min |
5 min |
analogue RX |
EZ |
5 min |
5 min |
analogue RZ |
EY |
5 min |
5 min |
analogue EY after correction with Quality-composit |
EH |
1 h |
hh:50 |
analogue RH (no preadjustment) 1 h summation of EY-composite |
EB |
1 h |
hh:50 |
analogue RB (with preadjustment) 1 h summation |
EW |
1 h |
hh:50 |
analogue RW (full adjustment) 1 h summation |
RADOLAN EX Product#
Load data from data source
[47]:
data, attrs = read_radolan("raa01-ex_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-ex_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-ex_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data and apply scale and offset
[48]:
data = np.ma.masked_equal(data, -9999) / 2 - 32.5
[49]:
plot_radolan(data, attrs, radolan_egrid_xy, clabel="dBZ")

RADOLAN EZ Product#
Load data from data source
[50]:
data, attrs = read_radolan("raa01-ez_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-ez_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-ez_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[51]:
data = np.ma.masked_equal(data, -9999)
[52]:
plot_radolan(data, attrs, radolan_egrid_xy, clabel="mm * 5min-1")

RADOLAN EY Product#
Load data from data source
[53]:
data, attrs = read_radolan("raa01-ey_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-ey_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-ey_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[54]:
data = np.ma.masked_equal(data, -9999)
[55]:
plot_radolan(data, attrs, radolan_egrid_xy, clabel="mm * 5min-1")

RADOLAN EH Product#
Load data from data source
[56]:
data, attrs = read_radolan("raa01-eh_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-eh_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-eh_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[57]:
data = np.ma.masked_equal(data, -9999)
[58]:
plot_radolan(data, attrs, radolan_egrid_xy, clabel="mm * h-1")

RADOLAN EB Product#
Load data from data source
[59]:
data, attrs = read_radolan("raa01-eb_10000-1408102050-dwd---bin.gz")
Downloading file 'radolan/showcase/raa01-eb_10000-1408102050-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/showcase/raa01-eb_10000-1408102050-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[60]:
data = np.ma.masked_equal(data, -9999)
[61]:
plot_radolan(data, attrs, radolan_egrid_xy, clabel="mm * h-1")

RADOLAN PG/PC Product#
The PG/PC product is a bit different from the normal RADOLAN formats. The header is actually the same, but the data is runlength encoded. Also the RADOLAN grid cells have 2km edge length (460x460 cells).
Load data from data source
[62]:
radfile = "raa00-pc_10015-1408030905-dwd---bin.gz"
radfile = wrl.util.get_wradlib_data_file("radolan/misc/" + radfile)
data, attrs = wrl.io.read_radolan_composite(radfile, missing=255)
radolan_grid_pc = wrl.georef.get_radolan_grid(460, 460)
Downloading file 'radolan/misc/raa00-pc_10015-1408030905-dwd---bin.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/misc/raa00-pc_10015-1408030905-dwd---bin.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Mask data
[63]:
data = np.ma.masked_equal(data, 255)
print(data.shape)
(460, 460)
[64]:
# plot the images side by side
pl.figure(figsize=(10, 8))
pl.subplot(111, aspect="equal")
X = radolan_grid_pc[:, :, 0]
Y = radolan_grid_pc[:, :, 1]
# color-scheme taken from DWD "legend_radar_products_pc.pdf"
colors = [
"lightgrey",
"yellow",
"lightblue",
"magenta",
"green",
"red",
"darkblue",
"darkred",
]
cmap = mpl.colors.ListedColormap(colors, name="DWD-pc-scheme")
bounds = np.arange(len(colors) + 1)
norm = mpl.colors.BoundaryNorm(bounds, cmap.N)
pl.pcolormesh(X, Y, data, cmap=cmap, norm=norm)
pl.xlim((X[0, 0], X[-1, -1]))
pl.ylim((Y[0, 0], Y[-1, -1]))
# add colorbar and do some magic for proper visualisation
cb = pl.colorbar(shrink=0.75, norm=norm, boundaries=bounds)
loc = bounds + 0.5
cb.set_ticks(loc[:-1])
labels = bounds[:-1]
cb.set_ticklabels(labels)
cl = cb.ax.get_yticklabels()
cl[-1].set_text("9")
cb.ax.set_yticklabels([elem.get_text() for elem in cl])
pl.title("RADOLAN PG Product \n" + attrs["datetime"].isoformat())
pl.grid(color="r")

RADOLAN-RADVOR#
Further new developments are the “RADVOR” products (Radar-based precipitation forecast in the shortest term). An opendata showcase is given below the following table.
ID |
INT |
avail |
Description |
---|---|---|---|
RV |
5 min |
5 min |
Analyzed and predicted precipitation, 1 analyse, 24 predictions (POLARA) |
RS |
60 min |
5 min |
Analyzed and predicted precipitation, 1 analyse, 24 predictions (POLARA) |
RQ |
60 min |
15 min |
Quantified precipitation analysis and prediction, 1 analyse, 2 predictions |
RE |
60 min |
5 min |
Analysis and prediction of the aggregate state and hail, 1 analyse, 24 predictions (POLARA) |
FS |
60 min |
15 min |
snow depth analysis and forecast, 1 analyse, 2 predictions |
FQ |
360 min |
15 min |
snow depth analysis and forecast, 1 analyse, 2 predictions |
RV-Product#
[65]:
fname = wrl.util.get_wradlib_data_file("radolan/radvor/DE1200_RV2210180700.tar.bz2")
fp = tarfile.open(fname)
fp.extractall()
names = fp.getnames()
buffer = [io.BytesIO(fp.extractfile(name).read()) for name in names]
for buf, name in zip(buffer, names):
buf.name = name
fp.close()
Downloading file 'radolan/radvor/DE1200_RV2210180700.tar.bz2' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/DE1200_RV2210180700.tar.bz2' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
[66]:
ds = wrl.io.open_radolan_mfdataset(buffer)
display(ds)
<xarray.Dataset> Dimensions: (prediction_time: 25, y: 1200, x: 1100, time: 1) Coordinates: * time (time) datetime64[ns] 2022-10-18T07:00:00 * y (y) float64 -4.822e+06 -4.821e+06 ... -3.624e+06 -3.623e+06 * x (x) float64 -5.432e+05 -5.422e+05 ... 5.548e+05 5.558e+05 * prediction_time (prediction_time) datetime64[ns] 2022-10-18T07:00:00 ...... Data variables: RV (prediction_time, y, x) float32 dask.array<chunksize=(1, 1200, 1100), meta=np.ndarray> Attributes: radarid: 10000 formatversion: 5 radolanversion: P300001H radarlocations: ['deasb', 'deboo', 'dedrs', 'deeis', 'deess', 'defbg', '... predictiontime: 0 moduleflag: 8
[67]:
ds.RV.plot(col="prediction_time", col_wrap=5, vmax=20)
[67]:
<xarray.plot.facetgrid.FacetGrid at 0x7f0d2aada490>

RE-Product#
This product isn’t implemented with all features, yet. Use with care!
[68]:
files = [
"radolan/radvor/RQ2210180700_060.gz",
"radolan/radvor/RQ2210180700_000.gz",
"radolan/radvor/RE2210180700_120.gz",
"radolan/radvor/RE2210180700_115.gz",
"radolan/radvor/RE2210180700_110.gz",
"radolan/radvor/RE2210180700_105.gz",
"radolan/radvor/RE2210180700_100.gz",
"radolan/radvor/RE2210180700_095.gz",
"radolan/radvor/RE2210180700_090.gz",
"radolan/radvor/RE2210180700_085.gz",
"radolan/radvor/RE2210180700_080.gz",
"radolan/radvor/RE2210180700_075.gz",
"radolan/radvor/RE2210180700_070.gz",
"radolan/radvor/RE2210180700_065.gz",
"radolan/radvor/RE2210180700_060.gz",
"radolan/radvor/RE2210180700_055.gz",
"radolan/radvor/RE2210180700_050.gz",
"radolan/radvor/RE2210180700_045.gz",
"radolan/radvor/RE2210180700_040.gz",
"radolan/radvor/RE2210180700_035.gz",
"radolan/radvor/RE2210180700_030.gz",
"radolan/radvor/RE2210180700_025.gz",
"radolan/radvor/RE2210180700_020.gz",
"radolan/radvor/RE2210180700_015.gz",
"radolan/radvor/RE2210180700_010.gz",
"radolan/radvor/RE2210180700_005.gz",
"radolan/radvor/RE2210180700_000.gz",
]
files = [wrl.util.get_wradlib_data_file(f) for f in files]
Downloading file 'radolan/radvor/RQ2210180700_060.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RQ2210180700_060.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RQ2210180700_000.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RQ2210180700_000.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_120.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_120.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_115.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_115.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_110.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_110.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_105.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_105.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_100.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_100.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_095.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_095.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_090.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_090.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_085.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_085.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_080.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_080.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_075.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_075.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_070.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_070.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_065.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_065.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_060.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_060.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_055.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_055.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_050.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_050.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_045.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_045.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_040.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_040.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_035.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_035.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_030.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_030.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_025.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_025.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_020.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_020.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_015.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_015.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_010.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_010.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_005.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_005.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
Downloading file 'radolan/radvor/RE2210180700_000.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RE2210180700_000.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
[69]:
ds = wrl.io.open_radolan_mfdataset(files)
display(ds)
<xarray.Dataset> Dimensions: (time: 1, y: 900, x: 900, prediction_time: 25) Coordinates: * time (time) datetime64[ns] 2022-10-18T07:00:00 * y (y) float64 -4.672e+06 -4.671e+06 ... -3.774e+06 -3.773e+06 * x (x) float64 -5.232e+05 -5.222e+05 ... 3.748e+05 3.758e+05 * prediction_time (prediction_time) datetime64[ns] 2022-10-18T07:00:00 ...... Data variables: RE (prediction_time, y, x) float32 dask.array<chunksize=(1, 900, 900), meta=np.ndarray> RQ (prediction_time, y, x) float32 dask.array<chunksize=(12, 900, 900), meta=np.ndarray> Attributes: radarid: 10000 formatversion: 5 radolanversion: P300001H radarlocations: ['deasb', 'deboo', 'dedrs', 'deeis', 'deess', 'defbg', '... predictiontime: 0 moduleflag: 8 quantification: 16
[70]:
ds.RE.plot(col="prediction_time", col_wrap=5, vmax=2)
[70]:
<xarray.plot.facetgrid.FacetGrid at 0x7f0d15af9950>

RQ-Product#
[71]:
files = [
"radolan/radvor/RQ2210180700_120.gz",
"radolan/radvor/RQ2210180700_060.gz",
"radolan/radvor/RQ2210180700_000.gz",
]
files = [wrl.util.get_wradlib_data_file(f) for f in files]
Downloading file 'radolan/radvor/RQ2210180700_120.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/radolan/radvor/RQ2210180700_120.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.
[72]:
ds = wrl.io.open_radolan_mfdataset(files)
display(ds)
<xarray.Dataset> Dimensions: (prediction_time: 3, y: 900, x: 900, time: 1) Coordinates: * time (time) datetime64[ns] 2022-10-18T07:00:00 * y (y) float64 -4.672e+06 -4.671e+06 ... -3.774e+06 -3.773e+06 * x (x) float64 -5.232e+05 -5.222e+05 ... 3.748e+05 3.758e+05 * prediction_time (prediction_time) datetime64[ns] 2022-10-18T07:00:00 ...... Data variables: RQ (prediction_time, y, x) float32 dask.array<chunksize=(1, 900, 900), meta=np.ndarray> Attributes: radarid: 10000 formatversion: 5 radolanversion: 2.29.1 radarlocations: ['asb', 'boo', 'drs', 'eis', 'ess', 'fbg', 'fld', 'hnr',... predictiontime: 0 moduleflag: 8 quantification: 0
[73]:
ds.RQ.plot(col="prediction_time")
[73]:
<xarray.plot.facetgrid.FacetGrid at 0x7f0d158f0250>
