Seven-parameter statistical model for BRDF in the UV band

A new semi-empirical seven-parameter BRDF model is developed in the UV band using experimentally measured data. The model is based on the five-parameter model of Wu and the fourteen-parameter model of Renhorn and Boreman. Surface scatter, bulk scatter and retroreflection scatter are considered. An optimizing modeling method, the artificial immune network genetic algorithm, is used to fit the BRDF measurement data over a wide range of incident angles. The calculation time and accuracy of the fiveand seven-parameter models are compared. After fixing the seven parameters, the model can well describe scattering data in the UV band. ©2012 Optical Society of America OCIS codes: (290.1483) BSDF, BRDF, and BTDF; (290.5820) Scattering measurements; (290.5880) Scattering, rough surfaces. References and links 1. C. Lavigne, G. Durand, and A. Roblin, “Ultraviolet light propagation under low visibility atmospheric conditions and its application to aircraft landing aid,” Appl. Opt. 45(36), 9140–9150 (2006). 2. M. A. Velazco-Roa and S. N. Thennadil, “Estimation of complex refractive index of polydisperse particulate systems from multiple-scattered ultraviolet-visible-near-infrared measurements,” Appl. Opt. 46(18), 3730–3735 (2007). 3. M. Minnaert, “The reciprocity principle of linear photometry,” Astrophys. J. 93, 403–410 (1941). 4. J. Stover, Optical Scattering, Measurement and Analysis (SPIE Press, 1995). 5. C. L. Walthall, J. M. Norman, J. M. Welles, G. Campbell, and B. L. Blad, “Simple equation to approximate the bidirectional reflectance from vegetative canopies and bare soil surfaces,” Appl. Opt. 24(3), 383–387 (1985). 6. P. Beckman and A. Spizzichino, The Scattering of Electromagnetic Waves from Rough Surfaces (Pergamon, 1963). 7. J. L. Roujean, M. Leroy, and P. Y. Deschamps, “A bidirectional reflectance model of the earth’s surface for the correction of remote sensing data,” J. Geophys. Res. 97(20), 455–468 (1992). 8. X. Li and A. H. Strahler, “Geometric-optical bidirectional reflectance modeling of the discrete crown vegetation canopy: effect of crown shape and mutual shadowing,” IEEE Trans. Geosci. Rem. Sens. 30(2), 276–292 (1992). 9. I. G. E. Renhorn and G. D. Boreman, “Analytical fitting model for rough-surface BRDF,” Opt. Express 16(17), 12892–12898 (2008). 10. K. E. Torrance and E. M. Sparrow, “Theory for off-specular reflection from roughened surfaces,” J. Opt. Soc. Am. 57(9), 1105–1114 (1967). 11. Z. S. Wu, D. H. Xie, P. H. Xie, and Q. N. Wei, “Modeling reflectance function from rough surface and algorithms,” Acta Opt. Sin. 22, 897–901 (2002). 12. K. J. Dana, B. van Ginneken, S. K. Nayar, and J. J. Koenderink, “Reflectance and texture of real-world surfaces,” ACM Trans. Graph. 18(1), 1–34 (1999). 13. Y. Barnes and J. J. Hsia, “UV bidirectional reflectance distribution function measurements for diffusers,” Proc. SPIE 1764, 285–288 (1993). 14. M. P. Newell, L. A. Whitlock, and R. A. M. Keski-Kuha, “Extreme ultraviolet scatter from particulate contaminated mirrors,” Proc. SPIE 2541, 174–185 (1995). 15. M. P. Newell and R. A. M. Keski-Kuha, “Bidirectional reflectance distribution function of diffuse extreme ultraviolet scatterers and extreme ultraviolet baffle materials,” Appl. Opt. 36(22), 5471–5475 (1997). 16. T. Zurbuchen, P. Bochsler, and F. Scholze, “Reflection of ultraviolet light at 121.6 nm from rough surfaces,” Opt. Eng. 34(5), 1303–1315 (1995). 17. M. P. Newell and R. A. M. Keski-Kuha, “Extreme ultraviolet BRDF measurements: instrumentation and results,” Proc. SPIE 2864, 453–464 (1996). 18. C. Amra, “Light scattering from multilayer optics. I. Tools of investigation,” J. Opt. Soc. Am. A 11(1), 197–210 (1994). 19. C. Amra, “Light scattering from multilayer optics. II. Application to experiment,” J. Opt. Soc. Am. A 11(1), 211–226 (1994). #165850 $15.00 USD Received 30 Mar 2012; revised 24 Apr 2012; accepted 4 May 2012; published 11 May 2012 (C) 2012 OSA 21 May 2012 / Vol. 20, No. 11 / OPTICS EXPRESS 12085 20. C. Amra, “From light scattering to the microstructure of thin-film multilayers,” Appl. Opt. 32(28), 5481–5491 (1993). 21. J. M. Elson, J. P. Rahn, and J. M. Bennett, “Light scattering from multilayer optics: comparison of theory and experiment,” Appl. Opt. 19(5), 669–679 (1980). 22. C. Amra, C. Grèzes-Besset, and L. Bruel, “Comparison of surface and bulk scattering in optical multilayers,” Appl. Opt. 32(28), 5492–5503 (1993). 23. H. L. Zhang, Z. S. Wu, Y. H. Cao, and G. Zhang, “Measurement and statistical modeling of BRDF of various samples,” Opt. Appl. 40, 197–208 (2010).


