NDPointIndex

This uses one of the new indexes xarray is providing. See xarray-indexes for more in-depth details on that.

Needs latest xarray ‘2025.7.1’

import cmweather  # noqa
import numpy as np
import xarray as xr
from open_radar_data import DATASETS

xr.__version__

'2026.4.0'

Step-by-step guide

Load radar data

fname = DATASETS.fetch("DWD-Vol-2_99999_20180601054047_00.h5")
ds = xr.open_dataset(fname, engine="gamic", group="sweep_9")
display(ds)

<xarray.Dataset> Size: 17MB
Dimensions:            (azimuth: 360, range: 1000)
Coordinates:
  * azimuth            (azimuth) float64 3kB 0.5109 1.519 2.519 ... 358.5 359.5
    elevation          (azimuth) float64 3kB ...
    time               (azimuth) datetime64[ns] 3kB ...
  * range              (range) float32 4kB 75.0 225.0 ... 1.498e+05 1.499e+05
    longitude          float64 8B ...
    latitude           float64 8B ...
    altitude           float64 8B ...
Data variables: (12/18)
    DBZH               (azimuth, range) float32 1MB ...
    DBZV               (azimuth, range) float32 1MB ...
    KDP                (azimuth, range) float32 1MB ...
    RHOHV              (azimuth, range) float32 1MB ...
    DBTH               (azimuth, range) float32 1MB ...
    DBTV               (azimuth, range) float32 1MB ...
    ...                 ...
    sweep_mode         <U20 80B ...
    sweep_number       int64 8B ...
    prt_mode           <U7 28B ...
    follow_mode        <U7 28B ...
    sweep_fixed_angle  float64 8B ...
    nyquist_velocity   float64 8B ...
Attributes:
    source:         gamic
    ant_gain_h:     43
    ant_gain_v:     43
    noise_power_h:  -3.40133
    noise_power_v:  -3.13518
    rx_loss_h:      3
    rx_loss_v:      3

xarray.Dataset

• Dimensions:
  • azimuth: 360
  • range: 1000
• Coordinates: (7)
  • azimuth
    (azimuth)
    float64
    0.5109 1.519 2.519 ... 358.5 359.5

    standard_name :

    ray_azimuth_angle

    long_name :

    azimuth_angle_from_true_north

    units :

    degrees

    axis :

    radial_azimuth_coordinate

    array([  0.510864,   1.51886 ,   2.518616, ..., 357.511597, 358.511353,
           359.505615], shape=(360,))

  • elevation
    (azimuth)
    float64
    ...

    standard_name :

    ray_elevation_angle

    long_name :

    elevation_angle_from_horizontal_plane

    units :

    degrees

    axis :

    radial_elevation_coordinate

    [360 values with dtype=float64]

  • time
    (azimuth)
    datetime64[ns]
    ...

    standard_name :

    time

    [360 values with dtype=datetime64[ns]]

  • range
    (range)
    float32
    75.0 225.0 ... 1.498e+05 1.499e+05

    units :

    meters

    standard_name :

    projection_range_coordinate

    long_name :

    range_to_measurement_volume

    axis :

    radial_range_coordinate

    meters_between_gates :

    150.0

    spacing_is_constant :

    true

    meters_to_center_of_first_gate :

    75.0

    array([7.50000e+01, 2.25000e+02, 3.75000e+02, ..., 1.49625e+05, 1.49775e+05,
           1.49925e+05], shape=(1000,), dtype=float32)

  • longitude
    ()
    float64
    ...

    long_name :

    longitude

    units :

    degrees_east

    standard_name :

    longitude

    [1 values with dtype=float64]

  • latitude
    ()
    float64
    ...

    long_name :

    latitude

    units :

    degrees_north

    positive :

    up

    standard_name :

    latitude

    [1 values with dtype=float64]

  • altitude
    ()
    float64
    ...

    long_name :

    altitude

    units :

    meters

    standard_name :

    altitude

    [1 values with dtype=float64]

• Data variables: (18)
  • DBZH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_h

    long_name :

    Equivalent reflectivity factor H

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • DBZV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_v

    long_name :

    Equivalent reflectivity factor V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • KDP
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    deg/km

    units :

    degrees per kilometer

    standard_name :

    radar_specific_differential_phase_hv

    long_name :

    Specific differential phase HV

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • RHOHV
    (azimuth, range)
    float32
    ...

    format :

    UV16

    is_dft :

    0

    units :

    unitless

    standard_name :

    radar_correlation_coefficient_hv

    long_name :

    Correlation coefficient HV

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • DBTH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_h

    long_name :

    Total power H (uncorrected reflectivity)

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • DBTV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_v

    long_name :

    Total power V (uncorrected reflectivity)

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • ZDR
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dB

    units :

    dB

    standard_name :

    radar_differential_reflectivity_hv

    long_name :

    Log differential reflectivity H/V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • VRADH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per seconds

    standard_name :

    radial_velocity_of_scatterers_away_from_instrument_h

    long_name :

    Radial velocity of scatterers away from instrument H

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • VRADV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per second

    standard_name :

    radial_velocity_of_scatterers_away_from_instrument_v

    long_name :

    Radial velocity of scatterers away from instrument V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • WRADH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per seconds

    standard_name :

    radar_doppler_spectrum_width_h

    long_name :

    Doppler spectrum width H

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • WRADV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per second

    standard_name :

    radar_doppler_spectrum_width_v

    long_name :

    Doppler spectrum width V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • PHIDP
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    deg

    units :

    degrees

    standard_name :

    radar_differential_phase_hv

    long_name :

    Differential phase HV

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • sweep_mode
    ()
    <U20
    ...

    [1 values with dtype=<U20]

  • 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
    ()
    float64
    ...

    standard_name :

    nyquist_velocity

    units :

    m s-1

    [1 values with dtype=float64]

• Attributes: (7)

  source :

  gamic

  ant_gain_h :

  43

  ant_gain_v :

  43

  noise_power_h :

  -3.40133

  noise_power_v :

  -3.13518

  rx_loss_h :

  3

  rx_loss_v :

  3

Georeference

Add x,y,z - 2D-coordinates.

ds = ds.xradar.georeference()
display(ds)

<xarray.Dataset> Size: 26MB
Dimensions:            (azimuth: 360, range: 1000)
Coordinates:
  * azimuth            (azimuth) float64 3kB 0.5109 1.519 2.519 ... 358.5 359.5
    elevation          (azimuth) float64 3kB 0.5988 0.5988 ... 0.5988 0.5988
    time               (azimuth) datetime64[ns] 3kB ...
  * range              (range) float32 4kB 75.0 225.0 ... 1.498e+05 1.499e+05
    x                  (azimuth, range) float64 3MB 0.6686 2.006 ... -1.293e+03
    y                  (azimuth, range) float64 3MB 74.99 225.0 ... 1.499e+05
    z                  (azimuth, range) float64 3MB 310.8 312.4 ... 3.2e+03
    longitude          float64 8B 6.457
    latitude           float64 8B 50.93
    altitude           float64 8B 310.0
    crs_wkt            int64 8B 0
Data variables: (12/18)
    DBZH               (azimuth, range) float32 1MB ...
    DBZV               (azimuth, range) float32 1MB ...
    KDP                (azimuth, range) float32 1MB ...
    RHOHV              (azimuth, range) float32 1MB ...
    DBTH               (azimuth, range) float32 1MB ...
    DBTV               (azimuth, range) float32 1MB ...
    ...                 ...
    sweep_mode         <U20 80B ...
    sweep_number       int64 8B ...
    prt_mode           <U7 28B ...
    follow_mode        <U7 28B ...
    sweep_fixed_angle  float64 8B ...
    nyquist_velocity   float64 8B ...
Attributes:
    source:         gamic
    ant_gain_h:     43
    ant_gain_v:     43
    noise_power_h:  -3.40133
    noise_power_v:  -3.13518
    rx_loss_h:      3
    rx_loss_v:      3

xarray.Dataset

• Dimensions:
  • azimuth: 360
  • range: 1000
