sandbox/extra_functions.R

In andrewpbray/UrbanClimate:

slice_time <- function(x, start_time, end_time, ref_start = "2005-01", ref_end = "2100-12") {
  start_period <- as.period(ymd(paste0(start_time, "-01")) - ymd(paste0(ref_start, "-01")))
  start_ind <- round(day(start_period)/31) + 1
  end_period <- as.period(ymd(paste0(end_time, "-01")) - ymd(paste0(ref_start, "-01")))
  end_ind <- round(day(end_period)/31) + 1
  return(map(x, ~ .x[ , , start_ind:end_ind]))
}

deseason_ts <- function(timeseries) {
  len <- length(timeseries)
  st <- rep(NA, len)
  for (i in 1:12) {
    st[i, 1] <- mean(timeseries[seq(from = i, to = len - 12 + i, by = 12)], na.rm = T)
  }
  st <- rep(st, len/12)
  ds_Ts <- timeseries - st
}