Next Article in Journal
Species Discrimination among Three Kinds of Puffer Fish Using an Electronic Nose Combined with Olfactory Sensory Evaluation
Previous Article in Journal
Fiber-Optical Sensors: Basics and Applications in Multiphase Reactors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sunglint Detection for Unmanned and Automated Platforms

by
Shungudzemwoyo Pascal Garaba
1,*,
Jan Schulz
1,
Marcel Robert Wernand
2 and
Oliver Zielinski
1
1
Institute for Chemistry and Biology of the Marine Environment—Terramare, Carl von Ossietzky University of Oldenburg, Schleusenstraβe 1, Wilhelmshaven 26382, Germany
2
Royal Netherlands Institute for Sea Research, Physical Oceanography, Marine Optics & Remote Sensing, P.O. Box 59, Den Burg, Texel 1790AB, The Netherlands
*
Author to whom correspondence should be addressed.
Sensors 2012, 12(9), 12545-12561; https://doi.org/10.3390/s120912545
Submission received: 14 July 2012 / Revised: 20 August 2012 / Accepted: 29 August 2012 / Published: 13 September 2012
(This article belongs to the Section Remote Sensors)

Abstract

: We present an empirical quality control protocol for above-water radiometric sampling focussing on identifying sunglint situations. Using hyperspectral radiometers, measurements were taken on an automated and unmanned seaborne platform in northwest European shelf seas. In parallel, a camera system was used to capture sea surface and sky images of the investigated points. The quality control consists of meteorological flags, to mask dusk, dawn, precipitation and low light conditions, utilizing incoming solar irradiance (ES) spectra. Using 629 from a total of 3,121 spectral measurements that passed the test conditions of the meteorological flagging, a new sunglint flag was developed. To predict sunglint conspicuous in the simultaneously available sea surface images a sunglint image detection algorithm was developed and implemented. Applying this algorithm, two sets of data, one with (having too much or detectable white pixels or sunglint) and one without sunglint (having least visible/detectable white pixel or sunglint), were derived. To identify the most effective sunglint flagging criteria we evaluated the spectral characteristics of these two data sets using water leaving radiance (LW) and remote sensing reflectance (RRS). Spectral conditions satisfying ‘mean LW (700–950 nm) < 2 mW·m−2·nm−1·Sr−1’ or alternatively ‘minimum RRS (700–950 nm) < 0.010 Sr−1’, mask most measurements affected by sunglint, providing an efficient empirical flagging of sunglint in automated quality control.

1. Introduction

Automated and unmanned remote sensing from above-water, airborne and satellite platforms is a non-invasive approach in observing marine biochemical and geophysical characteristics on a regional or global scale [13]. Inevitably remote sensing measurements from these platforms are prone to meteorological conditions and sunglint contamination. Sunglint is a transient anomaly that occurs when sunlight is reflected from the seawater surface directly into the down looking optical sensor [4,5]. It is a product of Fresnel reflection from a number of ‘dancing facets’ on a wind disturbed seawater surface. Sunglint is influenced by the position of the sun, viewing angle of the optical sensor, water refractive index, cloud cover, wind direction, and speed [68].

Kay et al. [9] assessed sunglint correction models for optical measurements in marine environments identifying that most of the models rely partly on the black pixel assumption [10] and tend to some extent correct glint pixels. The black pixel assumption postulates that water leaving radiance is insignificant in the near infra-red spectrum. However, several reports have shown that in coastal and turbid waters this assumption is not valid. Hence correction models employing this assumption have a high probability of over- or underestimation of apparent and inherent optical properties [5,1012]. Thus, development of a further sunglint flag for masking obviously contaminated measurements is desirable, hence minimising the probability of errors likely to occur when using correction models.

Wernand [13] described a meteorological quality flagging method to optimise automated hyperspectral measurements for coastal and shelf seas. It masks incoming solar irradiance (ES) taken during low light, dusk, dawn and under rainfall. Additionally, to reduce measurements affected by sunglint he suggests the use of two optical sensors looking in different azimuthal directions to measure water surface leaving radiance. The lowest water surface leaving radiance for each measurement is then assumed to have the least sunglint. While this setup is useful and minimises sunglint effects on measurements, it requires additional sensors to be installed and there is a possible risk of sunglint affecting both sensors.

In this report we aim to derive sunglint quality control flags for automated and unmanned above-water hyperspectral measurements to complement meteorological flags reported by Wernand [13]. For validation we use sky and sea surface images from a simultaneously operated camera system. The automatic sunglint detection algorithm is implemented to provide an objective way of distinguishing sunglint from non-sunglint situations.

2. Data and Methods

Above-water hyperspectral radiometric measurements were conducted aboard the R/V Heincke cruise HE302, between 21 April and 14 May 2009 in the northwest European shelf seas (Figure 1). The campaign was within the scope of the North Sea Coast Harmful Algal Bloom (NORCOHAB II) field campaign.

2.1. Instrumentation

A RAMSES-ACC hyperspectral cosine irradiance meter (TriOS, Germany) was used to measure incoming solar radiation, ES (λ). Two RAMSES-ARC hyperspectral radiance meters (7° field-of-view in air), were used to detect the sea surface radiance Lsfcsfc, Φ, λ) and sky radiance Lskysky, Φ, λ). A frame (see Figure 2) designed to hold the irradiance sensor facing upwards, with the sky and sea surface radiance sensors at zenith angles θsfc = 45° and θsky = 135°, was fixed to the mast of the ship facing starboard, 12 m above sea surface. These spectral measurements were automatically collected at 15 min intervals over a spectral range λ = 320–950 nm in steps of 5 nm.

