Editing Phong reflection model

{{Short description|Shading algorithm in computer graphics}}
{{distinguish|Phong shading}}
[[Image:Phong components version 4.png|thumb|700px|Visual illustration of the Phong equation: here the light is white, the ambient and diffuse colors are both blue, and the specular color is white, reflecting a small part of the light hitting the surface, but only in very narrow highlights. The intensity of the diffuse component varies with the direction of the surface, and the ambient component is uniform (independent of direction).]]

The '''Phong reflection model''' (also called '''Phong illumination''' or '''Phong lighting''') is an [[Empirical modelling|empirical model]] of the [[illumination model|local illumination]] of points on a [[Surface (mathematics)|surface]] designed by the computer graphics researcher [[Bui Tuong Phong]]. In [[3D computer graphics]], it is sometimes referred to as <!--DO NOT LINK-->"Phong shading"<!--DO NOT LINK-->, particularly if the model is used with the [[Phong shading|interpolation method of the same name]] and in the context of [[pixel shader]]s or other places where a lighting calculation can be referred to as “[[shading]]”.

== History ==

The Phong reflection model was developed by [[Bui Tuong Phong]] at the [[University of Utah]], who published it in his 1975 Ph.D. dissertation.<ref>Bui Tuong Phong, [http://www.cs.northwestern.edu/~ago820/cs395/Papers/Phong_1975.pdf Illumination for computer generated pictures], Communications of ACM 18 (1975), no. 6, 311–317.</ref><ref>University of Utah School of Computing, http://www.cs.utah.edu/school/history/#phong-ref</ref> It was published in conjunction with a method for interpolating the calculation for each individual [[pixel]] that is rasterized from a polygonal surface model; the interpolation technique is known as [[Phong shading]], even when it is used with a reflection model other than Phong's. Phong's methods were considered radical at the time of their introduction, but have since become the de facto baseline shading method for many rendering applications. Phong's methods have proven popular due to their generally efficient use of computation time per rendered pixel.

== Concepts ==

Phong reflection is an empirical model of local illumination. It describes the way a surface reflects light as a combination of the [[diffuse reflection]] of rough surfaces with the [[specular reflection]] of shiny surfaces. It is based on Phong's informal observation that shiny surfaces have small intense [[specular highlight]]s, while dull surfaces have large highlights that fall off more gradually. The model also includes an ''ambient'' term to account for the small amount of light that is scattered about the entire scene.

For each light source in the scene, components <math>i_\text{s}</math> and <math>i_\text{d}</math> are defined as the intensities (often as [[RGB]] values) of the specular and diffuse components of the light sources, respectively. A single term <math>i_\text{a}</math> controls the ambient lighting; it is sometimes computed as a sum of contributions from all light sources.

For each ''material'' in the scene, the following parameters are defined:

:<math>k_\text{s}</math>, which is a specular reflection constant, the ratio of reflection of the specular term of incoming light,
:<math>k_\text{d}</math>, which is a diffuse reflection constant, the ratio of reflection of the diffuse term of incoming light ([[Lambertian reflectance]]),
:<math>k_\text{a}</math>, which is an ambient reflection constant, the ratio of reflection of the ambient term present in all points in the scene rendered, and
:<math>\alpha</math>, which is a ''shininess'' constant for this material, which is larger for surfaces that are smoother and more mirror-like.  When this constant is large the specular highlight is small.

[[File:Blinn_Vectors.svg|thumb|Vectors for calculating Phong and Blinn–Phong shading]]

Furthermore, there is

:<math>\text{lights}</math>, which is the [[set (mathematics)|set]] of all light sources,
:<math>\hat{L}_m</math>, which is the direction vector from the point on the surface toward each light source (<math>m</math> specifies the light source),
:<math>\hat{N}</math>, which is the [[Surface normal|normal]] at this point on the surface,
:<math>\hat{R}_m</math>, which is the direction that a perfectly reflected ray of light would take from this point on the surface, and
:<math>\hat{V}</math>, which is the direction pointing towards the viewer (such as a virtual camera).

Then the Phong reflection model provides an equation for computing the illumination of each surface point <math>I_\text{p}</math>:

:<math>I_\text{p} = k_\text{a} i_\text{a} + \sum_{m\;\in\;\text{lights}} (k_\text{d} (\hat{L}_m \cdot \hat{N}) i_{m,\text{d}} + k_\text{s} (\hat{R}_m \cdot \hat{V})^{\alpha}i_{m,\text{s}}).</math>
<!--  WHICH variable is the index of summation in this sum, running through the list of all members of the set called "lights"?? -->
<!-- SEVERAL variables depend on the light:  direction vectors L and R, intensities i_d and i_s.  All of these could be subscripted with a light index if desired, but it's not needed. -->

where the direction vector <math>\hat{R}_m</math> is calculated as the [[Specular reflection#Direction of reflection|reflection]] of <math>\hat{L}_m</math> on the surface characterized by the surface normal <math>\hat{N}</math> using

:<math>\hat{R}_m = 2(\hat{L}_m\cdot \hat{N})\hat{N} - \hat{L}_m</math>

The hats indicate that the vectors are [[Unit vector|normalized]]. The diffuse term is not affected by the viewer direction (<math>\hat{V}</math>).  The specular term is large only when the viewer direction (<math>\hat{V}</math>) is aligned with the reflection direction <math>\hat{R}_m</math>.  Their alignment is measured by the <math>\alpha</math> power of the cosine of the angle between them.  The cosine of the angle between the normalized vectors <math>\hat{R}_m</math> and <math>\hat{V}</math> is equal to their [[dot product]].  When <math>\alpha</math> is large, in the case of a nearly mirror-like reflection, the specular highlight will be small, because any viewpoint not aligned with the reflection will have a cosine less than one which rapidly approaches zero when raised to a high power.

Although the above formulation is the common way of presenting the Phong reflection model, each term should only be included if the term's dot product is positive. (Additionally, the specular term should only be included if the dot product of the diffuse term is positive.)

When the color is represented as [[RGB]] values, as often is the case in [[computer graphics]], this equation is typically modeled separately for R, G and B intensities, allowing different reflection constants <math>k_\text{a},</math> <math>k_\text{d}</math> and <math>k_\text{s}</math> for the different [[Channel (digital image)|color channels]].

When implementing the Phong reflection model, there are a number of methods for approximating the model, rather than implementing the exact formulas, which can speed up the calculation; for example, the [[Blinn–Phong reflection model]] is a modification of the Phong reflection model, which is more efficient if the viewer and the light source are treated to be at infinity.

Another approximation<ref>{{cite web |url=http://dicklyon.com/tech/Graphics/Phong_TR-Lyon.pdf |title=Phong Shading Reformulation for Hardware Renderer Simplification |last1=Lyon |first1=Richard F. |date= August 2, 1993 |access-date=7 March 2011}}</ref> that addresses the calculation of the exponentiation in the specular term is the following: Considering that the specular term should be taken into account only if its dot product is positive, it can be approximated as

:<math>\max(0, \hat{R}_m \cdot \hat{V})^\alpha = \max(0, 1-\lambda)^{\beta \gamma} = \left(\max(0,1-\lambda)^\beta\right)^\gamma \approx \max(0, 1 - \beta \lambda)^\gamma </math>

where <math>\lambda = 1 - \hat{R}_m \cdot \hat{V}</math>, and <math>\beta = \alpha / \gamma\,</math> is a real number which doesn't have to be an integer. If <math>\gamma</math> is chosen to be a power of 2, i.e. <math>\gamma = 2^n</math> where <math>n</math> is an integer, then the expression <math>(1 - \beta\lambda)^\gamma</math> can be more efficiently calculated by squaring <math>(1 - \beta\lambda)\ n</math> times, i.e.

:<math>(1 - \beta\lambda)^\gamma \,=\, (1 - \beta\lambda)^{2^n} \,=\, (1 - \beta\lambda)^{\overbrace{\scriptstyle 2\,\cdot\,2\,\cdot\,\dots\,\cdot\,2}^n} \,=\, (\dots((1 - \beta\lambda)\overbrace{^2)^2\dots)^2}^n.</math>

This approximation of the specular term holds for a sufficiently large integer <math>\gamma</math> (typically, 4 or 8 will be enough).

Furthermore, the value <math>\lambda</math> can be approximated as <math>\lambda = (\hat{R}_m - \hat{V})\cdot(\hat{R}_m - \hat{V}) / 2</math>, or as <math>\lambda = (\hat{R}_m \times \hat{V})\cdot(\hat{R}_m \times \hat{V}) / 2.</math> The latter is much less sensitive to normalization errors in <math>\hat{R}_m</math> and <math>\hat{V}</math> than Phong's dot-product-based <math>\lambda = 1 - \hat{R}_m \cdot \hat{V}</math> is{{Citation needed|reason=Lyon-1993 does not reference this anywhere|date=April 2022}}, and practically doesn't require <math>\hat{R}_m</math> and <math>\hat{V}</math> to be normalized{{Citation needed|reason=Lyon-1993 does not reference this anywhere|date=April 2022}} except for very low-resolved triangle meshes.

This method substitutes a few multiplications for a variable exponentiation, and removes the need for an accurate reciprocal-square-root-based vector normalization.

== Inverse model ==

The Phong reflection model in combination with [[Phong shading]] is an approximation of shading of objects in real life. This means that the Phong equation can relate the shading seen in a [[photograph]] with the surface normals of the visible object. Inverse refers to the wish to estimate the surface normals given a rendered image, natural or computer-made.

The Phong reflection model contains many parameters, such as the surface diffuse reflection parameter ([[albedo]]) which may vary within the object. Thus the normals of an object in a photograph can only be determined, by introducing additional information such as the number of lights, light directions and reflection parameters.

For example, we have a cylindrical object, for instance a finger, and wish to compute the normal <math>N=[N_x, N_z]</math> on a line on the object. We assume only one light, no specular reflection, and uniform known (approximated) reflection parameters. We can then simplify the Phong equation to:

:<math>I_p(x) = C_a +  C_d (L(x) \cdot N(x)) </math>

With <math>C_a</math> a constant equal to the ambient light and <math>C_d</math> a constant equal to the diffusion reflection. We can re-write the equation to:

:<math>( I_p(x)- C_a ) / C_d  = L(x) \cdot N(x) </math>

Which can be rewritten for a line through the cylindrical object as:

:<math>( I_p- C_a ) / C_d  = L_x N_x + L_z N_z </math>

For instance if the light direction is 45 degrees above the object <math>L=[0.71, 0.71]</math> we get two equations with two unknowns.

:<math>( I_p- C_a ) / C_d  = 0.71 N_x + 0.71 N_z </math>
:<math>1 = \sqrt{(N_x^2 + N_z^2)}</math>

Because of the powers of two in the equation there are two possible solutions for the normal direction. Thus some prior information of the geometry is needed to define the correct normal direction. The normals are directly related to angles of inclination of the line on the object surface. Thus the normals allow the calculation of the relative surface heights of the line on the object using a line integral, if we assume a continuous surface.

If the object is not cylindrical, we have three unknown normal values <math>N=[N_x, N_y, N_z]</math>. Then the two equations still allow the normal to rotate around the view vector, thus additional constraints are needed from prior geometric information. For instance in [[facial recognition system|face recognition]] those geometric constraints can be obtained using [[principal component analysis]] (PCA) on a database of depth-maps of faces, allowing only surface normals solutions which are found in a normal population.<ref>{{cite book| title=Model-Based Illumination Correction for Face Images in Uncontrolled Scenarios|date=September 2009|last1=Boom |first1= B.J. |last2=Spreeuwers|first2=  L.J. |last3=Veldhuis|first3=  R.N.J.| volume=5702| pages=33–40| doi=10.1007/978-3-642-03767-2 | issue=2009| series=Lecture Notes in Computer Science| editor1-last=Jiang| editor1-first=Xiaoyi| editor2-last=Petkov| editor-link2= Nicolai Petkov| editor2-first=Nicolai| isbn=978-3-642-03766-5| url=https://halshs.archives-ouvertes.fr/halshs-00420059/document| bibcode=2009LNCS.5702.....J| hdl=11693/26732| hdl-access=free}}</ref>

== Applications ==

The Phong reflection model is often used together with [[Phong shading]] to shade surfaces in [[3D computer graphics]] software. Apart from this, it may also be used for other purposes. For example, it has been used to model the reflection of [[thermal radiation]] from the [[Pioneer program|Pioneer probes]] in an attempt to explain the [[Pioneer anomaly]].<ref>{{cite journal|title=Modelling the reflective thermal contribution to the acceleration of the Pioneer spacecraft|author1=F. Francisco |author2=O. Bertolami |author3=P. J. S. Gil |author4=J. Páramos |doi=10.1016/j.asr.2011.10.016|year=2012|journal=Advances in Space Research|volume=49|issue=3|pages=337–346|arxiv=1103.5222|bibcode=2012AdSpR..49..579S}}</ref>

== See also ==

* {{anl|Bidirectional reflectance distribution function}}
* {{anl|Blinn–Phong shading model}}
* {{anl|List of common shading algorithms}}
* {{anl|Gamma correction}}
* {{anl|Phong shading}}
* {{anl|Specular highlight}}

== References ==

{{refs}}

==External links==

* [https://web.archive.org/web/20160525135133/http://michal.is/projects/phong-reflection-model-matlab/ Phong reflection model in Matlab]
* {{usurped|1=[https://web.archive.org/web/20180816064924/http://www.sunandblackcat.com/tipFullView.php?l=eng&topicid=30 Phong reflection model in GLSL]}}

[[Category:Computer graphics algorithms]]
[[Category:Shading]]
[[Category:Vietnamese inventions]]