NURBS曲线曲面

$NURBS$曲线曲面

曲线定义：

$NURBS$曲线通常称为非均匀有理B样条，其数学定义如下：
$$P(t)=\frac{{\sum }_{i=0}^n{N}_{i,k}(t){\omega }_i{P}_i}{{\sum }_{i=0}^n{N}_{i,k}(t){\omega }_i},t\in [0,1]$$
式中，${\omega }_i,i=0,1,...,n$称为权因子，分别与控制点$P_i,i=0,1,...,n$相联系。首末权因子${\omega }_0,{\omega }_n>0$,其余${\omega }_i\ge 0$,且顺序k个权因子不同时为零。把首末权因子都等于1的$NURBS$曲线称为标准型$NURBS$曲线，否则称为非标准型$NURBS$曲线。$N_{i,k}(t)$是由节点向量$T=[t_0,t_1,...,t_{n+k+1}]$上按照$deBoor-Cox$递推公式决定的k次规范B样条基函数。对于$NURBS$开曲线，通常将两端节点的重复度取为$k+1$,即$t_0=t_1=...=t_k,t_{n+1}=t_{n+2}=...=t_{n+k+1}$。

权因子：

权因子对$NURBS$曲线的影响是：

• 当保持控制顶点与其他权因子不变，减小或增加权因子时，起到把曲线推离和拉向相应顶点的作用。

曲面定义

$$p(u,v)=\frac{{\sum }_{i=0}^m{\sum }_{j=0}^n{N}_{i,p}(u)N_{j,q}(v){\omega }_{ij}{P}_{ij}}{{\sum }_{i=0}^m{\sum }_{j=0}^n{N}_{i,p}(u)N_{j,q}(v){\omega }_{ij}},(u,v)\in [0,1]\times [0,1]$$

式中，$P_{ij},i=0,1,...,m,j=0,1,...,n$构成控制网格。${\omega }_{ij}$是与定点$P_{ij}$联系的权因子，规定四角定点处用正权因子即${\omega }_{00},{\omega }_{m0},{\omega }_{0n},{\omega }_{mn}>0$,其余${\omega }_{ij}\ge 0$且顺序$p\times q$个权因子不同是为零。$N_{i,p}(u),i=0,1,...,m$和$N_{j,q}(v),j=0,1,...,n$分别为定义在节点矢量U和V上的非有理B样条基函数。节点矢量U=$[u_0,u_1,...,u_{m+p+1}]$与V=$[v_0,v_1,...,v_{n+q+1}]$按Hartley-Judd(哈特利-贾德)方法决定（第一个节点和最后一个节点必须是重复度为k+1。这就产生了所谓的clamped B-样条曲线，产生的曲线的开始和结尾连接在一起形成了一个闭环）。代码如下所示：

void HJ_calculate_knot(void)
{
	//两端点取K+1个重复度
	for(int i=0; i<=k; i++)
		knot[i]=0.0;                 
	for(int i=n+1; i<=n+k+1; i++)
		knot[i]=1.0;
	
	//计算n-k个内节点
	for(int i=k+1; i<=n; i++)
	{
		double sum=0.0;

		for(int j=k+1; j<=i; j++)
		{
			//计算分子项
			double numerator=0.0;
			for(int loop=j-k; loop<=j-1; loop++)
			{
				numerator+=sqrt((P[loop].x-P[loop-1].x)*(P[loop].x-P[loop-1].x)
		                      +(P[loop].y-P[loop-1].y)*(P[loop].y-P[loop-1].y));
			}

			//计算分母项
			double denominator=0.0;
			for(int loop1=k+1; loop1<=n+1; loop1++)
			{
				for(int loop2=loop1-k; loop2<=loop1-1; loop2++)
				{
					denominator+=sqrt((P[loop2].x-P[loop2-1].x)*(P[loop2].x-P[loop2-1].x)
		                           +(P[loop2].y-P[loop2-1].y)*(P[loop2].y-P[loop2-1].y));
				}
			}
			sum+=numerator/denominator;
		}
		knot[i]=sum;
	}
}