Improvement of coated aluminum sheet pBRDF model based on scattering and phase function optimization

The pBRDF model is able to relate the properties of target materials to the polarization information of incident and reflected light, and is an important basis for obtaining polarization information of targets in space. It is an important basis for obtaining target polarization information and polarization detection of space targets. P-G model is the first strictly pBRDF model officially released, but there are still deficiencies. In this paper, we first analyze the assumption framework of the P-G model, derive the imperfections in the framework through the analysis of the assumption framework, and add scattering and phase function to enhance the existing model. On the basis of the existing P-G model and parameter inversion, the output results of the model are compared with the experimental data through simulation, and the results show that the relative error of the target's linear polarizability is reduced under the improved model, which proves the accuracy and precision of the improved model.


PG model
The material rough surface is combination of the infinity specular micro-mirror, and the material interior is the scattering of atoms and molecules, Fig. 1 is a description for them.When the light hits on the material surface, three types of photons are generated.The first type is the photon that directly reflects on the current mirror surface, as shown in A. This type of photon obeys the Fresnel reflection law and generates polarization.The second type is the photon that selectively absorbed and scattered by material atoms and molecules and reappear on the current surface, as shown in B. Those are considered to be unpolarized.The third type is considered as unpolarized.These photons which are not emitted from the current surface have experienced reflections or transmissions more than 2 times, as shown in C.
Based on the assumptions, the first type of photons needs to solve the S-wave and P-wave components of specular reflected light or solve the Mueller matrix.At the same time, since the surface is a micro-mirror combination, the orientation of the mirror normal will form a statistical distribution on the macroscopic level.It is necessary to establish the probability distribution function of the surface normal of the specular element.Since there will be obscuration and shadow between the mirror and mirror, it is also important to construct the obscuration and shadowing function.For the second type of photons, due to the limited number of collisions, the scattered light has different strengths in different directions, that is, volume scattering occurs.As to the third type of photons, no obvious directionality exists due to the enough collisions.It is equivalent to the direct-current background which is called diffuse scattering.
The PG model 7 based on the assumptions includes specular reflection, diffuse scattering and volume scattering, as shown in Eq. 1. (1)

