三次B样条曲线拟合算法

1 三次B样条曲线方程

B样条曲线分为近似拟合和插值拟合，所谓近似拟合就是不过特征点，而插值拟合就是通过特征点，但是插值拟合需要经过反算得到控制点再拟合出过特征点的B样条曲线方程。这里会一次介绍两种拟合算法。首先介绍B样条的曲线方程。

B样条曲线的总方程为：

其中P _i 是控制曲线的特征点，F_i,k​(u)则是K阶B样条基函数。

1.1 三次B样条曲线方程中基函数为：

其中 表示阶乘。化成看的明白的式子就是（以四个控制点为例）：

将图片上的基函数代入到方程（1）中，就是：

方程（3）就是三次B样条曲线方程。上式(2)的j是什么意思，其实j就是控制点的索引值。这里我把书上的公式摘抄下来，可能看起来更为清晰。参考书籍：《计算机图形学 第3版》何援军 第13章

三次B样条曲线计算

2 三次B样条曲线近似拟合 

近似拟合很简单。不需要求控制点，求得， 由上述方程（3），代入，就可以得到由这四个点近似拟合的一段三次B样条曲线，起始点在P ₀，终点在P ₁，对于闭合轮廓，最后一段可以取前两点做辅助，拟合实验结果我最后一块给出。这种近似拟合曲线光滑，但是最大不足就是不过特征点，也就是不过P _i ​ ，需要过点需要反求控制点再拟合。

3 三次B样条插值拟合

插值拟合较为复杂。其实也不算是很复杂，找资料过程和理解过程是一个复杂的过程。不过有了前面大神做工作，我们只是借用别人的成果写代码就好了。我给大家看一篇论文，大家可以百度或者去知网搜索，闭合 B 样条曲线控制点的快速求解算法及应用。文章讲解了反求控制点的具体步骤，写的非常详细，基本上贴近代码的那种。大家可以根据这篇论文反求控制点，拟合出来的三次B样条曲线是经过P i 。代码就不放了，很多，可以根据我给的那篇论文直接编写相应代码，有问题可以私信我，知无不言。

4 拟合结果

 

原轮廓

近似拟合轮廓。可以看到没过黑色特征点，只是近似拟合

 

 插值拟。可以看到曲线经过黑色特征点，不过有一些不足之处。

 5 总结

 三次B样条曲线拟合轮廓效果还是可以，较之Beizer（可以参考我博客三次Beizer曲线拟合算法），B样条将一些细节描述的很好，很多细节之处都贴近原轮廓，但是有一些不足之处，可以看到对直线拟合效果不是很好。两篇博客都是关于闭合轮廓的拟合，对于非闭合或者只是一段曲线拟合，还有一种曲线是很好的，《数值分析》提到过，叫三次样条插值拟合，拟合效果很好，我做过拟合一元三次方程曲线，拟合效果跟原曲线非常贴近，不过过程中需要用到追赶法，而追赶法需要满足一个条件，对于闭合曲线三次样条插值是不满足这个条件的，所以我没去深研究，大家可以去试一试。谢谢大家！

 

原文链接：https://blog.csdn.net/liumangmao1314/article/details/54588155