library(hydrotoolbox)
#> Useful information sources:
#> + help(package = "hydrotoolbox")
#> + vignette(package = "hydrotoolbox")
#> + https://gitlab.com/ezetoum27/hydrotoolbox
#> *******************************************
#> Don't forget to cite:
#> + citation(package = "hydrotoolbox")La Autoridad Interjurisdiccional de las Cuencas de los ríos Limay, Neuquén y Negro (AIC), es un organismo que tiene por objeto entender en todo lo relativo a la administración, control, uso y preservación de las cuencas de los ríos mencionados. Entre otras tareas, mantiene operativas una gran cantidad de estaciones nivo-meteorológicas. Estos datos sólo son accesibles mediante acuerdo legal con el AIC.
El paquete hydrotoolbox ofrece la posibilidad de
leer estos archivos (formato .xls) de manera automática
mediante la función read_aic(). Al hacerlo, se cargará al
Global Environment de R un
data.frame con los datos del archivo original. Cabe
destacar que esta función rellena automáticamente los vacíos existentes
entre registros con NA_real_.
# if you want to get the sheet's names
read_aic(path = path, get_sheet = TRUE)
# loading the series in a data.frame
my_data <-
read_aic(path = path,
by = 'day',
sheet = 'CERRO NEVADO')Si bien esta función resulta de gran utilidad, a medida que la cantidad de variables a analizar crece, cargar estas tablas, ordenarlas y modificarlas, se vuelve tarea complicada. La solución que ofrece hydrotoolbox es la de trabajar con los objetos y métodos que el paquete provee. En las siguientes secciones muestro cómo usarlos.
Como menciono en los principios de diseño de este paquete
(vignette('package_overview', package = 'hydrotoolbox')),
los datos que se registran en las estaciones deben almacenarse en un
mismo objeto. Por ello primero habrá que crear dicho objeto (o estación
hidro-meteorológica) y luego usar hm_build(), un método que
permite cargar automáticamente al objeto todas las variables que la
estación real registra.
# set path and file name
my_path <- '/home/.../my_folder'
my_file <- 'AIC_data.xls'
# build the station
nevado <-
hm_create() %>%
hm_build(bureau = 'aic',
path = my_path,
file_name = my_file,
slot_name = c('precip', 'tmax', 'tmin',
'tmean', 'rh', 'wspd',
'wdir', 'swe'),
by = 'day',
sheet = 'CERRO NEVADO')Dado que la función constructora es la única que difiere de lo
desarrollado para los datos del SNIH, recomiendo (re)visitar esta viñeta
(vignette('snih_arg', package = 'hydrotoolbox')).