ToGrid

In this notebook we show the production of a reflectivity composite from 3 neighboring radars to a common cartesian grid, using the sampling volume size as a quality criterion.

import tempfile
import warnings

import cmweather
import numpy as np
import pyproj
import xarray as xr
import matplotlib.pyplot as plt

import wradlib as wrl
import wradlib_data

warnings.filterwarnings("ignore")

/home/docs/checkouts/readthedocs.org/user_builds/wradlib/conda/latest/lib/python3.13/site-packages/tqdm/auto.py:21: 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 radar data

First, we import measurements from three belgian radars. This is done using xarray.open_dataset using xradar.io.backends.odim.OdimBackendEntrypoint.

filenames = ["bejab.pvol.hdf", "bewid.pvol.hdf", "behel.pvol.hdf"]
paths = {f.split(".")[0][2:]: wradlib_data.DATASETS.fetch(f"hdf5/{f}") for f in filenames}
ctree = xr.DataTree()
for radar, filename in paths.items():
    ctree[radar] = xr.open_dataset(filename, engine="odim", group="sweep_0").chunk().wrl.georef.georeference().set_coords("sweep_mode")
display(ctree)

<xarray.DataTree>
Group: /
├── Group: /jab
│       Dimensions:            (azimuth: 360, range: 598)
│       Coordinates: (12/15)
│         * azimuth            (azimuth) float32 1kB 0.5 1.5 2.5 ... 357.5 358.5 359.5
│           elevation          (azimuth) float64 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│           time               (azimuth) datetime64[ns] 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│         * range              (range) float32 2kB 250.0 750.0 ... 2.982e+05 2.988e+05
│           x                  (azimuth, range) float64 2MB 2.182 6.545 ... -2.605e+03
│           y                  (azimuth, range) float64 2MB 250.0 750.0 ... 2.986e+05
│           ...                 ...
│           bins               (azimuth, range) float32 861kB 250.0 750.0 ... 2.988e+05
│           sweep_mode         <U20 80B 'azimuth_surveillance'
│           longitude          float64 8B 3.064
│           latitude           float64 8B 51.19
│           altitude           float64 8B 50.0
│           crs_wkt            int64 8B 0
│       Data variables:
│           DBZH               (azimuth, range) float64 2MB dask.array<chunksize=(360, 598), meta=np.ndarray>
│           sweep_number       int64 8B ...
│           prt_mode           <U7 28B ...
│           follow_mode        <U7 28B ...
│           sweep_fixed_angle  float64 8B ...
│           nyquist_velocity   object 8B ...
│       Attributes:
│           Conventions:  ODIM_H5/V2_2
├── Group: /wid
│       Dimensions:            (azimuth: 360, range: 1000)
│       Coordinates: (12/15)
│         * azimuth            (azimuth) float32 1kB 0.5 1.5 2.5 ... 357.5 358.5 359.5
│           elevation          (azimuth) float64 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│           time               (azimuth) datetime64[ns] 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│         * range              (range) float32 4kB 125.0 375.0 ... 2.496e+05 2.499e+05
│           x                  (azimuth, range) float64 3MB 1.091 3.272 ... -2.179e+03
│           y                  (azimuth, range) float64 3MB 125.0 375.0 ... 2.497e+05
│           ...                 ...
│           bins               (azimuth, range) float32 1MB 125.0 375.0 ... 2.499e+05
│           sweep_mode         <U20 80B 'azimuth_surveillance'
│           longitude          float64 8B 5.506
│           latitude           float64 8B 49.91
│           altitude           float64 8B 590.0
│           crs_wkt            int64 8B 0
│       Data variables:
│           DBZH               (azimuth, range) float64 3MB dask.array<chunksize=(360, 1000), meta=np.ndarray>
│           sweep_number       int64 8B ...
│           prt_mode           <U7 28B ...
│           follow_mode        <U7 28B ...
│           sweep_fixed_angle  float64 8B ...
│           nyquist_velocity   object 8B ...
│       Attributes:
│           Conventions:  ODIM_H5/V2_2
└── Group: /hel
        Dimensions:            (azimuth: 360, range: 800)
        Coordinates: (12/15)
          * azimuth            (azimuth) float32 1kB 0.5 1.5 2.5 ... 357.5 358.5 359.5
            elevation          (azimuth) float64 3kB dask.array<chunksize=(360,), meta=np.ndarray>
            time               (azimuth) datetime64[ns] 3kB dask.array<chunksize=(360,), meta=np.ndarray>
          * range              (range) float32 3kB 125.0 375.0 ... 1.996e+05 1.999e+05
            x                  (azimuth, range) float64 2MB 1.091 3.272 ... -1.744e+03
            y                  (azimuth, range) float64 2MB 125.0 375.0 ... 1.998e+05
            ...                 ...
            bins               (azimuth, range) float32 1MB 125.0 375.0 ... 1.999e+05
            sweep_mode         <U20 80B 'azimuth_surveillance'
            longitude          float64 8B 5.406
            latitude           float64 8B 51.07
            altitude           float64 8B 140.0
            crs_wkt            int64 8B 0
        Data variables:
            DBZH               (azimuth, range) float64 2MB dask.array<chunksize=(360, 800), meta=np.ndarray>
            sweep_number       int64 8B ...
            prt_mode           <U7 28B ...
            follow_mode        <U7 28B ...
            sweep_fixed_angle  float64 8B ...
            nyquist_velocity   object 8B ...
        Attributes:
            Conventions:  ODIM_H5/V2_2

xarray.DataTree

