在R中将每日气象周转换为标准气象周

Question

我在SO中看到了很多关于使用xts，zoo或lubridate包将每日数据转换为每周数据的问题。没有一个答案适合我的问题。我尝试了以下代码

library(zoo)
library(lubridate)
library(xts)
library(tidyverse)
#Calculation for multistation
set.seed(123)
df <- data.frame("date"= seq(from = as.Date("1970-1-1"), to = as.Date("2000-12-31"), by = "day"),
"Station1" = runif(length(seq.Date(as.Date("1970-1-1"), as.Date("2000-12-31"), "days")), 10, 30),
"Station2" = runif(length(seq.Date(as.Date("1970-1-1"), as.Date("2000-12-31"), "days")), 11, 29),
"Station3" = runif(length(seq.Date(as.Date("1970-1-1"), as.Date("2000-12-31"), "days")), 9, 28))
head(df)
# Aggregate over week
df %>% 
mutate(Week = week(ymd(date)),
Year = year(ymd(date))) %>% 
pivot_longer(-c(Week, date, Year), values_to = "value", names_to = "Station") %>% 
group_by(Year, Week, Station) %>% 
summarise(Weekly = mean(value)) %>% 
arrange(Station) %>% 
print(n = 55)

从输出中你可以看到 1970 年包含 53 周，我不想要。我想从每年的第一天开始一周，第 52 周在非闰年应该有 8 天，如果是闰年，第 9 周和第 52 周应该有 8 天，这样每年只包含 52 周。如何在 R 中做到这一点？

Answer 1

为什么不写一个函数，根据你给出的定义给出气象周？包lubridate将为您提供带有yday的一年中的某一天，它可以充当正确周标签向量的索引。这些很容易用简单的模块化数学和串联来构造。

然后，您只需要弄清楚您是否处于闰年，这也可以使用lubridate::leap_year.将它们组合成一个ifelse，你就有一个易于使用的功能：

met_week <- function(dates)
{
normal_year <- c((0:363 %/% 7 + 1), 52)
leap_year   <- c(normal_year[1:59], 9, normal_year[60:365])
year_day    <- lubridate::yday(dates)
return(ifelse(lubridate::leap_year(dates), leap_year[year_day], normal_year[year_day]))
}

你可以做到

df %>% mutate(week = met_week(date))

Answer 2

您可以在一年中的某一天手动执行此操作，不确定是否已经为此构建了函数。

df %>% 
mutate(Week = pmin(52, ceiling(yday(date) / 7)),
Year = year(ymd(date)))