sp::gridded的sf等效实现方案技术咨询

一、无需手动指定分辨率的sf等效流程 #

你可以复用原栅格的空间参数（范围、分辨率、CRS），结合sf和terra实现无缝转换，完全不用手动输入分辨率。核心逻辑是保留栅格单元格ID，转回时基于原栅格模板匹配赋值。

完整代码示例 #

# 加载依赖包
library(terra)
library(sf)

# 加载示例栅格（替换为你的文件路径）
r <- rast("SR_50M_alaska_nad.tif")

# 1. 栅格转带单元格ID的sf对象
# 提取坐标、属性及单元格ID
df <- as.data.frame(r, xy = TRUE, cells = TRUE)
# 转为带CRS的sf点对象
sf_points <- st_as_sf(df, coords = c("x", "y"), crs = crs(r))

# 2. 对sf数据框做自定义处理（示例：修改属性值）
sf_points$SR_50M_alaska_nad <- sf_points$SR_50M_alaska_nad + 10

# 3. sf转回栅格（复用原栅格参数）
# 创建与原栅格完全一致的空模板
r_new <- terra::rast(r)
# 通过单元格ID匹配，将处理后的值赋值给新栅格
values(r_new) <- sf_points$SR_50M_alaska_nad[match(cellFromXY(r_new, st_coordinates(sf_points)), sf_points$cell)]

二、关于是否必须替换sp的建议 #

• 短期：sp目前仍可正常运行，若现有代码稳定，无需紧急替换。
• 长期：r-spatial社区已明确将sf、terra、stars作为未来核心技术栈，sp不再新增功能，仅做必要bug修复。从代码可维护性、新包兼容性及未来功能支持角度，建议逐步替换sp相关实现。
• 过渡技巧：优先在新开发代码中采用sf/terra方案，再逐步迁移旧代码，降低重构风险。

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