• /jab(22)

  • Dimensions:
    • azimuth: 360
    • range: 598
  • Coordinates: (15)
    • azimuth
      (azimuth)
      float32
      0.5 1.5 2.5 ... 357.5 358.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5],shape=(360,), dtype=float32)

    • elevation
      (azimuth)
      float64
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      ray_elevation_angle

      long_name :

      elevation_angle_from_horizontal_plane

      units :

      degrees

      axis :

      radial_elevation_coordinate

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • time
      (azimuth)
      datetime64[ns]
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      time

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type datetime64[ns] numpy.ndarray

    • range
      (range)
      float32
      250.0 750.0 ... 2.982e+05 2.988e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      500.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      250.0

      array([2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05], shape=(598,), dtype=float32)

    • x
      (azimuth, range)
      float64
      2.182 6.545 ... -2.605e+03

      axis :

      X

      long_name :

      Easting

      standard_name :

      projection_x_coordinate

      units :

      metre

      array([[ 2.18159727e+00,  6.54478976e+00,  1.09079795e+01, ...,2.59674093e+03,  2.60109715e+03,  2.60545336e+03],[ 6.54411939e+00,  1.96323521e+01,  3.27205766e+01, ...,7.78942244e+03,  7.80248976e+03,  7.81555703e+03],[ 1.09046481e+01,  3.27139340e+01,  5.45232063e+01, ...,1.29797311e+04,  1.30015056e+04,  1.30232799e+04],...,[-1.09046863e+01, -3.27140486e+01, -5.45233972e+01, ...,-1.29797766e+04, -1.30015511e+04, -1.30233255e+04],[-6.54406817e+00, -1.96321984e+01, -3.27203204e+01, ...,-7.78936147e+03, -7.80242869e+03, -7.81549585e+03],[-2.18157538e+00, -6.54472410e+00, -1.09078701e+01, ...,-2.59671488e+03, -2.60107106e+03, -2.60542722e+03]],shape=(360, 598))

    • y
      (azimuth, range)
      float64
      250.0 750.0 ... 2.981e+05 2.986e+05

      axis :

      Y

      long_name :

      Northing

      standard_name :

      projection_y_coordinate

      units :

      metre

      array([[2.49985540e+02, 7.49956387e+02, 1.24992692e+03, ...,2.97556150e+05, 2.98055323e+05, 2.98554492e+05],[2.49909397e+02, 7.49727957e+02, 1.24954620e+03, ...,2.97465517e+05, 2.97964537e+05, 2.98463555e+05],[2.49757125e+02, 7.49271141e+02, 1.24878484e+03, ...,2.97284269e+05, 2.97782985e+05, 2.98281699e+05],...,[2.49757125e+02, 7.49271141e+02, 1.24878484e+03, ...,2.97284269e+05, 2.97782985e+05, 2.98281699e+05],[2.49909397e+02, 7.49727957e+02, 1.24954620e+03, ...,2.97465517e+05, 2.97964537e+05, 2.98463555e+05],[2.49985540e+02, 7.49956387e+02, 1.24992692e+03, ...,2.97556150e+05, 2.98055323e+05, 2.98554492e+05]], shape=(360, 598))

    • z
      (azimuth, range)
      float64
      51.31 53.96 ... 6.849e+03 6.87e+03

      standard_name :

      height_above_ground

      units :

      metre

      array([[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],...,[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573]], shape=(360, 598))

    • gr
      (azimuth, range)
      float64
      250.0 750.0 ... 2.981e+05 2.986e+05

      standard_name :

      distance_from_radar

      units :

      metre

      array([[2.49995061e+02, 7.49984946e+02, 1.24997452e+03, ...,2.97567481e+05, 2.98066672e+05, 2.98565861e+05],[2.49995066e+02, 7.49984961e+02, 1.24997454e+03, ...,2.97567487e+05, 2.98066678e+05, 2.98565867e+05],[2.49995068e+02, 7.49984965e+02, 1.24997455e+03, ...,2.97567488e+05, 2.98066679e+05, 2.98565868e+05],...,[2.49995069e+02, 7.49984970e+02, 1.24997456e+03, ...,2.97567490e+05, 2.98066681e+05, 2.98565870e+05],[2.49995065e+02, 7.49984957e+02, 1.24997454e+03, ...,2.97567485e+05, 2.98066676e+05, 2.98565865e+05],[2.49995061e+02, 7.49984946e+02, 1.24997452e+03, ...,2.97567481e+05, 2.98066672e+05, 2.98565861e+05]], shape=(360, 598))

    • rays
      (azimuth, range)
      float32
      0.5 0.5 0.5 ... 359.5 359.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],[  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],[  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],...,[357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],[358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],[359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]],shape=(360, 598), dtype=float32)

    • bins
      (azimuth, range)
      float32
      250.0 750.0 ... 2.982e+05 2.988e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      500.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      250.0

      array([[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],...,[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05]], shape=(360, 598), dtype=float32)

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      3.064

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(3.0642)

    • latitude
      ()
      float64
      51.19

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(51.1917)

    • altitude
      ()
      float64
      50.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(50.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["AEQD centered at 51.1917N 3.0642E",BASEGEOGCRS["unknown",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]]],CONVERSION["unknown",METHOD["Azimuthal Equidistant",ID["EPSG",1125]],PARAMETER["Latitude of natural origin",51.1917,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",3.0642,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      unknown

      horizontal_datum_name :

      World Geodetic System 1984

      projected_crs_name :

      AEQD centered at 51.1917N 3.0642E

      grid_mapping_name :

      azimuthal_equidistant

      latitude_of_projection_origin :

      51.1917

      longitude_of_projection_origin :

      3.0642

      false_easting :

      0.0

      false_northing :

      0.0

      array(0)

  • Data variables: (6)
    • DBZH
      (azimuth, range)
      float64
      dask.array<chunksize=(360, 598), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 1.64 MiB 1.64 MiB Shape (360, 598) (360, 598) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

  • Attributes: (1)

    Conventions :

    ODIM_H5/V2_2

  /wid(22)

  • Dimensions:
    • azimuth: 360
    • range: 1000
  • Coordinates: (15)
    • azimuth
      (azimuth)
      float32
      0.5 1.5 2.5 ... 357.5 358.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5],shape=(360,), dtype=float32)

    • elevation
      (azimuth)
      float64
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      ray_elevation_angle

      long_name :

      elevation_angle_from_horizontal_plane

      units :

      degrees

      axis :

      radial_elevation_coordinate

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • time
      (azimuth)
      datetime64[ns]
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      time

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type datetime64[ns] numpy.ndarray

    • range
      (range)
      float32
      125.0 375.0 ... 2.496e+05 2.499e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05, 2.49625e+05,2.49875e+05], shape=(1000,), dtype=float32)

    • x
      (azimuth, range)
      float64
      1.091 3.272 ... -2.179e+03

      axis :

      X

      long_name :

      Easting

      standard_name :

      projection_x_coordinate

      units :

      metre

      array([[ 1.09072932e+00,  3.27218746e+00,  5.45364492e+00, ...,2.17504505e+03,  2.17722396e+03,  2.17940286e+03],[ 3.27185179e+00,  9.81555385e+00,  1.63592539e+01, ...,6.52446477e+03,  6.53100082e+03,  6.53753685e+03],[ 5.45197759e+00,  1.63559302e+01,  2.72598795e+01, ...,1.08718970e+04,  1.08827882e+04,  1.08936794e+04],...,[-5.45199668e+00, -1.63559875e+01, -2.72599750e+01, ...,-1.08719351e+04, -1.08828263e+04, -1.08937175e+04],[-3.27182618e+00, -9.81547702e+00, -1.63591258e+01, ...,-6.52441371e+03, -6.53094970e+03, -6.53748568e+03],[-1.09071838e+00, -3.27215464e+00, -5.45359021e+00, ...,-2.17502323e+03, -2.17720212e+03, -2.17938099e+03]],shape=(360, 1000))

    • y
      (azimuth, range)
      float64
      125.0 375.0 ... 2.495e+05 2.497e+05

      axis :

      Y

      long_name :

      Northing

      standard_name :

      projection_y_coordinate

      units :

      metre

      array([[1.24984828e+02, 3.74954426e+02, 6.24923946e+02, ...,2.49234733e+05, 2.49484410e+05, 2.49734087e+05],[1.24946759e+02, 3.74840218e+02, 6.24733599e+02, ...,2.49158818e+05, 2.49408419e+05, 2.49658020e+05],[1.24870627e+02, 3.74611824e+02, 6.24352944e+02, ...,2.49007004e+05, 2.49256452e+05, 2.49505901e+05],...,[1.24870627e+02, 3.74611824e+02, 6.24352944e+02, ...,2.49007004e+05, 2.49256452e+05, 2.49505901e+05],[1.24946759e+02, 3.74840218e+02, 6.24733599e+02, ...,2.49158818e+05, 2.49408419e+05, 2.49658020e+05],[1.24984828e+02, 3.74954426e+02, 6.24923946e+02, ...,2.49234733e+05, 2.49484410e+05, 2.49734087e+05]],shape=(360, 1000))

    • z
      (azimuth, range)
      float64
      590.7 592.0 ... 5.566e+03 5.575e+03

      standard_name :

      height_above_ground

      units :

      metre

      array([[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],...,[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531]], shape=(360, 1000))

    • gr
      (azimuth, range)
      float64
      125.0 375.0 ... 2.495e+05 2.497e+05

      standard_name :

      distance_from_radar

      units :

      metre

      array([[1.24989588e+02, 3.74968705e+02, 6.24947743e+02, ...,2.49244223e+05, 2.49493910e+05, 2.49743596e+05],[1.24989591e+02, 3.74968712e+02, 6.24947755e+02, ...,2.49244228e+05, 2.49493915e+05, 2.49743601e+05],[1.24989591e+02, 3.74968714e+02, 6.24947758e+02, ...,2.49244230e+05, 2.49493916e+05, 2.49743602e+05],...,[1.24989592e+02, 3.74968716e+02, 6.24947763e+02, ...,2.49244231e+05, 2.49493918e+05, 2.49743604e+05],[1.24989590e+02, 3.74968710e+02, 6.24947752e+02, ...,2.49244227e+05, 2.49493913e+05, 2.49743600e+05],[1.24989588e+02, 3.74968704e+02, 6.24947743e+02, ...,2.49244223e+05, 2.49493910e+05, 2.49743596e+05]],shape=(360, 1000))

    • rays
      (azimuth, range)
      float32
      0.5 0.5 0.5 ... 359.5 359.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],[  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],[  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],...,[357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],[358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],[359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]],shape=(360, 1000), dtype=float32)

    • bins
      (azimuth, range)
      float32
      125.0 375.0 ... 2.496e+05 2.499e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],...,[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05]], shape=(360, 1000), dtype=float32)

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      5.506

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(5.5056)

    • latitude
      ()
      float64
      49.91

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(49.9143)

    • altitude
      ()
      float64
      590.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(590.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["AEQD centered at 49.9143N 5.5056E",BASEGEOGCRS["unknown",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]]],CONVERSION["unknown",METHOD["Azimuthal Equidistant",ID["EPSG",1125]],PARAMETER["Latitude of natural origin",49.9143,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",5.5056,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      unknown

      horizontal_datum_name :

      World Geodetic System 1984

      projected_crs_name :

      AEQD centered at 49.9143N 5.5056E

      grid_mapping_name :

      azimuthal_equidistant

      latitude_of_projection_origin :

      49.9143

      longitude_of_projection_origin :

      5.5056

      false_easting :

      0.0

      false_northing :

      0.0

      array(0)

  • Data variables: (6)
    • DBZH
      (azimuth, range)
      float64
      dask.array<chunksize=(360, 1000), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 2.75 MiB 2.75 MiB Shape (360, 1000) (360, 1000) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

  • Attributes: (1)

    Conventions :

    ODIM_H5/V2_2

  /hel(22)

  • Dimensions:
    • azimuth: 360
    • range: 800
  • Coordinates: (15)
    • azimuth
      (azimuth)
      float32
      0.5 1.5 2.5 ... 357.5 358.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5],shape=(360,), dtype=float32)

    • elevation
      (azimuth)
      float64
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      ray_elevation_angle

      long_name :

      elevation_angle_from_horizontal_plane

      units :

      degrees

      axis :

      radial_elevation_coordinate

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • time
      (azimuth)
      datetime64[ns]
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      time

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type datetime64[ns] numpy.ndarray

    • range
      (range)
      float32
      125.0 375.0 ... 1.996e+05 1.999e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05, 1.99625e+05,1.99875e+05], shape=(800,), dtype=float32)

    • x
      (azimuth, range)
      float64
      1.091 3.272 ... -1.744e+03

      axis :

      X

      long_name :

      Easting

      standard_name :

      projection_x_coordinate

      units :

      metre

      array([[ 1.09078715e+00,  3.27236094e+00,  5.45393406e+00, ...,1.73927186e+03,  1.74145170e+03,  1.74363153e+03],[ 3.27202525e+00,  9.81607424e+00,  1.63601212e+01, ...,5.21727952e+03,  5.22381835e+03,  5.23035717e+03],[ 5.45226664e+00,  1.63567974e+01,  2.72613247e+01, ...,8.69369790e+03,  8.70459373e+03,  8.71548955e+03],...,[-5.45228573e+00, -1.63568547e+01, -2.72614202e+01, ...,-8.69372834e+03, -8.70462421e+03, -8.71552006e+03],[-3.27199964e+00, -9.81599741e+00, -1.63599931e+01, ...,-5.21723869e+03, -5.22377746e+03, -5.23031623e+03],[-1.09077621e+00, -3.27232811e+00, -5.45387934e+00, ...,-1.73925442e+03, -1.74143423e+03, -1.74361404e+03]],shape=(360, 800))

    • y
      (azimuth, range)
      float64
      125.0 375.0 ... 1.996e+05 1.998e+05

      axis :

      Y

      long_name :

      Northing

      standard_name :

      projection_y_coordinate

      units :

      metre

      array([[1.24991454e+02, 3.74974305e+02, 6.24957077e+02, ...,1.99300221e+05, 1.99550004e+05, 1.99799787e+05],[1.24953383e+02, 3.74860091e+02, 6.24766721e+02, ...,1.99239515e+05, 1.99489223e+05, 1.99738930e+05],[1.24877248e+02, 3.74631685e+02, 6.24386045e+02, ...,1.99118117e+05, 1.99367672e+05, 1.99617227e+05],...,[1.24877248e+02, 3.74631685e+02, 6.24386045e+02, ...,1.99118117e+05, 1.99367672e+05, 1.99617227e+05],[1.24953383e+02, 3.74860091e+02, 6.24766721e+02, ...,1.99239515e+05, 1.99489223e+05, 1.99738930e+05],[1.24991454e+02, 3.74974305e+02, 6.24957077e+02, ...,1.99300221e+05, 1.99550004e+05, 1.99799787e+05]], shape=(360, 800))

    • z
      (azimuth, range)
      float64
      140.7 142.0 ... 3.532e+03 3.539e+03

      standard_name :

      height_above_ground

      units :

      metre

      array([[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],...,[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903]], shape=(360, 800))

    • gr
      (azimuth, range)
      float64
      125.0 375.0 ... 1.996e+05 1.998e+05

      standard_name :

      distance_from_radar

      units :

      metre

      array([[1.24996215e+02, 3.74988584e+02, 6.24980876e+02, ...,1.99307810e+05, 1.99557603e+05, 1.99807395e+05],[1.24996217e+02, 3.74988591e+02, 6.24980888e+02, ...,1.99307813e+05, 1.99557606e+05, 1.99807399e+05],[1.24996218e+02, 3.74988593e+02, 6.24980891e+02, ...,1.99307814e+05, 1.99557607e+05, 1.99807400e+05],...,[1.24996219e+02, 3.74988596e+02, 6.24980895e+02, ...,1.99307816e+05, 1.99557609e+05, 1.99807401e+05],[1.24996217e+02, 3.74988589e+02, 6.24980885e+02, ...,1.99307812e+05, 1.99557605e+05, 1.99807398e+05],[1.24996215e+02, 3.74988584e+02, 6.24980875e+02, ...,1.99307809e+05, 1.99557602e+05, 1.99807395e+05]], shape=(360, 800))

    • rays
      (azimuth, range)
      float32
      0.5 0.5 0.5 ... 359.5 359.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],[  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],[  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],...,[357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],[358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],[359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]],shape=(360, 800), dtype=float32)

    • bins
      (azimuth, range)
      float32
      125.0 375.0 ... 1.996e+05 1.999e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],...,[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05]], shape=(360, 800), dtype=float32)

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      5.406

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(5.4064)

    • latitude
      ()
      float64
      51.07

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(51.069072)

    • altitude
      ()
      float64
      140.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(140.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["AEQD centered at 51.069072N 5.4064E",BASEGEOGCRS["unknown",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]]],CONVERSION["unknown",METHOD["Azimuthal Equidistant",ID["EPSG",1125]],PARAMETER["Latitude of natural origin",51.069072,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",5.4064,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      unknown

      horizontal_datum_name :

      World Geodetic System 1984

      projected_crs_name :

      AEQD centered at 51.069072N 5.4064E

      grid_mapping_name :

      azimuthal_equidistant

      latitude_of_projection_origin :

      51.069072

      longitude_of_projection_origin :

      5.4064

      false_easting :

      0.0

      false_northing :

      0.0

      array(0)

  • Data variables: (6)
    • DBZH
      (azimuth, range)
      float64
      dask.array<chunksize=(360, 800), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 2.20 MiB 2.20 MiB Shape (360, 800) (360, 800) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

  • Attributes: (1)

    Conventions :

    ODIM_H5/V2_2

Plot Overview

fig, axs = plt.subplots(1, 3, figsize=(18, 5), sharey=True)
ax = axs.flat[0]
ctree["jab"].ds.DBZH.wrl.vis.plot(ax=ax, vmin=0, vmax=60)
ax.set_aspect("equal")
ax.set_title("Radar Jabbeke")
ax = axs.flat[1]
ctree["wid"].ds.DBZH.wrl.vis.plot(ax=ax, vmin=0, vmax=60)
ax.set_aspect("equal")
ax.set_title("Radar Wideumont")
ax = axs.flat[2]
ctree["hel"].ds.DBZH.wrl.vis.plot(ax=ax, vmin=0, vmax=60)
ax.set_aspect("equal")
ax.set_title("Radar Helchteren")
fig.tight_layout()

Georeference UTM

proj_utm = pyproj.CRS.from_epsg(32632)
kwargs = dict(crs=proj_utm)
for radar in ctree.children:
    ctree[radar].ds = ctree[radar].ds.wrl.georef.georeference(crs=proj_utm)
display(ctree)

<xarray.DataTree>
Group: /
├── Group: /jab
│       Dimensions:            (azimuth: 360, range: 598)
│       Coordinates: (12/15)
│         * azimuth            (azimuth) float32 1kB 0.5 1.5 2.5 ... 357.5 358.5 359.5
│           elevation          (azimuth) float64 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│           time               (azimuth) datetime64[ns] 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│         * range              (range) float32 2kB 250.0 750.0 ... 2.982e+05 2.988e+05
│           x                  (azimuth, range) float64 2MB 8.539e+04 ... 1.074e+05
│           y                  (azimuth, range) float64 2MB 5.688e+06 ... 5.986e+06
│           ...                 ...
│           bins               (azimuth, range) float32 861kB 250.0 750.0 ... 2.988e+05
│           sweep_mode         <U20 80B 'azimuth_surveillance'
│           longitude          float64 8B 3.064
│           latitude           float64 8B 51.19
│           altitude           float64 8B 50.0
│           crs_wkt            int64 8B 0
│       Data variables:
│           DBZH               (azimuth, range) float64 2MB dask.array<chunksize=(360, 598), meta=np.ndarray>
│           sweep_number       int64 8B ...
│           prt_mode           <U7 28B ...
│           follow_mode        <U7 28B ...
│           sweep_fixed_angle  float64 8B ...
│           nyquist_velocity   object 8B ...
│       Attributes:
│           Conventions:  ODIM_H5/V2_2
├── Group: /wid
│       Dimensions:            (azimuth: 360, range: 1000)
│       Coordinates: (12/15)
│         * azimuth            (azimuth) float32 1kB 0.5 1.5 2.5 ... 357.5 358.5 359.5
│           elevation          (azimuth) float64 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│           time               (azimuth) datetime64[ns] 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│         * range              (range) float32 4kB 125.0 375.0 ... 2.496e+05 2.499e+05
│           x                  (azimuth, range) float64 3MB 2.492e+05 ... 2.588e+05
│           y                  (azimuth, range) float64 3MB 5.535e+06 ... 5.785e+06
│           ...                 ...
│           bins               (azimuth, range) float32 1MB 125.0 375.0 ... 2.499e+05
│           sweep_mode         <U20 80B 'azimuth_surveillance'
│           longitude          float64 8B 5.506
│           latitude           float64 8B 49.91
│           altitude           float64 8B 590.0
│           crs_wkt            int64 8B 0
│       Data variables:
│           DBZH               (azimuth, range) float64 3MB dask.array<chunksize=(360, 1000), meta=np.ndarray>
│           sweep_number       int64 8B ...
│           prt_mode           <U7 28B ...
│           follow_mode        <U7 28B ...
│           sweep_fixed_angle  float64 8B ...
│           nyquist_velocity   object 8B ...
│       Attributes:
│           Conventions:  ODIM_H5/V2_2
└── Group: /hel
        Dimensions:            (azimuth: 360, range: 800)
        Coordinates: (12/15)
          * azimuth            (azimuth) float32 1kB 0.5 1.5 2.5 ... 357.5 358.5 359.5
            elevation          (azimuth) float64 3kB dask.array<chunksize=(360,), meta=np.ndarray>
            time               (azimuth) datetime64[ns] 3kB dask.array<chunksize=(360,), meta=np.ndarray>
          * range              (range) float32 3kB 125.0 375.0 ... 1.996e+05 1.999e+05
            x                  (azimuth, range) float64 2MB 2.483e+05 ... 2.564e+05
            y                  (azimuth, range) float64 2MB 5.664e+06 ... 5.863e+06
            ...                 ...
            bins               (azimuth, range) float32 1MB 125.0 375.0 ... 1.999e+05
            sweep_mode         <U20 80B 'azimuth_surveillance'
            longitude          float64 8B 5.406
            latitude           float64 8B 51.07
            altitude           float64 8B 140.0
            crs_wkt            int64 8B 0
        Data variables:
            DBZH               (azimuth, range) float64 2MB dask.array<chunksize=(360, 800), meta=np.ndarray>
            sweep_number       int64 8B ...
            prt_mode           <U7 28B ...
            follow_mode        <U7 28B ...
            sweep_fixed_angle  float64 8B ...
            nyquist_velocity   object 8B ...
        Attributes:
            Conventions:  ODIM_H5/V2_2

xarray.DataTree

• /jab(22)

  • Dimensions:
    • azimuth: 360
    • range: 598
  • Coordinates: (15)
    • azimuth
      (azimuth)
      float32
      0.5 1.5 2.5 ... 357.5 358.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5],shape=(360,), dtype=float32)

    • elevation
      (azimuth)
      float64
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      ray_elevation_angle

      long_name :

      elevation_angle_from_horizontal_plane

      units :

      degrees

      axis :

      radial_elevation_coordinate

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • time
      (azimuth)
      datetime64[ns]
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      time

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type datetime64[ns] numpy.ndarray

    • range
      (range)
      float32
      250.0 750.0 ... 2.982e+05 2.988e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      500.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      250.0

      array([2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05], shape=(598,), dtype=float32)

    • x
      (azimuth, range)
      float64
      8.539e+04 8.543e+04 ... 1.074e+05

      axis :

      X

      long_name :

      Easting

      standard_name :

      projection_x_coordinate

      units :

      metre

      array([[ 85387.99614664,  85432.79968212,  85477.6056995 , ...,112463.48754843, 112509.65992423, 112555.83441459],[ 85392.34570106,  85445.84832822,  85499.35341433, ...,117635.84397661, 117690.68210988, 117745.52226877],[ 85396.68710641,  85458.87252485,  85521.06037625, ...,122798.00995288, 122861.49600322, 122924.98398809],...,[ 85374.91180587,  85393.54669751,  85412.1841216 , ...,96900.90081408,  96920.99478672,  96941.0911264 ],[ 85379.2779378 ,  85406.64508304,  85434.01474695, ...,102094.5248043 , 102123.32262749, 102152.12273611],[ 85383.63978279,  85419.73060544,  85455.82392991, ...,107282.524279  , 107320.01568308, 107357.50928837]],shape=(360, 598))

    • y
      (azimuth, range)
      float64
      5.688e+06 5.689e+06 ... 5.986e+06

      axis :

      Y

      long_name :

      Northing

      standard_name :

      projection_y_coordinate

      units :

      metre

      array([[5688155.89845245, 5688654.73588166, 5689153.57254225, ...,5984708.56943385, 5985206.34421309, 5985704.11619082],[5688155.46951505, 5688653.44900436, 5689151.4276383 , ...,5984183.31118806, 5984680.18138038, 5985177.04870138],[5688154.96472923, 5688651.93458238, 5689148.90349414, ...,5983568.01394866, 5984063.82912611, 5984559.64136594],...,[5688156.72902829, 5688657.22780708, 5689157.72608105, ...,5985742.57902385, 5986242.16218336, 5986741.742771  ],[5688156.52830379, 5688656.62556724, 5689156.72223746, ...,5985488.33309761, 5985987.46460367, 5986486.59345801],[5688156.25141791, 5688655.79484361, 5689155.33758812, ...,5985143.6251097 , 5985642.15374975, 5986140.6796615 ]],shape=(360, 598))

    • z
      (azimuth, range)
      float64
      51.31 53.96 ... 6.849e+03 6.87e+03

      standard_name :

      height_above_ground

      units :

      metre

      array([[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],...,[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573]], shape=(360, 598))

    • gr
      (azimuth, range)
      float64
      250.4 751.3 ... 2.985e+05 2.99e+05

      standard_name :

      distance_from_radar

      units :

      metre

      array([[2.50422917e+02, 7.51268339e+02, 1.25211322e+03, ...,2.98036535e+05, 2.98536446e+05, 2.99036355e+05],[2.50422916e+02, 7.51268304e+02, 1.25211310e+03, ...,2.98029076e+05, 2.98528962e+05, 2.99028846e+05],[2.50422912e+02, 7.51268258e+02, 1.25211297e+03, ...,2.98021692e+05, 2.98521554e+05, 2.99021414e+05],...,[2.50422941e+02, 7.51268513e+02, 1.25211367e+03, ...,2.98059399e+05, 2.98559386e+05, 2.99059370e+05],[2.50422931e+02, 7.51268450e+02, 1.25211351e+03, ...,2.98051701e+05, 2.98551662e+05, 2.99051621e+05],[2.50422922e+02, 7.51268389e+02, 1.25211336e+03, ...,2.98044078e+05, 2.98544014e+05, 2.99043947e+05]], shape=(360, 598))

    • rays
      (azimuth, range)
      float32
      0.5 0.5 0.5 ... 359.5 359.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],[  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],[  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],...,[357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],[358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],[359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]],shape=(360, 598), dtype=float32)

    • bins
      (azimuth, range)
      float32
      250.0 750.0 ... 2.982e+05 2.988e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      500.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      250.0

      array([[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],...,[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05]], shape=(360, 598), dtype=float32)

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      3.064

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(3.0642)

    • latitude
      ()
      float64
      51.19

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(51.1917)

    • altitude
      ()
      float64
      50.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(50.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      WGS 84

      horizontal_datum_name :

      World Geodetic System 1984 ensemble

      projected_crs_name :

      WGS 84 / UTM zone 32N

      grid_mapping_name :

      transverse_mercator

      latitude_of_projection_origin :

      0.0

      longitude_of_central_meridian :

      9.0

      false_easting :

      500000.0

      false_northing :

      0.0

      scale_factor_at_central_meridian :

      0.9996

      array(0)

  • Data variables: (6)
    • DBZH
      (azimuth, range)
      float64
      dask.array<chunksize=(360, 598), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 1.64 MiB 1.64 MiB Shape (360, 598) (360, 598) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

  • Attributes: (1)

    Conventions :

    ODIM_H5/V2_2

  /wid(22)

  • Dimensions:
    • azimuth: 360
    • range: 1000
  • Coordinates: (15)
    • azimuth
      (azimuth)
      float32
      0.5 1.5 2.5 ... 357.5 358.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5],shape=(360,), dtype=float32)

    • elevation
      (azimuth)
      float64
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      ray_elevation_angle

      long_name :

      elevation_angle_from_horizontal_plane

      units :

      degrees

      axis :

      radial_elevation_coordinate

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • time
      (azimuth)
      datetime64[ns]
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      time

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type datetime64[ns] numpy.ndarray

    • range
      (range)
      float32
      125.0 375.0 ... 2.496e+05 2.499e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05, 2.49625e+05,2.49875e+05], shape=(1000,), dtype=float32)

    • x
      (azimuth, range)
      float64
      2.492e+05 2.492e+05 ... 2.588e+05

      axis :

      X

      long_name :

      Easting

      standard_name :

      projection_x_coordinate

      units :

      metre

      array([[249155.06781167, 249168.91818746, 249182.76894182, ...,263144.17662993, 263158.38117039, 263172.58604245],[249157.24559267, 249175.45152815, 249193.65783911, ...,267484.99165246, 267503.5416776 , 267522.09201729],[249159.42060071, 249181.97654959, 249204.53287042, ...,271820.06491209, 271842.95447603, 271865.84433713],...,[249148.52444577, 249149.28809458, 249150.05212842, ...,250100.41962382, 250101.56510125, 250102.71095908],[249150.70660723, 249155.83457771, 249160.96293153, ...,254450.63731808, 254456.13833851, 254461.63972348],[249152.88792831, 249162.37853931, 249171.86953152, ...,258798.95174924, 258808.80618375, 258818.66096652]],shape=(360, 1000))

    • y
      (azimuth, range)
      float64
      5.535e+06 5.535e+06 ... 5.785e+06

      axis :

      Y

      long_name :

      Northing

      standard_name :

      projection_y_coordinate

      units :

      metre

      array([[5535082.98834073, 5535332.67681461, 5535582.36516813, ...,5783891.60265555, 5784140.95708747, 5784390.31097351],[5535082.84846902, 5535332.25718955, 5535581.66577647, ...,5783606.34399779, 5783855.40659879, 5784104.46864295],[5535082.67060929, 5535331.72360051, 5535580.77644495, ...,5783245.44215016, 5783494.13732825, 5783742.8319392 ],...,[5535083.17968139, 5535333.25086672, 5535583.32197181, ...,5784292.75569952, 5784542.53124173, 5784792.30627485],[5535083.15397266, 5535333.17373047, 5535583.1933946 , ...,5784234.86733203, 5784484.57830042, 5784734.28874682],[5535083.09018563, 5535332.98235931, 5535582.87442595, ...,5784101.13355945, 5784350.70413927, 5784600.27418483]],shape=(360, 1000))

    • z
      (azimuth, range)
      float64
      590.7 592.0 ... 5.566e+03 5.575e+03

      standard_name :

      height_above_ground

      units :

      metre

      array([[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],...,[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531]], shape=(360, 1000))

    • gr
      (azimuth, range)
      float64
      125.0 375.1 ... 2.496e+05 2.498e+05

      standard_name :

      distance_from_radar

      units :

      metre

      array([[1.25036200e+02, 3.75108523e+02, 6.25180747e+02, ...,2.49326605e+05, 2.49576363e+05, 2.49826121e+05],[1.25036202e+02, 3.75108523e+02, 6.25180739e+02, ...,2.49323419e+05, 2.49573172e+05, 2.49822923e+05],[1.25036201e+02, 3.75108517e+02, 6.25180721e+02, ...,2.49320273e+05, 2.49570019e+05, 2.49819765e+05],...,[1.25036206e+02, 3.75108558e+02, 6.25180830e+02, ...,2.49336430e+05, 2.49586208e+05, 2.49835986e+05],[1.25036203e+02, 3.75108544e+02, 6.25180798e+02, ...,2.49333113e+05, 2.49582884e+05, 2.49832655e+05],[1.25036201e+02, 3.75108531e+02, 6.25180768e+02, ...,2.49329836e+05, 2.49579601e+05, 2.49829366e+05]],shape=(360, 1000))

    • rays
      (azimuth, range)
      float32
      0.5 0.5 0.5 ... 359.5 359.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],[  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],[  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],...,[357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],[358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],[359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]],shape=(360, 1000), dtype=float32)

    • bins
      (azimuth, range)
      float32
      125.0 375.0 ... 2.496e+05 2.499e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],...,[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05]], shape=(360, 1000), dtype=float32)

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      5.506

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(5.5056)

    • latitude
      ()
      float64
      49.91

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(49.9143)

    • altitude
      ()
      float64
      590.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(590.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      WGS 84

      horizontal_datum_name :

      World Geodetic System 1984 ensemble

      projected_crs_name :

      WGS 84 / UTM zone 32N

      grid_mapping_name :

      transverse_mercator

      latitude_of_projection_origin :

      0.0

      longitude_of_central_meridian :

      9.0

      false_easting :

      500000.0

      false_northing :

      0.0

      scale_factor_at_central_meridian :

      0.9996

      array(0)

  • Data variables: (6)
    • DBZH
      (azimuth, range)
      float64
      dask.array<chunksize=(360, 1000), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 2.75 MiB 2.75 MiB Shape (360, 1000) (360, 1000) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

  • Attributes: (1)

    Conventions :

    ODIM_H5/V2_2

  /hel(22)

  • Dimensions:
    • azimuth: 360
    • range: 800
  • Coordinates: (15)
    • azimuth
      (azimuth)
      float32
      0.5 1.5 2.5 ... 357.5 358.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5],shape=(360,), dtype=float32)

    • elevation
      (azimuth)
      float64
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      ray_elevation_angle

      long_name :

      elevation_angle_from_horizontal_plane

      units :

      degrees

      axis :

      radial_elevation_coordinate

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • time
      (azimuth)
      datetime64[ns]
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      time

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type datetime64[ns] numpy.ndarray

    • range
      (range)
      float32
      125.0 375.0 ... 1.996e+05 1.999e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05, 1.99625e+05,1.99875e+05], shape=(800,), dtype=float32)

    • x
      (azimuth, range)
      float64
      2.483e+05 2.483e+05 ... 2.564e+05

      axis :

      X

      long_name :

      Easting

      standard_name :

      projection_x_coordinate

      units :

      metre

      array([[248254.94425219, 248269.32729914, 248283.71072562, ...,259834.56946303, 259849.23963983, 259863.9101571 ],[248257.1218574 , 248275.8601124 , 248294.59874379, ...,263305.66315785, 263324.68156188, 263343.70029198],[248259.29660856, 248282.38436319, 248305.47249058, ...,266772.06787742, 266795.42847086, 266818.7893756 ],...,[248248.40092683, 248249.69732801, 248250.99411538, ...,249403.67337631, 249405.27576479, 249406.87853456],[248250.58315594, 248256.24401404, 248261.90525657, ...,252882.57112738, 252888.53204489, 252894.49333051],[248252.7644635 , 248262.78793505, 248272.81178882, ...,256359.85310797, 256370.17034974, 256380.48794601]],shape=(360, 800))

    • y
      (azimuth, range)
      float64
      5.664e+06 5.664e+06 ... 5.863e+06

      axis :

      Y

      long_name :

      Northing

      standard_name :

      projection_y_coordinate

      units :

      metre

      array([[5663775.98675203, 5664025.6597276 , 5664275.332581  , ...,5862690.33808305, 5862939.77692429, 5863189.21529733],[5663775.8422333 , 5664025.22616145, 5664274.60995416, ...,5862455.80508824, 5862704.94503751, 5862954.08450717],[5663775.65973034, 5664024.67864266, 5664273.69740635, ...,5862160.76677905, 5862409.53212979, 5862658.29699009],...,[5663776.1920423 , 5664026.27562859, 5664276.35913298, ...,5863030.2826498 , 5863280.16309272, 5863530.04310489],[5663776.16168362, 5664026.18454247, 5664276.20730594, ...,5862977.6249208 , 5863227.43407482, 5863477.24278506],[5663776.09324629, 5664025.97922041, 5664275.86508573, ...,5862864.2969032 , 5863113.95883823, 5863363.620317  ]],shape=(360, 800))

    • z
      (azimuth, range)
      float64
      140.7 142.0 ... 3.532e+03 3.539e+03

      standard_name :

      height_above_ground

      units :

      metre

      array([[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],...,[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903]], shape=(360, 800))

    • gr
      (azimuth, range)
      float64
      125.0 375.1 ... 1.996e+05 1.999e+05

      standard_name :

      distance_from_radar

      units :

      metre

      array([[1.25043498e+02, 3.75130416e+02, 6.25217235e+02, ...,1.99376159e+05, 1.99626028e+05, 1.99875898e+05],[1.25043499e+02, 3.75130416e+02, 6.25217226e+02, ...,1.99374099e+05, 1.99623964e+05, 1.99873828e+05],[1.25043499e+02, 3.75130410e+02, 6.25217208e+02, ...,1.99372060e+05, 1.99621919e+05, 1.99871779e+05],...,[1.25043504e+02, 3.75130451e+02, 6.25217317e+02, ...,1.99382483e+05, 1.99632369e+05, 1.99882254e+05],[1.25043501e+02, 3.75130437e+02, 6.25217285e+02, ...,1.99380353e+05, 1.99630233e+05, 1.99880113e+05],[1.25043499e+02, 3.75130424e+02, 6.25217255e+02, ...,1.99378244e+05, 1.99628118e+05, 1.99877993e+05]], shape=(360, 800))

    • rays
      (azimuth, range)
      float32
      0.5 0.5 0.5 ... 359.5 359.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],[  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],[  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],...,[357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],[358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],[359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]],shape=(360, 800), dtype=float32)

    • bins
      (azimuth, range)
      float32
      125.0 375.0 ... 1.996e+05 1.999e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],...,[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05]], shape=(360, 800), dtype=float32)

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      5.406

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(5.4064)

    • latitude
      ()
      float64
      51.07

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(51.069072)

    • altitude
      ()
      float64
      140.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(140.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      WGS 84

      horizontal_datum_name :

      World Geodetic System 1984 ensemble

      projected_crs_name :

      WGS 84 / UTM zone 32N

      grid_mapping_name :

      transverse_mercator

      latitude_of_projection_origin :

      0.0

      longitude_of_central_meridian :

      9.0

      false_easting :

      500000.0

      false_northing :

      0.0

      scale_factor_at_central_meridian :

      0.9996

      array(0)

  • Data variables: (6)
    • DBZH
      (azimuth, range)
      float64
      dask.array<chunksize=(360, 800), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 2.20 MiB 2.20 MiB Shape (360, 800) (360, 800) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

  • Attributes: (1)

    Conventions :

    ODIM_H5/V2_2

fig, axs = plt.subplots(1, 3, figsize=(18, 5), sharey=True)
ax = axs.flat[0]
ctree["jab"].ds.DBZH.wrl.vis.plot(ax=ax, vmin=0, vmax=60)
ax.set_aspect("equal")
ax.set_title("Radar Jabbeke")
ax = axs.flat[1]
ctree["wid"].ds.DBZH.wrl.vis.plot(ax=ax, vmin=0, vmax=60)
ax.set_aspect("equal")
ax.set_title("Radar Wideumont")
ax = axs.flat[2]
ctree["hel"].ds.DBZH.wrl.vis.plot(ax=ax, vmin=0, vmax=60)
ax.set_aspect("equal")
ax.set_title("Radar Helchteren")
fig.tight_layout()

Calculate Quality

for radar in ctree.children:
    rng = ctree[radar].ds.range
    azimuth = ctree[radar].ds.azimuth
    rscale = rng.diff("range").median().values
    qual = rng.wrl.qual.pulse_volume(rscale, 1.0).expand_dims(azimuth=ctree[radar].ds.azimuth)
    ctree[radar].ds = ctree[radar].ds.assign(QUAL=qual)
display(ctree)

<xarray.DataTree>
Group: /
├── Group: /jab
│       Dimensions:            (azimuth: 360, range: 598)
│       Coordinates: (12/15)
│         * azimuth            (azimuth) float32 1kB 0.5 1.5 2.5 ... 357.5 358.5 359.5
│           elevation          (azimuth) float64 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│           time               (azimuth) datetime64[ns] 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│         * range              (range) float32 2kB 250.0 750.0 ... 2.982e+05 2.988e+05
│           x                  (azimuth, range) float64 2MB 8.539e+04 ... 1.074e+05
│           y                  (azimuth, range) float64 2MB 5.688e+06 ... 5.986e+06
│           ...                 ...
│           bins               (azimuth, range) float32 861kB 250.0 750.0 ... 2.988e+05
│           sweep_mode         <U20 80B 'azimuth_surveillance'
│           longitude          float64 8B 3.064
│           latitude           float64 8B 51.19
│           altitude           float64 8B 50.0
│           crs_wkt            int64 8B 0
│       Data variables:
│           DBZH               (azimuth, range) float64 2MB dask.array<chunksize=(360, 598), meta=np.ndarray>
│           sweep_number       int64 8B ...
│           prt_mode           <U7 28B ...
│           follow_mode        <U7 28B ...
│           sweep_fixed_angle  float64 8B ...
│           nyquist_velocity   object 8B ...
│           QUAL               (azimuth, range) float64 2MB 7.477e+03 ... 1.068e+10
│       Attributes:
│           Conventions:  ODIM_H5/V2_2
├── Group: /wid
│       Dimensions:            (azimuth: 360, range: 1000)
│       Coordinates: (12/15)
│         * azimuth            (azimuth) float32 1kB 0.5 1.5 2.5 ... 357.5 358.5 359.5
│           elevation          (azimuth) float64 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│           time               (azimuth) datetime64[ns] 3kB dask.array<chunksize=(360,), meta=np.ndarray>
│         * range              (range) float32 4kB 125.0 375.0 ... 2.496e+05 2.499e+05
│           x                  (azimuth, range) float64 3MB 2.492e+05 ... 2.588e+05
│           y                  (azimuth, range) float64 3MB 5.535e+06 ... 5.785e+06
│           ...                 ...
│           bins               (azimuth, range) float32 1MB 125.0 375.0 ... 2.499e+05
│           sweep_mode         <U20 80B 'azimuth_surveillance'
│           longitude          float64 8B 5.506
│           latitude           float64 8B 49.91
│           altitude           float64 8B 590.0
│           crs_wkt            int64 8B 0
│       Data variables:
│           DBZH               (azimuth, range) float64 3MB dask.array<chunksize=(360, 1000), meta=np.ndarray>
│           sweep_number       int64 8B ...
│           prt_mode           <U7 28B ...
│           follow_mode        <U7 28B ...
│           sweep_fixed_angle  float64 8B ...
│           nyquist_velocity   object 8B ...
│           QUAL               (azimuth, range) float64 3MB 934.6 ... 3.735e+09
│       Attributes:
│           Conventions:  ODIM_H5/V2_2
└── Group: /hel
        Dimensions:            (azimuth: 360, range: 800)
        Coordinates: (12/15)
          * azimuth            (azimuth) float32 1kB 0.5 1.5 2.5 ... 357.5 358.5 359.5
            elevation          (azimuth) float64 3kB dask.array<chunksize=(360,), meta=np.ndarray>
            time               (azimuth) datetime64[ns] 3kB dask.array<chunksize=(360,), meta=np.ndarray>
          * range              (range) float32 3kB 125.0 375.0 ... 1.996e+05 1.999e+05
            x                  (azimuth, range) float64 2MB 2.483e+05 ... 2.564e+05
            y                  (azimuth, range) float64 2MB 5.664e+06 ... 5.863e+06
            ...                 ...
            bins               (azimuth, range) float32 1MB 125.0 375.0 ... 1.999e+05
            sweep_mode         <U20 80B 'azimuth_surveillance'
            longitude          float64 8B 5.406
            latitude           float64 8B 51.07
            altitude           float64 8B 140.0
            crs_wkt            int64 8B 0
        Data variables:
            DBZH               (azimuth, range) float64 2MB dask.array<chunksize=(360, 800), meta=np.ndarray>
            sweep_number       int64 8B ...
            prt_mode           <U7 28B ...
            follow_mode        <U7 28B ...
            sweep_fixed_angle  float64 8B ...
            nyquist_velocity   object 8B ...
            QUAL               (azimuth, range) float64 2MB 934.6 8.411e+03 ... 2.39e+09
        Attributes:
            Conventions:  ODIM_H5/V2_2

xarray.DataTree

• /jab(23)

  • Dimensions:
    • azimuth: 360
    • range: 598
  • Coordinates: (15)
    • azimuth
      (azimuth)
      float32
      0.5 1.5 2.5 ... 357.5 358.5 359.5

      array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5],shape=(360,), dtype=float32)

    • elevation
      (azimuth)
      float64
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      ray_elevation_angle

      long_name :

      elevation_angle_from_horizontal_plane

      units :

      degrees

      axis :

      radial_elevation_coordinate

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • time
      (azimuth)
      datetime64[ns]
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      time

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type datetime64[ns] numpy.ndarray

    • range
      (range)
      float32
      250.0 750.0 ... 2.982e+05 2.988e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      500.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      250.0

      array([2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05], shape=(598,), dtype=float32)

    • x
      (azimuth, range)
      float64
      8.539e+04 8.543e+04 ... 1.074e+05

      axis :

      X

      long_name :

      Easting

      standard_name :

      projection_x_coordinate

      units :

      metre

      array([[ 85387.99614664,  85432.79968212,  85477.6056995 , ...,112463.48754843, 112509.65992423, 112555.83441459],[ 85392.34570106,  85445.84832822,  85499.35341433, ...,117635.84397661, 117690.68210988, 117745.52226877],[ 85396.68710641,  85458.87252485,  85521.06037625, ...,122798.00995288, 122861.49600322, 122924.98398809],...,[ 85374.91180587,  85393.54669751,  85412.1841216 , ...,96900.90081408,  96920.99478672,  96941.0911264 ],[ 85379.2779378 ,  85406.64508304,  85434.01474695, ...,102094.5248043 , 102123.32262749, 102152.12273611],[ 85383.63978279,  85419.73060544,  85455.82392991, ...,107282.524279  , 107320.01568308, 107357.50928837]],shape=(360, 598))

    • y
      (azimuth, range)
      float64
      5.688e+06 5.689e+06 ... 5.986e+06

      axis :

      Y

      long_name :

      Northing

      standard_name :

      projection_y_coordinate

      units :

      metre

      array([[5688155.89845245, 5688654.73588166, 5689153.57254225, ...,5984708.56943385, 5985206.34421309, 5985704.11619082],[5688155.46951505, 5688653.44900436, 5689151.4276383 , ...,5984183.31118806, 5984680.18138038, 5985177.04870138],[5688154.96472923, 5688651.93458238, 5689148.90349414, ...,5983568.01394866, 5984063.82912611, 5984559.64136594],...,[5688156.72902829, 5688657.22780708, 5689157.72608105, ...,5985742.57902385, 5986242.16218336, 5986741.742771  ],[5688156.52830379, 5688656.62556724, 5689156.72223746, ...,5985488.33309761, 5985987.46460367, 5986486.59345801],[5688156.25141791, 5688655.79484361, 5689155.33758812, ...,5985143.6251097 , 5985642.15374975, 5986140.6796615 ]],shape=(360, 598))

    • z
      (azimuth, range)
      float64
      51.31 53.96 ... 6.849e+03 6.87e+03

      standard_name :

      height_above_ground

      units :

      metre

      array([[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],...,[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573],[  51.31267297,   53.96011096,   56.63700501, ..., 6829.30705617,6849.46524672, 6869.65280573]], shape=(360, 598))

    • gr
      (azimuth, range)
      float64
      250.4 751.3 ... 2.985e+05 2.99e+05

      standard_name :

      distance_from_radar

      units :

      metre

      array([[2.50422917e+02, 7.51268339e+02, 1.25211322e+03, ...,2.98036535e+05, 2.98536446e+05, 2.99036355e+05],[2.50422916e+02, 7.51268304e+02, 1.25211310e+03, ...,2.98029076e+05, 2.98528962e+05, 2.99028846e+05],[2.50422912e+02, 7.51268258e+02, 1.25211297e+03, ...,2.98021692e+05, 2.98521554e+05, 2.99021414e+05],...,[2.50422941e+02, 7.51268513e+02, 1.25211367e+03, ...,2.98059399e+05, 2.98559386e+05, 2.99059370e+05],[2.50422931e+02, 7.51268450e+02, 1.25211351e+03, ...,2.98051701e+05, 2.98551662e+05, 2.99051621e+05],[2.50422922e+02, 7.51268389e+02, 1.25211336e+03, ...,2.98044078e+05, 2.98544014e+05, 2.99043947e+05]], shape=(360, 598))

    • rays
      (azimuth, range)
      float32
      0.5 0.5 0.5 ... 359.5 359.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],[  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],[  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],...,[357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],[358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],[359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]],shape=(360, 598), dtype=float32)

    • bins
      (azimuth, range)
      float32
      250.0 750.0 ... 2.982e+05 2.988e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      500.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      250.0

      array([[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],...,[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05],[2.5000e+02, 7.5000e+02, 1.2500e+03, ..., 2.9775e+05, 2.9825e+05,2.9875e+05]], shape=(360, 598), dtype=float32)

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      3.064

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(3.0642)

    • latitude
      ()
      float64
      51.19

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(51.1917)

    • altitude
      ()
      float64
      50.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(50.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      WGS 84

      horizontal_datum_name :

      World Geodetic System 1984 ensemble

      projected_crs_name :

      WGS 84 / UTM zone 32N

      grid_mapping_name :

      transverse_mercator

      latitude_of_projection_origin :

      0.0

      longitude_of_central_meridian :

      9.0

      false_easting :

      500000.0

      false_northing :

      0.0

      scale_factor_at_central_meridian :

      0.9996

      array(0)

  • Data variables: (7)
    • DBZH
      (azimuth, range)
      float64
      dask.array<chunksize=(360, 598), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 1.64 MiB 1.64 MiB Shape (360, 598) (360, 598) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

    • QUAL
      (azimuth, range)
      float64
      7.477e+03 6.729e+04 ... 1.068e+10

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      500.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      250.0

      array([[7.47681591e+03, 6.72913432e+04, 1.86920398e+05, ...,1.06057211e+10, 1.06413706e+10, 1.06770799e+10],[7.47681591e+03, 6.72913432e+04, 1.86920398e+05, ...,1.06057211e+10, 1.06413706e+10, 1.06770799e+10],[7.47681591e+03, 6.72913432e+04, 1.86920398e+05, ...,1.06057211e+10, 1.06413706e+10, 1.06770799e+10],...,[7.47681591e+03, 6.72913432e+04, 1.86920398e+05, ...,1.06057211e+10, 1.06413706e+10, 1.06770799e+10],[7.47681591e+03, 6.72913432e+04, 1.86920398e+05, ...,1.06057211e+10, 1.06413706e+10, 1.06770799e+10],[7.47681591e+03, 6.72913432e+04, 1.86920398e+05, ...,1.06057211e+10, 1.06413706e+10, 1.06770799e+10]], shape=(360, 598))

  • Attributes: (1)

    Conventions :

    ODIM_H5/V2_2

  /wid(23)

  • Dimensions:
    • azimuth: 360
    • range: 1000
  • Coordinates: (15)
    • azimuth
      (azimuth)
      float32
      0.5 1.5 2.5 ... 357.5 358.5 359.5

      array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5],shape=(360,), dtype=float32)

    • elevation
      (azimuth)
      float64
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      ray_elevation_angle

      long_name :

      elevation_angle_from_horizontal_plane

      units :

      degrees

      axis :

      radial_elevation_coordinate

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • time
      (azimuth)
      datetime64[ns]
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      time

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type datetime64[ns] numpy.ndarray

    • range
      (range)
      float32
      125.0 375.0 ... 2.496e+05 2.499e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05, 2.49625e+05,2.49875e+05], shape=(1000,), dtype=float32)

    • x
      (azimuth, range)
      float64
      2.492e+05 2.492e+05 ... 2.588e+05

      axis :

      X

      long_name :

      Easting

      standard_name :

      projection_x_coordinate

      units :

      metre

      array([[249155.06781167, 249168.91818746, 249182.76894182, ...,263144.17662993, 263158.38117039, 263172.58604245],[249157.24559267, 249175.45152815, 249193.65783911, ...,267484.99165246, 267503.5416776 , 267522.09201729],[249159.42060071, 249181.97654959, 249204.53287042, ...,271820.06491209, 271842.95447603, 271865.84433713],...,[249148.52444577, 249149.28809458, 249150.05212842, ...,250100.41962382, 250101.56510125, 250102.71095908],[249150.70660723, 249155.83457771, 249160.96293153, ...,254450.63731808, 254456.13833851, 254461.63972348],[249152.88792831, 249162.37853931, 249171.86953152, ...,258798.95174924, 258808.80618375, 258818.66096652]],shape=(360, 1000))

    • y
      (azimuth, range)
      float64
      5.535e+06 5.535e+06 ... 5.785e+06

      axis :

      Y

      long_name :

      Northing

      standard_name :

      projection_y_coordinate

      units :

      metre

      array([[5535082.98834073, 5535332.67681461, 5535582.36516813, ...,5783891.60265555, 5784140.95708747, 5784390.31097351],[5535082.84846902, 5535332.25718955, 5535581.66577647, ...,5783606.34399779, 5783855.40659879, 5784104.46864295],[5535082.67060929, 5535331.72360051, 5535580.77644495, ...,5783245.44215016, 5783494.13732825, 5783742.8319392 ],...,[5535083.17968139, 5535333.25086672, 5535583.32197181, ...,5784292.75569952, 5784542.53124173, 5784792.30627485],[5535083.15397266, 5535333.17373047, 5535583.1933946 , ...,5784234.86733203, 5784484.57830042, 5784734.28874682],[5535083.09018563, 5535332.98235931, 5535582.87442595, ...,5784101.13355945, 5784350.70413927, 5784600.27418483]],shape=(360, 1000))

    • z
      (azimuth, range)
      float64
      590.7 592.0 ... 5.566e+03 5.575e+03

      standard_name :

      height_above_ground

      units :

      metre

      array([[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],...,[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531],[ 590.65541586,  591.97176982,  593.2954868 , ..., 5557.4869679 ,5566.13948897, 5574.79923531]], shape=(360, 1000))

    • gr
      (azimuth, range)
      float64
      125.0 375.1 ... 2.496e+05 2.498e+05

      standard_name :

      distance_from_radar

      units :

      metre

      array([[1.25036200e+02, 3.75108523e+02, 6.25180747e+02, ...,2.49326605e+05, 2.49576363e+05, 2.49826121e+05],[1.25036202e+02, 3.75108523e+02, 6.25180739e+02, ...,2.49323419e+05, 2.49573172e+05, 2.49822923e+05],[1.25036201e+02, 3.75108517e+02, 6.25180721e+02, ...,2.49320273e+05, 2.49570019e+05, 2.49819765e+05],...,[1.25036206e+02, 3.75108558e+02, 6.25180830e+02, ...,2.49336430e+05, 2.49586208e+05, 2.49835986e+05],[1.25036203e+02, 3.75108544e+02, 6.25180798e+02, ...,2.49333113e+05, 2.49582884e+05, 2.49832655e+05],[1.25036201e+02, 3.75108531e+02, 6.25180768e+02, ...,2.49329836e+05, 2.49579601e+05, 2.49829366e+05]],shape=(360, 1000))

    • rays
      (azimuth, range)
      float32
      0.5 0.5 0.5 ... 359.5 359.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],[  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],[  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],...,[357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],[358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],[359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]],shape=(360, 1000), dtype=float32)

    • bins
      (azimuth, range)
      float32
      125.0 375.0 ... 2.496e+05 2.499e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],...,[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 2.49375e+05,2.49625e+05, 2.49875e+05]], shape=(360, 1000), dtype=float32)

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      5.506

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(5.5056)

    • latitude
      ()
      float64
      49.91

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(49.9143)

    • altitude
      ()
      float64
      590.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(590.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      WGS 84

      horizontal_datum_name :

      World Geodetic System 1984 ensemble

      projected_crs_name :

      WGS 84 / UTM zone 32N

      grid_mapping_name :

      transverse_mercator

      latitude_of_projection_origin :

      0.0

      longitude_of_central_meridian :

      9.0

      false_easting :

      500000.0

      false_northing :

      0.0

      scale_factor_at_central_meridian :

      0.9996

      array(0)

  • Data variables: (7)
    • DBZH
      (azimuth, range)
      float64
      dask.array<chunksize=(360, 1000), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 2.75 MiB 2.75 MiB Shape (360, 1000) (360, 1000) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

    • QUAL
      (azimuth, range)
      float64
      934.6 8.411e+03 ... 3.735e+09

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,3.71973916e+09, 3.72720113e+09, 3.73467046e+09],[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,3.71973916e+09, 3.72720113e+09, 3.73467046e+09],[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,3.71973916e+09, 3.72720113e+09, 3.73467046e+09],...,[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,3.71973916e+09, 3.72720113e+09, 3.73467046e+09],[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,3.71973916e+09, 3.72720113e+09, 3.73467046e+09],[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,3.71973916e+09, 3.72720113e+09, 3.73467046e+09]],shape=(360, 1000))

  • Attributes: (1)

    Conventions :

    ODIM_H5/V2_2

  /hel(23)

  • Dimensions:
    • azimuth: 360
    • range: 800
  • Coordinates: (15)
    • azimuth
      (azimuth)
      float32
      0.5 1.5 2.5 ... 357.5 358.5 359.5

      array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5],shape=(360,), dtype=float32)

    • elevation
      (azimuth)
      float64
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      ray_elevation_angle

      long_name :

      elevation_angle_from_horizontal_plane

      units :

      degrees

      axis :

      radial_elevation_coordinate

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • time
      (azimuth)
      datetime64[ns]
      dask.array<chunksize=(360,), meta=np.ndarray>

      standard_name :

      time

      Array Chunk Bytes 2.81 kiB 2.81 kiB Shape (360,) (360,) Dask graph 1 chunks in 2 graph layers Data type datetime64[ns] numpy.ndarray

    • range
      (range)
      float32
      125.0 375.0 ... 1.996e+05 1.999e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05, 1.99625e+05,1.99875e+05], shape=(800,), dtype=float32)

    • x
      (azimuth, range)
      float64
      2.483e+05 2.483e+05 ... 2.564e+05

      axis :

      X

      long_name :

      Easting

      standard_name :

      projection_x_coordinate

      units :

      metre

      array([[248254.94425219, 248269.32729914, 248283.71072562, ...,259834.56946303, 259849.23963983, 259863.9101571 ],[248257.1218574 , 248275.8601124 , 248294.59874379, ...,263305.66315785, 263324.68156188, 263343.70029198],[248259.29660856, 248282.38436319, 248305.47249058, ...,266772.06787742, 266795.42847086, 266818.7893756 ],...,[248248.40092683, 248249.69732801, 248250.99411538, ...,249403.67337631, 249405.27576479, 249406.87853456],[248250.58315594, 248256.24401404, 248261.90525657, ...,252882.57112738, 252888.53204489, 252894.49333051],[248252.7644635 , 248262.78793505, 248272.81178882, ...,256359.85310797, 256370.17034974, 256380.48794601]],shape=(360, 800))

    • y
      (azimuth, range)
      float64
      5.664e+06 5.664e+06 ... 5.863e+06

      axis :

      Y

      long_name :

      Northing

      standard_name :

      projection_y_coordinate

      units :

      metre

      array([[5663775.98675203, 5664025.6597276 , 5664275.332581  , ...,5862690.33808305, 5862939.77692429, 5863189.21529733],[5663775.8422333 , 5664025.22616145, 5664274.60995416, ...,5862455.80508824, 5862704.94503751, 5862954.08450717],[5663775.65973034, 5664024.67864266, 5664273.69740635, ...,5862160.76677905, 5862409.53212979, 5862658.29699009],...,[5663776.1920423 , 5664026.27562859, 5664276.35913298, ...,5863030.2826498 , 5863280.16309272, 5863530.04310489],[5663776.16168362, 5664026.18454247, 5664276.20730594, ...,5862977.6249208 , 5863227.43407482, 5863477.24278506],[5663776.09324629, 5664025.97922041, 5664275.86508573, ...,5862864.2969032 , 5863113.95883823, 5863363.620317  ]],shape=(360, 800))

    • z
      (azimuth, range)
      float64
      140.7 142.0 ... 3.532e+03 3.539e+03

      standard_name :

      height_above_ground

      units :

      metre

      array([[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],...,[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903],[ 140.65541597,  141.97177081,  143.29548955, ..., 3525.05038896,3532.23305328, 3539.42318903]], shape=(360, 800))

    • gr
      (azimuth, range)
      float64
      125.0 375.1 ... 1.996e+05 1.999e+05

      standard_name :

      distance_from_radar

      units :

      metre

      array([[1.25043498e+02, 3.75130416e+02, 6.25217235e+02, ...,1.99376159e+05, 1.99626028e+05, 1.99875898e+05],[1.25043499e+02, 3.75130416e+02, 6.25217226e+02, ...,1.99374099e+05, 1.99623964e+05, 1.99873828e+05],[1.25043499e+02, 3.75130410e+02, 6.25217208e+02, ...,1.99372060e+05, 1.99621919e+05, 1.99871779e+05],...,[1.25043504e+02, 3.75130451e+02, 6.25217317e+02, ...,1.99382483e+05, 1.99632369e+05, 1.99882254e+05],[1.25043501e+02, 3.75130437e+02, 6.25217285e+02, ...,1.99380353e+05, 1.99630233e+05, 1.99880113e+05],[1.25043499e+02, 3.75130424e+02, 6.25217255e+02, ...,1.99378244e+05, 1.99628118e+05, 1.99877993e+05]], shape=(360, 800))

    • rays
      (azimuth, range)
      float32
      0.5 0.5 0.5 ... 359.5 359.5 359.5

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],[  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],[  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],...,[357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],[358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],[359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]],shape=(360, 800), dtype=float32)

    • bins
      (azimuth, range)
      float32
      125.0 375.0 ... 1.996e+05 1.999e+05

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],...,[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05],[1.25000e+02, 3.75000e+02, 6.25000e+02, ..., 1.99375e+05,1.99625e+05, 1.99875e+05]], shape=(360, 800), dtype=float32)

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      5.406

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(5.4064)

    • latitude
      ()
      float64
      51.07

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(51.069072)

    • altitude
      ()
      float64
      140.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(140.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      WGS 84

      horizontal_datum_name :

      World Geodetic System 1984 ensemble

      projected_crs_name :

      WGS 84 / UTM zone 32N

      grid_mapping_name :

      transverse_mercator

      latitude_of_projection_origin :

      0.0

      longitude_of_central_meridian :

      9.0

      false_easting :

      500000.0

      false_northing :

      0.0

      scale_factor_at_central_meridian :

      0.9996

      array(0)

  • Data variables: (7)
    • DBZH
      (azimuth, range)
      float64
      dask.array<chunksize=(360, 800), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 2.20 MiB 2.20 MiB Shape (360, 800) (360, 800) Dask graph 1 chunks in 2 graph layers Data type float64 numpy.ndarray

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

    • QUAL
      (azimuth, range)
      float64
      934.6 8.411e+03 ... 2.39e+09

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,2.37765077e+09, 2.38361726e+09, 2.38959135e+09],[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,2.37765077e+09, 2.38361726e+09, 2.38959135e+09],[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,2.37765077e+09, 2.38361726e+09, 2.38959135e+09],...,[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,2.37765077e+09, 2.38361726e+09, 2.38959135e+09],[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,2.37765077e+09, 2.38361726e+09, 2.38959135e+09],[9.34601988e+02, 8.41141789e+03, 2.33650497e+04, ...,2.37765077e+09, 2.38361726e+09, 2.38959135e+09]], shape=(360, 800))

  • Attributes: (1)

    Conventions :

    ODIM_H5/V2_2

