Next Article in Journal
Proactive and Reactive Transmission Power Control for Energy-Efficient On-Body Communications
Next Article in Special Issue
Integration of GPS Precise Point Positioning and MEMS-Based INS Using Unscented Particle Filter
Previous Article in Journal
Study of the Interaction of Trastuzumab and SKOV3 Epithelial Cancer Cells Using a Quartz Crystal Microbalance Sensor
Previous Article in Special Issue
Tightly Coupled Integration of Ionosphere-Constrained Precise Point Positioning and Inertial Navigation Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Evaluation of Skylight Polarization Patterns for Navigation

College of Mechatronic Engineering and Automation, National University of Defense Technology, Changsha 410073, Hunan, China
*
Author to whom correspondence should be addressed.
Sensors 2015, 15(3), 5895-5913; https://doi.org/10.3390/s150305895
Submission received: 23 October 2014 / Revised: 15 December 2014 / Accepted: 28 February 2015 / Published: 10 March 2015
(This article belongs to the Special Issue Inertial Sensors and Systems)

Abstract

:
Skylight polarization provides a significant navigation cue for certain polarization-sensitive animals. However, the precision of the angle of polarization (AOP) of skylight for vehicle orientation is not clear. An evaluation of AOP must be performed before it is utilized. This paper reports an evaluation of AOP of skylight by measuring the skylight polarization patterns of clear and cloudy skies using a full-sky imaging polarimetry system. AOP measurements of skylight are compared with the pattern calculated by the single-scattering Rayleigh model and these differences are quantified. The relationship between the degree of polarization (DOP) and the deviation of AOP of skylight is thoroughly studied. Based on these, a solar meridian extracted method is presented. The results of experiments reveal that the DOP is a key parameter to indicate the accuracy of AOP measurements, and all the output solar meridian orientations extracted by our method in both clear and cloudy skies can achieve a high accuracy for vehicle orientation.

1. Introduction

Navigation is a critical and difficult issue in many applications for autonomous vehicles. Many animals can utilize the natural polarization patterns for navigation. It has provided inspiration for us. The desert ant Cataglyphis bicolor, for example, is able to utilize the skylight polarization pattern as a compass to forage for meters and return back to its nest on a straight line [1]. A range of marine animals, such as fishes and mollusks, can perceive the polarization of light, and use it for orientation and navigation [2]. Furthermore, dung-beetles can perceive the moonlight polarization patterns [3]. “Learning from nature” is a sagacious idea to solve navigation problems [4]. To date, there is rich literature about bio-inspired navigation sensors detecting the skylight polarization pattern [4,5,6,7,8,9,10]. In ground vehicle orientation, for example, a polarized light sensor measures the angle of polarization (AOP) of skylight which is perpendicular to the plane of scattering determined by the observer, the celestial point observed and the sun, according to the single-scattering Rayleigh model [11,12]. Then, we can obtain the vehicle direction with respect to the position of the sun. Combining this information with the azimuth of the sun, the vehicle heading can be determined. However, depending on the solar elevation and the cloud cover, the skylight polarization pattern differs more or less from the theoretical pattern calculated by the single-scattering Rayleigh atmosphere [7,11,12,13,14]. The error can significantly reduce the accuracy of the heading solution. Thus, the celestial polarization pattern must be evaluated prior to using of the AOP measurement of skylight.
In practice, sky polarization patterns are not straightforward because of the complex and varying atmospheric conditions. Numerous theoretical and experimental investigations have focused on the subject [7,13,14,15,16,17]. Additionally, the underwater polarization patterns of skylight have been researched in many literatures [18,19,20,21,22,23]. Clear sky polarization patterns are generally consistent with the results calculated by the single-scattering Rayleigh theory apart from regions near the sun, anti-sun and neutral points [11,13]. Yet, in cloudy skies, the measured degree of polarization (DOP) can be largely reduced by the clouds compared with that of clear skies, and the AOP of skylight will also suffer some distortions when clouds are present [11,12,13]. In this case, the DOP pattern differs greatly from the single-scattering Rayleigh theory, and it is unreliable as orientation cue compared with the AOP pattern [13,24,25]; the AOP may be either parallel to the scattering plane (i.e., 90° from the Rayleigh sky case), in the same direction as the clear sky AOP, or undefined [12]. Furthermore, when the sun is hidden by thicker clouds, both DOP and AOP exhibit departures from that of the clear sky [11,12]. Consequently, patterns of AOP of skylight are not always perpendicular to the scattering plane, which is a thorny point for navigation. We must quantify these differences between AOP measurements and theoretical (single-scattering Rayleigh) patterns. However, quantitative analysis of these differences is scarce for vehicle orientation.
It is noticed that the aforementioned AOP patterns share a common characteristic that a deviated AOP is often accompanied with a low DOP. Per definition, the DOP is the measure of the proportion of light which is actually polarized to the total incident light [6]. Later on, we will prove that the lower the DOP is, the larger the error of the AOP measurement will be. Furthermore, crickets can perceive the polarized stimuli down to a DOP level of less than 7% [26], but crickets cannot perceive the skylight polarization patterns if DOP is lower than the threshold value 5% [11].
In this paper, we measured the polarization patterns of clear and cloudy skies by a CCD camera with a fish-eye lens, and compared the results with the corresponding theoretical patterns quantitatively. Especially, we conducted an in-depth study of the relationship between the value of the DOP and the departure of AOP. Based on the evaluation results, a solar meridian extracted method is presented. The experiment results show that DOP is a good parameter to indicate the accuracy of AOP measurements of skylight, and the maximum error of solar meridian orientations extracted by our method is less than 0.5° even on cloudy days.

