Rendering equation

Last updated January 10, 2024

In computer graphics, the rendering equation is an integral equation in which the equilibrium radiance leaving a point is given as the sum of emitted plus reflected radiance under a geometric optics approximation. It was simultaneously introduced into computer graphics by David Immel et al.^[1] and James Kajiya ^[2] in 1986. The various realistic rendering techniques in computer graphics attempt to solve this equation.

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

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)

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

$L_{\text{o}}(\mathbf {x} ,\omega _{\text{o}},\lambda ,t)$ is the total spectral radiance of wavelength $\lambda$ directed outward along direction $\omega _{\text{o}}$ at time $t$ , from a particular position $\mathbf {x}$
$\mathbf {x}$ is the location in space
$\omega _{\text{o}}$ is the direction of the outgoing light
$\lambda$ is a particular wavelength of light
$t$ is time
$L_{\text{e}}(\mathbf {x} ,\omega _{\text{o}},\lambda ,t)$ is emitted spectral radiance
$L_{\text{r}}(\mathbf {x} ,\omega _{\text{o}},\lambda ,t)$ is reflected spectral radiance
$\int _{\Omega }\dots \operatorname {d} \omega _{\text{i}}$ is an integral over $\Omega$
$\Omega$ is the unit hemisphere centered around $\mathbf {n}$ containing all possible values for $\omega _{\text{i}}$ where $\omega _{\text{i}}\cdot \mathbf {n} >0$
$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 $\omega _{\text{i}}$ to $\omega _{\text{o}}$ at position $\mathbf {x}$ , time $t$ , and at wavelength $\lambda$
$\omega _{\text{i}}$ is the negative direction of the incoming light
$L_{\text{i}}(\mathbf {x} ,\omega _{\text{i}},\lambda ,t)$ is spectral radiance of wavelength $\lambda$ coming inward toward $\mathbf {x}$ from direction $\omega _{\text{i}}$ at time $t$
$\mathbf {n}$ is the surface normal at $\mathbf {x}$
$\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 $\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 — $L_{\text{o}}$ may be sampled at or integrated over sections of the visible spectrum to obtain, for example, a trichromatic color sample. A pixel value for a single frame in an animation may be obtained by fixing $t;$ motion blur can be produced by averaging $L_{\text{o}}$ over some given time interval (by integrating over the time interval and dividing by the length of the interval).^[4]

Note that a solution to the rendering equation is the function $L_{\text{o}}$ . The function $L_{\text{i}}$ is related to $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 realistic rendering. One approach to solving the equation is based on finite element methods, leading to the radiosity algorithm. Another approach using Monte Carlo methods has led to many different algorithms including path tracing, photon mapping, and Metropolis light transport, among others.

Limitations

Although the equation is very general, it does not capture every aspect of light reflection. Some missing aspects include the following:

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,
Non-linear effects, where very intense light can increase the energy level of an electron with more energy than that of a single photon (this can occur if the electron is hit by two photons at the same time), and emission of light with higher frequency than the frequency of the light that hit the surface suddenly becomes possible, and
Doppler effect, where light that bounces off an object moving at a very high speed will get its wavelength changed: if the light bounces off an object that is moving towards it, the light will be 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.

Related Research Articles

$<span class="mw-page-title-main">Diffraction</span> Phenomenon of the motion of waves$

Diffraction is the interference or bending of waves around the corners of an obstacle or through an aperture into the region of geometrical shadow of the obstacle/aperture. The diffracting object or aperture effectively becomes a secondary source of the propagating wave. Italian scientist Francesco Maria Grimaldi coined the word diffraction and was the first to record accurate observations of the phenomenon in 1660.

<span class="mw-page-title-main">Reflectance</span> Capacity of an object to reflect light

The reflectance of the surface of a material is its effectiveness in reflecting radiant energy. It is the fraction of incident electromagnetic power that is reflected at the boundary. Reflectance is a component of the response of the electronic structure of the material to the electromagnetic field of light, and is in general a function of the frequency, or wavelength, of the light, its polarization, and the angle of incidence. The dependence of reflectance on the wavelength is called a reflectance spectrum or spectral reflectance curve.