Gridding

First, we create a Datset with cartesian coordinates, which can hold our three radars. We are using wradlib.comp.CompMethods.togrid to interpolate our sweep data to the cartesian domain. As an interpolator we use wradlib.ipol.Nearest.

xmin, xmax, ymin, ymax = ctree.wrl.util.bbox()
x = np.linspace(xmin, xmax + 1000.0, 1000)
y = np.linspace(ymin, ymax + 1000.0, 1000)
cart = xr.Dataset(coords={"x": (["x"], x), "y": (["y"], y)})
display(cart)

<xarray.Dataset> Size: 16kB
Dimensions:  (x: 1000, y: 1000)
Coordinates:
  * x        (x) float64 8kB -2.143e+05 -2.136e+05 ... 4.991e+05 4.999e+05
  * y        (y) float64 8kB 5.285e+06 5.286e+06 ... 5.987e+06 5.988e+06
Data variables:
    *empty*

xarray.Dataset

• Dimensions:
  • x: 1000
  • y: 1000
• Coordinates: (2)
  • x
    (x)
    float64
    -2.143e+05 -2.136e+05 ... 4.999e+05

    array([-214275.862723, -213561.016705, -212846.170687, ...,  498425.617206,
            499140.463224,  499855.309242], shape=(1000,))

  • y
    (y)
    float64
    5.285e+06 5.286e+06 ... 5.988e+06

    array([5285121.537223, 5285825.098598, 5286528.659973, ..., 5986572.228301,
           5987275.789676, 5987979.351051], shape=(1000,))

