if (FALSE) {
# you can download the data from:
# https://gitlab.com/ezetoum27/hydrotoolbox/-/tree/master/my_data
# set the data path
my_path <- "./home/my_folder/my_data"
#///////////////////////////////
# Rectangular data
# txt, csv, csv2 and others.
# See readr package.
#///////////////////////////////
library(readr)
#* Case 1: single file - many numeric variables
hm_create() %>%
hm_build_generic(path = my_path,
file_name = "ianigla_cuevas.csv",
slot_name = c("tair", "rh", "patm",
"precip", "wspd", "wdir",
"kin", "hsnow", "tsoil"),
by = c("hour"),
FUN = read_csv,
col_select = !Est & !YJday & !hh.mm.ss & !bat.Volts
) %>%
hm_show()
#* Case 2: multiple files (one per observation)
hm_create() %>%
hm_build_generic(path = my_path,
file_name = c("h_relativa_cuevas.csv",
"p_atm_cuevas.csv",
"precip_total_cuevas.csv",
"temp_aire_cuevas.csv",
"vel_viento_cuevas.csv"),
slot_name = c("rh", "patm", "precip",
"tair", "wspd"),
by = c("hour", "45 min", "30 min", "1 hour", "15 min"),
FUN = read_csv ) %>%
hm_show()
#///////////////////////////////
# Excel files
# Recommended package => readxl
#///////////////////////////////
library(readxl)
#* Case 1: single file - one sheet - many numeric variables
hm_create() %>%
hm_build_generic(path = my_path,
file_name = "mnemos_guido.xlsx",
slot_name = c("qd"),
by = c("day"),
FUN = read_excel,
sheet = 1L,
skip = 3,
out_name = list("q_m3/s")
) %>% hm_show()
#* Case 2: single file - multiple sheets (one per variable)
hm_create() %>%
hm_build_generic(path = my_path,
file_name = "mnemos_guido.xlsx",
slot_name = c("qd", "evap", "tair",
"tmax", "tmin"),
by = c(q = "day", evap = "day", tair = "6 hour",
tmax = "day", tmin = "day"),
FUN = read_excel,
sheet = c(1L:5L),
skip = 3,
out_name = list( c("q_m3/s", "flag"),
c("evap_mm", "flag"),
c("tair", "flag"),
c("tmax", "flag"),
c("tmin", "flag")
)
) %>%
hm_show()
#* Case 3: multiple files - one sheet per file
hm_create() %>%
hm_build_generic(path = my_path,
file_name = c("discharge_daily.xlsx",
"air_teperature_subdaily.xlsx"),
slot_name = c("qd", "tair"),
by = c(q = "day", tair = "6 hour"),
FUN = read_excel,
sheet = c(1L, 1L),
skip = 3,
out_name = list( c("q_m3/s", "flag"),
c("tair", "flag"))
) %>%
hm_show()
}
Run the code above in your browser using DataLab