Next Article in Journal
Inferring Changes in Arctic Sea Ice through a Spatio-Temporal Logistic Autoregression Fitted to Remote-Sensing Data
Previous Article in Journal
A Review of Artificial Intelligence and Remote Sensing for Archaeological Research
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Degree of Polarization Calculation for Laser Backscattering from Typical Geometric Rough Surfaces at Long Distance

1
State Key Laboratory of Pulsed Power Laser Technology, National University of Defense Technology, Hefei 230000, China
2
Anhui Province Key Laboratory of Electronic Restriction, National University of Defense Technology, Hefei 230037, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(23), 6001; https://doi.org/10.3390/rs14236001
Submission received: 20 October 2022 / Revised: 19 November 2022 / Accepted: 22 November 2022 / Published: 26 November 2022
(This article belongs to the Section Engineering Remote Sensing)

Abstract

:
Measurement of the degree of polarization of backscattering light from rough surfaces plays an important role in targets-detection applications. The polarization bidirectional reflectance function is the key to establish the relation between the polarization states of incidence and backscattering light. For the purpose of obtaining a polarized bidirectional reflectance distribution function (pBRDF) of a realistic, complicated target, it is decomposed as typical geometric surfaces and analytically calculated as the degree of polarization of the backscattering light, using a microfacet model, under conditions in which the scale of the target is far less than the target distance. In an experiment testing several typical geometric models, the results coincided with the theoretical calculation. The degree of polarization varied substantially as the rotation angle of the target changed, but showed little dependence on the size of target. The results have potential in applications discriminating between targets at different spatial orientations.

1. Introduction

Rough surfaces scattering light can carry varied information about the target characteristics, such as the surface material, the spatial form, or the movement status. Measuring the parameters, such as the echo waveform, the intensity, the phase, or the polarization states of scattering light to extract these target characteristics facilitates varied application in the realm of target discrimination [1,2,3,4,5], vegetation detection [6,7], the earth surface investigation [8,9,10], ocean water monitoring [11,12,13], atmospheric remote sensing [14,15], light detection and ranging (LiDAR) [16,17,18,19,20], hyperspectral data analysis [21] and topography analysis [22].
Polarization is the mode of spatial orientations in the electric field. An arbitrary polarization state can be described by four Stokes parameters. More explicitly, the Stokes parameter S 0 represents the total light intensity, while S 1 , S 2 and S 3 represent the intensity difference between horizontal and vertical polarized light, 45 and 45 polarized light, and right and left circular polarized light, respectively. Polarization detection can obtain these Stokes parameters as well as the degree of polarization (DOP) [13,23], or an angle of polarization (AOP) [24,25]. It has the following major advantages against intensity detection. Firstly, it can discriminate between different surface materials since the surfaces with different roughness parameters or refractive indices have significant differences in DOP [3]. Secondly, atmospheric turbulence affects polarization much less than the light intensity or the phase in question. Therefore, polarization detection has higher visibility and thus a longer imaging distance in foggy or hazed atmospheres. The polarization detection can thus be used in dehazing imaging technology [25,26].
The theoretical model for rough surface scattering is generally based on the bidirectional reflectance distribution function (BRDF) [27,28,29], which describes the distribution characteristics of scattered light in the upper hemisphere space of a rough surface. Such a function includes many complex physical processes such as specular reflection, diffuse reflection, bulk scattering, and others involving many parameters that affect the characteristic of scattering light, making the exact model for rough surfaces scattering impractical in the analysis or calculations. Various approximate models are proposed to simplify the numerical calculation in practical applications according to the scale of rough surfaces parameters such as the standard deviation of roughness σ h and correlation length l. When wavelength λ σ h , l or λ σ h , l , the Kirchhoff–Beckmann’s approach [30,31] or perturbation method [32], respectively, are suitable approximation theories. The most practical one for objective scattering measurement is the Torrance–Sparrow model, also known as microfacet theory [28,33], treating macro rough surfaces as innumerable small perfectly smooth facets, the size of which is sufficiently large compared with the wavelength, such that reflection on them obeys Fresnel’s law. The normal directions of microfacets are random, but are usually assumed to obey normal distribution.
On the other hand, the polarization characteristic of scattering light is described by the polarization BRDF (pBRDF), connecting the Stokes vectors of incident and scattering light, this idea was firstly proposed by D. Flynn [34]. G. Priest [35] extended the unpolarized BRDF to the polarization version using the Mueller matrix, giving an explicit method to obtain the Stokes vector of backscattered light from rough surfaces. Priest’s model considered only the specular component of scattered light, which failed to describe the depolarization effect. Various modified models are proposed to take the contribution of diffuse scattering into account [23,36,37]. M. Hyde introduced the directional hemispherical reflectance function and treated the rough surfaces as ideal Lambertian [36], which approximately coincides the experiment result. I. Renhorn presented a four-parameter pBRDF model, accurately fitting the s and p components of polarization [38,39]. H. Liu proposed another empirical model, considering not only the reflectance, but also refraction on rough surfaces, to describe the subsurface scattering according to Fresnel’s law [40].
Measuring the polarization states and DOP from them is necessary for the remote polarization detection of realistic targets. However, in application, the current literature is mainly focused on polarization imaging [41] or discrimination between different kinds of surface materials [42,43] based on the degree of polarization. In the theoretical research area, pBRDF models are concerned with dependence on the incidence or reflection angles [36,44,45,46]. This is meaningful and essential in improving the accuracy of pBRDF. However, pBRDF models rely on so many variables, making it difficult to perform calculations in the practical applications. For example, for a given satellite model, the general method for obtaining the pBRDF is measuring the polarization states of backscattering light at different incidence and reflection angles. This is not feasible when the target is a long-distance target.
Therefore, it is necessary to give a numerical or even analytical calculation for the pBRDF of targets. The numerical calculation of the Stokes vector or Mueller matrix could be complicated due to the complication of the geometry of the target surfaces. However, the irregular objectives can be simplified to be the combination of several typical regular geometric surfaces, which easily calculate the backscattering Stokes vector and thus the DOP. For example, Figure 1 shows a model of LIDAR imaging satellite, which consists of typical geometric shapes: two rectangle solar panels (blue parts), two cuboid or prism main bodies (yellow parts), and an antenna (gray part). Each color corresponds to different surface-coating materials. For detecting the polarization of backscattering light from the remote satellite, it is useful to calculate from the individual parts. If the DOP and polarization states of these regular geometric surfaces are already known, then the discrimination for complicated targets can be simplified, since fewer pulses are required for transmitting between laser and target.
In this paper, based on the specular and diffuse components pBRDF model presented in Ref. [36], a novel method to access the Stokes parameter and DOP of backscattering light from rough surfaces is presented. By giving analytical mathematical formulas of pBRDF for typical geometric surfaces, the Stokes parameters and DOP can thus be calculated as functions of geometric and spatial parameters. Firstly, a practical fitting guidance is given for the diffuse component of BRDF as a function of the incident angle. The fitting coefficients depend only on the surface material, which can be easily measured in application scenarios. Then, the analytical calculations follow to show the dependence on backscattering light DOP, of geometric parameters, such as target size and spatial orientation. Finally, except for the case that the rotation angle of a planar surface is greater than around 50 , the experiment results of DOP of backscattering light from squares and discs agreed with the theoretical calculation within the experiment precision. The work is expected to pave a way for discrimination between the remote target based on the geometric size or spatial orientation, instead of materials.

2. Method

2.1. General Theories of pBRDF

