在R中栅格化线要素：如何为栅格像元分配最大值？

问题背景 #

需要栅格化空间线要素，为每个被线接触/相交的栅格像元分配对应所有线的最大值。使用terra包处理SpatVect格式的线要素，优先采用terra::rasterize方案，也接受stars::st_rasterize方案。目前发现terra::rasterize的fun参数仅支持点矢量数据，对线要素不生效，且不想使用基于sp对象的老旧方法raster::rasterize。

问题复现示例 #

library("terra")

# 示例数据
f <- system.file("ex/lux.shp", package="terra")
v <- vect(f)
lns <- as.lines(v)
r <- rast(v, res=.2)

# 尝试用terra::rasterize的max/mean函数
x_max <- rasterize(lns, r, fun="max", field = "POP")
x_mean <- rasterize(lns, r, fun="mean", field = "POP")

# 绘图查看结果
fontsize <- 0.7
par(mfrow=c(1,3))
plot(lns, y = "POP", main = "带原始值的线要素")
text(lns, "POP", cex = fontsize)
plot(x_max, main = "使用fun='max'栅格化结果")
text(x_max, cex = fontsize)
plot(lns, add = T)
plot(x_mean, main = "使用fun='mean'栅格化结果")
text(x_mean, cex = fontsize)
plot(lns, add = T)

（注：运行上述代码后，x_max和x_mean未按预期计算每个像元的最大值/平均值，仅保留了单条线的属性值）

可行解决方案 #

方案1：基于terra包的提取+聚合赋值 #

利用terra::extract获取每个栅格像元对应的所有线要素属性值，再按像元计算最大值，最后赋值到目标栅格：

library(terra)

# 1. 获取每个栅格像元对应的线要素POP值，同时记录像元ID
extracted <- extract(lns, r, field = "POP", cells = TRUE)

# 2. 按像元ID分组，计算每组的最大值
max_vals <- tapply(extracted$POP, extracted$cell, max)

# 3. 创建结果栅格并赋值
r_max <- r
r_max[as.integer(names(max_vals))] <- max_vals

# 查看结果
plot(r_max, main = "像元对应线要素的POP最大值")
plot(lns, add = TRUE)

方案2：基于stars包的直接栅格化 #

使用stars包的st_rasterize，原生支持对线要素栅格化时指定聚合函数：

library(stars)
library(sf)

# 将SpatVect转换为sf对象
lns_sf <- st_as_sf(lns)

# 创建目标栅格模板
stars_r <- st_as_stars(st_bbox(lns_sf), dx = 0.2, dy = 0.2)

# 栅格化并计算每个像元的最大值
rasterized_max <- st_rasterize(lns_sf["POP"], stars_r, fun = max)

# 查看结果
plot(rasterized_max, main = "stars栅格化的POP最大值")
plot(st_geometry(lns_sf), add = TRUE)

说明 #

• terra::rasterize目前对线要素确实不支持fun参数的聚合计算，方案1通过先提取后聚合的方式绕开了这一限制，完全基于terra生态实现。
• 方案2利用stars包的原生支持，代码更简洁，直接在栅格化步骤完成聚合计算。

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