使用Terra::buffer创建点缓冲区并保留原属性及合并方法咨询

一、创建保留原属性的点缓冲区 #

Terra的buffer()函数默认会保留输入矢量的全部属性，无需额外配置，直接调用即可完成需求。

举个实操例子：

1. 先构造带属性的点矢量数据：

library(terra)

# 创建点坐标
pts <- vect(cbind(c(10, 20, 30), c(40, 50, 60)), crs = "EPSG:4326")
# 添加自定义属性字段
pts$ID <- c(1, 2, 3)
pts$Name <- c("PointA", "PointB", "PointC")

2. 生成指定米数的缓冲区（注意：如果用WGS84地理坐标系（EPSG:4326），需先转换为投影坐标系（如UTM），因为buffer的宽度单位是米，地理坐标系单位为度，无法直接计算米级缓冲区）：

# 转换为对应区域的UTM投影（示例用32N带号，需根据你的点位置调整）
pts_utm <- project(pts, "EPSG:32632")
# 创建1000米半径的缓冲区
buf_with_attr <- buffer(pts_utm, width = 1000)
# 若需要转回地理坐标系，执行以下代码
buf_with_attr_wgs <- project(buf_with_attr, "EPSG:4326")

查看缓冲区属性，会发现原点点的ID、Name等字段完全保留：

head(values(buf_with_attr))

二、无属性缓冲区与原数据集按位置合并 #

如果已经得到无属性的缓冲区，可通过terra::join()函数实现类似GIS“按位置合并”的功能，基于空间交集匹配原数据的属性。

示例步骤：

1. 先模拟无属性缓冲区（比如删除现有缓冲区的属性字段）：

# 模拟无属性缓冲区：清空所有属性
buf_no_attr <- buf_with_attr
values(buf_no_attr) <- NULL

2. 按空间位置合并，将原点点的属性关联到缓冲区上：

# 基于空间交集匹配，保留缓冲区几何，添加原点点的属性
buf_joined <- join(buf_no_attr, pts_utm, type = "intersect")

参数说明：

• type = "intersect"：表示只要缓冲区与原点点空间相交（此处每个缓冲区对应一个点，属于完全包含关系），就将点的属性匹配到缓冲区
• 若存在一对多的匹配场景，可通过how参数控制匹配逻辑，比如how = "left"会保留所有缓冲区，即使没有匹配到对应点

查看合并结果，原属性已成功关联：

head(values(buf_joined))

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