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_mappingCreate 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>
fig = plt.figure(figsize=(5, 5))
pm = swp.DBZ.wrl.vis.plot(crs={"latmin": 3e3}, fig=fig)
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 ...