算法产生背景
Bezier曲线有以下几个不足点,所以导致出现了B-spline算法:
- 一旦确定特征多边形,就确定了曲线的阶次
- Bezier曲线拼接复杂(需要满足几何连续性,参数连续性等)
- Bezier曲线不能作局部修改(只能整体修改)
B-spline算法是整条曲线用一段一段的曲线连接而成,采用分段连续多段式生成
B-spline曲线定义
B-spline曲线定义为:
其中
是特征多边形的顶点;
称为k阶(k-1次)基函数,B-spline算法阶数是次数加1,这是和Bezier算法的一个不同之处;定义域的解释之后会给出,先给出基函数算法。
B-spline基函数的求出算法应用最广泛的是deBoor-cox递推算法:
规定0/0=0
B-spline曲线的定义域为
。设U为所有节点矢量的集合,显见节点表个数为n+k+1个。举例说明,当n=4,k=4时,有
。第一项为
,由deBoor-cox算法可知,其涉及到
到
五个点;第二项同理,设计到
到
五个点;其余同理;所以可以画出区间对应坐标轴:
区间合法所需要的条件为:区间内必须有足够基函数与顶点对应,也即区间中基函数覆盖较多的区间才是一个合法区间。所以上例中对应的合法区间为
也就是
B-spline的类型划分
有均匀B样条曲线,准均匀B样条曲线,分段Bezier曲线,非均匀B样条曲线这几种方法,这里采用均匀B样条曲线的方法绘制曲线。
均匀B样条曲线
定义:当节点沿参数轴均匀等距分布,即
时,为均匀B样条函数,比如:
B-spline基函数(
)的递推公式计算需要用到这里的思想。这里还使用上面当n=4,k=4的例子,从上面的定义可知,我们可以把节点集合
写成这样
也即
,当然U也可以写成其他集合,因为计算的时候是一个比例,所以对计算结果并没有影响,综上就可以简单的算出对应基函数的值。
代码实现
def