Introduction
Interest in the study of scattering and reflective properties in the UV band (200 nm to 400 nm) has recently increased.The applications of such properties include UV space object detection and subsurface defect short wavelength monitoring in the semiconductor industry, among others [1,2].
The bidirectional reflectance distribution function (BRDF) is often used to describe the directional dependence of the scattering properties of a surface.The BRDF has been extensively studied and surveyed in various fields [3][4][5][6][7][8][9][10].Many BRDF models have been developed, and they can be classified into two categories, namely, purely empirical and analytical models.
Both models have their advantages and disadvantages.The purely empirical BRDF models are simple and useful but have no physical basis, such as the Minnaert BRDF model [3], Walthall model [5], etc.The analytical BRDF models are derived from more complex physical theory by simplifying some assumptions and approximations.They are needed in many applications, especially in simulation software, but they usually require the determination of many parameters.Examples include the Roujean [7], LiSparse-Dense BRDF [8], and Renhorn and Boreman 14-parameter [9] models.
Most BRDF models are developed over the visible and near-infrared spectrum.There are numerous papers that have discussed the measurement or model method in the visible to infrared band [7][8][9][10][11][12].Although there are some related papers reporting UV scatter in rough surfaces [13][14][15][16][17], a suitable BRDF model in the UV band is still lacking.In the UV band, the refraction light spectrum is often not easy to be measured, and we should consider more complex model to describe the measurement data.Bulk or volume scattering are usually caused by the in-homogeneities of materials.Relative studies about volume scattering provide a powerful tool to investigate the scattering properties of roughness surfaces [18][19][20][21][22].
The present paper aims to provide a model of reflection by a rough surface that can successfully predict the experimental findings in the UV band.In section one, the schematic diagram of the instrument used to performing scatter measurements is presented.In section two, a novel seven-parameter model is developed.There are three terms in this model: surface scatter, bulk scatter, and retro-reflection.An optimizing modeling method is used to model the BRDF measurement data of typical samples in the UV band.In the last section, the calculation results using the new seven-parameter model are compared with the fiveparameter model.

