Algorithme de de casteljau

Algorithme de de casteljau: Algorithme de De Casteljau

L'algorithme de De Casteljau est un algorithme récursif trouvé par Paul de Casteljau pour approximer efficacement les polynômes écrit dans la base de Bernstein.

Cet algorithme peut être utilisé pour dessiner des courbes et des surfaces de Bézier. L'idée principale dans ce cas repose sur le fait qu'une restriction d'une courbe de Bézier est aussi une courbe de Bézier. L'algorithme calcule de manière efficace le point de parametres $t = T$ et les points de contrôle des courbes de la restriction à $t \leqslant T$ et à $t \geqslant T$ . On applique alors de nouveau l'algorithme sur les deux restrictions jusqu'à réaliser un critère donné (celui-ci peut être par exemple que la précision soit inférieure au pixel).

Cet algorithme semble ne plus être le plus efficace car il ne permettrait pas d'utiliser l'antialiasing étant donné qu'il travaille pixel par pixel et ne donne pas d'information sur la tangente.

Sommaire

1 Historique

2 L'algorithme de calcul d'un point

2.1 Principe

2.2 Algorithme

2.3 Eléments de preuve de l'algorithme

2.4 Complexité

3 L'algorithme dans le cas de courbes de Bézier

3.1 Principe

3.2 Algorithme

3.3 Complexité

4 L'algorithme dans le cas de surfaces de Bézier

5 Voir aussi

Historique

Historiquement, c'est avec cet algorithme que les travaux de M. De Casteljau commençaient en 1959 chez Citroën. Ils étaient publiés comme des rapports techniques, tenus très au secret par Citroën.

Ces travaux restèrent inconnus jusqu'en 1975 quand W. Böhm en a pris connaissance et les a rendu public. Cet algorithme a été très utile pour l'informatique qui utilise les courbes de bézier dans de nombreux cas (logiciels de dessin, de modélisation, ...), et sans lequel le développement de l'utilisation des courbes de Pierre Bezier n'aurait pas pu se faire.

L'algorithme de calcul d'un point

Principe

Considérons une courbe de Bézier définie par les points de contrôles $P_0,\dots,P_N$ , où les $P i$ sont des points de $\R^m, m \geqslant 2$ . Ici on souhaite simplement calculer le point de paramètre $t \in [0,1]$ .

Principe de l'algorithme de De Casteljau

Comme on peut le voir sur l'image, en calculant les barycentres de paramètres ${t,1 - t}$ des points de contrôle consécutifs de la courbe, puis les barycentres de même paramètres de ces barycentres et ainsi de suite itérativement, on définit de cette manière une suite de listes de points que l'on va indexer $P^j_i$ , où $P^j_i$ est le barycentre de $\{P^{j-1}_i,P^{j-1}_{i+1}\}$ . La dernière liste ne contient en fait qu'un point, qui est le point de la courbe de paramètre $t$ .

Algorithme

En pseudo-code, ceci donne:

// Calcul des points intermédiaires Pour j de 1 à N faire | | Pour i de 0 à N-j faire | | | | T[i][j] = t*T[i+1][j-1] + (1-t)*T[i][j-1] | | | Afficher T[0][N] // Afficher (ou stocker) le point

Eléments de preuve de l'algorithme

Pour démontrer l'algorithme, il faut prouver par récurrence limitée que

$\forall i \in [0,n], \forall t \in [0,1], B(t) = \sum_{k=0}^{n} P^{0}_{k} B^{n}_{k}(t) = \sum_{k=0}^{n-i} P^{i}_{k} B^{n-i}_{k}(t)$

Et d'appliquer la formule en $i = n$ , ce qui nous donne le résultat directement.

La récurrence se démontre facilement en utilisant la propriété de construction des points $P_i^j = (1-t)P_i^{j-1} + tP_{i+1}^{j-1}$ pour séparer la somme en deux, puis en faisant un changement d'indice et en utilisant une propriété des polynômes de Bernstein $B_i^n(t) = (1-t)B_i^{n-1}(t) + t B_{i+1}^{n-1}(t)$ pour regrouper les deux sommes. Pour réintégrer les cas particuliers, il faut utiliser deux autres propriétés des polynômes de Bernstein: $B^n_n (t) = t^n = t B^{n-1}_{n-1}$ et $B^n_0 (t) = (1-t)^n = (1-t) B^{n-1}_{0}$

Complexité

L'éxécution d'une étape de l'algorithme est quadratique en nombre de points de contrôle de la courbe (la boucle imbriquée donne lieu à $\mathcal{O}(N^2)$ opérations).

L'algorithme dans le cas de courbes de Bézier

Considérons encore une courbe de Bézier définie par les points de contrôles $P_0,\dots,P_N$ , où les $P i$ sont des points de $\R^m, m>=2$ .