Specific detectivity, or D*, for a photodetector is a figure of merit used to characterize performance, equal to the reciprocal of noise-equivalent power (NEP), normalized per square root of the sensor's area and frequency bandwidth.

Wien's displacement law states that the black-body radiation curve for different temperatures will peak at different wavelengths that are inversely proportional to the temperature. The shift of that peak is a direct consequence of the Planck radiation law, which describes the spectral brightness or intensity of black-body radiation as a function of wavelength at any given temperature. However, it had been discovered by Wilhelm Wien several years before Max Planck developed that more general equation, and describes the entire shift of the spectrum of black-body radiation toward shorter wavelengths as temperature increases.

<span class="mw-page-title-main">Bremsstrahlung</span> Electromagnetic radiation due to deceleration of charged particles

In particle physics, bremsstrahlung is electromagnetic radiation produced by the deceleration of a charged particle when deflected by another charged particle, typically an electron by an atomic nucleus. The moving particle loses kinetic energy, which is converted into radiation, thus satisfying the law of conservation of energy. The term is also used to refer to the process of producing the radiation. Bremsstrahlung has a continuous spectrum, which becomes more intense and whose peak intensity shifts toward higher frequencies as the change of the energy of the decelerated particles increases.

In physics, Planck's law describes the spectral density of electromagnetic radiation emitted by a black body in thermal equilibrium at a given temperature $T$ , when there is no net flow of matter or energy between the body and its environment.

In optical physics, transmittance of the surface of a material is its effectiveness in transmitting radiant energy. It is the fraction of incident electromagnetic power that is transmitted through a sample, in contrast to the transmission coefficient, which is the ratio of the transmitted to incident electric field.

In radiometry, radiance is the radiant flux emitted, reflected, transmitted or received by a given surface, per unit solid angle per unit projected area. Radiance is used to characterize diffuse emission and reflection of electromagnetic radiation, and to quantify emission of neutrinos and other particles. The SI unit of radiance is the watt per steradian per square metre. It is a directional quantity: the radiance of a surface depends on the direction from which it is being observed.

In physics, a wave vector is a vector used in describing a wave, with a typical unit being cycle per metre. It has a magnitude and direction. Its magnitude is the wavenumber of the wave, and its direction is perpendicular to the wavefront. In isotropic media, this is also the direction of wave propagation.

In radiometry, radiant intensity is the radiant flux emitted, reflected, transmitted or received, per unit solid angle, and spectral intensity is the radiant intensity per unit frequency or wavelength, depending on whether the spectrum is taken as a function of frequency or of wavelength. These are directional quantities. The SI unit of radiant intensity is the watt per steradian, while that of spectral intensity in frequency is the watt per steradian per hertz and that of spectral intensity in wavelength is the watt per steradian per metre —commonly the watt per steradian per nanometre. Radiant intensity is distinct from irradiance and radiant exitance, which are often called intensity in branches of physics other than radiometry. In radio-frequency engineering, radiant intensity is sometimes called radiation intensity.

The Compton wavelength is a quantum mechanical property of a particle, defined as the wavelength of a photon the energy of which is the same as the rest energy of that particle. It was introduced by Arthur Compton in 1923 in his explanation of the scattering of photons by electrons.

The bidirectional reflectance distribution function (BRDF), symbol $, is a function of four real variables that defines how light from a source is reflected off an opaque surface. It is employed in the optics of real-world light, in computer graphics algorithms, and in computer vision algorithms. The function takes an incoming light direction,, and outgoing direction,, and returns the ratio of reflected radiance exiting along to the irradiance incident on the surface from direction . Each direction is itself parameterized by azimuth angle and zenith angle, therefore the BRDF as a whole is a function of 4 variables. The BRDF has units sr -1, with steradians (sr) being a unit of solid angle.$

In radiometry, radiant exitance or radiant emittance is the radiant flux emitted by a surface per unit area, whereas spectral exitance or spectral emittance is the radiant exitance of a surface per unit frequency or wavelength, depending on whether the spectrum is taken as a function of frequency or of wavelength. This is the emitted component of radiosity. The SI unit of radiant exitance is the watt per square metre, while that of spectral exitance in frequency is the watt per square metre per hertz (W·m⁻²·Hz⁻¹) and that of spectral exitance in wavelength is the watt per square metre per metre (W·m⁻³)—commonly the watt per square metre per nanometre. The CGS unit erg per square centimeter per second is often used in astronomy. Radiant exitance is often called "intensity" in branches of physics other than radiometry, but in radiometry this usage leads to confusion with radiant intensity.

The theoretical and experimental justification for the Schrödinger equation motivates the discovery of the Schrödinger equation, the equation that describes the dynamics of nonrelativistic particles. The motivation uses photons, which are relativistic particles with dynamics described by Maxwell's equations, as an analogue for all types of particles.

The linear attenuation coefficient, attenuation coefficient, or narrow-beam attenuation coefficient characterizes how easily a volume of material can be penetrated by a beam of light, sound, particles, or other energy or matter. A coefficient value that is large represents a beam becoming 'attenuated' as it passes through a given medium, while a small value represents that the medium had little effect on loss. The SI unit of attenuation coefficient is the reciprocal metre (m⁻¹). Extinction coefficient is another term for this quantity, often used in meteorology and climatology. Most commonly, the quantity measures the exponential decay of intensity, that is, the value of downward e-folding distance of the original intensity as the energy of the intensity passes through a unit thickness of material, so that an attenuation coefficient of 1 m⁻¹ means that after passing through 1 metre, the radiation will be reduced by a factor of e, and for material with a coefficient of 2 m⁻¹, it will be reduced twice by e, or e². Other measures may use a different factor than e, such as the decadic attenuation coefficient below. The broad-beam attenuation coefficient counts forward-scattered radiation as transmitted rather than attenuated, and is more applicable to radiation shielding. The mass attenuation coefficient is the attenuation coefficient normalized by the density of the material.

Kramers' law is a formula for the spectral distribution of X-rays produced by an electron hitting a solid target. The formula concerns only bremsstrahlung radiation, not the element specific characteristic radiation. It is named after its discoverer, the Dutch physicist Hendrik Anthony Kramers.

In the study of heat transfer, absorptance of the surface of a material is its effectiveness in absorbing radiant energy. It is the ratio of the absorbed to the incident radiant power.

In physics and engineering, the envelope of an oscillating signal is a smooth curve outlining its extremes. The envelope thus generalizes the concept of a constant amplitude into an instantaneous amplitude. The figure illustrates a modulated sine wave varying between an upper envelope and a lower envelope. The envelope function may be a function of time, space, angle, or indeed of any variable.

Heat transfer physics describes the kinetics of energy storage, transport, and energy transformation by principal energy carriers: phonons, electrons, fluid particles, and photons. Heat is energy stored in temperature-dependent motion of particles including electrons, atomic nuclei, individual atoms, and molecules. Heat is transferred to and from matter by the principal energy carriers. The state of energy stored within matter, or transported by the carriers, is described by a combination of classical and quantum statistical mechanics. The energy is different made (converted) among various carriers. The heat transfer processes are governed by the rates at which various related physical phenomena occur, such as the rate of particle collisions in classical mechanics. These various states and kinetics determine the heat transfer, i.e., the net rate of energy storage or transport. Governing these process from the atomic level to macroscale are the laws of thermodynamics, including conservation of energy.

Stimulated Raman spectroscopy, also referred to as stimulated Raman scattering (SRS) is a form of spectroscopy employed in physics, chemistry, biology, and other fields. The basic mechanism resembles that of spontaneous Raman spectroscopy: a pump photon, of the angular frequency $, which is scattered by a molecule has some small probability of inducing some vibrational transition, as opposed to inducing a simple Rayleigh transition. This makes the molecule emit a photon at a shifted frequency. However, SRS, as opposed to spontaneous Raman spectroscopy, is a third-order non-linear phenomenon involving a second photon—the Stokes photon of angular frequency —which stimulates a specific transition. When the difference in frequency between both photons resembles that of a specific vibrational transition the occurrence of this transition is resonantly enhanced. In SRS, the signal is equivalent to changes in the intensity of the pump and Stokes beams. The signals are typically rather low, of the order of a part in 10^5, thus calling for modulation-transfer techniques: one beam is modulated in amplitude and the signal is detected on the other beam via a lock-in amplifier. Employing a pump laser beam of a constant frequency and a Stokes laser beam of a scanned frequency allows for the unraveling of the spectral fingerprint of the molecule. This spectral fingerprint differs from those obtained by other spectroscopy methods such as Rayleigh scattering as the Raman transitions confer to different exclusion rules than those that apply for Rayleigh transitions.$

References

↑ Immel, David S.; Cohen, Michael F.; Greenberg, Donald P. (1986). "A radiosity method for non-diffuse environments" (PDF). In David C. Evans; RussellJ. Athay (eds.). SIGGRAPH '86. Proceedings of the 13th annual conference on Computer graphics and interactive techniques. pp. 133–142. doi:10.1145/15922.15901. ISBN 978-0-89791-196-2. S2CID 7384510.
↑ Kajiya, James T. (1986). "The rendering equation" (PDF). In David C. Evans; RussellJ. Athay (eds.). SIGGRAPH '86. Proceedings of the 13th annual conference on Computer graphics and interactive techniques. pp. 143–150. doi:10.1145/15922.15902. ISBN 978-0-89791-196-2. S2CID 9226468.
↑ Watt, Alan; Watt, Mark (1992). "12.2.1 The path tracing solution to the rendering equation". Advanced Animation and Rendering Techniques: Theory and Practice . Addison-Wesley Professional. p. 293. ISBN 978-0-201-54412-1.
↑ Owen, Scott (September 5, 1999). "Reflection: Theory and Mathematical Formulation" . Retrieved 2008-06-22.
↑ Kajiya, James T.; Von Herzen, Brian P. (1984), "Ray tracing volume densities", ACM SIGGRAPH Computer Graphics, 18 (3): 165–174, CiteSeerX 10.1.1.128.3394 , doi:10.1145/964965.808594
↑ Smith, Adam M.; Skorupski, James; Davis, James (2008). Transient Rendering (PDF) (Technical report). UC Santa Cruz. UCSC-SOE-08-26.

External links

Lecture notes from Stanford University course CS 348B, Computer Graphics: Image Synthesis Techniques

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Immel, David S.; Cohen, Michael F.; Greenberg, Donald P. (1986). "A radiosity method for non-diffuse environments" (PDF). In David C. Evans; RussellJ. Athay (eds.). SIGGRAPH '86. Proceedings of the 13th annual conference on Computer graphics and interactive techniques. pp. 133–142. doi:10.1145/15922.15901. ISBN 978-0-89791-196-2. S2CID 7384510.

[2] Kajiya, James T. (1986). "The rendering equation" (PDF). In David C. Evans; RussellJ. Athay (eds.). SIGGRAPH '86. Proceedings of the 13th annual conference on Computer graphics and interactive techniques. pp. 143–150. doi:10.1145/15922.15902. ISBN 978-0-89791-196-2. S2CID 9226468.

[3] Watt, Alan; Watt, Mark (1992). "12.2.1 The path tracing solution to the rendering equation". Advanced Animation and Rendering Techniques: Theory and Practice . Addison-Wesley Professional. p. 293. ISBN 978-0-201-54412-1.

[4] Owen, Scott (September 5, 1999). "Reflection: Theory and Mathematical Formulation" . Retrieved 2008-06-22.

[5] Kajiya, James T.; Von Herzen, Brian P. (1984), "Ray tracing volume densities", ACM SIGGRAPH Computer Graphics, 18 (3): 165–174, CiteSeerX 10.1.1.128.3394 , doi:10.1145/964965.808594

[6] Smith, Adam M.; Skorupski, James; Davis, James (2008). Transient Rendering (PDF) (Technical report). UC Santa Cruz. UCSC-SOE-08-26.

[1]

[2]

[3]

[4]

[5]

[6]