A DualDome D12 (Mobotix AG, Langmeil, Germany) camera system with field-of-view set to 45°, was used to capture sky and sea surface images simultaneous to hyperspectral measurements, as illustrated in Figure 2. Positioning height (∼12 m above sea surface) of camera and optical sensors proved to be unaffected by sea spray. The camera's field-of-view was set congruent with the area observed by the radiometers, Lsky and Lsfc. Ship's position and heading were recorded by a Differential Global Position System (DGPS) and sampling times were logged in Coordinated Universal Time (UTC).

2.2. Methods

The first quality control step involved implementing the meteorological flagging [13] on ES (λ) measurements using MATLAB 2010a (The MathWorks, GmbH, Ismaning, Germany). The three meteorological flag conditions are:

  • ES (λ = 480 nm) > 20 mW·m−2·nm−1 setting a threshold for which significant ES (λ) can be measured,

  • ES (λ = 470 nm)/ES (λ = 680 nm) < 1 masking spectra affected by dawn/dusk radiation,

  • ES (λ = 940 nm)/ES (λ = 370 nm) < 0.25 masking spectra affected by rainfall and high humidity.

ES (λ) that passed this meteorological flagging and corresponding Lsfcsfc, Φ, λ), Lskysky, Φ, λ) measurements were used to derive water leaving radiance, LWsfc, Φ, λ), and remote sensing reflectance, RRS (θ, Φ, λ) according to Equation (1) [14]:

R R S = L W E S = L sfc ( ρ air sea L sky ) E S
where ρair-sea is the air-sea interface reflectance factor assumed to be ρair-sea = 0.0256 correcting for sky reflected radiation. We tested the performance of sunglint correction models by Gould et al. [5] Mobley [8], Lee et al. [4] and Ruddick et al. [15] using its cloud cover and wind components of Ruddick et al. [15] (see Supplementary Material). Optical measurements and computed parameters can be freely accessed via the PANGAEA database of the World Data Center for Marine Environmental Sciences (WDC-Mare, http://www.pangaea.de): http://www.doi.pangaea.de/10.1594/PANGAEA.759690.

2.2.1. Automated Sunglint Image Detection Algorithm

A well exposed digital greyscale image of the sea surface normally shows several peaks scattered across the grey level histogram. Between peaks the histogram shows plateaus and has an averaged centred and balanced density distribution. When overexposed, images show a significant higher count of bright colours with white dominating. Consequently, the density distribution in the histogram shifts characteristically towards the white end where mostly one prominent peak is observed. The proposed automated sunglint detection algorithm is based on simple histogram arithmetic to evaluate the gray value ratio between dominating dark and bright counts (see Figure 3).

Processing of each full colour (24 bit colour depth) image starts with cropping image borders for 5% on all sides. As the camera is mounted stationary on the mast, field of view is prone to ship movements. Cropping has the advantage to reduce wash and white caps close to the vessel's hull, the hull itself or atmospheric interferences, depending on displacement from the vertical. Additionally it focuses the evaluated part of the image to the reading spot of the hyperspectral radiometer on the sea surface.

The cropped colour image is converted to greyscale and colour depth is reduced to b bits (with b < 24). In our case we use b = 8, resulting in a normal greyscale image with 256 levels. The histogram is computed for the greyscale image, normalised and multiplied by 2b.

Prior to determining the bright/dark ratio a lower (TL) and upper (TU) threshold is defined, with constraints 0 ≤ TL < TU < 2b. Threshold used here have been identified empirically and are in good congruence with results obtained by a human investigator. A value TL = 2b/2 is a suitable first, and TU = 2b × 0.75 second threshold (respectively values 128 and 192 for a 28 level greyscale image). In the histogram local maxima are identified within the intervals 0 ≤ × ≤ TL and TU ≤ × < 2b. Values between TL < × < TU are excluded.

Between the two maxima a line segment is created and its slope calculated. The perpendicular to the line segment is calculated, intersecting it at the half of the length of the line segment. A special case appears if the slope of the line segment is zero, as it results in a division by zero when calculating the slope of the perpendicular. In this case a vertical line is created instead. With the previously defined constraints the slope cannot become infinity.

The intersection of the finally determined straight line with the abscissa is used as indicator for overexposure. If this value is larger than TU the image is tagged as potentially overexposed.

2.2.2. Spectral Analysis for Sunglint Flag

The empirical sunglint flag, based on spectral information, was developed on the premise that open seawater is assumed to absorb all light in the NIR. Thus, any light signal measured by an optical sensor would be sea surface reflectance or atmospheric scattered radiance. However, in turbid waters multiple scattering influenced by optically active seawater constituents also contributes to this light signal measurable by an optical sensor [5,12,16]. Removing a fraction ρair-sea of the sky radiance Lsky has become standard procedure to minimise sunglint in water leaving radiance LW from, above-water remote sensing. However, this approach like any correction model does not fully guarantee sunglint free measurements [4,8].

Our objective is to provide a reliable sunglint flagging with the advantage of sea-surface images to validate an empirical spectral approach to eliminating measurements affected by sunglint. In Figure 4 a simplified activity diagram illustrates the steps that were implemented in this sunglint flag investigation and evaluation:

  • The automated sunglint image detection algorithm (Section 2.2.1) was applied to sea surface images that are matching the unmasked spectra validated with the meteorological flagging [13]. Based on the computed probability of sunglint contamination in the images, they were classified into Nns–image set without sunglint or Ns–sunglint-affected image set.

  • Spectra analysis of the two sets Nns and Ns was used to identify unique characteristics with respect to LW (λ) and RRS (λ). This included the investigation of each individual spectra and mean spectra for the two sets over the whole spectrum range (λ = 320–950 nm), to obtain a general overview on typical spectra for the sets Nns and Ns;

  • To obtain distinguishing features for spectra data sets in Ns and Nns utilising findings from step 2 (i.e., spectrum shape, behaviour and magnitude variations of both LW (λ) and RRS (λ) in the VIS and NIR) a combination of inequality equations and band ratios was implemented. These combinations include:

    • using the NIR mean spectra to obtain threshold values unique in set Nns and Ns. This procedure was repeated using the NIR minimum spectra obtained from the statistical analysis of the Nns and Ns sets respectively,

    • performing spectral band ratioing on characteristic spectral bands both in VIS and NIR, here λ = 400 nm, 460 nm, 760 nm, 940 nm [9,13,17].

  • To test the performance of the above mentioned procedures a reanalysis of the original united data set of all valid spectra from step 1 was run. Based on the comparison of this computation with the independent image algorithm (Section 2.2.1.) effective sunglint flagging criteria were identified.

3. Results and Discussion

A total of 629 from 3,121 spectral measurements passed the test conditions of the meteorological flagging. The image assessment of these unmasked 629 spectra came up with 501 images free of sunglint, Nns (having least visible/detectable white pixel or sunglint), and 128 sunglint-affected images, Ns (having too much or detectable white pixels or sunglint).

3.1. Sunglint Image Analysis

Automated and unmanned optical measurements from a seaborne platform are challenging. It is difficult to adhere to recommended sensor setups for θsfc, θsky, Φ; thus it is inevitable to collect measurements affected by sunglint, whitecap and foam [1820]. Figure 5 demonstrates typical situations also noted during the image inspection in step 2 [Section 2.2.1 and Figure 3(a,b)].

Using the automated sunglint image detection algorithm, the output was either an image with a probability of sunglint or without. Sunglint (Figure 5(A,B1,C)) as well as whitecaps and foam (Figure 5(B2,C)) were identified as main sources of contamination, with the latter two presumably resulting in similar sunlight influenced spectral patterns. This image algorithm detects overexposure which is a result of sunglint and/or whitecaps and foam; we however assume that sunglint plays a major role.

3.2. Sunglint Flag

A spectral assessment of the sample sets indicated that RRS (NIR) and LW (NIR) are significantly enhanced for Ns, compared to Nns. The enhanced spectra in set Ns were assumed to be a result of sunglint, whitecaps, and foam. In Figure 6 normalised mean spectral shapes for Nns (501 spectra) and Ns (128 spectra) illustrate these findings. Data normalisation was applied to simplify the visual comparison of spectra, dividing each LW (λ) and RRS (λ) measurement by the maximum value for each measurement. However, for determining the flagging criteria, the actual computed spectral measurements LW (λ) and RRS (λ) were used.

Figure 6 shows how the mean spectra, with respect to normLW and normRRS values, are enhanced in the presence of sunglint in both the VIS and NIR compared to the non-sunglint situation. To further investigate sunglint spectral characteristics, 13 of the 128 sea surface images in Ns were identified to be highly affected by sunglint (see Figure 7). Their spectral shapes also reveal the same trend of enhanced water leaving and reflectance signal over the measured spectrum.

Spectral band ratios, a conventional approach in remote sensing algorithms, were implemented to identify differences in Nns and Ns. In previous reports related to sunglint correction, the following characteristic spectral bands were used mainly due to their physical properties: (a) oxygen absorption band (λ ≈ 760 nm) which has been used in a prior sunglint correction model [21]; (b) water and precipitable water vapour absorption band (λ ≈ 940 nm) known to interact with solar radiation in the NIR [17,22]; and (c) coloured dissolved organic matter absorption from UV to visible, here investigated at (λ ≈ 400 nm) and (λ ≈ 460 nm), contributing to sunglint through multiple scattering [23,24].

Several tests were implemented to obtain the most effective sunglint flagging. Test results are summarized in Table 1. A ranking was introduced grouping the methods according to their performance in the detection of sunglint affected data.

The method of thresholds was chosen because it is a widely used technique in developing flagging and validation algorithms [13,25,26]. Another benefit of using a threshold is the possibility of adjusting them to improve sensitivity of algorithms or models [27]. The threshold values were obtained by repetitive testing of possible threshold values i.e., after statistically computing the min, max and mean for spectra in the NIR. The computed statistics provided a starting test point and altered to obtain the best threshold values aimed at; (i) masking/eliminating as many measurements in the sunglint-affected set Ns; and (ii) unmasking/keeping as many measurements in the sunglint free set Nns. The performance test summarized in Table 1, revealed the most effective sunglint flagging criteria if the effectiveness of both tasks is taken into account;

  • ‘mean (LW)NIR < 2 mW·m−2·nm−1·Sr−1’, being 97% effective in sets Nns and 91% effective in Ns (91%); Or alternatively

  • ‘minimum (RRS)NIR < 0.010 Sr−1’, with an effectiveness of 95% in sets Nns and 94% in Ns.

3.3. Remote Sensing Reflectance

In the previous section, a sunglint flagging was identified which can be implemented using either LW or RRS. Both, the meteorological flagging conditions [13] and herein proposed sunglint flag criteria are presented in Table 2. The first three flags rely on the incoming radiation, Es, thus masking measurements taken during dusk or with too low incoming solar radiation (Flag 1), during dawn (Flag 2), or under rainfall (Flag 3). Equation 1 is then applied to derive the water leaving radiance, LW and remote sensing reflectance, RRS, followed by the sunglint flag validation. In this study the sunglint flag was implemented using LW (Flag 4a) because it is the first product of Equation 1 but can be replaced by RRS (Flag 4b) as summarised in Table 2. In order not to limit flagging to one remote sensing product, e.g., LW or RRS only, Flag 4a or Flag 4b were proposed.

4. Conclusions and Outlook

In this study we developed two alternative sunglint flagging criteria as a quality control procedure for unmanned and automated above-water platforms measuring water reflectance. To provide an independent evaluation method for sunglint situations, an automated sunglint image detection algorithm was successfully implemented. The Flag 1 (the minimal ES accepted) is a bit arbitrary and certainly depends on the geographical position. Flag 2 can be discussed, although this flag does not depend on the region but more or less on sun rise and fall. We choose for a more or less average shaped spectrum (comparable with the form of the Es spectrum at high noon, only amplitude differs. Flag 4 is depends on the region as we are using a empirical approach and therefore makes the use of the camera important.

We assumed sunglint to be the main cause of error in the collected measurements after applying the meteorological flagging [13]. However, the influence of whitecaps and foam has been reported to cause both: (i) a decrease in reflectance in the NIR due to radiation absorption by large air bubbles [28,29], or physical coolness of residual foam [30]; (ii) enhanced reflectance occurring as soon as waves break generating thick strong reflecting foam [20]. The image analysis revealed that sunglint was also present in sea surface images influenced by whitecaps/foam, and it was not possible to distinguish how each of them contributes to the spectra. It was therefore assumed that for the available measurements sunglint and whitecaps/foam led to contaminated measurements. It is for this reason that in future studies the contribution of whitecaps and foam be specifically investigated with respect to sunglint flagging, possibly utilising novel glint measurement apparatus e.g., [31,32] integrated with sunglint and wave models such as [4,5,7,15,33].

The sunglint image detection algorithm introduced here has shown to be a valuable additional tool for detecting contaminated sensor readings. However, the used values for TL and TU are empirical. Results with the given values have shown a good congruence with visual inspection by a human investigator. It might useful to adjust thresholds to specific situations or locations. Higher counts in the bright part of the histogram result in the situation that the ‘target’ line intersects the x-axis at higher values. Thus, increased numbers of bright pixels simultaneously emphasize the shape of the histogram on the ‘bright side’. Using different values for TL and TU impacts the sensitivity of the function for darker/brighter pixels.With respects to sensor setup we suggest that prior to automated and unmanned optical measurements the Solar Position Algorithm [34] be utilised to predict or identify optimal relative angle of optical sensor setup to the sun's position despite the limiting factors of R/V pitch, roll and yaw motions [14,18,35]. Here the idea is not to necessarily change the cruise track but to perform underway optical measurements avoiding sunglint. Alternative approaches to minimise sunglint, but increasing technical requirements, would be (a) to turn the sensors by some automatic device e.g., RFlex ( http://www.sourceforge.net/p/rflex/wiki/Home/) or (b) to have more sensors and cameras looking at different azimuthal directions. However for fixed platforms, such as piles, these alternative approaches can be avoided if SPA is utilised beforehand to limit sunglint influenced measurements. Furthermore, we provide a set of step that can be followed to perform sunglint correction on above-water automated measuring platforms (Appendix A).

Recently operational oceanographic observatories are becoming more prominent and at the same time hyperspectral radiance sensor technology becomes increasingly affordable [2,3]. Therefore the application of reflectance measurements above the water surface, from stationary and moving platforms alike, is expected to significantly grow in numbers. Given this enormous amount of data, favourably processed in real-time, effective quality control procedures like the ones discussed here are more than supporting tools, they are a crucial prerequisite for trustworthy and manageable information.

Acknowledgments

The authors extend their gratitude to the master and crew of R/V Heincke, and to R. H. Henkel, B. Krock, B. Saworski and D. Voβ for their support during the field campaign. This work was supported by Institute of Marine Resources (IMARE) GmbH subsidised by European Regional Development Fund (ERDF), University of Applied Science Bremerhaven and Coastal Observation System for Northern and Arctic Seas (COSYNA).

Appendix A

The following recommendations assume automated and unmanned platform above-water radiometry is performed using the setup explained in this report. For a detailed explanation of the correction models refer to the respective protocols.

Prerequisite information:

-

Wind speed [m/s]

-

Cloud cover

Lsky(750)/Es(750) [15]

<0.05 means clear sky

≥0.05 means cloudy sky

-

Ship heading

-

Ship GPS data

-

Sea surface images

Step 1

Apply three meteorological flags:

  • ES (λ = 480 nm) > 20 mW·m−2·nm−1 setting a threshold for which significant ES (λ) can be measured.

  • ES (λ = 470 nm)/ES (λ = 680 nm) < 1 masking spectra affected by dawn/dusk radiation.

  • ES (λ = 940 nm)/ES (λ = 370 nm) < 0.25 masking spectra affected by rainfall and high humidity [13].

Step 2

Alternative 1

Usage of the collected sea surface images and implementation of the sunglint detection algorithm. It eliminates images too contaminated with sunglint. Respective matching spectra of images successfully passing the image detection algorithm are used in Step 3.

Alternative 2

Determine LW and RRS and apply the spectral conditions supplied for sunglint flagging. To validate and verify flagging performance cross check with the image algorithm.

If not satisfactory, follow the procedures described in Section 2.2.2. ‘Spectral Analysis for Sunglint Flag’ and adjust the sunglint flag conditions appropriately.

Alternative 3

Determine the relative azimuthal angle of the sensor to the sun using e.g., Solar Position Algorithm, SPA [34]. Using the ship's position the SPA will compute the sun's azimuthal and zenith angle at a given space-time spot.

Extract spectra collected in the optimal relative azimuthal angle of sensor to sun 90° ≤ Φ ≤ 135°.

In this case if required, verify measurements with the image sunglint detection algorithm to further eliminate images with too much glint or above average white pixels.

Proceed with next step.

Step 3—Sunglint Correction

To determine which coefficient optimally removes sea surface reflectance/sunglint a number of methods are available. It is important to test them and decide on which correction model to apply from these:

Lee et al. [4]

Lee et al., recommends a spectral optimization where ρair-sea = 0.022 = Fresnel Reflectance. Additionally a bias delta (Δ), obtained from comparing modeled RRS and in-situ RRS, is then used as residual sea surface reflectance. To calculate RRS use Equation (2):

R R S = L sfc E S ρ air sea L sky E S Δ

Ruddick et al. [15]

Ruddick et al., provides a model based on cloud cover and wind speed, W that is aimed at correcting sunglint for measurements from turbid to very turbid water, to calculate RRS use Equation (3).

-

Lsky(750)/Ed(750) test for presence of cloud cover in Lsky field of view [15]

<0.05 means clear sky ρair-sea = 0.0256 + 0.00039 W + 0.000034 W2

≥0.05, cloudy sky ρair-sea = 0.0256

R R S = L W E S = L sfc ( ρ air sea L sky ) E S

Gould et al. [5]

According to Gould et al. [5] correction can be performed for sunglint by assuming that sunglint consists of (a) spectrally variable sky radiance ρair-sea (0.021) × Lsky and (b) spectrally-flat sunglint and cloud reflected radiance, B. To calculate RRS use Equation (4):

R R S = L sfc E S ρ air sea L sky E S B

Mobley [8]

Mobley provides a number of simulations and recommends ρair-sea = 0.028 for clear skies and wind speed < 5 m/s. For higher wind speeds change adjust ρair-sea accordingly. For cloudy skies at all wind speeds he recommends ρair-sea = 0.028. However, take note that the higher the ρair-sea is you might over correct and get negative reflectance in the NIR. To calculate RRS use Equation (5):

R R S = L W E S = L sfc ( ρ air sea L sky ) E S

References

  1. International Ocean Colour Coordinating Group (IOCCG). Remote Sensing of Ocean Colour in Coastal, and Other Optically-Complex, Waters; Sathyendranath, S., Ed.; IOCCG: Dartmouth, NS, Canada, 2000. [Google Scholar]
  2. Moore, C.; Barnard, A.; Fietzek, P.; Lewis, M.R.; Sosik, H.M.; White, S.; Zielinski, O. Optical tools for ocean monitoring and research. Ocean Sci. 2009, 5, 661–684. [Google Scholar]
  3. Glenn, S.; Schofield, O.; Dickey, T.D.; Chant, R.; Kohut, T.; Barrier, H.; Bosch, J.; Bowers, L.; Creed, E.; Halderman, C.; et al. The expanding role of ocean color and optics in the changing field of operational oceanography. Oceanography 2004, 17, 86–95. [Google Scholar]
  4. Lee, Z.; Ahn, Y.-H.; Mobley, C.; Arnone, R. Removal of surface-reflected light for the measurement of remote-sensing reflectance from an above-surface platform. Opt Express 2010, 18, 26313–26324. [Google Scholar]
  5. Gould, R.W.; Arnone, R.A.; Sydor, M. Absorption, scattering, and, remote-sensing reflectance relationships in coastal waters: Testing a new inversion algorithm. J. Coast. Res. 2001, 17, 328–341. [Google Scholar]
  6. Zhang, H.; Wang, M. Evaluation of sun glint models using MODIS measurements. J. Quant. Spectrosc. Radiat. Transf. 2010, 111, 492–506. [Google Scholar]
  7. Cox, C.; Munk, W. Measurement of the roughness of the sea surface from photographs of the sun's glitter. J. Opt. Soc. Am. 1954, 44, 838–850. [Google Scholar]
  8. Mobley, C.D. Estimation of the remote-sensing reflectance from above-surface measurements. Appl. Opt. 1999, 38, 7442–7455. [Google Scholar]
  9. Kay, S.; Hedley, J.; Lavender, S. Sun glint correction of high and low spatial resolution images of aquatic scenes: A review of methods for visible and near-infrared wavelengths. Remote Sens. 2009, 1, 697–730. [Google Scholar]
  10. Siegel, D.A.; Wang, M.; Maritorena, S.; Robinson, W. Atmospheric correction of satellite ocean color imagery: The black pixel assumption. Appl. Opt. 2000, 39, 3582–3591. [Google Scholar]
  11. Shi, W.; Wang, M. An assessment of the black ocean pixel assumption for MODIS SWIR bands. Remote Sens. Environ. 2009, 113, 1587–1597. [Google Scholar]
  12. Jamet, C.; Loisel, H.; Kuchinke, C.P.; Ruddick, K.; Zibordi, G.; Feng, H. Comparison of three SeaWiFS atmospheric correction algorithms for turbid waters using AERONET-OC measurements. Remote Sens. Environ. 2011, 115, 1955–1965. [Google Scholar]
  13. Wernand, M.R. Guidelines for (Ship-Borne) Auto-Monitoring of Coastal and Ocean Colour. Proceedings of Ocean Optics XVI, Santa Fe, NM, USA, 18– 22 November 2002; Ackleson, S.G., Trees, C., Eds.; p. 13.
  14. Mueller, J.L.; Davis, C.; Arnone, R.; Frouin, R.; Carder, K.; Lee, Z.P.; Steward, R.G.; Hooker, S.; Mobley, C.D.; McLean, S. Above-Water Radiance and Remote Sensing Reflectance Measurement and Analysis Protocols; Goddard Space Flight Space Center: Greenbelt, MD, USA, 2003. [Google Scholar]
  15. Ruddick, K.G.; de Cauwer, V.; Park, Y.J.; Moore, G. Seaborne measurements of near infrared water-leaving reflectance: The similarity spectrum for turbid waters. Limnol. Oceanogr. 2006, 51, 1167–1179. [Google Scholar]
  16. Doxaran, D.; Babin, M.; Leymarie, E. Near-Infrared light scattering by particles in coastal waters. Opt Express 2007, 15, 12834–12849. [Google Scholar]
  17. Halthore, R.N.; Eck, T.F.; Holben, B.N.; Markham, B.L. Sun photometric measurements of atmospheric water vapor column abundance in the 940-nm band. J. Geophys. Res. 1997, 102, 4343–4352. [Google Scholar]
  18. Hooker, S.B.; Morel, A. Platform and environmental effects on above-water determinations of water-leaving radiances. J. Atmos. Ocean. Technol. 2003, 20, 187–205. [Google Scholar]
  19. Fougnie, B.; Frouin, R.; Lecomte, P.; Deschamps, P.-Y. Reduction of skylight reflection effects in the above-water measurement of diffuse marine reflectance. Appl. Opt. 1999, 38, 3844–3856. [Google Scholar]
  20. Moore, K.D.; Voss, K.J.; Gordon, H.R. Spectral reflectance of whitecaps: Instrumentation, calibration, and performance in coastal waters. J. Atmos. Ocean. Technol. 1998, 15, 496–509. [Google Scholar]
  21. Kutser, T.; Vahtmäe, E.; Praks, J. A sun glint correction method for hyperspectral imagery containing areas with non-negligible water leaving NIR signal. Remote Sens. Environ. 2009, 113, 2267–2274. [Google Scholar]
  22. Kleidman, R.G.; Kaufman, Y.J.; Gao, B.C.; Remer, L.A.; Brackett, V.G.; Ferrare, R.A.; Browell, E.V.; Ismail, S. Remote sensing of total precipitable water vapor in the near-IR over ocean glint. Geophys. Res. Lett. 2000, 27, 2657–2660. [Google Scholar]
  23. Bricaud, A.; Morel, A.; Prieur, L. Absorption by dissolved organic matter of the sea (yellow substance) in the UV and visible domains. Limnol. Oceanogr. 1981, 26, 43–53. [Google Scholar]
  24. Gallegos, C.L.; Correll, D.L.; Pierce, J.W. Modeling spectral diffuse attenuation, absorption, and scattering coefficients in a turbid estuary. Limnol. Oceanogr. 1990, 35, 1486–1502. [Google Scholar]
  25. Wang, W.; Jiang, W. Study on the seasonal variation of the suspended sediment distribution and transportation in the East China Seas based on SeaWiFS data. J. Ocean Univ China (English Edition) 2008, 7, 385–392. [Google Scholar]
  26. Lavender, S.J.; Pinkerton, M.H.; Moore, G.F.; Aiken, J.; Blondeau-Patissier, D. Modification to the atmospheric correction of SeaWiFS ocean colour images over turbid waters. Cont. Shelf Res. 2005, 25, 539–555. [Google Scholar]
  27. Freeman, E.A.; Moisen, G.G. A comparison of the performance of threshold criteria for binary classification in terms of predicted prevalence and kappa. Ecol. Model. 2008, 217, 48–58. [Google Scholar]
  28. Whitlock, C.H.; Bartlett, D.S.; Gurganus, E.A. Sea foam reflectance and influence on optimum wavelength for remote sensing of ocean aerosols. Geophys. Res. Lett. 1982, 9, 719–722. [Google Scholar]
  29. Frouin, R.; Schwindling, M.; Deschamps, P.-Y. Spectral reflectance of sea foam in the visible and near-infrared: In situ measurements and remote sensing implications. J. Geophys. Res. 1996, 101, 14361–14371. [Google Scholar]
  30. Marmorino, G.O.; Smith, G.B. Bright and dark ocean whitecaps observed in the infrared. Geophys. Res. Lett. 2005, 32, L11604. [Google Scholar]
  31. Ottaviani, M.; Stamnes, K.; Koskulics, J.; Eide, H.; Long, S.R.; Su, W.; Wiscombe, W. Light reflection from water waves: Suitable setup for a polarimetric investigation under controlled laboratory conditions. J. Atmos. Ocean. Technol. 2008, 25, 715–728. [Google Scholar]
  32. Ottaviani, M.; Merck, C.; Long, S.; Koskulics, J.; Stamnes, K.; Su, W.; Wiscombe, W. Time-resolved polarimetry over water waves: Relating glints and surface statistics. Appl. Opt. 2008, 47, 1638–1648. [Google Scholar]
  33. Munk, W. An inconvenient sea truth: Spread, steepness, and skewness of surface slopes. Annu. Rev. Marine Sci. 2009, 1, 377–415. [Google Scholar]
  34. Reda, I.; Andreas, A. Solar position algorithm for solar radiation applications. Sol Energy 2004, 76, 577–589. [Google Scholar]
  35. Hooker, S.B.; Zibordi, G. Platform perturbations in above-water radiometry. Appl. Opt. 2005, 44, 553–567. [Google Scholar]
Figure 1. R/V Heincke HE302 cruise track (red line) where above-water hyperspectral optical measurements were taken between 21 April and 14 May 2009.
Figure 1. R/V Heincke HE302 cruise track (red line) where above-water hyperspectral optical measurements were taken between 21 April and 14 May 2009.
Sensors 12 12545f1 1024
Figure 2. Optical sensor and camera system setup on the mast of R/V Heincke. Highlighted in red on the left side is the RAMSES hyperspectral radiometers setup, and on the right side is the DualDome D12 camera system with a sample set of captured sky and sea surface images, not to scale.
Figure 2. Optical sensor and camera system setup on the mast of R/V Heincke. Highlighted in red on the left side is the RAMSES hyperspectral radiometers setup, and on the right side is the DualDome D12 camera system with a sample set of captured sky and sea surface images, not to scale.
Sensors 12 12545f2 1024
Figure 3. Schematic diagram of the automated sunglint image detection algorithm. (A) Cropping original sea surface image; (B) Converting cropped image to greyscale; (C) Extract grey level histogram of cropped and converted image; (D) Define lower (TL) and upper (TU) threshold; (E) Find position and magnitude of local maxima of darker (xTL, yTL) and brighter (xTU, yTU) colours and calculate slope of line segment; (FG) Calculate midpoint (midx, midy) of line segment between (xTL, yTL) and (xTU, yTU); (H) Calculate slope and intercept of the target line being perpendicular to the previously calculated slope line crossing the midpoint (midx, midy); (I) Calculate target line's crossing with abscissa. If and only if this value is larger than TU, the image is tagged as probably overexposed.
Figure 3. Schematic diagram of the automated sunglint image detection algorithm. (A) Cropping original sea surface image; (B) Converting cropped image to greyscale; (C) Extract grey level histogram of cropped and converted image; (D) Define lower (TL) and upper (TU) threshold; (E) Find position and magnitude of local maxima of darker (xTL, yTL) and brighter (xTU, yTU) colours and calculate slope of line segment; (FG) Calculate midpoint (midx, midy) of line segment between (xTL, yTL) and (xTU, yTU); (H) Calculate slope and intercept of the target line being perpendicular to the previously calculated slope line crossing the midpoint (midx, midy); (I) Calculate target line's crossing with abscissa. If and only if this value is larger than TU, the image is tagged as probably overexposed.
Sensors 12 12545f3 1024
Figure 4. A simplified activity diagram showing the four steps involved in generating the new sunglint flag. (1) Automated sunglint image detection and sorting of spectra into Nns–image set without sunglint or Ns–sunglint-affected image set; (2) Analysis of the two spectra sets (λ = 320–950 nm); (3) Specific spectra analysis in visible λ = 320–700 nm and near infra-red λ = 700–950 nm range (4) Performance test and identification of effective sunglint flagging criteria.
Figure 4. A simplified activity diagram showing the four steps involved in generating the new sunglint flag. (1) Automated sunglint image detection and sorting of spectra into Nns–image set without sunglint or Ns–sunglint-affected image set; (2) Analysis of the two spectra sets (λ = 320–950 nm); (3) Specific spectra analysis in visible λ = 320–700 nm and near infra-red λ = 700–950 nm range (4) Performance test and identification of effective sunglint flagging criteria.
Sensors 12 12545f4 1024
Figure 5. Starboard side sky (top) and sea surface (bottom) images captured during R/V Heincke field campaign HE302 showing sources of erroneous measurements; (A)—sunglint; (B1)—sunglint and (B2)—whitecap or foam; and (C)—a combination of sunglint, whitecaps and foam.
Figure 5. Starboard side sky (top) and sea surface (bottom) images captured during R/V Heincke field campaign HE302 showing sources of erroneous measurements; (A)—sunglint; (B1)—sunglint and (B2)—whitecap or foam; and (C)—a combination of sunglint, whitecaps and foam.
Sensors 12 12545f5 1024
Figure 6. Normalised mean spectral shapes for the sunglint free set Nns (top) and sunglint set Ns (bottom). The green line highlights the spectral limits for the VIS (λ = 320–700 nm) and NIR (λ = 700–950 nm).
Figure 6. Normalised mean spectral shapes for the sunglint free set Nns (top) and sunglint set Ns (bottom). The green line highlights the spectral limits for the VIS (λ = 320–700 nm) and NIR (λ = 700–950 nm).
Sensors 12 12545f6 1024
Figure 7. Thirteen samples of spectral shapes for sea surface images strongly affected by sunglint from the sunglint image set Ns. On the left is the normalized water leaving radiance, normLW and on the right is the normalized remote sensing reflectance, normRRS. It is observed that sunglint enhances the measured signal especially in the NIR for both water leaving and remote sensing reflectance. The red spectrum shows the mean spectral shape and the green line highlights the VIS and NIR spectral ranges.
Figure 7. Thirteen samples of spectral shapes for sea surface images strongly affected by sunglint from the sunglint image set Ns. On the left is the normalized water leaving radiance, normLW and on the right is the normalized remote sensing reflectance, normRRS. It is observed that sunglint enhances the measured signal especially in the NIR for both water leaving and remote sensing reflectance. The red spectrum shows the mean spectral shape and the green line highlights the VIS and NIR spectral ranges.
Sensors 12 12545f7 1024
Table 1. A summarized evaluation of the most effective sunglint flagging ranked according to their performance in masking sunglint spectra in a reanalysis. To check for effectiveness, the percentage E (Ns) % and E (Nns) % was derived by dividing the number of spectra masked or unmasked by each condition with the actually number of spectra in sets (sunglint set Ns—128 spectra and non sunglint set Nns—501 spectra).
Table 1. A summarized evaluation of the most effective sunglint flagging ranked according to their performance in masking sunglint spectra in a reanalysis. To check for effectiveness, the percentage E (Ns) % and E (Nns) % was derived by dividing the number of spectra masked or unmasked by each condition with the actually number of spectra in sets (sunglint set Ns—128 spectra and non sunglint set Nns—501 spectra).
Sunglint Flag Test ConditionSunglint Affected Observations MaskedE (Ns) %Sunglint Free Observations UnmaskedE (Nns) %
Minimum (LW)NIR < 0.3 mW/m2·nm·Sr1239645491
Mean (RRS)NIR < 0.010 Sr−11239645390
Minimum (RRS)NIR < 0.010 Sr−11219547194
Minimum (LW)NIR < 0.4 mW/m2·nm·Sr1199347895
Mean (LW)NIR < 2.000 mW/m2·nm·Sr1179148797
Minimum (LW)NIR < 0.5 mW/m2·nm·Sr1159048697
Minimum (RRS)NIR < 0.012 Sr−11159047996
Mean (RRS)NIR < 0.015 Sr−11128847595
LW (940)/LW (400) < 0.1451118748096
Mean (LW)NIR < 2.500 mW/m2·nm·Sr1108649298
LW (765)/LW (460) < 0.2901108647795
Minimum(RRS)NIR < 0.015 Sr−11098549098
RRS (765)/RRS (400) < 0.7001098546893
RRS (760)/RRS (400) < 0.7001088447094
LW (940)/LW (460) < 0.0901088448897
RRS (760)/RRS (460) < 0.6501078447394
Table 2. Summarised meteorological and sunglint flag conditions. The meteorological criteria [13] are represented by Flags 1–3. The sunglint flag can be implemented either as Flag 4a or Flag 4b depending on availability of measured water leaving radiance, LW or remote sensing reflectance, RRS.
Table 2. Summarised meteorological and sunglint flag conditions. The meteorological criteria [13] are represented by Flags 1–3. The sunglint flag can be implemented either as Flag 4a or Flag 4b depending on availability of measured water leaving radiance, LW or remote sensing reflectance, RRS.
Flag NamePurposeTest Conditions
Flag 1The ‘minimal flag’ sets the lower limit for which significant incoming solar radiation can be measured.Es (480 nm) > 20 mW/m2 nm
Flag 2The ‘shape flag’ will mask optical measurements influenced by dusk ‘red colouring of the sky’ or dawn radiation.Es (470 nm)/Es (680 nm) < 1
Flag 3The ‘rainfall flag’ will mask optical measurements influenced by precipitation or high humidity.Es (940 nm)/Es (370 nm) > 0.25
Flag 4aThe ‘sunglint flag’ will mask optical measurements influence by sunglint based on LW.Mean Lw (700–950 nm) < 2 m W/m2 nm Sr
Flag 4bAlternative ‘sunglint flag’ based on RRS.Minimum RRS (700–950 nm) < 0.010 Sr−1

Share and Cite

MDPI and ACS Style

Garaba, S.P.; Schulz, J.; Wernand, M.R.; Zielinski, O. Sunglint Detection for Unmanned and Automated Platforms. Sensors 2012, 12, 12545-12561. https://doi.org/10.3390/s120912545

AMA Style

Garaba SP, Schulz J, Wernand MR, Zielinski O. Sunglint Detection for Unmanned and Automated Platforms. Sensors. 2012; 12(9):12545-12561. https://doi.org/10.3390/s120912545

Chicago/Turabian Style

Garaba, Shungudzemwoyo Pascal, Jan Schulz, Marcel Robert Wernand, and Oliver Zielinski. 2012. "Sunglint Detection for Unmanned and Automated Platforms" Sensors 12, no. 9: 12545-12561. https://doi.org/10.3390/s120912545

APA Style

Garaba, S. P., Schulz, J., Wernand, M. R., & Zielinski, O. (2012). Sunglint Detection for Unmanned and Automated Platforms. Sensors, 12(9), 12545-12561. https://doi.org/10.3390/s120912545

Article Metrics

Back to TopTop