This subsection explains the general theory for the Stokes parameters and microfacets model for pBRDF calculations. The polarization states are described by Stokes vectors,
S = S 0 S 1 S 2 S 3 T
where the four parameters can be defined by the electric field components E x = a 1 e i δ 1 , E y = a 2 e i δ 2 , as:
S 0 = a 1 2 + a 2 2 S 1 = a 1 2 a 2 2 S 2 = 2 a 1 a 2 cos ( δ 1 δ 2 ) S 2 = 2 a 1 a 2 sin ( δ 1 δ 2 )
According to this, the degree of polarization is defined as:
D O P = S 1 2 + S 2 2 + S 3 2 S 0
The input and output Stokes vectors can then be connected by the pBRDF matrix of the rough surfaces, whose elements are F j k , as follows:
F j k ( θ i , θ r , ϕ r ϕ i ) = d S o u t , j S i n , k cos θ i d Ω i
which is expressed by the differential Stokes parameter d S o u t , j and differential solid angle d Ω i , or equivalently, the integral form:
S o u t = Ω i F ( θ i , θ r , ϕ r ϕ i ) · S i n cos θ i d Ω i
where S o u t , j and S i n , k are the jth and kth components of output and input Stokes parameters, respectively, θ i , θ r are incidence and reflectance angles, respectively, ϕ i , ϕ r are azimuth angles of incidence and reflectance directions, respectively, and the integral is over the solid angle Ω i subtended by the target surface at the observation point. Furthermore, the reversibility of light requires the symmetry for F j k such that
F j k ( θ i , θ r , ϕ ) = F j k ( θ r , θ i , ϕ ) = F j k ( θ i , θ r , ϕ )
where ϕ denotes ϕ r ϕ i for simplicity.
To express the explicit form of pBRDF matrix F , it is necessary to describe the scattering behavior of light on the rough surface. Each single reflection on a microfacet obeys Fresnel’s law, thus, for a fixed-incident angle θ i , the reflection angle θ r can vary from 0 to π , since the random distribution of the normal of the microfacet, and vice versa. According to the reversibility of light, it is expected that the non-polarized BRDF is unity if the integrate area Ω r is spread across the space, if there is no energy loss due to absorption of the surface material [35],
0 2 π d ϕ 0 π f ( θ i , θ r , ϕ ) cos θ r sin θ r d θ r = 1
However, for the sake of backscattering light, the reflection angle is confined from 0 to π / 2 only, then the physical meaning of that π / 2 θ r π is that the light experiences multiple reflections on the microfacets. Therefore, the pBRDF matrix can be decomposed as two terms: the single reflection term and multiple reflection term [36], the former is the product of BRDF f s ( θ i , θ r , ϕ ) and Mueller matrix M ,
F = F s + F d = f s ( θ i , θ r , ϕ ) M ( θ i , θ r , ϕ ) + F d
It is of course quite difficult to evaluate the multiple reflection BRDF; however, as a good approximation, for an ideal Lambertian, it can be treated as isotropic diffusion in every direction, which implies f d is a constant for a fixed incident angle. As a result, the diffusion term of BRDF can be expressed as [36]
f d ( θ i ) = 1 π 1 D H R ( θ i )
where the directional hemispherical reflectance (DHR) is defined as
D H R ( θ i ) = 0 2 π d ϕ 0 π / 2 f s ( θ i , θ r , ϕ ) cos θ r sin θ r d θ r
which depends only on the incident angle θ i and the roughness parameters. Then, the corresponding diffuse component of pBRDF is
F 00 d = f d ( θ i ) M 00 ( θ i , θ r , ϕ )
while other matrix elements vanish.
In the long-distance backscattering-detection scenario, where the “long-distance” condition is defined as: the ratio of target size a to light transmission distance L is much smaller than unity, for example, a / L 10 2 , and it is proper to assume that the laser source and detector are coaxial. In this case, the incident and reflection angles are treated exactly the same, θ i = θ r and ϕ = 0 , resulting in the independence of the Mueller matrix element of an incident angle,
M = M 00 0 0 0 0 M 00 0 0 0 0 M 22 M 23 0 0 M 23 M 22
where the elements M 00 , M 22 and M 23 can be determined only by the surface roughness parameters, more explicitly, as the standard deviation of roughness σ h , correlation length l and the complex refractive index, respectively.
To sum up, the pBRDF matrix in case of coaxial system has the following form:
F = ( f s + f d ) M 00 0 0 0 0 f s M 00 0 0 0 0 f s M 22 f s M 23 0 0 f s M 23 f s M 22
If the input light is horizontally linear polarized, then
S i n = I z 1 1 0 0 T
where I z is light intensity at distance L along the propagating direction, by assuming that the laser beam profile is Gaussian, then
I z = I w 0 k L e x 2 + y 2 k 2 L 2
where w 0 is the beam waist and k is the full divergence angle of the laser, and x and y are the coordinates perpendicular to the light propagating direction. Thus, k L is the beam radius at the distance of L.
S o u t , 0 = Ω i F 00 ( θ i , θ r , ϕ ) I z cos θ i d Ω i S o u t , 1 = Ω i F 11 ( θ i , θ r , ϕ ) I z cos θ i d Ω i S o u t , 2 = S o u t , 3 = 0
The degree of polarization can thus be calculated according Equation (3):
D O P = 1 + Ω i f d ( θ i ) I z cos θ i d Ω i Ω i f s ( θ i , θ i , 0 ) I z cos θ i d Ω i 1
which is independent of the Muller matrix elements. Denote these two surface integrals in Equation (17) as I d and I s . Then, the task of this paper is to calculate the degree of polarization of different typical geometric surfaces at a given incident direction to the target objects.

2.1.1. Specular Component

The incident (reflection) angle θ i ( θ r ) is the polar angle between the macro surface’s normal and irradiation (radiation) light. Meanwhile, the angle between the normal of microfacet and macro surface is θ , which is randomly oriented, and is supposed to have a Gaussian distribution probability density function [28,35],
p ( θ ) = 1 2 π σ 2 cos 3 θ e tan 2 θ 2 σ 2
where σ = 2 σ h / l . According to the microfacets model, the specular component of BRDF is [36]
f s ( θ i , θ r , ϕ ) = G ( θ i , θ r , ϕ ) e tan 2 θ 2 σ 2 8 π σ 2 cos 4 θ cos θ i cos θ r
where
G ( θ i , θ r , ϕ ) = m i n 1 , 2 cos θ i cos θ cos β , 2 cos θ r cos θ cos β
is called the shadowing/masking function, which could keep the BRDF finite while the incident or reflection angle is close to π / 2 . Each single reflection on a microfacet obeys Fresnel’s law, thus these geometric relations connecting all angles are defined:
cos 2 β = cos θ i cos θ r + sin θ i sin θ r cos ϕ cos θ = cos θ i + cos θ r 2 cos β
In the practical case of the coaxial system, ϕ i = ϕ r , θ i = θ r = θ , and β = 0 , thus Equation (19) can be simplified for the integration of I s , as follows:
I s = 1 8 π σ 2 Ω G ( θ , θ , 0 ) e tan 2 θ 2 σ 2 I z cos 5 θ d Ω
In our previous work [47], the explicit specular component of pBRDF was calculated analytically or numerically on varied surfaces, such as planar discs, squares, cones and cylinders.

2.1.2. Diffuse Component

