Furuno

import xarray as xr
import xradar as xd
from open_radar_data import DATASETS

Download

Fetching Furuno radar data file from open-radar-data repository.

filename_scnx = DATASETS.fetch("2006_20220324_000000_000.scnx.gz")
filename_scn = DATASETS.fetch("0080_20210730_160000_01_02.scn.gz")

Downloading file '2006_20220324_000000_000.scnx.gz' from 'https://github.com/openradar/open-radar-data/raw/main/data/2006_20220324_000000_000.scnx.gz' to '/home/docs/.cache/open-radar-data'.

Downloading file '0080_20210730_160000_01_02.scn.gz' from 'https://github.com/openradar/open-radar-data/raw/main/data/0080_20210730_160000_01_02.scn.gz' to '/home/docs/.cache/open-radar-data'.

xr.open_dataset

Making use of the xarray furuno backend.

scn format

ds = xr.open_dataset(filename_scn, engine="furuno")
display(ds)

<xarray.Dataset> Size: 28MB
Dimensions:            (azimuth: 1376, range: 602)
Coordinates:
    elevation          (azimuth) float32 6kB ...
  * range              (range) float32 2kB 25.0 75.0 ... 3.002e+04 3.008e+04
    time               (azimuth) datetime64[ns] 11kB ...
    longitude          float64 8B ...
    latitude           float64 8B ...
    altitude           float64 8B ...
  * azimuth            (azimuth) float32 6kB 0.21 0.47 0.74 ... 359.7 359.9
Data variables: (12/14)
    RATE               (azimuth, range) float32 3MB ...
    DBZH               (azimuth, range) float32 3MB ...
    VRADH              (azimuth, range) float32 3MB ...
    ZDR                (azimuth, range) float32 3MB ...
    KDP                (azimuth, range) float32 3MB ...
    PHIDP              (azimuth, range) float32 3MB ...
    ...                 ...
    QUAL               (azimuth, range) uint16 2MB ...
    sweep_mode         <U20 80B ...
    sweep_number       int64 8B ...
    prt_mode           <U7 28B ...
    follow_mode        <U7 28B ...
    sweep_fixed_angle  float64 8B ...

xarray.Dataset

• Dimensions:
  • azimuth: 1376
  • range: 602
• Coordinates: (7)
  • elevation
    (azimuth)
    float32
    ...

    standard_name :

    ray_elevation_angle

    long_name :

    elevation_angle_from_horizontal_plane

    units :

    degrees

    axis :

    radial_elevation_coordinate

    [1376 values with dtype=float32]

  • range
    (range)
    float32
    25.0 75.0 ... 3.002e+04 3.008e+04

    units :

    meters

    standard_name :

    projection_range_coordinate

    long_name :

    range_to_measurement_volume

    axis :

    radial_range_coordinate

    meters_between_gates :

    50.0

    spacing_is_constant :

    true

    meters_to_center_of_first_gate :

    25.0

    array([2.5000e+01, 7.5000e+01, 1.2500e+02, ..., 2.9975e+04, 3.0025e+04,
           3.0075e+04], dtype=float32)

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

    standard_name :

    time

    [1376 values with dtype=datetime64[ns]]

  • 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]

  • azimuth
    (azimuth)
    float32
    0.21 0.47 0.74 ... 359.7 359.9

    standard_name :

    ray_azimuth_angle

    long_name :

    azimuth_angle_from_true_north

    units :

    degrees

    axis :

    radial_azimuth_coordinate

    array([2.1000e-01, 4.7000e-01, 7.4000e-01, ..., 3.5942e+02, 3.5968e+02,
           3.5994e+02], dtype=float32)