Measurement of BRDF in the UV band
A BRDF is actually an angle-resolved energy distribution.As shown in Fig. 1, an element of surface dA is illuminated by a source of incident wave with wave vector i k .Letting r k be the reflected direction wave vector, the symbol ẑ denotes the normal of the mean surface of the z axis and n denotes the normal direction of the micro-facet dA .Letting α be the angle between n and ẑ , where γ is the angle between i k and n , angles α and γ satisfy the following relationships ( ) where θ and φ are the zenith and azimuthal angles, respectively.The subscripts i and r denote the incident and reflected angles, respectively.The definition of BRDF, r f , can be expressed as the differential radiance dL r scattered by a uniformly illuminated, homogeneous material per unit differential incident irradiance dE i , , ,  , , , , Another equivalent definition of the BRDF is ( ) where s Ω is the solid angle through which the scattered power s P is collected.It is then normalized with respect to the total incident power i P .On the other hand, cos r θ is the factor that can be thought of as a correction to adjust the illuminated area to its apparent area when viewed from the direction of the scatter.Figure 2 shows the schematic diagram of the instrument used to perform angle-resolved optical scatter measurements.This BRDF measurement installation system was designed and constructed by the Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences.This measurement system has three parts, i.e., the source, sample holder, and receiver.The design details of each part are described elsewhere [23].
The UV light source used in these investigations was a 150 mW-power deuterium lamp, which can generate continuous emission spectra.After emerging from the exit slit, the beam is quasi-collimated by passing through a convex lens.After passing through two mirrors, the beam is reflected to the slit.By passing through the slit, the beam enters the sample and is then be collected by the detector.The system can measure three-dimensional BRDF data by rotating the sample and light source using electromotors A, B, and C in different threedimensional directions.The detector used was a spectrometer made by Ocean Optic Co. in USA.The sample can be rotated about a vertical axis to allow different angles of incidence, and the detector can be rotated independently about the same axis to allow the measurement of angular scatter distribution.The step of the zenith angle is 1°, and the step of the azimuth angle is 5°.Near the specular direction, the sample angle is 1° to guarantee sufficient precision.The accuracy of the zenith and azimuth angles is 0.1°, and the relative error of the measured BRDF of this system is less than 5%.The incident UV wavelength varies from 250.45288 nm to 368.89652 nm.The total 1024 output voltages s V for the sample are measured within this wavelength range.The voltages of reflectance standard plate ref V are also measured within the corresponding wavelengths.Then, the ratio is as follows ( , ; , , ) ( , ; , , ) where ( , ; , , ) f θ ϕ θ ϕ λ and ( , ; , , ) f θ ϕ θ ϕ λ are the BRDFs of the sample and reflectance standard plate at a given incident wavelengths λ , respectively.
Notably, V ref of the reference standard plate must be determined under the same measurement conditions used to detect the samples, and all measurements must be performed correspondingly.In our experiments, the reference standard plate is the pressed polytetrafluoroethylene plate, which can be considered as a perfect Lambertian plate.

Development of a seven-parameter BRDF model
Torrance and Sparrow [10] have deduced a useful light-scattering model, which assumes that the surface consists of small, randomly disposed mirror-like facets.An implied criterion of this model is that the root-mean-square (RMS) surface roughness is greater than the wavelength of the incident radiant energy.There are three parameters in the Torrance-Sparrow model, based on which Wu [11] has proposed a five-parameter BRDF model, whose expression can be expressed as ( , , , ) , , , cos exp 1 cos cos cos cos 1 1 cos where For isotropic materials at a given wavelength, the BRDF can be safely assumed to be invariant to azimuthal rotations of the incident and 0  UV wavelengths are shorter than visible and infrared bands.Thus, a new model that suited the experimental measurement data in the UV band is developed.Our model partly adopted the model of Wu and partly adopted the Renhorn-Boreman fourteen-parameter analytical The corresponding bulk scatter and retro-reflection terms are introduced to improve the accuracy of the simulation results of the five-parameter model.
According to the Renhorn-Boreman model [9], a Gaussian surface auto-covariance function ( , ) g x y is given by ( , ) exp ( ) where g σ is the RMS surface roughness and ρ is the inverse of the surface correlation length.Its Fourier transform corresponding to the surface power spectrum is ( , ) G ξ η .The parameters ξ and η are direction cosines defined as 0 sin( ) An exponential auto-covariance function is given by 2 ( , ) exp ( ) The Fourier transformation results of Eq. ( 9) can be obtained as follows using Mathematic under the condition 0 1, 0 1 In the Renhorn-Boreman model, the exponential statistical distribution is combined with the Gaussian distribution.Taking the angle of incidence into account, the Renhorn-Boreman model yields the two-dimensional Lorentzian BRDF by ( ) where σ is the integrated reflectance.
The Renhorn-Boreman BRDF model, which considers the surface scattering and includes the bulk and retro-reflection scattering sections, is partly adopted.The total BRDF is given by three separate parts as follows The surface scattering BRDF part is described by Eq. ( 6), whereas the bulk scattering and retro-reflection scattering BRDF are described as follows where 1 ρ and 2 ρ are the parameters related to the bulk and retro-reflection scatterings, respectively.Given that the proportions of the bulk and retro-reflection scattering sections are very small, only the masking and shadowing effects are considered in the surface scattering section.For the bulk and retro-reflection scattering sections, only the two-dimensional Lorentzian BRDF is considered.