In the current work, for calculating the degree of polarization, the remain part is to give the expression for diffuse component I d , and firstly the result of diffuse BRDF f d in Equation (9), based on Equations (10) and (19). It is impossible to obtain an analytic form for f d because of the complexity of the relations of the angles. However, there is an approximate approach to give a relation between f d and cos θ i . Since f d is finite as cos θ i varies in the interval 0 , 1 , one can expand f d as the Taylor series of cos θ i , as follows:
f d ( θ i ) cos θ i = m = 1 a m cos m θ i
where the coefficients a m are determined only by the roughness parameter σ . The numerical simulation of f d ( θ i ) cos θ i is plotted as Figure 2.
Polynomial fitting of these curves shows that it is accurate enough to maintain the first 4 terms of series Equation (23), as shown in Figure 3.
Assume that the equation of the macro surface is F ( x , y , z ) = 0 , and that the light source is at point Q = ( 0 , 0 , 1 ) , then for an arbitrary point P on the surface, the incident angle θ i can be expressed by the coordinates ( x , y , z ) :
cos θ i = n · P Q | P Q |
where n = F x , F y , F z is the unit normal vector of the surface, assuming that it is normalized. The element solid angle d Ω i is
d Ω i = n · P Q | P Q | 3 d S
Substituting Equations (23)–(25) into Equation (17), the approximate expression for I d is:
I d I w 0 k m = 1 4 a m ( n · P Q ) m + 1 e x 2 + y 2 k 2 ( 1 z ) 2 ( 1 z ) | P Q | m + 3 d S
To summarize, if the surface equation, laser parameters and roughness σ are given, one can firstly numerically calculate the diffuse component of BRDF, f d cos θ i , and perform polynomial fitting to obtain the series coefficients. Then, the diffuse component of pBRDF can be calculated through Equation (26).

3. Theoretical Calculations

In this section, the diffuse component of pBRDF is calculated, I d , and thus the degree of polarization for specific typical geometric rough surfaces. The calculation of the specular component can be found in our previous work, Ref. [47].

3.1. General Planar Surfaces

The most representative regular planar surfaces are of course disc and rectangle. Without loss of generality, suppose the surface equation is A x + B y + C z + D = 0 , and the unit normal vector is n = A , B , C . Therefore, according to the geometry,
n · P Q = C + D | P Q | = x 2 + y 2 + ( 1 z ) 2 1 / 2
Substitute Equation (27) into Equation (26), one can obtain
I d = I w 0 k m = 1 4 a m ( C + D ) m + 1 e x 2 + y 2 k 2 ( 1 z ) 2 ( 1 z ) x 2 + y 2 + ( 1 z ) 2 m + 3 2 d S
In the practical applications, the distance between target and the laser source will be much larger than the target size, which implies that x , y , z 1 . Thus, the quadratic or higher terms in the expansion of denominator in Equation (28) can be omitted, which becomes:
( 1 z ) 1 x 2 + y 2 + ( 1 z ) 2 m + 3 2 1 + ( m + 4 ) z = 1 ( m + 4 ) ( D + A x + B y ) C
and
e x 2 + y 2 k 2 ( 1 z ) 2 e x 2 + y 2 k 2

3.2. Disc

As the simplest example of planar surface, the disc is common and widely used in practical remote detection. Suppose that the disc center locates in the origin, and the unit normal vector and radius of the disc are n = 0 , sin α , cos α and R, respectively. The planar equation is y sin α + z cos α = 0 . To integrate Equation (28) analytically, a coordinates transformation was performed,
x y z = 1 0 0 0 cos α sin α 0 sin α cos α x y z
where ( x , y , z ) are the local coordinates on the disc. In this coordinates system, the surface equation is x 2 + y 2 R 2 and z = 0 , or equivalently, the cylindrical coordinates version ρ R and z = 0 . Then, Equation (28) becomes
I d I w 0 k m = 1 4 a m cos m + 2 α e x 2 + y 2 cos 2 α k 2 1 + ( m + 4 ) y sin α d x d y
The integral area is the disc, x 2 + y 2 R 2 . Note that the second term of the integral kernel is odd about y , and thus vanishes through the integration. Thus, Equation (32) becomes
I d = I w 0 k m = 1 4 a m cos m + 2 α 0 R e ρ 2 k 2 ρ d ρ 0 2 π e ρ 2 sin 2 α sin 2 ϕ k 2 d ϕ
Using the well-known Poisson’s integral expression for the Bessel function of order ν ,
J ν ( z ) = ( z / 2 ) ν π Γ ( ν + 1 / 2 ) 0 π e i z cos ϕ sin 2 ν ϕ d ϕ
the second integral of Equation (33) can be explicitly expressed by Bessel function:
0 2 π e ρ 2 sin 2 α sin 2 ϕ k 2 d ϕ = 2 π I 0 ρ 2 sin 2 α 2 k 2 e ρ 2 sin 2 α 2 k 2
where Γ ( x ) is the gamma function, I 0 ( x ) is the modified Bessel function of order 0 of the first kind. Let ξ = ρ 2 / k 2 , finally, the diffuse component of backscattering light is
I d = 2 π I w 0 k m = 1 4 a m cos m + 2 α 0 R I 0 ρ 2 sin 2 α 2 k 2 e ρ 2 k 2 1 sin 2 α 2 ρ d ρ = π I k w 0 f d ( α ) cos 3 α 0 R 2 / k 2 I 0 sin 2 α 2 ξ e 1 sin 2 α 2 ξ d ξ
which in general cannot be simplified further.
From Equation (36), it is easy to see the dependence of I d on all parameters. Firstly, the coefficients a m are dependent on surface roughness σ . Secondly, the laser divergence angle k contributes the upper bound of the integral and overall value. Finally, the geometric parameter, disc radius R, also determines the integral upper bound, while rotation angle α contributes the integral kernel.
Figure 4 shows the curves of diffuse component I d as a function of the disc rotation angle α and the radius R, at four different surface roughness parameters σ , according to Equation (36). As expected, at a fixed geometric configuration ( α , R fixed), as σ increases, the diffuse backscattering light intensity will become stronger, this is in contrast to the specular component, see Ref. [47]. From the physical explanation, σ is the standard variance of the tangent of microfacets normal vector direction, θ . When σ increases, the probability of a large angle of θ is higher, thus the multi-reflection becomes substantial. In Figure 4b, the diffuse component backscattering light intensity tends to be asymptotic as the disc radius increases close to the beam diameter, which coincides with that of the specular component [47].
The degree of polarization of planar discs is plotted in Figure 5. The DOP is high when the rotate angle α is small, and decreases monotonically and rapidly as α increases when the surface is “smooth”, i.e., σ is small, as shown with the blue line in Figure 5a. This is because, in this condition, the specular component of pBRDF dominates. However, when the surface becomes “rougher”, the diffuse component of pBRDF dominates. However, when the surface becomes “rougher”, the diffuse component dominates, and thus DOP becomes low when α is small. The other three lines in Figure 5a imply that DOP first increases before a maximum value then drops rapidly to 0 as α increases. The physical explanation of this phenomenon is that the specular component is more sensitive to the rotation angle than the diffuse component. Therefore, the former contributes to DOP at a relative narrow range of α than the latter. When α is larger than about 60 , I s rapidly vanishes, so that DOP tends to 0.
On the other hand, note that the DOP is calculated in the approximation condition of long distance, R 1 , for example, R 10 2 . The relative error of DOP is less than 10 2 orders of magnitude. Within this accuracy, the DOP is essentially a constant at a certain roughness σ , as Figure 5b suggested. Figure 6 plots the DOP increases with disc radius. The relative change rate of DOP is only about 0.239066 / 0.239050 1 6.7 × 10 5 , as the radius varies from 0 to 0.1, which is 3 orders of magnitudes lower than that of R 2 0.01 . Thus, in this accuracy condition, the contribution of higher order, O ( x 2 ) and O ( y 2 ) , in Equation (32) must be taken into account.

3.3. Rectangle