2. Materials and Methods

2.1. Instrument Description

The full-sky imaging polarimetry system is shown in Figure 1. The technique used in this paper is similar to that of [13]. The full-sky imaging polarimetry mainly consists of a CCD camera (product model: KTJ-125M) equipped with a fish-eye lens (F-number = 1.6, focal length = 1.8 mm) and a manually rotated linear polarizer (product model: LPVISE200-A, diameter = 50.8 mm) which is in front of the fish-eye lens. As a result, the angle of view decreases from 180°–108°. All instruments are mounted on an optical flat (product model: OTSB33-1), which ensure the optical axes of the CCD camera, fish-eye lens and linear polarizer parallel to each other, and vertical pointing toward the zenith during the measuring process. Because the angle of view of the lens is colossal, the four-image technique is indispensable to improve the estimation accuracy of the skylight polarization patterns [27]. In this paper, we rotated the linear polarizer to four different relative orientations (0°, 45°, 90° and 135° clockwise from the top view of the camera) to obtain four digital images for the calculation of polarization patterns. The linear polarizer is fixed in an inner ring of a lens holder. The inner ring has a scale. Hence, we can rotate it to a certain position in the outer ring freely and accurately.
Figure 1. Full-sky imaging polarimetry system.
Figure 1. Full-sky imaging polarimetry system.
Sensors 15 05895 g001

2.2. Instrument Calibration

In this paper, we focus on the evaluation of AOP and DOP patterns of skylight, and the following two characteristics of the polarimetry system are taken into consideration in the calibration: (1) the measured reduced Mueller matrix of the fish-eye lens by the method in [25]; (2) the distortion of the fish-eye lens according to the calibration toolbox for Matlab [28,29].
It has been shown that the fish-eye lens has little effect on the light polarization patterns, which can be ignored in the measurement [13]. However, the fish-eye lens causes serious angle distortion on the incident light [29]. Consequently, we mainly calibrate the angle distortion caused by the fish-eye lens and the calibration results are listed in Table 1. The image size is 1280 × 960 pixels.
Table 1. The calibration results.
Table 1. The calibration results.
Calibration Parameters (Pixel)Values
Focal length fc[475.87278, 475.91533]
Principal point Cc[643.71206, 477.27085]
Skew coefficient αc0.00047
Distortion coefficients kc[−0.26688, 0.07423, 0.00011, 0.00023, −0.00940]

2.3. Polarization Patterns Estimation and Analysis

In order to analyze the skylight polarization patterns, we calculate the AOP and DOP based on the Stokes polarization parameters. In general, light in the atmosphere is mostly linearly polarized, and the circularly polarized light is usually negligibly small [15,30], so the Stokes’ famous intensity formula can be simplified as Equation (1) [31].
I ( α ) = 1 2 [ I + Q cos 2 α + U sin 2 α ]
where I , Q and U are Stokes parameters. I ( α ) is the intensity measured at a polarizer orientation angle α. To improve the estimation accuracy of the skylight polarization patterns, the three Stokes parameters are calculated by rotating the transmission axis of the polarizer to the angles α = 0°, 45°, 90° and 135°, respectively [16]. They can be obtained by Equation (2).
[ I Q U ] = [ ( I ( 0 ° ) + I ( 45 ° ) + I ( 90 ° ) + I ( 135 ° ) ) / 2 I ( 0 ° ) I ( 90 ° ) I ( 45 ° ) I ( 135 ° ) ]
Then, the resulting degree of (linear) polarization P and its angle ψ are computed as Equations (3) and (4) [31]:
In general, the measurement errors of AOP patterns result from sensor manufacturing (scale factors, and offsets, etc.) and calibration procedure. However, we shall mathematically prove that the error of the AOP is also affected by the value of DOP. We therefore perform the variation operation on both sides of Equation (4). Then, we can get Equation (5):
P = Q 2 + U 2 I
ψ = 1 2 arctan | U Q |
δψ = U 2 ( Q 2 + U 2 ) δ Q + Q 2 ( Q 2 + U 2 ) δ U
Let U = IP·sin 2ψ and, Q = IP·cos 2ψ so a simple form of Equation (5) is given by Equation (6).
δψ = sin 2 ψ 2 I P δ Q + cos 2 ψ 2 I P δ U
Note that the amplitudes of the coefficients of δ Q and δU decrease monotonically with the value of P for a given ψ and I; that is, other things being equal, a low DOP will result in a poor AOP measurement.

2.4. Using the Angle of Polarization (AOP) Patterns for Orientation

