基于时空匹配从SpatRaster提取点数据的问题与优化

一、去除data.table匹配后的重复列 #

用data.table做最近时间匹配后出现重复列，本质是匹配过程中保留了源数据（栅格+点）的重复键列（比如空间坐标、时间列）。处理方式很直接：

• 手动删除冗余列：如果重复列是x/y（栅格单元格坐标和点的匹配坐标），或者time（点的时间和栅格的匹配时间），直接用data.table的:=操作符删除不需要的列。
• 提前重命名列：在匹配前给栅格的时间/坐标列重命名（比如把栅格时间改成raster_time），避免和点数据的列名冲突，后续按需保留或删除。

二、更简便的时空匹配方法 #

不用手动写data.table滚动匹配，直接用stars包的st_extract并设置time_match="nearest"，就能自动完成空间+最近时间的匹配，一步到位且不会产生重复列。另外，terra包的extract也支持时空匹配，适合习惯terra生态的用户。

完整示例代码 #

# 加载依赖包
library(terra)
library(data.table)
library(stars)
library(sf)

# ----------------------
# 1. 构造模拟数据
# ----------------------
# 创建带时间维度的SpatRaster（含u、v流速分量）
r <- rast(nrow=10, ncol=10, xmin=0, xmax=10, ymin=0, ymax=10)
times <- seq(as.POSIXct("2023-01-01 00:00:00"), as.POSIXct("2023-01-01 12:00:00"), by="2 hour")
r <- lapply(times, function(t) setValues(r, runif(ncell(r)))) %>% rast() %>% setNames("u")
r$v <- r * 0.8  # 添加v分量
time(r) <- times # 绑定时间维度

# 创建带时间属性的点数据
points_df <- data.frame(
  id = 1:5,
  lon = c(2.5, 4.5, 6.5, 8.5, 1.5),
  lat = c(3.5, 5.5, 7.5, 1.5, 9.5),
  time = seq(as.POSIXct("2023-01-01 01:00:00"), as.POSIXct("2023-01-01 11:00:00"), by="2 hour")
)

# ----------------------
# 方法1：data.table匹配+去重
# ----------------------
# 栅格转长格式data.table，重命名时间列避免冲突
r_dt <- as.data.table(r, xy=TRUE, na.rm=FALSE)
setnames(r_dt, "time", "raster_time")

# 点数据转data.table，获取对应栅格单元格的中心坐标
points_dt <- as.data.table(points_df)
points_dt[, c("x", "y") := as.data.table(xyFromCell(r, cellFromXY(r, .SD[, .(lon, lat)])))]

# 设置匹配键，滚动匹配最近时间
setkey(r_dt, x, y, raster_time)
setkey(points_dt, x, y, time)
matched_dt <- r_dt[points_dt, roll="nearest"]

# 删除重复的空间坐标列，按需保留时间列
matched_dt[, c("x", "y") := NULL]
# 若不需要栅格匹配时间，可删除：matched_dt[, raster_time := NULL]
cat("data.table匹配去重结果:\n")
print(matched_dt)

# ----------------------
# 方法2：stars一键时空匹配（推荐）
# ----------------------
# 点数据转带时间维度的stars对象
points_stars <- st_as_stars(
  st_as_sf(points_df, coords=c("lon", "lat"), crs=crs(r)),
  time = ~time
)

# 提取时指定时间匹配为最近邻
extracted_stars <- st_extract(r, points_stars, time_match = "nearest")
extracted_dt <- as.data.table(extracted_stars)

cat("\nstars一键匹配结果:\n")
print(extracted_dt)

关键说明 #

• stars::st_extract的time_match="nearest"参数是核心：默认精确匹配时间，不匹配则返回NA；设置该参数后自动匹配最近的时间层，无需手动处理时间匹配逻辑。
• 如果习惯用terra，可以用terra::extract结合time参数指定点的时间，通过method="nearest"实现最近时间匹配，逻辑类似stars。

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