Another typical planar surface is a rectangle. Suppose the side lengths are a and b, respectively, then the surface equation is the same as the disc discussed in the previous subsection. According to Equation (32), the diffuse component of backscattering light intensity is
I d = I w 0 k m = 1 4 a m cos m + 2 α a 2 a 2 d x b cos α 2 b cos α 2 e x 2 + y 2 cos 2 α k 2 1 + ( m + 4 ) y sin α d y = π I w 0 k f d ( α ) cos 2 α Φ a 2 k Φ b cos α 2 k
where
Φ ( x ) = 2 π 0 x e t 2 d t
is the error function.
It is still worthy considering a spacial case that a = b , a square. Figure 7 shows the diffuse component of backscattering light intensity, I d , as a function of the planar rotation angle α and square length a, respectively. The square length in Figure 7a is the same as the disc diameter in Figure 4a. The behavior of these curves are consistent with that of discs, based on the same physical principles, except the absolute intensity differs, which results from the irradiance areas of the two surfaces being difference.
Additionally, comparing the DOP of square and disc with the same length and diameter at different rotation angle, it is clear that the relative difference is at an order of magnitude of 0.01, as shown in Figure 8. Considering that the calculation accuracy is of 10 2 orders of magnitude, one can conclude that within this precision, the DOP of the disc and square of the same sizes are equal. The reasonable physical meaning is that the DOP of a planar surface is determined essentially by the material (the roughness σ ) and incident angle α , rather than the geometric size or shape.

3.4. Typical Curved Surfaces

For the typical curved surfaces, the most common ones are cone, hemisphere and cylinder. For calculation of the former one, first consider a general situation that the arbitrary revolved surface described by equation F ( ρ , z ) = 0 . Assuming that ρ F 2 + z F 2 = 1 , according to Equation (26), the diffuse component backscattering light intensity is
I d = I w 0 k m = 1 4 a m z F ρ ρ F + z z F m + 1 ( 1 z ) ρ 2 + ( 1 z ) 2 ( m + 3 ) / 2 e ρ 2 k 2 ( 1 z ) 2 d S
where the element area is
d S = 1 | ρ F | ρ d ϕ d z
For a cone with a half vertex angle α ( 0 , π ) and bottom radius R, or height H, the surface equation is ρ cos α + z sin α = 0 . Under this definition, the vertex of the cone is upward when α ( 0 , π / 2 ) and downward otherwise. Then, according to geometry relation of this configuration, the element area d S = | sec α | ρ d ϕ d z = csc α ρ d ρ d ϕ . Thus,
I d = 2 π I w 0 k m = 1 4 a m sin m α 0 R e ρ 2 k 2 ( 1 + ρ cot α ) 2 ( 1 + ρ cot α ) ( 1 + ρ 2 csc 2 α + 2 ρ cot α ) m + 3 2 ρ d ρ
Very carefully expand the denominator of Equation (41) as series. It can be carried out only under the long-distance condition that both R 1 and H 1 . When R H , or equivalently, | cot α | 1 ,
I d 2 π I w 0 k m = 1 4 a m sin m α 0 R e ρ 2 k 2 1 ( m + 4 ) ρ cot α ρ d ρ = I 1 d + I 2 d
where
I 1 d = π I w 0 k m = 1 4 a m sin m α 1 e R 2 / k 2
I 2 d = π I w 0 k cot α m = 1 4 ( m + 4 ) a m sin m α R e R 2 / k 2 k π 2 Φ R k
are the first and second orders of approximation for I d , respectively. As Figure 9 shows, the former is the diffuse component of backscattering light from the projection on O x y plane of the cone. If α = π / 2 , the cone is reduced to a disc of radius R, which can be easily verified by comparing the results of Equation (43) with α = π / 2 and Equation (36) with α = 0 . In the second order, Equation (44) is odd around α = π / 2 , which is the consequence to the non-planar property.
Similarly, when R H , or equivalently, | tan α | 1 , Equation (42) can be expanded as
I d 2 π I w 0 k m = 1 4 a m sin m α 0 H e z 2 tan 2 α k 2 tan 2 α 1 + ( m + 4 ) z z d z = I 1 d + I 2 d
where now
I 1 d = π I w 0 k m = 1 4 a m sin m α 1 e H 2 tan 2 α / k 2
I 2 d = π I w 0 k cos α m = 1 4 ( m + 4 ) a m sin m 1 α k π 2 Φ H tan α k H tan α e H 2 tan 2 α k 2
Figure 10 shows the degree of polarization of backscattering light from the cone at different roughness varies as functions of the half-vertex angle and bottom radius, respectively. The DOP versus α exhibits a similar behavior as that of a disc, Figure 5a, since the major contribution portion is I 1 d . However, the DOP is not symmetrical around α = π / 2 due to the presence of I 2 d . It must be emphasized that in the realistic target detection, the difference between D O P ( π / 2 + α ) and D O P ( π / 2 α ) is expected to be larger than the result given in this article. The reason is when the cone half-vertex angle π / 2 < α < π , the multi-reflection occurs not only between the microfacets, but also the macro-internal surface of the cone. This effect will be studied in the future work. Figure 10b shows the dependence of DOP on the bottom radius of the cone. For a specific roughness σ , it exhibits a monotonic increase as R increases, until R reaches the beam radius k ( 1 z ) k , resulting from the increasing irradiated area by the light source.

4. Results

Figure 11 shows the experimental setup for testing the polarization state of backscattering light from typical geometric targets. A 532 nm pulsed laser of 10 μ J output energy and nanosecond pulse width is modulated by a signal generator at 10 Hz repetition frequency. An optical attenuator consists by a half-wave plate and two polarized beam splitters (PBSs), after which the light is horizontally polarized. The extinction ratio of each PBS is about 30 dB, then the total extinction ratio is 60 dB. The convex lens diverges the beam to about a 30 cm radius, which completely covers the target surface, located 15 m away from the laser. The beam divergence k, appearing in Equation (15), thus is k = 30 cm/15 m = 0.02. The laser angle of pitch was adjusted to be vertical to the reference board. The pedestal can be rotated to control the incident angle. The backscattering beam is collected by a optical camera with a photon detector integrated inside. Finally the wave pulses are recorded by an oscilloscope.
The Stokes parameter S 0 , S 1 and S 2 can be measured by rotating the polarizer before the camera. Since the incidence light is horizontally polarized and the target surface material is non-conductive, the refractive index is real. According to Equation (13), the nonvanishing components of the backscattering Stokes parameter can only be S 0 and S 1 . When rotating the polarizer, the detector response has maximum and minimum values, I m a x and I m i n , which can be expressed by Stokes parameters as
I m a x = S 0 + 1 2 S 0 S 1 I m i n = 1 2 S 0 S 1
Thus, the degree of polarization is
D O P = I m a x I m i n I m a x + I m i n
The targets models are made by hardened white card paper, with a roughness parameter σ estimated to be at 0.1∼1 order of magnitude. The geometric parameters are listed in Table 1. The first two columns are planar targets, for verifying DOP as either the size or rotating angle, corresponding to the theoretical calculation result of Figure 5. The vertex angles of cones are fixed at α = 63.4 to verify the influence of the bottom radius on the degree of polarization.

5. Discussion

