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 plt
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:
    plt.ion()
import numpy as np

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 = plt.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)
    plt.xlabel("x [km]")
    plt.ylabel("y [km]")
    plt.title(
        "{0} Product\n{1}".format(attrs["producttype"], attrs["datetime"].isoformat())
    )
    plt.xlim((x[0, 0], x[-1, -1]))
    plt.ylim((y[0, 0], y[-1, -1]))
    plt.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")

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

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

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

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
plt.figure(figsize=(10, 8))
plt.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)
plt.pcolormesh(X, Y, data, cmap=cmap, norm=norm)
plt.xlim((X[0, 0], X[-1, -1]))
plt.ylim((Y[0, 0], Y[-1, -1]))

# add colorbar and do some magic for proper visualisation
cb = plt.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])
plt.title("RADOLAN PG Product \n" + attrs["datetime"].isoformat())
plt.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 0x7fbfff93a910>
../../../_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 0x7fbfeaa724d0>
../../../_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 0x7fc0003abb10>
../../../_images/notebooks_fileio_radolan_radolan_showcase_150_1.png