The polarization patterns of skylight can be measured by two methods: point-source based and image based [32]. The image-based method can obtain the polarization information of a region of the sky. Accordingly, we can use the AOP patterns of this region of the sky for orientation. According to the single-scattering Rayleigh model, one important feature of polarization information of a region with a center at the zenith is that the polarization pattern of skylight is symmetrical along the solar-antisolar meridian [13,25,33]. Hence, we can estimate the orientation of the solar meridian by detecting the symmetry axis of the polarization pattern of skylight. Combining the solar meridian orientation with the azimuth of the sun, we can robustly determine the vehicle heading.
In the orientation process, the image based polarization sensor, such as our full-sky imaging polarimetry system, is aligned with the body frame (b-frame) of the vehicle. We assume that the polarization sensor frame (s-frame) coincides with b-frame. The s-frame consists of three orthogonal axes where the xs-axis is parallel with the row, ys-axis is parallel with the column of the pixels of the CCD sensor, respectively, and zs-axis is the optical axes pointing in the upward direction. The b-frame consists of three orthogonal axes where the x-axis points in the forward direction, ys-axis is in the direction of transverse motion, and z-axis is the upward direction. We let the x and y axis be parallel with the xs and ys axis, respectively. During the measuring process, we keep the xy plane horizontal; that is, the z-axis vertical points toward the zenith.
Then, the angle between geographic north and the solar meridian in the xy plane denotes the azimuth of the sun ϕ s , the angle between the positive x axis and the solar meridian in the xy plane denotes the solar meridian orientation θ in the b-frame, and the angle between geographic north and the positive x axis denotes the vehicle heading φ Figure 2 shows the relationship among the three angles in the xy plane of the b-frame.
Therefore, the vehicle heading φ can be determined by the following equation. We perform the variation operation on both sides of Equation (7), and we can get Equation (8).
φ = ϕ s θ
where ϕ s can be obtained by the solar position calculation algorithm in [34] with an accuracy of 0.0083°, and θ can be solved by the AOP patterns measured by the polarization sensor.
δφ = δ ϕ s δθ
Neglecting the small value δΦs, Equation (8) can be rewritten as Equation (9).
δφ δθ
Note that the accuracy of vehicle heading φ strongly depends on the solution precision of the solar meridian orientation θ. Hence, it is necessary to solve θ with a high accuracy.
Figure 2. The relationship between the azimuth of the sun, the solar meridian orientation, and the vehicle heading in the horizontal plane.
Figure 2. The relationship between the azimuth of the sun, the solar meridian orientation, and the vehicle heading in the horizontal plane.
Sensors 15 05895 g002

2.5. Solar Meridian Extracted from the Angle of Polarization (AOP) Patterns

Let ψ(x, y) be the AOP patterns of skylight measured by our polarimetry system, defined on a circle of radius L centered at the origin. Suppose t is the axis of symmetry of ψ(x, y) which passes through the origin and makes an angle θ with the positive x axis (Figure 3). It can be calculated as Equation (10).
Figure 3. The coordinate system for the symmetry measurement.
Figure 3. The coordinate system for the symmetry measurement.
Sensors 15 05895 g003
t = x cosθ y sinθ s = x sinθ + y cosθ
The measurements are considered to be reflectionally symmetric if ψ(t, s) = ψ(t, −s) and antisymmetric if ψ(t, s) = −ψ(t, −s) [35]. Evidently, the AOP measurements ψ(x, y) is antisymmetric as expected from Rayleigh theory [29].
A measure of symmetry with respect to the t axis is computed as Equation (11) [35]:
S θ { ψ } = 1 2 ψ ( t , s ) ψ ( t , s ) d s d t ψ 2 ( t , s ) d s d t + 1 2
where S θ { ψ } [ 0 , 1 ] . Finding the minimum of S θ { ψ } parameterized by the axis t, we will get the only axis of symmetry of ψ(x, y), and this axis of symmetry is the orientation of solar and antisolar meridian.
Note that the AOP patterns may suffer some distortions, for instance, when clouds appear. This has negative effects on the calculation of the measure of symmetry Sθ {ψ}. In order to solve this problem, we set all the AOP patterns to zero when the corresponding AOP deviations are large. The experimental results in Section 4.1 will show that a deviated AOP is often accompanied with a low DOP. Hence, we consider the AOP deviations are large when the corresponding DOP pattern is less than a DOP threshold. Pthreshold Note that crickets cannot perceive the skylight polarization patterns if DOP is lower than the threshold value 5% [11]. We use Pthreshold = 5% for our solar meridian extraction algorithm. Furthermore, the AOP measurements with the same distance from the origin as these points have also been set to zero.
According to the Rayleigh theory, the absolute values of AOP measurements on the solar and antisolar meridian are equal to 90° [29]. Calculating the average slope of those points that satisfy the condition | ψ ( x , y ) | 90 ° , we can obtain the rough value of the solar meridian orientation θ0. The solar meridian extraction method in [36] only uses the AOP measurements on the solar and antisolar meridian. If those points are distorted by the clouds, it will fail to obtain the solar meridian orientation.
In order to obtain a more stable solar meridian orientation, a refinement method is presented in this paper. The procedures of our method are described as follows:
Step 1. Calculate the AOP measurements ψ(x, y) and the DOP measurements P(x, y) using Equations 2 and 3.
Step 2. Set the AOP measurements to zero when the corresponding AOP deviations are large. The experimental results in Section 4.1 will show that a deviated AOP is often accompanied with a low DOP. Hence, we consider the AOP deviations are large when the corresponding DOP pattern is less than a DOP threshold Pthreshold. We use Pthreshold = 5% in the experiment.
Step 3. Set the AOP measurements with the same distance from the origin as these points in Step 2 to zero.
Step 4. Extract solar–antisolar meridian points in the AOP measurements ψ ( x , y ) , where the absolute values of AOP are very close to 90°. Calculate the average slope of these points as the rough value of the solar meridian orientation θ0.
Step 5. Calculate S θ { ψ } for different θ values (θ1, θ2, … , θn) using Equation (11), where the value range of θ is set from θ 0 ε to θ 0 + ε . ε is a small positive number and n is a positive integer. We use ε = 2° and n = 10 in the experiment. Then, we can obtain a S θ { ψ } series ( S θ 1 { ψ } , S θ 2 { ψ } , , S θ n { ψ } ) .
Step 6. θ i , i [ 1 , n ] can be selected as the accurate solar meridian orientation if and only if S θ i { ψ } is the minimum value in the series ( S θ 1 { ψ } , S θ 2 { ψ } , , S θ n { ψ } ) .
For convenience, we refer to this solar meridian extraction method as “our method”, and refer to the method not containing steps 2 and 3 as “our method without pretreatment” for distinguishing between them.

