pydaymet 0.2.0

Access the climate data from the Daymet database through ORNL's RESTful service.

Package	Description	Status
Hydrodata	Access NWIS, HCDN 2009, NLCD, and SSEBop databases
PyGeoOGC	Send queries to any ArcGIS RESTful-, WMS-, and WFS-based services
PyGeoUtils	Convert responses from PyGeoOGC’s supported web services to datasets
PyNHD	Navigate and subset NHDPlus (MR and HR) using web services
Py3DEP	Access topographic data through National Map’s 3DEP web service
PyDaymet	Access Daymet for daily climate data both single pixel and gridded

PyDaymet: Daily climate data through Daymet

🚨 This package is under heavy development and breaking changes are likely to happen. 🚨

Features

PyDaymet is a part of Hydrodata software stack and provides an interface to access to daily climate data through the Daymet RESTful service. Both single pixel and gridded data can be requested which are returned as pandas.DataFrame for single pixel requests and xarray.Dataset for gridded data requests. Additionally, it can compute Potential EvapoTranspiration (PET) using UN-FAO 56 paper method for both single pixel and gridded data.

You can try using PyDaymet without installing it on you system by clicking on the binder badge below the PyDaymet banner. A Jupyter notebook instance with the Hydrodata software stack pre-installed will be launched in your web browser and you can start coding!

Moreover, requests for additional functionalities can be submitted via issue tracker.

Installation

You can install PyDaymet using pip after installing libgdal on your system (for example, in Ubuntu run sudo apt install libgdal-dev):

$ pip install pydaymet

Alternatively, PyDaymet can be installed from the conda-forge repository using Conda:

$ conda install -c conda-forge pydaymet

Quick start

PyDaymet offers two functions for getting climate data; get_byloc and get_bygeom. The arguments of these functions are identical except the first argument where the latter should be polygon and the former should be a coordinate (a tuple of length two as in (x, y)). The input geometry or coordinate can be in any valid CRS (defaults to EPSG:4326). The dates argument can be either a tuple of length two like (start_str, end_str) or a list of years like [2000, 2005]. It is noted that both functions have a pet flag for computing PET.

from pynhd import NLDI
import pydaymet as daymet

dates = ("2000-01-01", "2000-06-12")
variables = ["prcp", "tmin"]

geometry = NLDI.getfeature_byid("nwissite", "USGS-01031500", basin=True).geometry[0]
clm_g = daymet.get_bygeom(geometry, dates, variables=variables, pet=True)

coords = (-1431147.7928, 318483.4618)
crs = "epsg:3542"
clm_p = daymet.get_byloc(coords, dates, crs=crs, variables=variables, pet=True)

Some example plots are shown below:

Contributing

Contributions are very welcomed. Please read CONTRIBUTING.rst file for instructions.

Credits

Credits to Koen Hufkens for his implementation of accessing the Daymet RESTful service, daymetpy.