https://github.com/csiro-hydroinformatics/efts
Handling ensemble forecast time series in hydrology, meteorology and possibly other domains
https://github.com/csiro-hydroinformatics/efts
ensemble-prediction hydrologic-database time-series-forecast
Last synced: 5 months ago
JSON representation
Handling ensemble forecast time series in hydrology, meteorology and possibly other domains
- Host: GitHub
- URL: https://github.com/csiro-hydroinformatics/efts
- Owner: csiro-hydroinformatics
- Created: 2018-02-13T05:51:31.000Z (over 8 years ago)
- Default Branch: master
- Last Pushed: 2024-08-30T00:04:04.000Z (almost 2 years ago)
- Last Synced: 2025-09-10T06:43:19.423Z (9 months ago)
- Topics: ensemble-prediction, hydrologic-database, time-series-forecast
- Language: R
- Size: 241 KB
- Stars: 1
- Watchers: 1
- Forks: 1
- Open Issues: 1
-
Metadata Files:
- Readme: README.Rmd
Awesome Lists containing this project
README
---
output: github_document
---
```{r, echo = FALSE}
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.path = "man/figures/README-"
)
```
# efts
Reading and writing Ensemble Forecast Time Series in netCDF files.
## Overview
Plain text files are not well suited to storing the large volumes of data generated for and by ensemble streamflow forecasts with numerical weather prediction models. netCDF is a binary file format developed primarily for climate, ocean and meteorological data. netCDF has traditionally been used to store time slices of gridded data, rather than complete time series of point data. **efts** is for handling the latter.
**efts** is designed to handle netCDF data following the [NetCDF for Water Forecasting Conventions v2.0](https://github.com/jmp75/efts/blob/master/docs/netcdf_for_water_forecasting.md)
## Installation
```{r, eval = FALSE}
# From CRAN:
install.packages("efts")
# Or the the development version from GitHub:
# install.packages("devtools")
devtools::install_github("jmp75/efts")
```
## Example use
This is an extract from one of the package vignettes.
```{r}
library(efts)
ext_data <- system.file('extdata', package='efts')
rain_file <- file.path(ext_data, 'Upper_Murray_sample_rain.nc')
stopifnot(file.exists(rain_file))
rain_dat <- open_efts(rain_file)
class(rain_dat)
```
`rain_dat` has methods to discover and retrieve data in the file.
```{r}
cat(sprintf("This rainfall data set has data for %s stations, the lead time dimension is '%s' because this is not forecast data\n", rain_dat$get_station_count(),
rain_dat$get_lead_time_count()))
rain_dat$get_variable_names()
rain_dat$get_variable_dim_names("rain_der")
```
`rain_der` in this instance has two dimensions, but even if it had been defined as a 3 or 4 dimension data, or in different orders, the method `get_all_series` just does the low-level processing to present a meaninful multivariate `xts` series.
```{r}
d <- rain_dat$get_all_series(variable_name = 'rain_der')
head(d)
```
`xts` may insist on warning that the "timezone of object (UTC) is different than current timezone ().". This is normal and the series is exactly as it should be.
```{r obs_rainfall_plot}
plot(d[1:48], main="Interpolated rainfall (mm/h)")
```
You should use the `close` method once you are done with accessing the data set object and its underlying netCDF file.
```{r}
rain_dat$close()
```
## Related work
Placeholder section, see whether there is an intersect with:
* https://www.r-pkg.org/pkg/ensemblepp
* https://www.r-pkg.org/pkg/tsensembler
* https://www.r-pkg.org/pkg/hyfo