Comparison of the seven-and five-parameter BRDF models
At this point, a new BRDF model has been developed.There are seven parameters to be determined in this model, namely, { } k k k a b ρ ρ .The artificial immune network algorithm is used to fit the parameters in the model.The optimization method is used to determine the optimum values of the fitting parameters of the BRDF to minimize the squared error between the measured data and the model.
( with the fitted BRDF derived from the seven-parameter (solid curves in blue) and five-parameter (short dash curves in red) models.The sample used in this simulation is a fine grinding rough aluminum plate.The zenith angles of incidence are 20°, 30°, 45°, and 60°.The scattering angle of scatter is over a range of angles from -70° to 70°.In the given wavelengths, six groups of measured data at various angles of incidence are usually used to fit the parameters to be determined (only four groups of measured data are shown in the Figs.[3][4][5][6].For convenience, the azimuthal angles of incidence and scatter are both set to 0°.The wavelengths of incidence are 266nm   Another example is shown in Figs. 5 and 6.Table 2 lists the model parameters, RMS error, and computing time of sample #2.The sample used in this case is a white painted surface.
The mirror-reflective section of sample #2 is more obvious than that of sample #1.Thus, sample #2 is relatively smoother than sample #1.The seven-parameter model still has less RMS error for different incidence wavelengths in this case.

Conclusion
A new seven-parameter BRDF model of rough surfaces in the UV band is established.An angle-resolved BRDF measurement system operates on some typical samples at different incident angles.The measured data are compared with the results of the seven-parameter model.Both for different wavelengths of the same sample or for different samples, the new seven-parameter model meets the experimental data better than the five-parameter one.It can significantly improve the calculation accuracy, although it costs a little more calculation time.
Both seven-and five-parameter model has its advantages and disadvantages.Optimization method will become more difficult when parameter in a model is so many.And the calculation time will be longer.If the precision requirement is not high, for most engineeringoriented application fields, the five-parameter model is a good enough model.But for some This paper offers a new developed BRDF model of rough surfaces in the UV band.The scattering properties of materials in the UV band are elucidated.And this kind of researches is helpful to understand the scattering properties of materials in the UV band.

iϕ=
. There are two components in Eq.(6); one is the specular reflection from mirror-like surface facets, and the other is a diffuse component.In the modified Torrance-Sparrow model of Wu[11], there are five parameters a , b , b k , d k , and r k to be determined.b k is the mirror-direction component, d k is related with the diffuse reflection component, and r k is related to the distribution of the #165850 -$15.00USD Received 30 Mar 2012; revised 24 Apr 2012; accepted 4 May 2012; published 11 May 2012 (C) 2012 OSA subsurface dA (determined by the slope distribution of the surface).
function of the sub-surface, exp[ (1 cos ) ] description of the Fresnel reflectance function, and ( ) ϕ λ is the masking and shadowing effect.In this model, the exponential function in the Torrance-Sparrow model is substituted by an elliptical function to describe the distribution of the normal of the facets.The exponential function with two parameters is used to substitute the Fresnel reflectance function to avoid the calculation of many trigonometric functions.The model can be used to describe isotropic surfaces with non-polarization incident light.The selection criteria of these five parameters are to minimize the RMS errors between the simulation and experimental data.
Fourier transform of the surface auto-covariance function is

Fig. 4 .(
Fig. 4. Comparison of the seven-and five-parameter models of sample #1 for various angles of incidence when wavelength of incidence is =369nm λ

Fig. 5 .-Fig. 6 .
Fig. 5. Comparison of the seven-and five-parameter models of sample 2# for various angles of incidence when wavelength of incidence is =266nm λ

#
165850 -$15.00USD Received 30 Mar 2012; revised 24 Apr 2012; accepted 4 May 2012; published 11 May 2012 (C) 2012 OSA special use such as in the UV band, it is difficult to detect the scattering signal.If you want to improve the calculation precision, more precision model should be used.
To ensure that the analytical two-dimensional integral of the BRDF is unitary, the parameter N gives #165850 -$15.00USD Received 30 Mar 2012; revised 24 Apr 2012; accepted 4 May 2012; published 11 May 2012 ) Measured BRDF data of sample #1: Fine grinding rough aluminum plate.