第四章 曲线和曲面.ppt

void bez_to_points(int n,int npoints,double P[],double points[]) // P为控制点坐标 points为采用几何作图算法生成的Bezier曲线上的离散点序列 离散点序列points的个数为npoints 控制点P的个数为n +1 { double t,delt; delt=1.0/(double)npoints;//将参数t npoints等分 t=0.0; for(int i=0;i=npoints;i++) { points[i]=decas(n, P, t); //分别求出npoints+1个离散点points的坐标 t+=delt; } } double decas(int n,double P[],double t) { int m,i; double *R, *Q, P0; R = new double[n +1]; Q = new double[n +1]; for(i=0;i=n;i++) R[i]= P [i]; //将控制点坐标P保存于R中 //需要作n次外部循环，方能产生最终Bezier曲线在点t的值 for(m=n;m0;m--) { //n次Bezier曲线在点t的值，可由两条n-1次Bezier曲线 //在点t的值通过线性组合而求得。 Q[i]= R [i]+t*( R [i+1]- R [i]); } for(i=0;i= m -1;i++) R[i]= Q [i]; } P0=R[0]; delete R; delete Q; return (P0); } 设给出四点的坐标是(1，1)，(2，3)，(4，3)，(3，1)，求所确定三次Bezier曲线在t=1/3时的值P(1/3)，算法的计算过程 Bezier几何作图算法计算过程 设控制点序列P0，P1，…，Pn确定的n次Bezier曲线是P(t)，用如下递归方式计算另一组点集： 如果令Pa(s)和Pb(s)分别是以控制点序列 和 确定的Bezier曲线，其中0≤s≤1，那么就有： 己知四点P0，P1，P2，P3，确定了一条三次Bezier曲线P(t)，可写出下式， 分裂法中的递归计算 分裂法的示意图 可通过计算验证Bezier曲线由前后两段构成。现以P0的系数为例，验证两端它的系数是相等的。 左端显然就是 。先看右端。若0≤t≤ ，这时就用前半段的表达式，观察分裂计算图注意到 中有1份P0， 中有 份， 中是 份， 中是 份，因此全部P0的系数是： 右端对 ≤t≤1，注意到仅 中有 份的P0，知P0的系数是： 设己知三次Bezier曲线P(t)的控制顶点是P0，P1，P2，P3，在P( )处将曲线分为两段，求出前半段的控制顶点Q0，Ql，Q2，Q3和后半段的控制顶点R0，R1，R2，R3，。有算法如下 void split_Bezier(Point P[]) { Point R[4],Q[4]; int i,j; for(i=0;i=3;i++) R[i]=P[i]; for(i=0;i=2;i++) { Q[i]=R[0]; for(j=0;j=2-i;j++) { R[j].x=(R[j].x+R[j+1].x)/2; //分别对相邻两控制点间的线段进行分裂 R[j].y=(R[j].y+R[j+1].y)/2; } } Q[3]=R[0]; } 分裂算法的计算 根据Bezier曲线的凸包性质，知道曲线上任意一点到线段P0P3的距离，小于P1和P2到线段P0P3距离中的较大者，即有： 任意事先给定的对画出曲线近似程度的要求ε0，可以取max(d(P1，P0P3)，d(P2，P0P3))ε为分裂停止的条件。 void new_split_Bezier(Point P[]) { Point R[4],Q[4]; int i,j; const double epsilon=0.01; if(maxdistance(P)epsilon) /*maxdistance(P) 为求max(d(P1,P0P3)，d(P2,P0P3))的函数*/ { MoveTo(P[0].x,P[0].y); LineTo(P[3].x,P[3