如何避免R语言sf包st_buffer生成的缓冲区边缘出现粗糙锯齿

• st_buffer函数默认的圆弧拟合精度不足：默认nQuadSegs参数值为8，即每四分之一圆弧仅用8段直线近似，当缓冲区半径较大（本次用到100km）时，边缘会出现明显的多边形锯齿
• 坐标系选择不当：如果输入的质心和几何对象使用地理坐标系（如WGS84，单位为经纬度），球面转平面的投影变形会导致缓冲区形状扭曲，同时距离单位也会出现计算偏差

1. 提高缓冲区圆弧拟合精度 #

调整st_buffer的nQuadSegs参数，增大每四分之一圆的拟合线段数，常规场景设置为30~100即可得到光滑的圆形边界，数值越高边缘越平滑，仅会小幅增加计算量。

2. 转换为适配的投影坐标系 #

执行缓冲区运算前，先将所有几何对象转换为平面投影坐标系（单位为米），避免地理坐标系的投影变形问题，可优先选用研究区域对应的UTM投影、地方坐标系等适配性更高的坐标系统。

library(sf)
library(ggplot2)

# 转换为适配的投影坐标系，示例为UTM 50N（EPSG:32650，适合中国东部区域，可替换为你所在区域对应的EPSG编码）
centroid_proj <- st_transform(centroid, crs = 32650)
geometry_proj <- st_transform(geometry, crs = 32650)

# 生成高精度缓冲区，设置nQuadSegs参数提升平滑度
circle <- st_buffer(centroid_proj, dist = 100000, nQuadSegs = 60)

# 执行相交运算
intersection <- st_intersection(geometry_proj, circle)

# 可视化
ggplot() + 
  geom_sf(data = circle, fill = "yellow", alpha = 0.5) +
  geom_sf(data = intersection, fill = "blue", alpha = 0.7)

内容的提问来源于stack exchange，提问作者Murilo Donin de Souza