Introduction à l'algorithme du lancer de rayons

2.a. Principe optique

La formation d'une image par un système optique (œil, caméra, etc) obéit aux lois de la perspective, expliquées dans Projection perspective. Le système est constitué d'un centre optique C et d'un plan de projection (ou plan image) Π.

La lumière provenant d'un point M d'un objet suit le rayon MC pour parvenir au plan image.

Figure pleine page

La lumière qui part de M pour parvenir au plan image provient des sources de lumière. Elle peut en provenir directement, ou par l'intermédiaire de la diffusion et de la réflexion effectuée par les autres objets.

Le lancer de rayons consiste à considérer le trajet de la lumière dans le sens inverse du sens réel, en partant du centre optique C vers le point M. L'éclairement (et la couleur) du point image P (intersection de ce rayon avec le plan image) est égal au flux lumineux le long du rayon MP. Ce flux lumineux provient de

• la lumière provenant directement des sources,
• la lumière diffusée par d'autres objets.

On fait une distinction entre la réflexion spéculaire, qui suit la loi de Descartes de la réflexion, et la réflexion diffuse. Cette dernière est de loin la plus difficile à simuler, en raison du grand nombre de rayons qu'elle implique (en principe une infinité). Dans cette introduction, nous nous limitons à la première cause d'éclairement. Pour obtenir un rendu réaliste, il faut tenir compte de la lumière diffusée, qui est un phénomène très important dans la réalité. On suppose de plus que les sources de lumière sont ponctuelles.

Un autre phénomène à prendre en compte est la réflexion et la transmission éventuelles de la lumière au point M. Nous nous limiterons à la réflexion, mais la transmission se traite de manière analogue.

Le point de départ de l'algorithme est le lancer d'un rayon dans la direction CP, où P est le point de l'image dont on cherche l'éclairement. On commence par déterminer le point M d'intersection de ce rayon avec le premier objet rencontré. C'est ce point qui apparaît sur l'image en P.

Pour calculer le flux lumineux provenant directement des sources, on lance des rayons du point M vers les sources, qui sont ici ponctuelles. Ces rayons sont appelés rayons d'éclairement (représentés en bleu sur la figure). La figure montre deux rayons d'éclairement : le premier, lancé vers la source 1, l'atteint effectivement. Ce rayon contribue donc à l'éclairement de la surface au point M. Le second, lancé vers la source 2, ne l'atteint pas à cause de la présence d'un objet sur son chemin. La surface de l'objet au point M peut être plus ou moins réfléchissante. Pour prendre cela en compte, on lance à partir du point M un rayon réfléchi. Si ce rayon réfléchi rencontre un objet en un point M', la lumière provenant de M' parvient au point P par la réfléxion en M : l'éclairement du point M' contribue donc à l'éclairement du point P sur le plan image.

On voit bien la récursivité de l'algorithme : le flux lumineux qui quitte M' en direction de M se calcule de la même manière, en lançant des rayons vers les sources, et en lançant un rayon réfléchi si la surface est réfléchissante au point M'.

Dans certains cas, le tracé récursif des rayons réfléchis peut conduire à une boucle sans fin, ou du moins limitée par le nombre de récursions autorisé. Il faut donc limiter volontairement le nombre de récursions, ce qui revient à limiter le nombre de réflexions multiples.

2.b. Modèle de surface

Voyons comment une source de lumière ponctuelle contribue au flux lumineux. La figure suivante montre le rayon incident, le rayon envoyé vers la source et la normale à la surface au point M.

Figure pleine page

Il faut bien noter que le sens réel de la lumière est inverse de celui représenté par les flèches, qui correspond aux sens des rayons calculés. On cherche le flux lumineux envoyé à partir du point M dans la direction du rayon incident. En général, ce flux dépend des angles entre les deux rayons et de la normale n à la surface au point M.

Le modèle le plus simple pour représenter les propriétés optiques d'une surface est le modèle de Lambert. Il consiste à attribuer le facteur cos(θ) pour l'éclairement du point M, et à supposer que le flux lumineux envoyé dans la direction du rayon incident est indépendant de l'angle qu'il fait avec la normale. On introduit aussi un coefficient de diffusion pour la surface, et le flux s'écrit donc :

$$\Phi =k_{d}Icos\theta \text{(1)}$$

où I est l'intensité de la source.

Les grandeurs Φ, I et k_d dépendent en fait de la longueur d'onde λ. La couleur de l'objet dépend de la fonction k_d(λ). La source elle-même peut déliver une intensité I(λ) fonction de la longueur d'onde. Une manière simple de traiter l'influence de la longueur d'onde est de considérer trois composantes RVB (rouge,vert, bleu) pour ces grandeurs. Par exemple, la surface d'un objet rouge aura une composante rouge de son coefficient de diffusion plus grande que les composantes bleue et verte. Ce mode de calcul permet d'obtenir directement les composantes RVB du point P sur l'image finale.

Dans la suite de cette introduction, on se contente de calculer des flux monochromes, sans tenir compte de l'influence de la longueur d'onde.

En ce qui concerne le rayon réfléchi, on doit définir un coefficient de réflexion. Physiquement, ce coefficient de réflexion dépend de l'angle entre le rayon incident et la normale. Pour simplifier, on considèrera un coefficient de réflexion indépendant de l'angle.