Model improvement
In the field of scientific research, the integration of the scattering and phase functions in the polarized bidirectional reflectance distribution function (pBRDF) is fundamental to the comprehensive simulation of the complex interactions of light with material surfaces.pBRDF, a widely used method in computer graphics and image rendering, provides a sophisticated method for simulating realistic illumination and material reflection properties.The symbiotic application of scattering and phase functions is of critical importance by promoting a holistic understanding of how incident light interacts with and scatters from material surfaces.Scattering functions, including the normal distribution function (NDF) and the micro surface shading function, depict the distribution of micro surface normal and account for light occlusion due to microgeometry.The phase function takes into account the light transmission from the light source to the camera and considers the surface geometry and the spatial relationship between the light source and the camera.The synergistic utilization of the scattering and phase functions in the pBRDF accurately describes the process of light propagation and interaction on the surface of the material, which improves the realism and fidelity of the final rendering.
Butler and Nauyoks et al. analyzed both the linear optical model based on the assumption of specular microelement, and the diffractive optical model based on the assumption of surface continuous cosine.The results show that there are some differences between the two 16 .The new model combines the basic methods of linear optics and diffractive optics.
Assuming that the particles are irregular in shape and have no statistically significant directionality, the rough surface is the combination of an infinite number of hemispherical particles, shown as Fig. 2. For the reflected light, at the medium interface, there is a phase shift along the cross-section direction, and it will penetrate a certain depth in the direction perpendicular to the interface, that is, an evanescent wave.In the near field of surface fine structures, the evanescent waves effect cannot be ignored.
The evanescent wave penetration depth is a function of the angle of incidence and the complex refractive index.Defined as the position where the amplitude drops to 1/e.If there exist a utility surface, it represents a virtual reflection surface, and the influence of phase shift and penetration depth on the light intensity distribution is ignored.Since the solution of the electric field on the fine structure is still reasonable, the reflection of the material surface can be regarded as spherical particle scattering according to the symmetry of the reflection.At this time, the material surface reflected light can be equivalent to the mirror light reflection from the inside, and the scattering is generated.
The material scattering inside can be considered as a combination of shallow scattering and deep scattering.Unlike class B photons type in the PG model, shallow scattering produces polarization.For there are enough sufficient photon collisions, the deep scattering is unpolarized, which is similar to the PG model.
Compared with the micro-mirror combination, there is no clear reflection between the internal atoms and molecules and the external space on the surface of a real material, which is essentially the result of the internal atoms and molecules arrangement on the surface.Therefore, it can be modeled with a unified structural assumption.
The sphere of the material surface is approximated as a tile, so it can be equivalent to a scatter body.The scattering effect is a combination of single-particle scattering, and the equivalent scattering effect can be seen www.nature.com/scientificreports/as the Henyey-Greenstein (HG) phase function 17 The principal optical axis is equivalent to the material normal direction.The structure of the HG phase function is shown as Eq. 2.
where is the asymmetry factor.In the equation, θ r is the included angle of the reflected light deviating from the principal optical axis.
The HG phase function is the intensity at the angle of the main optical axis.Its polarization still needs the Mueller matrix, as shown in Eq. 7.
In Eq. 4, R is the reflectivity, Ψ is the tangent angle formed by the ratio of the amplitudes of the p-wave and the s-wave, and Δ is the phase difference between the p-wave and the s-wave, all of which can be solved by the Fresnel equation.
The phase function P represents the normalized probability distribution of the reflected light intensity for the reflection, so the bidirectional reflection distribution of the reflection intensity can be expressed as Eq. 5.
For the interior scattering, an assumption of a thick isotropic non-absorbing medium is put forward in the model.In the upper hemisphere space of the material, the bidirectional reflection distribution contains two parts, shallow and deep scattering.For shallow scattering, the Mueller matrix can be directly replaced by the coefficient, since most matter is non-transparent and scattering is a smaller magnitude than reflection, as shown in Eq. 6.
In the equation, ρ s is the shallow scattering coefficient, which is equivalent to the volume scatter coefficient, g s is the asymmetric factor of scattering, θ i and θ r are the pitch angles of the incident and emission light, θ s is the angle between the scattered light and the main optical axis.
The deep scattering is theoretically expressed by Eq. 7. Since this value is also small in most cases, it may be equal to ρ d Then, the intensity BRDF can be expressed as Eq. 8.Only shallow scattering can produce polarization.For the degree of polarization, the proportional coefficient η can be set to calculate the degree of polarization of the shallow layer, since only part of the shallow scattering can produce polarized light, The overall linear polarization is shown in Eq. 9.
In this equation, Ψ is the tangent angle by the amplitude ratio of the p-wave and the s-wave.Compared with the PG model, the new model is mainly reflected in the following two aspects.A new distribution function replaces the original probability distribution function of the surface normal.The parameter B n is theoretically 1 but materially not 1 is removed, the obscuration and shadowing function removed as well.Shallow scatter model is refined.It is possible to generate partially polarization, and the intensity can increase (or decrease) as the phase increases, which is consistent with the actual situation.

