Courbe d'intersection

En géométrie, une courbe d'intersection est une courbe commune à deux objets géométriques. Dans le cas le plus simple, l'intersection de deux plans non parallèles dans l'espace euclidien est une droite. En général, une courbe d'intersection est constituée des points communs à deux surfaces qui se croisent transversalement, ce qui signifie qu'en tout point commun les normales à la surface ne sont pas colinéaires. Cette restriction exclut les cas où les surfaces se touchent ou ont des parties de surface communes.

Droite rouge comme intersection de deux plans, vert et bleu.

La détermination analytique de la courbe d'intersection de deux surfaces n'est aisée que dans des cas simples ; par exemple :

• l'intersection de deux plans
• la section plane d'une quadrique (sphère, cylindre, cône, etc.)
• l'intersection de deux quadriques dans des cas particuliers.

Pour le cas général, la littérature fournit des algorithmes de calcul des points d'intersection de deux surfaces.

Droite d'intersection de deux plans

Soient deux plans $\varepsilon _{i}:\quad {\vec {n}}_{i}\cdot {\vec {x}}=d_{i},\quad i=1,2,\quad {\vec {n}}_{1},{\vec {n}}_{2}$ linéairement indépendants, c'est-à-dire que les plans ne sont pas parallèles.

On recherche une représentation paramétrique ${\vec {x}}={\vec {p}}+t{\vec {r}}$ de la droite d'intersection.

Un vecteur directeur de la droite peut être obtenu à partir du produit vectoriel des vecteurs normaux : ${\vec {r}}={\vec {n}}_{1}\times {\vec {n}}_{2}$ .

Un point $P:{\vec {p}}$ de la droite d'intersection peut être déterminé en coupant les plans donnés $\varepsilon _{1},\varepsilon _{2}$ avec le plan$\varepsilon _{3}:{\vec {x}}=s_{1}{\vec {n}}_{1}+s_{2}{\vec {n}}_{2}$, qui est orthogonal à $\varepsilon _{1}$ et $\varepsilon _{2}$ . En injectant la représentation paramétrique de $\varepsilon _{3}$ dans les équations de $\varepsilon _{1}$ et $\varepsilon _{2}$ on obtient les paramètres $s_{1}$ et $s_{2}$ .

$P:{\vec {p}}={\frac {d_{1}({\vec {n}}_{2}\cdot {\vec {n}}_{2})-d_{2}({\vec {n}}_{1}\cdot {\vec {n}}_{2})}{({\vec {n}}_{1}\cdot {\vec {n}}_{1})({\vec {n}}_{2}\cdot {\vec {n}}_{2})-({\vec {n}}_{1}\cdot {\vec {n}}_{2})^{2}}}{\vec {n}}_{1}+{\frac {d_{2}({\vec {n}}_{1}\cdot {\vec {n}}_{1})-d_{1}({\vec {n}}_{1}\cdot {\vec {n}}_{2})}{({\vec {n}}_{1}\cdot {\vec {n}}_{1})({\vec {n}}_{2}\cdot {\vec {n}}_{2})-({\vec {n}}_{1}\cdot {\vec {n}}_{2})^{2}}}{\vec {n}}_{2}\ .$

Exemple

on cherche la droite ${\mathcal {D}}$, intersection de $\varepsilon _{1}:\ x+2y+z=1,\quad \varepsilon _{2}:\ 2x-3y+2z=2\ .$

Les vecteurs normaux sont ${\vec {n}}_{1}=(1,2,1)^{\top },\ {\vec {n}}_{2}=(2,-3,2)^{\top }$ et la direction de la droite d'intersection est donc ${\vec {r}}={\vec {n}}_{1}\times {\vec {n}}_{2}=(7,0,-7)^{\top }$ . Pour le point $P:{\vec {p}}$, on obtient de la formule ci-dessus ${\vec {p}}={\tfrac {1}{2}}(1,0,1)^{\top }\ .$ Ainsi

