Rendering equation

In

The physical basis for the rendering equation is the law of conservation of energy. Assuming that L denotes radiance, we have that at each particular position and direction, the outgoing light (L_o) is the sum of the emitted light (L_e) and the reflected light (L_r). The reflected light itself is the sum from all directions of the incoming light (L_i) multiplied by the surface reflection and cosine of the incident angle.

Equation form

The rendering equation may be written in the form

${\displaystyle L_{\text{o}}(\mathbf {x} ,\omega _{\text{o}},\lambda ,t)=L_{\text{e}}(\mathbf {x} ,\omega _{\text{o}},\lambda ,t)+L_{\text{r}}(\mathbf {x} ,\omega _{\text{o}},\lambda ,t)}$

${\displaystyle L_{\text{r}}(\mathbf {x} ,\omega _{\text{o}},\lambda ,t)=\int _{\Omega }f_{\text{r}}(\mathbf {x} ,\omega _{\text{i}},\omega _{\text{o}},\lambda ,t)L_{\text{i}}(\mathbf {x} ,\omega _{\text{i}},\lambda ,t)(\omega _{\text{i}}\cdot \mathbf {n} )\operatorname {d} \omega _{\text{i}}}$

where

• ${\displaystyle L_{\text{o}}(\mathbf {x} ,\omega _{\text{o}},\lambda ,t)}$ is the total spectral radiance of wavelength ${\displaystyle \lambda }$ directed outward along direction ${\displaystyle \omega _{\text{o}}}$ at time ${\displaystyle t}$, from a particular position ${\displaystyle \mathbf {x} }$
• ${\displaystyle \mathbf {x} }$ is the location in space
• ${\displaystyle \omega _{\text{o}}}$ is the direction of the outgoing light
• ${\displaystyle \lambda }$ is a particular wavelength of light
• ${\displaystyle t}$ is time
• ${\displaystyle L_{\text{e}}(\mathbf {x} ,\omega _{\text{o}},\lambda ,t)}$ is emitted spectral radiance
• ${\displaystyle L_{\text{r}}(\mathbf {x} ,\omega _{\text{o}},\lambda ,t)}$ is reflected spectral radiance
• ${\displaystyle \int _{\Omega }\dots \operatorname {d} \omega _{\text{i}}}$ is an integral over ${\displaystyle \Omega }$
• ${\displaystyle \Omega }$ is the unit hemisphere centered around ${\displaystyle \mathbf {n} }$ containing all possible values for ${\displaystyle \omega _{\text{i}}}$ where ${\displaystyle \omega _{\text{i}}\cdot \mathbf {n} >0}$
• ${\displaystyle f_{\text{r}}(\mathbf {x} ,\omega _{\text{i}},\omega _{\text{o}},\lambda ,t)}$ is the bidirectional reflectance distribution function, the proportion of light reflected from ${\displaystyle \omega _{\text{i}}}$ to ${\displaystyle \omega _{\text{o}}}$ at position ${\displaystyle \mathbf {x} }$, time ${\displaystyle t}$, and at wavelength ${\displaystyle \lambda }$
• ${\displaystyle \omega _{\text{i}}}$ is the negative direction of the incoming light
• ${\displaystyle L_{\text{i}}(\mathbf {x} ,\omega _{\text{i}},\lambda ,t)}$ is spectral radiance of wavelength ${\displaystyle \lambda }$ coming inward toward ${\displaystyle \mathbf {x} }$ from direction ${\displaystyle \omega _{\text{i}}}$ at time ${\displaystyle t}$
• ${\displaystyle \mathbf {n} }$ is the surface normal at ${\displaystyle \mathbf {x} }$
• ${\displaystyle \omega _{\text{i}}\cdot \mathbf {n} }$ is the weakening factor of outward irradiance due to incident angle, as the light flux is smeared across a surface whose area is larger than the projected area perpendicular to the ray. This is often written as ${\displaystyle \cos \theta _{i}}$.

Two noteworthy features are: its linearity—it is composed only of multiplications and additions, and its spatial homogeneity—it is the same in all positions and orientations. These mean a wide range of factorings and rearrangements of the equation are possible. It is a Fredholm integral equation of the second kind, similar to those that arise in quantum field theory.^[3]

Note this equation's spectral and time dependence — ${\displaystyle L_{\text{o}}}$ may be sampled at or integrated over sections of the

${\displaystyle t;}$

${\displaystyle L_{\text{o}}}$

Note that a solution to the rendering equation is the function ${\displaystyle L_{\text{o}}}$. The function ${\displaystyle L_{\text{i}}}$ is related to ${\displaystyle L_{\text{o}}}$ via a ray-tracing operation: The incoming radiance from some direction at one point is the outgoing radiance at some other point in the opposite direction.

Applications

Solving the rendering equation for any given scene is the primary challenge in

• Transmission, which occurs when light is transmitted through the surface, such as when it hits a glass object or a water
  surface,
• Subsurface scattering, where the spatial locations for incoming and departing light are different. Surfaces rendered without accounting for subsurface scattering may appear unnaturally opaque — however, it is not necessary to account for this if transmission is included in the equation, since that will effectively include also light scattered under the surface,
• Polarization, where different light polarizations will sometimes have different reflection distributions, for example when light bounces at a water surface,
• Phosphorescence, which occurs when light or other electromagnetic radiation is absorbed at one moment and emitted at a later moment, usually with a longer wavelength (unless the absorbed electromagnetic radiation is very intense),
• Interference
  , where the wave properties of light are exhibited,
• Fluorescence, where the absorbed and emitted light have different wavelengths,
• emission
  of light with higher frequency than the frequency of the light that hit the surface suddenly becomes possible, and
• blueshifted and the photons will be packed more closely so the photon flux will be increased; if it bounces off an object moving away from it, it will be redshifted and the photon flux will be decreased. This effect becomes apparent only at speeds comparable to the speed of light
  , which is not the case for most rendering applications.

For scenes that are either not composed of simple surfaces in a vacuum or for which the travel time for light is an important factor, researchers have generalized the rendering equation to produce a volume rendering equation[5] suitable for volume rendering and a transient rendering equation^[6] for use with data from a time-of-flight camera.

References