0

我正在生成具有 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;
    }
}
}
4

1 回答 1

1

当 p != 0 但 KN[i+p] == KN[i] 时，在 Basis() 中会出现除以零错误，这就是为什么当结相等时它不起作用的原因，并且 KN 应该包含 14 个项目，而不是 11。

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

于 2014-04-28T01:19:58.540 回答