• Data variables: (14)
  • RATE
    (azimuth, range)
    float32
    ...

    standard_name :

    rainfall_rate

    units :

    mm h-1

    long_name :

    rainfall_rate

    [828352 values with dtype=float32]

  • DBZH
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_equivalent_reflectivity_factor_h

    units :

    dBZ

    long_name :

    Equivalent reflectivity factor H

    [828352 values with dtype=float32]

  • VRADH
    (azimuth, range)
    float32
    ...

    standard_name :

    radial_velocity_of_scatterers_away_from_instrument_h

    units :

    meters per seconds

    long_name :

    Radial velocity of scatterers away from instrument H

    [828352 values with dtype=float32]

  • ZDR
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_differential_reflectivity_hv

    units :

    dB

    long_name :

    Log differential reflectivity H/V

    [828352 values with dtype=float32]

  • KDP
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_specific_differential_phase_hv

    units :

    degrees per kilometer

    long_name :

    Specific differential phase HV

    [828352 values with dtype=float32]

  • PHIDP
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_differential_phase_hv

    units :

    degrees

    long_name :

    Differential phase HV

    [828352 values with dtype=float32]

  • RHOHV
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_correlation_coefficient_hv

    units :

    unitless

    long_name :

    Correlation coefficient HV

    [828352 values with dtype=float32]

  • WRADH
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_doppler_spectrum_width_h

    units :

    meters per seconds

    long_name :

    Doppler spectrum width H

    [828352 values with dtype=float32]

  • QUAL
    (azimuth, range)
    uint16
    ...

    [828352 values with dtype=uint16]

  • 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]

• Indexes: (2)
  • range
    PandasIndex

    PandasIndex(Index([   25.0,    75.0,   125.0,   175.0,   225.0,   275.0,   325.0,   375.0,
             425.0,   475.0,
           ...
           29625.0, 29675.0, 29725.0, 29775.0, 29825.0, 29875.0, 29925.0, 29975.0,
           30025.0, 30075.0],
          dtype='float32', name='range', length=602))

  • azimuth
    PandasIndex

    PandasIndex(Index([0.20999999344348907,  0.4699999988079071,  0.7400000095367432,
                           1.0,  1.2599999904632568,  1.5299999713897705,
            1.7899999618530273,   2.049999952316284,   2.319999933242798,
            2.5799999237060547,
           ...
             357.6199951171875,   357.8800048828125,   358.1399841308594,
             358.4100036621094,  358.66998291015625,  358.92999267578125,
             359.1999816894531,  359.41998291015625,  359.67999267578125,
            359.94000244140625],
          dtype='float32', name='azimuth', length=1376))

• Attributes: (0)

Plot Time vs. Azimuth

ds.azimuth.plot(y="time")

[<matplotlib.lines.Line2D at 0x7fbd56004650>]

Plot Range vs. Time

We need to sort by time and specify the y-coordinate.

ds.DBZH.sortby("time").plot(y="time")

<matplotlib.collections.QuadMesh at 0x7fbd5611b150>

Plot Range vs. Azimuth

ds.DBZH.sortby("azimuth").plot(y="azimuth")

<matplotlib.collections.QuadMesh at 0x7fbd4df46f10>

scnx format

ds = xr.open_dataset(filename_scnx, engine="furuno")
display(ds)

<xarray.Dataset> Size: 23MB
Dimensions:            (azimuth: 722, range: 936)
Coordinates:
    elevation          (azimuth) float32 3kB ...
  * range              (range) float32 4kB 37.5 112.5 ... 7.009e+04 7.016e+04
    time               (azimuth) datetime64[ns] 6kB ...
    longitude          float64 8B ...
    latitude           float64 8B ...
    altitude           float64 8B ...
  * azimuth            (azimuth) float32 3kB 0.19 0.68 1.16 ... 359.2 359.7
Data variables: (12/14)
    RATE               (azimuth, range) float32 3MB ...
    DBZH               (azimuth, range) float32 3MB ...
    VRADH              (azimuth, range) float32 3MB ...
    ZDR                (azimuth, range) float32 3MB ...
    KDP                (azimuth, range) float32 3MB ...
    PHIDP              (azimuth, range) float32 3MB ...
    ...                 ...
    QUAL               (azimuth, range) uint16 1MB ...
    sweep_mode         <U20 80B ...
    sweep_number       int64 8B ...
    prt_mode           <U7 28B ...
    follow_mode        <U7 28B ...
    sweep_fixed_angle  float64 8B ...

xarray.Dataset

• Dimensions:
  • azimuth: 722
  • range: 936
