Imaging Components, Systems, and Processing

Relating the statistics of the angle of linear polarization to measurement uncertainty of the Stokes vector

[+] Author Affiliations
Meredith Kupinski

University of Arizona, College of Optical Sciences, Tucson, Arizona 85721, United States

University of Arizona, Department of Medical Imaging, Tucson, Arizona 85721, United States

Russell Chipman

University of Arizona, College of Optical Sciences, Tucson, Arizona 85721, United States

Eric Clarkson

University of Arizona, College of Optical Sciences, Tucson, Arizona 85721, United States

University of Arizona, Department of Medical Imaging, Tucson, Arizona 85721, United States

Opt. Eng. 53(11), 113108 (Nov 12, 2014). doi:10.1117/1.OE.53.11.113108
History: Received July 2, 2014; Revised October 2, 2014; Accepted October 10, 2014
Text Size: A A A

Open Access Open Access

Abstract.  This work shows an analytic solution to the central moments of the angle of linear polarization (AoLP) when the linear Stokes parameters are independent and Gaussian distributed with different means but equal variance. Such a result is useful for distinguishing AoLP features from noise in polarimetry. When the DoLP is high relative to the measurement uncertainty of the linear Stokes vector, AoLP statistics have been shown to be well approximated by a Gaussian distribution. When the DoLP is zero, AoLP values are uniformly distributed. In general, the probability density function (PDF) of AoLP does not have a closed-form solution and this is the first report, to our knowledge, on an exact analytic form for the central moments of the AoLP. This analytic form will be useful when the AoLP is of interest even when the DoLP is low and the corresponding PDF on the AoLP is in between the extreme cases of a Gaussian or a uniform distribution. We also show that a simple propagation of error (PE) analysis underestimates the AoLP variance at extremely low DoLP but is verified for cases of DoLP that are high relative to the Stokes measurement uncertainty. An example use of the AoLP variance in imaging polarimetry is presented.

Figures in this Article

The Stokes vector that is limited to the three first components is commonly referred to as the linear Stokes vector1 and includes intensity and two linear Stokes parameters. From its measurement, two primary parameters of interest are calculated: the degree of linear polarization (DoLP) and the angle of linear polarization (AoLP). The AoLP and DoLP are nonlinear transforms of the two linear Stokes parameters from Cartesian to polar coordinates. The AoLP expresses the orientation of the linear polarization state. The DoLP expresses the magnitude of the linear polarization state. This work presents an exact analytic solution for the central moments of the AoLP when the linear Stokes parameters are drawn from uncorrelated Gaussian distributions of different means but equal variance. Fluctuations due to a change in the object’s polarization state during data acquisition (e.g., atmospheric turbulence) or fluctuations due to system artifacts (e.g., error from mechanical motion of optical components) could potentially introduce significant correlations. To test the validity of the uncorrelated assumption both the polarimetric instrumentation and the postprocessing calculations must be considered. If in postprocessing, different raw measurements are used to calculate each linear Stokes parameter then the measurement statistics will be uncorrelated owing to the independence of each measurement. Correlations would be introduced by reusing an individual measurement in the calculation of both linear Stokes parameters. Division of amplitude, aperture, or wavefront polarimeters traditionally use independent measurements to estimate each linear Stokes parameter. Traditional wire-grid polarizers also produce linear Stokes parameters with measurement noise that is uncorrelated since a different pair of independent measurements are used to calculate each.

If the Stokes parameters are calculated by combining a weighted sum of numerous independent measurements then, owing to the central limit theorem, a Gaussian probability density function (PDF) is expected with increasing confidence as the number of measurements used in the calculation increases.2 For example, modulated polarimeters collect a series of measurements while simultaneously changing the instrument’s polarimetric response. Using this strategy, the object’s polarization can be encoded temporally, spatially, and/or spectrally in the measurements. Then, in postprocessing numerous independent measurements are combined, usually by fitting to a model, to form an estimate of the object’s polarization.

Physical models of polarization are often specified in polar coordinates (e.g., astronomy and material scattering properties), therefore, characterizing AoLP statistics is of interest.3,4 Although an integral form of the PDF of the AoLP is known, this is the first report, to our knowledge, on an exact analytic expression for the central moments of the AoLP and their relation to the measurement uncertainty of the linear Stokes vector. When the DoLP is high relative to the measurement uncertainty of the linear Stokes vector, AoLP statistics have been shown to be well approximated by a Gaussian distribution;5 this is consistent with the results presented in this paper. As the DoLP approaches the measurement uncertainty of the linear Stoke parameters, the distribution on AoLP is no longer well approximated by a Gaussian distribution; this is also consistent with the results presented in this paper. This low DoLP limit is of interest and importance to researchers in astrophysics,5 reconnaissance,6 and atmospheric aerosol retrieval.7

Consider the two linear Stokes parameters q and u as independent and normally distributed random variables with equal variance, pr(q)=N(q¯,σu,q2) and pr(u)=N(u¯,σu,q2). Here, all quantities are normalized by the intensity (i.e., the first linear Stokes component), as in q=Q/I and u=U/I where I denotes the intensity. In practice, the PDF on q and u will depend upon how these quantities are measured. If in postprocessing the Stokes parameters are calculated from a sum, or weighted sum, of raw measurements then we can argue from the central limit theorem that the pr(u) and pr(u) will be approximately Gaussian if the raw measurements are sufficiently independent. From pr(u) and pr(u), a transformation of random variables is used to obtain pr(θ) as a function of the means and variance of the two linear Stokes parameters. Beginning with the joint PDF Display Formula

prθ,d(θ,d)=prq,u[dcos(2θ),dsin(2θ)]|det[cos(2θ)2dsin(2θ)sin(2θ)2dcos(2θ)]|=2d2πσu,q2exp{[dcos(2θ)q¯]2+[dsin(2θ)u¯]22σu,q2},(1)
where θ is the AoLP and d is the DoLP. Here, the transformations q=dcos(2θ) and u=dsin(2θ) have been used to relate the linear Stokes parameters to the AoLP and the DoLP. Marginalizing over AoLP yields the Rice distribution on DoLP.8 Marginalizing over the DoLP to obtain the PDF of AoLP Display Formula
pr(θ)=01dpprθ,d(θ,p)=1πed¯220dp˜p˜ep˜22ep˜[q¯cos(2θ)+u¯sin(2θ)],(2)
where all parameters denoted with the tilde have been normalized by σu,q, e.g., u¯˜=(u¯/σu,q) and d¯˜2=u¯˜2+q¯˜2. This integral form is the moment generating function of the Rayleigh distribution; it is not closed-form but can be expressed using error functions as in 9.

The definition of the α-order central moment of the AoLP is Display Formula

(θϕ¯)α=ϕ¯π/2ϕ¯+π/2dθpr(θ)(θϕ¯)α.(3)
Here, ϕ¯=1/2atan (u¯/q¯) is half the angle between the mean of prq,u(q,u) and the q-axis, as depicted in Fig. 1. The distribution pr(θ) is symmetric about ϕ¯ and the expectation over all possible values of θ is evaluated over a π interval centered on ϕ¯. This definition of pr(θ) avoids the spurious values in calculations of angular quantities if wrapping of the directional statistics is not accounted for.

Graphic Jump LocationF1 :

In (a) and (b), the yellow points are samples from prq,u(q,u) where the spread of this distribution is characterized by d¯˜=u¯2+q¯2/σu,q and both low and high examples of this value are provided. Two points are labeled: the mean r¯˜ and the n’th sample r˜n. ϕ¯ is the angle of linear polarization (AoLP) of r¯˜ and θn is the AoLP of r˜n. The AoLP’s central moments are related to the angular spread from ϕ¯. A histogram of (θϕ¯) is given in (c) and (d); the red line is a least-squares fit to a Gaussian distribution.

In the 1, we show that the α-order central moment of the AoLP can be equivalently expressed as Display Formula

(θϕ¯)α=12α+1πn=Bnαfn(d¯˜),(4)
where the scalar Bnα and the function fn are defined below. As expected, the central moments are circularly symmetric with respect to the average value of the AoLP, in other words, only the DoLP and the uncertainty of the linear Stokes parameters σu,q are in this expression; ϕ¯ is not. The term that can be factored as independent of the DoLP and only dependent upon the order of the moment is Display Formula
Bnα=12πππβαeinβdβ,(5)
where β=2θ2ϕ¯ describes the angular deviation as depicted in Fig. 1. For the first and second moments, which are calculated in Secs. 3 and 4Display Formula
Bn1={(1)ninn00otherwise,(6)
Display Formula
Bn2={2(1)nn2n0.1cmπ23otherwise.(7)
The other term depends on the DoLP and the linear Stokes uncertainty, but not the α moment Display Formula
fn(d¯˜)={π32ed¯24d¯[In12(d¯24)+In+12(d¯24)]n02πotherwise.(8)
Here, I is a modified Bessel function of the first kind. 1 contains a detailed derivation of Eq. (4) from Eq. (3).

The analytic expression for the AoLP central moments is an infinite sum which, in practice, must be truncated, therefore, the computed value depends upon the number of terms retained Display Formula

(θϕ¯)α(θϕ¯)α^(K)=12α+1πn=KKBnαfn(d¯˜).(9)
When performing this computation, it is useful to notice that fn=fn. The analytic form for the square root of the second central moment, i.e., the standard deviation, of the AoLP is shown in Fig. 2(b) for K=500.

Graphic Jump LocationF2 :

The AoLP standard deviation is reported in degrees and calculated from (a) sample statistics [Eq. (12)], (b) the new analytic form presented in this work [Eq. (9)], and (c) the propagation of error approximation [Eq. (11)]. The maximum differences between methods (b) and (c) with respect to method (a) are 0.03 deg and 4.72 deg, respectively. This agreement between an asymptotically unbiased estimate and the analytic form is expected and verifies the precision of the analytic expression and its derivation presented in the 1.

To verify the precision of the new analytic expression for the central moments of the AoLP, the square root of the second central moment (i.e., α=2) is compared to two other numerical methods for calculating the AoLP standard deviation.

When a random variable of unknown variance has a nonlinear relationship to random variables of known variance the propagation of error (PE) method approximates the unknown variance by a truncated Taylor series as described in 10. When the linear Stokes parameters are uncorrelated, the PE approximation of the AoLP variance is Display Formula

σθ2σu2(θu)2+σq2(θq)2.(10)
By applying the assumption of equal variances for each linear Stokes parameter, the PE approximation becomes Display Formula
σ^PEθ=12d˜=σu,q2d,(11)
which is displayed in Fig. 2(c).

The AoLP variance can also be calculated directly from Nsamp linear Stokes parameter samples Display Formula

σ^θ2(Nsamp)=1Nsampn=1Nsamp(θnϕ¯)2=1Nsampn=1Nsamp[12acos(r¯·rn|r¯||rn|)]2,(12)
where rn is the n’th sample generated by sampling from the Gaussian distributions pr(q) and pr(u), the mean values of these distributions are r¯=[u¯,q¯], and the dot product r¯·rn is used to calculate half of the angle between the vectors (θnϕ¯) (depicted in Fig. 1). This sample variance is asymptotically unbiased, therefore, it will approach the true value of the AoLP variance as Nsamp increases. To verify the precision of the analytic form, a true value for the AoLP variance is needed. To provide this comparison, the sample variance is calculated using a very large number of samples. The calculation using Nsamp=106 is displayed in Fig. 2(a). The analytic form is closer to the sample variance than the PE approximation; the maximum difference is given as a text box in Fig. 2 and is 0.03 deg for the analytic form and 4.72 deg for the PE approximation.

To investigate the sensitivity between number of terms retained in Eq. (9) and the AoLP variance, a comparison is shown in Fig. 3(a). The number of terms required for convergence depends upon the AoLP variance. Accuracy at low AoLP variances (i.e., high DoLP) requires retaining more terms. For low DoLP environments and σθ>10deg, convergence is achieved at even 20 terms. This convergence behavior is advantageous since the analytic expression for AoLP variance is most useful in low DoLP environments. In relatively high DoLP environments, the analytic expression is not necessary because simpler approximations for the variance of AoLP, such as the PE method, are valid.

Graphic Jump LocationF3 :

(a) To show the effect of truncating the infinite sum in Eq. (4), the calculation is compared for different numbers of terms retained. Convergence with respect to K is slowest at low values of σθ. The results become indistinguishable between K=250 and K=500. The black line is σ^PEθ and does not agree with the sample statistics or the analytic form at high values of σθ. (b) The under estimation of σθ by σ^PEθ is shown to occur at DoLP values below 2% when σu,q=0.005.

Figure 3(b) is a comparison between the PE estimate and the analytic expression for the AoLP variance at low DoLP. For σu,q=0.005, the two methods agree above about 2% DoLP. For DoLP lower than 2%, the methods begin to disagree and the PE analysis predicts lower values for σθ. At 1% DoLP, the difference between the two methods is approximately 5 deg. This difference is potentially significant because statistical testing becomes most useful when the DoLP is low.

Figure 4 shows a cloud measurement from the Ground-based Multiangle SpectroPolarimetric Imager (GroundMSPI).11 This is an ideal application for our analytic AoLP variance since 47% of the pixels in the scene are at or below 2% DoLP; see Fig. 4(b). AoLP statistical analysis can be used to test GroundMSPI’s sensitivity in distinguishing single from multiple scattering events. In the case of single scattering, the AoLP of the cloud will either be the same as the AoLP of the sky or 90 deg from it as described in Refs. 12 and 13. Given GroundMSPI’s measurement uncertainty in the Stokes parameters of ±0.005,14 the analytic form presented in this paper is used to calculate the standard deviation of the AoLP measured at each image pixel. Pixels are selected from the GroundMSPI image based on three different σAoLP thresholds of: 3 deg, 6 deg, and 9 deg. These thresholds result in including 40%, 48%, and 60% of the image pixels; respectively. In Fig. 5, these AoLP values are plotted versus DoLP to separate sky from cloud pixels. The AoLP is plotted with error bars ±2σθ from the measured value. The ±2σθ interval has been shown in 5 to produce a 95.45% confidence level on AoLP. The AoLP uncertainty is lowest for the sky measurements since these have a higher DoLP value than the cloud [see Figs. 4(b) and 4(d)]. As the σAoLP threshold is increased, more low DoLP cloud pixels are included and a majority of these are ±90deg from the sky’s AoLP, which indicates single scattering. Multiple scattering would further decrease the DoLP, and at the lowest σAoLP threshold, more low DoLP points that are intermediate to sky and ±90deg sky AoLP appear.

Graphic Jump LocationF4 :

An 865 nm GroundMSPI cloud measurement acquired on August 16, 2013 13:27 (PST) at 32°N, 110°W. Here, the solar angle is [23 deg, 40 deg] and the view angle ranged from [37 deg to 46 deg, 172 deg to 188 deg] within the field of view (FOV). Angles are expressed as [zenith, azimuth]. The resulting scattering angle over the FOV was 145 deg to 159 deg. The measured linear Stokes components are used to calculate (b) DoLP and (c) AoLP. The assumption σq=σu=0.005 is used to calculate (d), the standard deviation of AoLP.

Graphic Jump LocationF5 :

AoLP of GroundMSPI measurements in Fig. 4 is plotted versus DoLP. This separates the sky pixels, which are more polarized, from the cloud pixels. The sky’s average AoLP is 90 deg in scattering coordinates as expected and denoted by a red line. Each AoLP point includes an error bar of ±2σAoLP. A single event scattering within the cloud results in an AoLP that equals the sky’s AoLP or differs by ±90deg.

To further investigate the presence of multiple scattering, Fig. 6(a) shows all AoLP measurements with a ±2σAoLP confidence interval that does not overlap the sky’s AoLP within a 12 deg boundary. This thresholding results in including 4% of the pixels and the location of these are displayed in Fig. 6(b). A higher density of intermediate AoLP values are seen as the DoLP decreases which would be expected for the AoLP of multiple scattering events.13 This analysis suggests the possibility that GroundMSPI is capable of detecting multiple scattering events and that a very low number of multiple scattering events is present in this scene but is not comprehensive enough to conclude that the intermediate population of AoLP values is caused by multiple scattering. Instead this analysis motivates further statistical AoLP studies on more polarimetric cloud images, as well as inspection of instrument calibration data to check whether the Gaussian noise assumption is valid. This cursory statistical analysis of MSPI cloud images is presented to motivate the usefulness of an analytic form for the AoLP variance when the linear Stokes parameters are uncorrelated and normally distributed with equal variance.

Graphic Jump LocationF6 :

(a) To distinguish multiple scattering, pixels with a ±2σAoLP confidence interval that do not coincide with the sky’s AoLP are displayed. The sky’s AoLP is defined as 90 deg or 180 deg within a 12 deg range. These intervals are denoted by the red horizontal lines. This thresholding results in including 4% of the pixels in the image. A higher density of intermediate AoLP values is measured as the DoLP decreases, which would be expected for the AoLP of multiple scattering events. AoLP measurements that are below the sky AoLP (plotted in blue) are more dense than above the sky AoLP (plotted in cyan). (b) The locations of intermediate AoLP values are displayed as an overlay of dots on the intensity image; blue or cyan denotes the measurements below or above the sky AoLP.

Here, we have provided an analytic form which relates the central moments of the AoLP to the Gaussian statistics of the linear Stokes vector. This expression yields a functional relationship between the measurement uncertainty in AoLP and the measurement uncertainty in DoLP. This is especially useful for high fidelity imaging polarimetry measurements of low DoLP scenes. In such an environment, simple Gaussian assumptions on the statistics of the AoLP are not valid because the high measurement precision of the AoLP statistics may contain useful information about the scene. Additionally, we have shown that a simple PE analysis underestimates the AoLP variance at extremely low DoLP values, but is accurate at DoLP values that are high relative to the Stokes measurement uncertainty. An example use of the analytic form on AoLP variance is presented for an example in imaging polarimetry.

Substituting Eq. (2) into Eq. (3) yields an expression for the AoLP central moments involving two integrals Display Formula

(θϕ¯)α=1πed¯˜220dp˜p˜ep˜22ϕ¯π/2ϕ¯+π/2dθep˜(q¯˜cos(2θ)+u¯˜sin(2θ))(θϕ¯)α.(13)
From Fig. 1, it can be seen that q¯˜=d¯˜cos(2ϕ¯) and u¯˜=d¯˜sin(2ϕ¯). Using the equality cos(A)cos(B)+sin(A)sin(B)=cos(BA) and the change of variables β=2θ2ϕ¯ allows the angular integral from Eq. (13) to be expressed as Display Formula
12α+1ππdββαep˜d¯˜cos(β).(14)
To represent βα on the [π,π] interval, the function can be expanded in a Fourier series Display Formula
βα=Bnαeinβ,(15)
where the Fourier coefficients are defined by Display Formula
Bnα=12πππdββαeinβ.(16)
This integral can be represented recursively for n0 by integrating by parts Display Formula
Bnα=(1)nπα1i2n[1(1)α]+αinBnα1(17)
and the integral can be solved for n=0Display Formula
B0α=πα+12π(α+1)[1(1)α+1].(18)
When α=0 the integral in Eq. (16) simplifies to Display Formula
Bn0=12πππdβeinβ=δn0,(19)
where δ is the Kronecker delta function. Since the Fourier coefficients have closed-form solutions, substituting Eq. (15) into Eq. (14) is helpful because it leads to an expression Display Formula
12α+1n=Bnαππdβeinβep˜d¯˜cos(β)(20)
that can be related to a modified Bessel function Display Formula
In(x)=12πππexcos(θ)einθdθ.(21)
Now the integral expression over the angle is replaced by an infinite sum involving modified Bessel functions Display Formula
2π2α+1n=BnαIn(p˜d¯˜).(22)
By substituting these results into Eq. (13), the AoLP central moments are Display Formula
(θϕ¯)α=12αed¯˜22n=Bnα0dp˜p˜ep˜22In(p˜d¯˜).(23)
Noting the recurrence relations for derivatives of modified Bessel functions Display Formula
ddxIn(x)=12In1(x)+12In+1(x)(24)
and using integration by parts from page 259 of 15 leads to the equality Display Formula
0dp˜p˜ep˜22In(p˜d¯˜)=d¯˜2π2ed¯˜24[In12(d¯˜24)+In+12(d¯˜24)]n0.(25)

For n=0, relating the modified Bessel function to a Bessel function with a purely complex argument allows a solution using integration by parts Display Formula

0dp˜p˜ep˜22I0(p˜d¯˜)=0dp˜p˜ep˜22J0(ip˜d¯˜)=ed¯˜22(26)
Substituting these integral solutions into Eq. (23) leads to Eq. (4) which is presented in the body of this paper. This analytic expression offers a solution for the central moments of the AoLP that involves an infinite series of modified Bessel functions to avoid evaluating the expected value over two integrals. This infinite series has been shown to be computationally useful because it converges relatively quickly as presented in Fig. 3(a).

This work is supported by NSF/CHE-1313892. Christine Bradley provided the GroundMSPI images and brought the application of AoLP analysis to distinguish single and multiple scattering to our attention. Discussions with Michael Garay, Anthony Davis, and David Diner of the Jet Propulsion Laboratory, California Institute of Technology improved the interpretation of GroundMSPI data.

Chipman  R. A., “Polarimetry,” in OSA Handbook of Optics. , pp. 22.21 –22.35,  McGraw-Hill ,  New York  (1995).
DasGupta  A., Asymptotic Theory of Statistics and Probability. ,  Springer ,  New York  (2008).
Clark  D., Stewart  B. G., “Statistical methods of stellar polarimetry,” Vistas Astron.. 29, (1 ), 27 –51 (1986). 0083-6656 CrossRef
Goudail  F., Bénière  A., “Estimation precision of the degree of linear polarization and of the angle of polarization in the presence of different sources of noise,” Appl. Opt.. 49, (4 ), 683 –693 (2010). 0003-6935 CrossRef
Naghizadeh-Khouei  J., Clarke  D., “On the statistical behavior of the position angle of linear polarization,” Astron. Astrophys.. 274, , 968  (1993). 0004-6361 
Coulson  K. L., “Effects of reflection properties of natural surfaces in aerial reconnaissance,” Appl. Opt.. 5, (6 ), 905 –917 (1966). 0003-6935 CrossRef
Mishchenko  M., Travis  L., “Satellite retrieval of aerosol properties over the ocean using polarization as well as intensity of reflected sunlight,” J. Geophys. Res. Atmos.. 102, (D14 ), 16989 –17013 (1997). 0148-0227 CrossRef
Rice  S. O., “Mathematical analysis of random noise,” Bell Syst. Tech. J.. 24, (1 ), 46 –156 (1945). 0005-8580 CrossRef
Quinn  J. L., “Bayesian analysis of polarization measurements,” Astron. Astrophys.. 538, , A65  (2012). 0004-6361 CrossRef
Bevington  P. R., Robinson  D. K., Data Reduction and Error Analysis for the Physical Sciences. , Vol. 2,  McGraw-Hill ,  New York  (1969).
Diner  D. et al., “Exploration of a polarized surface bidirectional reflectance model using the ground-based multiangle spectropolarimetric imager,” Atmosphere. 3, (4 ), 591 –619 (2012). 2073-4433 CrossRef
Pust  N., Shaw  J., “Digital all-sky polarization imaging of partly cloudy skies,” Appl. Opt.. 47, (34 ), H190 –H198 (2008). 0003-6935 CrossRef
Kokhanovsky  A., “Optical properties of terrestrial clouds,” Earth Sci. Rev.. 64, (3 ), 189 –241 (2004). 0012-8252 CrossRef
Mahler  A., Diner  D., Chipman  R., “Analysis of static and time-varying polarization errors in the multiangle spectropolarimetric imager,” Appl. Opt.. 50, (14 ), 2080 –2087 (2011). 0003-6935 CrossRef
Olver  F. et al., NIST Handbook of Mathematical Functions. , Vol. 256,  Cambridge University Press ,  Cambridge  (2010).

Meredith Kupinski earned a BS degree in imaging technology from the Rochester Institute of Technology and both MS and PhD degrees in optical science from University of Arizona. She is an assistant research professor at the University of Arizona’s College of Optical Sciences working on statistical optics, imaging science, and probabilistic decision theory. She is the recipient of a Science, Engineering, and Education for Sustainability (SEES) fellowship from the National Science Foundation.

Russell Chipman received his BS degree from MIT and his MS/PhD degree in optical science from the University of Arizona. He is the co-PI of the Jet Propulsion Labs Multi-Spectral Polarimetric Imager, a candidate instrument for the global monitoring of aerosols from space to support measurements of the global energy balance for climate studies. He is a fellow of both OSA and SPIE and the 2007 recipient of SPIE’s Stokes award for polarimetry.

Eric Clarkson obtained a BA degree in mathematics, physics, and philosophy from Rice University, an MS degree in physics and a PhD degree in mathematics from Arizona State University, and an MS degree in optical sciences from the University of Arizona. He is a professor at the University of Arizona with joint appointments in the Department of Medical Imaging and the College of Optical Sciences. His research interests include statistical decision theory, statistical estimation, information theory, and image analysis.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Meredith Kupinski ; Russell Chipman and Eric Clarkson
"Relating the statistics of the angle of linear polarization to measurement uncertainty of the Stokes vector", Opt. Eng. 53(11), 113108 (Nov 12, 2014). ; http://dx.doi.org/10.1117/1.OE.53.11.113108


Figures

Graphic Jump LocationF1 :

In (a) and (b), the yellow points are samples from prq,u(q,u) where the spread of this distribution is characterized by d¯˜=u¯2+q¯2/σu,q and both low and high examples of this value are provided. Two points are labeled: the mean r¯˜ and the n’th sample r˜n. ϕ¯ is the angle of linear polarization (AoLP) of r¯˜ and θn is the AoLP of r˜n. The AoLP’s central moments are related to the angular spread from ϕ¯. A histogram of (θϕ¯) is given in (c) and (d); the red line is a least-squares fit to a Gaussian distribution.

Graphic Jump LocationF2 :

The AoLP standard deviation is reported in degrees and calculated from (a) sample statistics [Eq. (12)], (b) the new analytic form presented in this work [Eq. (9)], and (c) the propagation of error approximation [Eq. (11)]. The maximum differences between methods (b) and (c) with respect to method (a) are 0.03 deg and 4.72 deg, respectively. This agreement between an asymptotically unbiased estimate and the analytic form is expected and verifies the precision of the analytic expression and its derivation presented in the 1.

Graphic Jump LocationF3 :

(a) To show the effect of truncating the infinite sum in Eq. (4), the calculation is compared for different numbers of terms retained. Convergence with respect to K is slowest at low values of σθ. The results become indistinguishable between K=250 and K=500. The black line is σ^PEθ and does not agree with the sample statistics or the analytic form at high values of σθ. (b) The under estimation of σθ by σ^PEθ is shown to occur at DoLP values below 2% when σu,q=0.005.

Graphic Jump LocationF4 :

An 865 nm GroundMSPI cloud measurement acquired on August 16, 2013 13:27 (PST) at 32°N, 110°W. Here, the solar angle is [23 deg, 40 deg] and the view angle ranged from [37 deg to 46 deg, 172 deg to 188 deg] within the field of view (FOV). Angles are expressed as [zenith, azimuth]. The resulting scattering angle over the FOV was 145 deg to 159 deg. The measured linear Stokes components are used to calculate (b) DoLP and (c) AoLP. The assumption σq=σu=0.005 is used to calculate (d), the standard deviation of AoLP.

Graphic Jump LocationF5 :

AoLP of GroundMSPI measurements in Fig. 4 is plotted versus DoLP. This separates the sky pixels, which are more polarized, from the cloud pixels. The sky’s average AoLP is 90 deg in scattering coordinates as expected and denoted by a red line. Each AoLP point includes an error bar of ±2σAoLP. A single event scattering within the cloud results in an AoLP that equals the sky’s AoLP or differs by ±90deg.

Graphic Jump LocationF6 :

(a) To distinguish multiple scattering, pixels with a ±2σAoLP confidence interval that do not coincide with the sky’s AoLP are displayed. The sky’s AoLP is defined as 90 deg or 180 deg within a 12 deg range. These intervals are denoted by the red horizontal lines. This thresholding results in including 4% of the pixels in the image. A higher density of intermediate AoLP values is measured as the DoLP decreases, which would be expected for the AoLP of multiple scattering events. AoLP measurements that are below the sky AoLP (plotted in blue) are more dense than above the sky AoLP (plotted in cyan). (b) The locations of intermediate AoLP values are displayed as an overlay of dots on the intensity image; blue or cyan denotes the measurements below or above the sky AoLP.

Tables

References

Chipman  R. A., “Polarimetry,” in OSA Handbook of Optics. , pp. 22.21 –22.35,  McGraw-Hill ,  New York  (1995).
DasGupta  A., Asymptotic Theory of Statistics and Probability. ,  Springer ,  New York  (2008).
Clark  D., Stewart  B. G., “Statistical methods of stellar polarimetry,” Vistas Astron.. 29, (1 ), 27 –51 (1986). 0083-6656 CrossRef
Goudail  F., Bénière  A., “Estimation precision of the degree of linear polarization and of the angle of polarization in the presence of different sources of noise,” Appl. Opt.. 49, (4 ), 683 –693 (2010). 0003-6935 CrossRef
Naghizadeh-Khouei  J., Clarke  D., “On the statistical behavior of the position angle of linear polarization,” Astron. Astrophys.. 274, , 968  (1993). 0004-6361 
Coulson  K. L., “Effects of reflection properties of natural surfaces in aerial reconnaissance,” Appl. Opt.. 5, (6 ), 905 –917 (1966). 0003-6935 CrossRef
Mishchenko  M., Travis  L., “Satellite retrieval of aerosol properties over the ocean using polarization as well as intensity of reflected sunlight,” J. Geophys. Res. Atmos.. 102, (D14 ), 16989 –17013 (1997). 0148-0227 CrossRef
Rice  S. O., “Mathematical analysis of random noise,” Bell Syst. Tech. J.. 24, (1 ), 46 –156 (1945). 0005-8580 CrossRef
Quinn  J. L., “Bayesian analysis of polarization measurements,” Astron. Astrophys.. 538, , A65  (2012). 0004-6361 CrossRef
Bevington  P. R., Robinson  D. K., Data Reduction and Error Analysis for the Physical Sciences. , Vol. 2,  McGraw-Hill ,  New York  (1969).
Diner  D. et al., “Exploration of a polarized surface bidirectional reflectance model using the ground-based multiangle spectropolarimetric imager,” Atmosphere. 3, (4 ), 591 –619 (2012). 2073-4433 CrossRef
Pust  N., Shaw  J., “Digital all-sky polarization imaging of partly cloudy skies,” Appl. Opt.. 47, (34 ), H190 –H198 (2008). 0003-6935 CrossRef
Kokhanovsky  A., “Optical properties of terrestrial clouds,” Earth Sci. Rev.. 64, (3 ), 189 –241 (2004). 0012-8252 CrossRef
Mahler  A., Diner  D., Chipman  R., “Analysis of static and time-varying polarization errors in the multiangle spectropolarimetric imager,” Appl. Opt.. 50, (14 ), 2080 –2087 (2011). 0003-6935 CrossRef
Olver  F. et al., NIST Handbook of Mathematical Functions. , Vol. 256,  Cambridge University Press ,  Cambridge  (2010).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

PubMed Articles
Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.