使用R语言terra包从数据框创建栅格时如何保留x/y坐标？

核心问题原因 #

terra::setValues()是**按栅格单元格的固有顺序（行优先，从上到下、从左到右）**赋值的，完全不参考数据框里的x/y坐标。如果你的数据框行顺序和栅格单元格顺序不匹配，就会出现赋值错位的问题。

两种正确的坐标匹配赋值方法 #

方法1：用rasterize()（推荐）

这个方法直接基于坐标将点数据的属性值匹配到对应栅格单元格，是最稳妥的方式：

# 1. 将数据框转为空间点矢量对象
library(terra)
pts <- vect(你的数据框名, geom = c("x", "y"), crs = crs(你的空栅格名))

# 2. 将点的预测值栅格化到空栅格中
filled_rast <- rasterize(pts, 你的空栅格名, field = "你的预测值列名")

方法2：先匹配单元格索引再赋值

如果需要严格对应栅格中心坐标，可以先计算每个坐标对应的单元格索引，再按索引顺序赋值：

# 1. 获取空栅格所有单元格的中心坐标
cell_coords <- xyFromCell(你的空栅格名, 1:ncell(你的空栅格名))

# 2. 匹配数据框坐标对应的单元格索引
你的数据框名$cell_idx <- match(
  paste(你的数据框名$x, 你的数据框名$y),
  paste(cell_coords[,1], cell_coords[,2])
)

# 3. 按单元格顺序排序后赋值
filled_rast <- setValues(
  你的空栅格名,
  你的数据框名$你的预测值列名[order(你的数据框名$cell_idx)]
)

注意：这个方法要求数据框的x/y和栅格中心坐标完全一致，否则需要用near()或distance()函数找最近的单元格。

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