3. Simulation Studies

In this section, simulations have been made to assess the performance of the proposed solar meridian extraction methods and compare them with the solar meridian extraction method in [36]. For convenience, we call this method the “fitting based method” in the following text. Assume that the solar elevation is 50° and the solar meridian orientation in the b-frame is −30°. Each map contains approximately 125,600 pixels. In order to simulate the noise involved in the AOP measurements, zero-mean white Gaussian noise is added to the simulation data. The standard variance of noise σ is set to 0, 1, 2 and 4 degrees, respectively. Through Monte Carlo simulations, we evaluate the performance of the algorithm by computing the mean error.
Figure 4 contains the plots of the estimates for different methods when the standard variance of noise σ is 4 degree. The simulation results are summarized in Table 2. It can be seen from the table that our method performs slightly better than the fitting based method. All the methods can obtain stable estimates even under the condition that the standard variance of noise σ is 4 degrees in the AOP measurements. It is noticed that the fitting based method only use the AOP measurements on the solar and antisolar meridian. If those points are distorted by the clouds, the solar meridian orientation will not be obtained.
Figure 4. The solar meridian orientations extracted by (a) fitting based method; (b) our method without pretreatment and (c) our method when the standard variance of noise σ is 4 degrees.
Figure 4. The solar meridian orientations extracted by (a) fitting based method; (b) our method without pretreatment and (c) our method when the standard variance of noise σ is 4 degrees.
Sensors 15 05895 g004aSensors 15 05895 g004b
Table 2. The error of the solar meridian orientation extracted by the three methods
Table 2. The error of the solar meridian orientation extracted by the three methods
Measurement NoiseCalculation MethodsMean Error (°)
No noise Fitting based method0.18
Our method without pretreatment0.15
Our method0.06
σ = 1°Fitting based method0.20
Our method without pretreatment0.32
Our method0.16
σ = 2°Fitting based method0.23
Our method without pretreatment0.36
Our method0.19
σ = 4°Fitting based method0.35
Our method without pretreatment0.37
Our method0.20

4. Experiment Results

4.1. Skylight Polarization Measurements and Evaluation