Principe

On va ici appliquer l'algorithme précédent pour trouver le point de paramètre $1 / 2$ , et conserver les barycentres intermédiaires. En effet, la courbe de Bézier de points de contrôle $P^i_0,i\in [0,N]$ est égale à la restriction de la courbe originale à $t\leqslant 1/2$ , et la courbe de Bézier de points de contrôle $P^{N-i}_i,i\in [0,N]$ est égale à la restriction de la courbe originale à $t\geqslant 1/2$ .

On affiche ou mémorise le point de paramètre $1 / 2$ (qui est $P^N_0$ ) et applique récursivement l'algorithme sur les deux courbes (de même nombre de points de contrôle que l'originale). On s'arrête quand un critère à choisir est vérifié (typiquement la distance entre les points inférieure à une limite).

Plutot que le paramètre $1 / 2$ , on pourrait prendre un paramètre quelconque et l'algorithme fonctionnerait encore, mais le paramètre $1 / 2$ est celui qui converge en moyenne le plus rapidement.

Algorithme

Initialisation: affecter le tableau des points de contrôle dans les $P^0_i$

Voir que le critère d'arrêt n'est pas vérifié: il y a plusieurs possibilités dont:

Si on fait tous les calculs avec des nombres entiers, un choix peut être de s'arrêter lorsque tous les points sont confondus. (c'est-à-dire que la courbe n'est représentée que par un seul pixel)

Si le calcul n'est pas fait sur des nombres entiers on peut s'arrêter quand les points sont distants d'une distance inférieure à une valeur choisie

On effectue M itérations puis on relie les points obtenus (c'est-à-dire que l'on trace le polygone de Bézier), M étant déterminé empiriquement

Calcul les points intermédiaires de l'algorithme : il y a deux méthodes possibles qui donnent finalement les mêmes points

Méthode constructive (en construisant une suite de milieux): On définit itérativement les points $P^1_0,\dots,P^1_{N-1}, P^2_0,\dots,P^2_{N-2},\dots,P^{N-1}_0, P^{N-1}_1, P^N_0$ tels que $P^j_i$ soit le milieu de $[P^{j-1}_iP^{j-1}_{i+1}]$

Méthode matricielle (en construisant directement les points comme barycentre, les coefficients étant donnés par les matrices de De Casteljau) :

$\begin{pmatrix} P_0^0 \\ \vdots \\ P_0^N \end{pmatrix} = D_0^N * \begin{pmatrix} P_0^0 \\ \vdots \\ P_N^0 \end{pmatrix}$ et $\begin{pmatrix} P_0^N \\ \vdots \\ P_N^0 \end{pmatrix} = D_1^N * \begin{pmatrix} P_0^0 \\ \vdots \\ P_N^0 \end{pmatrix}$

Mémorisation: on mémorise le point $P^N_0$ qui est sur la courbe

Appel récursif: on appelle l'algorithme sur les deux courbes de Bézier intermédiaires définies par les points $P^i_0,i\in [0,n]$ d'une part, $P^{n-i}_i,i\in [0,n]$ d'autre part.

Voici un exemple d'implémentation de l'algorithme en pseudo-code avec pour critère d'arrêt l'égalité des points (on travaille donc sur des entiers) et la construction constructives pour calculer les points intermédiaires:

Entrée : tableau T[0][0...N] des coordonnées des points de contrôle.

Si T[0][0] = T[0][1] = ... = T[0][N] //si le critère d'arrêt est vérifié on s'arrête alors | | fin | Sinon //pour dessiner | | // Calcul des points intermédiaires | Pour i de 1 à N faire | | | | Pour j de 0 à N-i faire | | | | | | T[i][j] = milieu de T[i-1][j] T[i-1][j+1] | | Afficher T[N][0] // Afficher (ou stocker) le point milieu | | // Construction des courbes restreintes | Pour i de 0 à N faire | | | | T'[i] = T[i][0] | | T"[i] = T[N-i][i] | | // Appel récursif | de_Casteljau(T') | de_Casteljau(T")

Complexité

Si on prend comme critère d'arrêt un nombre d'appels récursifs constant, comme le nombre de points de contrôle est constant pendant les appels récursifs reste constant et qu'à chaque étape de récursion on double le nombre de courbes étudiées, la complexité de l'algorithme est en $\mathcal{O}(N^2 * 2^M)$ avec $M$ le nombre de récursions (mais il est linéaire en nombre de points calculés car pour $M$ itérations il y a $2 M$ points calculés).

L'algorithme dans le cas de surfaces de Bézier

Une surface de Bézier est définie par une double somme de points: $P(t_1,t_2) = \sum_{i=1}^m \sum_{j=1}^n P_{i,j} B_m^i (t_1) B_n^j (t_2)$

