在B样条(BSpline)中维护上向量(Upvector)的技术求助

问题原因 #

当前代码直接用普通线性插值（LerpN）处理控制点的上向量，但这种插值方式没考虑向量的方向连续性——当相邻两个控制点的上向量方向接近相反时，线性插值会直接穿过向量空间的“反面”，导致180度跳转的突变。虽然后面用叉乘重新对齐切线，但插值阶段的方向断层已经埋下了问题。

修正方案 #

1. 用球面线性插值（SLERP）替代普通Lerp #

SLERP是专门针对单位向量的平滑插值方法，它会沿着球面的最短路径过渡两个向量，从根本上避免方向突变。

2. 确保插值前的向量走最短路径 #

在做SLERP前，先计算两个上向量的点积，如果点积为负，说明两个向量方向接近相反，此时反转其中一个向量，让插值走最短路径。

3. 修正后的代码 #

utype GetUpVector(float thePos)
{
    vtype aTan = GetTangent(thePos).Normal(); // 先把切线归一化，避免后续计算误差

    float aPos = thePos - floorf(thePos);
    int aIPos = (int)floorf(thePos);

    // 取出相邻控制点的上向量并归一化
    utype up0 = mPoints[_clamp(0, aIPos, mPoints.Size()-1)].mUp.Normal();
    utype up1 = mPoints[_clamp(0, aIPos+1, mPoints.Size()-1)].mUp.Normal();

    utype aUp;
    // 检查两个向量的点积，判断是否需要反转其中一个来走最短插值路径
    float dotProduct = up0.Dot(up1);
    if (dotProduct < 0.0f)
    {
        up1 = -up1;
        dotProduct = -dotProduct;
    }

    // 当两个向量几乎同向时，用Lerp替代SLERP，避免数值不稳定
    if (dotProduct > 0.9995f)
    {
        LerpN(aUp, up0, up1, aPos);
    }
    else
    {
        // 计算SLERP的权重
        float theta = acosf(dotProduct);
        float sinTheta = sinf(theta);
        float weight0 = sinf((1.0f - aPos) * theta) / sinTheta;
        float weight1 = sinf(aPos * theta) / sinTheta;
        aUp = up0 * weight0 + up1 * weight1;
    }

    // 最后正交化，确保上向量和切线严格垂直
    vtype aLeft = aTan.Cross(aUp);
    vtype aResult = aLeft.Cross(aTan);
    return aResult.Normal();
}

额外注意事项 #

• 提前把所有控制点的上向量都归一化，减少每次插值时的重复计算。
• 如果你的样条是闭合的，要额外处理首尾控制点的上向量连续性，保证循环处的插值不会出现突变。

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