• Data variables: (0)

gtree = xr.DataTree()
for radar in ctree.children:
    gridded = ctree[radar].ds.wrl.comp.togrid(cart)
    gtree[radar] = gridded
display(gtree)

<xarray.DataTree>
Group: /
├── Group: /jab
│       Dimensions:            (y: 1000, x: 1000)
│       Coordinates:
│         * y                  (y) float64 8kB 5.285e+06 5.286e+06 ... 5.988e+06
│         * x                  (x) float64 8kB -2.143e+05 -2.136e+05 ... 4.999e+05
│           sweep_mode         <U20 80B 'azimuth_surveillance'
│           longitude          float64 8B 3.064
│           latitude           float64 8B 51.19
│           altitude           float64 8B 50.0
│           crs_wkt            int64 8B 0
│       Data variables:
│           DBZH               (y, x) float64 8MB dask.array<chunksize=(1000, 1000), meta=np.ndarray>
│           QUAL               (y, x) float64 8MB nan nan nan nan ... nan nan nan nan
│           sweep_number       int64 8B ...
│           prt_mode           <U7 28B ...
│           follow_mode        <U7 28B ...
│           sweep_fixed_angle  float64 8B ...
│           nyquist_velocity   object 8B ...
├── Group: /wid
│       Dimensions:            (y: 1000, x: 1000)
│       Coordinates:
│         * y                  (y) float64 8kB 5.285e+06 5.286e+06 ... 5.988e+06
│         * x                  (x) float64 8kB -2.143e+05 -2.136e+05 ... 4.999e+05
│           sweep_mode         <U20 80B 'azimuth_surveillance'
│           longitude          float64 8B 5.506
│           latitude           float64 8B 49.91
│           altitude           float64 8B 590.0
│           crs_wkt            int64 8B 0
│       Data variables:
│           DBZH               (y, x) float64 8MB dask.array<chunksize=(1000, 1000), meta=np.ndarray>
│           QUAL               (y, x) float64 8MB nan nan nan nan ... nan nan nan nan
│           sweep_number       int64 8B ...
│           prt_mode           <U7 28B ...
│           follow_mode        <U7 28B ...
│           sweep_fixed_angle  float64 8B ...
│           nyquist_velocity   object 8B ...
└── Group: /hel
        Dimensions:            (y: 1000, x: 1000)
        Coordinates:
          * y                  (y) float64 8kB 5.285e+06 5.286e+06 ... 5.988e+06
          * x                  (x) float64 8kB -2.143e+05 -2.136e+05 ... 4.999e+05
            sweep_mode         <U20 80B 'azimuth_surveillance'
            longitude          float64 8B 5.406
            latitude           float64 8B 51.07
            altitude           float64 8B 140.0
            crs_wkt            int64 8B 0
        Data variables:
            DBZH               (y, x) float64 8MB dask.array<chunksize=(1000, 1000), meta=np.ndarray>
            QUAL               (y, x) float64 8MB nan nan nan nan ... nan nan nan nan
            sweep_number       int64 8B ...
            prt_mode           <U7 28B ...
            follow_mode        <U7 28B ...
            sweep_fixed_angle  float64 8B ...
            nyquist_velocity   object 8B ...

