## Abstract

A popular class of BRDF models is the microfacet models, where geometric optics is assumed. In contrast, more complex physical optics models may more accurately predict the BRDF, but the calculation is more resource intensive. These seemingly disparate approaches are compared in detail for the rough and smooth surface approximations of the modified Beckmann-Kirchhoff BRDF model, assuming Gaussian surface statistics. An approximation relating standard Fresnel reflection with the semi-rough surface polarization term, *Q*, is presented for unpolarized light. For rough surfaces, the angular dependence of direction cosine space is shown to be identical to the angular dependence in the microfacet distribution function. For polished surfaces, the same comparison shows a breakdown in the microfacet models. Similarities and differences between microfacet BRDF models and the modified Beckmann-Kirchhoff model are identified. The rationale for the original Beckmann-Kirchhoff
${F}_{bk}^{2}$ geometric term relative to both microfacet models and generalized Harvey-Shack model is presented. A modification to the geometric
${F}_{bk}^{2}$ term in original Beckmann-Kirchhoff BRDF theory is proposed.

© 2015 Optical Society of America

## 1. Introduction

The bidirectional reflectance distribution function (BRDF) defines the ratio of the reflected radiance to the incident irradiance for materials, and is commonly used in remote sensing [1,2], scene generation [3–5], and computer graphics [6–9]. Nicodemus formally defined the BRDF as [10]

*θ*(representing elevation angle) and

*ϕ*(representing azimuthal angle) as measured from the surface normal; see Fig. 1(a).

*λ*represents the wavelength of light.

We are interested in using the BRDF to render scenes which are both radiometrically and spectrally accurate, yet are computationally inexpensive to implement. Therefore, in this paper, two different approaches to modeling BRDF will be compared: microfacet model BRDFs and scalar wave optics model BRDFs. The former assumes geometric optics only, but results in several easy-to-use closed-form expressions for the BRDF that can be written in a common form [11]. The latter uses linear systems theory to develop a scalar BRDF based on wave optics, but is much more cumbersome to apply to resource-constrained problems. The purpose of this paper is to understand how the physics of a more complete scalar wave optics model may be related to microfacet closed-form BRDF models. All components of microfacet models are directly related to components of a scalar wave optics model through the development of novel relationships between differing elements in these models. Surface reflection is fundamentally a physical optics problem, but understanding how a closed-form geometric optics model relates to a more complete, but more computationally complex, physical optics model may lead to development of models that are more radiometrically and spectrally accurate for applications that require a closed-form approximation, as is desirable for resource-constrained applications of the BRDF such as remote sensing.

To accomplish this task, first these models are introduced. Then, the rough and smooth surface limits of the wave optics model are compared to the common form used by many microfacet models. Relationships are then developed between the different components of the models, suggesting a different approach to modeling the geometric attenuation (shadowing and masking) term in microfacet models, as well as explaining limitations in microfacet models by quantifying differences in the angular terms in the models. Finally, the relationships developed in this paper are used to explain the role of the original Beckmann-Kirchhoff geometric term as compared to both models, and to explain experimental results in [12]. No BRDF model can be reasonably simple and simultaneously perfect, but comparing the microfacet class of BRDF models to a scalar wave optics model and developing relationships between all components of these models may lead to future development of closed-form BRDF models that correct for deficiencies in microfacet models.

## 2. Background

This section discusses two classes of BRDF models: the microfacet class of models and the Modified Beckmann-Kirchhoff (MBK) scalar wave optics model. These models form the basis of comparison in the remainder of the paper. Prior work in [11–15] has set the stage for comparing the simplicity of a closed-form approximation given by a geometric optics BRDF model with the accuracy of a more robust scalar wave optics model such as MBK BRDF theory. A summary of these key results is presented as background.

#### 2.1. Microfacet BRDF models

The microfacet class of BRDF models includes several popular models, such as Torrance-Sparrow [3], Blinn-Phong [6], Cook-Torrance [7], Ward-Duer [16, 17], Hyde [18], Priest-Germer [19], and Beard-Maxwell [4]. These microfacet models for isotropic surfaces can all be written in a general form as [11]

*V*, and the third to Lambertian scattering.

*ρ*,

_{s}*ρ*, and