Experiments were carried out in an outdoor environment of Changsha. Clear sky polarization patterns were performed on 8 July 2014 at 19:13, and partly cloudy sky polarization patterns were performed on 10 July 2014 at 18:50 and 19:10 in the same location. The solar elevations are 2.27°, 6.87°, and 2.81°, respectively. They were calculated by the solar position calculation algorithm in [34]. We refer to the three sky conditions simply as “clear sky”, “cloudy sky I”, and “cloudy sky II”. Patterns of DOP and AOP of both clear and cloudy skies were determined and visualized as color-coded circular maps by our system only in visible blue (470 ± 30 nm) part of the spectrum (Figure 5). Each map contains approximately 1,286,729 pixels.
We calculated polarization patterns of the Rayleigh sky at the same positions of the sun as those in clear and cloudy skies. Patterns of Rayleigh sky corresponding to the clear sky are shown in Figure 6. Using the measured AOP pattern of the sky in the visible blue part of the spectrum range and at three different sun positions, we calculated the AOP deviation Δψ = | ψ meas ψ Rayleigh | at every celestial point between the real and single-scattering Rayleigh skies, where ψ meas is the AOP of skylight measured by our full-sky imaging polarimetry, and ψ Rayleigh is the AOP of skylight calculated on the basis of the single-scattering Rayleigh model. The AOP deviations of three sky conditions are shown in row (iv) of Figure 5.
Comparison of the measured AOP patterns (row (ii) of Figure 5) and theoretical patterns (Figure 6a) indicate that, apart from regions near the sun, anti-sun, and underneath clouds, the single-scattering Rayleigh model describe the characteristics of the sky polarization patterns relatively well. The measured AOP patterns of the clear sky follow almost the same patterns of the Rayleigh model (row (iv) of Figure 5a).
It is noticed that the degree of polarization under cloudy skies are less than that under the clear sky (row (iii) of Figure 5). The clouds largely distort the normal distribution of the DOP as it occurs in the clear sky. They result in very low DOP patterns (row (iii) of Figure 5b,c). Furthermore, the AOP patterns in cloudy skies are also suffered distortion when clouds appear (row (ii) of Figure 5b,c, and row (iv) of Figure 5b,c.
The cumulative percentages of the AOP deviations of three sky conditions are shown in Figure 7. In the clear sky, 84% of the deviations are less than 4°, and 60% of them are less than 2°; In the cloudy sky I, 49% of the deviations are less than 4°, and 27% of them are less than 2°; In the cloudy sky II, 64% of the deviations are less than 4°, and 38% of them are less than 2°. These results indicate that even if the clouds are present, the percentage of the sky that follows the Rayleigh model with an accuracy of Δ ψ 4 ° is large. We can utilize the AOP patterns of skylight for orientation even on cloudy days.
Figure 5. Color photo, AOP, and DOP were measured by our instrument in blue part of the spectrum in (a) Clear sky; (b) Cloudy sky I; and (c) Cloudy sky II. Row (iv) is the AOP deviation Δψ = | ψ meas ψ Rayleigh | of the three conditions.
Figure 5. Color photo, AOP, and DOP were measured by our instrument in blue part of the spectrum in (a) Clear sky; (b) Cloudy sky I; and (c) Cloudy sky II. Row (iv) is the AOP deviation Δψ = | ψ meas ψ Rayleigh | of the three conditions.
Sensors 15 05895 g005aSensors 15 05895 g005b
Figure 6. Patterns of skylight calculated by the single-scattering Rayleigh model (the position of the sun is the same as that in the clear sky). (a) AOP of Rayleigh sky; (b) DOP of Rayleigh sky.
Figure 6. Patterns of skylight calculated by the single-scattering Rayleigh model (the position of the sun is the same as that in the clear sky). (a) AOP of Rayleigh sky; (b) DOP of Rayleigh sky.
Sensors 15 05895 g006
However, it is undeniable that a certain percentage of data points are “non-Rayleigh” points, such as the deviations of AOP are greater than 4°. Using the “non-Rayleigh” points for orientation can significantly decrease the accuracy of the heading solution.
It is noteworthy that the most striking AOP differences between the actual and the theoretical patterns are often present in the regions near the sun, anti-sun, neutral points and underneath clouds. These regions share a common characteristic that the degree of polarization is very low.
Figure 7. The cumulative percentages of the AOP deviations of the three sky conditions.
Figure 7. The cumulative percentages of the AOP deviations of the three sky conditions.
Sensors 15 05895 g007
Figure 8 shows the distribution of the AOP deviations in different intervals of DOP patterns. The proportions of the AOP deviations Δ ψ 4 ° in each DOP interval are given on the top of the bars in the bar graphs. Analyzing the figure, we can establish the following trends:
(1) The percentage of the sky that follows the Rayleigh model with accuracy of Δ ψ 4 ° or Δ ψ 2 ° in high DOP intervals is larger than that in low DOP intervals.
(2) In various DOP intervals, the higher the DOP, the larger the percentage of the sky polarization pattern that follows the Rayleigh model in this DOP interval will be. In the clear sky, the proportion increases from 26%–100% as the DOP increases from 0 to its highest value. In the cloudy sky I, the proportion increases from 22%–90% as the DOP increases from 0 to its highest value, and from 28%–92% in the cloudy sky II as the DOP increases from 0 to its highest value. So, it will be more likely to detect the accurate AOP patterns of skylight in the region with high DOP patterns than that in the region with low DOP patterns.
Figure 8. The distribution of the AOP deviations in different DOP intervals. (a) Clear sky; (b) Cloudy sky I; (c) Cloudy sky II.
Figure 8. The distribution of the AOP deviations in different DOP intervals. (a) Clear sky; (b) Cloudy sky I; (c) Cloudy sky II.
Sensors 15 05895 g008aSensors 15 05895 g008b
It can be seen from Figure 9 that the AOP deviations decrease evidently as the DOP increase from 0 to its highest value for all the three sky conditions. In the clear sky, the mean deviations of AOP can reach to 1.5° when the value of DOP is greater than 40%, and the cumulative percentage is 42%. In the cloudy sky I, the mean deviations of AOP can reach to 2.4° when the value of DOP is greater than 21%, and the cumulative percentage is 18%. In the cloudy sky II, the mean deviations of AOP are less than 2.0° when the value of DOP is greater than 26%, and the cumulative percentage is 26%.
These results show that even if in cloudy days, we could use the AOP patterns of skylight for orientation with a high accuracy by selectively utilizing the AOP patterns in the region with high DOP patterns.
Figure 9. The relationships between deviations of AOP and the values of DOP (row (ii) is the drawing of partial enlargement of row (i) in the red boxes).
Figure 9. The relationships between deviations of AOP and the values of DOP (row (ii) is the drawing of partial enlargement of row (i) in the red boxes).
Sensors 15 05895 g009aSensors 15 05895 g009b

4.2. Solar Meridian Orientation Calculation

We experimentally evaluate the performance of the solar meridian extraction algorithm with the AOP patterns of all the three sky conditions. Each AOP pattern contains approximately 1,286,729 pixels. The data is too large to be handled. Thus, we select the data in a circle of radius 100 pixels centered at the origin for the extraction algorithm.
The experimental results in Section 4.1 have showed that a deviated AOP is often accompanied with a low DOP. Hence, we set the AOP patterns to zero when the corresponding DOP pattern is less than a DOP threshold P threshold. Also, the AOP measurements at the points with the same distance from the origin as those points have been set to zero. We use P threshold. = 5% for our solar meridian extraction algorithm.
Both our method without pretreatment and the fitting based method are used for comparison. The reference heading φ is determined by an inertial navigation system (INS). It is able to yield the azimuth up to 0.01° accuracy in 10 minutes. Then, the reference solar meridian orientation θ in the b-frame can be obtained by θ   =   ϕ s φ . In the experiment, the reference solar meridian orientations of the three sky conditions are −1.50°, −1.62° and −1.92°, respectively.
The output solar meridian orientations extracted by the three methods of the three sky conditions are shown in Figure 10. Table 3 shows the result of comparison of the solar meridian orientation errors of the three methods. We can see that our method can take full advantage of all the useful AOP measurements for the determination of the solar meridian orientations, and performs slightly better than the other two methods. All the methods can obtain stable estimates even when the clouds have great influence on the AOP patterns.
Figure 10. The AOP and DOP patterns in a circle of radius 100 pixels centered at the origin. (a) Clear sky, (b) Cloudy sky I, and (c) Cloudy sky II. (The white lines in row (ii) are the solar meridian orientations extracted by fitting based method; the white lines in row (iii) are the solar meridian orientations extracted by our method without pretreatment; the white lines in row (iv) are the solar meridian orientations extracted by our method; the black crosses in row (iv) indicate the points set to zero).
Figure 10. The AOP and DOP patterns in a circle of radius 100 pixels centered at the origin. (a) Clear sky, (b) Cloudy sky I, and (c) Cloudy sky II. (The white lines in row (ii) are the solar meridian orientations extracted by fitting based method; the white lines in row (iii) are the solar meridian orientations extracted by our method without pretreatment; the white lines in row (iv) are the solar meridian orientations extracted by our method; the black crosses in row (iv) indicate the points set to zero).
Sensors 15 05895 g010aSensors 15 05895 g010b
Table 3. The error of the solar meridian orientation extracted by the three methods
Table 3. The error of the solar meridian orientation extracted by the three methods
AlgorithmsClear SkyCloudy Sky ICloudy Sky II
Fitting based method0.33°0.48°0.47°
Our method without pretreatment0.27°0.38°1.86°
Our method0.27°0.38°0.41°
It takes about 0.5 seconds in a Matlab 7.8 implementation running on a standard computer with an Intel Core 2.93 GHz processor for our solar meridian extraction method.

5. Discussion

Experiments in this paper were performed in three sky conditions. The results show that even if the clouds are present, the percentage of the sky that follows the Rayleigh model with a high accuracy is large. Consequently, we can utilize the AOP patterns of skylight for orientation even in cloudy days. However, the Mie and geometric pattern would be the dominant one under the condition that the atmosphere with relatively high concentrations of micro and above sized aerosols. It is noticed that the Mie scattering are dramatically different from that for the Rayleigh and geometric scattering [12,21,23]. The polarized portion of the Rayleigh scattered components are always polarized perpendicular to the scattering plane [11,12], while the polarized component of the Mie scattering results can potentially be polarized either parallel or perpendicular to the scattering plane [12,21,23]. If some of the AOP measurements are parallel to the scattering plane which is determined by the observer, the celestial point observed and the sun, the relationship between the AOP measurement and the azimuth of the sun would be uncertain. Then, it will fail to use the AOP measurement combining with the solar azimuth for orientation. Fortunately, it would be Rayleigh instead of Mie scattering even in turbid atmosphere when we use sensors to receive the polarized light that are sensitive at wavelengths longer than that which is visible [37]. It is worth researching for our future work.
Actually, the experiments on studying the efficiency and accuracy of utilizing the AOP patterns for vehicle orientation in the paper are static tests. We need to manually rotate the linear polarizer of the polarimetry system to four different relative orientations to obtain four digital images for the calculation of polarization patterns, and cannot collect these data in real-time. Accordingly, we cannot perform the dynamic experiments by this polarimetry system. In order to resolve this problem, many real-time image based polarization sensors with different structures and/or different work principles have been developed, such as the multichannel camera-based polarization sensor [7,10,32], the CCD polarization imaging sensor [38], and the CMOS image sensor [6]. Some outdoor experiments have been done by these sensors [10,32]. Utilizing the real-time image based polarization sensor for dynamic vehicle orientation would be desirable for our future research.

6. Conclusions

An evaluation of skylight AOP has been presented that compares the measurements of full-sky imaging polarimetry system with the theoretical results for vehicle orientation. Its accuracy is quantitatively evaluated in both clear and cloudy skies. The primary results demonstrate that we can utilize the angle of polarization patterns of skylight for vehicle orientation even in cloudy days, and the DOP is a key parameter to indicate the accuracy of the AOP measurement. Based on the evaluation results of the skylight polarization patterns, a solar meridian extracted method is presented. It can achieve a high accuracy for vehicle orientation on both clear and cloudy days.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No. 61104201), and Research Project of National University of Defense Technology (Grant No. JC14-03-04). The authors are also grateful to the anonymous reviewers for their helpful comments.