Le principe de l'algorithme est de réécrire la formule sous la forme:

$P(t_1,t_2) = \sum_{i=1}^m B_m^i (t_1) \left( \sum_{j=1}^n P_{i,j} B_n^j (t_2) \right)$

et en renommant $Q_{i} (t_2) = \sum_{j=1}^n P_{i,j} B_n^j (t_2)$ , on obtient

$P(t_1,t_2) = \sum_{i=1}^m Q_{i} (t_2) B_m^i (t_1)$

En remarquant que les $Q i (t 2)$ sont des points de courbes de Bézier, le principe de l'algortithme arrive. A chaque itération de l'algorithme

calculer les points $Q i (1 / 2)$ et les restrictions des courbes de Bézier pour chaque $i$ par l'algorithme de De Casteljau sur les courbes

calculer puis afficher/stocker les points de la courbe de Bézier de points de contrôle les $Q i (1 / 2)$ par l'algorithme de De Casteljau sur les courbes

appliquer récursivement l'algorithme sur les deux surfaces obtenues en regroupant les restrictions pour $t_2 \leqslant 1/2$ et $t_2 \geqslant 1/2$

Voir aussi

Courbe de Bézier en C qui utilise cet algorithme.

Ce document provient de « Algorithme de De Casteljau ».

Catégorie : Algorithme

Contenu soumis à la licence CC-BY-SA. Source : Article Algorithme de de casteljau de Wikipédia en français (auteurs)

Игры ⚽ Нужно сделать НИР?

Regardez d'autres dictionnaires:

Algorithme De De Casteljau — L algorithme de De Casteljau est un algorithme récursif trouvé par Paul de Casteljau pour approximer efficacement les polynômes écrit dans la base de Bernstein. Cet algorithme peut être utilisé pour dessiner des courbes et des surfaces de Bézier … Wikipédia en Français
Algorithme de de Casteljau — L algorithme de De Casteljau est un algorithme récursif trouvé par Paul de Casteljau pour approximer efficacement les polynômes écrit dans la base de Bernstein. Cet algorithme peut être utilisé pour dessiner des courbes et des surfaces de Bézier … Wikipédia en Français
Algorithme de De Casteljau — L algorithme de De Casteljau est un algorithme récursif trouvé par Paul de Casteljau pour approximer efficacement les polynômes écrit dans la base de Bernstein. Cet algorithme peut être utilisé pour dessiner des courbes et des surfaces de Bézier … Wikipédia en Français
Matrice De De Casteljau — Les matrices de de Casteljau sont des matrices de Markov triangulaires (ou leurs transposées suivant les conventions) principalement utilisées dans l algorithme de De Casteljau. Pour une taille N fixée, il y deux matrices D0 et D1 définies par … Wikipédia en Français
Matrice de de Casteljau — Les matrices de de Casteljau sont des matrices de Markov triangulaires (ou leurs transposées suivant les conventions) principalement utilisées dans l algorithme de De Casteljau. Pour une taille N fixée, il y deux matrices D0 et D1 définies par … Wikipédia en Français
Matrice de de casteljau — Les matrices de de Casteljau sont des matrices de Markov triangulaires (ou leurs transposées suivant les conventions) principalement utilisées dans l algorithme de De Casteljau. Pour une taille N fixée, il y deux matrices D0 et D1 définies par … Wikipédia en Français
Matrice de De Casteljau — Les matrices de de Casteljau sont des matrices de Markov triangulaires (ou leurs transposées suivant les conventions) principalement utilisées dans l algorithme de De Casteljau. Pour une taille N fixée, il y deux matrices D0 et D1 définies par où … Wikipédia en Français
De casteljau — Paul de Faget de Casteljau Paul de Faget de Casteljau (né en 1930 à Besançon) est un mathématicien et physicien français. Il a fait ses études à l École normale supérieure de Paris. Il est connu pour sa découverte des formes à pôles en 1959 et l… … Wikipédia en Français
Paul Faget de Casteljau — Paul de Faget de Casteljau Paul de Faget de Casteljau (né en 1930 à Besançon) est un mathématicien et physicien français. Il a fait ses études à l École normale supérieure de Paris. Il est connu pour sa découverte des formes à pôles en 1959 et l… … Wikipédia en Français
Paul de Casteljau — Paul de Faget de Casteljau Paul de Faget de Casteljau (né en 1930 à Besançon) est un mathématicien et physicien français. Il a fait ses études à l École normale supérieure de Paris. Il est connu pour sa découverte des formes à pôles en 1959 et l… … Wikipédia en Français

Mark and share
Search through all dictionaries
Translate…
Search Internet

Share the article and excerpts