• Coordinates: (7)
  • elevation
    (azimuth)
    float32
    ...

    standard_name :

    ray_elevation_angle

    long_name :

    elevation_angle_from_horizontal_plane

    units :

    degrees

    axis :

    radial_elevation_coordinate

    [722 values with dtype=float32]

  • range
    (range)
    float32
    37.5 112.5 ... 7.009e+04 7.016e+04

    units :

    meters

    standard_name :

    projection_range_coordinate

    long_name :

    range_to_measurement_volume

    axis :

    radial_range_coordinate

    meters_between_gates :

    75.0

    spacing_is_constant :

    true

    meters_to_center_of_first_gate :

    37.5

    array([3.75000e+01, 1.12500e+02, 1.87500e+02, ..., 7.00125e+04, 7.00875e+04,
           7.01625e+04], dtype=float32)

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

    standard_name :

    time

    [722 values with dtype=datetime64[ns]]

  • 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]

  • azimuth
    (azimuth)
    float32
    0.19 0.68 1.16 ... 359.2 359.7

    standard_name :

    ray_azimuth_angle

    long_name :

    azimuth_angle_from_true_north

    units :

    degrees

    axis :

    radial_azimuth_coordinate

    array([1.9000e-01, 6.8000e-01, 1.1600e+00, ..., 3.5866e+02, 3.5918e+02,
           3.5967e+02], dtype=float32)

• Data variables: (14)
  • RATE
    (azimuth, range)
    float32
    ...

    standard_name :

    rainfall_rate

    units :

    mm h-1

    long_name :

    rainfall_rate

    [675792 values with dtype=float32]

  • DBZH
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_equivalent_reflectivity_factor_h

    units :

    dBZ

    long_name :

    Equivalent reflectivity factor H

    [675792 values with dtype=float32]

  • VRADH
    (azimuth, range)
    float32
    ...

    standard_name :

    radial_velocity_of_scatterers_away_from_instrument_h

    units :

    meters per seconds

    long_name :

    Radial velocity of scatterers away from instrument H

    [675792 values with dtype=float32]

  • ZDR
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_differential_reflectivity_hv

    units :

    dB

    long_name :

    Log differential reflectivity H/V

    [675792 values with dtype=float32]

  • KDP
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_specific_differential_phase_hv

    units :

    degrees per kilometer

    long_name :

    Specific differential phase HV

    [675792 values with dtype=float32]

  • PHIDP
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_differential_phase_hv

    units :

    degrees

    long_name :

    Differential phase HV

    [675792 values with dtype=float32]

  • RHOHV
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_correlation_coefficient_hv

    units :

    unitless

    long_name :

    Correlation coefficient HV

    [675792 values with dtype=float32]

  • WRADH
    (azimuth, range)
    float32
    ...

    standard_name :

    radar_doppler_spectrum_width_h

    units :

    meters per seconds

    long_name :

    Doppler spectrum width H

    [675792 values with dtype=float32]

  • QUAL
    (azimuth, range)
    uint16
    ...

    [675792 values with dtype=uint16]

  • 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]

• Indexes: (2)
  • range
    PandasIndex

    PandasIndex(Index([   37.5,   112.5,   187.5,   262.5,   337.5,   412.5,   487.5,   562.5,
             637.5,   712.5,
           ...
           69487.5, 69562.5, 69637.5, 69712.5, 69787.5, 69862.5, 69937.5, 70012.5,
           70087.5, 70162.5],
          dtype='float32', name='range', length=936))

  • azimuth
    PandasIndex

    PandasIndex(Index([0.1899999976158142, 0.6800000071525574,  1.159999966621399,
           1.6899999380111694,  2.169999837875366, 2.6499998569488525,
            3.179999828338623,  3.669999837875366,  4.150000095367432,
            4.679999828338623,
           ...
           355.17999267578125, 355.66998291015625, 356.19000244140625,
           356.67999267578125,  357.1600036621094, 357.69000244140625,
           358.16998291015625,  358.6600036621094, 359.17999267578125,
           359.66998291015625],
          dtype='float32', name='azimuth', length=722))

• Attributes: (0)

Plot Time vs. Azimuth

ds.azimuth.plot(y="time")

[<matplotlib.lines.Line2D at 0x7fbd4c5bf750>]

