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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_20_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_26_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_32_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_38_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_44_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_50_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_56_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_62_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_68_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_74_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_80_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_86_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_92_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_97_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_105_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_111_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_117_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_123_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_129_0.png

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")

../../../_images/notebooks_fileio_radolan_radolan_showcase_136_0.png

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>

../../../_images/notebooks_fileio_radolan_radolan_showcase_142_1.png

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>

../../../_images/notebooks_fileio_radolan_radolan_showcase_146_1.png

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>

../../../_images/notebooks_fileio_radolan_radolan_showcase_150_1.png

RADOLAN Product Showcase

Contents

RADOLAN Product Showcase#

Setup Environment#

RADOLAN Composite#

National Composites#

RADOLAN RX Product#

RADOLAN RZ Product#

RADOLAN RY Product#

RADOLAN RH Product#

RADOLAN RB Product#

RADOLAN RL Product#

RADOLAN RW Product#

RADOLAN RU Product#

RADOLAN SQ Product#

RADOLAN SH Product#

RADOLAN SF Product#

RADOLAN W1 Product#

RADOLAN W2 Product#

RADOLAN WX Product#

Extended RADOLAN Composites#

RADOLAN EX Product#

RADOLAN EZ Product#

RADOLAN EY Product#

RADOLAN EH Product#

RADOLAN EB Product#

RADOLAN PG/PC Product#

RADOLAN-RADVOR#

RV-Product#

RE-Product#

RQ-Product#