R/intervalIntersectv2.R

In fjbaron/accelerator: A Tidy Framework For Accelerometry

#' Intersect the intervals given by two dataframes
#' Generates the intersection dataframe resulting from the intersection of the intervals represented by two dataframes of intervals
#'
#' @param interval1 first dataframe
#' @param interval2 second dataframe

#' @return A dataframe of intervals representing the inersection
#'
#' @export
intervalIntersectv2=function(interval1,interval2){
  if(is.null(interval1)) return(interval2 %>% mutate(fromNew=from,toNew=to,day=Sys.Date()) %>% filter(FALSE))
  limits=c(first(interval2$from),last(interval2$to))
  interval1=interval1 %>% mutate(.dummyColumn=1) 
  if(!is.na(limits[1])) interval1=interval1 %>% filter(to>=limits[1],from<=limits[2]) 

  interval2=interval2 %>% mutate(.fromCut=from,.toCut=to) %>% mutate(.dummyColumn=1) %>% select(.fromCut,.toCut,.dummyColumn,everything(),-from,-to)
  resultado=interval1 %>% full_join(interval2,by=".dummyColumn") %>%
    mutate(solapan=as.integer(as.period(lubridate::intersect(interval(from,to),interval(.fromCut,.toCut)),"seconds"))) %>%
    filter(!is.na(solapan)) %>% as_tibble() %>%
    mutate(fromNew=pmax(from,.fromCut),toNew=pmin(to,.toCut)) %>%
    select(from,to,fromNew,toNew,everything(),-.dummyColumn, -.fromCut, -.toCut,-solapan,) %>% 
    filter(fromNew<toNew)
  resultado
}