xarray.DataTree

• /jab(14)

  • Dimensions:
    • y: 1000
    • x: 1000
  • Coordinates: (7)
    • y
      (y)
      float64
      5.285e+06 5.286e+06 ... 5.988e+06

      array([5285121.537223, 5285825.098598, 5286528.659973, ..., 5986572.228301,5987275.789676, 5987979.351051], shape=(1000,))

    • x
      (x)
      float64
      -2.143e+05 -2.136e+05 ... 4.999e+05

      array([-214275.862723, -213561.016705, -212846.170687, ...,  498425.617206,499140.463224,  499855.309242], shape=(1000,))

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      3.064

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(3.0642)

    • latitude
      ()
      float64
      51.19

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(51.1917)

    • altitude
      ()
      float64
      50.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(50.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      WGS 84

      horizontal_datum_name :

      World Geodetic System 1984 ensemble

      projected_crs_name :

      WGS 84 / UTM zone 32N

      grid_mapping_name :

      transverse_mercator

      latitude_of_projection_origin :

      0.0

      longitude_of_central_meridian :

      9.0

      false_easting :

      500000.0

      false_northing :

      0.0

      scale_factor_at_central_meridian :

      0.9996

      array(0)

  • Data variables: (7)
    • DBZH
      (y, x)
      float64
      dask.array<chunksize=(1000, 1000), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 7.63 MiB 7.63 MiB Shape (1000, 1000) (1000, 1000) Dask graph 1 chunks in 9 graph layers Data type float64 numpy.ndarray

    • QUAL
      (y, x)
      float64
      nan nan nan nan ... nan nan nan nan

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      500.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      250.0

      array([[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],...,[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan]], shape=(1000, 1000))

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

  /wid(14)

  • Dimensions:
    • y: 1000
    • x: 1000
  • Coordinates: (7)
    • y
      (y)
      float64
      5.285e+06 5.286e+06 ... 5.988e+06

      array([5285121.537223, 5285825.098598, 5286528.659973, ..., 5986572.228301,5987275.789676, 5987979.351051], shape=(1000,))

    • x
      (x)
      float64
      -2.143e+05 -2.136e+05 ... 4.999e+05

      array([-214275.862723, -213561.016705, -212846.170687, ...,  498425.617206,499140.463224,  499855.309242], shape=(1000,))

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      5.506

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(5.5056)

    • latitude
      ()
      float64
      49.91

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(49.9143)

    • altitude
      ()
      float64
      590.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(590.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      WGS 84

      horizontal_datum_name :

      World Geodetic System 1984 ensemble

      projected_crs_name :

      WGS 84 / UTM zone 32N

      grid_mapping_name :

      transverse_mercator

      latitude_of_projection_origin :

      0.0

      longitude_of_central_meridian :

      9.0

      false_easting :

      500000.0

      false_northing :

      0.0

      scale_factor_at_central_meridian :

      0.9996

      array(0)

  • Data variables: (7)
    • DBZH
      (y, x)
      float64
      dask.array<chunksize=(1000, 1000), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 7.63 MiB 7.63 MiB Shape (1000, 1000) (1000, 1000) Dask graph 1 chunks in 9 graph layers Data type float64 numpy.ndarray

    • QUAL
      (y, x)
      float64
      nan nan nan nan ... nan nan nan nan

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],...,[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan]], shape=(1000, 1000))

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

  /hel(14)

  • Dimensions:
    • y: 1000
    • x: 1000
  • Coordinates: (7)
    • y
      (y)
      float64
      5.285e+06 5.286e+06 ... 5.988e+06

      array([5285121.537223, 5285825.098598, 5286528.659973, ..., 5986572.228301,5987275.789676, 5987979.351051], shape=(1000,))

    • x
      (x)
      float64
      -2.143e+05 -2.136e+05 ... 4.999e+05

      array([-214275.862723, -213561.016705, -212846.170687, ...,  498425.617206,499140.463224,  499855.309242], shape=(1000,))

    • sweep_mode
      ()
      <U20
      'azimuth_surveillance'

      array('azimuth_surveillance', dtype='<U20')

    • longitude
      ()
      float64
      5.406

      long_name :

      longitude

      units :

      degrees_east

      standard_name :

      longitude

      array(5.4064)

    • latitude
      ()
      float64
      51.07

      long_name :

      latitude

      units :

      degrees_north

      positive :

      up

      standard_name :

      latitude

      array(51.069072)

    • altitude
      ()
      float64
      140.0

      long_name :

      altitude

      units :

      meters

      standard_name :

      altitude

      array(140.)

    • crs_wkt
      ()
      int64
      0

      crs_wkt :

      PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

      semi_major_axis :

      6378137.0

      semi_minor_axis :

      6356752.314245179

      inverse_flattening :

      298.257223563

      reference_ellipsoid_name :

      WGS 84

      longitude_of_prime_meridian :

      0.0

      prime_meridian_name :

      Greenwich

      geographic_crs_name :

      WGS 84

      horizontal_datum_name :

      World Geodetic System 1984 ensemble

      projected_crs_name :

      WGS 84 / UTM zone 32N

      grid_mapping_name :

      transverse_mercator

      latitude_of_projection_origin :

      0.0

      longitude_of_central_meridian :

      9.0

      false_easting :

      500000.0

      false_northing :

      0.0

      scale_factor_at_central_meridian :

      0.9996

      array(0)

  • Data variables: (7)
    • DBZH
      (y, x)
      float64
      dask.array<chunksize=(1000, 1000), meta=np.ndarray>

      _Undetect :

      0.0

      long_name :

      Equivalent reflectivity factor H

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      Array Chunk Bytes 7.63 MiB 7.63 MiB Shape (1000, 1000) (1000, 1000) Dask graph 1 chunks in 9 graph layers Data type float64 numpy.ndarray

    • QUAL
      (y, x)
      float64
      nan nan nan nan ... nan nan nan nan

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      250.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      125.0

      array([[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],...,[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan]], shape=(1000, 1000))

    • sweep_number
      ()
      int64
      ...

      [1 values with dtype=int64]

    • prt_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • follow_mode
      ()
      <U7
      ...

      [1 values with dtype=<U7]

    • sweep_fixed_angle
      ()
      float64
      ...

      [1 values with dtype=float64]

    • nyquist_velocity
      ()
      object
      ...

      standard_name :

      nyquist_velocity

      units :

      m s-1

      [1 values with dtype=object]

fig, axs = plt.subplots(1, 3, figsize=(18, 5), sharey=True)
ax = axs.flat[0]
gtree["jab"].ds.DBZH.plot(ax=ax, vmin=0, vmax=60, cmap="HomeyerRainbow")
ax.set_aspect("equal")
ax.set_title("Radar Jabbeke")
ax = axs.flat[1]
gtree["wid"].ds.DBZH.plot(ax=ax, vmin=0, vmax=60, cmap="HomeyerRainbow")
ax.set_aspect("equal")
ax.set_title("Radar Wideumont")
ax = axs.flat[2]
gtree["hel"].ds.DBZH.plot(ax=ax, vmin=0, vmax=60, cmap="HomeyerRainbow")
ax.set_aspect("equal")
ax.set_title("Radar Helchteren")
fig.tight_layout()

fig, axs = plt.subplots(1, 3, figsize=(18, 5), sharey=True)
ax = axs.flat[0]
gtree["jab"].ds.QUAL.plot(ax=ax)
ax.set_aspect("equal")
ax.set_title("Radar Jabbeke")
ax = axs.flat[1]
gtree["wid"].ds.QUAL.plot(ax=ax)
ax.set_aspect("equal")
ax.set_title("Radar Wideumont")
ax = axs.flat[2]
gtree["hel"].ds.QUAL.plot(ax=ax)
ax.set_aspect("equal")
ax.set_title("Radar Helchteren")
fig.tight_layout()

Compositing

Before compositing we combine the three radar grids as well as the quality grids into one Dataset, respectively.

radars = xr.DataArray(gtree.children, dims="radar")
radargrids = xr.concat([gtree[radar].ds.DBZH for radar in gtree.children], dim=radars)
# normalizing Quality between 1. and 0.
qualitygrids = xr.concat([1. - (gtree[radar].ds.QUAL / gtree[radar].ds.QUAL.max()) for radar in gtree.children], dim=radars)

display(radargrids)
display(qualitygrids)

<xarray.DataArray 'DBZH' (radar: 3, y: 1000, x: 1000)> Size: 24MB
dask.array<concatenate, shape=(3, 1000, 1000), dtype=float64, chunksize=(1, 1000, 1000), chunktype=numpy.ndarray>
Coordinates:
  * radar       (radar) <U3 36B 'jab' 'wid' 'hel'
    longitude   (radar) float64 24B 3.064 5.506 5.406
    latitude    (radar) float64 24B 51.19 49.91 51.07
    altitude    (radar) float64 24B 50.0 590.0 140.0
  * y           (y) float64 8kB 5.285e+06 5.286e+06 ... 5.987e+06 5.988e+06
  * x           (x) float64 8kB -2.143e+05 -2.136e+05 ... 4.991e+05 4.999e+05
    sweep_mode  <U20 80B 'azimuth_surveillance'
    crs_wkt     int64 8B 0
Attributes:
    _Undetect:      0.0
    long_name:      Equivalent reflectivity factor H
    standard_name:  radar_equivalent_reflectivity_factor_h
    units:          dBZ

xarray.DataArray

'DBZH'

• radar: 3
• y: 1000
• x: 1000

• dask.array<chunksize=(1, 1000, 1000), meta=np.ndarray>

  Array Chunk Bytes 22.89 MiB 7.63 MiB Shape (3, 1000, 1000) (1, 1000, 1000) Dask graph 3 chunks in 29 graph layers Data type float64 numpy.ndarray

• Coordinates: (8)
  • radar
    (radar)
    <U3
    'jab' 'wid' 'hel'

    array(['jab', 'wid', 'hel'], dtype='<U3')

  • longitude
    (radar)
    float64
    3.064 5.506 5.406

    long_name :

    longitude

    units :

    degrees_east

    standard_name :

    longitude

    array([3.0642, 5.5056, 5.4064])

  • latitude
    (radar)
    float64
    51.19 49.91 51.07

    long_name :

    latitude

    units :

    degrees_north

    positive :

    up

    standard_name :

    latitude

    array([51.1917  , 49.9143  , 51.069072])

  • altitude
    (radar)
    float64
    50.0 590.0 140.0

    long_name :

    altitude

    units :

    meters

    standard_name :

    altitude

    array([ 50., 590., 140.])

  • y
    (y)
    float64
    5.285e+06 5.286e+06 ... 5.988e+06

    array([5285121.537223, 5285825.098598, 5286528.659973, ..., 5986572.228301,
           5987275.789676, 5987979.351051], shape=(1000,))

  • x
    (x)
    float64
    -2.143e+05 -2.136e+05 ... 4.999e+05

    array([-214275.862723, -213561.016705, -212846.170687, ...,  498425.617206,
            499140.463224,  499855.309242], shape=(1000,))

  • sweep_mode
    ()
    <U20
    'azimuth_surveillance'

    array('azimuth_surveillance', dtype='<U20')

  • crs_wkt
    ()
    int64
    0

    crs_wkt :

    PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

    semi_major_axis :

    6378137.0

    semi_minor_axis :

    6356752.314245179

    inverse_flattening :

    298.257223563

    reference_ellipsoid_name :

    WGS 84

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    WGS 84

    horizontal_datum_name :

    World Geodetic System 1984 ensemble

    projected_crs_name :

    WGS 84 / UTM zone 32N

    grid_mapping_name :

    transverse_mercator

    latitude_of_projection_origin :

    0.0

    longitude_of_central_meridian :

    9.0

    false_easting :

    500000.0

    false_northing :

    0.0

    scale_factor_at_central_meridian :

    0.9996

    array(0)

• Attributes: (4)

  _Undetect :

  0.0

  long_name :

  Equivalent reflectivity factor H

  standard_name :

  radar_equivalent_reflectivity_factor_h

  units :

  dBZ

<xarray.DataArray 'QUAL' (radar: 3, y: 1000, x: 1000)> Size: 24MB
array([[[nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan],
        ...,
        [nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan]],

       [[nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan],
        ...,
        [nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan]],

       [[nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan],
        ...,
        [nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan],
        [nan, nan, nan, ..., nan, nan, nan]]], shape=(3, 1000, 1000))
Coordinates:
  * radar       (radar) <U3 36B 'jab' 'wid' 'hel'
    longitude   (radar) float64 24B 3.064 5.506 5.406
    latitude    (radar) float64 24B 51.19 49.91 51.07
    altitude    (radar) float64 24B 50.0 590.0 140.0
  * y           (y) float64 8kB 5.285e+06 5.286e+06 ... 5.987e+06 5.988e+06
  * x           (x) float64 8kB -2.143e+05 -2.136e+05 ... 4.991e+05 4.999e+05
    sweep_mode  <U20 80B 'azimuth_surveillance'
    crs_wkt     int64 8B 0
Attributes:
    units:                           meters
    standard_name:                   projection_range_coordinate
    long_name:                       range_to_measurement_volume
    axis:                            radial_range_coordinate
    meters_between_gates:            500.0
    spacing_is_constant:             true
    meters_to_center_of_first_gate:  250.0

xarray.DataArray

'QUAL'

• radar: 3
• y: 1000
• x: 1000

• nan nan nan nan nan nan nan nan ... nan nan nan nan nan nan nan nan

  array([[[nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan],
          ...,
          [nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan]],
  
         [[nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan],
          ...,
          [nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan]],
  
         [[nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan],
          ...,
          [nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan],
          [nan, nan, nan, ..., nan, nan, nan]]], shape=(3, 1000, 1000))

• Coordinates: (8)
  • radar
    (radar)
    <U3
    'jab' 'wid' 'hel'

    array(['jab', 'wid', 'hel'], dtype='<U3')

  • longitude
    (radar)
    float64
    3.064 5.506 5.406

    long_name :

    longitude

    units :

    degrees_east

    standard_name :

    longitude

    array([3.0642, 5.5056, 5.4064])

  • latitude
    (radar)
    float64
    51.19 49.91 51.07

    long_name :

    latitude

    units :

    degrees_north

    positive :

    up

    standard_name :

    latitude

    array([51.1917  , 49.9143  , 51.069072])

  • altitude
    (radar)
    float64
    50.0 590.0 140.0

    long_name :

    altitude

    units :

    meters

    standard_name :

    altitude

    array([ 50., 590., 140.])

  • y
    (y)
    float64
    5.285e+06 5.286e+06 ... 5.988e+06

    array([5285121.537223, 5285825.098598, 5286528.659973, ..., 5986572.228301,
           5987275.789676, 5987979.351051], shape=(1000,))

  • x
    (x)
    float64
    -2.143e+05 -2.136e+05 ... 4.999e+05

    array([-214275.862723, -213561.016705, -212846.170687, ...,  498425.617206,
            499140.463224,  499855.309242], shape=(1000,))

  • sweep_mode
    ()
    <U20
    'azimuth_surveillance'

    array('azimuth_surveillance', dtype='<U20')

  • crs_wkt
    ()
    int64
    0

    crs_wkt :

    PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

    semi_major_axis :

    6378137.0

    semi_minor_axis :

    6356752.314245179

    inverse_flattening :

    298.257223563

    reference_ellipsoid_name :

    WGS 84

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    WGS 84

    horizontal_datum_name :

    World Geodetic System 1984 ensemble

    projected_crs_name :

    WGS 84 / UTM zone 32N

    grid_mapping_name :

    transverse_mercator

    latitude_of_projection_origin :

    0.0

    longitude_of_central_meridian :

    9.0

    false_easting :

    500000.0

    false_northing :

    0.0

    scale_factor_at_central_meridian :

    0.9996

    array(0)

• Attributes: (7)

  units :

  meters

  standard_name :

  projection_range_coordinate

  long_name :

  range_to_measurement_volume

  axis :

  radial_range_coordinate

  meters_between_gates :

  500.0

  spacing_is_constant :

  true

  meters_to_center_of_first_gate :

  250.0

Then we finally can call wradlib.comp.CompMethods.compose_weighted to create the final output.

composite = radargrids.wrl.comp.compose_weighted(qualitygrids)
display(composite)

<xarray.DataArray 'DBZH' (y: 1000, x: 1000)> Size: 8MB
dask.array<where, shape=(1000, 1000), dtype=float64, chunksize=(1000, 1000), chunktype=numpy.ndarray>
Coordinates:
  * y           (y) float64 8kB 5.285e+06 5.286e+06 ... 5.987e+06 5.988e+06
  * x           (x) float64 8kB -2.143e+05 -2.136e+05 ... 4.991e+05 4.999e+05
    sweep_mode  <U20 80B 'azimuth_surveillance'
    crs_wkt     int64 8B 0
Attributes:
    _Undetect:      0.0
    long_name:      Equivalent reflectivity factor H
    standard_name:  radar_equivalent_reflectivity_factor_h
    units:          dBZ

xarray.DataArray

'DBZH'

• y: 1000
• x: 1000

• dask.array<chunksize=(1000, 1000), meta=np.ndarray>

  Array Chunk Bytes 7.63 MiB 7.63 MiB Shape (1000, 1000) (1000, 1000) Dask graph 1 chunks in 51 graph layers Data type float64 numpy.ndarray

• Coordinates: (4)
  • y
    (y)
    float64
    5.285e+06 5.286e+06 ... 5.988e+06

    array([5285121.537223, 5285825.098598, 5286528.659973, ..., 5986572.228301,
           5987275.789676, 5987979.351051], shape=(1000,))

  • x
    (x)
    float64
    -2.143e+05 -2.136e+05 ... 4.999e+05

    array([-214275.862723, -213561.016705, -212846.170687, ...,  498425.617206,
            499140.463224,  499855.309242], shape=(1000,))

  • sweep_mode
    ()
    <U20
    'azimuth_surveillance'

    array('azimuth_surveillance', dtype='<U20')

  • crs_wkt
    ()
    int64
    0

    crs_wkt :

    PROJCRS["WGS 84 / UTM zone 32N",BASEGEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],MEMBER["World Geodetic System 1984 (G2296)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 32N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Navigation and medium accuracy spatial referencing."],AREA["Between 6°E and 12°E, northern hemisphere between equator and 84°N, onshore and offshore. Algeria. Austria. Cameroon. Denmark. Equatorial Guinea. France. Gabon. Germany. Italy. Libya. Liechtenstein. Monaco. Netherlands. Niger. Nigeria. Norway. Sao Tome and Principe. Svalbard. Sweden. Switzerland. Tunisia. Vatican City State."],BBOX[0,6,84,12]],ID["EPSG",32632]]

    semi_major_axis :

    6378137.0

    semi_minor_axis :

    6356752.314245179

    inverse_flattening :

    298.257223563

    reference_ellipsoid_name :

    WGS 84

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    WGS 84

    horizontal_datum_name :

    World Geodetic System 1984 ensemble

    projected_crs_name :

    WGS 84 / UTM zone 32N

    grid_mapping_name :

    transverse_mercator

    latitude_of_projection_origin :

    0.0

    longitude_of_central_meridian :

    9.0

    false_easting :

    500000.0

    false_northing :

    0.0

    scale_factor_at_central_meridian :

    0.9996

    array(0)

• Attributes: (4)

  _Undetect :

  0.0

  long_name :

  Equivalent reflectivity factor H

  standard_name :

  radar_equivalent_reflectivity_factor_h

  units :

  dBZ

Plot Result

composite.where(composite>0.1).plot(cmap="HomeyerRainbow", vmin=0, vmax=60)

<matplotlib.collections.QuadMesh at 0x79012c35bed0>