_{v}*ρ*are fit parameters that are allowed to vary with wavelength and signify the amount of surface, volume, and diffuse (Lambertian) reflectance for a material. For surface-reflecting materials in a single bounce model,

_{d}*ρ*= 0 and

_{v}*ρ*= 0.

_{d}*F*is the Fresnel reflectance as derived in many optics and electromagnetism textbooks such as [20], with the usual parameters of complex refractive index

*ñ*=

*n*+

*iκ*.

*D*is the microfacet surface normal statistical distribution.

_{μ}*G*is a geometric attenuation (shadowing and masking) term.

*σ*is a conversion from scattering cross section to BRDF, defined as

Overviews of popular microfacet models are given in [11, 21]. In this paper, all surfaces are assumed to be isotropic; *ϕ _{i}* = 0 is chosen as the orientation of the

*ϕ*axis. Microfacet coordinates

*θ*(representing the half angle, or angle of the bisector) and

_{h}*θ*(representing the incident angle in the microfacet’s frame of reference) have been used since the Torrance-Sparrow BRDF in 1967 [3], but were only recently formally defined as [22]

_{d}*is a rotation about some axis*

_{a}*a*. This formal definition has the advantage of specifying the

*ϕ*dependence of ${\widehat{\omega}}_{h}$ and ${\widehat{\omega}}_{d}$ in addition to the

*θ*dependence, which is of particular importance for anisotropic surfaces. In the isotropic case, only the

*θ*and

_{h}*θ*angles are computed using equations that appear in many papers in microfacet literature [3, 4, 7, 16, 18, 19, 23], and are given as

_{d}The Gaussian distribution *D _{g}* is used in the Torrance-Sparrow [3], Cook-Torrance [7], Ward-Duer [16, 17] (in the small angle approximation), Priest-Germer [19], and Hyde [18] BRDF models as well as being approximated by the popular cosine lobe distribution used in Blinn-Phong [6], and is given by [11]

*σ*is commonly thought of as statistically describing the microfacet orientation, but is actually a fit parameter which will be shown here to vary with wavelength, indicating a breakdown of the geometric optics model. From geometric optics,

_{g}*σ*represents the probability distribution of microsurface normals. Mathematically, the microfacet surface normal probability distribution can be calculated from the surface height profile by calculating the normal to the derivative of the surface height profile. The angular dependence of the distribution function is primarily driven by tan

_{g}^{2}

*θ*. The remaining terms, and in particular the 1/cos

_{h}^{4}

*θ*term, normalize the Gaussian distribution to ensure it is a proper probability distribution function in projected area (as Trowbridge indicates is the proper normalization in [24]); that is,

_{h}Unlike linear systems models, the closed-form expression for microfacet models presented in Eq. (2) is used for both very rough and smooth surfaces.

#### 2.2. Modified Beckmann-Kirchhoff model

Several different physical optics models have been developed, deriving from the Kirchhoff tangent plane approximation, small perturbation approximation, and other unifying methods; for a 2004 survey of over 30 of these models, see [25]. Elfouhaily noted all of these models had deficiencies. Krywonos then developed Generalized Harvey-Shack (GHS) theory and MBK theory to correct some of these deficiencies.

More recently, a method for estimating surface statistics using the tangent plane approximation and assuming a Gaussian distribution of surface properties was developed in [26]. Although this method did not include out-of-plane scatter, its success in predicting surface statistics lends credence to using a physical optics model based on the tangent plane approximation, which is used by the Beckmann-Kirchhoff (BK) scattering model.

The BK scattering model was formulated in 1963 [27]. Similar to [26], the BK model uses the tangent plane approximation and assumes a Gaussian distribution of surface properties. However, the original BK theory was formulated prior to formal definition of the BRDF in 1977 [10]. This BK model was modified when the linear systems diffraction BRDF theory was generalized to all incident and scattered angles, as well as for a variety of surfaces, by Krywonos [15]. Linear systems diffraction BRDF theory was originally formulated for paraxial reflection in the Harvey-Shack BRDF model, then was extended to arbitrary angles in the GHS BRDF model [14, 15, 28]. The GHS BRDF model developed by Krywonos does not possess a known closed-form solution in general. GHS theory is more physically rigorous, but is more computationally intensive and unfortunately, this method fails to distinguish between scattering due to surface reflection versus scattering due to volumetric effects. In both cases, the surface is treated as adding an overall phase offset, regardless of whether the scatter is due to the surface or due to impurities in the material. Krywonos also proposed MBK BRDF theory in [15,28,29] that is equal to his GHS theory in the limit of either a polished or very rough surface [28, 29]. This formulation assumes Gaussian statistics, but results in an infinite series expression for the BRDF that more readily enables comparison to the microfacet class of models without sacrificing the accuracy of the GHS model, from which the modifications to the original BK theory were derived.