${\mathcal {D}}={\begin{cases}x&={\tfrac {1}{2}}+t\\y&=0\\z&={\tfrac {1}{2}}-t\end{cases}},\ t\in \mathbb {R}$

est une représentation paramétrique de la droite d'intersection.

Remarques :

1. Dans des cas particuliers, la détermination de la droite d'intersection par le pivot de Gauss peut être plus rapide.
2. Si l'un (ou les deux) des plans est donné par sa représentation paramétrique telle que ${\vec {x}}={\vec {p}}+s{\vec {v}}+t{\vec {w}}$, on obtient ${\vec {n}}={\vec {v}}\times {\vec {w}}$ comme vecteur normal et l'équation est : ${\vec {n}}\cdot {\vec {x}}={\vec {n}}\cdot {\vec {p}}$ .

Courbe d'intersection d'un plan et d'une quadrique

Dans tous les cas, la courbe d'intersection d'un plan et d'une quadrique (sphère, cylindre, cône...) est une conique. Une application importante des sections planes de quadriques est les lignes de contour de quadriques. Dans tous les cas (projection parallèle ou centrale), les courbes de niveau des quadriques sont des sections coniques.

Courbe d'intersection d'un cylindre ou d'un cône et d'une quadrique

Il est facile de déterminer les points d'intersection d'une droite avec une quadrique (c'est-à-dire une intersection droite-sphère) ; il suffit de résoudre une équation quadratique. Ainsi, toute courbe d'intersection d'un cône ou d'un cylindre (qui sont deux surfaces réglées) avec une quadrique est constituée des points d'intersection des droites et de la quadrique.

Les images montrent les possibilités qui se présentent lors de l'intersection d'un cylindre et d'une sphère :

1. la courbe a une seule boucle
2. la courbe a deux boucles non sécantes
3. la sphère et le cylindre se touchent en un point singulier. La courbe d'intersection est auto-sécante. La courbe est alors une fenêtre de Viviani
4. si le cylindre et la sphère ont le même rayon et que le milieu de la sphère est situé sur l'axe du cylindre, alors la courbe d'intersection se compose uniquement de points singuliers (formant un cercle). Cette définition peut servir à définir le contour apparent de la sphère.

• 
  Intersection d'une sphère et d'un cylindre à une boucle
• 
  Intersection d'une sphère et d'un cylindre à deux boucles

• 
  Intersection d'une sphère et d'un cylindre avec un point singulier
• 
  Intersection d'une sphère et d'un cylindre de même rayon

Cas général : méthode de marche

Courbe d'intersection : principe de l'algorithme de marche

En général, il n'y a pas de fonctionnalités particulières à exploiter. Une possibilité pour déterminer un polygone de points de la courbe d'intersection de deux surfaces est la méthode de marche. Il se compose de deux parties essentielles:

1. La première partie est l'algorithme de point de courbe, qui détermine à partir d'un point de départ au voisinage des deux surfaces un point sur la courbe d'intersection. L'algorithme dépend essentiellement de la représentation des surfaces données. La situation la plus simple est celle où les deux surfaces sont implicitement données par des équations $f_{1}(x,y,z)=0,\ f_{2}(x,y,z)=0$, car les fonctions fournissent des informations sur les distances aux surfaces et indiquent via les gradients le chemin vers les surfaces. Si une ou les deux surfaces sont données par leurs représentations paramétriques, les avantages du cas implicite n'existent pas. Dans ce cas, l'algorithme de point de courbe utilise des procédures chronophages comme la détermination du point de pied d'une perpendiculaire sur une surface.
2. La deuxième partie de la méthode de marche commence par un premier point sur la courbe d'intersection, détermine la direction de la courbe d'intersection à l'aide des normales de surface, puis fait un pas avec une longueur de pas donnée dans la direction de la ligne tangente, afin d'obtenir un point de départ pour un deuxième point de courbe...

L'algorithme de marche produit pour tout point de départ un polygone sur la courbe d'intersection. Si la courbe d'intersection se compose de deux parties, l'algorithme doit être exécuté en utilisant un deuxième point de départ pratique. L'algorithme est assez robuste. Habituellement, les points singuliers ne posent aucun problème, car la chance de rencontrer exactement un point singulier est très faible (voir image : intersection d'un cylindre et de la surface $x^{4}+y^{4}+z^{4}=1$ ).

• 
  Intersection de la surface d'équation $x^{4}+y^{4}+z^{4}=1$ avec un cylindre : cas de deux boucles disjointes
• 
  Intersection de la surface d'équation $x^{4}+y^{4}+z^{4}=1$ avec un cylindre : cas à une boucle
• 
  Intersection de la surface d'équation $x^{4}+y^{4}+z^{4}=1$ avec un cylindre : cas avec un point singulier

Application : ligne de contour

Un point $(x,y,z)$ de la courbe de niveau d'une surface implicite d'équation $f(x,y,z)=0$ et projection parallèle de direction ${\vec {v}}$ doit remplir la condition $g(x,y,z)=\nabla f(x,y,z)\cdot {\vec {v}}=0$, parce que ${\vec {v}}$ doit être un vecteur tangent, ce qui signifie que tout point de contour est un point de la courbe d'intersection des deux surfaces implicites

$f(x,y,z)=0,\ g(x,y,z)=0$ .

Pour les quadriques, $g$ est toujours une fonction linéaire. Ainsi la courbe de niveau d'une quadrique est toujours une section plane (c'est-à-dire une section conique).

La courbe de niveau de la surface $f(x,y,z)=x^{4}+y^{4}+z^{4}-1=0$ (voir photo) a été tracé par la méthode de la marche.

Remarque

La détermination d'un polygone de contour d'une surface paramétrique ${\vec {x}}={\vec {x}}(s,t)$ nécessite de tracer une courbe implicite dans le plan des paramètres.

Condition pour les points de contour : $g(s,t)=({\vec {x}}_{s}(s,t)\times {\vec {x}}_{t}(s,t))\cdot {\vec {v}}=0$ .

Courbe d'intersection de deux polyèdres

Courbe d'intersection entre polyèdres : trois maisons

Intersection de polyèdres : deux tores

La courbe d'intersection de deux polyèdres est un polygone (voir intersection de trois maisons). L'affichage d'une surface définie paramétriquement se fait généralement en cartographiant un réseau rectangulaire dans l'espace. Les quadrilatères spatiaux sont presque plats. Ainsi, pour l'intersection de deux surfaces définies paramétriquement, l'algorithme d'intersection de deux polyèdres peut être utilisé.

Articles connexes

• Théorie des intersections
• Exemples de courbes d'intersection
  • Conique
  • Section torique
  • Ovale de Cassini
  • Courbe de la crêpe
  • Fenêtre de Viviani

Références

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Intersection curve » (voir la liste des auteurs).

Bibliographie

• (en) C. L. Bajaj, C. M. Hoffmann et R.E. Lynch, « Tracing surface intersections », Computer Aided Geometric Design, vol. 5, n^o 4,‎ 1988, p. 285-307 (DOI 10.1016/0167-8396(88)90010-6).
• (en) R. E. Barnhill et S. N. Kersey, « A marching method for parametric surface/surface intersection », Computer Aided Geometric Design, vol. 7, n^os 1-4,‎ 1990, p. 257-280 (DOI 10.1016/0167-8396(90)90035-P).
• (en) R. Barnhill, G. Farin, M. Jordan et B. Piper, « Surface/Surface intersection », Computer Aided Geometric Design, vol. 4, n^os 1-2,‎ 1987, p. p 3-16 (DOI 10.1016/0167-8396(87)90020-3).
• (en) Changhe Tu, Wenping Wang et Jiaye Wang, « Classifying the nonsingular intersection curve of two quadric surfaces », Geometric Modeling and Processing. Theory and Applications. GMP 2002. Proceedings, Wako, Japan,‎ 2002, p. 23-32 (DOI 10.1109/GMAP.2002.1027493).

Liens externes

• COURBE ALGÉBRIQUE 3D sur MathCurve