To verify the dependency of DOP on the incident angle, the DOP of squares and discs were tested. Since the energy of pulses fluctuated, the measured pulse signal introduced errors, which could reduce SNR and affect the integration to obtain total pulse energy. To suppress the measurement error caused by laser power fluctuation, the experiment repeated 10 times for each geometric surface, the square and disc. Within each repetition, the pedestal angle α varies from 0 to 60 in every 5 , while the polarizer angle ψ independently varies from 0 to 90 in every 5 . In every angle set ( α , ψ ) , 20 averaged pulses were measured to suppress the detector electronic noise. Figure 12 shows the measurement results of DOP of these two kinds of planar surfaces, which exhibits basically the same behavior that increases monotonically as the rotation angle increases. Comparing the theoretical plot and measured data in Figure 12, the tendency of DOP curves coincides with the theoretical calculation. Specifically, when α is less than about 20 , the mean measured results coincide with the theoretical calculation, while when α is larger than about 25 , the mean measured results are a bit less than the theoretical results. Additionally, with the rotation angle α becoming larger, the measured data (the dots in Figure 12) spread and variate from the mean value. There are two major reasons. One is that when the rotation angle is large, the total intensity of the backscattering light is very weak, while the environment light from the reference board (see Figure 11) is inevitable. While α increases, the intensity from target falls, then the component of the environment light dominates, which depolarizes the backscattering light substantially. This can be qualitatively analyzed as follows. If the unpolarized environment light intensity is I e n , then the measurement results are
I m a x = S 0 + 1 2 S 0 S 1 + 1 2 I e n I m i n = 1 2 S 0 S 1 + 1 2 I e n
respectively. Thus, the realistic DOP is
D O P m e a s u r e = I m a x I m i n I m a x + I m i n + I e n < D O P t h e o r e
On the other hand, the total light intensity rapidly decreases with α increasing, and the signal-to-noise ratio falls, making it difficult to extract the signal pulse from the backscattering light waveform. Thus, the measurement error is relatively large compared with small α , which explains the spread of the dots of the measured data in Figure 12.
To verify the dependence of DOP on geometric size, five squares and five discs with different lengths were first tested. The pedestal angle was fixed such that the light was of normal incidence and measured backscattered light pulses in every 5 of the polarizer angle. Each measurement result also contains 20 averaged laser pulses backscattered from the target. Furthermore, the experiment also repeated 10 times for each target. The experiment results are shown in Figure 13a,b. The experiment results agreed well with the theoretical calculations. However, when the size of squares or discs is small, the measured data show more fluctuations than that of the large-sized squares or discs. The reason is that, when the target size is small, the total light intensity that comes from the target surface is weak, leading to a low SNR. Thus, the extracted waveform from the signal pulse has a relatively larger deviation, as the dots suggested at a = 8 cm or r = 4 cm in Figure 13a and Figure 13b, respectively.
On the other hand, the degree of polarization of five cones with different bottom radii and fixed vertex angles were measured, as shown in Figure 13c. The DOP behavior is similar to that of the former two geometric shapes and basically coincides with the theoretical calculation. However, for the cones, there is another source of error. The cones are cylindrically symmetrical when the light is of normal incidence. The backscattered light from an arbitrary generatrix at the azimuth angle of ϕ is equal to that from a plane whose normal direction is a function of ϕ and α . Thus, the DOP contribution is expected to be axial symmetric around the polarization direction of incident light. However, the beam profile deviates from Gaussian, destroying such axial symmetry. Therefore, the DOP of backscattered light depends on the beam profile, resulting in extra measurement errors for different sizes of cones.
The goal of the theoretical calculation of backscattering Stokes parameters and DOP is to discriminate remote targets by using the LIDAR. There are two major ways to give the pBRDF and thus the Stokes parameters in the current literature. One way is the measuring method [17,48,49,50,51]. For a specific target, for example, a plane with a coating material whose surface parameters like roughness or refractive index are unknown, measure the polarization states of backscattering light at different incidence and reflectance angles ( θ i , θ r , ϕ ) . Therefore, according to Equation (4), the pBRDF as a function of these angles is given. However, the measurement method suffers from a disadvantage of un-universality. If the structures or surface parameters change, the measured pBRDF is invalid.
The other way to obtain the polarization state from the target is numerical simulation calculation by dividing the rough surface into smaller but finite facets [14,37,52,53]. There are two simulation methods, one is the Method of Moments (MoM) [36], which is to numerically solve integral equations about electromagnetic fields. It is an accurate solution to the electric field vector, however, it is confined to give the polarization state of backscattering light from a small area surface, as Ref. [36], millimeter level, inadequate for the purpose of realistic target detection. The other one is directly summing the backscattering light from all small planar facets. However, the parameters like normal directions or roughness are different in each facets. However, the numerical simulation method costs much computational time and resources due to the complexity of solving integral equations or calculating the 16 Mueller matrix elements for each small facet, as the number of facets increase, especially when the realistic targets are complex surfaces, while the simulation accuracy increases in contrast.
Fortunately, typical geometric surfaces have analytic pBRDF formulas. The previous work [47], as well as this current article, derived the mathematical formulas of pBRDF for some typical geometric surfaces that provide another possible way to access the pBRDF and Stokes parameter of the backscattering light with less computational time required. What is more essential, one can discriminate the target geometric size and spatial orientation. More explicitly, according to both theoretical and experimental results in this article, the DOP is dependent substantially on the rotation angle of the target, but insensitive to the size of the target. On the other hand, according to the results of our previous work [47], the Stokes parameter S 0 and S 1 stand in contrast. Based on this fact, in the further work of remote sensing applications, assuming that the targets are rigid bodies, and the geometric parameters, denoted by ( q 1 , q 2 , ) , such as radius, heights, lengths, and vertex angles, are fixed for a given target, but the gesture parameters, denoted by ( p 1 , p 2 , ) , such as spatial positions x, y, the rotation angles, can vary, then the Stokes parameters of backscattering light are the function of these parameters:
S = S ( q 1 , q 2 , , p 1 , p 2 , )
where the explicit form of the function is analytically calculated by the formulas shown in both current and previous works [47] for typical geometric surfaces. The DOP can be calculated by Equation (3) if the four Stokes parameters are measured. As the target moves, the gesture parameters vary, thus the measured Stokes parameters consequently change,
S = S ( q 1 , q 2 , , p 1 , p 2 , ) S = S ( q 1 , q 2 , , p 1 , p 2 , ) S = S ( q 1 , q 2 , , p 1 , p 2 , )
Then, according to the specific form of the function, one can recover the gesture parameters, expressed by the measured Stokes and geometric parameters, leading to the discrimination of target sizes and spatial status.
However, the mathematical formula method has its limitations as well. Firstly, the Stokes parameters of realistic targets with complicated geometric structures can usually be calculated by summing up the Stokes parameters from all typical geometric surface parts, as shown in Figure 1. The multi-reflection between different parts could further depolarize the backscattering light while the shadowing effect could also change the polarization, which is worth considering in the future work. Secondly, the accuracy of the pBRDF model, especially the shadowing and masking function, G ( θ i , θ r , ϕ ) model, plays a critical role in the theoretical calculation of the polarization states of the backscattering light, including DOP. Therefore, the development of a higher-precision pBRDF model is urgently needed to aid in further calculations, assessing more complicated targets, and for applications in remote sensing.

6. Conclusions

To summarize, this work provides an alternative method to access the pBRDF of a rough surface, and thus the Stokes parameters of backscattering light, by giving analytical formulas of typical geometric surfaces. Compared to the current ways, the direct measurement or numerical simulation, the method presented here has the advantage of universality and simplicity. For theoretical calculations, the model of diffuse components of BRDF was firstly simplified, giving a polynomial fitting, which is practical for the theoretical analysis of backscattering Stokes parameters. Then, by using this fitting model of diffuse component BRDF f d ( cos θ i ) , the analytical calculation of DOP of backscattering light from several kinds of typical geometric rough surfaces is established. The theoretical calculations for different sizes of planar and cone targets agreed well with experiment results within our measurement accuracy, the latter represents a current method to obtain pBRDF. The results for planar targets at different rotation angles basically agreed with the theoretical calculations when the rotation angle is less than about 50 , but deviate when it increased. The imperfection of results is caused by the non-Gaussian distribution of the beam profile; the accuracy of the optical element, such as the polarizer; and the unpolarized environmental reflection light from background objects, such as the reference board, the pedestal or the connecting stick, which decreases the DOP. The agreement of the mathematical formulas and measurement results provides the reliability for the potential application in remote targets discrimination. However, the mathematical formula method of backscattering Stokes parameters still needs to be improved in the future practical application scenarios, considering the combination of typical geometric surfaces to the realistic complex target.

