寻求在JavaScript中实现scipy的splprep与splev的方案建议

1. 手动生成适配的节点数组（解决b-spline包的痛点） #

你用b-spline包时缺合适的knot值，本质是因为splprep的核心逻辑之一就是根据输入点自动计算参数化节点，而很多基础B样条库只提供插值/求值能力。可以自己实现这两步：

• 参数化处理：用splprep默认的累积弦长参数化，计算每个点到起点的距离累加，再归一到[0,1]得到参数数组t。
• 生成节点数组：假设用3阶（三次）样条，节点数组需要满足：前3个节点重复参数起点（0），后3个重复参数终点（1），中间节点用t的内部值（去掉首尾的参数点）。示例代码：

// 计算累积弦长参数化
function getParametricPoints(points) {
  let t = [0];
  let total = 0;
  for (let i = 1; i < points.length; i++) {
    const dx = points[i][0] - points[i-1][0];
    const dy = points[i][1] - points[i-1][1];
    total += Math.hypot(dx, dy);
    t.push(total);
  }
  return t.map(v => v / total);
}

// 生成3阶样条的节点数组
function generateKnots(t) {
  const k = 3;
  const knots = [];
  // 前k个节点为0
  for (let i = 0; i < k; i++) knots.push(0);
  // 中间节点用参数数组的内部值
  for (let i = 1; i < t.length - 1; i++) knots.push(t[i]);
  // 后k个节点为1
  for (let i = 0; i < k; i++) knots.push(1);
  return knots;
}

生成knots后，就可以传给b-spline包构建样条，再实现类似splev的求值逻辑（遍历基函数计算加权和）。

2. 用WebAssembly复用Fitpack核心逻辑 #

如果不想自己造轮子，直接复用scipy依赖的Fitpack（Fortran代码）是最准确的方案：

• 找到Fitpack的C语言移植版本，用Emscripten编译成WebAssembly模块。
• 在JS里调用wasm导出的splprep和splev函数，完全复刻scipy的行为，包括平滑参数s的设置、多维度曲线拟合等。

3. 实现最小二乘B样条拟合（处理非插值场景） #

如果需要splprep的拟合功能（而非严格插值，即允许曲线偏离原始点以平滑），可以自己实现最小二乘求解：

• 按步骤1生成节点数组。
• 对每个原始点的参数t_i，计算所有B样条基函数的值，构建矩阵A。
• 用QR分解或伪逆求解线性方程组A * c = y（y是原始点的x/y坐标），得到样条系数c。
• splev的逻辑就是对输入的参数u，计算基函数值并和系数加权求和。

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