2.c. Algorithme

L'algorithme est principalement constitué d'une fonction lancer_rayon(rayon) qui cherche le point de la scène rencontré par un rayon et calcule le flux lumineux le long de ce rayon. La fonction s'appelle elle-même récursivement dans le cas où le point rencontré est sur une surface réfléchissante.

Voici la séquence des opérations effectuées par cette fonction :

• Calcul des intersections avec les tous les objets de la scène. Lorsqu'il y a une intersection, on stocke la distance t (analogue à un temps de parcours) entre le point de départ et le point d'intersection.
• Si aucun objet n'est rencontré, renvoi de l'éclairement du fond.
• L'intersection la plus proche est retenue, celle dont la valeur de t est la plus petite.
• Le flux est initialisé à zéro.
• Pour chaque source ponctuelle, envoi d'un rayon vers la source. Les intersections de ce rayon avec les objets de la scène sont recherchées. Si aucune intersection n'est trouvée, la contribution de cette source au flux lumineux est ajoutée, en utilisant un modèle de surface.
• Si la surface est réfléchissante, calcul du rayon réfléchi en utilisant la loi de Descarte de la refléxion spéculaire. La fonction lancer_rayon(rayon) est appelée (appel récursif) et le résultat est ajouté au flux, avec un facteur de multiplication (coefficient de réflexion).
• Renvoi du flux.

Pour effectuer la synthèse d'une image complète, il faut échantillonner le plan image, par exemple avec un quadrillage régulier, puis appeler la fonction lancer_rayon(rayon) pour chaque rayon CP et attribuer l'éclairement obtenu au point P du plan image. Une manière simple de le faire est de lancer un rayon pour chaque pixel de l'image finale.

2.d. Exemple de calcul d'intersection

Le calcul de l'intersection d'un rayon avec un objet de la scène dépend évidemment de la forme de l'objet. Nous allons considérer le cas simple d'un objet sphérique.

Pour chaque rayon lancé, il faut rechercher l'intersection avec les N objets de la scène. Pour les scènes très complexes, cela fait beaucoup de calculs, d'autant plus qu'il y a au moins autant de rayons que de pixels dans l'image. On est alors amené à utiliser une méthode d'accélération, qui consiste à éliminer le plus rapidement possible les objets qui ne sont pas coupés par le rayon. Nous n'aborderons pas ce sujet dans cette introduction. On note seulement qu'une méthode d'accélération (parmis d'autres) consiste à définir pour chaque objet une sphère englobante. On commence par rechercher l'intersection du rayon avec cette sphère. Si cette intersection existe, on peut aborder le calcul de l'intersection avec l'objet, qui peut être beaucoup plus complexe.

Un rayon est défini par un point initial P et un vecteur unitaire U. Une sphère est définie par son centre C et son rayon r. Voyons comment déterminer l'intersection éventuelle du rayon avec la sphère. Une première idée est de considérer l'équation paramétrique du rayon et de rechercher le paramètre t pour l'intersection avec la sphère : cela donne une équation du second degré pour t. L'inconvénient de cette méthode est qu'elle nécessite beaucoup de calculs (en particulier le calcul d'une racine carrée), même lorsque le rayon ne coupe pas la sphère, ce qui est finalement le cas le plus fréquent.

Une méthode géométrique ([1]) permet d'éliminer rapidement le cas du rayon ne rencontrant pas la sphère, et de calculer le paramètre t seulement lorsque cela est nécessaire.

Figure pleine page

On commence par comparer la distance PC avec le rayon r pour savoir si le point de départ du rayon se trouve dans ou hors la sphère. Soit H le projeté orthogonal de C sur le rayon. La distance algébrique de P à H est

$$t_{PH}=\overrightarrow{PC}\cdot \overrightarrow{U}\text{(2)}$$

Si cette distance est négative et si P est en dehors de la sphère, alors le rayon ne coupe pas la sphère.

Le carré de la distance entre le centre de la sphère et le point H est :

$$D^2=C{H}^2=P{C}^2-t_{PH}^2\text{(3)}$$

Le carré de la distance MH est :

$$t_{MH}^2=M{H}^2=r^2-D^2\text{(4)}$$

Si ce carré est négatif, c'est que le point H est hors de la sphère et donc que le rayon ne coupe pas la sphère. On peut finalement calculer la distance t=PM. Si le point P est hors de la sphère, on a :

$$t=t_{PH}-\sqrt{t_{MH}^2}\text{(5)}$$

Si le point P est dans la sphère, on a :

$$t=t_{PH}+\sqrt{t_{MH}^2}\text{(6)}$$

Si la valeur de t est nulle, c'est que le point initial du rayon se trouve déjà sur la surface. Dans ce cas, on considère qu'il n'y a pas d'intersection. Ce cas doit en principe se produire lorsqu'on lance un nouveau rayon à partir de M (rayon d'éclairement ou rayon réfléchi). Cependant, des erreurs d'arrondis dans les calculs peuvent placer le point M légèrement à l'intérieur de la sphère, ce qui conduira à l'obtention d'une intersection avec la sphère, qui n'a physiquement pas lieu puisque le point se trouve en fait sur la surface. Pour éviter ce problème, on modifie légèrement la valeur de t afin que le point M soit de manière certaine hors de la sphère, par exemple en multipliant t par 0.999.