Author Contributions

Tao Ma and Xiaoping Hu conceived and designed the experiments; Yujie Wang and Zhiwen Xian performed the experiments; Lilian Zhang, Junxiang Lian and Xiaofeng He analyzed the data; Tao Ma wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Collett, M.; Collett, T.; Bisch, S.; Wehner, R. Local and global vectors in desert ant navigation. Nature 1998, 394, 269–272. [Google Scholar] [CrossRef]
  2. Shashar, N.; Johnsen, S.; Lerner, A.; Sabbah, S.; Chiao, C.-C.; Mathger, L.M.; Hanlon, R.T. Underwater linear polarization: Physical limitations to biological functions. Phil. Trans. R. Soc. B 2011, 366, 649–654. [Google Scholar] [CrossRef] [PubMed]
  3. Horváth, G.; Varju, D. Polarized light in animal vision: Polarization patterns in nature. Springer-Verlag: Berlin, Germany, 2004. [Google Scholar]
  4. Lambrinos, D.; Möller, R.; Labhart, T.; Pfeifer, R.; Wehner, R. A mobile robot employing insect strategies for navigation. Robot. Auton. Syst. 2000, 30, 39–64. [Google Scholar] [CrossRef]
  5. Chu, J.; Zhao, K.; Zhang, Q.; Wang, Tichang. Construction and performance test of a novel polarization sensor for navigation. Sensor. Actuat. A-Phys. 2008, 148, 75–82. [Google Scholar] [CrossRef]
  6. Sarkar, M.; Theuwissen, A. A Biologically Inspired cmos Image Sensor; Springer-Verlag Berlin Heidelberg: New York, NY, USA, 2013. [Google Scholar]
  7. Horvath, G.; Barta, A.; Gal, J.; Suhai, B.; Haiman, O. Ground-based full-sky imaging polarimetry of rapidly changing skies and its use for polarimetric cloud detection. Appl. Opt. 2002, 41, 543–559. [Google Scholar] [CrossRef] [PubMed]
  8. Chahl, J.; Mizutani, A. Biomimetic attitude and orientation sensors. IEEE Sensor. J. 2012, 12, 289–297. [Google Scholar] [CrossRef]
  9. Karman, S.B.; Diah, S.Z.M.; Gebeshuber, I.C. Bio-inspired polarized skylight-based navigation sensors: A review. Sensors 2012, 12, 14232–14261. [Google Scholar] [CrossRef] [PubMed]
  10. Sturzl, W.; Carey, N. A Fisheye Camera System for Polarisation Detection on Uavs. In Proceedings of Computer Vision-ECCV 2012 Workshops and Demonstrations, Florence, Italy, 7–13 October 2012; Volume 7584, pp. 431–440.
  11. Pomozi, I.; Horvath, G.; Wehner, R. How the clear-sky angle of polarization pattern continues underneath clouds: Full-sky measurements and implications for animal orientation. J. Exp. Biol. 2001, 204, 2933–2942. [Google Scholar] [PubMed]
  12. Pust, N.J.; Shaw, J.A. Digital all-sky polarization imaging of partly cloudy skies. Appl. Opt. 2008, 47, H190–H198. [Google Scholar] [CrossRef] [PubMed]
  13. Miyazaki, D.; Ammar, M.; Kawakami, R.; Ikeuchi, K. Estimating sunlight polarization using a fish-eye lens. IPSJ J. 2008, 49, 1234–1246. [Google Scholar]
  14. Suhai, B.; Horvath, G. How well does the rayleigh model describe the e-vector distribution of skylight in clear and cloudy conditions? A full-sky polarimetric study. J. Opt. Soc. Am. A 2004, 21, 1669–1676. [Google Scholar] [CrossRef]
  15. Kreuter, A.; Zangerl, M.; Schwarzmann, M.; Blumthaler, M. All-sky imaging: A simple, versatile system for atmospheric research. Appl. Opt. 2009, 48, 1091–1097. [Google Scholar] [CrossRef] [PubMed]
  16. Raymond, L.; Lee, J.; Samudio, O.R. Spectral polarization of clear and hazy coastal skies. Appl. Opt. 2012, 51, 7499–7508. [Google Scholar] [CrossRef] [PubMed]
  17. Kreuter, A.; Blumthaler, M. Feasibility of polarized all-sky imaging for aerosol characterization. Atmos. Meas. Tech. 2013, 6, 1845–1854. [Google Scholar] [CrossRef]
  18. Horvath, G.; Varju, D. Underwater refraction-polarization patterns of skylight perceived by aquatic animals through snell's wiindow of the flat water surface. Vision Res. 1995, 35, 1651–1666. [Google Scholar] [CrossRef] [PubMed]
  19. Sabbah, S.; Lerner, A.; Erlick, C.; Shashar, N. Under water polarization vision-a physical examination. Recent Res. Devel. Exp. Theor. Biol. 2005, 1, 123–176. [Google Scholar]
  20. Sabbah, S.; Barta, A.; Gál, J.; Horváth, G.; Shashar, N. Experimental and theoretical study of skylight polarization transmitted through snell’s window of a flat water surface. J. Opt. Soc. Am. A 2006, 23, 1978–1988. [Google Scholar] [CrossRef]
  21. Lerner, A.; Sabbah, S.; Erlick, C.; Shashar, N. Navigation by light polarization in clear and turbid waters. Philos. Tran. R. Soc. B: Biol. Sci. 2011, 366, 671–679. [Google Scholar] [CrossRef]
  22. You, Y.; Tonizzo, A.; Gilerson, A.A.; Cummings, M.E.; Brady, P.; Sullivan, J.M.; Twardowski, M.S.; Dierssen, H.M.; Ahmed, S.A.; Kattawar, G.W. Measurements and simulations of polarization states of underwater light in clear oceanic waters. Appl. Opt. 2011, 50, 4873–4893. [Google Scholar] [CrossRef] [PubMed]
  23. Lerner, A.; Shashar, N.; Haspel, C. Sensitivity study on the effects of hydrosol size and composition on linear polarization in absorbing and nonabsorbing clear and semi-turbid waters. J. Opt. Soc. Am. A 2012, 29, 2394–2405. [Google Scholar] [CrossRef]
  24. Brines, M.L.; Gould, J.L. Skylight polarization patterns and animal orientation. J. Exp. Biol. 1982, 96, 69–91. [Google Scholar]
  25. Gal, J.; Horvath, G.; Meyer-Rochow, V.B.; Wehner, R. Polarization patterns of the summer sky and its neutral points measured by full-sky imaging polarimetry in finnish lapland north of the arctic circle. Proc. R. Soc. Lond. A 2001, 457, 1385–1399. [Google Scholar] [CrossRef]
  26. Henze, M.J.; Labhart, T. Haze, clouds and limited sky visibility: Polarotactic orientation of crickets under difficult stimulus conditions. J. Exp. Biol. 2007, 210, 3266–3276. [Google Scholar] [CrossRef] [PubMed]
  27. Raymond, L.; Lee, J. Digital imaging of clear-sky polarization. Appl. Opt. 1998, 37, 1465–1476. [Google Scholar]
  28. Calibration Toolbox for Matlab. Available online: http://www.vision.caltech.edu/bouguetj/calib_doc/index.html (accessed on 4 March 2015).
  29. Wang, Y.; Hu, X.; Lian, J.; Zhang, L.; Xian, Z.; Ma, T. Design of a device for skylight polarization measurements. Sensors 2014, 14, 14916–14931. [Google Scholar] [CrossRef] [PubMed]
  30. Voss, K.J.; Liu, Y. Polarized radiance distribution measurements of skylight. I. System description and characterization. Appl. Opt. 1997, 36, 6083–6094. [Google Scholar] [CrossRef] [PubMed]
  31. Goldstein, D. Polarized Light; Marcel Dekker, Inc.: New York, NY, USA, 2003. [Google Scholar]
  32. Wang, D.; Liang, H.; Zhu, H.; Zhang, S. A bionic camera-based polarization navigation sensor. Sensors 2014, 14, 13006–13023. [Google Scholar] [CrossRef] [PubMed]
  33. Berry, M.V.; Dennis, M.R.; Lee, R.L., Jr. Polarization singularities in the clear sky. New J. Phys. 2004, 6, 1–14. [Google Scholar] [CrossRef]
  34. Blanco-Muriel, M.; Alarcon-Padilla, D.C.; Lopezmoratalla, T.; Lara-coira, M. Computing the solar vector. Sol. Energy 2001, 70, 431–441. [Google Scholar] [CrossRef]
  35. Kiryati, N.; Gofman, Y. Detecting symmetry in grey level images: The global optimization approach. Int. J. Comput. Vision 1998, 29, 29–45. [Google Scholar] [CrossRef]
  36. Shabayek, A.E.R. Combining Omnidirectional Vision with Polarization Vision for Robot Navigation. Ph.D. Thesis, Universite de Bourgogne, Dijon, France, 2012. [Google Scholar]
  37. Bohren, C.F.; Huffman, D.R. Absorption and Scattering of Light by Smalll Particles; John Wiley & Sons, Inc.: New York, NY, USA, 1983. [Google Scholar]
  38. Gruev, V.; Perkins, R.; York, T. CCD polarization imaging sensor with aluminum nanowire optical filters. Opt. Express 2010, 18, 19087–19094. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Ma, T.; Hu, X.; Zhang, L.; Lian, J.; He, X.; Wang, Y.; Xian, Z. An Evaluation of Skylight Polarization Patterns for Navigation. Sensors 2015, 15, 5895-5913. https://doi.org/10.3390/s150305895

AMA Style

Ma T, Hu X, Zhang L, Lian J, He X, Wang Y, Xian Z. An Evaluation of Skylight Polarization Patterns for Navigation. Sensors. 2015; 15(3):5895-5913. https://doi.org/10.3390/s150305895

Chicago/Turabian Style

Ma, Tao, Xiaoping Hu, Lilian Zhang, Junxiang Lian, Xiaofeng He, Yujie Wang, and Zhiwen Xian. 2015. "An Evaluation of Skylight Polarization Patterns for Navigation" Sensors 15, no. 3: 5895-5913. https://doi.org/10.3390/s150305895

APA Style

Ma, T., Hu, X., Zhang, L., Lian, J., He, X., Wang, Y., & Xian, Z. (2015). An Evaluation of Skylight Polarization Patterns for Navigation. Sensors, 15(3), 5895-5913. https://doi.org/10.3390/s150305895

Article Metrics

Back to TopTop