In GHS and MBK theory, the BRDF is calculated in direction cosine space (*α, β*) instead of spherical space (*θ*, *ϕ*). This space is defined as

*π*terms here arise from a difference in defining the

*ϕ*= 0 location in the microfacet model as compared to the linear systems model. For in-plane forward scatter,

*ϕ*=

_{s}*π*in both models as presented here (backscatter is

*ϕ*= 0); however, in-plane forward scatter data is represented by

_{s}*ϕ*= 0 in some linear systems papers [14, 15]. For isotropic samples, the angular dependence simplifies to

_{s}*ν*is the notation used by Beckmann [27].

_{xy}The MBK formulation presented in Eq. (11) is the Angle Spread Function (ASF), *f _{a}*, which is related to the radiance

*L*; the ASF is equal to the BRDF if

*F*= 1 (total Fresnel reflectance). Assuming the surface statistics follow a Gaussian autocovariance, and taking into account the 1

*/λ*

^{2}term present in MBK as shown in [29], the ASF is calculated as

*K*is a renormalization term that ensures conservation of energy,

*l*is the correlation length of the surface, and

_{c}*g*is related to the surface height

*σ*and given as [28]

_{s}When written in terms of scattered radiance *L _{s}*,

*K*is defined as [15, 28, 29]

*L _{s}* ≥ 0 everywhere, so

*K*≥1 for ensures all incident angles. Physically,

*K*is a unitless quantity that energy is not lost to evanescent waves by redistributing energy that may fall outside the unit circle of real space defined by the direction cosines. For some BRDF applications, this term can be problematic to compute because it depends on knowing the surface distribution shape, and involves computation of a double integral expression that does not possess an analytic solution for most surface distributions.

For a smooth surface, *σ _{s} << λ* and as such only the

*m*= 1 term of the sum in Eq. (11) is significant. In [29], Krywonos defines a smooth surface as

*g*(

*θ*) < 0.025, which results in less than 1% error compared to the infinite series summation when

_{i}, θ_{s}*θ*70

_{i}≤*°*. By using Eqs. (10)–(12) and multiplying by the polarization factor

*Q*to convert the ASF to a BRDF, this polished approximation

*f*can be written as

_{p}*g*)

*≈*1 since

*σ*for a polished surface, and

_{s}<< λ*K ≈*1 for a polished surface, as shown in [15].

At the other extreme, if a surface is very rough, the BRDF *f _{vr}* as given in [29] can be rewritten using Eqs. (10)–(12) as

*Q*to convert from ASF to BRDF. In [29], Krywonos defines very rough as

*g*(

*θ*) > 800, to result in less than 1% maximum error as compared to the infinite summation when

_{i}, θ_{s}*θ*70

_{i}≤*°*.

The appropriate Fresnel term to convert ASF to BRDF is the polarization factor *Q* as described by Stover in [30]. The more popular unpolarized Fresnel reflectance term *F* derived in many electromagnetism and optics texts assumes plane wave incidence and an infinitely smooth surface. The derivation of the polarization factor *Q* in [30] relaxes the infinite plane wave and perfectly smooth assumptions in the derivation of *F*. For *s*-polarized light incident on a surface, *Q _{s}* =

*Q*+

_{ss}*Q*. For

_{sp}*p*-polarized light incident on a surface,

*Q*=

_{p}*Q*+

_{ps}*Q*. For unpolarized light incident on a surface,

_{pp}*Q*=

*Q*+

_{s}*Q*. If

_{p}*θ*=

_{i}*θ*and

_{s}*ϕ*= 180

_{s}*°*(specular reflection), the equations for

*Q*and

_{s}*Q*reduce to the standard Fresnel equations

_{p}*F*and

_{s}*F*.

_{p}In [15, 31], careful measurements of material surface statistics were made by Stover. The measured surface statistics were used to compute the Power Spectral Density (PSD) of each surface. Stover then measured the BRDF of these surfaces. The GHS BRDF model was used to compute the PSD from the BRDF data. The computed PSD using GHS theory closely matched the measured surface statistics of the surfaces. Similarly, in [32], Schröder *et al.* used GHS theory to predict the scatter off an arbitrary (but known) surface of arbitrary roughness, both rough and smooth. Since GHS matches MBK in the very rough or very smooth approximation, *f _{p}* and

*f*from MBK theory are considered to be accurate physical models in the limit of a polished or very rough surface, respectively.

_{vr}## 3. Analysis

To compare these BRDF models, the ratio of microfacet BRDFs to the MBK BRDF is taken in both the very rough approximation and the smooth surface approximation, assuming an isotropic sample. We also assume surface statistics (surface height and autocorrelation length) follow a Gaussian distribution.

Since MBK theory can be used to predict accurate surface statistics, the ratios *f _{μ}/f_{p}* and

*f*would equal 1 if the microfacet model also modeled physical truth. Differences from 1 indicate deviations in the microfacet model from the more physical MBK model, and may suggest how to alter the microfacet model to improve its accuracy while still maintaining a closed-form approximation to the BRDF.

_{μ}/f_{vr}#### 3.1. Very rough surface comparison

The ratio of microfacet BRDFs *f _{μ}* to the very rough surface MBK BRDF

*f*is

_{vr}The exponential term will be analyzed first. The angular dependence is contained in two terms: tan ^{2}*θ _{h}* from the microfacet model and
${\eta}_{r}^{2}/{\left(\mathrm{cos}{\theta}_{i}+\mathrm{cos}{\theta}_{s}\right)}^{2}$ from MBK. Although these terms were derived with a completely different physical interpretation (geometric optics versus physical optics), it will be shown here that these terms are equal at all incident and scattered angles, in-plane and out-of-plane. Starting with Eq. (6) and using basic trigonometry identities to rewrite cos

*θ*in terms of tan

_{h}*θ*:

_{h}The denominator already matches the denominator of
${\eta}_{r}^{2}/{\left(\mathrm{cos}{\theta}_{i}+\mathrm{cos}{\theta}_{s}\right)}^{2}$, so we only need to show that the numerator equals
${\eta}_{r}^{2}$ as given by Eq. (10). This is performed by use of the double angle formula for cos(2*θ _{d}*) and adding zero:

The angular term in the exponential from the scalar wave optics model,
${\eta}_{r}^{2}/{\left(\mathrm{cos}{\theta}_{i}+\mathrm{cos}{\theta}_{s}\right)}^{2}$, ranges from 0 to ∞. This is in contrast to the result obtained later in this paper when examining a polished surface. (The angular term in the exponential from the microfacet model, tan^{2}*θ _{h}*, always ranges from 0 to ∞.)

From geometric optics *σ _{g}* is commonly thought of as the Gaussian width of the probability distribution of microsurface normals. Mathematically, the microfacet surface normals can be calculated from the surface height profile and is related to the normal to the derivative of the surface height profile; however, only the probability distribution of the surface profile is specified. Instead, from Eq. (16), we solve for

*σ*for a very rough surface using the remaining terms in the exponential, obtaining

_{g}Recall *σ _{s}* represents the width of the Gaussian in the vertical (height) direction and

*l*represents the width of the Gaussian in the horizontal (correlation length) direction. This result is identical to the result obtained in the Hyde polarized BRDF model’s parameterization of the Gaussian distribution given in [18] that was derived using the Method of Moments, suggesting the methodology employed in this comparison is correct. However, Eq. (20) only holds for the very rough surface approximation. A different result for the polished surface analysis is obtained later in the paper in Eq. (25), showing this interpretation of

_{c}*σ*is only valid for very rough surfaces, and is not true in general. This point will be discussed in more depth after comparing the polished surface approximation to the microfacet model.

_{g}Eq. (16) has now simplified to

Direct comparisons between *Q* and *F* vary significantly in magnitude. Stover explains the *s* polarization, but does not suggest an approximation valid at all angles and indices of refraction for the *p* polarization (and thus the unpolarized term *Q* = *Q _{s}* +

*Q*), although he does present approximations for a few limiting cases [30]. To improve on this understanding, the polarization factor

_{p}*Q*was analyzed from the perspective of comparing it to not only Fresnel reflectance, but also to the remaining angular terms from the microfacet model, since the microfacet model uses

*F*instead of

*Q*.

In [13], we discovered that, for in-plane angles, there is an approximate relationship *S* between the unpolarized Fresnel reflectance and the polarization factor, *Q*, given as

*Q*=

*Q*+

_{s}*Q*=

_{p}*F*+

_{s}*F*= 2

_{p}*F*if

*θ*=

_{i}*θ*and

_{s}*ϕ*= 180

_{s}*°*(standard Fresnel reflection). However, in [13], there was not a solid basis for this approximation. In this paper, the rationale for Eq. (22) is now evident, as this relationship explicitly appears in Eq. (21).

In [13], we also only examined the in-plane relationship. In this paper, a more complete analysis is performed by computing the relative difference between 2*F/Q* and *S*, including at out-of-plane angles, where the relative difference is defined as

Figure 2 was created at *θ _{i}* = (15°, 30°, 45°, 60

*°*) for five different indices of refraction. For each angle, the five different indices of refraction were chosen to represent a wide range of materials: high

*n*and high

*κ*:

*ñ*= 4 + 10

*i*; low

*n*and low (but nonzero)

*κ*:

*ñ*= 1.5 +

*i*, low

*n*and zero

*κ*:

*ñ*= 1.4, moderate

*n*and

*κ*:

*ñ*= 1.7+5

*i*, and

*n*< 1:

*ñ*= 0.25+3

*i*. When

*θ*=

_{i}*θ*and

_{s}*ϕ*= 180

_{s}*°*(specular reflection),

*R*= 0. In general, even out-of-plane, the error is relatively small if

_{d}*θ*and

_{i}*θ*are both small, but increases when either

_{s}*θ*or

_{i}*θ*is large. This result suggests that the cross section conversion term

_{s}*σ*(

*θ*) in Eq. (3) and the microfacet distribution normalization 1/cos

_{i},θ_{s}^{4}

*θ*in Eq. (7) both arise from using

_{h}*F*instead of

*Q*.

This result also suggests a novel interpretation for *G*. In geometric optics, *G* = 1 when the incident and scattered angles are both small–that is, there is no shadowing and masking. This is consistent with the low relative difference observed in Fig. 2 at small incident and scattered angles. However, the relative difference at large angles does not agree with a geometric optics interpretation. Since the BRDF is fundamentally a wave optics problem, better approximations for *G* may be derived by solving 2*FG/Q* = *S* for *G* instead of from the geometric optics approach detailed in [33]; such a solution is beyond the scope of this paper. Alternatively, *Q* does have a closed-form solution (although it is substantially more complex than *F*). For increased accuracy, *Q/*2 could be used in place of *F* using this observed relationship.

#### 3.2. Polished surface comparison

Much of the detail for a similar analysis of *f _{μ}/f_{p}* in the polished surface approximation is contained in [13], where the microfacet model was compared to GHS. However, GHS was shown to reduce to MBK in the polished surface approximation [15]. Key points from that analysis are presented here.

The ratio *f _{μ}* /

*f*results in

_{p}Compared to the very rough surface approximation, the angular dependence of the exponential terms has changed, so the microfacet angular dependence is no longer equal to the wave optics angular dependence. From Eq. (19), instead of the exponential from scalar wave optics having an angular dependence of tan^{2} *θ _{h}*, the angular dependence is
${\eta}_{r}^{2}=\left({\mathrm{tan}}^{2}{\theta}_{h}\right){\left(\mathrm{cos}{\theta}_{i}+\mathrm{cos}{\theta}_{s}\right)}^{2}$. Unlike with the very rough surface, the polished surface angular term in the exponential from scalar wave optics,

*η*, only ranges from 0 to 4, not from 0 to ∞, but the microfacet model angular term, tan

_{r}^{2}

*θ*, still ranges from 0 to ∞. Additionally, a different relationship for

_{h}*σ*is observed, given as

_{g}When making this substitution for *σ _{g}*, it is observed that

*ρ*remains a unitless quantity, but varies with surface height and wavelength; that is, ${\rho}_{s}\propto {\sigma}_{s}^{2}/{\lambda}^{2}$. Additionally, by applying Eq. (22) to substitute for 2

_{s}*F/Q*, it becomes apparent that there remains a (cos

*θ*+ cos

_{i}*θ*)

_{s}^{−4}term in the ratio, suggesting the microfacet model would be multiplied by a (cos

*θ*+ cos

_{i}*θ*)

_{s}^{4}term.

Since there is now a wavelength dependence both in the distribution function and in *ρ _{s}*, it becomes apparent that a geometric optics model is not suited for this regime. Furthermore, the angular dependence is no longer simply tan

^{2}

*θ*, but now is

_{h}*η*. Equation (19) does relate the two expressions, but since the two are not equal, the microfacet model does not appear to be suitable for fitting to polished surfaces, unless corrections for the differences noted here were present.

_{r}#### 3.3. Original Beckmann-Kirchhoff modification

In the original BK theory,
${F}_{bk}^{2}$ was used instead of *K*, and was described as a geometrical term [27]. Although Krywonos used *K* in place of
${F}_{bk}^{2}$ in MBK,
${F}_{bk}^{2}$ is much simpler to compute and may be desirable in some applications of BRDF if a more complete understanding of its role were developed. That purpose is explored in this section. This geometrical factor
${F}_{bk}^{2}$ was given as [27]

*F _{bk}* is not to be confused with Fresnel reflectance; it represents a geometric term that does not depend on index of refraction. Its connection to the microfacet model was not clear from prior work on MBK in [27, 28]; in this section, we develop an explanation and a minor modification for this term. Rewriting
${F}_{bk}^{2}$ in microfacet coordinates

*θ*and

_{h}*θ*:

_{d}Using basic trigonometric manipulations to put *S* in microfacet coordinates, this term can be rewritten as

These equations suggest the purpose of
${F}_{bk}^{2}$ is to approximate *Q/*2 so that when BK theory is multiplied by the Fresnel equation *F*, it scales approximately correctly. However, there is an issue with
${F}_{bk}^{2}$ as presented. When *θ _{i}* ≠

*θ*, ${F}_{bk}^{2}$ does not obey Helmholtz reciprocity and differs from this approximation for

_{s}*Q/*2

*F*. For this reason, we propose using 1

*/S*rather than the original ${F}_{bk}^{2}$ term; using 1

*/S*makes the term functionally equivalent to the cross section conversion and microfacet distribution normalization, with the modified angular term that arises from

*g*in MBK.

#### 3.4. BRDF scaling by 1/σ

In [12], experimentally measured BRDF data was presented in direction cosine space, which aligns BRDF data at all incident angles, other than differing in height. We observed that multiplying measured BRDF data by a term inspired by the microfacet cross section conversion, cos*θ _{i}* cos

*θ*, resulted in an approximate height alignment. With the 1

_{s}*/S*term written in micro-facet coordinates above, the scaling in direction cosine space that was observed experimentally in [12] can be explained. When at the specular peak, cos

*θ*= 1 (

_{h}*θ*= 0). Near the specular peak,

_{h}*θ*is small and thus cos

_{h}*θ*1. Thus,

_{h}≈*S ≈*cos

*θ*cos

_{i}*θ*as was observed in [12] to be the scaling term that aligned the BRDF magnitude data in direction cosine space.

_{s}## 4. Conclusion

The terms of the MBK BRDF model are now compared to microfacet models, enabling the MBK model to be rewritten in microfacet coordinates (assuming Gaussian statistics) for a very rough surface (*g* > 800) as

*g*< 0.025),

The infinite summation can also be written in microfacet coordinates as

*g*is defined in Eq. (12). In this form, a direct comparison between the microfacet class of BRDF models given in Eq. (2) is now possible, which is expected to lead to more physical closed-form BRDF models for applications such as remote sensing and computer graphics.

To obtain the necessary relationships that were used to obtain the above form for the MBK model in microfacet coordinates, and to better comprehend what the microfacet model represents, the very rough surface approximation of MBK was compared to the general form of the microfacet model. In the rough surface approximation, the angular dependence of MBK (and thus also of GHS) was shown to be exactly equal to the angular dependence in the microfacet distribution function. The remaining angular terms were compared to the ratio of the standard Fresnel term with the polarization factor *Q*. This led to a novel approximation for *Q*, and suggested where certain terms present in the microfacet BRDF models arise. In particular, part of the microfacet distribution function normalization and the cross section conversion term were found to result from using standard unpolarized Fresnel reflectance *F* instead of the polarization factor *Q*. The approximation was found to be relatively accurate, except when *θ _{i}* or

*θ*is large; in the microfacet model, this region is where there would be significant shadowing and masking. Since the BRDF is fundamentally a physical optics problem, this observation suggests that instead of attempting to derive the shadowing and masking term from geometric optics, better results may be obtained by developing a modification to the ratio 2

_{s}*F/Q*, particularly at large incident or scattered angles. Alternatively, it is now possible to modify the microfacet model to use

*Q*instead of

*F*using this relationship.

Next, key differences between the polished surface analysis that were detailed in [13] and the rough surface analysis developed in this paper are presented. This resulted in a modification to the angular and wavelength dependencies in the microfacet model for a polished surface, and showed the breakdown of a geometric model for BRDF.

Additionally, the original BK geometric term ${F}_{bk}^{2}$ was analyzed, developing an explanation for this term relative to both the microfacet and linear systems models. A modification to ${F}_{bk}^{2}$ in the BK model was also proposed. The physical insight developed in this paper is expected to lead to future development of better closed-form approximations to the BRDF that preserve the relative simplicity of the microfacet model, while increasing the fidelity of the microfacet BRDF approximation.

## Acknowledgments

We would like to thank the United States Air Force Research Laboratory, Sensors Directorate, for supporting this research. The views expressed in this paper are those of the authors and do not necessarily reflect the official policy or position of the United States Air Force, the U. S. Department of Defense, or the U. S. government.

## References and links

**1. **J. R. Schott, *Fundamentals of Polarimetric Remote Sensing* (SPIE, 2009). [CrossRef]

**2. **M. T. Eismann, *Hyperspectral Remote Sensing* (SPIE, 2012). [CrossRef]

**3. **K. E. Torrance and E. M. Sparrow, “Theory of off-specular reflection from roughened surfaces,” J. Opt. Soc. Am. **57**, 1105–1114 (1967). [CrossRef]

**4. **J. R. Maxwell, J. Beard, S. Weiner, D. Ladd, and S. Ladd, “Bidirectional reflectance model validation and utilization,” Environmental Research Institute of Michigan (ERIM) Technical Report AFAL-TR-73-303 (1973).

**5. **D. R. Crow, C. F. Coker, D. L. Garbo, and E. M. Olson, “Closed-loop real-time infrared scene generator,” Proc. SPIE **3368**, 342 (1998). [CrossRef]

**6. **J. F. Blinn, “Models of light reflection for computer synthesized pictures,” in Proceedings of the 4th Annual Conference on Computer Graphics and Interactive Techniques, (ACM, 1977), pp. 192–198.

**7. **R. L. Cook and K. E. Torrance, “A reflectance model for computer graphics,” ACM Trans. Graphics **1**, 7–24 (1982). [CrossRef]

**8. **W. Matusik, H. Pfister, M. Brand, and L. McMillan, “A data-driven reflectance model,” ACM Trans. Graphics **22**, 759–769 (2003). [CrossRef]

**9. **A. Ngan, F. Durand, and W. Matusik, “Experimental analysis of BRDF models,” in Proceedings of the Euro-graphics Symposium on Rendering, (Eurographics Association) pp. 117–226 (2005).

**10. **F. E. Nicodemus, J. C. Richmond, J. J. Hsia, I. W. Ginsberg, and T. Limperis, “Geometrical considerations and nomenclature for reflectance,” *National Bureau of Standards Monograph 160*, Department of Commerce (1977).

**11. **S. D. Butler and M. A. Marciniak, “Robust categorization of microfacet BRDF models to enable flexible application-specific BRDF adaptation,” Proc. SPIE **9205**, 920506 (2014). [CrossRef]

**12. **S. D. Butler, S. E. Nauyoks, and M. A. Marciniak, “Experimental analysis of bidirectional reflectance distribution function cross section conversion term in direction cosine space,” Opt. Lett. **40**, 2445–2448 (2015). [CrossRef] [PubMed]

