R语言实现xy坐标定义的任意折线单向偏移的方法

方案一：基于sf包直接生成单侧偏移（推荐） #

该方案无需拆分双向偏移结果，直接指定偏移方向生成目标曲线，不受原始曲线点间距影响，通用性最强。
核心原理是调用sf::st_buffer()的singleSide参数：

• singleSide = TRUE生成曲线右侧偏移线
• singleSide = FALSE生成曲线左侧偏移线
• dist参数控制偏移距离，正负值可配合方向调整偏移朝向

代码示例：

# 加载依赖包
library(sf)

# 1. 生成示例曲线
t <- seq(10, 0, by = -0.05)
curve_df <- data.frame(
  x = t * cos(t), 
  y = t * sin(t)
)

# 2. 转换为sf线要素对象
curve_sf <- st_sfc(st_linestring(as.matrix(curve_df)), crs = NA)

# 3. 直接生成单侧偏移线
# 右侧偏移0.5单位
offset_right <- st_buffer(
  curve_sf, 
  dist = 0.5, 
  singleSide = TRUE, 
  endCapStyle = "FLAT", 
  joinStyle = "ROUND"
)
# 左侧偏移0.5单位
offset_left <- st_buffer(
  curve_sf, 
  dist = 0.5, 
  singleSide = FALSE, 
  endCapStyle = "FLAT", 
  joinStyle = "ROUND"
)

# 4. 绘图验证
plot(st_geometry(curve_sf), col = "black", lwd = 2)
plot(st_geometry(offset_right), col = "red", lwd = 2, add = TRUE)
plot(st_geometry(offset_left), col = "blue", lwd = 2, add = TRUE)
legend("topleft", 
       legend = c("原曲线", "右侧偏移", "左侧偏移"),
       col = c("black", "red", "blue"),
       lwd = 2, bty = "n")

方案二：基于polyclip双向偏移结果拆分 #

如果需要沿用原有polyclip的计算结果，可通过向量叉乘判断偏移点的相对位置完成拆分，不受点间距限制。
核心原理是计算每个偏移点与原曲线最近线段的叉乘值，根据叉乘正负判断点在曲线的左侧还是右侧，按方向归类即可。

代码示例：

library(polyclip)

# 1. 生成点间距较大的示例曲线（验证拆分效果）
t <- seq(10, 0, by = -0.5)
curve_df <- data.frame(
  x = t * cos(t), 
  y = t * sin(t)
)

# 2. 生成双向偏移结果
offset_raw <- polylineoffset(
  curve_df, 
  delta = 0.5, 
  jointype = "round", 
  endtype = "openbutt"
)[[1]]
offset_df <- as.data.frame(offset_raw)

# 3. 定义叉乘函数判断点相对线段的朝向
point_side <- function(px, py, x1, y1, x2, y2) {
  return((x2 - x1) * (py - y1) - (y2 - y1) * (px - x1))
}

# 4. 给每个偏移点标注所属方向
offset_df$side <- NA
for (i in 1:nrow(offset_df)) {
  px <- offset_df$x[i]
  py <- offset_df$y[i]
  min_dist <- Inf
  curr_side <- 0
  # 遍历原曲线所有线段找到最近的线段
  for (j in 1:(nrow(curve_df)-1)) {
    x1 <- curve_df$x[j]
    y1 <- curve_df$y[j]
    x2 <- curve_df$x[j+1]
    y2 <- curve_df$y[j+1]
    # 计算点到线段的距离
    a <- x2 - x1
    b <- y2 - y1
    c <- px - x1
    d <- py - y1
    dot <- c*a + d*b
    len_sq <- a^2 + b^2
    param <- ifelse(len_sq == 0, 0, max(0, min(1, dot / len_sq)))
    xx <- x1 + param * a
    yy <- y1 + param * b
    dist <- sqrt((px - xx)^2 + (py - yy)^2)
    # 更新最近线段对应的方向
    if (dist < min_dist) {
      min_dist <- dist
      curr_side <- point_side(px, py, x1, y1, x2, y2)
    }
  }
  offset_df$side[i] <- ifelse(curr_side > 0, "left", "right")
}

# 5. 拆分左右偏移线
offset_left <- offset_df[offset_df$side == "left", ]
offset_right <- offset_df[offset_df$side == "right", ]

# 6. 绘图验证
plot(curve_df, type = "l", lwd = 2)
lines(offset_left$x, offset_left$y, col = "blue", lwd = 2)
lines(offset_right$x, offset_right$y, col = "red", lwd = 2)

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