二次b样条曲线插值方法.docxVIP

二次b样条曲线插值方法 1 动态参数化方法 插值是曲线曲线形状和逆向工程中基本而重要的技术。b样条是cagd中曲线曲线形状应用广泛的工具。插值b样条是指将一组数据点指定给一系列数据点，并构建b样条曲线以使序列位于这些数据点。对于指定的一组数据点和端点的限制，插值曲线的确定取决于数据点的参数化和节点向量的选择。在确定数量后，不同的数据点参数化方法和不同的节点向量导致不同形状的插值曲线。因此，要形成一个良好的形状（例如光顺）或b样条的预期形状，必须建立成本的参数化方法。在主要参数化方法中，等距参数化方法使数据点的参数值在参数间隔内均匀分布。输入钩子的参数化方法基于钩子的轮廓，因为钩子的长度类似于曲线的圆弧。参数化方法是基于钩子轮廓的概念。参数化方法是改进累积绳长参数的方法。采用foley参数化方法，充分考虑了绳长与相邻约束之间的角度。节点向量的选择通常符合参数化，但也存在不一致的情况。例如，在文献中，节点向量是通过上述参数化方法进行选择的，相应的b样条基函数最大值的参数值作为数据点的参数值。在这项工作中，我们希望参数化和节点向量的选择一致。 以上提到的参数化方法主要是根据数据点的分布及其几何性质预先确定各数据点的参数值.但由于数据点的分布与插值曲线几何性质之间的关系十分复杂,还未得到很好的研究,所以这些参数化方法时常不能直观地控制B样条插值曲线的几何形状,难使插值曲线具预期的几何性质,导致在有些情况下,这些参数化方法失效,插值曲线形状完全失控. 参数化为构造B样条插值曲线提供了自由度,但在以往的研究中,这些自由度并未得到充分利用.如何充分利用这些自由度构造出具有满意形状的B样条插值曲线是一个困难的问题.本文从相对简单的二次B样条曲线入手,基于插值曲线本身的几何性质来进行参数化,以使插值曲线具有预期的几何性质.不同于以往的参数化方法,本文的方法不是预先选定数据点的参数值,而是在构造插值曲线的过程中根据B样条插值曲线预期的几何约束条件,递推地确定各数据点的参数值,因此称之为动态参数化方法.实验结果表明,本文给出的参数化方法十分有效. 本文的第2节给出二次B样条曲线的有关公式;第3节分别以曲线在各数据点处的切向和曲线各段的相对高度作为插值曲线的几何约束条件,用动态参数化构造相应的二次B样条插值曲线,并给出相应的算法;第4节给出本文方法的应用例子;第5节小结本文的结果. 2 节点向量ttn 二次B样条曲线表示如下: 其中,d0,d1,…,dn为曲线的控制顶点,Ni,2(t)为二次规范B样条基函数,由节点向量T={t0,t1,…,tn+3}确定,t0≤t1≤…≤tn+3.当t0=t1=t2,tn+1=tn+2=tn+3时,相应的二次B样条曲线为端点插值的,即c(t2)=d0,c(tn+1)=dn.本文讨论C1连续的二次B样条曲线,为此下面总是设节点向量T在[t2,tn+1]中无重节点,即t2t3…tn+1. 根据k次B样条基函数的deBoor-Cox公式, 规定 和式(1)直接推算得,二次B样条曲线c(t)在节点ti和参数值处的位置向量可分别表示为 其中,;c(t)的一阶导向量(切向量)可表示为 ti≤tti+1,i=2,3,…,n(4) c(t)在节点ti处的一阶导向量可表示为 3 动态参数化构造顺序 给定平面上n个互异的数据点p1,p2,…,pn,本节将利用动态参数化方法构造顺序通过这些数据点且满足一定的几何约束条件的二次B样条插值曲线. 3.1 次b样条插值曲线 设已知数据点pi处的单位切向量vi,i=1,2,…,n.若不能直接获得vi,可根据数据点预先计算预期的vi,相关的计算方法见文献.记αi为vi到向量pi+1-pi的有向夹角,θi为vi到vi+1的有向夹角,i=1,2,…,n-1(如图1). 定义1. 设u,v为向量,u到v的有向夹角记作〈u,v〉.当u和v方向相反时,〈u,v〉=π;否则,〈u,v〉指从u的正向到v的正向的角度绝对值π的夹角,夹角的符号取逆时针方向为正,顺时针方向为负. 下面要构造如式(1)表示的二次B样条插值曲线c(t),使得 根据B样条曲线的性质和式(2)、式(5),在区间[ti,ti+1]上c(t)为两端分别切于控制多边形的边di-1-di-2和di-di-1且包含在三角形Δdi-2di-1di中的凸曲线段.因此,为了保证满足式(6)、式(7)的二次B样条插值曲线c(t)存在,对单位切向量vi作如下限定是合理的: 构造在各数据点具指定切向的二次B样条插值曲线的方法由下面定理1给出. 定理1. 给定平面上n个互异的数据点p1,p2,…,pn以及n个满足式(8)的单位切向量v1,v2,…,vn,取控制顶点序列d0,d1,…,dn如下: 其中,li为过点pi以vi为方向的直线,取节点向量T如下: 其中,为满