R/import_Adata.R

In lvaudor/woody: Use wood occurrence data and discharge data for wood flux prediction

#' Imports data regarding annotation times
#' @param path the path towards the annotation times file
#' @param site the name of the site considered
#' @importFrom magrittr %>%
#' @export
#' @return Adata
#' @examples
#' Adata=import_Adata("data-raw/annot_times_Ain.csv",site="Ain")
import_Adata=function(path, site){
  Adata=readr::read_csv(path) %>%
    dplyr::mutate(start=lubridate::dmy_hms(start),
                  end=lubridate::dmy_hms(end))

  cut_obs_times=function(start,end){
    Time=lubridate::floor_date(start,"hour")
    Timend=lubridate::floor_date(end,"hour")
    if(Time!=Timend){
      seqtimes=seq(from=Time+lubridate::hours(1),to=Timend,by="hours")
      start=c(start,seqtimes)
      end=c(seqtimes,end)
    }
    result=tibble::tibble(start=start,
                          end=end) %>%
      dplyr::mutate(Time=lubridate::floor_date(start,"hour"))
    return(result)
  }
  result=purrr::map2_df(Adata$start,
                        Adata$end,
                        cut_obs_times) %>%
    dplyr::mutate(P_obs=as.numeric(difftime(end,start,units="hours"))) %>%
    dplyr::group_by(Time) %>%
    dplyr::summarise(P_obs=sum(P_obs),
                     .groups="drop")
return(result)
}