求适用于非规整XYZ点集的Java平滑3D多项式插值方案

核心思路：分量式参数化插值 #

不管是2D还是3D曲线插值，核心逻辑都是用统一的参数t来驱动三个坐标分量的独立插值，完全不需要要求x/y/z坐标严格递增——参数t只和点的排列顺序挂钩，和坐标值无关。这完美解决你提到的y/z坐标不单调的问题。

1. Catmull-Rom样条（最推荐） #

你之前了解的Catmull-Rom完全适配3D场景，实现起来也简单：

• 把你的点按预期的曲线顺序排列成(x0,y0,z0), (x1,y1,z1), ..., (xn,yn,zn)
• 生成参数序列t：可以用均匀参数（t从0到n，每个点对应t=i），或者更平滑的弧长参数化（t的增量等于相邻点之间的欧氏距离，让曲线的速度更均匀）
• 对x、y、z三个数组分别应用Catmull-Rom插值，用同一个t值计算三个分量的插值结果，组合起来就是3D曲线上的点

Java实现选项：

• JOML库：直接提供CatmullRomSpline类，支持3D向量插值。示例代码片段：

import org.joml.Vector3f;
import org.joml.CatmullRomSpline;

// 初始化点数组
Vector3f[] points = new Vector3f[] {
    new Vector3f(x0, y0, z0),
    new Vector3f(x1, y1, z1),
    // ... 其他点
};
// 创建样条，设置是否闭合曲线（false为开放）
CatmullRomSpline spline = new CatmullRomSpline(points, false);

// 插值获取t位置的点（t范围0到1，对应从第一个点到最后一个点）
Vector3f interpolated = new Vector3f();
spline.interpolate(t, interpolated);
float x = interpolated.x;
float y = interpolated.y;
float z = interpolated.z;

• 轻量手动实现：如果不想依赖库，核心插值公式很简单，针对单个分量（比如x）的Catmull-Rom计算：

public static float catmullRom(float t, float p0, float p1, float p2, float p3) {
    float t2 = t * t;
    float t3 = t2 * t;
    return 0.5f * ((2 * p1) +
                   (-p0 + p2) * t +
                   (2 * p0 - 5 * p1 + 4 * p2 - p3) * t2 +
                   (-p0 + 3 * p1 - 3 * p2 + p3) * t3);
}

对x、y、z分别调用这个函数，传入对应点的分量即可。

2. 三次样条插值的3D扩展 #

你之前用的2D三次样条可以直接扩展到3D，同样采用参数化思路：

• 不用x作为参数，而是用基于点顺序的参数t（和Catmull-Rom的参数逻辑一致）
• 对x、y、z三个分量分别使用三次样条插值，用同一个t值计算结果

Java实现选项：

• Apache Commons Math：用SplineInterpolator类分别处理x、y、z数组：

import org.apache.commons.math3.analysis.interpolation.SplineInterpolator;
import org.apache.commons.math3.analysis.polynomials.PolynomialSplineFunction;

// 生成参数t数组（这里用均匀参数）
double[] t = new double[points.length];
for (int i = 0; i < t.length; i++) {
    t[i] = i;
}

// 分别创建三个分量的样条函数
SplineInterpolator interpolator = new SplineInterpolator();
PolynomialSplineFunction xSpline = interpolator.interpolate(t, xArray);
PolynomialSplineFunction ySpline = interpolator.interpolate(t, yArray);
PolynomialSplineFunction zSpline = interpolator.interpolate(t, zArray);

// 插值获取t位置的点
double interpolatedX = xSpline.value(tVal);
double interpolatedY = ySpline.value(tVal);
double interpolatedZ = zSpline.value(tVal);

关于Tricubic插值的说明 #

你提到的Tricubic插值是用于3D网格体数据的插值（比如三维空间中的标量场），完全不适合单条3D曲线的场景，直接排除即可。

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