在R语言中匹配数据框日期与日期向量中最近日期的方法

我有一个包含日期列的数据框，以及一个向量形式的参考日期列表。已生成可复现的数据框与向量（r对象仅用于创建数据框，可忽略）：

i=500;j=500;n=50
r <- do.call(stack, replicate(n,raster(matrix(runif(i*j), i, j))))
( dates <- seq(as.Date("2020/10/23"), by = 16, length.out = 20) ) 

s <- as(sampleRandom(r[[1]], 100, sp=TRUE), "sf")
df <-  as.data.frame(s)
df$date <- sample(dates, nrow(s), replace=TRUE)

需求：为数据框每一行找到参考日期向量中与之最近的日期，并将该日期添加为数据框的新列。

补充说明：此日期匹配将用于把向量日期对应到栅格堆栈中的相应栅格，进而基于日期提取数据——这也是s最初为sf类型、保留r对象的原因。

方法一：Base R 原生实现 #

无需额外包，通过计算日期差的绝对值找到最近的参考日期：

# 确保date列是标准Date类型（若已为Date类型可跳过）
df$date <- as.Date(df$date)

# 添加最近参考日期列
df$nearest_date <- dates[apply(df$date, 1, function(x) {
  which.min(abs(difftime(x, dates, units = "days")))
})]

方法二：用lubridate包简化计算 #

如果已安装lubridate包，日期运算会更直观：

library(lubridate)

df$nearest_date <- dates[apply(df$date, 1, function(x) {
  which.min(abs(x - dates))
})]

针对栅格提取的优化建议 #

考虑到你最终要从栅格堆栈r中提取数据，可以直接保存匹配的日期索引，避免重复计算：

# 获取每行对应的参考日期在dates向量中的索引
df$date_index <- apply(df$date, 1, function(x) {
  which.min(abs(difftime(x, dates, units = "days")))
})

# 利用索引直接从栅格堆栈提取对应栅格的像素值
df$raster_value <- extract(r, s, index = df$date_index)

内容的提问来源于stack exchange，提问作者jwell123