Learn R Programming

⚠️There's a newer version (0.7.1) of this package.Take me there.

tidyhydat

ß

What does tidyhydat do?

  • Provides functions (hy_*) that access hydrometric data from the HYDAT database, a national archive of Canadian hydrometric data and return tidy data.
  • Provides functions (realtime_*) that access Environment and Climate Change Canada’s real-time hydrometric data source.
  • Provides functions (search_*) that can search through the approximately 7000 stations in the database and aid in generating station vectors
  • Keep functions as simple as possible. For example, for daily flows, the hy_daily_flows() function queries the database, tidies the data and returns a tibble of daily flows.

Installation

You can install tidyhydat from CRAN:

install.packages("tidyhydat")

To install the development version of the tidyhydat package, you can install directly from the rOpenSci development server:

install.packages("tidyhydat", repos = "https://dev.ropensci.org")

Usage

More documentation on tidyhydat can found at the rOpenSci doc page: https://docs.ropensci.org/tidyhydat/

When you install tidyhydat, several other packages will be installed as well. One of those packages, dplyr, is useful for data manipulations and is used regularly here. To use actually use dplyr in a session you must explicitly load it. A helpful dplyr tutorial can be found here.

library(tidyhydat)
library(dplyr)

HYDAT download

To use many of the functions in the tidyhydat package you will need to download a version of the HYDAT database, Environment and Climate Change Canada’s database of historical hydrometric data then tell R where to find the database. Conveniently tidyhydat does all this for you via:

download_hydat()

This downloads (with your permission) the most recent version of HYDAT and then saves it in a location on your computer where tidyhydat’s function will look for it. Do be patient though as this can take a long time! To see where HYDAT was saved you can run hy_default_db(). Now that you have HYDAT downloaded and ready to go, you are all set to begin looking at Canadian hydrometric data.

Real-time

To download real-time data using the datamart we can use approximately the same conventions discussed above. Using realtime_dd() we can easily select specific stations by supplying a station of interest:

realtime_dd(station_number = "08MF005")
#>   Queried on: 2024-01-09 05:08:33.18974 (UTC)
#>   Date range: 2023-12-10 to 2024-01-09 
#> # A tibble: 10,972 × 8
#>    STATION_NUMBER PROV_TERR_STATE_LOC Date               
#>    <chr>          <chr>               <dttm>             
#>  1 08MF005        BC                  2023-12-10 08:00:00
#>  2 08MF005        BC                  2023-12-10 08:05:00
#>  3 08MF005        BC                  2023-12-10 08:10:00
#>  4 08MF005        BC                  2023-12-10 08:15:00
#>  5 08MF005        BC                  2023-12-10 08:20:00
#>  6 08MF005        BC                  2023-12-10 08:25:00
#>  7 08MF005        BC                  2023-12-10 08:30:00
#>  8 08MF005        BC                  2023-12-10 08:35:00
#>  9 08MF005        BC                  2023-12-10 08:40:00
#> 10 08MF005        BC                  2023-12-10 08:45:00
#> # ℹ 10,962 more rows
#> # ℹ 5 more variables: Parameter <chr>, Value <dbl>,
#> #   Grade <chr>, Symbol <chr>, Code <chr>

Or we can use realtime_ws:

realtime_ws(
  station_number = "08MF005",
  parameters = c(46, 5), ## see param_id for a list of codes
  start_date = Sys.Date() - 14,
  end_date = Sys.Date()
)
#> All station successfully retrieved
#> All parameters successfully retrieved
#> # A tibble: 1,323 × 10
#>    STATION_NUMBER Date                Name_En    Value Unit 
#>    <chr>          <dttm>              <chr>      <dbl> <chr>
#>  1 08MF005        2023-12-25 00:00:00 Water tem…  5.14 °C   
#>  2 08MF005        2023-12-25 01:00:00 Water tem…  5.12 °C   
#>  3 08MF005        2023-12-25 02:00:00 Water tem…  5.11 °C   
#>  4 08MF005        2023-12-25 03:00:00 Water tem…  5.1  °C   
#>  5 08MF005        2023-12-25 04:00:00 Water tem…  5.1  °C   
#>  6 08MF005        2023-12-25 05:00:00 Water tem…  5.07 °C   
#>  7 08MF005        2023-12-25 06:00:00 Water tem…  5.06 °C   
#>  8 08MF005        2023-12-25 07:00:00 Water tem…  5.04 °C   
#>  9 08MF005        2023-12-25 08:00:00 Water tem…  5.02 °C   
#> 10 08MF005        2023-12-25 09:00:00 Water tem…  5    °C   
#> # ℹ 1,313 more rows
#> # ℹ 5 more variables: Grade <chr>, Symbol <chr>,
#> #   Approval <chr>, Parameter <dbl>, Code <chr>

Compare realtime_ws and realtime_dd

tidyhydat provides two methods to download realtime data. realtime_dd() provides a function to import .csv files from here. realtime_ws() is an client for a web service hosted by ECCC. realtime_ws() has several difference to realtime_dd(). These include:

  • Speed: The realtime_ws() is much faster for larger queries (i.e. many stations). For single station queries to realtime_dd() is more appropriate.
  • Length of record: realtime_ws() records goes back further in time.
  • Type of parameters: realtime_dd() are restricted to river flow (either flow and level) data. In contrast realtime_ws() can download several different parameters depending on what is available for that station. See data("param_id") for a list and explanation of the parameters.
  • Date/Time filtering: realtime_ws() provides argument to select a date range. Selecting a data range with realtime_dd() is not possible until after all files have been downloaded.