Author Contributions

Conceptualization, S.S.; Data curation, X.Z. and Y.L.; Supervision, S.X. and J.F.; Project administration, Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 61871389), the Advanced Laser Technology Laboratory of Anhui Province, China (No. KY21C608), the Research Plan Project of the National University of Defense Technology (No. ZK18-01-02) and the Post-graduate Scientic Research Innovation Project of Hunan Province (No. CX20220007).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ratliff, B.M.; LeMaster, D.A.; Mack, R.T.; Villeneuve, P.V.; Weinheimer, J.J.; Middendorf, J.R. Detection and tracking of RC model aircraft in LWIR microgrid polarimeter data. In Polarization Science and Remote Sensing V; International Society for Optics and Photonics: Bellingham, WA, USA, 2011; Volume 8160, p. 816002. [Google Scholar]
  2. Gavrilov, N.; Kshevetskii, S. Verifications of the nonlinear numerical model and polarization relations of atmospheric acoustic-gravity waves. Geosci. Model Dev. Discuss. 2014, 7, 7805–7822. [Google Scholar]
  3. Lavigne, D.A.; Breton, M.; Fournier, G.; Charette, J.F.; Pichette, M.; Rivet, V.; Bernier, A.P. Target discrimination of man-made objects using passive polarimetric signatures acquired in the visible and infrared spectral bands. In Polarization Science and Remote Sensing V; International Society for Optics and Photonics: Bellingham, WA, USA, 2011; Volume 8160, p. 816007. [Google Scholar]
  4. Forssell, G.; Hedborg-Karlsson, E. Measurements of polarization properties of camouflaged objects and of the denial of surfaces covered with cenospheres. In Targets and Backgrounds IX: Characterization and Representation; International Society for Optics and Photonics: Bellingham, WA, USA, 2003; Volume 5075, pp. 246–258. [Google Scholar]
  5. Li, X.; Zhang, L.; Qi, P.; Zhu, Z.; Xu, J.; Liu, T.; Zhai, J.; Hu, H. Are Indices of Polarimetric Purity Excellent Metrics for Object Identification in Scattering Media? Remote Sens. 2022, 14, 4148. [Google Scholar] [CrossRef]
  6. Duggin, M.J.; Kinn, G.J. Vegetative target enhancement in natural scenes using multiband polarization methods. In Polarization Analysis and Measurement IV; SPIE: Bellingham, WA, USA, 2002; Volume 4481, pp. 281–291. [Google Scholar]
  7. Li, S.; Han, X.; Weng, F. Monitoring Land Vegetation from Geostationary Satellite Advanced Himawari Imager (AHI). Remote Sens. 2022, 14, 3817. [Google Scholar] [CrossRef]
  8. Zhang, Z.; Yan, L.; Jiang, X.; Ding, J.; Zhang, F.; Jiang, K.; Shang, K. Exploring the Potential of Optical Polarization Remote Sensing for Oil Spill Detection: A Case Study of Deepwater Horizon. Remote Sens. 2022, 14, 2398. [Google Scholar] [CrossRef]
  9. Qu, Y.; Liang, S.; Liu, Q.; He, T.; Liu, S.; Li, X. Mapping surface broadband albedo from satellite observations: A review of literatures on algorithms and products. Remote Sens. 2015, 7, 990–1020. [Google Scholar] [CrossRef] [Green Version]
  10. Ding, A.; Jiao, Z.; Dong, Y.; Zhang, X.; Peltoniemi, J.I.; Mei, L.; Guo, J.; Yin, S.; Cui, L.; Chang, Y.; et al. Evaluation of the snow albedo retrieved from the snow kernel improved the Ross-Roujean BRDF model. Remote Sens. 2019, 11, 1611. [Google Scholar] [CrossRef] [Green Version]
  11. Shaw, J.A. Degree of linear polarization in spectral radiances from water-viewing infrared radiometers. Appl. Opt. 1999, 38, 3157–3165. [Google Scholar] [CrossRef]
  12. Hieronymi, M. Polarized reflectance and transmittance distribution functions of the ocean surface. Opt. Express 2016, 24, A1045–A1068. [Google Scholar] [CrossRef]
  13. Touzi, R.; Hurley, J.; Vachon, P.W. Optimization of the degree of polarization for enhanced ship detection using polarimetric RADARSAT-2. IEEE Trans. Geosci. Remote Sens. 2015, 53, 5403–5424. [Google Scholar] [CrossRef]
  14. Zhang, Y.; Xuan, J.; Zhao, H.; Song, P.; Zhang, Y.; Xu, W. Improved atmospheric effect elimination method for the roughness estimation of painted surfaces. Opt. Lett. 2018, 43, 1079–1082. [Google Scholar] [CrossRef]
  15. Li, J.; Qiu, S.; Zhang, Y.; Yang, B.; Gao, C.; Qian, Y.; Liu, Y.; Zhao, Y. Assessment of BRDF Impact on VIIRS DNB from Observed Top-of-Atmosphere Reflectance over Dome C in Nighttime. Remote Sens. 2021, 13, 301. [Google Scholar] [CrossRef]
  16. Guo, R.; Jiang, Z.; Jin, Z.; Zhang, Z.; Zhang, X.; Guo, L.; Hu, Y. Reflective Tomography Lidar Image Reconstruction for Long Distance Non-Cooperative Target. Remote Sens. 2022, 14, 3310. [Google Scholar] [CrossRef]
  17. Brown, A.J.; Michaels, T.I.; Byrne, S.; Sun, W.; Titus, T.N.; Colaprete, A.; Wolff, M.J.; Videen, G.; Grund, C.J. The case for a modern multiwavelength, polarization-sensitive LIDAR in orbit around Mars. J. Quant. Spectrosc. Radiat. Transf. 2015, 153, 131–143. [Google Scholar] [CrossRef] [Green Version]
  18. Brown, A.J. Equivalence relations and symmetries for laboratory, LIDAR, and planetary Müeller matrix scattering geometries. JOSA A 2014, 31, 2789–2794. [Google Scholar] [CrossRef] [PubMed]
  19. Dong, Q.; Huang, Z.; Li, W.; Li, Z.; Song, X.; Liu, W.; Wang, T.; Bi, J.; Shi, J. Polarization Lidar Measurements of Dust Optical Properties at the Junction of the Taklimakan Desert–Tibetan Plateau. Remote Sens. 2022, 14, 558. [Google Scholar] [CrossRef]
  20. Kong, Z.; Yin, Z.; Cheng, Y.; Li, Y.; Zhang, Z.; Mei, L. Modeling and evaluation of the systematic errors for the polarization-sensitive imaging lidar technique. Remote Sens. 2020, 12, 3309. [Google Scholar] [CrossRef]
  21. Brown, A.J. Spectral curve fitting for automatic hyperspectral data analysis. IEEE Trans. Geosci. Remote Sens. 2006, 44, 1601–1608. [Google Scholar] [CrossRef] [Green Version]
  22. Wen, N.; Zeng, F.; Dai, K.; Li, T.; Zhang, X.; Pirasteh, S.; Liu, C.; Xu, Q. Evaluating and Analyzing the Potential of the Gaofen-3 SAR Satellite for Landslide Monitoring. Remote Sens. 2022, 14, 4425. [Google Scholar] [CrossRef]
  23. Yang, M.; Xu, W.; Sun, Z.; Jia, A.; Xiu, P.; Chen, W.; Li, L.; Zheng, C.; Li, J. Degree of polarization modeling based on modified microfacet pBRDF model for material surface. Opt. Commun. 2019, 453, 124390. [Google Scholar] [CrossRef]
  24. Shen, Y.; Chen, B.; He, C.; He, H.; Guo, J.; Wu, J.; Elson, D.S.; Ma, H. Polarization Aberrations in High-Numerical-Aperture Lens Systems and Their Effects on Vectorial-Information Sensing. Remote Sens. 2022, 14, 1932. [Google Scholar] [CrossRef]
  25. Liang, J.; Ren, L.; Ju, H.; Zhang, W.; Qu, E. Polarimetric dehazing method for dense haze removal based on distribution analysis of angle of polarization. Opt. Express 2015, 23, 26146–26157. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, X.; Hu, T.; Li, D.; Guo, K.; Gao, J.; Guo, Z. Performances of polarization-retrieve imaging in stratified dispersion media. Remote Sens. 2020, 12, 2895. [Google Scholar] [CrossRef]
  27. Nicodemus, F.E. Directional reflectance and emissivity of an opaque surface. Appl. Opt. 1965, 4, 767–775. [Google Scholar] [CrossRef]
  28. Torrance, K.E.; Sparrow, E.M. Theory for off-specular reflection from roughened surfaces. Josa 1967, 57, 1105–1114. [Google Scholar] [CrossRef]
  29. Nicodemus, F.E.; Richmond, J.C.; Hsia, J.J.; Ginsberg, I.; Limperis, T. Geometrical considerations and nomenclature for reflectance. NBS Monogr. 1992, 160, 4. [Google Scholar]
  30. Leader, J. Bidirectional scattering of electromagnetic waves from rough surfaces. J. Appl. Phys. 1971, 42, 4808–4816. [Google Scholar] [CrossRef]
  31. Beckmann, P.; Spizzichino, A. The Scattering of Electromagnetic Waves from Rough Surfaces; Artech House, Inc.: Norwood, MA, USA, 1987. [Google Scholar]
  32. Valenzuela, G. Depolarization of EM waves by slightly rough surfaces. IEEE Trans. Antennas Propag. 1967, 15, 552–557. [Google Scholar] [CrossRef]
  33. Barrick, D. Rough surface scattering based on the specular point theory. IEEE Trans. Antennas Propag. 1968, 16, 449–454. [Google Scholar] [CrossRef]
  34. Flynn, D.S.; Alexander, C. Polarized surface scattering expressed in terms of a bidirectional reflectance distribution function matrix. Opt. Eng. 1995, 34, 1646–1650. [Google Scholar]
  35. Priest, R.G.; Meier, S.R. Polarimetric microfacet scattering theory with applications to absorptive and reflective surfaces. Opt. Eng. 2002, 41, 988–993. [Google Scholar] [CrossRef]
  36. Hyde, M.W., IV; Schmidt, J.D.; Havrilla, M.J. A geometrical optics polarimetric bidirectional reflectance distribution function for dielectric and metallic surfaces. Opt. Express 2009, 17, 22138–22153. [Google Scholar] [CrossRef] [PubMed]
  37. Zhan, H.; Voelz, D.G. Modified polarimetric bidirectional reflectance distribution function with diffuse scattering: Surface parameter estimation. Opt. Eng. 2016, 55, 123103. [Google Scholar] [CrossRef]
  38. Renhorn, I.G.; Hallberg, T.; Bergström, D.; Boreman, G.D. Four-parameter model for polarization-resolved rough-surface BRDF. Opt. Express 2011, 19, 1027–1036. [Google Scholar] [CrossRef] [PubMed]
  39. Renhorn, I.G.; Hallberg, T.; Boreman, G.D. Efficient polarimetric BRDF model. Opt. Express 2015, 23, 31253–31273. [Google Scholar] [CrossRef] [PubMed]
  40. Liu, H.; Zhu, J.; Wang, K.; Xu, R. Polarized BRDF for coatings based on three-component assumption. Opt. Commun. 2017, 384, 118–124. [Google Scholar] [CrossRef]
  41. Rowe, M.; Pugh, E.; Tyo, J.S.; Engheta, N. Polarization-difference imaging: A biologically inspired technique for observation through scattering media. Opt. Lett. 1995, 20, 608–610. [Google Scholar] [CrossRef]
  42. Lavigne, D.A.; Breton, M.; Pichette, M.; Larochelle, V.; Simard, J.R. Evaluation of active and passive polarimetric electro-optic imagery for civilian and military targets discrimination. In Polarization: Measurement, Analysis, and Remote Sensing VIII; SPIE: Bellingham, WA, USA, 2008; Volume 6972, pp. 285–293. [Google Scholar]
  43. Vermeulen, A.; Devaux, C.; Herman, M. Retrieval of the scattering and microphysical properties of aerosols from ground-based optical measurements including polarization. I. Method. Appl. Opt. 2000, 39, 6207–6220. [Google Scholar] [CrossRef]
  44. Sun, Y. Statistical ray method for deriving reflection models of rough surfaces. JOSA A 2007, 24, 724–744. [Google Scholar] [CrossRef]
  45. Liu, H.; Zhu, J.; Wang, K. Modified polarized geometrical attenuation model for bidirectional reflection distribution function based on random surface microfacet theory. Opt. Express 2015, 23, 22788–22799. [Google Scholar] [CrossRef]
  46. Sun, L.; Zhao, F. Geometric attenuation factor based on scattering theory from randomly rough surface. Appl. Opt. 2021, 60, 476–483. [Google Scholar] [CrossRef]
  47. Shen, S.; Zhang, X.; Liu, Y.; Fang, J.; Xu, S.; Hu, Y. Calculation of Stokes vector of laser backscattering from typical geometric rough surfaces at a long distance. Appl. Opt. 2022, 61, 1766–1777. [Google Scholar] [CrossRef] [PubMed]
  48. Kalantari, E.; Molan, Y.E. Analytical BRDF model for rough surfaces. Optik 2016, 127, 1049–1055. [Google Scholar] [CrossRef]
  49. Zhang, Y.; Zhang, Y.; Zhao, H.; Wang, Z. Improved atmospheric effects elimination method for pBRDF models of painted surfaces. Opt. Express 2017, 25, 16458–16475. [Google Scholar] [CrossRef] [PubMed]
  50. Jiang, Y.; Li, Z. Mueller matrix of laser scattering by a two-dimensional randomly rough surface. J. Quant. Spectrosc. Radiat. Transf. 2022, 287, 108225. [Google Scholar] [CrossRef]
  51. Prokopenko, V.; Alekseev, S.; Matveev, N.; Popov, I. Simulation of the polarimetric bidirectional reflectance distribution function. Opt. Spectrosc. 2013, 114, 961–964. [Google Scholar] [CrossRef]
  52. Wang, S.; Xue, L.; Yan, K. Numerical calculation of light scattering from metal and dielectric randomly rough Gaussian surfaces using microfacet slope probability density function based method. J. Quant. Spectrosc. Radiat. Transf. 2017, 196, 183–200. [Google Scholar] [CrossRef]
  53. Letnes, P.A.; Maradudin, A.A.; Nordam, T.; Simonsen, I. Calculation of the Mueller matrix for scattering of light from two-dimensional rough surfaces. Phys. Rev. A 2012, 86, 031803. [Google Scholar] [CrossRef]
Figure 1. Model of a LIDAR imaging satellite, which consists of several typical geometric parts.
Figure 1. Model of a LIDAR imaging satellite, which consists of several typical geometric parts.
Remotesensing 14 06001 g001
Figure 2. Numerical simulation for f d ( θ i ) cos θ i as the incident angle θ i varies from 0 to 90 at different typical roughness parameter σ from 0.3 to 1.0.
Figure 2. Numerical simulation for f d ( θ i ) cos θ i as the incident angle θ i varies from 0 to 90 at different typical roughness parameter σ from 0.3 to 1.0.
Remotesensing 14 06001 g002
Figure 3. Polynomial fitting of the f d cos θ i as cos θ i at σ = 0.8 , keeping only the lowest 4 orders in Equation (23). The coefficients are: a 1 = 0.0289 , a 2 = 0.125 , a 3 = 0.607 and a 4 = 0.313 .
Figure 3. Polynomial fitting of the f d cos θ i as cos θ i at σ = 0.8 , keeping only the lowest 4 orders in Equation (23). The coefficients are: a 1 = 0.0289 , a 2 = 0.125 , a 3 = 0.607 and a 4 = 0.313 .
Remotesensing 14 06001 g003
Figure 4. Diffuse component of backscattering light at different roughness σ , varies as (a): rotation angle α of the disc, at fixed radius R = 0.0053 ; and (b): disc radius R at fixed rotation angle α = 0 . The laser parameter k = 0.02 .
Figure 4. Diffuse component of backscattering light at different roughness σ , varies as (a): rotation angle α of the disc, at fixed radius R = 0.0053 ; and (b): disc radius R at fixed rotation angle α = 0 . The laser parameter k = 0.02 .
Remotesensing 14 06001 g004
Figure 5. Degree of polarization of backscattering light at different roughness σ , varies as (a) rotation angle α of the disc at fixed radius R = 0.0053 ; and (b) disc radius R at fixed rotation angle α = 0 . The laser parameter k = 0.02 . The non-smoothness of the curves at about α π / 4 results from the shadowing/masking function Equation (20).
Figure 5. Degree of polarization of backscattering light at different roughness σ , varies as (a) rotation angle α of the disc at fixed radius R = 0.0053 ; and (b) disc radius R at fixed rotation angle α = 0 . The laser parameter k = 0.02 . The non-smoothness of the curves at about α π / 4 results from the shadowing/masking function Equation (20).
Remotesensing 14 06001 g005
Figure 6. DOP versus disc radius in a higher precision, at α = 0 , σ = 0.8 and k = 0.02 .
Figure 6. DOP versus disc radius in a higher precision, at α = 0 , σ = 0.8 and k = 0.02 .
Remotesensing 14 06001 g006
Figure 7. Diffuse component of backscattering light intensity at different roughness σ , varies as (a) rotation angle α of the square at fixed a = 0.01 ; and (b) square length a at fixed rotation angle α = 0 . The laser parameter k = 0.02 .
Figure 7. Diffuse component of backscattering light intensity at different roughness σ , varies as (a) rotation angle α of the square at fixed a = 0.01 ; and (b) square length a at fixed rotation angle α = 0 . The laser parameter k = 0.02 .
Remotesensing 14 06001 g007
Figure 8. The DOP relative difference between a square of length a and a disc of diameter a = 2 R , at α = 0 , σ = 0.8 and k = 0.02 .
Figure 8. The DOP relative difference between a square of length a and a disc of diameter a = 2 R , at α = 0 , σ = 0.8 and k = 0.02 .
Remotesensing 14 06001 g008
Figure 9. First two orders of approximation for the diffuse component of backscattering light intensity at different roughness σ = 0.8 , varies as (a) half vertex angle α of the cone at fixed radius R = 0.0067 ; and (b) bottom radius R at fixed half vertex angle α = π / 4 . The laser parameter k = 0.02 .
Figure 9. First two orders of approximation for the diffuse component of backscattering light intensity at different roughness σ = 0.8 , varies as (a) half vertex angle α of the cone at fixed radius R = 0.0067 ; and (b) bottom radius R at fixed half vertex angle α = π / 4 . The laser parameter k = 0.02 .
Remotesensing 14 06001 g009
Figure 10. DOP of backscattering light from the cone at different roughness varies as (a) half vertex angle α of the cone at fixed radius R = 0.0067 ; and (b) bottom radius R at fixed half vertex angle α = π / 4 . The laser parameter k = 0.02 .
Figure 10. DOP of backscattering light from the cone at different roughness varies as (a) half vertex angle α of the cone at fixed radius R = 0.0067 ; and (b) bottom radius R at fixed half vertex angle α = π / 4 . The laser parameter k = 0.02 .
Remotesensing 14 06001 g010
Figure 11. Experimental setup for testing the polarization state of backscattering light from typical geometric targets with rough surfaces.
Figure 11. Experimental setup for testing the polarization state of backscattering light from typical geometric targets with rough surfaces.
Remotesensing 14 06001 g011
Figure 12. Measured DOP of backscattering light from planar squares and discs rough surfaces.
Figure 12. Measured DOP of backscattering light from planar squares and discs rough surfaces.
Remotesensing 14 06001 g012
Figure 13. Measured DOP of backscattering light from different typical geometric rough surfaces. (a) Squares of length: 8 cm, 10 cm, 12 cm, 14 cm and 16 cm. (b) Discs of radius: 4 cm, 5 cm, 6 cm, 7 cm and 8 cm. (c) Cones of bottom radius: 4 cm, 5 cm, 6 cm, 7 cm and 8 cm.
Figure 13. Measured DOP of backscattering light from different typical geometric rough surfaces. (a) Squares of length: 8 cm, 10 cm, 12 cm, 14 cm and 16 cm. (b) Discs of radius: 4 cm, 5 cm, 6 cm, 7 cm and 8 cm. (c) Cones of bottom radius: 4 cm, 5 cm, 6 cm, 7 cm and 8 cm.
Remotesensing 14 06001 g013
Table 1. Geometric parameters of target models.
Table 1. Geometric parameters of target models.
SquaresDiscsCones
LengthRadiusRadiusHeight α
8 cm4 cm4 cm2 cm63
10 cm5 cm5 cm2.5 cm63
12 cm6 cm6 cm3 cm63
14 cm7 cm7 cm3.5 cm63
16 cm8 cm8 cm4 cm63
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shen, S.; Zhang, X.; Liu, Y.; Xu, S.; Fang, J.; Hu, Y. Degree of Polarization Calculation for Laser Backscattering from Typical Geometric Rough Surfaces at Long Distance. Remote Sens. 2022, 14, 6001. https://doi.org/10.3390/rs14236001

AMA Style

Shen S, Zhang X, Liu Y, Xu S, Fang J, Hu Y. Degree of Polarization Calculation for Laser Backscattering from Typical Geometric Rough Surfaces at Long Distance. Remote Sensing. 2022; 14(23):6001. https://doi.org/10.3390/rs14236001

Chicago/Turabian Style

Shen, Shiyang, Xinyuan Zhang, Yifan Liu, Shilong Xu, Jiajie Fang, and Yihua Hu. 2022. "Degree of Polarization Calculation for Laser Backscattering from Typical Geometric Rough Surfaces at Long Distance" Remote Sensing 14, no. 23: 6001. https://doi.org/10.3390/rs14236001

APA Style

Shen, S., Zhang, X., Liu, Y., Xu, S., Fang, J., & Hu, Y. (2022). Degree of Polarization Calculation for Laser Backscattering from Typical Geometric Rough Surfaces at Long Distance. Remote Sensing, 14(23), 6001. https://doi.org/10.3390/rs14236001

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop