## 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 $Fbk2$ geometric term relative to both microfacet models and generalized Harvey-Shack model is presented. A modification to the geometric $Fbk2$ 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]

$fr=dLs(ω^i,ω^s,λ)dEi(ω^i,λ),$
where $ω^i$ is the incident unit vector and $ω^s$ is the scattered (outgoing) unit vector, where in both cases $ω^=[1,θ,ϕ]T$ is a unit vector pointing in the incident (or scattered) direction, with spherical coordinates θ (representing elevation angle) and ϕ (representing azimuthal angle) as measured from the surface normal; see Fig. 1(a). λ represents the wavelength of light.

Fig. 1 Relevant coordinate systems for the microfacet (geometric) BRDF model. (a) BRDF coordinate system orientation relative to the surface normal $z^$. (b) Microfacet coordinates relative to the coordinates defined in (a). In all cases, $ω^$ has a θ and ϕ component in spherical coordinates, and is normalized to a magnitude of 1 [11].

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]

$fμ(ω^i,ω^s)=ρsF(θd)Dμ(θh)G(ω^i,ω^s)σ(θi,θs)+ρvV(ω^i,ω^s)+ρd/π,$
where the first group of terms is due to surface reflection, the second to volumetric scattering V, and the third to Lambertian scattering. ρs, ρv, and ρd 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, ρv = 0 and ρd = 0. 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 + . 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
$σ(θi,θs)=14cosθicosθs.$

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 θh (representing the half angle, or angle of the bisector) and θd (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]

$ω^h=ω^i+ω^s||ω^i+ω^s||,ω^d=ℝy(−θh)ℝz(−ϕh),$
where ℝa is a rotation about some axis a. This formal definition has the advantage of specifying the ϕ dependence of $ω^h$ and $ω^d$ in addition to the θ dependence, which is of particular importance for anisotropic surfaces. In the isotropic case, only the θh and θd angles are computed using equations that appear in many papers in microfacet literature [3, 4, 7, 16, 18, 19, 23], and are given as
$cos(2θd)=cosθicosθs+sinθisinθscosϕs,$
$cosθh=cosθi+cosθs2cosθd.$

The Gaussian distribution Dg 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]

$Dg(θh)=12πσg2cos4θhexp[−tan2θh2σg2],$
where σg 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 tan2 θh. The remaining terms, and in particular the 1/cos4 θh 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,
$∫02π∫0π/2Dg(θh)cosθhsinθhdθhdϕh=1.$

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

$Δα=αs−αi=λνy=sinθssin(ϕs−π)Δβ=βs−βi=λνx=sinθscos(ϕs−π)−sinθi,$
where the extra π 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 ϕs = 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
$ηr2=(Δα)2+(Δβ)2=sin2θi+sin2θs+2sinθisinθscosϕs=(λνxy2π)2,$
where νxy is the notation used by Beckmann [27].

The MBK formulation presented in Eq. (11) is the Angle Spread Function (ASF), fa, 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 12 term present in MBK as shown in [29], the ASF is calculated as

$fa=πKlc2λ2exp(−g)∑m=1∞gmm!mexp(−νxy2lc24m),$
where K is a renormalization term that ensures conservation of energy, lc is the correlation length of the surface, and g is related to the surface height σs and given as [28]
$g(θi,θs)=(2πσsλ)2(cosθi+cosθs)2.$

When written in terms of scattered radiance Ls, K is defined as [15, 28, 29]

$K(βi)=∫−∞∞∫−∞∞Ls(αs,βs−βi)dβsdαs∫−11∫−1−αs21−αs2Ls(αs,βs−βi)dβsdαs.$

Ls ≥ 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(θi, θs) < 0.025, which results in less than 1% error compared to the infinite series summation when θi70°. By using Eqs. (10)(12) and multiplying by the polarization factor Q to convert the ASF to a BRDF, this polished approximation fp can be written as

$fp=4π3lc2σs2Q(cosθi+cosθs)2λ4exp[−(πlcηrλ)2],$
where exp(−g) 1 since σs << λ for a polished surface, and K ≈ 1 for a polished surface, as shown in [15].

At the other extreme, if a surface is very rough, the BRDF fvr as given in [29] can be rewritten using Eqs. (10)(12) as

$fvr=KQlc24πσs2(cosθi+cosθs)2exp[−(lc2σs)2(ηrcosθi+cosθs)2],$
where again we have multiplied by the polarization factor Q to convert from ASF to BRDF. In [29], Krywonos defines very rough as g(θi, θs) > 800, to result in less than 1% maximum error as compared to the infinite summation when θi70°.

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, Qs = Qss + Qsp. For p-polarized light incident on a surface, Qp = Qps + Qpp. For unpolarized light incident on a surface, Q = Qs + Qp. If θi = θs and ϕs = 180° (specular reflection), the equations for Qs and Qp reduce to the standard Fresnel equations Fs and Fp.

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, fp and fvr from MBK theory are considered to be accurate physical models in the limit of a polished or very rough surface, respectively.

## 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μ/fp and fμ/fvr 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.

#### 3.1. Very rough surface comparison

The ratio of microfacet BRDFs fμ to the very rough surface MBK BRDF fvr is

$fμfvr=(ρs2K)(2FGQ)((cosθi+cosθs)24cosθicosθscos4θh)(σs2lcσg)2×exp[−(tan2θh2σg2−lc24σs2ηr2(cosθi+cosθs)2)].$

The exponential term will be analyzed first. The angular dependence is contained in two terms: tan 2θh from the microfacet model and $ηr2/(cosθi+cosθs)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 θh in terms of tan θh:

$tan2θh=4cos2θd−(cosθi+cosθs)2(cosθi+cosθs)2.$

The denominator already matches the denominator of $ηr2/(cosθi+cosθs)2$, so we only need to show that the numerator equals $ηr2$ as given by Eq. (10). This is performed by use of the double angle formula for cos(2θd) and adding zero:

$4cos2θd−(cosθi+cosθs)2=2+2(2cos2θd−1)−(cosθi+cosθs)2=2+2cos(2θd)−(cosθi+cosθs)2=sin2θi+cos2θi+sin2θs+cos2θs+2cos(2θd)−(cosθi+cosθs)2=sin2θi+sin2θs+2sinθisinθscosϕs=ηr2,$
where in the last line Eq. (5) was used to obtain $ηr2$. Thus, we have shown that linear systems direction cosine space for a very rough surface is exactly equal to the geometric microfacet model bisector space:
$tan2θh=(ηrcosθi+cosθs)2.$

The angular term in the exponential from the scalar wave optics model, $ηr2/(cosθi+cosθs)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, tan2θ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 σg for a very rough surface using the remaining terms in the exponential, obtaining

$σg,vr=σs2lc.$

Recall σs represents the width of the Gaussian in the vertical (height) direction and lc 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 σg 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.

Eq. (16) has now simplified to

$fμfvr=(ρs2K)(2FGQ)((cosθi+cosθs)24cosθicosθscos4θh).$

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 = Qs + Qp), although he does present approximations for a few limiting cases [30]. To improve on this understanding, the polarization factor 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

$S=4cosθicosθscos4θh(cosθi+cosθs)2≈2FQ,$
where the factor of 2 arises from Q = Qs +Qp = Fs +Fp = 2F if θi =θs and ϕs = 180° (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 2F/Q and S, including at out-of-plane angles, where the relative difference is defined as

$Rd=|2F/Q−S|2F/Q.$

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 + 10i; low n and low (but nonzero) κ: ñ = 1.5 + i, low n and zero κ: ñ = 1.4, moderate n and κ: ñ = 1.7+5i, and n < 1: ñ = 0.25+3i. When θi = θs and ϕs = 180° (specular reflection), Rd = 0. In general, even out-of-plane, the error is relatively small if θi and θs are both small, but increases when either θi or θs is large. This result suggests that the cross section conversion term σ (θis) in Eq. (3) and the microfacet distribution normalization 1/cos4 θh in Eq. (7) both arise from using F instead of Q.

Fig. 2 Surface plot of relative difference Rd for five different indices of refraction: (a-d) ñ = 4+10i, (e-h) ñ = 1.5+i, (i-l) ñ = 1.4, (m-p) ñ = 1.7+5i, (q-t) ñ = 0.25+3i. For each index, four different incident angles are plotted: (a,e,i,m,q) θi = 15°; (b,f,j,n,r) θi = 30°; (c,g,k,o,s) θi = 45°; (d,h,l,p,t) θi = 60°. The plots are all symmetric about ϕs = 180°. The dotted black line at ϕs = 180° on the right side of each plot represents in-plane scatter. In each case, Rd = 0 when θi = θs and ϕs = 180°; Rd is generally small when θi and θs are small.

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 2FG/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μ/fp 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μ / fp results in

$fμfp=(ρsλ28π2σs2)(2FGQ)(14cosθicosθscos4θh(cosθi+cosθs)2)×(λ22π2σg2lc2)exp[−(tan2θh2σg2−π2lc2ηr2λ2)].$

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 tan2 θh, the angular dependence is $ηr2=(tan2θh)(cosθi+cosθs)2$. Unlike with the very rough surface, the polished surface angular term in the exponential from scalar wave optics, ηr, only ranges from 0 to 4, not from 0 to ∞, but the microfacet model angular term, tan2θh, still ranges from 0 to ∞. Additionally, a different relationship for σg is observed, given as

$σg,p=λπlc2,$
which varies with wavelength, unlike in Eq. (20). This shows the microsurface normal interpretation in the microfacet model does not hold for a polished surface, since the microsurface normal distribution should not vary with wavelength.

When making this substitution for σg, it is observed that ρs remains a unitless quantity, but varies with surface height and wavelength; that is, $ρs∝σs2/λ2$. Additionally, by applying Eq. (22) to substitute for 2F/Q, it becomes apparent that there remains a (cosθi + cosθs)−4 term in the ratio, suggesting the microfacet model would be multiplied by a (cosθi + cosθ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 tan2θh, but now is ηr. 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.

#### 3.3. Original Beckmann-Kirchhoff modification

In the original BK theory, $Fbk2$ was used instead of K, and was described as a geometrical term [27]. Although Krywonos used K in place of $Fbk2$ in MBK, $Fbk2$ 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 $Fbk2$ was given as [27]

$Fbk2=(1+(cosθicosθs+sinθisinθscosϕs)cosθi(cosθi+cosθs))2.$

Fbk 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 $Fbk2$ in microfacet coordinates θh and θd:

$Fbk2=cos2θdcos2θicos2θh.$

Using basic trigonometric manipulations to put S in microfacet coordinates, this term can be rewritten as

$1S=cos2θdcosθicosθscos2θh≈Q2F.$

These equations suggest the purpose of $Fbk2$ 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 $Fbk2$ as presented. When θiθs, $Fbk2$ does not obey Helmholtz reciprocity and differs from this approximation for Q/2F. For this reason, we propose using 1/S rather than the original $Fbk2$ 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θs, resulted in an approximate height alignment. With the 1/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θh = 1 (θh = 0). Near the specular peak, θh is small and thus cosθh 1. Thus, S ≈ cosθi cosθs as was observed in [12] to be the scaling term that aligned the BRDF magnitude data in direction cosine space.

## 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

$fvr=KFG4cosθicosθs(lc22πσs2cos4θh)exp[−(lc2σs)2tan2θh],$
and for a polished surface (g < 0.025),
$fp=(2π)3FG(cosθi+cosθs)44cosθicosθscos4θh(σs2lc2λ4)×exp[−(πlcλ)2tan2θh(cosθi+cosθs)2].$

The infinite summation can also be written in microfacet coordinates as

$fMBK=πFG(cosθi+cosθs)2Klc22λ2cosθicosθscos4θhexp(−g)×∑m=1∞gmm!mexp(−lc24mλ2tan2θh(cosθi+cosθs)2),$
where 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 θs 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 2F/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 $Fbk2$ was analyzed, developing an explanation for this term relative to both the microfacet and linear systems models. A modification to $Fbk2$ 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.

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).

### References

• View by:
• |
• |
• |

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, 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).

#### 2015 (2)

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]

#### 2014 (2)

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]

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

#### 2013 (1)

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]

#### 2011 (2)

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]

#### 2007 (3)

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]

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]

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

#### 2006 (1)

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

#### 2004 (1)

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]

#### 2003 (1)

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

#### 2002 (1)

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]

#### 1998 (1)

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]

#### 1982 (1)

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

#### Beard, J.

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).

#### Beckmann, P.

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

#### Blinn, J. F.

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.

#### Brand, M.

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

#### Butler, S. D.

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]

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]

#### Chakrabarti, S.

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]

#### Choi, N.

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]

#### Coker, C. F.

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]

#### Cook, R. L.

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

#### Crow, D. R.

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]

#### Dorsey, J.

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

#### Duer, A.

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

#### Durand, F.

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).

#### Eismann, M. T.

M. T. Eismann, Hyperspectral Remote Sensing (SPIE, 2012).
[Crossref]

#### Elfouhaily, T. M.

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]

#### Garbo, D. L.

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]

#### Germer, T. A.

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.

#### Ginsberg, I. W.

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).

#### Guerin, C.-A.

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]

#### Harvey, J. E.

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]

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

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]

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]

#### Heitz, E.

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

#### Hsia, J. J.

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).

#### Jackson, J. D.

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

#### Krywonos, A.

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]

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

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]

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

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).

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).

#### Limperis, T.

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).

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]

#### Marciniak, M. A.

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]

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]

#### Matusik, W.

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

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).

#### Maxwell, J. R.

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).

#### McMillan, L.

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

#### Meier, S. R.

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]

#### Mendez, E. R.

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]

#### Nauyoks, S. E.

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]

#### Ngan, A.

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).

#### Nicodemus, F. E.

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).

#### Olson, E. M.

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]

#### Pfister, H.

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

#### Priest, R. G.

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]

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.

#### Richmond, J. C.

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).

#### Rushmeier, H.

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

#### Rusinkiewicz, S.

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

#### Schott, J. R.

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

#### Sillion, F.

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

#### Spizzichino, A.

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

#### Stover, J. C.

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]

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

#### Torrance, K. E.

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

#### Vernold, C. L.

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]

#### Ward, G. J.

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

#### Weiner, S.

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).

#### ACM Trans. Graphics (2)

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

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

#### J. Comp. Graph. Tech. (1)

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

#### J. Graphics, GPU, and Game Tools (1)

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

#### J. Opt. Sci. Am. A (1)

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]

#### Opt. Eng. (2)

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]

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]

#### Phys. Rev. A (1)

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]

#### Proc. SPIE (5)

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]

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]

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]

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

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]

#### Waves in Random Media (1)

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]

#### Other (14)

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

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

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

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

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.

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.

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

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

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

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).

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

M. T. Eismann, Hyperspectral Remote Sensing (SPIE, 2012).
[Crossref]

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).

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).

### Cited By

OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

### Figures (2)

Fig. 1 Relevant coordinate systems for the microfacet (geometric) BRDF model. (a) BRDF coordinate system orientation relative to the surface normal $z ^$. (b) Microfacet coordinates relative to the coordinates defined in (a). In all cases, $ω ^$ has a θ and ϕ component in spherical coordinates, and is normalized to a magnitude of 1 [11].
Fig. 2 Surface plot of relative difference Rd for five different indices of refraction: (a-d) ñ = 4+10i, (e-h) ñ = 1.5+i, (i-l) ñ = 1.4, (m-p) ñ = 1.7+5i, (q-t) ñ = 0.25+3i. For each index, four different incident angles are plotted: (a,e,i,m,q) θi = 15°; (b,f,j,n,r) θi = 30°; (c,g,k,o,s) θi = 45°; (d,h,l,p,t) θi = 60°. The plots are all symmetric about ϕs = 180°. The dotted black line at ϕs = 180° on the right side of each plot represents in-plane scatter. In each case, Rd = 0 when θi = θs and ϕs = 180°; Rd is generally small when θi and θs are small.

### Equations (31)

