sf包st_difference结果异常：无法生成指定范围环形缓冲区

你的问题出在浮点精度识别和多要素差集运算的逻辑上：当st_difference处理两个多要素sf对象时，默认逐要素对应计算差集，但如果sf无法精准识别内层缓冲区是对应外层缓冲区的子集，就会直接返回原外层要素。结合你的需求，以下是修正方案：

修正后的代码 #

library(sf)
set.seed(42)

# 直接用sf生成点，避免sp包类型转换的潜在问题
x_coords <- runif(n = 10, min = -100, max = 100)
y_coords <- runif(n = 10, min = -100, max = 100)
PhenomPoints_SF <- st_sf(geometry = st_sfc(st_multipoint(cbind(x_coords, y_coords))))

# 设置明确的笛卡尔坐标系，确保距离计算逻辑一致
st_crs(PhenomPoints_SF) <- 4978

# 生成高平滑度缓冲区，减少浮点精度误差
buffer_zones <- st_buffer(PhenomPoints_SF, dist = 20, nQuadSegs = 50)
buffer_zones_2 <- st_buffer(PhenomPoints_SF, dist = 19, nQuadSegs = 50)

# 用st_snap对齐顶点，确保sf识别内层缓冲区是外层的子集
buffer_zones_2_snapped <- st_snap(buffer_zones_2, buffer_zones, tolerance = 1e-6)
ring <- st_difference(buffer_zones, buffer_zones_2_snapped)

# 绘图验证
plot(buffer_zones, col = "Grey")
plot(buffer_zones_2, col = "Blue", add = TRUE)
plot(ring, col = "Yellow", add = TRUE)

关键修正点 #

• 移除全局赋值：原代码中的<<-会将变量存入全局环境，容易引发变量冲突，改用局部赋值<-更安全。
• 明确CRS设置：无CRS的sf对象可能导致距离计算逻辑模糊，设置笛卡尔坐标系（如EPSG:4978）确保距离计算的一致性。
• 提高缓冲区平滑度：将nQuadSegs从默认30提高到50，让缓冲区边缘更平滑，减少浮点精度带来的识别误差。
• 顶点对齐处理：用st_snap将内层缓冲区的顶点对齐到外层缓冲区的顶点，让sf能精准识别包含关系，从而正确计算差集。

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