import pandas as pd
from . import _csv, _TRACK, _netcdf, _tempestextremes, witrack
from . import ibtracs
rename_defaults = dict(
longitude="lon",
latitude="lat",
# Names for MIT netCDF
n_track="track_id",
lon_track="lon",
lat_track="lat",
# Names for CHAZ netCDF
stormID="track_id",
# Names for TRACK netCDF
TRACK_ID="track_id",
)
[docs]
def load(
filename=None,
source=None,
variable_names=None,
rename=dict(),
units=None,
baselon=None,
ibtracs_subset="wmo",
tempest_extremes_unstructured=False,
tempest_extremes_header_str="start",
track_calendar=None,
**kwargs,
):
"""Load track data
The optional parameters for different sources of tracks (currently **IBTrACS**,
**TRACK** and **TempestExtremes**) are named {source}_{parameter} (in lower case),
e.g. "ibtracs_subset".
Parameters
----------
filename : str, optional
The file to be loaded. If `source="ibtracs"`, this is not needed as the data is
either included in huracanpy or downloaded when called. If the filename is
provided for an online IBTrACS subset, then the raw downloaded data will be
saved there.
source : str, optional
If the file is not a CSV or NetCDF (identified by the file extension) then the
source needs to be specified to decide how to load the data
* **track**
* **track.tilt**
* **te**, **tempest**, **tempestextremes**, **uz** (For textual format, not csv)
* **witrack**
* **ibtracs**
* **csv**
* **netcdf**, **nc** (includes support for CHAZ & MIT-Open file provided appropriate variable_names)
variable_names : list of str, optional
When loading data from an ASCII file (TRACK or TempestExtremes), specify the
list of variables that have been added to the tracks. This does not include
variables which are in the track by default (e.g. time/lon/lat/vorticity for
TRACK). If a variable at multiple levels has been added to the tracks, then a
variable name for each level needs to be included (e.g. vorticity profiles
with TRACK). If the variable names are not given, then any additional
variables will be named variable_n, where n goes from 0 to the number of
variables
rename : dict, optional
A mapping of variable names to rename (for netCDF files only). Defaults are
* longitude -> lon
* latitude -> lat
To override any of these defaults, you can pass the same name, e.g.
>>> tracks = huracanpy.load(..., rename=dict(longitude="longitude"))
units : dict, optional
A mapping of variable names to units
>>> tracks = huracanpy.load(..., units=dict(vorticity="s**-1", slp="hPa"))
baselon : scalar, optional
Force the loaded longitudes into the range (baselon, baselon + 360). e.g.
(0, 360) or (-180, 180)
ibtracs_subset : str, default="wmo"
IBTrACS subset. Two offline versions are available:
* **wmo**: Data with the wmo_* variables. The data as reported by the WMO agency
responsible for each basin, so methods are not consistent across basins
* **usa** or **JTWC**: Data with the usa_* variables. The data as recorded by
the USA/Joint Typhoon Warning Centre. Methods are consistent across basins,
but may not be complete.
To download online data, the subsets are the different filesp rovided by IBTrACS.
* **ACTIVE**: TCs currently active
* **ALL**: Entire IBTrACS database
* Specific basins: **EP**, **NA**, **NI**, **SA**, **SI**, **SP**, **WP**
* **last3years**: self-explanatory
* **since1980**: Entire IBTrACS database since 1980 (advent of satellite era,
considered reliable from then on)
tempest_extremes_unstructured : bool, default=False
By default the first two columns in TempestExtremes files are the i, j indices
of the closest gridpoint, but for unstructured grids it is a single lookup index
so there is one less column
tempest_extremes_header_str : str, default="start"
This is an option in the Colin's load function, so I assume this can change
between files
track_calendar : str, optional
When loading data from a TRACK ASCII file, if the data uses a different
calendar to the default :class:`datetime.datetime`, then you can pass this
argument to load the times in as :class:`cftime.datetime` with the given
calendar instead
**kwargs
When loading tracks from a standard files these will be passed to the relevant
load function
* netCDF file - :func:`xarray.open_dataset`
* CSV file - :func:`pandas.read_csv`
* parquet file - :func:`pandas.read_parquet`
For CSV files pandas interprets "NA" as `nan` by default, which is overriden in
this function. To restore the pandas default behavious set
:code:`keep_default_NA=True` and :code:`na_values=[]`
Returns
-------
xarray.Dataset
"""
# Overwrite default arguments with explicit arguments passed to rename by putting
# "rename" second in this dictionary combination
rename = {**rename_defaults, **rename}
# If source is not given, try to derive the right function from the file extension
if source is None:
extension = filename.split(".")[-1]
if extension == "csv":
data = _csv.load(filename, **kwargs)
elif extension == "parquet":
data = _csv.load(filename, load_function=pd.read_parquet, **kwargs)
elif filename.split(".")[-1] == "nc":
data = _netcdf.load(filename, rename, **kwargs)
else:
raise ValueError("Source is set to None and file type is not detected")
# If source is given, use the relevant function
else:
source = source.lower()
if source == "track":
data = _TRACK.load(
filename, calendar=track_calendar, variable_names=variable_names
)
elif source == "track.tilt":
data = _TRACK.load_tilts(
filename,
calendar=track_calendar,
)
elif source in ["csv", "uz"]:
data = _csv.load(filename, **kwargs)
elif source in ["te", "tempest", "tempestextremes"]:
data = _tempestextremes.load(
filename,
variable_names,
tempest_extremes_unstructured,
tempest_extremes_header_str,
)
elif source == "witrack":
data = witrack.load(filename)
elif source == "ibtracs":
data = ibtracs.load(ibtracs_subset, filename, **kwargs)
elif source == "netcdf":
data = _netcdf.load(filename, rename, **kwargs)
else:
raise ValueError(f"Source {source} unsupported or misspelled")
# xarray.Dataset.rename only accepts keys that are actually in the dataset
rename = {key: rename[key] for key in rename if key in data}
if len(rename) > 0:
data = data.rename(rename)
if units is not None:
for varname in units:
data[varname].attrs["units"] = units[varname]
if baselon is not None:
data["lon"] = ((data.lon - baselon) % 360) + baselon
return data
