在R中计算可自定义插值点的2D样条路径的xspline替代方案

以下方案均无需触发绘图设备、支持自定义插值点数量，同时可覆盖更多样条类型调整需求。

方案1：grid包原生xsplinePoints()（零成本迁移） #

该方案完全兼容原有xspline()的shape参数逻辑，不需要修改原有参数配置即可直接替换，grid为R内置基础包，无需额外安装。

library(grid)
library(ggplot2)

# 原有控制点
x <- c(.1, .5, .7, .8)
y <- c(.9, .6, .5, .1)

# 无需调用plot.new()等绘图初始化函数
# 自定义插值点数量
n_interp <- 200
# 构造不绘图的xspline grob
xspline_grob <- xsplineGrob(x = x, y = y, shape = 1, n = n_interp)
# 提取插值坐标
path <- xsplinePoints(xspline_grob)
path_df <- data.frame(x = as.numeric(path$x), y = as.numeric(path$y))

# 绘图验证结果和原xspline一致
ggplot(data = NULL, aes(x, y)) +
  geom_point(data = path_df, size = 0.5) +
  geom_point(data = data.frame(x, y), size = 2, color = "blue")

优势：

• 完全兼容原有xspline的shape参数逻辑，迁移成本为0
• 无任何绘图设备调用要求
• 直接通过n参数自定义插值点数量

方案2：基础包splines实现多类型样条 #

如果需要尝试更多样条类型、精细调整路径形态，可使用R内置的splines包，支持B样条、自然样条、周期性样条等多种类型，可通过阶数、节点位置等参数灵活调整。

library(splines)
library(ggplot2)

# 控制点
x <- c(.1, .5, .7, .8)
y <- c(.9, .6, .5, .1)
# 以控制点序列为参数t，避免x值重复导致的插值异常
t <- seq_along(x)
# 自定义插值点数量
n_interp <- 200
t_interp <- seq(min(t), max(t), length.out = n_interp)

# 拟合3阶B样条，调整df参数可控制平滑度，也可直接指定knots自定义节点位置
spline_x <- lm(x ~ bs(t, df = 4))
spline_y <- lm(y ~ bs(t, df = 4))
# 生成插值坐标
path_df <- data.frame(
  x = predict(spline_x, newdata = data.frame(t = t_interp)),
  y = predict(spline_y, newdata = data.frame(t = t_interp))
)

# 绘图验证
ggplot(data = NULL, aes(x, y)) +
  geom_point(data = path_df, size = 0.5) +
  geom_point(data = data.frame(x, y), size = 2, color = "blue")

优势：

• 支持十余种样条类型，可满足不同的平滑度、边界约束需求
• 完全基于数值计算，无任何绘图相关依赖
• 调整维度丰富，可实现远多于单一shape参数的路径形态控制

方案3：bezier包实现贝塞尔样条 #

如果需要使用贝塞尔曲线类路径，可使用bezier包，适合需要精确控制曲率的场景。

# 首次使用先执行安装：install.packages("bezier")
library(bezier)
library(ggplot2)

x <- c(.1, .5, .7, .8)
y <- c(.9, .6, .5, .1)
ctrl_points <- cbind(x, y)
# 自定义插值点数量
n_interp <- 200
t_interp <- seq(0, 1, length.out = n_interp)
# 生成贝塞尔曲线插值点
path_mat <- bezier(t = t_interp, p = ctrl_points)
path_df <- as.data.frame(path_mat)
colnames(path_df) <- c("x", "y")

# 绘图验证
ggplot(data = NULL, aes(x, y)) +
  geom_point(data = path_df, size = 0.5) +
  geom_point(data = data.frame(x, y), size = 2, color = "blue")

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