使用Math.Net在三次样条曲线上查找参数t对应的点

在Math.Net Numerics中，没有直接的二维三次样条类，但可以通过分别对x坐标和y坐标基于参数t做一维三次样条插值，实现传入t返回对应二维点的功能。以下是具体实现步骤：

1. 准备基础数据 #

首先定义原始控制点，并生成参数t的数组（t可以是索引或累计弧长，后者能让曲线上的点分布更均匀）：

using MathNet.Numerics.Interpolation;
using System.Drawing; // WPF场景改用System.Windows.Point

// 原始控制点
var controlPoints = new List<PointF>
{
    new PointF(0, 0),
    new PointF(1, 3),
    new PointF(2, 1),
    new PointF(3, 4),
    new PointF(4, 2)
};

// 方式1：用索引作为参数t（简单直观）
var tValues = Enumerable.Range(0, controlPoints.Count).Select(i => (double)i).ToArray();

// 方式2：用累计弧长作为参数t（点分布更均匀）
// var tValues = new List<double> { 0 };
// double totalLength = 0;
// for (int i = 1; i < controlPoints.Count; i++)
// {
//     var dx = controlPoints[i].X - controlPoints[i-1].X;
//     var dy = controlPoints[i].Y - controlPoints[i-1].Y;
//     totalLength += Math.Sqrt(dx * dx + dy * dy);
//     tValues.Add(totalLength);
// }

// 提取x、y坐标数组
var xCoords = controlPoints.Select(p => (double)p.X).ToArray();
var yCoords = controlPoints.Select(p => (double)p.Y).ToArray();

2. 创建三次样条插值实例 #

使用CubicSpline分别对x和y坐标做自然样条插值（也可根据需求选择其他边界条件）：

// 创建x(t)和y(t)的三次样条
var splineX = CubicSpline.InterpolateNatural(tValues, xCoords);
var splineY = CubicSpline.InterpolateNatural(tValues, yCoords);

3. 实现传入t返回点的函数 #

通过样条的Evaluate方法计算对应t的x、y值，组合成二维点：

// 根据参数t获取曲线上的点
PointF GetSplinePoint(double t)
{
    double x = splineX.Evaluate(t);
    double y = splineY.Evaluate(t);
    return new PointF((float)x, (float)y);
}

4. 生成大量点用于绘制 #

通过遍历t的取值范围，生成足够多的点后用线段连接：

// 控制点的密度，step越小点越密集
double step = 0.01;
double startT = tValues[0];
double endT = tValues[^1];

var splineDrawPoints = new List<PointF>();
for (double t = startT; t <= endT; t += step)
{
    splineDrawPoints.Add(GetSplinePoint(t));
}

// 后续可通过GDI+/WPF等绘图库连接splineDrawPoints中的点完成绘制

注意：t的取值需在样条的定义域内（即tValues的最小值到最大值之间），超出范围会触发线性外插逻辑。

内容的提问来源于stack exchange，提问作者Samuel N Lewis