inst/doc/timezones.R

In ARUtools: Management and Processing of Autonomous Recording Unit (ARU) Data

## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----setup--------------------------------------------------------------------
library(ARUtools)
library(lubridate) # For working with date/times
library(dplyr) # For manipulating data

## -----------------------------------------------------------------------------
s <- example_sites_clean

# Force to a non-UTC timezones
s$date_time_start <- force_tz(s$date_time_start, "America/Toronto")
s$date_time_end <- force_tz(s$date_time_end, "America/Toronto")

## -----------------------------------------------------------------------------
m <- clean_metadata(project_files = example_files)
m <- add_sites(m, s)

## -----------------------------------------------------------------------------
tz(m$date_time)

## -----------------------------------------------------------------------------
m_est <- calc_sun(m, aru_tz = "America/Toronto")

## -----------------------------------------------------------------------------
m_local <- calc_sun(m, aru_tz = "local")

## -----------------------------------------------------------------------------
# Split by timezone
m1 <- filter(m, site_id %in% c("P06_1", "P09_1")) # Get P06_1 and P09_1
m2 <- filter(m, !site_id %in% c("P06_1", "P09_1")) # Get all except the above

# Calculate time to sunrise/sunset individually
m1_cst <- calc_sun(m1, aru_tz = "America/Winnipeg")
m2_est <- calc_sun(m2, aru_tz = "America/Toronto")

# Join them back in
m_joint <- bind_rows(m1_cst, m2_est)

## -----------------------------------------------------------------------------
m_mini <- filter(m, site_id %in% c("P01_1", "P06_1")) |>
  select(aru_id, site_id, longitude, latitude) |>
  distinct() |>
  cross_join(data.frame(date_time = c(
    "2020-05-02 05:00:00",
    "2020-05-02 06:00:00",
    "2020-05-02 07:00:00"
  ))) |>
  mutate(
    date = as_date(date_time),
    path = paste0(aru_id, "_", site_id, "_", hour(date_time), ".csv")
  )

## -----------------------------------------------------------------------------
calc_sun(m_mini, aru_tz = "America/Toronto") |>
  arrange(date_time)

## -----------------------------------------------------------------------------
calc_sun(m_mini, aru_tz = "local") |>
  arrange(date_time)