courbes de Bézier

Version imprimable

courbes de Bézier

Bonjour
Je suis sur la réalisation d'un logiciel de dessin vectoriel et je rencontre un pb avec les courbes de Bézier.
Le problème est je pense au niveau algo et pas directement lié au C, donc déplacez le topic si nécessaire...

Disons que ça marche a peu près mais c'est pas très précis.

voici ma fonction pour tracer une courbe:

formule pour un courbe de Bézier cubique (cf wikipédia):
B(t) = P0(1-t)^3 + 3*P1*t(1-t)^2 + 3*P2*t^2*(1-t) + P3*t^3
avec t compris entre 0 (début du traçé) et 1 (fin du traçé)
P0 = pt source, P1-P2 = pts de controles, P3 = pt cible

variables et fonctions:
this->aX, this->aY: position courante du point (donc P0 sur la formule)
ligne(floatx, float y); fonction pour tracer une segment du point courant (this->aX, this->aY) jusqu'à (x,y)
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 void Draw::bezierCubABS(float x1, float y1, float x2, float y2, float x, float y) { float ptx, pty; //point de destination de chaque segment du traçé float t; for(t=0; t<=1; t+=0.1) { ptx = this->aX*pow((1-t),3) + 3*x1*t*pow(1-t, 2) + 3*x2*pow(t, 2)*(1-t) + x*pow(t,3); pty = this->aY*pow((1-t),3) + 3*y1*t*pow(1-t, 2) + 3*y2*pow(t, 2)*(1-t) + y*pow(t,3); ligne(ptx, pty); } ligne(x,y); // finit le traçé au cas ou ça ne tombe pas pile poil }
La fonction me trace une espèce de patate qui ressemble à cercle, quand je diminue l'incrémentation de t ça ressemble plus à un cercle mais c'est du coup c'est segmenté.

Une idée?

17/12/2011, 19h10
natha31
résolu! xD
En ré-écrivant mon code pour ce topic je me suis aperçu d'une ânerie:
J'utilise mes this->aX, this->aY plusieurs fois alors qu'ils sont modifiés à chaque appel de ligne();
je l'ai ai donc stockés dans des variables provisoires:
Code:

1 2 float x0 = this->aX; float y0 = this->aY;

Avec l'algorithme de De Casteljau (C et bibliothèque SDL) :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
 
typedef struct _points{
        double x, y;
 }Points;
 
static void Casteljau(double t,Points* tab,int nbpoints,double *rx,double *ry)
{
	if (nbpoints == 1)
	{
		*rx = tab[0].x;
		*ry = tab[0].y;
		return;
	}
 
	Points subtab[nbpoints];
	int i;
	for(i=0;i<nbpoints-1;i++)
	{
		subtab[i].x = t*tab[i].x + (1-t)*tab[i+1].x;
		subtab[i].y = t*tab[i].y + (1-t)*tab[i+1].y;
	}
 
	Casteljau(t, subtab, nbpoints-1, rx, ry);
}
 
void ShowBezier(SDL_Surface* screen, Points* tab, double tstep, int nbPoints, Uint32 couleur)
{
	double t = 0.0;
	Points p, plast;
 
	plast.x=plast.y=0;
	double rx = 0.0;
	double ry = 0.0;
	int i = 0;
	while(t<=1.0)
	{
		Casteljau(t,tab,nbPoints,&rx,&ry);
		p.x = (int)rx;
		p.y = (int)ry;
		if (i!=0)
			ligne_Bresenham(screen,p,plast,couleur);
		plast.x = p.x;
		plast.y = p.y;
		t+=tstep;
		i++;
	}
}