Experimental verification and results
The light source used in the experiment was a 200W tungsten bromine lamp.The spectrometer is the USB2000 + XR of Ocean Optics, with a spectral range of 200-1050 nm and a spectral resolution of 1.7 nm.The fiber optic head is installed with a beam expanding optical path with a diameter of about 10 mm, and the polarization state can be changed by rotating the polarizer.The positional accuracy of the azimuth and pitch angle of the (5) www.nature.com/scientificreports/reflectance test carousel is better than ± 0.02°, and the platform can be adjusted up and down horizontally, and the measurement error can be controlled within 3% by the above experimental equipment's.The experimental equipment is shown in Fig. 3.
The test target was selected to be dark green coated aluminum plate, a color widely used to camouflage military armor materials.Dark green provides effective camouflage in many environments, allowing armored vehicles to better blend into natural environments such as forests, grasslands, and jungles, and reducing the likelihood of detection by the enemy.This is vital for military operations and concealment.Since the coating material is a mixture of complex material composition, the coating surface roughness is more likely to develop into scattering effect, so it is necessary to improve the polarization characteristic model of this coating material based on the scattering effect.The coating material is polymer resin, the thickness of the sample is 1 mm, and the thickness of the coating is 0.1 mm. Figure 4 shows the dark green coated aluminum plate used in the experiment.The general process of studying the pBRDF of rough surface of base material is as follows: measure the scattering characteristics of the base material under the irradiation of white light source, obtain the data of scattering characteristics of the base material, and then process the test data of the base material under the white light source to obtain the data of the corresponding polarization characteristics of the base material at other different observation azimuths and zenith angles, and establish the pBRDF model.Experiments are conducted on the experimental samples selected in this paper, firstly looking at the experimental test results of the light green paint surface.The same as the standard whiteboard experiments, the light source using tungsten bromine lamp white light source incidence, incident zenith angle of 20°, the detector using the lens before the addition of 400-675 nm filter SALSA camera, in the relative azimuth angle of 180° in the plane of the detection, the number of steps is 20°.The test results are shown in Fig. 5.
Figure 5 shows the test results of the green paint surface in the plane of 180° relative azimuth angle, specular reflection direction, i.e., when the observation zenith angle is 20°.From the figure can be initially understood: green paint surface in the direction of specular reflection and its vicinity, the reflection intensity is larger, and the intensity to both sides gradually decreasing; green paint surface backward scattering is not obvious, there is no more obvious non-specular reflection peak; green paint surface scattered light polarization with the change in the observation angle changes, and the magnitude of the change is larger, initially shows that the polarization detection, there is the best detection angle.
The test data of the total energy scattered from the surface of each experimental sample were normalized, and the results are shown in Figs. 6 and 7.
As can be seen from the Figs.6 and 7, the energy scattering characteristics of the coating material has a more consistent pattern of change, that is, in the direction of specular reflection there is a reflection peak in the other direction of scattering energy scattering energy gradually decreases, until a more uniform DC component intuitively illustrates the surface scattering contains at least two parts of the scattering behavior, the specular reflection of the microfacial elements and the surface of the many times the scattering The specular reflection of the microfacets and the multiple scattering from the surface.The specular reflection of microfacets is the main source of polarization information.
Figure 8 shows the polarizability plots of the PG model and the new model for dark green coated aluminum plates with wavelengths from 400 to 675 nm.The least squares method was used to fit the models to minimize the difference between the predicted values of the models and the actual observed values A total of 15 spectral curves of 440 nm, 500 nm, 550 nm, 600 nm, and 675 nm were selected and the statistics of the error results after parameter inversion were carried out, and the results are shown in the Table 1.

Discussion
The traditional pBRDF model is improved based on the scattering and phase functions, and the target parameters are inverted by the least squares method, and the inversion results are substituted into the model, which results in better fitting of the simulation curves and measurements.By comparing the relative errors of line polarization between the improved model and the traditional P-G model, it is concluded that the relative errors of line polarization of dark green coated aluminum plate are reduced by 77.31%, 3%, 60.17%, 11.6% and 15.63% under the improved model, respectively.It is proved that the proposed pBRDF model has better accuracy.In this paper, the internal intrinsic properties affecting the polarization characteristics of the target are investigated through

Figure 2 .
Figure 2. The basic assumption of the new model.

Figure 5 .
Figure 5. Test results of pBRDF in the direction of specular reflection on the surface of green paint at an angle of incidence of 20°.

Figure 6 .
Figure 6.Variation curves of strength of dark green paint-coated aluminum plate with detection zenith angle in different bands.

Figure 7 .
Figure 7. Variation curves of polarizability of dark green paint-coated aluminum plate with detection zenith angle in different bands.

Figure 8 .
Figure 8. Polarizability plots of PG model and new model of dark green coated aluminum plate from 400 to 675 nm.

Table 1 .
Comparison of errors between the PG model and the new model.