Voir la source de Ombrage de Phong

Le terme '''''ombrage de Phong''''' désigne à la fois le '''modèle d'illumination de Phong''' et l{{'}}'''interpolation de Phong''', deux [[algorithme]]s de traitement 3D en [[infographie]]. Tous les deux furent développés par [[Bùi Tường Phong]] et publiés en 1973.

== Modèle d'illumination de Phong ==
=== Description du modèle ===
L'illumination de Phong est un modèle local, c'est-à-dire que le calcul se fait en chaque point. Ce [[Loi empirique|modèle empirique]] n'a rien d'exact, mais permet de calculer de manière crédible la lumière réfléchie par le point étudié, pour cela il combine trois éléments : la lumière ambiante, la lumière diffuse ([[Loi de Lambert|modèle lambertien]]) et la [[lumière spéculaire]] (voir [[Réflexion (optique)#Réflexion spéculaire|réflexion spéculaire]]).

Le but est de calculer l'intensité lumineuse qui va être émise par réflexion par le point étudié, éclairé par une source supposée ponctuelle, dans une direction précise (celle de l'observateur).
Pour cela, la lumière est séparée en trois composantes :
# La ''composante ambiante'' représente les parasites provenant d'autre chose que la source considérée, la lumière réfléchie par d'autres points par exemple. La lumière ambiante est supposée égale en tout point de l'espace ;
# La lumière incidente est réfléchie dans toutes les directions. La ''composante diffuse'' indique l'intensité qui repart en tenant compte de l'inclinaison avec laquelle la lumière incidente arrive sur la surface, mais en supposant que l'intensité est la même quelle que soit la direction que prend le rayon réfléchi ;
# Malgré tout, il y a plus de lumière renvoyée dans la direction de la réflexion géométrique (celle dans laquelle repartirait le rayon en arrivant sur un miroir). Le rôle de la ''composante spéculaire'' est d'en tenir compte.

[[Image:Phong components version 4.png|655px|center|thumb|Illustration du modèle de Phong : l'intensité totale est la somme des trois composantes décrites ci-dessus.]]

=== Paramètres ===

[[Image:Reflexion_a_deux_directions.svg|right|frame|Définition des différents vecteurs.]]
* On définit pour chaque matériau des constantes caractéristiques :
** <math>k_a \in [0,1]</math> : constante liée à la composante ambiante, la proportion de lumière renvoyée ;
** <math>k_d \in [0,1]</math> : constante liée à la composante diffuse ;
** <math>k_s \in [0,1]</math> : constante liée à la composante spéculaire ;
** <math>\alpha \gg 1</math> : constante liée au ''[[Brillant et mat|brillant]]'' du matériau : plus <math>\alpha</math> est grand, plus la surface est brillante. Cette constante peut prendre des valeurs ''élevées'' : 10, 100 ou plus.
(Les valeurs indiquées sont seulement un ordre de grandeur.)

* On appelle <math>i_a</math>, <math>i_d</math> et <math>i_s</math> l'intensité des lumières incidentes ambiante, diffuse et spéculaire. <math>I_a</math>, <math>I_d</math>, <math>I_s</math> et <math>I</math> les intensités réfléchies, <math>I</math> étant le total.

* On définit les vecteurs suivants : <math>\vec L</math> pour la lumière, <math>\vec N</math> pour la normale à la surface, <math>\vec V</math> pour la direction de vue de l'observateur et <math>\vec R</math> pour la direction dans laquelle serait réfléchie la lumière sur un miroir.<br /><math>\vec R</math> se déduit par la relation <math>\vec R = 2(\vec N \cdot \vec L)\vec N - \vec L = 2\cos \theta \vec N - \vec L</math>.
Tous ces vecteurs doivent être normalisés pour que les produits scalaires donnent simplement le cosinus de l'angle entre les vecteurs.

=== Formules ===
La composante ambiante, <math>I_a</math> est simplement donnée par :

:<math>I_a = i_a.k_a</math>

La composante diffuse, <math>I_d</math> est donnée par :

:<math>I_d = i_d k_d (\vec L \cdot \vec N) = i_d k_d \cos \theta</math>
:De la même manière que le Soleil chauffe plus au zénith, <math>I_d</math> est maximale lorsque la lumière arrive sur la surface selon la direction normale, c'est-à-dire  <math>\vec L = \vec N</math>

La composante spéculaire, <math>I_s</math> est donnée par : 

:<math>I_s = i_s k_s (\vec R \cdot \vec V)^\alpha = i_s k_s \cos^\alpha \Omega</math>
:Si <math>\vec V = \vec R</math>, <math>I_s</math> est maximale. Plus <math>\alpha</math> est grand plus les taches lumineuses seront petites.

L'intensité réfléchie totale est ainsi :

:<math>I = I_a + I_d + I_s</math>

Pour plusieurs sources lumineuses, on obtient comme formule complète :

:<math>I = i_a k_a + \sum_{n \in \{sources\} } \Big( i_{d,n} k_d (\vec L_n \cdot \vec N) + i_{s,n} k_s (\vec R_n \cdot \vec V)^\alpha \Big) </math>

=== Application à l'infographie ===
Lorsqu'on souhaite appliquer cet algorithme en infographie, on sépare les composantes rouge/vert/bleu de la couleur de la texture pour le point considéré. Puis on applique la formule à <math>i_a = i_d = i_{texture}</math> pour les composantes ambiante et diffuse.

== Les failles du modèle de Phong ==
La composante spéculaire se base sur deux vecteurs directionnels, l'observateur et la lumière, et interdit toute [[radiosité]] du modèle. Ce modèle est empirique et n'est fondé sur aucune théorie physique, mais uniquement sur les observations de Phong. Ce modèle ne prévoit pas la diffusion de la lumière avec la distance.

== L'interpolation de Phong ==
[[Image:Phong shading.png|right|frame|Une sphère interpolée et rendue selon Phong.]]
Cette méthode, à mettre en parallèle avec l'[[ombrage de Gouraud|interpolation de Gouraud]], produit de bons résultats de rendus, souvent plus réalistes que son prédécesseur.

Le principal problème de l'ombrage de Gouraud est qu'il ne calcule que les sommets des polygones : une source lumineuse spéculaire placée au centre d'un triangle n'apparaîtra pas. Ce problème est réglé avec l'interpolation de Phong.

Il est à noter que, dans un triangle, les trois sommets ont la même normale, ainsi d'ailleurs qu'en tout point du triangle. Pour un rendu plus lisse, la normale de chaque sommet d'un triangle peut être calculée en faisant la moyenne de toutes les surfaces utilisant ce sommet. On se retrouve alors avec trois normales qui peuvent être différentes.

Soient trois sommets distincts : ''v''<sub>1</sub>, ''v''<sub>2</sub> et ''v''<sub>3</sub>, ayant pour vecteurs unitaires normaux ''n''<sub>1</sub>, ''n''<sub>2</sub> et ''n''<sub>3</sub>. Comme pour l'interpolation de Gouraud, celle-ci se fait de façon linéaire sur toute la surface du triangle <math>v_1v_2v_3</math>, seulement elle se fait depuis les trois vecteurs normaux des sommets, c'est-à-dire que nous interpolons en fait les vecteurs normaux au lieu des couleurs.

À la différence de l'interpolation de Gouraud cependant, le calcul ne se fait pas sur trois points par surface, mais pour tous les points d'une surface {{Incise|ou plus raisonnablement, sur plusieurs subdivisions de points. Cette méthode, bien plus lente, est parfois traitée directement par le matériel, via les ''[[shader]]s''|stop}}.

== Voir aussi ==
===Articles connexes===
* [[Ombrage de Gouraud]]
* [[Bùi Tường Phong]]

=== Bibliographie ===
* {{en}} [[Bùi Tường Phong|Bui Tuong Phong]], « {{Langue|en|Illumination for Computer Generated Pictures}} », dans ''[[Communications of the ACM]]'', {{vol.}}18, {{numéro}}6, {{date-|juin 1975}} {{doi|10.1145/360825.360839}}.

== Liens externes ==
{{Liens}}

{{Palette|Types d'ombrage}}
{{Portail|imagerie numérique}}

[[Catégorie:Imagerie numérique]]
[[Catégorie:3D]]