如何用R的sf包围绕单点生成6个等距100米的新simple features点？

完全可以用sf包轻松实现，核心思路是通过极坐标转笛卡尔坐标计算出6个等距点的位置，再结合sf的几何操作完成转换。以下是完整示例：

步骤说明 #

• 先确保你的点对象使用投影坐标系（单位为米，比如UTM系列），如果是WGS84（EPSG:4326）这类地理坐标系，需要先转换，否则距离计算会出错。
• 6个等距点的角度间隔为60°（即π/3弧度），从0°开始依次生成。

代码示例 #

# 加载sf包
library(sf)

# 1. 创建示例中心点（替换成你自己的sf对象）
# 这里先创建WGS84坐标系的点，再转成UTM 50N（单位米）
center_point <- st_sfc(st_point(c(116.397, 39.908)), crs = 4326) %>% 
  st_transform(32650)

# 2. 生成6个等距角度（弧度）
angles <- seq(0, 2*pi - pi/3, by = pi/3)

# 3. 提取中心点的坐标
center_coords <- st_coordinates(center_point)
x_center <- center_coords[1, "X"]
y_center <- center_coords[1, "Y"]

# 4. 计算6个新点的坐标
new_coords <- data.frame(
  X = x_center + 100 * cos(angles),
  Y = y_center + 100 * sin(angles)
)

# 5. 转换为sf对象，保持和中心点相同的坐标系
new_points <- st_as_sf(new_coords, coords = c("X", "Y"), crs = st_crs(center_point))

# 6. 可选：如果需要转回WGS84地理坐标系
new_points_wgs84 <- st_transform(new_points, 4326)

# 查看结果
print(new_points)

关键细节 #

• 如果你的原始点已经是投影坐标系（米单位），可以跳过st_transform步骤。
• R的三角函数默认使用弧度，seq(0, 2*pi - pi/3, by = pi/3)会生成0, π/3, 2π/3,...,5π/3共6个角度，确保点均匀分布在圆周上。
• 最终的new_points就是围绕中心点、距离100米、等距分布的6个新点的sf对象。

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