**13. **S. D. Butler, S. E. Nauyoks, and M. A. Marciniak, “Comparison of microfacet BRDF model elements to diffraction BRDF model elements,” Proc. SPIE **9472**, 94720C (2015). [CrossRef]

**14. **J. E. Harvey and A. Krywonos, “Unified scatter model for rough surfaces at large incident and scatter angles,” Proc. SPIE **6672**, 66720C (2007). [CrossRef]

**15. **A. Krywonos, J. E. Harvey, and N. Choi, “Linear systems formulation of scattering theory for rough surfaces with arbitrary incident and scattering angles,” J. Opt. Sci. Am. A **28**, 1121–1138 (2011). [CrossRef]

**16. **G. J. Ward, “Measuring and modeling anisotropic surfaces,” in Proceedings of SIGGRAPH Computer Graphics ’92, (ACM), 265–272. (1992).

**17. **A. Duer, “An improved normalization for the Ward reflectance model,” J. Graphics, GPU, and Game Tools **11**, 51–59 (2006). [CrossRef]

**18. **M. W. Hyde IV, J. D. Schmidt, and M. J. Havrilla, “A geometrical optics polarimetric bidirectional reflectance distribution function for dielectric and metallic surfaces,” Opt. Express **17**, 22138–22153 (2009). [CrossRef] [PubMed]

**19. **R. G. Priest and T. A. Germer, “Polarimetric BRDF in the microfacet model: Theory and measurements,” inProceedings of 2000 Meeting of the MSS Specialty Sensors Group on Passive Sensors, (Naval Research Lab, 2000), pp. 169–182.

**20. **J. D. Jackson, *Classical Electrodynamics*, 3rd ed. (John Wiley and Sons, Inc., 1999).

**21. **J. Dorsey, H. Rushmeier, and F. Sillion, *Digital Modeling of Material Appearance* (Morgan Kaufmann, 2007).

**22. **S. Rusinkiewicz, “A new change of variables for efficient BRDF representation,” in *Rendering Techniques ’98*, (Springer Vienna, 1998), pp. 11–22. [CrossRef]

**23. **R. G. Priest and S. R. Meier, “Polarimetric microfacet scattering theory with applications to absorptive and reflective surfaces,” Opt. Eng. **41**(5), 988–993 (2002). [CrossRef]

**24. **T. S. Trowbridge and K. P. Reitz, “Average irregularity reprentation of a rough surface for ray reflection,” J. Opt. Soc. Am. **65**, 531–536 (1975). [CrossRef]

**25. **T. M. Elfouhaily and C.-A. Guerin, “A critical survey of approximate scattering wave theories from random rough surfaces,” Waves in Random Media **14**, R1–R40 (2004). [CrossRef]

**26. **S. Chakrabarti, A. A. Maradudin, and E. R. Mendez, “Reconstruction of the surface-height autocorrelation function of a randomly rough dielectric surface from incoherent light scattering,” Phys. Rev. A **88**, 013812 (2013). [CrossRef]

**27. **P. Beckmann and A. Spizzichino, *The Scattering of Electromagnetic Waves from Rough Surfaces* (MacMillan, 1963).

**28. **J. E. Harvey, A. Krywonos, and C. L. Vernold, “Modified Beckmann-Kirchhoff scattering model for rough surfaces with large incident and scattering angles,” Opt. Eng. **46**(7), 078002 (2007). [CrossRef]

**29. **A. Krywonos, “Predicting surface scatter using a linear systems formulation of non-paraxial scalar diffraction,” Ph.D. Dissertation, University of Central Florida (2006).

**30. **J. C. Stover, *Optical Scattering: Measurement and Analysis*, 3rd ed. (SPIE, 2012).

**31. **J. C. Stover and J. E. Harvey, “Unified scatter model for rough surfaces at large incident and scatter angles,” Proc. SPIE **6672**, 66720B (2007). [CrossRef]

**32. **S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, and J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Express **19**, 9820–9835 (2011). [CrossRef] [PubMed]

**33. **E. Heitz, “Understanding the masking-shadowing function in microfacet-based BRDFs,” J. Comp. Graph. Tech. **3**, 32–91 (2014).