Question

我正在生成具有10个控制点的3阶B样条曲线。我的结矢量有14个条目，范围从0到1.由于某种原因，我的曲线始终从原点开始和结束。如果任何结等于另一个结，它也无法绘制曲线。

这是我的基本功能代码：

if( p == 0 ) {
    if( (KN[i] <= u) && (u < KN[i+1]) ){
        return 1;
    } 
    return 0;
}

else{
    u1 = (u - KN[i])/(KN[i+p]-KN[i]);
    u2 = (KN[i+p+1]-u)/(KN[i+p+1]-KN[i+1]);

    return u1*Basis(i,(p-1),u)+u2*Basis(i+1,p-1,u);
}

其中i是当前点（0到9），p是度，3，u是我的域[0,1]中的位置。我的曲线制作函数调用此基函数，如下所示：

void makeCurve (int n, int p){
    G_rgb(0,1,0);

    double u, x, y;
    int i;

    for(u = 0; u <= 1; u += 0.0001){
        x = 0;
        y = 0;

        for(i = 0; i < n; i ++){
            x += Basis(i,p,u) * PX[i];
            y += Basis(i,p,u) * PY[i];
         }

         G_point(x,y);
    }
}

其中n是点数，10，p是度，3。

此处参考该计划的全部内容：

double PX[20], PY[20], KN[11];

int factorial( int n )
{
int res, i;
    if( n <= 0 ){
        return 1 ; // safeguard 0 and -ve
    }

    res = n ;

    for(i = n-1; i > 1; i --){
        res *= i;
    }

    return res;
}

double Basis(int i, int p, double u){
double u1, u2;

if( p == 0 ) {
    if( (KN[i] <= u) && (u < KN[i+1]) ){
        return 1;
    } 
    return 0;
}
else{
    u1 = (u - KN[i])/(KN[i+p]-KN[i]);
    u2 = (KN[i+p+1]-u)/(KN[i+p+1]-KN[i+1]);

    return u1*Basis(i,(p-1),u)+u2*Basis(i+1,p-1,u);
}
}

void makeCurve (int n, int p){

G_rgb(0,1,0);

double u, x, y;
int i;

for(u = 0; u <= 1; u += 0.0001){
    x = 0;
    y = 0;

    for(i = 0; i < n; i ++){
        x += Basis(i,p,u) * PX[i];
        y += Basis(i,p,u) * PY[i];
    }

    G_point(x,y);
}
}


void makeLines(int n){
    int i;

    G_rgb(1,1,1);

for(i = 0; i < n - 1; i ++){
    G_line(PX[i],PY[i],PX[i+1],PY[i+1]);
}
}

void makeDots(int n){

int i;

G_rgb(1,0,0);

for(i = 0; i < n; i ++){
    G_fill_circle(PX[i],PY[i],5);
}

}

void main(){

G_init_graphics(800,800);
G_rgb(0,0,0);
G_clear();

double p[2], u;

int i, n, m;
n = 0;

for(i = 0; i < 10; i ++){

    G_wait_click(p);

    PX[i] = p[0];
    PY[i] = p[1];

    G_rgb(1,0,0);
    G_fill_circle(p[0],p[1],5);

    if(i != 0){
        G_rgb(1,1,1);
        G_line(p[0],p[1],PX[i-1],PY[i-1]);
    }
}

// so n = 9 and p = 6 so m = 13 and there are 14 knots
KN[0] = 0;
KN[1] = 0.03;
KN[2] = 0.1;
KN[3] = 0.12;
KN[4] = 0.14;
KN[5] = 0.28;
KN[6] = 0.42;
KN[7] = 0.57;
KN[8] = 0.71;
KN[9] = 0.85;
KN[10] = 0.9;
KN[11] = 0.99;
KN[12] = 0.999;
KN[13] = 1;

makeCurve(10, 3);

while(1 == 1){
    if(G_wait_key() == 'q') {
        break;
    }
}
}

Answer 1

当p！= 0但是KN [i + p] == KN [i]时，你在Basis（）中得到除零误差，这就是为什么当结节相等时它不起作用的原因，以及KN应该包含14个项目，而不是11个。

gcc和clang都会自动警告你KN的大小错误，所以你应该考虑更新你的编译器。

B样条曲线结束并从原点开始？

1 个答案: