Atmospheric visibility sensor based on backscattering using correlation coding method

An new atmospheric visibility sensor is described which uses an empirical relationship between atmospheric visibility and lidar backscatter signals. Speciﬁcally, autocorrelation analysis is performed for lidar backscatter signals and the autocorrelation coefﬁcients are shown to be a function of the correlation distance. The slope of this function is then con-nected to the atmospheric visibility under different weather conditions. It is experimentally demonstrated that this new sensor can reliably estimate the atmospheric visibility in the range of 0–1,200m.


Introduction
A clear visibility of the surrounding environment is crucial in various types of transportation means. Dense fog or rain developed alongside roadways may considerably decrease the safety of motorists, thus automated detectors and monitoring of the visibility, fog, rain, snow or dust is of the upmost importance. As a feasible variable, atmospheric visibility is difficult to assess. Weather, sun angle, light intensity, darkness adaptation, availability of appropriate visibility targets and individual physical abilities determine the level of human's perception quality of atmospheric conditions. Automated visibility sensors determine the visibility range by measuring light being scattered or by utilization of extinction coefficients, locally over a short range (Siikamaki 2004;Babari et al. 2011). Optical sensors such as transmissometers or scatterometers were developed to address atmospheric visibility and visual range evaluation. The transmissometer was the first instrument developed to provide a standard for visibility measurement. The device estimates the attenuation of a light beam emitted from a source to a receiver at a certain path length and correlates obtained result with human observations. The light attenuation phenomena is due to scattering or refraction, depending on the utilized medium. However, transmissometers' physical dimensions and long baselines create a multiplicity of disadvantages when performance is considered. Demanding maintenance, high manufacturing costs and small dynamic range with insufficient accuracy discard transmissiometers from being abundantly popular. In case of the scatterometer, measurements are based on the amount of light being scattered by aerosols in an optical volume, observed within a small solid angle. The scatterometer assesses the dispersion of a light beam (Kazovsky 1985;Kemp et al. 2012), nevertheless reliability of a typical scatterometer is reported as questionable (Babari et al. 2011). Atmospheric visibility assessment is also performed with the use of light backscattering phenomenon (Taillade et al. 2008). Therein discussed is the thoroughly performed examination of the analytical model supporting backscattered luminance investigation (with an optical device geometry modeling), and as authors report, a good performance in agreement with experimental results. Another approach towards cruise supporting systems (Hautiere et al. 2008), with respect to the vehicular transportation domain, shows several visual data analysis approaches to enable estimation of the atmospheric visibility in mobile applications, with low margin of error.
In this paper, a brand new method for atmospheric visibility estimation is depicted. Presented sensor setup employs light backscattering effect with use of correlation coding technique. Evaluation of atmospheric visibility range is performed by signal correlation analysis.

Atmospheric visibility estimation based on backscattering using correlation coding technique
The principle of operation of the designed sensor is based on the analysis of backscattered optical signal. Device's transmitter emits an optical signal into the air where the light is subjected to backward scattering by air particles (for instance: rain drops, snowflakes or fog) and arrives at the receiver where it falls onto a photodetector. The backscattered optical signals within a range from 0 to 250 m are being analyzed. The principle of operation of the presented sensor is depicted on the schematic diagram in the Fig. 1. Consequently, the setup configuration of the described device is shown in the Fig. 2. Concerning architectural layout of the presented sensor, it is composed of transmitting and receiving parts. The transmitting module generates pulsed laser light with the wavelength of 808 nm and the average optical power of 250 mW (the Bob 808-2W laser is utilized). The optical signal generated consists of a 7.5 Mbit/s binary stream with 2 9 − 1 length of Pseudo Random Bit Sequence (PRBS). As presented on the Fig. 1, characteristic of emitted radiation is described as nearly parallel as it concentrates in a form of a regular beam with the diameter equal to 50 mm. Transmitting beam intersects with the receiving one at a distance of 250 m (what further implies the zone taken under consideration in autocorrelation analysis) from the surface of both active parts. In terms of the divergence, the receiving beam parameter is equal to approximately 1 of the angular spread, while the receiver area diameter is 120 mm. At the receiving end, the avalanche photodiode Hamamatsu C5331-02 (APD module on the Fig. 1) and 808 nm narrow bandpass optical filter are used. In order to sustain reliable measurements, it is of the upmost importance to prevent the receiver from falling into a saturation mode. Sensor parameters (i.e.: the laser optical power, photodiode active area size, analog to digital converter dynamic range) have been precisely selected, as it is crucial to avoid saturation state of the receiving part. To complete the description of the designed device, the transmitting part consists of the push-pull power driver and synchronization branch towards the wave generator. Receiving section utilizes HYPRES HA103-35 signal amplifier and 16-bit analog to digital converter-ADC 08L060. To organize the workflow, PIC 16F630 microchip is applied to support the control functions, as well as communication with outer devices by means of Universal Serial Bus (USB).
Further, concerning methodology insights, as it is presented on the plots within of the autocorrelation coefficient, as a function of considered distance, reflecting the outdoor atmospheric conditions (as presented by the stacked images to the right). One can see the variety of autocorrelation peaks and respective slopes of resulting falling edges. Achieved results imply relevant relationship between the slope of the autocorrelation peak (graphical assessment presented on the Fig. 4) and atmospheric visibility. Elaborating on the correlation technique itself, signals obtained at the receiving part are processed in the modulo operation implemented on the mesh of correlators (in count and length matching PRBS sequence, however single bit shift operation is performed), executed on the Xilinx Spartan family chip. Consequently, having calculated the value of correlation distance reflecting the maximum level of the autocorrelation (initial position of the autocorrelation peak-L), one has to determine corresponding values of the autocorrelation index for a set of corresponding distances: L+dL, L+2dL, L+3dL and L+4dL; where dL is the distance step of 10 m. Next, approaching the acquisition of the slope of the autocorrelation peak, data fitting algorithm is applied to the values represented by a sequence of points (A, B, C, D, E), in order to achieve the required Fig. 4 Calculation of the slope of the autocorrelation peak common axis. Assumption is made, that the resulting slope of the fitting line represents the overall slope of the autocorrelation peak (as Fig. 4 depicts, vertical axis indicates the correlation coefficient value, whereas the horizontal axis informs what correlation distance is considered onwards). Finally, the obtained values are compared to the reference atmospheric visibility factors, as to calibrate the device.

Reference atmospheric visibility calculation
The standard visibility V n is defined in relation to the distance R for which the threshold of the human discrimination is estimated for the contrast value of k min = 0.02 (Frank et al. 2004). According to Koschmieder (1924), the contrast formula (1) is as follows: Components of the equation state, the parameter b is dependent on k min ; solution for k min = 0.02 and R = V n results in parameter b equal to 3.912. The Eq.
Concerning assessment of a reference atmospheric visibility, utilized method for obtaining of the actual k[dB] parameter, is the analysis of visual data (i.e. still frame from a digital camera (Babari et al. 2011;Hautiere et al. 2008)). The irradiance of the captured picture is determined and its distribution along any of the preferred profiles (crosshair selection) is represented on the same picture (consider output presented on the Fig. 5). Analyzing further, Fig. 5 Contrast estimation for a specific foggy scenery the contrast level of the dark object, located over the known distance R is estimated (Holejko et al. 2013). Figure 5 depicts exemplary scenario, where the input parameters are: R = 1.2 km and k ≈ 0.18, thus the actual standard visibility is approximated on the level of V n ≈ 2.7 km.

Results and discussion
For the final clause, concerning results discussion, as plotted on the Fig. 6, relationship between reference atmospheric visibility and autocorrelation peak is depicted. The vertical axis represents the slope coefficient (i.e. the angular tilt) indicating atmospheric visibility assessed by autocorrelation analysis. As one can observe, the steepness of the characteristic resides upon the approximate range of V n = 100 m towards 1,200 m. Therefore, the plausible accuracy and measurement performance of the device might be reported on such defined distance interval. Next, in order to assess the link between the reference atmospheric visibility (V n ) and the relevant slope of the autocorrelation peak (ρ) one can position following equation: The numerical approach of obtaining coefficients in (4) is given by the optimal curve fitting algorithm of the dual-term exponential approximation (MathWorks Inc. MATLAB 2013). As presented on the Fig. 6, where positioning of data fitting curve and measurement outcomes is achieved, the theoretical range of the atmospheric visibility assessment of the designed sensor is limited towards the vicinity of the 3,500 m marker. However, the reliable interval of estimation is, as previously mentioned, spanned approximately on the range from 100 to 1,200 m.

Conclusion
Presented design is based on the detection of backscattered optical signals, that further utilizes correlation coding technique to assess the atmospheric visibility. By series of experimental measurements, as demonstrated in the article, evaluation of the performance of the designed device is stated. Investigation indicates the overall superiority of the correlation coding method applied in evaluation of discussed parameter (V n ). Hence, the possibility of assessing the atmospheric visibility over the range (as the numerical spatial parameter, not only as a single point figure) is achieved. What is more, pioneered solution is able to determine the V n parameter to the extent of 1,200 m, what is obtained with real-time processing speed and appropriate accuracy. The designed device is leading its kind, as it utilizes continuous wave, PRBS sequence signal generator, that if necessary might support extended range of measurements. Further, presented gauge might be characterized as a low complexity, cost effective, lightweight and compact device, suitable for continuous measurements and functional improvements of contemporarily deployed systems. Utilizing only a single light source and an avalanche photodiode receiver, the simplicity in design is achieved (measurement system does not imply additional optical switches, amplifiers nor photodiodes). Additionally, while estimation of the atmospheric visibility is performed, the described apparatus might be simultaneously utilized in the area of cloud base elevation measurement (ceilometer).
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.