Plotting

Plot methods are also provided to quickly visualize realtime data:

realtime_ex <- realtime_dd(station_number = "08MF005")

plot(realtime_ex)

and also historical data:

hy_ex <- hy_daily_flows(station_number = "08MF005", start_date = "2013-01-01")

plot(hy_ex)

Getting Help or Reporting an Issue

To report bugs/issues/feature requests, please file an issue.

These are very welcome!

How to Contribute

If you would like to contribute to the package, please see our CONTRIBUTING guidelines.

Please note that this project is released with a Contributor Code of Conduct. By participating in this project you agree to abide by its terms.

Citation

Get citation information for tidyhydat in R by running:

To cite package 'tidyhydat' in publications use:

  Albers S (2017). "tidyhydat: Extract and Tidy
  Canadian Hydrometric Data." _The Journal of Open
  Source Software_, *2*(20). doi:10.21105/joss.00511
  <https://doi.org/10.21105/joss.00511>,
  <http://dx.doi.org/10.21105/joss.00511>.

A BibTeX entry for LaTeX users is

  @Article{,
    title = {tidyhydat: Extract and Tidy Canadian Hydrometric Data},
    author = {Sam Albers},
    doi = {10.21105/joss.00511},
    url = {http://dx.doi.org/10.21105/joss.00511},
    year = {2017},
    publisher = {The Open Journal},
    volume = {2},
    number = {20},
    journal = {The Journal of Open Source Software},
  }

License

Copyright 2017 Province of British Columbia

Licensed under the Apache License, Version 2.0 (the “License”); you may not use this file except in compliance with the License. You may obtain a copy of the License at

https://www.apache.org/licenses/LICENSE-2.0

Unless required by applicable law or agreed to in writing, software distributed under the License is distributed on an “AS IS” BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the License for the specific language governing permissions and limitations under the License.

Copy Link

Version

Install

install.packages('tidyhydat')

Monthly Downloads

712

Version

0.6.1

License

Apache License (== 2.0) | file LICENSE

Issues

Pull Requests

Stars

Forks

Maintainer

Sam Albers

Last Published

January 11th, 2024

Functions in tidyhydat (0.6.1)

hy_stn_data_coll

Extract station data collection from HYDAT database
hy_sed_monthly_suscon

Extract monthly flows information from the HYDAT database
hy_data_types

DATA TYPES look-up table
hy_sed_samples

Extract instantaneous sediment sample information from the HYDAT database
hy_reg_office_list

Extract regional office list from HYDAT database
hy_stn_regulation

Extract station regulation from the HYDAT database
hy_stn_op_schedule

Extract station operation schedule from HYDAT database
hy_stn_datum_unrelated

Extract station datum unrelated from HYDAT database
pull_station_number

Convenience function to pull station number from tidyhydat functions
hy_stn_data_range

Extract station data range from HYDAT database
hy_stn_datum_conv

Extract station datum conversions from HYDAT database
plot

Plot historical and realtime data
param_id

Parameter ID
hy_sed_samples_psd

Extract instantaneous sediment sample particle size distribution information from the HYDAT database
realtime_add_local_datetime

Add local datetime column to realtime tibble
hy_version

Extract version number from HYDAT database
realtime_stations

Download a tibble of active realtime stations
reexports

Objects exported from other packages
search_stn_name

A search function for hydrometric station name or number
hy_set_default_db

Set the default database path
hy_dir

Output OS-independent path to the HYDAT sqlite database
realtime_daily_mean

Calculate daily means from higher resolution realtime data
hy_sed_daily_suscon

Extract daily suspended sediment concentration information from the HYDAT database
tidyhydat-package

tidyhydat: Extract and Tidy Canadian 'Hydrometric' Data
realtime_ws

Download realtime data from the ECCC web service
realtime_dd

Download a tibble of realtime river data from the last 30 days from the Meteorological Service of Canada datamart
hy_stn_remarks

Extract station remarks from HYDAT database
hy_sed_monthly_loads

Extract monthly flows information from the HYDAT database
station_choice

Function to chose a station based on consistent arguments for hydat functions.
realtime_plot

Convenience function to plot realtime data
hy_stations

Extract station information from the HYDAT database
hy_src

Open a connection to the HYDAT database
hy_test_db

Get the location of the HYDAT database
tidyeval

Tidy eval helpers
hy_daily_flows

Extract daily flows information from the HYDAT database
download_hydat

Download and set the path to HYDAT
allstations

All Canadian stations
hy_agency_list

hy_agency_list function
hy_annual_stats

Extract annual statistics information from the HYDAT database
hy_data_symbols

DATA SYMBOLS look-up table
hy_daily

Extract all daily water level and flow measurements
hy_daily_levels

Extract daily levels information from the HYDAT database
hy_annual_instant_peaks

Extract annual max/min instantaneous flows and water levels from HYDAT database
hy_remote

Get the version date of HYDAT that is current on the ECCC website
hy_monthly_flows

Extract monthly flows information from the HYDAT database
hy_monthly_levels

Extract monthly levels information from the HYDAT database
hy_plot

This function is deprecated in favour of generic plot methods
hy_sed_daily_loads

Extract daily sediment load information from the HYDAT database
hy_datum_list

Extract datum list from HYDAT database