Plot Range vs. Time

We need to sort by time and specify the y-coordinate.

ds.DBZH.sortby("time").plot(y="time")

<matplotlib.collections.QuadMesh at 0x7fbd4c654b50>

Plot Range vs. Azimuth

ds.DBZH.sortby("azimuth").plot(y="azimuth")

<matplotlib.collections.QuadMesh at 0x7fbd56fded90>

open_furuno_datatree

Furuno scn/scnx files consist only of one sweep. But we might load and combine several sweeps into one DataTree.

dtree = xd.io.open_furuno_datatree(filename_scn)
display(dtree)

<xarray.DatasetView> Size: 232B
Dimensions:              ()
Data variables:
    volume_number        int64 8B 0
    platform_type        <U5 20B 'fixed'
    instrument_type      <U5 20B 'radar'
    time_coverage_start  <U20 80B '2021-07-30T16:00:00Z'
    time_coverage_end    <U20 80B '2021-07-30T16:00:14Z'
    longitude            float64 8B 15.45
    altitude             float64 8B 407.9
    latitude             float64 8B 47.08
Attributes:
    Conventions:      None
    instrument_name:  None
    version:          None
    title:            None
    institution:      None
    references:       None
    source:           None
    history:          None
    comment:          im/exported using xradar

datatree.DataTree

• Groups: (1)
  <xarray.DatasetView> Size: 28MB
  Dimensions:            (azimuth: 1376, range: 602)
  Coordinates:
      elevation          (azimuth) float32 6kB ...
    * range              (range) float32 2kB 25.0 75.0 ... 3.002e+04 3.008e+04
      time               (azimuth) datetime64[ns] 11kB 2021-07-30T16:00:06.2777...
      longitude          float64 8B ...
      latitude           float64 8B ...
      altitude           float64 8B ...
    * azimuth            (azimuth) float32 6kB 0.21 0.47 0.74 ... 359.7 359.9
  Data variables: (12/14)
      RATE               (azimuth, range) float32 3MB ...
      DBZH               (azimuth, range) float32 3MB ...
      VRADH              (azimuth, range) float32 3MB ...
      ZDR                (azimuth, range) float32 3MB ...
      KDP                (azimuth, range) float32 3MB ...
      PHIDP              (azimuth, range) float32 3MB ...
      ...                 ...
      QUAL               (azimuth, range) uint16 2MB ...
      sweep_mode         <U20 80B ...
      sweep_number       int64 8B ...
      prt_mode           <U7 28B ...
      follow_mode        <U7 28B ...
      sweep_fixed_angle  float64 8B ...

  sweep_0

  • Groups: (0)
  • Dimensions:
    • azimuth: 1376
    • range: 602
  • Coordinates: (7)
    • elevation
      (azimuth)
      float32
      ...

      standard_name :

      ray_elevation_angle

      long_name :

      elevation_angle_from_horizontal_plane

      units :

      degrees

      axis :

      radial_elevation_coordinate

      [1376 values with dtype=float32]

    • range
      (range)
      float32
      25.0 75.0 ... 3.002e+04 3.008e+04

      units :

      meters

      standard_name :

      projection_range_coordinate

      long_name :

      range_to_measurement_volume

      axis :

      radial_range_coordinate

      meters_between_gates :

      50.0

      spacing_is_constant :

      true

      meters_to_center_of_first_gate :

      25.0

      array([2.5000e+01, 7.5000e+01, 1.2500e+02, ..., 2.9975e+04, 3.0025e+04,
             3.0075e+04], dtype=float32)

    • time
      (azimuth)
      datetime64[ns]
      2021-07-30T16:00:06.277723500 .....

      standard_name :

      time

      array(['2021-07-30T16:00:06.277723500', '2021-07-30T16:00:06.288556500',
             '2021-07-30T16:00:06.299389500', ..., '2021-07-30T16:00:06.245224500',
             '2021-07-30T16:00:06.256057500', '2021-07-30T16:00:06.266890500'],
            dtype='datetime64[ns]')

    • 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]

    • azimuth
      (azimuth)
      float32
      0.21 0.47 0.74 ... 359.7 359.9

      standard_name :

      ray_azimuth_angle

      long_name :

      azimuth_angle_from_true_north

      units :

      degrees

      axis :

      radial_azimuth_coordinate

      array([2.1000e-01, 4.7000e-01, 7.4000e-01, ..., 3.5942e+02, 3.5968e+02,
             3.5994e+02], dtype=float32)

  • Data variables: (14)
    • RATE
      (azimuth, range)
      float32
      ...

      standard_name :

      rainfall_rate

      units :

      mm h-1

      long_name :

      rainfall_rate

      [828352 values with dtype=float32]

    • DBZH
      (azimuth, range)
      float32
      ...

      standard_name :

      radar_equivalent_reflectivity_factor_h

      units :

      dBZ

      long_name :

      Equivalent reflectivity factor H

      [828352 values with dtype=float32]

    • VRADH
      (azimuth, range)
      float32
      ...

      standard_name :

      radial_velocity_of_scatterers_away_from_instrument_h

      units :

      meters per seconds

      long_name :

      Radial velocity of scatterers away from instrument H

      [828352 values with dtype=float32]

    • ZDR
      (azimuth, range)
      float32
      ...

      standard_name :

      radar_differential_reflectivity_hv

      units :

      dB

      long_name :

      Log differential reflectivity H/V

      [828352 values with dtype=float32]

    • KDP
      (azimuth, range)
      float32
      ...

      standard_name :

      radar_specific_differential_phase_hv

      units :

      degrees per kilometer

      long_name :

      Specific differential phase HV

      [828352 values with dtype=float32]

    • PHIDP
      (azimuth, range)
      float32
      ...

      standard_name :

      radar_differential_phase_hv

      units :

      degrees

      long_name :

      Differential phase HV

      [828352 values with dtype=float32]

    • RHOHV
      (azimuth, range)
      float32
      ...

      standard_name :

      radar_correlation_coefficient_hv

      units :

      unitless

      long_name :

      Correlation coefficient HV

      [828352 values with dtype=float32]

    • WRADH
      (azimuth, range)
      float32
      ...

      standard_name :

      radar_doppler_spectrum_width_h

      units :

      meters per seconds

      long_name :

      Doppler spectrum width H

      [828352 values with dtype=float32]

    • QUAL
      (azimuth, range)
      uint16
      ...

      [828352 values with dtype=uint16]

    • 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]

  • Attributes: (0)

