xarray CfRadial1 backend#

In this example, we read CfRadial1 data files using the xradar cfradial1 backend.

Data is also exported to ODIM_H5 and CfRadial2.

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

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

warnings.filterwarnings("ignore")

Load CfRadial1 Volume Data#

We use the functionality provided now by xradar to read the CfRadial1 data into a DataTree.

fpath = "netcdf/cfrad.20080604_002217_000_SPOL_v36_SUR.nc"
f = wradlib_data.DATASETS.fetch(fpath)
vol = xd.io.open_cfradial1_datatree(f)
Downloading file 'netcdf/cfrad.20080604_002217_000_SPOL_v36_SUR.nc' from 'https://github.com/wradlib/wradlib-data/raw/main/data/netcdf/cfrad.20080604_002217_000_SPOL_v36_SUR.nc' to '/home/docs/.cache/wradlib-data'.

Inspect RadarVolume#

display(vol)
<xarray.DataTree>
Group: /
β”‚   Dimensions:              (sweep: 9)
β”‚   Coordinates:
β”‚       latitude             float64 8B ...
β”‚       longitude            float64 8B ...
β”‚       altitude             float64 8B ...
β”‚   Dimensions without coordinates: sweep
β”‚   Data variables:
β”‚       volume_number        int32 4B ...
β”‚       platform_type        |S32 32B ...
β”‚       primary_axis         |S32 32B ...
β”‚       status_str           |S1 1B ...
β”‚       instrument_type      |S32 32B ...
β”‚       time_coverage_start  |S32 32B ...
β”‚       time_coverage_end    |S32 32B ...
β”‚       sweep_group_name     (sweep) <U7 252B 'sweep_0' 'sweep_1' ... 'sweep_8'
β”‚       sweep_fixed_angle    (sweep) float32 36B ...
β”‚   Attributes: (12/13)
β”‚       Conventions:         CF/Radial instrument_parameters radar_parameters rad...
β”‚       version:             1.2
β”‚       title:               TIMREX
β”‚       institution:         
β”‚       references:          
β”‚       source:              
β”‚       ...                  ...
β”‚       comment:             
β”‚       instrument_name:     SPOLRVP8
β”‚       site_name:           
β”‚       scan_name:           
β”‚       scan_id:             0
β”‚       platform_is_mobile:  false
β”œβ”€β”€ Group: /sweep_0
β”‚       Dimensions:                    (azimuth: 483, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:15:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 483B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_1
β”‚       Dimensions:                    (azimuth: 483, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:16:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 483B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_2
β”‚       Dimensions:                    (azimuth: 482, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:17:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 482B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_3
β”‚       Dimensions:                    (azimuth: 483, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:17:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 483B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_4
β”‚       Dimensions:                    (azimuth: 481, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:18:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 481B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_5
β”‚       Dimensions:                    (azimuth: 482, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:19:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 482B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_6
β”‚       Dimensions:                    (azimuth: 482, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:20:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 482B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_7
β”‚       Dimensions:                    (azimuth: 484, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:21:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 484B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
└── Group: /sweep_8
        Dimensions:                    (azimuth: 483, range: 996)
        Coordinates:
          * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
            elevation                  (azimuth) float32 2kB ...
            time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:21:...
          * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
        Data variables: (12/18)
            sweep_number               int32 4B ...
            sweep_mode                 <U20 80B 'azimuth_surveillance'
            prt_mode                   |S32 32B ...
            follow_mode                |S32 32B ...
            sweep_fixed_angle          float32 4B ...
            pulse_width                (azimuth) float32 2kB ...
            ...                         ...
            r_calib_index              (azimuth) int8 483B ...
            measured_transmit_power_h  (azimuth) float32 2kB ...
            measured_transmit_power_v  (azimuth) float32 2kB ...
            scan_rate                  (azimuth) float32 2kB ...
            DBZ                        (azimuth, range) float32 2MB ...
            VR                         (azimuth, range) float32 2MB ...

Inspect root group#

The sweep dimension contains the number of scans in this radar volume. Further the dataset consists of variables (location coordinates, time_coverage) and attributes (Conventions, metadata).

vol.root
<xarray.DataTree>
Group: /
β”‚   Dimensions:              (sweep: 9)
β”‚   Coordinates:
β”‚       latitude             float64 8B ...
β”‚       longitude            float64 8B ...
β”‚       altitude             float64 8B ...
β”‚   Dimensions without coordinates: sweep
β”‚   Data variables:
β”‚       volume_number        int32 4B ...
β”‚       platform_type        |S32 32B ...
β”‚       primary_axis         |S32 32B ...
β”‚       status_str           |S1 1B ...
β”‚       instrument_type      |S32 32B ...
β”‚       time_coverage_start  |S32 32B ...
β”‚       time_coverage_end    |S32 32B ...
β”‚       sweep_group_name     (sweep) <U7 252B 'sweep_0' 'sweep_1' ... 'sweep_8'
β”‚       sweep_fixed_angle    (sweep) float32 36B ...
β”‚   Attributes: (12/13)
β”‚       Conventions:         CF/Radial instrument_parameters radar_parameters rad...
β”‚       version:             1.2
β”‚       title:               TIMREX
β”‚       institution:         
β”‚       references:          
β”‚       source:              
β”‚       ...                  ...
β”‚       comment:             
β”‚       instrument_name:     SPOLRVP8
β”‚       site_name:           
β”‚       scan_name:           
β”‚       scan_id:             0
β”‚       platform_is_mobile:  false
β”œβ”€β”€ Group: /sweep_0
β”‚       Dimensions:                    (azimuth: 483, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:15:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 483B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_1
β”‚       Dimensions:                    (azimuth: 483, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:16:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 483B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_2
β”‚       Dimensions:                    (azimuth: 482, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:17:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 482B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_3
β”‚       Dimensions:                    (azimuth: 483, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:17:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 483B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_4
β”‚       Dimensions:                    (azimuth: 481, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:18:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 481B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_5
β”‚       Dimensions:                    (azimuth: 482, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:19:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 482B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_6
β”‚       Dimensions:                    (azimuth: 482, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:20:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 482B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
β”œβ”€β”€ Group: /sweep_7
β”‚       Dimensions:                    (azimuth: 484, range: 996)
β”‚       Coordinates:
β”‚         * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
β”‚           elevation                  (azimuth) float32 2kB ...
β”‚           time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:21:...
β”‚         * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
β”‚       Data variables: (12/18)
β”‚           sweep_number               int32 4B ...
β”‚           sweep_mode                 <U20 80B 'azimuth_surveillance'
β”‚           prt_mode                   |S32 32B ...
β”‚           follow_mode                |S32 32B ...
β”‚           sweep_fixed_angle          float32 4B ...
β”‚           pulse_width                (azimuth) float32 2kB ...
β”‚           ...                         ...
β”‚           r_calib_index              (azimuth) int8 484B ...
β”‚           measured_transmit_power_h  (azimuth) float32 2kB ...
β”‚           measured_transmit_power_v  (azimuth) float32 2kB ...
β”‚           scan_rate                  (azimuth) float32 2kB ...
β”‚           DBZ                        (azimuth, range) float32 2MB ...
β”‚           VR                         (azimuth, range) float32 2MB ...
└── Group: /sweep_8
        Dimensions:                    (azimuth: 483, range: 996)
        Coordinates:
          * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
            elevation                  (azimuth) float32 2kB ...
            time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:21:...
          * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
        Data variables: (12/18)
            sweep_number               int32 4B ...
            sweep_mode                 <U20 80B 'azimuth_surveillance'
            prt_mode                   |S32 32B ...
            follow_mode                |S32 32B ...
            sweep_fixed_angle          float32 4B ...
            pulse_width                (azimuth) float32 2kB ...
            ...                         ...
            r_calib_index              (azimuth) int8 483B ...
            measured_transmit_power_h  (azimuth) float32 2kB ...
            measured_transmit_power_v  (azimuth) float32 2kB ...
            scan_rate                  (azimuth) float32 2kB ...
            DBZ                        (azimuth, range) float32 2MB ...
            VR                         (azimuth, range) float32 2MB ...

Inspect sweep group(s)#

The sweep-groups can be accessed via their respective keys. The dimensions consist of range and time with added coordinates azimuth, elevation, range and time. There will be variables like radar moments (DBZH etc.) and sweep-dependent metadata (like fixed_angle, sweep_mode etc.).

display(vol["sweep_0"])
<xarray.DataTree 'sweep_0'>
Group: /sweep_0
    Dimensions:                    (sweep: 9, azimuth: 483, range: 996)
    Coordinates:
      * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
        elevation                  (azimuth) float32 2kB ...
        time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:15:...
      * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
    Dimensions without coordinates: sweep
    Data variables: (12/18)
        sweep_number               int32 4B ...
        sweep_mode                 <U20 80B 'azimuth_surveillance'
        prt_mode                   |S32 32B ...
        follow_mode                |S32 32B ...
        sweep_fixed_angle          float32 4B ...
        pulse_width                (azimuth) float32 2kB ...
        ...                         ...
        r_calib_index              (azimuth) int8 483B ...
        measured_transmit_power_h  (azimuth) float32 2kB ...
        measured_transmit_power_v  (azimuth) float32 2kB ...
        scan_rate                  (azimuth) float32 2kB ...
        DBZ                        (azimuth, range) float32 2MB ...
        VR                         (azimuth, range) float32 2MB ...

Georeferencing#

sweep_mode is assigned coordinate, as we need it available on the DataArray.

swp = vol["sweep_0"].to_dataset(inherit="all_coords")
swp = swp.assign_coords(sweep_mode=swp.sweep_mode)
swp = swp.wrl.georef.georeference()
display(swp)
<xarray.Dataset> Size: 23MB
Dimensions:                    (azimuth: 483, range: 996)
Coordinates: (12/15)
  * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
    elevation                  (azimuth) float32 2kB 0.5164 0.5219 ... 0.5219
    time                       (azimuth) datetime64[ns] 4kB 2008-06-04T00:15:...
  * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
    x                          (azimuth, range) float64 4MB -6.556e-06 ... -1...
    y                          (azimuth, range) float64 4MB 150.0 ... 1.493e+05
    ...                         ...
    bins                       (azimuth, range) float32 2MB 150.0 ... 1.494e+05
    sweep_mode                 <U20 80B 'azimuth_surveillance'
    latitude                   float64 8B 22.53
    longitude                  float64 8B 120.4
    altitude                   float64 8B 45.0
    crs_wkt                    int64 8B 0
Data variables: (12/17)
    sweep_number               int32 4B ...
    prt_mode                   |S32 32B ...
    follow_mode                |S32 32B ...
    sweep_fixed_angle          float32 4B ...
    pulse_width                (azimuth) float32 2kB ...
    prt                        (azimuth) float32 2kB ...
    ...                         ...
    r_calib_index              (azimuth) int8 483B ...
    measured_transmit_power_h  (azimuth) float32 2kB ...
    measured_transmit_power_v  (azimuth) float32 2kB ...
    scan_rate                  (azimuth) float32 2kB ...
    DBZ                        (azimuth, range) float32 2MB ...
    VR                         (azimuth, range) float32 2MB ...

Inspect radar moments#

The DataArrays can be accessed by key or by attribute. Each DataArray has dimensions and coordinates of its parent dataset. There are attributes connected which are defined by Cf/Radial standard.

display(swp.DBZ)
<xarray.DataArray 'DBZ' (azimuth: 483, range: 996)> Size: 2MB
[481068 values with dtype=float32]
Coordinates: (12/15)
  * azimuth     (azimuth) float32 2kB 0.0 0.75 1.5 2.25 ... 357.8 358.5 359.2
    elevation   (azimuth) float32 2kB 0.5164 0.5219 0.5164 ... 0.5219 0.5219
    time        (azimuth) datetime64[ns] 4kB 2008-06-04T00:15:34 ... 2008-06-...
  * range       (range) float32 4kB 150.0 300.0 450.0 ... 1.492e+05 1.494e+05
    x           (azimuth, range) float64 4MB -6.556e-06 ... -1.955e+03
    y           (azimuth, range) float64 4MB 150.0 300.0 ... 1.492e+05 1.493e+05
    ...          ...
    bins        (azimuth, range) float32 2MB 150.0 300.0 ... 1.492e+05 1.494e+05
    sweep_mode  <U20 80B 'azimuth_surveillance'
    latitude    float64 8B 22.53
    longitude   float64 8B 120.4
    altitude    float64 8B 45.0
    crs_wkt     int64 8B 0
Attributes:
    long_name:             Computed Horizontal Co-polar Reflectivit
    standard_name:         equivalent_reflectivity_factor
    units:                 dBZ
    threshold_field_name:  
    threshold_value:       -9999.0
    sampling_ratio:        1.0
    grid_mapping:          grid_mapping

Create simple plot#

Using xarray features a simple plot can be created like this. Note the sortby('time') method, which sorts the radials by time.

For more details on plotting radar data see under Visualization.

swp.DBZ.sortby("time").plot(x="range", y="time", add_labels=False)
<matplotlib.collections.QuadMesh at 0x72db00af4d70>
../../../_images/5f8c34267d5ebf17d89c27faf1878cbcac64d0e52d7045e3fa5ae0f63a6948c6.png
fig = plt.figure(figsize=(5, 5))
pm = swp.DBZ.wrl.vis.plot(crs={"latmin": 3e3}, fig=fig)
../../../_images/9450f18123b6e36cde04476fc80dfb5aea92f5ccefa3ab070bbb1860da27c319.png

Use xr.open_dataset to retrieve explicit group#

Warning

Since \(\omega radlib\) version 2.0 all xarray backend related functionality is imported from xradar-package.

swp_b = xr.open_dataset(
    f, engine="cfradial1", group="sweep_1", backend_kwargs=dict(reindex_angle=False)
)
display(swp_b)
<xarray.Dataset> Size: 4MB
Dimensions:                    (azimuth: 483, range: 996)
Coordinates:
  * azimuth                    (azimuth) float32 2kB 0.0 0.75 ... 358.5 359.2
    elevation                  (azimuth) float32 2kB ...
    time                       (azimuth) datetime64[ns] 4kB ...
  * range                      (range) float32 4kB 150.0 300.0 ... 1.494e+05
    latitude                   float64 8B ...
    longitude                  float64 8B ...
    altitude                   float64 8B ...
Data variables: (12/18)
    sweep_number               int32 4B ...
    sweep_mode                 <U20 80B ...
    prt_mode                   |S32 32B ...
    follow_mode                |S32 32B ...
    sweep_fixed_angle          float32 4B ...
    pulse_width                (azimuth) float32 2kB ...
    ...                         ...
    r_calib_index              (azimuth) int8 483B ...
    measured_transmit_power_h  (azimuth) float32 2kB ...
    measured_transmit_power_v  (azimuth) float32 2kB ...
    scan_rate                  (azimuth) float32 2kB ...
    DBZ                        (azimuth, range) float32 2MB ...
    VR                         (azimuth, range) float32 2MB ...