$f r = d L s ( ω ^ i , ω ^ s , λ ) d E i ( ω ^ i , λ ) ,$
$f μ ( ω ^ i , ω ^ s ) = ρ s F ( θ d ) D μ ( θ h ) G ( ω ^ i , ω ^ s ) σ ( θ i , θ s ) + ρ v V ( ω ^ i , ω ^ s ) + ρ d / π ,$
$σ ( θ i , θ s ) = 1 4 cos θ i cos θ s .$
$ω ^ h = ω ^ i + ω ^ s | | ω ^ i + ω ^ s | | , ω ^ d = ℝ y ( − θ h ) ℝ z ( − ϕ h ) ,$
$cos ( 2 θ d ) = cos θ i cos θ s + sin θ i sin θ s cos ϕ s ,$
$cos θ h = cos θ i + cos θ s 2 cos θ d .$
$D g ( θ h ) = 1 2 π σ g 2 cos 4 θ h exp [ − tan 2 θ h 2 σ g 2 ] ,$
$∫ 0 2 π ∫ 0 π / 2 D g ( θ h ) cos θ h sin θ h d θ h d ϕ h = 1 .$
$Δ α = α s − α i = λ ν y = sin θ s sin ( ϕ s − π ) Δ β = β s − β i = λ ν x = sin θ s cos ( ϕ s − π ) − sin θ i ,$
$η r 2 = ( Δ α ) 2 + ( Δ β ) 2 = sin 2 θ i + sin 2 θ s + 2 sin θ i sin θ s cos ϕ s = ( λ ν x y 2 π ) 2 ,$
$f a = π K l c 2 λ 2 exp ( − g ) ∑ m = 1 ∞ g m m ! m exp ( − ν x y 2 l c 2 4 m ) ,$
$g ( θ i , θ s ) = ( 2 π σ s λ ) 2 ( cos θ i + cos θ s ) 2 .$
$K ( β i ) = ∫ − ∞ ∞ ∫ − ∞ ∞ L s ( α s , β s − β i ) d β s d α s ∫ − 1 1 ∫ − 1 − α s 2 1 − α s 2 L s ( α s , β s − β i ) d β s d α s .$
$f p = 4 π 3 l c 2 σ s 2 Q ( cos θ i + cos θ s ) 2 λ 4 exp [ − ( π l c η r λ ) 2 ] ,$
$f v r = K Q l c 2 4 π σ s 2 ( cos θ i + cos θ s ) 2 exp [ − ( l c 2 σ s ) 2 ( η r cos θ i + cos θ s ) 2 ] ,$
$f μ f v r = ( ρ s 2 K ) ( 2 F G Q ) ( ( cos θ i + cos θ s ) 2 4 cos θ i cos θ s cos 4 θ h ) ( σ s 2 l c σ g ) 2 × exp [ − ( tan 2 θ h 2 σ g 2 − l c 2 4 σ s 2 η r 2 ( cos θ i + cos θ s ) 2 ) ] .$
$tan 2 θ h = 4 cos 2 θ d − ( cos θ i + cos θ s ) 2 ( cos θ i + cos θ s ) 2 .$
$4 cos 2 θ d − ( cos θ i + cos θ s ) 2 = 2 + 2 ( 2 cos 2 θ d − 1 ) − ( cos θ i + cos θ s ) 2 = 2 + 2 cos ( 2 θ d ) − ( cos θ i + cos θ s ) 2 = sin 2 θ i + cos 2 θ i + sin 2 θ s + cos 2 θ s + 2 cos ( 2 θ d ) − ( cos θ i + cos θ s ) 2 = sin 2 θ i + sin 2 θ s + 2 sin θ i sin θ s cos ϕ s = η r 2 ,$
$tan 2 θ h = ( η r cos θ i + cos θ s ) 2 .$
$σ g , v r = σ s 2 l c .$
$f μ f v r = ( ρ s 2 K ) ( 2 F G Q ) ( ( cos θ i + cos θ s ) 2 4 cos θ i cos θ s cos 4 θ h ) .$
$S = 4 cos θ i cos θ s cos 4 θ h ( cos θ i + cos θ s ) 2 ≈ 2 F Q ,$
$R d = | 2 F / Q − S | 2 F / Q .$
$f μ f p = ( ρ s λ 2 8 π 2 σ s 2 ) ( 2 F G Q ) ( 1 4 cos θ i cos θ s cos 4 θ h ( cos θ i + cos θ s ) 2 ) × ( λ 2 2 π 2 σ g 2 l c 2 ) exp [ − ( tan 2 θ h 2 σ g 2 − π 2 l c 2 η r 2 λ 2 ) ] .$
$σ g , p = λ π l c 2 ,$
$F b k 2 = ( 1 + ( cos θ i cos θ s + sin θ i sin θ s cos ϕ s ) cos θ i ( cos θ i + cos θ s ) ) 2 .$
$F b k 2 = cos 2 θ d cos 2 θ i cos 2 θ h .$
$1 S = cos 2 θ d cos θ i cos θ s cos 2 θ h ≈ Q 2 F .$
$f v r = K F G 4 cos θ i cos θ s ( l c 2 2 π σ s 2 cos 4 θ h ) exp [ − ( l c 2 σ s ) 2 tan 2 θ h ] ,$
$f p = ( 2 π ) 3 F G ( cos θ i + cos θ s ) 4 4 cos θ i cos θ s cos 4 θ h ( σ s 2 l c 2 λ 4 ) × exp [ − ( π l c λ ) 2 tan 2 θ h ( cos θ i + cos θ s ) 2 ] .$
$f M B K = π F G ( cos θ i + cos θ s ) 2 K l c 2 2 λ 2 cos θ i cos θ s cos 4 θ h exp ( − g ) × ∑ m = 1 ∞ g m m ! m exp ( − l c 2 4 m λ 2 tan 2 θ h ( cos θ i + cos θ s ) 2 ) ,$