• Coordinates: (11)
  • azimuth
    (azimuth)
    float64
    0.5109 1.519 2.519 ... 358.5 359.5

    standard_name :

    ray_azimuth_angle

    long_name :

    azimuth_angle_from_true_north

    units :

    degrees

    axis :

    radial_azimuth_coordinate

    array([  0.510864,   1.51886 ,   2.518616, ..., 357.511597, 358.511353,
           359.505615], shape=(360,))

  • elevation
    (azimuth)
    float64
    0.5988 0.5988 ... 0.5988 0.5988

    standard_name :

    ray_elevation_angle

    long_name :

    elevation_angle_from_horizontal_plane

    units :

    degrees

    axis :

    radial_elevation_coordinate

    array([0.598755, 0.598755, 0.598755, ..., 0.598755, 0.598755, 0.598755],
          shape=(360,))

  • time
    (azimuth)
    datetime64[ns]
    ...

    standard_name :

    time

    [360 values with dtype=datetime64[ns]]

  • range
    (range)
    float32
    75.0 225.0 ... 1.498e+05 1.499e+05

    units :

    meters

    standard_name :

    projection_range_coordinate

    long_name :

    range_to_measurement_volume

    axis :

    radial_range_coordinate

    meters_between_gates :

    150.0

    spacing_is_constant :

    true

    meters_to_center_of_first_gate :

    75.0

    array([7.50000e+01, 2.25000e+02, 3.75000e+02, ..., 1.49625e+05, 1.49775e+05,
           1.49925e+05], shape=(1000,), dtype=float32)

  • x
    (azimuth, range)
    float64
    0.6686 2.006 ... -1.293e+03

    standard_name :

    easting

    long_name :

    Easting

    units :

    meters

    axis :

    X

    array([[ 6.68649890e-01,  2.00594930e+00,  3.34324822e+00, ...,
             1.33357288e+03,  1.33490927e+03,  1.33624566e+03],
           [ 1.98776872e+00,  5.96330506e+00,  9.93883993e+00, ...,
             3.96445807e+03,  3.96843090e+03,  3.97240374e+03],
           [ 3.29549797e+00,  9.88649209e+00,  1.64774838e+01, ...,
             6.57262758e+03,  6.57921410e+03,  6.58580061e+03],
           ...,
           [-3.25599138e+00, -9.76797234e+00, -1.62799509e+01, ...,
            -6.49383460e+03, -6.50034216e+03, -6.50684971e+03],
           [-1.94823801e+00, -5.84471294e+00, -9.74118643e+00, ...,
            -3.88561697e+03, -3.88951080e+03, -3.89340463e+03],
           [-6.47081083e-01, -1.94124289e+00, -3.23540422e+00, ...,
            -1.29055548e+03, -1.29184876e+03, -1.29314205e+03]],
          shape=(360, 1000))

  • y
    (azimuth, range)
    float64
    74.99 225.0 ... 1.497e+05 1.499e+05

    standard_name :

    northing

    long_name :

    Northing

    units :

    meters

    axis :

    Y

    array([[7.49901777e+01, 2.24970491e+02, 3.74950750e+02, ...,
            1.49562377e+05, 1.49712255e+05, 1.49862134e+05],
           [7.49668101e+01, 2.24900389e+02, 3.74833912e+02, ...,
            1.49515772e+05, 1.49665604e+05, 1.49815435e+05],
           [7.49207150e+01, 2.24762103e+02, 3.74603436e+02, ...,
            1.49423839e+05, 1.49573578e+05, 1.49723318e+05],
           ...,
           [7.49224423e+01, 2.24767285e+02, 3.74612073e+02, ...,
            1.49427284e+05, 1.49577027e+05, 1.49726770e+05],
           [7.49678478e+01, 2.24903502e+02, 3.74839100e+02, ...,
            1.49517842e+05, 1.49667676e+05, 1.49817509e+05],
           [7.49903669e+01, 2.24971059e+02, 3.74951696e+02, ...,
            1.49562754e+05, 1.49712633e+05, 1.49862512e+05]],
          shape=(360, 1000))

  • z
    (azimuth, range)
    float64
    310.8 312.4 ... 3.196e+03 3.2e+03

    meters_between_gates :

    150.0

    spacing_is_constant :

    true

    meters_to_center_of_first_gate :

    75.0

    units :

    meters

    standard_name :

    projection_range_coordinate

    long_name :

    range_to_measurement_volume

    axis :

    radial_range_coordinate

    array([[ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201],
           [ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201],
           [ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201],
           ...,
           [ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201],
           [ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201],
           [ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201]], shape=(360, 1000))

  • longitude
    ()
    float64
    6.457

    long_name :

    longitude

    units :

    degrees_east

    standard_name :

    longitude

    array(6.456949)

  • latitude
    ()
    float64
    50.93

    long_name :

    latitude

    units :

    degrees_north

    positive :

    up

    standard_name :

    latitude

    array(50.928727)

  • altitude
    ()
    float64
    310.0

    long_name :

    altitude

    units :

    meters

    standard_name :

    altitude

    array(310.)

  • crs_wkt
    ()
    int64
    0

    crs_wkt :

    PROJCRS["unknown",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",50.9287272,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",6.4569489,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 :

    unknown

    grid_mapping_name :

    azimuthal_equidistant

    latitude_of_projection_origin :

    50.9287272

    longitude_of_projection_origin :

    6.4569489

    false_easting :

    0.0

    false_northing :

    0.0

    array(0)

• Data variables: (18)
  • DBZH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_h

    long_name :

    Equivalent reflectivity factor H

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • DBZV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_v

    long_name :

    Equivalent reflectivity factor V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • KDP
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    deg/km

    units :

    degrees per kilometer

    standard_name :

    radar_specific_differential_phase_hv

    long_name :

    Specific differential phase HV

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • RHOHV
    (azimuth, range)
    float32
    ...

    format :

    UV16

    is_dft :

    0

    units :

    unitless

    standard_name :

    radar_correlation_coefficient_hv

    long_name :

    Correlation coefficient HV

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • DBTH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_h

    long_name :

    Total power H (uncorrected reflectivity)

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • DBTV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_v

    long_name :

    Total power V (uncorrected reflectivity)

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • ZDR
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dB

    units :

    dB

    standard_name :

    radar_differential_reflectivity_hv

    long_name :

    Log differential reflectivity H/V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • VRADH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per seconds

    standard_name :

    radial_velocity_of_scatterers_away_from_instrument_h

    long_name :

    Radial velocity of scatterers away from instrument H

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • VRADV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per second

    standard_name :

    radial_velocity_of_scatterers_away_from_instrument_v

    long_name :

    Radial velocity of scatterers away from instrument V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • WRADH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per seconds

    standard_name :

    radar_doppler_spectrum_width_h

    long_name :

    Doppler spectrum width H

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • WRADV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per second

    standard_name :

    radar_doppler_spectrum_width_v

    long_name :

    Doppler spectrum width V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • PHIDP
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    deg

    units :

    degrees

    standard_name :

    radar_differential_phase_hv

    long_name :

    Differential phase HV

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • sweep_mode
    ()
    <U20
    ...

    [1 values with dtype=<U20]

  • 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
    ()
    float64
    ...

    standard_name :

    nyquist_velocity

    units :

    m s-1

    [1 values with dtype=float64]

• Attributes: (7)

  source :

  gamic

  ant_gain_h :

  43

  ant_gain_v :

  43

  noise_power_h :

  -3.40133

  noise_power_v :

  -3.13518

  rx_loss_h :

  3

  rx_loss_v :

  3

Add NDPointIndex

This uses scipy.KDTree under the hood. See also Indexes-Section in the html-repr below.

ds = ds.set_xindex(("x", "y"), xr.indexes.NDPointIndex)
display(ds)

<xarray.Dataset> Size: 26MB
Dimensions:            (azimuth: 360, range: 1000)
Coordinates:
  * azimuth            (azimuth) float64 3kB 0.5109 1.519 2.519 ... 358.5 359.5
    elevation          (azimuth) float64 3kB 0.5988 0.5988 ... 0.5988 0.5988
    time               (azimuth) datetime64[ns] 3kB ...
  * range              (range) float32 4kB 75.0 225.0 ... 1.498e+05 1.499e+05
    z                  (azimuth, range) float64 3MB 310.8 312.4 ... 3.2e+03
  * x                  (azimuth, range) float64 3MB 0.6686 2.006 ... -1.293e+03
  * y                  (azimuth, range) float64 3MB 74.99 225.0 ... 1.499e+05
    longitude          float64 8B 6.457
    latitude           float64 8B 50.93
    altitude           float64 8B 310.0
    crs_wkt            int64 8B 0
Data variables: (12/18)
    DBZH               (azimuth, range) float32 1MB ...
    DBZV               (azimuth, range) float32 1MB ...
    KDP                (azimuth, range) float32 1MB ...
    RHOHV              (azimuth, range) float32 1MB ...
    DBTH               (azimuth, range) float32 1MB ...
    DBTV               (azimuth, range) float32 1MB ...
    ...                 ...
    sweep_mode         <U20 80B ...
    sweep_number       int64 8B ...
    prt_mode           <U7 28B ...
    follow_mode        <U7 28B ...
    sweep_fixed_angle  float64 8B ...
    nyquist_velocity   float64 8B ...
Indexes:
  ┌ x        NDPointIndex (ScipyKDTreeAdapter)
  └ y
Attributes:
    source:         gamic
    ant_gain_h:     43
    ant_gain_v:     43
    noise_power_h:  -3.40133
    noise_power_v:  -3.13518
    rx_loss_h:      3
    rx_loss_v:      3

xarray.Dataset

• Dimensions:
  • azimuth: 360
  • range: 1000
• Coordinates: (11)
  • azimuth
    (azimuth)
    float64
    0.5109 1.519 2.519 ... 358.5 359.5

    standard_name :

    ray_azimuth_angle

    long_name :

    azimuth_angle_from_true_north

    units :

    degrees

    axis :

    radial_azimuth_coordinate

    array([  0.510864,   1.51886 ,   2.518616, ..., 357.511597, 358.511353,
           359.505615], shape=(360,))

  • elevation
    (azimuth)
    float64
    0.5988 0.5988 ... 0.5988 0.5988

    standard_name :

    ray_elevation_angle

    long_name :

    elevation_angle_from_horizontal_plane

    units :

    degrees

    axis :

    radial_elevation_coordinate

    array([0.598755, 0.598755, 0.598755, ..., 0.598755, 0.598755, 0.598755],
          shape=(360,))

  • time
    (azimuth)
    datetime64[ns]
    ...

    standard_name :

    time

    [360 values with dtype=datetime64[ns]]

  • range
    (range)
    float32
    75.0 225.0 ... 1.498e+05 1.499e+05

    units :

    meters

    standard_name :

    projection_range_coordinate

    long_name :

    range_to_measurement_volume

    axis :

    radial_range_coordinate

    meters_between_gates :

    150.0

    spacing_is_constant :

    true

    meters_to_center_of_first_gate :

    75.0

    array([7.50000e+01, 2.25000e+02, 3.75000e+02, ..., 1.49625e+05, 1.49775e+05,
           1.49925e+05], shape=(1000,), dtype=float32)

  • z
    (azimuth, range)
    float64
    310.8 312.4 ... 3.196e+03 3.2e+03

    meters_between_gates :

    150.0

    spacing_is_constant :

    true

    meters_to_center_of_first_gate :

    75.0

    units :

    meters

    standard_name :

    projection_range_coordinate

    long_name :

    range_to_measurement_volume

    axis :

    radial_range_coordinate

    array([[ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201],
           [ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201],
           [ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201],
           ...,
           [ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201],
           [ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201],
           [ 310.78408538,  312.35424418,  313.92705368, ..., 3191.98063179,
            3196.19244125, 3200.40690201]], shape=(360, 1000))

  • x
    (azimuth, range)
    float64
    0.6686 2.006 ... -1.293e+03

    standard_name :

    easting

    long_name :

    Easting

    units :

    meters

    axis :

    X

    array([[ 6.68649890e-01,  2.00594930e+00,  3.34324822e+00, ...,
             1.33357288e+03,  1.33490927e+03,  1.33624566e+03],
           [ 1.98776872e+00,  5.96330506e+00,  9.93883993e+00, ...,
             3.96445807e+03,  3.96843090e+03,  3.97240374e+03],
           [ 3.29549797e+00,  9.88649209e+00,  1.64774838e+01, ...,
             6.57262758e+03,  6.57921410e+03,  6.58580061e+03],
           ...,
           [-3.25599138e+00, -9.76797234e+00, -1.62799509e+01, ...,
            -6.49383460e+03, -6.50034216e+03, -6.50684971e+03],
           [-1.94823801e+00, -5.84471294e+00, -9.74118643e+00, ...,
            -3.88561697e+03, -3.88951080e+03, -3.89340463e+03],
           [-6.47081083e-01, -1.94124289e+00, -3.23540422e+00, ...,
            -1.29055548e+03, -1.29184876e+03, -1.29314205e+03]],
          shape=(360, 1000))

  • y
    (azimuth, range)
    float64
    74.99 225.0 ... 1.497e+05 1.499e+05

    standard_name :

    northing

    long_name :

    Northing

    units :

    meters

    axis :

    Y

    array([[7.49901777e+01, 2.24970491e+02, 3.74950750e+02, ...,
            1.49562377e+05, 1.49712255e+05, 1.49862134e+05],
           [7.49668101e+01, 2.24900389e+02, 3.74833912e+02, ...,
            1.49515772e+05, 1.49665604e+05, 1.49815435e+05],
           [7.49207150e+01, 2.24762103e+02, 3.74603436e+02, ...,
            1.49423839e+05, 1.49573578e+05, 1.49723318e+05],
           ...,
           [7.49224423e+01, 2.24767285e+02, 3.74612073e+02, ...,
            1.49427284e+05, 1.49577027e+05, 1.49726770e+05],
           [7.49678478e+01, 2.24903502e+02, 3.74839100e+02, ...,
            1.49517842e+05, 1.49667676e+05, 1.49817509e+05],
           [7.49903669e+01, 2.24971059e+02, 3.74951696e+02, ...,
            1.49562754e+05, 1.49712633e+05, 1.49862512e+05]],
          shape=(360, 1000))

  • longitude
    ()
    float64
    6.457

    long_name :

    longitude

    units :

    degrees_east

    standard_name :

    longitude

    array(6.456949)

  • latitude
    ()
    float64
    50.93

    long_name :

    latitude

    units :

    degrees_north

    positive :

    up

    standard_name :

    latitude

    array(50.928727)

  • altitude
    ()
    float64
    310.0

    long_name :

    altitude

    units :

    meters

    standard_name :

    altitude

    array(310.)

  • crs_wkt
    ()
    int64
    0

    crs_wkt :

    PROJCRS["unknown",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",50.9287272,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",6.4569489,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 :

    unknown

    grid_mapping_name :

    azimuthal_equidistant

    latitude_of_projection_origin :

    50.9287272

    longitude_of_projection_origin :

    6.4569489

    false_easting :

    0.0

    false_northing :

    0.0

    array(0)

• Data variables: (18)
  • DBZH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_h

    long_name :

    Equivalent reflectivity factor H

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • DBZV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_v

    long_name :

    Equivalent reflectivity factor V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • KDP
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    deg/km

    units :

    degrees per kilometer

    standard_name :

    radar_specific_differential_phase_hv

    long_name :

    Specific differential phase HV

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • RHOHV
    (azimuth, range)
    float32
    ...

    format :

    UV16

    is_dft :

    0

    units :

    unitless

    standard_name :

    radar_correlation_coefficient_hv

    long_name :

    Correlation coefficient HV

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • DBTH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_h

    long_name :

    Total power H (uncorrected reflectivity)

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • DBTV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_v

    long_name :

    Total power V (uncorrected reflectivity)

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • ZDR
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dB

    units :

    dB

    standard_name :

    radar_differential_reflectivity_hv

    long_name :

    Log differential reflectivity H/V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • VRADH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per seconds

    standard_name :

    radial_velocity_of_scatterers_away_from_instrument_h

    long_name :

    Radial velocity of scatterers away from instrument H

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • VRADV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per second

    standard_name :

    radial_velocity_of_scatterers_away_from_instrument_v

    long_name :

    Radial velocity of scatterers away from instrument V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • WRADH
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per seconds

    standard_name :

    radar_doppler_spectrum_width_h

    long_name :

    Doppler spectrum width H

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • WRADV
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per second

    standard_name :

    radar_doppler_spectrum_width_v

    long_name :

    Doppler spectrum width V

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • PHIDP
    (azimuth, range)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    deg

    units :

    degrees

    standard_name :

    radar_differential_phase_hv

    long_name :

    Differential phase HV

    _Undetect :

    0.0

    [360000 values with dtype=float32]

  • sweep_mode
    ()
    <U20
    ...

    [1 values with dtype=<U20]

  • 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
    ()
    float64
    ...

    standard_name :

    nyquist_velocity

    units :

    m s-1

    [1 values with dtype=float64]

• Indexes: (1)
  • x
    y
    NDPointIndex (ScipyKDTreeAdapter)

    <xarray.indexes.nd_point_index.NDPointIndex object at 0x7959c575f230>

• Attributes: (7)

  source :

  gamic

  ant_gain_h :

  43

  ant_gain_v :

  43

  noise_power_h :

  -3.40133

  noise_power_v :

  -3.13518

  rx_loss_h :

  3

  rx_loss_v :

  3

Plot

This works as usual with the new NDPointIndex.

ds.DBZH.plot(
    x="x", y="y", xlim=(-10e3, 10e3), ylim=(-10e3, 10e3), cmap="HomeyerRainbow", vmin=0
)

<matplotlib.collections.QuadMesh at 0x7959c49b8980>

Nearest neighbour interpolation with NDPointIndex

Create 1D DataArrays for x and y selection

y = xr.DataArray(np.arange(-100e3, 100e3, 500), dims="y", name="y", attrs=ds.y.attrs)
x = xr.DataArray(np.arange(-100e3, 100e3, 500), dims="x", name="x", attrs=ds.x.attrs)

Select with above 1D DataArrays

actual = ds.sel(y=y, x=x, method="nearest")

Assign the 1D DataArrays

actual = actual.assign(x=x, y=y)

display(actual)

<xarray.Dataset> Size: 13MB
Dimensions:            (y: 400, x: 400)
Coordinates:
  * y                  (y) float64 3kB -1e+05 -9.95e+04 ... 9.9e+04 9.95e+04
  * x                  (x) float64 3kB -1e+05 -9.95e+04 ... 9.9e+04 9.95e+04
    elevation          (y, x) float64 1MB 0.5988 0.5988 0.5988 ... 0.5988 0.5988
    time               (y, x) datetime64[ns] 1MB ...
    range              (y, x) float32 640kB 1.415e+05 1.411e+05 ... 1.408e+05
    azimuth            (y, x) float64 1MB 224.5 224.5 224.5 ... 44.51 44.51
    z                  (y, x) float64 1MB 2.968e+03 2.956e+03 ... 2.948e+03
    longitude          float64 8B 6.457
    latitude           float64 8B 50.93
    altitude           float64 8B 310.0
    crs_wkt            int64 8B 0
Data variables: (12/18)
    DBZH               (y, x) float32 640kB ...
    DBZV               (y, x) float32 640kB ...
    KDP                (y, x) float32 640kB ...
    RHOHV              (y, x) float32 640kB ...
    DBTH               (y, x) float32 640kB ...
    DBTV               (y, x) float32 640kB ...
    ...                 ...
    sweep_mode         <U20 80B ...
    sweep_number       int64 8B ...
    prt_mode           <U7 28B ...
    follow_mode        <U7 28B ...
    sweep_fixed_angle  float64 8B ...
    nyquist_velocity   float64 8B ...
Attributes:
    source:         gamic
    ant_gain_h:     43
    ant_gain_v:     43
    noise_power_h:  -3.40133
    noise_power_v:  -3.13518
    rx_loss_h:      3
    rx_loss_v:      3

xarray.Dataset

• Dimensions:
  • y: 400
  • x: 400
• Coordinates: (11)
  • y
    (y)
    float64
    -1e+05 -9.95e+04 ... 9.95e+04

    standard_name :

    northing

    long_name :

    Northing

    units :

    meters

    axis :

    Y

    array([-100000.,  -99500.,  -99000., ...,   98500.,   99000.,   99500.],
          shape=(400,))

  • x
    (x)
    float64
    -1e+05 -9.95e+04 ... 9.95e+04

    standard_name :

    easting

    long_name :

    Easting

    units :

    meters

    axis :

    X

    array([-100000.,  -99500.,  -99000., ...,   98500.,   99000.,   99500.],
          shape=(400,))

  • elevation
    (y, x)
    float64
    0.5988 0.5988 ... 0.5988 0.5988

    standard_name :

    ray_elevation_angle

    long_name :

    elevation_angle_from_horizontal_plane

    units :

    degrees

    axis :

    radial_elevation_coordinate

    array([[0.598755, 0.598755, 0.598755, ..., 0.598755, 0.598755, 0.598755],
           [0.598755, 0.598755, 0.598755, ..., 0.598755, 0.598755, 0.598755],
           [0.598755, 0.598755, 0.598755, ..., 0.598755, 0.598755, 0.598755],
           ...,
           [0.598755, 0.598755, 0.598755, ..., 0.598755, 0.598755, 0.598755],
           [0.598755, 0.598755, 0.598755, ..., 0.598755, 0.598755, 0.598755],
           [0.598755, 0.598755, 0.598755, ..., 0.598755, 0.598755, 0.598755]],
          shape=(400, 400))

  • time
    (y, x)
    datetime64[ns]
    ...

    standard_name :

    time

    [160000 values with dtype=datetime64[ns]]

  • range
    (y, x)
    float32
    1.415e+05 1.411e+05 ... 1.408e+05

    units :

    meters

    standard_name :

    projection_range_coordinate

    long_name :

    range_to_measurement_volume

    axis :

    radial_range_coordinate

    meters_between_gates :

    150.0

    spacing_is_constant :

    true

    meters_to_center_of_first_gate :

    75.0

    array([[141525., 141075., 140775., ..., 140475., 140775., 141075.],
           [141075., 140775., 140475., ..., 140025., 140475., 140775.],
           [140775., 140475., 140025., ..., 139725., 140025., 140475.],
           ...,
           [140475., 140025., 139725., ..., 139275., 139725., 140025.],
           [140775., 140475., 140025., ..., 139725., 140025., 140475.],
           [141075., 140775., 140475., ..., 140025., 140475., 140775.]],
          shape=(400, 400), dtype=float32)

  • azimuth
    (y, x)
    float64
    224.5 224.5 224.5 ... 44.51 44.51

    standard_name :

    ray_azimuth_angle

    long_name :

    azimuth_angle_from_true_north

    units :

    degrees

    axis :

    radial_azimuth_coordinate

    array([[224.51385498, 224.51385498, 224.51385498, ..., 135.51361084,
            135.51361084, 135.51361084],
           [225.519104  , 224.51385498, 224.51385498, ..., 135.51361084,
            135.51361084, 134.51385498],
           [225.519104  , 225.519104  , 224.51385498, ..., 135.51361084,
            134.51385498, 134.51385498],
           ...,
           [314.50836182, 314.50836182, 314.50836182, ...,  44.50836182,
             45.51635742,  45.51635742],
           [314.50836182, 314.50836182, 314.50836182, ...,  44.50836182,
             44.50836182,  45.51635742],
           [314.50836182, 314.50836182, 315.51086426, ...,  44.50836182,
             44.50836182,  44.50836182]], shape=(400, 400))

  • z
    (y, x)
    float64
    2.968e+03 2.956e+03 ... 2.948e+03

    meters_between_gates :

    150.0

    spacing_is_constant :

    true

    meters_to_center_of_first_gate :

    75.0

    units :

    meters

    standard_name :

    projection_range_coordinate

    long_name :

    range_to_measurement_volume

    axis :

    radial_range_coordinate

    array([[2968.47720336, 2956.28673058, 2948.17308678, ..., 2940.06992852,
            2948.17308678, 2956.28673058],
           [2956.28673058, 2948.17308678, 2940.06992852, ..., 2927.93513804,
            2940.06992852, 2948.17308678],
           [2948.17308678, 2940.06992852, 2927.93513804, ..., 2919.85849529,
            2927.93513804, 2940.06992852],
           ...,
           [2940.06992852, 2927.93513804, 2919.85849529, ..., 2907.7633576 ,
            2919.85849529, 2927.93513804],
           [2948.17308678, 2940.06992852, 2927.93513804, ..., 2919.85849529,
            2927.93513804, 2940.06992852],
           [2956.28673058, 2948.17308678, 2940.06992852, ..., 2927.93513804,
            2940.06992852, 2948.17308678]], shape=(400, 400))

  • longitude
    ()
    float64
    6.457

    long_name :

    longitude

    units :

    degrees_east

    standard_name :

    longitude

    array(6.456949)

  • latitude
    ()
    float64
    50.93

    long_name :

    latitude

    units :

    degrees_north

    positive :

    up

    standard_name :

    latitude

    array(50.928727)

  • altitude
    ()
    float64
    310.0

    long_name :

    altitude

    units :

    meters

    standard_name :

    altitude

    array(310.)

  • crs_wkt
    ()
    int64
    0

    crs_wkt :

    PROJCRS["unknown",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",50.9287272,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",6.4569489,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 :

    unknown

    grid_mapping_name :

    azimuthal_equidistant

    latitude_of_projection_origin :

    50.9287272

    longitude_of_projection_origin :

    6.4569489

    false_easting :

    0.0

    false_northing :

    0.0

    array(0)

• Data variables: (18)
  • DBZH
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_h

    long_name :

    Equivalent reflectivity factor H

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • DBZV
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_v

    long_name :

    Equivalent reflectivity factor V

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • KDP
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    deg/km

    units :

    degrees per kilometer

    standard_name :

    radar_specific_differential_phase_hv

    long_name :

    Specific differential phase HV

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • RHOHV
    (y, x)
    float32
    ...

    format :

    UV16

    is_dft :

    0

    units :

    unitless

    standard_name :

    radar_correlation_coefficient_hv

    long_name :

    Correlation coefficient HV

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • DBTH
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_h

    long_name :

    Total power H (uncorrected reflectivity)

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • DBTV
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dBZ

    units :

    dBZ

    standard_name :

    radar_equivalent_reflectivity_factor_v

    long_name :

    Total power V (uncorrected reflectivity)

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • ZDR
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    dB

    units :

    dB

    standard_name :

    radar_differential_reflectivity_hv

    long_name :

    Log differential reflectivity H/V

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • VRADH
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per seconds

    standard_name :

    radial_velocity_of_scatterers_away_from_instrument_h

    long_name :

    Radial velocity of scatterers away from instrument H

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • VRADV
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per second

    standard_name :

    radial_velocity_of_scatterers_away_from_instrument_v

    long_name :

    Radial velocity of scatterers away from instrument V

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • WRADH
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per seconds

    standard_name :

    radar_doppler_spectrum_width_h

    long_name :

    Doppler spectrum width H

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • WRADV
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    m/s

    units :

    meters per second

    standard_name :

    radar_doppler_spectrum_width_v

    long_name :

    Doppler spectrum width V

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • PHIDP
    (y, x)
    float32
    ...

    format :

    UV8

    is_dft :

    0

    unit :

    deg

    units :

    degrees

    standard_name :

    radar_differential_phase_hv

    long_name :

    Differential phase HV

    _Undetect :

    0.0

    [160000 values with dtype=float32]

  • sweep_mode
    ()
    <U20
    ...

    [1 values with dtype=<U20]

  • 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
    ()
    float64
    ...

    standard_name :

    nyquist_velocity

    units :

    m s-1

    [1 values with dtype=float64]

• Attributes: (7)

  source :

  gamic

  ant_gain_h :

  43

  ant_gain_v :

  43

  noise_power_h :

  -3.40133

  noise_power_v :

  -3.13518

  rx_loss_h :

  3

  rx_loss_v :

  3

Plot cartesian representation

actual.DBZH.plot(xlim=(-10e3, 10e3), ylim=(-10e3, 10e3), cmap="HomeyerRainbow", vmin=0)

<matplotlib.collections.QuadMesh at 0x7959c4714550>