• Dimensions:
• Coordinates: (0)
• Data variables: (8)
  • volume_number
    ()
    int64
    0

    array(0)

  • platform_type
    ()
    <U5
    'fixed'

    array('fixed', dtype='<U5')

  • instrument_type
    ()
    <U5
    'radar'

    array('radar', dtype='<U5')

  • time_coverage_start
    ()
    <U20
    '2021-07-30T16:00:00Z'

    array('2021-07-30T16:00:00Z', dtype='<U20')

  • time_coverage_end
    ()
    <U20
    '2021-07-30T16:00:14Z'

    array('2021-07-30T16:00:14Z', dtype='<U20')

  • longitude
    ()
    float64
    15.45

    long_name :

    longitude

    units :

    degrees_east

    standard_name :

    longitude

    array(15.44729)

  • altitude
    ()
    float64
    407.9

    long_name :

    altitude

    units :

    meters

    standard_name :

    altitude

    array(407.9)

  • latitude
    ()
    float64
    47.08

    long_name :

    latitude

    units :

    degrees_north

    positive :

    up

    standard_name :

    latitude

    array(47.07734)

• Attributes: (9)

  Conventions :

  None

  instrument_name :

  None

  version :

  None

  title :

  None

  institution :

  None

  references :

  None

  source :

  None

  history :

  None

  comment :

  im/exported using xradar

Plot Sweep Range vs. Time

dtree["sweep_0"].ds.DBZH.plot()

<matplotlib.collections.QuadMesh at 0x7fbd56f0a1d0>

Plot Sweep Range vs. Azimuth

dtree["sweep_0"].ds.DBZH.sortby("azimuth").plot(y="azimuth")

<matplotlib.collections.QuadMesh at 0x7fbd56ddb990>