Nurbs样条线算法推导和python实现.docxVIP

Nurbs样条线算法推导及python实现注：在现在的建模软件中Nurbs曲线还引入一个跨度的概念，在这里并不体现，这里是有几个控制点就是控制点数-1阶，最后好像到36阶曲线就发散了耶三个控制点：cbacbaa=(ax, ay)b=(bx, by)c=(cx, cy)0≤t≤1e为ab线上一点则：aee=(ex, ey)=(t*(bx-ax)+ax, t*(by-ay)+ay)同理，同一t值，f为bc线上一点为：f=(fx, fy)=(t*(cx-bx)+bx, t*(cy-by)+by)g点为动线ef上的动点为：g=(t*(fx-ex)+ex, t*(fy-ey)+ey)=(gx, gy)gx= t^2*(cx-2*bx+ax)+2*t*(bx-ax)+axgy= t^2*(cy-2*by+ay)+2*t*(by-ay)+ay动点g的运动方程入上图所示0:xx=ty=t四个控制点cdbacdba0≤t≤1a=(ax, ay) b=(bx, by) c=(cx, cy) d=(dx, dy) e=(ex, ey)=(t*(bx-ax)+ax, t*(by-ay)+ay) t?f = (fx , fy)=(t*(cx-bx)+bx, t*(cy-by)+by)t?g= (gx, gy)=(t*(dx-cx)+cx, t*(dy-cy)+cy)t?h=(hx, hy)=(t*(fx-ex)+ex, t*(fy-ey)+ey)tI = (ix , iy )=(t*(gx-fx)+fx, t*(gy-fy)+fy)tj = (jx, jy )=(t*(ix-hx)+hx, t*(iy-hy)+hy)ttN个控制点：def cv(N, t):n=[ Point(t*( N[i+1].x- N[i].x)+N[i].x, t*( N[i+1].y- N[i].y)+N[i].yfor i in range(len(N)-1) )]if len(n)==1: return n[0]else: return cv(n, t)0≤t≤1N=[p0, p1, p2…pn]print cv(N, t)python代码实现：（按分段数计算样条线上的点）def cv(N, seg):step=1/segsP=N[0]eP=N[-1]p=[sP]oldN=Nfor i in range(step, 1, step): t=i*stepcount= len(N)while count1:n=[ Point(t*( N[i+1].x- N[i].x)+N[i].x , t*( N[i+1].y- N[i].y)+N[i].y) for i in range(count-1) ] N=ncount-=1 N=oldN p.append(n[0])p.append(eP)return p用多项式的方法优化运算：各项常数确定：Cx, Cy因为y因为y和x的运算过程相同的故下面只求xxxxx?t?t?t?t??ttt?tt?t4C4?6?4?C?DDDD?EEE?FF?G?FEEDDDCCCCGFEDC0101234xxxx11,2,11,3,3,11,4,6,4,1各项系数增量scale:scalescaledef getScales(n): 各项常系数的缩放值 1,1 n=1 1,2,1 n=2 1,3,3,1 n=3 1,4,6,4,1 n=4 ... ... 返回第n行数据 n为控制点数目 if n2: return None n=n-2 s1=[1,1] scales=[1,1] for i in range(n): scales=[s1[j]+s1[j+1] for j in range(len(s1)-1)] scales.insert(0,1) scales.append(1) s1=scales return scales各项常系数C Cdef getCs(scales,points): 各项常系数 x_4 n=0 x_3,F_4 n=1 x_2,E_3,E_4 n=2 x_1,D_2,D_3,D_4 n=3 x_0,C_1,C_2,C_3,C_4