NEXRAD Level 2#

[1]:

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

import xradar as xd

Download#

Fetching NEXRAD Level2 radar data file from open-radar-data repository.

[2]:

filename = "../../../../Downloads/KDMX20200810_154746_V06"

[3]:

filename = DATASETS.fetch("KATX20130717_195021_V06")

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

open_nexradlevel2_datatree#

The same works analoguous with the datatree loader. But additionally we can provide a sweep string, number or list.

[12]:

help(xd.io.open_nexradlevel2_datatree)

Help on function open_nexradlevel2_datatree in module xradar.io.backends.nexrad_level2:

open_nexradlevel2_datatree(filename_or_obj, **kwargs)
    Open NEXRAD Level2 dataset as :py:class:`datatree.DataTree`.

    Parameters
    ----------
    filename_or_obj : str, Path, file-like or DataStore
        Strings and Path objects are interpreted as a path to a local or remote
        radar file

    Keyword Arguments
    -----------------
    sweep : int, list of int, optional
        Sweep number(s) to extract, default to first sweep. If None, all sweeps are
        extracted into a list.
    first_dim : str
        Can be ``time`` or ``auto`` first dimension. If set to ``auto``,
        first dimension will be either ``azimuth`` or ``elevation`` depending on
        type of sweep. Defaults to ``auto``.
    reindex_angle : bool or dict
        Defaults to False, no reindexing. Given dict should contain the kwargs to
        reindex_angle. Only invoked if `decode_coord=True`.
    fix_second_angle : bool
        If True, fixes erroneous second angle data. Defaults to ``False``.
    site_coords : bool
        Attach radar site-coordinates to Dataset, defaults to ``True``.
    kwargs : dict
        Additional kwargs are fed to :py:func:`xarray.open_dataset`.

    Returns
    -------
    dtree: datatree.DataTree
        DataTree

[13]:

dtree = xd.io.open_nexradlevel2_datatree(filename, sweep=4)
display(dtree)

Plot Sweep Range vs. Time#

[14]:

dtree["sweep_0"].ds.DBZH.sortby("time").plot(y="time", cmap="HomeyerRainbow")

[14]:

<matplotlib.collections.QuadMesh at 0x7f5c4105a590>

../_images/notebooks_NexradLevel2_22_1.png

Plot Sweep Range vs. Azimuth#

[15]:

dtree["sweep_0"].ds.DBZH.plot(cmap="HomeyerRainbow")

[15]:

<matplotlib.collections.QuadMesh at 0x7f5c2e047b10>

../_images/notebooks_NexradLevel2_24_1.png

[16]:

dtree = xd.io.open_nexradlevel2_datatree(filename, sweep="sweep_8")
display(dtree)

[17]:

dtree = xd.io.open_nexradlevel2_datatree(filename, sweep=[0, 1, 8])
display(dtree)

[18]:

dtree["sweep_0"]["sweep_fixed_angle"].values

[18]:

array(0.48339844)

[19]:

dtree["sweep_1"]["sweep_fixed_angle"].values

[19]:

array(0.48339844)

[20]:

dtree = xd.io.open_nexradlevel2_datatree(
    filename,
)
display(dtree)

[21]:

dtree["sweep_1"]

[ ]: