View source for Physically-based BRDF

== Introduction ==

Physically based rendering have been known for many years, but the "ad-hoc" rendering models (such as Phong) are still widely used  in game. These "ad-hoc" models require laborious tweaking to produce high-quality images. However, physically based, energy-conserving rendering models easily create materials that hold up under a variety of lighting environments.

'''Surprisingly''', physically based models are not more difficult to implement or evaluate than the traditional "ad-hoc" ones.

Since [[KlayGE]] 3.11.0, there is a calc_brdf() function in Lighting.fxml. It can be used for calculating the [[physically-based BRDF]] described by this article.

== Reflectance equation ==

The most common used rendering model in game describes only reflectance, not including terms such as SSS. The reflectance equation is:

<center><math>L_0(\mathbf{v})=\int_{\Omega} \rho(\mathbf{l},\mathbf{v}) \otimes L_i(\mathbf{l}) (\mathbf{n} \cdot \mathbf{l}) d\omega_i</math></center>

Here ρ('''l''', '''v''') is BRDF, L<sub>i</sub>('''l''') is the contribution from light source, ('''n'''·'''l''') is the angle between light and surface normal. This integration results the sum of all light sources contribute to a surface point.

== Diffuse term ==

The simplest BRDF is the Lambert. The well-known Lambertian BRDF in game is present as ('''n'''·'''l'''). However, it is part of reflectance equation, and lambertian term is actually a constant value:

<center><math>\rho_{lambert}(\mathbf{l},\mathbf{v})=\frac{\mathbf{c}_{diff}}{\pi}</math></center>

== Punctual Light Sources ==

These are classic point, directional, and spot lights in compute games. These local light source can be abstracted to a concept named "''punctual light sources''". They are infinitely small light with centain direction. Since what we want is the lighting reach a surface point, we don't need to consider the attenuation. A punctual light source are parameterized by the light color '''c<sub>light</sub>''' and the light direction vector '''l<sub>c</sub>'''. The light color '''c<sub>light</sub>''' is specified as the color a white Lambertian surface would have when illuminated by the light from a direction parallel to the surface normal ('''l<sub>c</sub>''' = '''n''').

How to calculate the contribution from this light source to a point? Here we will start by defining a tiny area light source centered on '''l<sub>c</sub>''', with a small angular extent ε. This tiny area light illuminates a shaded surface point with the incoming radiance function L<sub>tiny</sub>('''l'''). The incoming radiance function has the following two properties:

<center><math>\forall\mathbf{l}|\angle(\mathbf{l}, \mathbf{l_c})>\varepsilon, L_{tiny}(\mathbf{l})=0</math></center>

<center><math>if \mathbf{l_c}=\mathbf{n}, \mathbf{c}_{light}=\frac{1}{\pi}\int_{\Omega} L_{tiny}(\mathbf{l}) (\mathbf{n} \cdot \mathbf{l}) d\omega_i</math></center>

The first property says that no light is incoming for any light directions which form an angle greater than ε with '''l<sub>c</sub>'''. The second property follows from the definition of '''c<sub>light</sub>'''. Since the surface is white, '''c<sub>diff</sub>''' = 1. Applying reflectance equation and Lambert from last article, there's the second property. Because '''c<sub>light</sub>''' requires '''l<sub>c</sub>''' = '''n''', it is the limit as ε goes to 0:

<center><math>if \mathbf{l_c}=\mathbf{n}, \mathbf{c}_{light}=\lim_{\varepsilon \to 0}(\frac{1}{\pi}\int_{\Omega} L_{tiny}(\mathbf{l}) (\mathbf{n} \cdot \mathbf{l}) d\omega_i)</math></center>

Since '''l<sub>c</sub>''' = '''n''' and ε→0, we can assume ('''n'''·'''l''') = 1, which gives us:

<center><math>\mathbf{c}_{light}=\lim_{\varepsilon \to 0}(\frac{1}{\pi}\int_{\Omega} L_{tiny}(\mathbf{l}) d\omega_i)</math></center>

That is:

<center><math>\lim_{\varepsilon \to 0}(\int_{\Omega} L_{tiny}(\mathbf{l}) d\omega_i) = \pi \mathbf{c}_{light}</math></center>

Now we shall apply our tiny area light to a general BRDF, and look at its behavior in the limit as goes to 0:

<center><math>L_0(\mathbf{v})=\lim_{\varepsilon \to 0}(\int_{\Omega} \rho(\mathbf{l},\mathbf{v}) \otimes L_{tiny}(\mathbf{l}) (\mathbf{n} \cdot \mathbf{l}) d\omega_i)=\rho(\mathbf{l_c}, \mathbf{v}) \otimes \lim_{\varepsilon \to 0}(\int_{\Omega} L_{tiny}(\mathbf{l}) d\omega_i)(\mathbf{n} \cdot \mathbf{l_c})</math></center>

So

<center><math>L_0(\mathbf{v})=\pi \rho(\mathbf{l_c}, \mathbf{v}) \otimes \mathbf{c}_{light} (\mathbf{n} \cdot \mathbf{l_c})</math></center>

As you see, the tiny area light source term disappears. The remaining of our equation is what we familar with.

== Microfacet BRDF ==

The basis for most physically-based specular BRDF terms is ''microfacet theory''. This theory was developed to describe surface reflection from general (non-optically flat) surfaces. The basic assumption underlying microfacet theory is that the surface is composed of many microfacets, too small to be seen individually. Each microfacet is assumed to be optically flat.

Each of these microfacets reflects light from a given incoming direction into a single outgoing direction which depends on the orientation of the microfacet normal '''m'''. When evaluating a BRDF term, both the light direction '''l''' and the view direction '''v''' are specified. This means that of all the millions of microfacets on the surface, only those that happen to be angled just right to reflect '''l''' into '''v''' have any contribution to the BRDF value. In the figure below, we can see that these ''active microfacets'' have their surface normal '''m''' oriented exactly halfway between '''l''' and '''v'''. We will denote it as '''h'''.

[[File:microfacet.jpg|400px|thumb|center|Microfacet]]

Not all microfacets for which '''m''' = '''h''' will contribute to the reflection; some are blocked by other microfacets from the direction of '''l''' (shadowing), from the direction of '''v''' (masking), or from both. Microfacet theory assumes that all shadowed light is lost from the specular term. The microfacet specular BRDF term has the following form:

<center><math>\rho(\mathbf{l}, \mathbf{v})=\frac{F(\mathbf{l},\mathbf{h})G(\mathbf{l},\mathbf{v},\mathbf{h})D(\mathbf{h})}{4(\mathbf{n} \cdot \mathbf{l})(\mathbf{n} \cdot \mathbf{v})}</math></center>

F('''l''', '''h''') is the Fresnel reflectance of the ''active microfacets'' ('''m''' = '''h'''). G('''l''', '''v''', '''h''') is the proportion of active microfacets which are not shadowed or masked. D('''h''') is the microfacet normal distribution function evaluated at the active microfacet normal '''m''' = '''h'''; in other words, the concentration of microfacets with normals equal to '''h'''. Finally, the denominator 4('''n'''·'''l''')('''n'''·'''v''') is a correction factor which accounts for quantities being transformed between the local space of the microfacets and that of the overall macrosurface.

== Fresnel Term ==

The Fresnel reflectance term computes the fraction of light reflected from an optically flat surface. Its value depends on two things: the incoming angle (angle between light vector and surface normal) and the refractive index of the material. The Schlick approximation of Fresnel term is:

<center><math>F_{Schlick}(\mathbf{c}_{spec}, \mathbf{l}, \mathbf{h})=\mathbf{c}_{spec}+(1-\mathbf{c}_{spec})(1-(\mathbf{l} \cdot \mathbf{h})^5)</math></center>

which '''c'''<sub>spec</sub> is the specular color.

== Normal Distribution Term ==

In most surfaces, the microfacet's orientations are not uniformly distributed. Microfacet normals closer to the macroscopic surface normal tend to appear with smoother. The exact distribution is defined via the ''microfacet normal distribution function'' D('''m'''). The function D() determines the size, brightness, and shape of the specular highlight. Normal distribution functions have some kind of "roughness" or variance parameter (anisotropic functions typically have two variance parameters).

== Shadowing-Masking Term ==

The shadowing and masking term G('''l''', '''v''', '''h''') is also often called the ''geometry term''. In many cases, the shadowing-masking term partially cancels out the ('''n'''·'''l''')('''n'''·'''v''') denominator in microfacet BRDF equation, replacing it with some other expression such as max('''n'''·'''l''', '''n'''·'''v'''). The shadowing-masking term is essential for BRDF energy conservation--without such a term the BRDF can reflect arbitrarily more light energy than it receives. A key part of the microfacet BRDF is the ratio between the active area (combined area of the microfacets which reflect light energy from '''l''' to '''v''') and the total area (of the macroscopic surface). If shadowing and masking are not accounted for, then the active area may exceed the total area, an obvious impossibility which can lead to the BRDF not conserving energy, in some cases by a huge amount, as you see below:

[[File:shadowing_masking_term.jpg|400px|thumb|center|Shadowing and masking]]

On the top the flat macroscopic surface is shown in green, and the rugged microscopic surface is shown in blue. The facets for which '''m''' = '''h''' are marked in red. The projection of the macroscopic surface area onto the view direction is shown as a green line on the upper left. The projected areas of the individual red microfacets are shown as separate red lines. On the bottom left the areas of the red microfacets are added up without accounting for masking, resulting in an active area greater than the total area. So the BRDF illogically reflects more energy than it receives. On the right we see that the red areas are combined in a way that accounts for masking. The overlapping areas are no longer counted multiple times, and we see that the correct active area is smaller than the total area.

== Comparison ==

The reflectance equation common used in game look like this:

<center><math>L_{o}(\mathbf{v})=(\mathbf{c}_{diff} (\mathbf{n} \cdot \mathbf{l_c}) + \mathbf{c}_{spec}(\mathbf{r_v} \cdot \mathbf{l_c})^{\alpha}) \otimes \mathbf{c}_{light}</math></center>

According to Lambert equation, ('''n'''·'''l<sub>c</sub>''') is part of reflectance equation, not diffuse term. So it should be:

<center><math>L_{o}(\mathbf{v})=(\mathbf{c}_{diff} + \mathbf{c}_{spec}(\mathbf{r_v} \cdot \mathbf{l_c})^{\alpha}) \otimes \mathbf{c}_{light} (\mathbf{n} \cdot \mathbf{l_c})</math></center>

where '''r<sub>v</sub>''' is the view vector reflected about the normal. It's from Phong model. It doesn't seem to correspond to anything from microfacet theory before. The Blinn-Phong model is very similar to this equation, but it uses the more physically meaningful half-vector '''h''', which is microfacet's normal '''m'''. Changing from Phong to Blinn-Phong gives us the following model:

<center><math>L_{o}(\mathbf{v})=(\mathbf{c}_{diff} + \mathbf{c}_{spec}(\mathbf{n} \cdot \mathbf{h})^{\alpha}) \otimes \mathbf{c}_{light} (\mathbf{n} \cdot \mathbf{l_c})</math></center>

Recall the Lambert term, the diffuse term now is the same with physically-based one:

<center><math>\mathbf{c}_{diff} = \pi \frac {\mathbf{c}_{diff}} {\pi}</math></center>

So, now we focus on specular term. If we look at the specular term only in the equation before:

<center><math>L_{o}(\mathbf{v})=(\mathbf{n} \cdot \mathbf{h})^{\alpha} \mathbf{c}_{spec} \otimes \mathbf{c}_{light} (\mathbf{n} \cdot \mathbf{l_c})</math></center>

Compare to microfacet specular:

<center><math>L_{o}(\mathbf{v})=\pi \frac{D(\mathbf{h})G(\mathbf{l_c},\mathbf{v},\mathbf{h})}{4(\mathbf{n} \cdot \mathbf{l_c})(\mathbf{n} \cdot \mathbf{v})}F(\mathbf{c}_{spec}, \mathbf{l_c},\mathbf{h}) \otimes \mathbf{c}_{light} (\mathbf{n} \cdot \mathbf{l_c})</math></center>

Firstly, we can see from these two equations, the '''c'''<sub>spec</sub> in game can be replaced by Fresnel's F('''c'''<sub>spec</sub>, '''l<sub>c</sub>''', '''h''').

Secondly, I mentioned before that D('''h''') is related to roughness. So D('''h''') and Blinn-Phong's ('''n'''·'''h''')<sup>α</sup> have similar functions. However, to convert the ('''n'''·'''h''')<sup>α</sup> term into a microfacet normal distribution function it must be correctly normalized. It means the sum of the signed projected areas of the microfacets needs to equal the signed projected area of the macroscopic surface; this must hold true for any viewing direction. Mathematically, this means that the function must fulfil this equation for any '''v''':

<center><math>(\mathbf{v} \cdot \mathbf{n}) = \int_{\Theta} D(\mathbf{m})(\mathbf{v} \cdot \mathbf{m}) d\omega_m</math></center>

This equation holds for any viewing direction. In the special case, '''v''' = '''n''':

<center><math>1 = \int_{\Theta} D(\mathbf{m})(\mathbf{n} \cdot \mathbf{m}) d\omega_m</math></center>

The Blinn-Phong cosine power term can be made to obey this equation by multiplying it with a simple normalization factor:

<center><math>D_{BP}=\frac {\alpha + 2} {2\pi}(\mathbf{n} \cdot \mathbf{m})^{\alpha}</math></center>

The only term left is:

<center><math>\frac {G(\mathbf{l_c},\mathbf{v},\mathbf{h})}{(\mathbf{n} \cdot \mathbf{l_c})(\mathbf{n} \cdot \mathbf{v})}</math></center>

It's shadowing/masking, or geometry term divided by the foreshortening factors. The physical meaning of this term is visibility. Since our modified specular model has no visibility term, we will simply set it to 1.  That is:

<center><math>G_{implicit}(\mathbf{l_c},\mathbf{v},\mathbf{h}) = (\mathbf{n} \cdot \mathbf{l_c})(\mathbf{n} \cdot \mathbf{v})</math></center>

This is actually a plausible geometry term for a heightfield microsurface (which is what the Blinn-Phong normal distribution function corresponds to, since it is zero for all backfacing microfacets). G<sub>implicit</sub>() is equal to 1 when '''l''' = '''n''' and '''v''' = '''n''', which is correct for a heightfield (no microfacets are occluded from the direction of the macrosurface normal). It goes to 0 for either glancing view angles or glancing light angles, which again is correct (the probability of a microfacet being occluded by other microfacets increases with viewing angle).

Until now, several terms from microfacet specular BRDF are all explained, leaving us this rendering model:

<center><math>L_{o}(\mathbf{v})=\frac {\alpha + 2} {8}(\mathbf{n} \cdot \mathbf{h})^{\alpha} F(\mathbf{c}_{spec}, \mathbf{l_c},\mathbf{h}) \otimes \mathbf{c}_{light} (\mathbf{n} \cdot \mathbf{l_c})</math></center>

The factor is because of ((α + 2) / 2π) * (π / 4) = (α + 2) / 8.

Add the diffuse term back, that is:

<center><math>L_{o}(\mathbf{v})=(\mathbf{c}_{diff} + \frac {\alpha + 2} {8}(\mathbf{n} \cdot \mathbf{h})^{\alpha} F(\mathbf{c}_{spec}, \mathbf{l_c},\mathbf{h})) \otimes \mathbf{c}_{light} (\mathbf{n} \cdot \mathbf{l_c})</math></center>

This is the final equation, which is '''physically-based rendering model''', derivated by this serial. Compare to the model common used in game, the calculated amount increases little. However, this equation not only ensure the energy-preserving, but also easier for artists to adjust parameters.

== Future works ==

A unsolved problem occurs in the context of very smooth surfaces with high specular powers. Such materials are important to model e.g., wet surfaces. However, the punctual light approximation breaks down in this case, yielding extremely intense highlights of subpixel size that are unrealistic and alias badly. What we would like to see is a sharp reflection of the shape of the light source, which requires some kind of area light approximation which is fast enough to use in games.

Another problem is that there are a variety of geometry terms in the literature. Do any of them provide a visual improvement over the "cheaper-than-free" implicit geometry function G<sub>implicit</sub> that is worth the extra cost? One candidate is the geometry factor proposed by Kelemen et. al. This is an approximation to the Cook-Torrance geometry factor but it is far cheaper to compute:

<center><math>\frac {G_{CT}(\mathbf{l_c}, \mathbf{v}, \mathbf{h})} {(\mathbf{n} \cdot \mathbf{l_c})(\mathbf{n} \cdot \mathbf{v})} \approx \frac {1} {(\mathbf{l_c} \cdot \mathbf{h})^2}</math></center>

note that ('''l<sub>c</sub>'''·'''h''') also exists in any case for Fresnel. It's a close approximation to the full Cook-Torrance geometry factor  divided by the foreshortening terms.

The third problem is the Fresnel of diffuse term. It should be noted that there is no such thing as "a surface without specular". In reality such materials have '''c'''<sub>spec</sub> values around 0.03-0.06, and very low values of α (around 0.1-2.0). At glancing angles, even the most "matte" surfaces have noticeable specular appearance; the lack of this effect is another reason why so many game environments appear unrealistic.

[[Category:Rendering]]
[[Category:Knowledge base]]

[[zh:基于物理的BRDF]]