Prerpints.org logo
Preprint
Article

An Evaluation of Sun-Glint Correction Methods for UAV-Derived Secchi Depth Estimations in Inland Water Bodies

Altmetrics

Downloads

120

Views

42

Comments

0

A peer-reviewed article of this preprint also exists.

supplementary.docx (88.97KB )

Submitted:

24 July 2023

Posted:

25 July 2023

You are already at the latest version

Alerts
Abstract
This study investigates the application of Unmanned Aerial Vehicles (UAVs) equipped with a Micasense RedEdge-MX multispectral camera for the estimation of Secchi Depth (SD) in inland water bodies. The research analyzed and compared five sun-glint correction methodologies — Hedley, Goodman, Lyzenga, Joyce, and threshold removed glint to model the SD values derived from UAV multispectral imagery, highlighting the role of reflectance accuracy and algorithmic precision in SD modeling. While Goodman's method showed a higher correlation (0.92) with in situ SD measurements, Hedley's method exhibited the smallest average deviation (0.65 m), suggesting its potential in water resource management, environmental monitoring, and ecological modeling. The study also underscored the Quasi-Analytical Algorithm (QAA) potential in estimating SD due to its flexibility to process data from various sensors without requiring in situ measurements, offering scalability for large-scale water quality surveys. The accuracy of SD measures, calculated using QAA, was related to variability in water constituents of coloured dissolved organic matter and the solar zenith angle. A practical workflow for SD acquisition using UAVs and multispectral data was proposed for monitoring inland water bodies.
Keywords: 
Subject: Environmental and Earth Sciences  -   Remote Sensing

1. Introduction

Secchi depth (SD), an essential measure of water transparency in aquatic ecosystems, provides a critical indication of water quality and ecological health [1,2,3]. In Europe, SD is important for following the water quality rules set by the European Water Framework Directive (WFD), which mandates member states to uphold ecological standards [4]. SD thus helps identify impairments affecting water quality and aids in devising and executing management plans to protect water bodies [5].
The Secchi disk is the main tool used for obtaining SD measurements, particularly in vast and complex aquatic ecosystems [6,7,8]. However, traditional methods – visual observations using a white Secchi disk, can be labour-intensive, time-consuming, and potentially influenced by observer bias [9,10]. They also typically offer point-level measurements, which may not reflect water transparency over larger areas [11,12].
The limitations of traditional methods have spurred interest in remote sensing techniques for estimating SD in complex aquatic ecosystems. Techniques ranging from satellite data to Unmanned Aerial Vehicles (UAVs) equipped with Light Detection and Ranging (LiDAR) systems and hyperspectral imagery have been explored for accuracy [13]. These techniques promise cost-effective and efficient high-resolution spatial data on water transparency over large areas, with multispectral imagery from UAVs emerging as a promising alternative [14].
SD is primarily influenced by three optical components: chlorophyll-a (Chl-a), coloured dissolved organic matter (CDOM), and total suspended matter (TSM) [15,16]. Various algorithms exist to calculate water parameters from multispectral data, including mechanistic models [17], artificial neural networks [18], and regression algorithms [19]. Each offers unique advantages and applicability, depending on the specific characteristics of the data and the aquatic ecosystem under study. For instance, a study by Chusnah and Chu [20] demonstrated the application of machine learning in estimating chlorophyll-a concentrations, which are commonly used as indicators for assessing the trophic level of lakes and the state of water quality. The study utilized machine learning to implement a band ratio algorithm and generate chlorophyll-a maps from Sentinel-2 and Sentinel-3 satellite images. However, mechanistic models, relying on physical laws and principles to simulate light attenuation and scattering processes in water, often provide the most reliable and accurate predictions, particularly when dealing with smaller datasets, that would not be sufficient enough for machine learning models [21].
Lee et al. [22] provided a foundational understanding of the optical properties influencing SD, which significantly contributed to the development of the Quasi-Analytical Algorithm (QAA) for more accurate and reliable water clarity estimations. The QAA, a commonly used mechanistic model for SD, has been utilized in various water bodies [21,23,24]. It has been applied to MODIS and MERIS satellite data, where it reduced the root mean square error (RMSE) of SD estimation from 1.5 m to 1.0 m. Furthermore, the QAA has been used to account for the residual error in reflectance data from MODIS satellite data, demonstrating its potential for remote sensing in monitoring and managing water resources [25]. The algorithm showed excellent results (R2 = 0.96, MAPD = 0.18) when validated with independent measurements covering oceanic, coastal, and lake waters [26]. However, previous studies have not delved into the potential benefits of utilizing multispectral cameras onboard UAVs. The integration of these cameras with UAVs offers potential improvements in spatial resolution and data availability, bridging the gap between in situ and satellite remote sensing measurements [27].
It is important to accurately account for reflected light from the water surface, more specifically the sun-glint effect, as it can lead to inaccuracies in further processing of water quality algorithms [27,28,29]. The most simple way of avoiding sun-glint is careful UAV flight time and direction planning, however, since the water surface is often uneven and it is hard to reduce sun-glint completely [30]. There are several methods of reducing the sun-glint effect in multispectral UAV images during post-processing, for example, M Muslim et al. [29] tested four methods proposed by Lyzenga et al. [31], Joyce [32], Hedley et al. [33], and Goodman et al. [34] and applied them to either the whole image or just the glinted area. Other studies used methods of Hochberg et al. [35], removing glinted pixels as NIR threshold, or using HydroLight simulations [27]. The application of these methods should align with the specific requirements of the study. For instance, if the primary objective entails bottom mapping, the method proposed by Lyzenga offers superior results, as demonstrated in M Muslim's et al. [29] study. On the other hand, if the analysis is focused on assessing surface water quality parameters, the Hedley method emerged as the preferred choice by Windle and Silsbe [27]. Thus, the selection of sun-glint correction techniques requires careful consideration of the objectives, ensuring the most effective and accurate outcomes in different scenarios.
This research aimed to evaluate the effectiveness of QAA for the determination of SD using multispectral cameras onboard UAVs. It was done by testing QAA on image datasets pre-processed using five different methods of sun-glint correction. Additionally, the study examined how water constituents and solar zenith angle affected the discrepancy between actual and modeled SD values. We hypothesise that the radiometric accuracy of calibrated UAS sensor should meet the required accuracy of 5 %, expected with ocean colour remote sensing when compared to in situ measures of Hooker et al. [36] and application of sun-glint correction methods will improve usability in SD modeling, by deviation of RMSD of the same 5 % from in situ measures, that would still provide practical results for water quality assessment.

2. Materials and Methods

The The study was performed in 43 water bodies in Lithuania (Figure 1), represented by high variability and proportion of optically active in-water components – turbidity, Chl-a and CDOM (Table S1), thus bolstering the robustness and generalizability of the findings. The research was conducted from May 2021 to May 2023.

2.1. UAV Data Collection and Calibration

Inspire 2 Unmanned Aerial Vehicle (UAV) was used to facilitate the acquisition of data, equipped with a RedEdge-MX camera [37]. The camera features five bands: blue (475 nm ± 16 nm), green (560 nm ± 13 nm), red (668 nm ± 8 nm), red edge (717 nm ± 6 nm), and near-infrared (NIR) (842 nm ± 28 nm). To ensure optimal image capture, the UAV was flown at a height of 60 m to optimize the balance between spatial resolution and the area covered in each image allowing for both relatively high-resolution imagery (from 3 to 4 cm/pixel) and a reasonable area coverage per flight (around 40 – 50 m2), with the camera programmed to capture images every three seconds. This frequency was necessitated by the inability to view the live feed and capture images in desired areas from the Micasense RedEdge-MX camera at a distance. The maximum distance between the SD measurement site, where GPS coordinates were recorded, and the image GPS coordinates were approximately 40 m, while the minimum distance was 0 when the measurement was above the boat. The operation of the UAV was coordinated by two individuals — one controlling the UAV from the shore and the other onboard a boat, both communicating via phone during the Secchi disk measurements. The onboard operator also concurrently measured other water quality parameters outlined in Section 2.2.
For image calibration, the original Micasense RedEdge calibration panel was used. According to Micasense's guidelines, a calibration panel image was captured at the shore before each flight, ensuring no shadows were cast on the reference panel or drone [37]. Three images were taken each time for assurance.

2.2. In Situ Data

Concurrent with the UAV image capturing, in situ SD measurements were performed using a 30-cm white Secchi disk, a widely accepted method for assessing water transparency in lakes, rivers, and oceans [38].
The parameters of Coloured Dissolved Organic Matter (CDOM) – the light-absorbing component of dissolved organic matter, Chlorophyll-a (Chl-a) – an indicator of phytoplankton biomass, and turbidity – a measure of the cloudiness or haziness of water were measured alongside SD to assess their potential influence on water transparency.
Water samples for Chl-a measurements were filtered through glass fiber GF/F filters with a nominal pore size of 0.7 μm and extracted into 90% acetone. Photosynthetic pigments were measured spectrophotometrically and estimated according to the trichromatic method [39,40]. CDOM was measured spectrophotometrically after filtration through 0.22 μm membrane filters. The CDOM absorption coefficient at 440 nm was derived according to Kirk [41]. SHIMADZU UV-2600 spectrophotometer was used for the analysis of Chl-a and CDOM. Turbidity was measured with a portable turbidity meter (Eutech InstrumentsTN-100, Landsmeer, The Netherlands) in the Nephelometric Turbidity Unit (NTU). The instrument has a light-emitting diode in the near-infrared range (Hach Lange at 860 nm and Eutech Instruments at 850 nm), and the detector measures the scatter at a 90° angle. This method is based on the International Standard Organisation (ISO) 7027.
The average SD of the lakes was 1.91 ± 1.53 m (± – standard deviation), ranging from 0.25 to 7.2 m. The mean CDOM was 2.95 ± 3.2 m−1 (min-max: 0.37 – 20.01 m−1), the mean Chl-a concentration – 26.59 ± 26.48 (1.13 – 113.23) and the mean turbidity – 6.76 ± 11.8 NTU (0.00 – 70.62 NTU).
In situ remote sensing reflectance Rrs was acquired to validate wavelengths from UAV observations of water surface reflectance. Rrs was measured in the spectral range of 400–800 nm by the simultaneous measurements of downwelling irradiance, upwelling radiance, and downwelling radiance, performed with a WISP-3 spectroradiometer [42]. Rrs was calculated according to equation 1:
R r s = L u − ρ L d E d
where Lu is the upwelling radiance, Ld – the downwelling radiance, Ed – the downwelling irradiance, ρ – a water surface reflectance factor equal to 0.028.
Central wavelengths of the Rededge MX camera were used (475 nm, 560 nm, 668 nm, 717 nm, 842 nm). NIR was excluded from comparison with in situ reflectances, which were not measured at 842 nm, however, NIR was still used for the sun-glint correction step section 1.4.

2.3. Pre-Processing the UAV Data

The collected UAV images underwent pre-processing of correcting atmospheric effects and standardizing the format suitable for analysis (Figure 2). The conversion of raw images to radiance involves several corrective steps to eliminate biases and errors, accounting for dark pixel offset, vignette effect, as well as aligning images due to distances between sensors, which could potentially affect the accuracy and reliability of the data. These steps were achieved using the Micasense Python (version 3.7, Python Software Foundation, 2018) workflow, mainly the function of raw_image_to_radiance, as described in the Micasense Github repository for users [37].
Since reflectance is relatively independent of illumination conditions [11,43], it is preferred for calculating SD using remote sensing images. The conversion of radiance to reflectance was performed using a reference panel of known reflectance to determine a scale factor between radiance and reflectance. This scale factor was applied to the entire image to obtain a reflectance image. The accuracy of the reflectance image was verified by extracting and checking the same reflectance panel region used to calculate the scale factor for any trends or inconsistencies [37].
Reflectance was further normalized by dividing it by π, assuming water as a Lambertian body, following the mathematical model (2).
Reflected radiance = ρ * (cos θ_i)/π
Here, cos θ_i represents the cosine of the angle between the incident light and the surface normal and ρ – a surface reflectance factor [43]. This normalization ensures that the reflectance values fall within a standardized range (0 to 1), facilitating consistent comparisons and calculations across different surfaces, lighting conditions, and measurement devices [44,45].

2.4. Image Masking and Sun-Glint Correction

From the initial set of water bodies, four were eliminated (Tūbausių Reservoir, Grabuostas Lake, Ilgis Lake, and Mušėjus Lake) from further analysis. One was excluded due to excessive cloud glint, while the other three measurements were performed too close to the shore or contained emerged macrophytes, leaving no sufficient area where the reflectance would not be affected (Figure 5a). One water body (Kruminių reservoir) was surveyed at a late time when the zenith angle was approximately 80 degrees, resulting in shadows covering a half portion of the area; however, the unshaded area could be used for the analyses (Figure 5b). In total, 39 water bodies were left after discarding unsuitable ones for statistical analysis.
Preprints 80381 g003
Figure 3. (a) RGB of an image affected by shadows, bottom reflectance and floating macrophytes after thresholding (Grabuostas Lake), (b) RGB of contaminated image (Kruminių reservoir) with shadow and shore but where values can still be used by selecting an area (red square) that was good, (c) RGB example of a relatively good image (Musia Lake) before removing sun-glint and boat.
Figure 3. (a) RGB of an image affected by shadows, bottom reflectance and floating macrophytes after thresholding (Grabuostas Lake), (b) RGB of contaminated image (Kruminių reservoir) with shadow and shore but where values can still be used by selecting an area (red square) that was good, (c) RGB example of a relatively good image (Musia Lake) before removing sun-glint and boat.
Preprints 80381 g003
To ensure accurate SD measurements, objects (boats, coast, macrophytes) and shadows potentially affecting the measurements' accuracy were removed. Firstly, the Normalized Difference Water Index (NDWI) was calculated to distinguish between water and non-water pixels (3). This index was then classified into binary parts using Yen thresholding [46,47], where non-water pixels were masked as NaN values. Yen method was used also on just the NIR band to classify areas where pixels were affected by sun-glint and masked as NaN in one of the sun-glint removal tests. This sun-glint removal method was later called the threshold-removed glint method.
N D W I = R r s G r e e n − R r s N i r R r s G r e e n + R r s N i r
In 11 images water waves caused a substantial amount of sun-glint (Figure 4a), despite images being captured early in the morning around 10 am GMT+3 when the sun reflection from calm water should not have reached the lens. Therefore, the sun-glint correction was performed using the algorithms developed by Hedley et al. [33], Goodman et al. [34], Lyzenga et al. [31], and Joyce [32], which assume that the water surface reflectance is a linear combination of water reflectance and sun-glint reflectance. The models were fitted using a set of training data consisting of image pixels where sun-glint was absent, calculated from an area with 10 % lowest value NIR pixels [27]. The model was then used to predict water reflectance for sun-glint-affected pixels, reducing reflectance values according to equations (4, 5, 6, 7).
The Hedley method (4) calculates the remote sensing reflectance (Rrs) for each pixel in each band. The method first finds the minimum NIR value and then calculates the slope for each band using linear regression. The Lyzenga and Joyce methods (5 and 6) are similar to the Hedley method but use respectively the mean and mode of the lowest 10% of NIR values instead of the minimum NIR value. The Goodman method (7) calculates the Rrs for each pixel in each band using a constant A (0.000019) and a factor B (0.1) that was multiplied by the difference between the red and NIR bands.
Hedley Rrsλ = λ + biλ ∗ (NIR - NIRmin)
Lyzenga Rrsλ = λ + bijλ ∗ (NIR - NIRmean)
Joyce Rrsλ = λ + biλ ∗ (NIR - NIRmode)
Goodman Rrsλ = λ – NIR + δ, where δ = A + B ∗ (Redλ - NIRλ)
where λ is the band of interest (blue, green, red, red edge, NIR), biλ is the slope of the band, bijλ is equal to covariance between λ and NIR divided by variance in NIR, and NIRmin NIRmode NIRmean is the minimum, mode or mean NIR value respectively. For Goodman’s method, δ is a constant offset across all wavelengths where A and B are constants (A = 0.000019 and B = 0.1). Additionally, removing only sun-glint affected areas (determined by the binary thresholding Yen algorithm) was compared as an alternative approach to sun-glint correction methods (Figure 4c).
The images of four water bodies (Vembutų, Stebuliškių, Pakapės, and Krūminių reservoirs) were also affected by either shore or cloud glint artefacts. However, most of these artefacts were successfully removed using binary thresholding, except for one image (Figure 3b) where the area had to be manually chosen because the thresholding method was unable to accurately separate the unwanted areas.

2.5. Secchi Depth Model

Quasi-analytical algorithm (QAA) proposed by Lee et al. [22], provides a robust framework for monitoring SD (8), particularly in scenarios where in situ measurements may be unavailable. Compared to the empirical approach, this semi-analytical method offers a significant advantage – it does not necessitate the recalibration of the retrieval algorithm with in situ data [23]. This enhances its utility in diverse settings. As such, it was considered more suitable for monitoring SD in various water bodies.
Lee et al. [22] paper introduced a mechanistic model that accounts for the effects of light attenuation, scattering, and reflection in the water column, as well as the properties of the Secchi disk itself, in determining the SD. The parameters that determine SD in this algorithm are the total absorption coefficient 'a and the total backscattering coefficient 'bb. From these parameters, the diffuse attenuation coefficient Kd was calculated, which is the main variable in the SD formula, besides the Rrs band.
K T _ K d = 1.04 ∗ 1 + 5.4 u 0.5 1 / 1 − s i n θ 2 R I 2 0.5
SD = 1 k t _ k d   ∗   m i n K d   l n 0.14 − m i n R r s 0.013
where minKd – minimum value Kd chosen out of Kd calculated with blue, green and red band, minRrs – above surface remote sensing reflectance of band that had the lowest Kd value. Kt_Kd – upwelling radiance diffuse attenuation coefficient and was used instead of a constant value of 1.5 as suggested by Jiang et al. [23], where θ is the solar zenith angle, RI is the refractive index value of water equal to 1.33 and u denotes the ratio of backscattering coefficient to the sum of absorption and backscattering coefficient.
The sun zenith angle was computed from in situ measurements using the time and location of the observation. This information was then used to calculate the solar position using the Python library pytz [48]. This approach provided a reliable means of determining the solar zenith angle, which is an important parameter in applications related to SD acquisition from Rrs.
Selecting the appropriate reference wavelength is of major importance for the final SD value [26]. The green band (560 nm) was used as a reference if the red band Rrs (668 nm) is < 0.0015 sr−1, otherwise the red band (668 nm) was used, with accordingly modified calculations of other parameters, as suggested by Lee et. Al. [22].
SD was calculated using the QAA method for images of 6 types: corrected after Hedley, Lyzenga, Joyce and Goodman sun-glint algorithms, threshold-removed images and images with no correction.

2.6. Validation and Interpret the Results

All of the final modeled SD values and reflectance images were averaged to ensure that images can be compared with in situ data. The final reflectances obtained after applying various correction methods and no correction (1.4 – 1.5 sections) were validated with in situ reflectances (WISP-3). The accuracy of these methods was evaluated by bias (9), while the root mean square deviation (RMSD) was used as an indicator of the QAA model‘s precision (10), and the Pearson’s correlation coefficient (r) described the relationship strength between the model‘s output and the real data values.
Bias = 1 N ∑ i = 1 N X e s t i m a t e d , i − X m e a s u r e d , i
RMSD = ∑ i = 1 N X e s t i m a t e d , i − X m e a s u r e d , i 2 N
The same accuracy and precision measures were applied to modeled SD values. The Generalized additive models (GAMs) were employed [49] to investigate the relationships between the difference in modeled and in situ SD (the response variable) and a set of independent variables: CDOM, Chl-a, turbidity, and Solar Zenith. The GAMs were chosen due to their flexibility in modeling non-linear relationships and their ability to handle interactions between predictors. Before analysis, the crosscorrelation (based on the Pearson correlation coefficient) between the independent variables was determined. The correlation was relatively high (r = 0.76) between Chl-a and turbidity, therefore turbidity was not included in the GAMs. Before the interpretation of the GAMs results, the residuals were visually inspected with diagnostic plots for normality and equal variance against the fitted values. F and p values were obtained to assess the relative importance and significance of the independent variables (p < 0.05 was considered as a statistically significant relationship). The fit of the model was evaluated using the explained deviance. A response plot was graphically represented to visualize the relationship between an explanatory variable and the response. As there were multiple predictors, each one was plotted separately with a smooth curve and the confidence interval of the effect.
To visualize the SD results, a subset of images with calculated SD values were exported as TIFF and mosaiced according to image metadata GPS coordinates in QGIS version 3.16 [50]. The result was then visualized on top of the RGB mosaic that was mosaiced using OpenDroneMap [51] photogrammetry software.

3. Results

3.1. Band Validation after Sun-Glint Correction

Across all correction methods and the in situ data, the green band (560 nm) consistently showed the peak mean value, while the blue band (475 nm) indicated the valley or lowest mean value. The general shape of the data appears to peak at the green band, with decreasing values on either side at the blue and red bands (668 nm), and a slight increase at the red edge band (717 nm). This pattern was consistent across all the correction methods and the in situ data, indicating the robustness of this spectral feature in the multispectral UAV image data.
The highest mean values across all bands were observed with the Lyzenga correction method, while the lowest was observed with the Goodman correction method (Figure 6). The in situ values were generally lower than the corrected values, with the exception of the Goodman correction method.
Preprints 80381 g005
Figure 5. Spectral scatterplots demonstrating the agreement between in situ measured and sun-glint corrected (or not corrected) reflectances in the UAV images. Statistics (r, RMSD, bias) and points are colour coded accordingly except for the red edge that is grey. A diagonal black line marks perfect agreement (1:1).
Figure 5. Spectral scatterplots demonstrating the agreement between in situ measured and sun-glint corrected (or not corrected) reflectances in the UAV images. Statistics (r, RMSD, bias) and points are colour coded accordingly except for the red edge that is grey. A diagonal black line marks perfect agreement (1:1).
Preprints 80381 g005
The validation results provided robust evidence regarding the performance of the correction methods and the reliability of the final reflectance values. The correlation for all methods (Figure 7) had the same trend where it was the lowest for the blue band (r = 0.38 ± 0.12) as well as for the red band (r = 0.56 ± 0.14) but higher for the green band (r = 0.75 ± 0.11) and red edge band (r = 0.80 ± 0.10). However, RMSD was relatively low for the blue band (RMSD = 0.0043 ± 0.0016) and red (RMSD = 0.0048 ± 0.0018) bands and slightly higher for the green (RMSD = 0.0060 ± 0.0013) and red edge (RMSD = 0.0053 ± 0.0014) bands it was relatively high. Bias for most methods followed a similar trend of larger underestimation for the blue band and green bands, then slight overestimation for the red band with Goodman, Hedley and threshold-removed glint methods but still undervaluation for other methods, and slightly lower undervaluation for the red edge band. Goodman‘s method was exceptional to these trends as green, red and red edge bands were overestimated.
Overall Goodman’s algorithm showed the highest correlation for the green and red edge bands, with r values of 0.90 and 0.92, respectively. However, it had low r values for the blue and red bands, with 0.32 and 0.71, respectively. The RMSD values for Goodman were relatively low across all bands, ranging from 0.002 to 0.005. Similarly, the Hedley algorithm consistently performed well across all bands, with r values ranging from 0.58 to 0.88. The RMSD values were similar to those of Goodman, ranging from 0.003 to 0.005.
The Joyce and Lyzenga algorithms showed matching performances, with lower r values and higher RMSD and bias values compared to Goodman and Hedley. The r values ranged from 0.27 to 0.69 for Joyce and from 0.26 to 0.67 for Lyzenga. The RMSD values ranged from 0.006 to 0.007 for both algorithms.
The control group with no correction applied showed moderate r values, ranging from 0.41 to 0.81. The RMSD values were similar to those of Goodman and Hedley, ranging from 0.005 to 0.006.

3.2. Validation of QAA SD Model

The performance of the QAA model‘s (Figure 8) ability to predict SD when compared with in situ measurements, showed a relatively high correlation across all methods with r ranging from 0.74 (threshold-removed glint) to 0.92 (Hedley glint correction). RMSD (from 0.65 to 1.05 m) and bias (from -0.78 to 0.58 m) showed acceptable results over all methods as well.
Comparing methods between themselves, the accuracy of all parameters had similar trends as accuracy for band comparison with in situ measurements (section 2.1), where the SD values were overestimated for the smaller SD values and underestimated for the larger ones, except for the Goodman’s method where the most of the SD values were overestimated.
Hedley’s sun-glint corrected images achieved the best results according to the RMSD measures (0.65 m), also r was just slightly smaller (0.91) than with the Goodman’s method (0.92) that achieved the best results based on the r value of 0.92. However, the RMSD value for Goodman’s method was relatively high (1.00 m) to the values of other methods. This method overestimated most of the values both for small and large SD values.
The Joyce and Lyzenga methods showed similar results between themselves with r values of 0.87 and 0.85, respectively. The RMSD values for Joyce and Lyzenga were 0.79 and 0.86 m, respectively and biases showed underestimation.
The worst performing method was when only sun-glint-affected pixels were removed and all other pixels were left unchanged (r = 0.75, RMDS = 1.05 m, bias = 0.13 m). The control group with no correction showed slightly better results: r value of 0.89, the RMSD value of 0.74 m and underestimation according to bias.

3.3. Relation with Water Constituents

The environmental factors in GAMs significantly explained (38.3 %) the variance of the difference in SD measurements. The interaction term of the sun zenith angle and CDOM was significant (F = 6.808, p-value < 0.05), suggesting that these factors together affect the accuracy of SD retrieval. The most important and statistically significant factor was CDOM (F = 10.47, p < 0.05), followed by the solar zenith (F = 4.84, p = 0.02). Chl-a values did not have a significant effect for the GAM model (F = 0.295, p = 0.59).
The modeled SD values were overestimated (to over 1.3 m) when CDOM values < 7.5 m−1 (Fig. 12), while underestimated (over 1 m) when CDOM was > 7.5 (although underestimation did not significantly change when CDOM >12). For the solar zenith, the SD values were overestimated (over 0.7 m) at the lower solar zenith angles (< 45 degrees), and then a relatively low effect (within ± 0.5 m) was between 45 and 75 degrees of the solar zenith angle and the underestimation sharply increased to over 1.5 m when the solar zenith angle was > 75 degrees.

4. Discussion

4.1. Advancements in SD Measurements

This research focused on evaluating the effectiveness of the QAA for determining SD using multispectral camera onboard UAVs. This novel approach expands the current understanding of measuring SD by introducing the potential use of UAVs, balancing the broad coverage provided by remote sensing methods with the high accuracy characteristic of in situ measurements, therefore serving as a more comprehensive alternative to traditional methods.
The experimental results did not uniformly support our hypothesis that the application of sun-glint correction methods would enhance the utility of multispectral images in SD modeling by reducing the RMSD by about 5 % and reducing the bias by a similar amount. The precision observed across all bands appears to correspond with the spectral band reflectance intensity, where bands demonstrating relatively lower reflectance intensity (e.g. blue band) typically displayed a lower compliance between the reflectance measured from the UAV and the in situ reflectance measurements and the reverse held true for bands with relatively higher reflectance intensity. Environmental factors like weak water surface signal and roughness of the water surface can introduce systematic and random errors, respectively, into water surface detection [13]. The correlation with the blue were identified as the least robust (Figure 7), which can be ascribed to their increased vulnerability to scattering, a trait inherent to these bands in the water environment [52].
Among the tested glint correction methods, the Hedley and Goodman methods emerged as the most effective across all bands, with Hedley demonstrating the lowest RMSD across all bands. For instance, a study conducted by M Muslim et al. [29] employed a similar methodology to our study, testing multiple sun-glint correction methods. Their findings indicated that the Lyzenga method yielded the most accurate results. However, it's important to note that the primary objective of their research was to map coral reefs, and in most of the study area, the bottom was visible, which may have influenced their results. In contrast, Windle and Silsbe [27] found that the Hedley method provided the highest correlation coefficient (r) and the lowest RMSD, corresponding closely with the results of our study.
The key distinction between the three glint correction methods (excluding the Goodman method) lies in the slope or covariance index of NIR band values. In the case of the Hedley method, the slope was the largest, as it takes the minimum value of the band as the starting point, as opposed to the mode (Joyce method) or the mean (Lyzenga method) [31,32,33]. This suggests that the images used in our study required a larger numeric correction due to the lowest initial reflectance values.
The results showed that the Goodman's method reflectance intensity values tend to be significantly different from other methods for the blue, red and red edge bands. This method, which was ranked as the second best performing in our study, adopted a different approach, using constant values to correct for the sun-glint effect [34]. This approach led to overall better r values between in situ and the green, red, and red edge bands on-board UAV. However, the performance of the blue band was significantly lower than when using the Hedley method. Despite this, the Goodman method demonstrated the lowest RMSD values and small overestimation bias, underscoring its effectiveness in applications where low reflectance intensity deviation from in situ radiometer is required.
The results also indicated that applying glint correction to the entire image rather than just the binary thresholded area can yield better results in terms of the overall accuracy of the reflectance values. This observation implies that the method of applying glint correction solely to the binary thresholded area might not be the most effective strategy, as it may overlook potential glint effects present beyond this area. Additional sun-glint correction research should be considered in the future for better generalization of reflectance correction [53].
The overall agreement between the in situ SD measurements and the modeled SD values might be connected to their handling of in-water constituents such as Coloured Dissolved Organic Matter (CDOM) and solar zenith angles, which significantly influenced the accuracy of the models. CDOM was found to predominantly influence the discrepancies between in situ measurements and modeled SD. Given that CDOM primarily absorbs light in the ultraviolet (UV) and blue regions of the spectrum, resulting in relatively low acquired signal by sensor [54], this agrees with our observation that the blue band demonstrated one of the weakest correlations between in situ and UAV-derived reflectance measurements [55]. It is plausible that Goodman's and Hedley’s methods more effectively managed this factor where Hedley’s had the highest r and Goodman’s had the lowest RMSD for blue band reflectance, resulting in a closer match with in situ values.
Similarly, high solar zenith angles above 70 o impacted the accuracy of modeled SD, leading to underestimations by up to 1.5 m. This is likely due to increased scattering and absorption of light at higher zenith angles, resulting in less light reaching the water's surface and thus larger differences between modeled and in situ SD values [56]. This reinforces the need for sophisticated algorithms that can accurately model these complex environmental factors in SD predictions.

4.2. Practical Applications

Accurate SD determination in large areas is particularly important given the increasing demand for high-resolution data on water transparency for applications including water resource management, environmental monitoring, and ecological modeling [57]. Measures of SD can reduce the need for boat measurements in lakes, also allowing for data collection at a higher frequency, surpassing traditional monthly monitoring intervals that may be inadequate for dynamic water bodies with recurring algal blooms, which in some cases can be inaccessible due to terrain or vegetation around the lake [58]. As well, UAVs could be used in shallow coastal waters (< 1.5 m), where a research vessel usually cannot access them, e.g. in the Curonian Lagoon [59].
The optimal fit model using Hedley sun-glint correction on the whole image was used to construct a mosaic from 45 individual images (Figure 9) of Kašučių Lake 2021 09 20. All of the pixels were left unmasked in order to show how shadows affect the final results, and therefore this aspect should be taken into consideration while planning the acquisitions.
This example shows the efficacy of the QAA best-fit model in transforming discrete data into coherent, large-scale representation of SD. The SD varied from 4 to 6 m, where the highest values were determined in the areas covered by tree shadows. Elsewhere SD values were more homogeneous (4.0 to 4.8 m), but some noise was still apparent in the orthomosaic (snow-like effect) as SD values are calculated for each pixel. For this reason, the mean value of the whole image was compared to in situ measurements in the area, instead of selecting a small square or point in the images, which would have resulted in lower accuracies.
Overall, UAVs can provide higher-resolution images than satellites, allowing for more precise measurements of water bodies, including small and shallow ponds [60,61], which are often not well captured by satellite imagery [61,62]. Furthermore, UAVs can be flown at a specific time and location, which is particularly important for SD measurements, as the measurement of water transparency can be influenced by various environmental factors, such as partial cloud cover, water turbidity, and the angle of the sun [27,63]. Although weather conditions influence the ability of UAV flight (with Inspire 2 UAV if winds exceed around 10 m/s), they can collect data during fully cloudy conditions, where optical satellite data are heavily affected by clouds. Additionally, they can be equipped with diverse sensors for capturing multi-spectral or hyperspectral images, thus providing valuable data on water quality and ecology [14,64]. Furthermore, UAVs can serve as ground truth to validate satellite data [65,66], thereby enhancing its accuracy and monitored area, e.g. one tile coverage of Sentinel-2 can reach from 100 000 ha, while UAVs cannot cover relatively large areas.
This method has potential applications in shaping policies and regulations related to water bodies. By yielding accurate, timely, and spatially inclusive data on water transparency, it equips decision-makers with invaluable resources, necessitating immediate intervention or remediation measures. Moreover, UAV technology presents a strategic advantage in regions where conventional methodologies struggle due to logistical complications. This is particularly evident when we consider the capabilities of UAVs for mapping vast areas. A single UAV flight, which takes approximately 25 min., can effectively map an area as large as 25 ha. To put this into context, Lithuania has around 6000 lakes, and about 340 of them are larger than 5 ha [67], which would be suitable for monitoring using Sentinel-2 satellite, assuming that the shape is not elongated. The rest of the smaller lakes could be monitored by applying methodology from this study and potentially could improve the accuracy and coverage of current monitoring, which covers around 80 lakes every 6 years in Lithuania. With the methods provided in this paper, processing workflow for one 5 ha lake should not take longer than 20 min. Moreover, this process can be fully automated, requiring only the supervision of the final results.

4.3. Future Research and Potential Limitations

The QAA model possesses the flexibility to process remote sensing data sourced from an array of sensors, including but not limited to MODIS, MERIS, OLCI, MSI and GOCI. By integrating a more diverse set of global in situ measurements corresponding to various water types, it's conceivable to further refine the precision and effectiveness of the QAA model.
While the research results demonstrate potential, it is essential to emphasize that the investigation was conducted across a broad range of unique aquatic environments and included measurements captured under variable sun zenith angles. Consequently, for future work, there's potential to modify the existing QAA used for SD estimations to better account for CDOM and sun zenith angle, that in this study have shown to be critical parameters. This modification could improve our understanding and predictions in the context of diverse and changeable aquatic conditions. There are several limitations that need to be addressed in future research. One of the main limitations is the assumption of linearity in the Hedley method. This method requires at least one dark pixel unaffected by glint [33] to be visible in the image, and if the whole image is affected by glint, this deglinting method will provide incorrect data. Understanding these limitations aids in contextualizing the results and fosters the development of more refined glint correction methodologies in the future. This is also relevant for a partial cloud glint, when clouds are reflected from the water surface in some, but not all areas of the image, for the analysis of multispectral drone imagery [33,68]. Most methodological studies for water surface mapping using UAVs suggest careful planning of flight time, preferably during conditions of clear skies and low sun-glint [30], however, this restricts one of the main advantages of UAVs – getting data on demand, especially if there is a need to visit several inland water bodies at the same day.
Many studies underscore the importance of choosing the appropriate quasi-analytical algorithm (QAA) for different optical water types [13,21,26]. The optical properties of coastal and inland waters, which are primarily influenced by the concentrations of suspended particulate matter, phytoplankton, and dissolved organic material, exhibit spatial and temporal variability, leading to diverse optical water types [69,70]. Applying a universal algorithm in these optically complex waters often results in significant uncertainties [69,70]. This insight highlights the importance of optical water classification to enhance retrieval accuracy, as indicated by numerous studies that developed class-specific algorithms for bio-optical parameters and achieved improvements by applying optical classification in the retrieval of these parameters [70,71,72,73]. In our study, relatively high concentrations of CDOM significantly affected the optical properties of the water and thereby introducing uncertainties in the results. One possible approach to improve the accuracy of the QAA is to calibrate the algorithm for the CDOM-dominated water bodies.
In addition to the aforementioned limitations, another significant challenge to consider is the interference of high vegetation pixels in the analysis of multispectral drone imagery [74]. Vegetation and their shadows, especially when in close proximity to water bodies, can skew the reflectance measurements and thus influence the accuracy of Secchi depth estimates. Consequently, implementing strategies to exclude these pixels during image analysis can enhance the reliability of the measurements.
The applicability of water parameters extends beyond SD, allowing for the calculation of additional parameters. Prior research has demonstrated the feasibility of employing multispectral UAVs for turbidity [75,76,77], Chl-a [75,77,78,79], CDOM [77], TSS [78,80], cyanobacteria [78] and macrophytes [46,61,81]. Given the inherent scalability of the UAV-based methodology, it stands as a promising tool for extensive SD assessments and water quality surveys, thereby facilitating large-scale studies focused on water transparency. This advancement brings us a step closer to exploiting the full potential of UAV-based remote sensing for assessing and monitoring aquatic environments.

5. Conclusions

This study revealed that the accuracy of SD measurements is profoundly influenced by sun-glint correction methods employed in UAV flights. There was a consistent agreement across all methods and the in situ radiometric data, particularly for the green band, emphasizing the robustness of multispectral UAV image data. Among the tested methods, Hedley's method demonstrated superior accuracy (RMSD = 0.65 m) and precision, thereby significantly contributing to the accuracy of the UAV-derived SD data.
Moreover, findings underscored the significant role of environmental factors, particularly the CDOM and solar zenith angle, causing inaccuracies in SD measurements: a solar zenith angle > 70 ° resulted in an underestimation of up to 1.5 m in modeled SD, while CDOM > 12 m−1 caused similar underestimations. Our research, therefore, supports the use of UAVs equipped with multispectral cameras as a viable method for SD determination in inland water bodies with SD of up to 7 m and lower than 12 m−1 CDOM. The results point towards an approach capable of reaching a correlation as high as 0.91 and reducing the RMSD by up to 12.85 % (Hedley’s method), thereby enhancing the versatility and reliability of SD measurements.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org. Table S1. Summary of investigated lakes and reservoirs in Lithuania. Descriptive statistics of in situ measured SD, chlorophyll-a, CDOM and turbidity.

Author Contributions

E.T.: Conceptualization, Methodology, Software, Validation, Formal Analysis, Investigation, Resources, Data Curation, Writing—Original Draft, Writing—Review and Editing, Visualization, Supervision; M.B.: Conceptualization, Methodology, Validation, Investigation, Resources, Data Curation, Writing—Review and Editing, Supervision; D.V.: Conceptualization, Validation, Investigation, Resources, Data Curation, Writing—Review and Editing, Supervision; J.G.: Conceptualization, Validation, Investigation, Resources, Writing—Review and Editing; I.B.: Resources, Data Curation, Funding Acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Doctorate scholarship program in Ecology and Environmental Sciences at Klaipeda University, Lithuania. The field campaigns were co-funded by Environmental Protection Agency contract “Assessment of the water condition and its more efficient management of the remote monitoring data collection, processing, use and storage system, to ensure accurate results (NUOTOLIS) [grant number 28T-2021-64/SUT-21P-20].

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cloern, J.E.; Abreu, P.C.; Carstensen, J.; Chauvaud, L.; Elmgren, R.; Grall, J.; Greening, H.; Johansson, J.O.R.; Kahru, M.; Sherwood, E.T.; et al. Human Activities and Climate Variability Drive Fast-Paced Change across the World’s Estuarine–Coastal Ecosystems. Global Change Biology 2016, 22, 513–529. [Google Scholar] [CrossRef] [PubMed]
  2. MacKay, M.D.; Neale, P.J.; Arp, C.D.; De Senerpont Domis, L.N.; Fang, X.; Gal, G.; Jöhnk, K.D.; Kirillin, G.; Lenters, J.D.; Litchman, E.; et al. Modeling Lakes and Reservoirs in the Climate System. Limnology and Oceanography 2009, 54, 2315–2329. [Google Scholar] [CrossRef]
  3. Paerl, H.W. Nuisance Phytoplankton Blooms in Coastal, Estuarine, and Inland Waters1. Limnology and Oceanography 1988, 33, 823–843. [Google Scholar] [CrossRef]
  4. European Parliament, Council of the European Union Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000 Establishing a Framework for Community Action in the Field of Water Policy; Luxembourg, 2000. 23 October 2000.
  5. Aas, E.; Høkedal, J.; Sørensen, K. Secchi Depth in the Oslofjord–Skagerrak Area: Theory, Experiments and Relationships to Other Quantities. Ocean Sci. 2014, 10, 177–199. [Google Scholar] [CrossRef]
  6. Ibelings, B.W.; Portielje, R.; Lammens, E.H.R.R.; Noordhuis, R.; Van Den Berg, M.S.; Joosse, W.; Meijer, M.L. Resilience of Alternative Stable States during the Recovery of Shallow Lakes from Eutrophication: Lake Veluwe as a Case Study. Ecosystems 2007, 10, 4–16. [Google Scholar] [CrossRef]
  7. Schultze-Lam, S.; Schultze-Lam, S.; Beveridge, T.J.; Des Marais, D.J. Whiting Events: Biogenic Origin Due to the Photosynthetic Activity of Cyanobacterial Picoplankton. Limnology and Oceanography 1997, 42, 133–141. [Google Scholar] [CrossRef] [PubMed]
  8. Spyrakos, E.; O’Donnell, R.; Hunter, P.D.; Miller, C.; Scott, M.; Simis, S.G.H.; Neil, C.; Barbosa, C.C.F.; Binding, C.E.; Bradt, S.; et al. Optical Types of Inland and Coastal Waters. Limnology and Oceanography 2018, 63, 846–870. [Google Scholar] [CrossRef]
  9. Jiang, H. Research on Transparency Derived from Remote Sensing and Its Spatial-Temporal Changes in Poyang Lake. China Rural Water and Hydropower 2012, 1, 30–33. [Google Scholar]
  10. Yu, D.F.; Xing, Q.G.; Lou, M.J.; Shi, P. Retrieval of Secchi Disk Depth in the Yellow Sea and East China Sea Using 8-Day MODIS Data. IOP Conf. Ser.: Earth Environ. Sci. 2014, 17, 012112. [Google Scholar] [CrossRef]
  11. Pham, T.N.; Ho, A.P.H.; Nguyen, T.V.; Nguyen, H.M.; Truong, N.H.; Huynh, N.D.; Nguyen, T.H.; Dung, L.T. Development of a Solar-Powered IoT-Based Instrument for Automatic Measurement of Water Clarity. Sensors 2020, 20, 2051. [Google Scholar] [CrossRef] [PubMed]
  12. Stock, A. Satellite Mapping of Baltic Sea Secchi Depth with Multiple Regression Models. International Journal of Applied Earth Observation and Geoinformation 2015, 40, 55–64. [Google Scholar] [CrossRef]
  13. Wang, D.; Xing, S.; He, Y.; Yu, J.; Xu, Q.; Li, P. Evaluation of a New Lightweight UAV-Borne Topo-Bathymetric LiDAR for Shallow Water Bathymetry and Object Detection. Sensors (Basel) 2022, 22, 1379. [Google Scholar] [CrossRef]
  14. Manfreda, S.; McCabe, M.F.; Miller, P.E.; Lucas, R.; Pajuelo Madrigal, V.; Mallinis, G.; Ben Dor, E.; Helman, D.; Estes, L.; Ciraolo, G.; et al. On the Use of Unmanned Aerial Systems for Environmental Monitoring. Remote Sensing 2018, 10, 641. [Google Scholar] [CrossRef]
  15. Castagna, A.; Amadei Martínez, L.; Bogorad, M.; Daveloose, I.; Dassevile, R.; Dierssen, H.M.; Beck, M.; Mortelmans, J.; Lavigne, H.; Dogliotti, A.; et al. Optical and Biogeochemical Properties of Belgian Inland and Coastal Waters; Biosphere – Biogeosciences. 2022. [Google Scholar]
  16. Garaba, S.P.; Friedrichs, A.; Voß, D.; Zielinski, O. Classifying Natural Waters with the Forel-Ule Colour Index System: Results, Applications, Correlations and Crowdsourcing. International journal of environmental research and public health 2015, 12, 16096–16109. [Google Scholar] [CrossRef] [PubMed]
  17. Giardino, C.; Bresciani, M.; Stroppiana, D.; Oggioni, A.; Morabito, G. Optical Remote Sensing of Lakes: An Overview on Lake Maggiore. J Limnol 2013, 73. [Google Scholar] [CrossRef]
  18. Anuar, N.N.M.; Fauzi, N.F.; Halim, H.; Khairudin, N.I.; Bakhtiar, N.S.A.; Shafii, N.H.B. Implementation of Long-Short Term Memory Neural Network (LSTM) for Predicting The Water Quality Parameters in Sungai Selangor. Journal of Computing Research and Innovation 2021. [Google Scholar] [CrossRef]
  19. Coskun, H.G.; Tanik, A.; Alganci, U.; Cigizoglu, H.K. Determination of Environmental Quality of a Drinking Water Reservoir by Remote Sensing, GIS and Regression Analysis. Water Air Soil Pollut 2008, 194, 275–285. [Google Scholar] [CrossRef]
  20. Chusnah, W.N.; Chu, H.-J. Estimating Chlorophyll-a Concentrations in Tropical Reservoirs from Band-Ratio Machine Learning Models. Remote Sensing Applications: Society and Environment 2022, 25, 100678. [Google Scholar] [CrossRef]
  21. Zhang, Y.; Zhang, Y.; Shi, K.; Zhou, Y.; Li, N. Remote Sensing Estimation of Water Clarity for Various Lakes in China. Water Research 2021, 192, 116844. [Google Scholar] [CrossRef]
  22. Lee, Z.; Shang, S.; Hu, C.; Du, K.; Weidemann, A.; Hou, W.; Lin, J.; Lin, G. Secchi Disk Depth: A New Theory and Mechanistic Model for Underwater Visibility. Remote Sensing of Environment 2015, 169, 139–149. [Google Scholar] [CrossRef]
  23. Jiang, D.; Matsushita, B.; Setiawan, F.; Vundo, A. An Improved Algorithm for Estimating the Secchi Disk Depth from Remote Sensing Data Based on the New Underwater Visibility Theory. ISPRS Journal of Photogrammetry and Remote Sensing 2019, 152. [Google Scholar] [CrossRef]
  24. Wei, L.; Wang, Z.; Huang, C.; Zhang, Y.; Wang, Z.; Xia, H.; Cao, L. Transparency Estimation of Narrow Rivers by UAV-Borne Hyperspectral Remote Sensing Imagery. IEEE Access 2020, 8, 168137–168153. [Google Scholar] [CrossRef]
  25. Chen, J.; Han, Q.; Chen, Y.; Li, Y. A Secchi Depth Algorithm Considering the Residual Error in Satellite Remote Sensing Reflectance Data. Remote Sensing 2019, 11, 1948. [Google Scholar] [CrossRef]
  26. Qing, S.; Cui, T.; Lai, Q.; Bao, Y.; Diao, R.; Yue, Y.; Hao, Y. Improving Remote Sensing Retrieval of Water Clarity in Complex Coastal and Inland Waters with Modified Absorption Estimation and Optical Water Classification Using Sentinel-2 MSI. International Journal of Applied Earth Observation and Geoinformation 2021, 102, 102377. [Google Scholar] [CrossRef]
  27. Windle, A.E.; Silsbe, G.M. Evaluation of Unoccupied Aircraft System (UAS) Remote Sensing Reflectance Retrievals for Water Quality Monitoring in Coastal Waters. Frontiers in Environmental Science 2021, 9. [Google Scholar] [CrossRef]
  28. De Keukelaere, L.; Moelans, R.; Knaeps, E.; Sterckx, S.; Reusen, I.; De Munck, D.; Simis, S.G.H.; Constantinescu, A.M.; Scrieciu, A.; Katsouras, G.; et al. Airborne Drones for Water Quality Mapping in Inland, Transitional and Coastal Waters—MapEO Water Data Processing and Validation. Remote Sensing 2023, 15, 1345. [Google Scholar] [CrossRef]
  29. M Muslim, A.; Wei Sheng, C.; Che Din, M.S.; Khalil, I.; Hossain, M. Coral Reef Mapping of UAV: A Comparison of Sun Glint Correction Methods. Remote Sensing 2019, 11, 2422. [Google Scholar] [CrossRef]
  30. Joyce, K.E.; Duce, S.; Leahy, S.M.; Leon, J.; Maier, S.W. Principles and Practice of Acquiring Drone-Based Image Data in Marine Environments. Mar. Freshwater Res. 2019, 70, 952. [Google Scholar] [CrossRef]
  31. Lyzenga, D.R.; Malinas, N.P.; Tanis, F.J. Multispectral Bathymetry Using a Simple Physically Based Algorithm. IEEE Trans. Geosci. Remote Sensing 2006, 44, 2251–2259. [Google Scholar] [CrossRef]
  32. Joyce, K.E. A Method for Mapping Live Coral Cover Using Remote Sensing. 2005.
  33. Hedley, J.D.; Harborne, A.R.; Mumby, P.J. Technical Note: Simple and Robust Removal of Sun Glint for Mapping Shallow-water Benthos. International Journal of Remote Sensing 2005, 26, 2107–2112. [Google Scholar] [CrossRef]
  34. Goodman, J.A.; Lee, Z.; Ustin, S.L. Influence of Atmospheric and Sea-Surface Corrections on Retrieval of Bottom Depth and Reflectance Using a Semi-Analytical Model: A Case Study in Kaneohe Bay, Hawaii. Appl. Opt. 2008, 47, F1. [Google Scholar] [CrossRef] [PubMed]
  35. Hochberg, E.J.; Andréfouët, S.; Tyler, M.R. Sea Surface Correction of High Spatial Resolution Ikonos Images to Improve Bottom Mapping in Near-Shore Environments. IEEE transactions on geoscience and remote sensing 2003, 41, 1724–1729. [Google Scholar] [CrossRef]
  36. Hooker, S.B.; Firestone, E.R.; Mcclain, C.R.; Esaias, W.E.; Barnes, W.; Guenther, B.; Endres, D.; Mitchell, B.G.; Barnes, R. SeaWiFS Calibration and Validation Plan, Volume 3; 1992. [Google Scholar]
  37. User Guide for MicaSense Sensors. Available online: https://support.micasense.com/hc/en-us/articles/360039671254-User-Guide-for-MicaSense-Sensors (accessed on 13 March 2023).
  38. Preisendorfer, R.W. Secchi Disk Science: Visual Optics of Natural Waters1. Limnology and Oceanography 1986, 31, 909–926. [Google Scholar] [CrossRef]
  39. Jeffrey, S.W.; Humphrey, G.F. New Spectrophotometric Equations for Determining Chlorophylls a, b, C1 and C2 in Higher Plants, Algae and Natural Phytoplankton. Biochemie und Physiologie der Pflanzen 1975, 167, 191–194. [Google Scholar] [CrossRef]
  40. Parsons, T.R.; Maita, Y.; Lalli, C.M. A Manual of Chemical and Biological Methods for Sea Water Analysis. Deep Sea Research Part A. Oceanographic Research Papers 1984, 31, 1523. [Google Scholar] [CrossRef]
  41. Kirk, J. Light and Photosynthesis in Aquatic Systems. Light and Photosynthesis in Aquatic Ecosystems, third edition 2010, VI, 1–651. [Google Scholar] [CrossRef]
  42. Hommersom, A.; Kratzer, S.; Laanen, M.; Ansko, I.; Ligi, M.; Bresciani, M.; Giardino, C.; Beltran-Abaunza, J.M.; Moore, G.; Wernand, M.R.; et al. Intercomparison in the Field between the New WISP-3 and Other Radiometers (TriOS Ramses, ASD FieldSpec, and TACCS). JARS 2012, 6, 063615. [Google Scholar] [CrossRef]
  43. Mobley, C. Estimation of the Remote-Sensing Reflectance from Above-Surface Measurements. Applied optics 2000, 38, 7442–7455. [Google Scholar] [CrossRef] [PubMed]
  44. Qiu, R.; Wang, S.; Shi, J.; Shen, W.; Zhang, W.; Zhang, F.; Li, J. Sentinel-2 MSI Observations of Water Clarity in Inland Waters across Hainan Island and Implications for SDG 6.3.2 Evaluation. Remote Sensing 2023, 15, 1600. [Google Scholar] [CrossRef]
  45. Wang, D.; Ma, R.; Xue, K.; Loiselle, S.A. The Assessment of Landsat-8 OLI Atmospheric Correction Algorithms for Inland Waters. Remote Sensing 2019, 11, 169. [Google Scholar] [CrossRef]
  46. Tiškus, E.; Vaičiūtė, D.; Bučas, M.; Gintauskas, J. Evaluation of Common Reed (Phragmites Australis) Bed Changes in the Context of Management Using Earth Observation and Automatic Threshold. European Journal of Remote Sensing 2023, 56, 2161070. [Google Scholar] [CrossRef]
  47. Yen, J.-C.; Chang, F.-J.; Chang, S. A New Criterion for Automatic Multilevel Thresholding. IEEE Transactions on Image Processing 1995, 4, 370–378. [Google Scholar] [CrossRef]
  48. Bishop Pytz Python Historical Timezone Library and Database. Available online: https://github.com/stub42/pytz (accessed on 3 April 2023).
  49. Wood, S.N. Fast Stable Restricted Maximum Likelihood and Marginal Likelihood Estimation of Semiparametric Generalized Linear Models. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 2011, 73, 3–36. [Google Scholar] [CrossRef]
  50. QGIS Version 3.16. Available online: https://qgis.org/en/site/ (accessed on 20 July 2023).
  51. Drone Mapping Software. Available online: https://www.opendronemap.org/ (accessed on 20 July 2023).
  52. Zaneveld, J.R.V.; Barnard, A.H.; Boss, E. Theoretical Derivation of the Depth Average of Remotely Sensed Optical Parameters. Opt. Express 2005, 13, 9052. [Google Scholar] [CrossRef] [PubMed]
  53. Overstreet, B.; Legleiter, C. Removing Sun Glint from Optical Remote Sensing Images of Shallow Rivers: Removing Sun Glint from River Imagery. Earth Surface Processes and Landforms 2016, 42. [Google Scholar] [CrossRef]
  54. Warren, M.; Simis, S.; Martinez-Vicente, V.; Poser, K.; Bresciani, M.; Alikas, K.; Spyrakos, E.; Giardino, C.; Ansper, A. Assessment of Atmospheric Correction Algorithms for the Sentinel-2A MultiSpectral Imager over Coastal and Inland Waters. Remote Sensing of Environment 2019, 225, 267–289. [Google Scholar] [CrossRef]
  55. Mamaghani, B.; Salvaggio, C. Multispectral Sensor Calibration and Characterization for SUAS Remote Sensing. Sensors 2019, 19, 4453. [Google Scholar] [CrossRef]
  56. Hashimoto, N.; Saito, Y.; Maki, M.; Homma, K. Simulation of Reflectance and Vegetation Indices for Unmanned Aerial Vehicle (UAV) Monitoring of Paddy Fields. Remote. Sens. 2019. [Google Scholar] [CrossRef]
  57. Soulignac, F.; Vinçon-Leite, B.; Lemaire, B.J.; Scarati Martins, J.R.; Bonhomme, C.; Dubois, P.; Mezemate, Y.; Tchiguirinskaia, I.; Schertzer, D.; Tassin, B. Performance Assessment of a 3D Hydrodynamic Model Using High Temporal Resolution Measurements in a Shallow Urban Lake. Environ Model Assess 2017, 22, 309–322. [Google Scholar] [CrossRef]
  58. Alikas, K.; Kratzer, S. Improved Retrieval of Secchi Depth for Optically-Complex Waters Using Remote Sensing Data. Ecological Indicators 2017, 77, 218–227. [Google Scholar] [CrossRef]
  59. Tyrimų laivas „Vėjūnas“ | Aplinkos apsaugos agentūra. Available online: https://www.gamta.lt/apie-agentura/tyrimu-laivas-vejunas/45 (accessed on 20 July 2023).
  60. Flynn, K.; Chapra, S. Remote Sensing of Submerged Aquatic Vegetation in a Shallow Non-Turbid River Using an Unmanned Aerial Vehicle. Remote Sensing 2014, 6, 12815–12836. [Google Scholar] [CrossRef]
  61. Müllerová, J.; Gago, X.; Bučas, M.; Company, J.; Estrany, J.; Fortesa, J.; Manfreda, S.; Michez, A.; Mokroš, M.; Paulus, G.; et al. Characterizing Vegetation Complexity with Unmanned Aerial Systems (UAS) – A Framework and Synthesis. Ecological Indicators 2021, 131, 108156. [Google Scholar] [CrossRef]
  62. Yang, S.; Wang, J.; Wang, P.; Gong, T.; Liu, H. Low Altitude Unmanned Aerial Vehicles (UAVs) and Satellite Remote Sensing Are Used to Calculated River Discharge Attenuation Coefficients of Ungauged Catchments in Arid Desert. Water 2019. [Google Scholar] [CrossRef]
  63. Yunlin, Z. Distribution, Seasonal Variation and Correlation Analysis of the Transparency in Taihu Lake. Transaction of Oceanology and Limnology 2003. [Google Scholar]
  64. Rahman, A.; Maulud, K.N.A.; Mohd, F.; Jaafar, O.; Tahar, K.N. Volumetric Calculation Using Low Cost Unmanned Aerial Vehicle (UAV) Approach. IOP Conference Series: Materials Science and Engineering 2017. [Google Scholar] [CrossRef]
  65. Otsu, K.; Pla, M.; Vayreda, J.; Brotóns, L. Calibrating the Severity of Forest Defoliation by Pine Processionary Moth with Landsat and UAV Imagery. Sensors (Basel, Switzerland) 2018. [Google Scholar] [CrossRef] [PubMed]
  66. Zhang, H.; Sun, Y.; Chang, L.-D.; Qin, Y.; Chen, J.; Qin, Y.; Du, J.; Yi, S.; Wang, Y. Estimation of Grassland Canopy Height and Aboveground Biomass at the Quadrat Scale Using Unmanned Aerial Vehicle. Remote. Sens. 2018. [Google Scholar] [CrossRef]
  67. Vaitekūnas, S.; Valančienė, E. Lietuvos Geografija (Geography of Lithuania); Alma litera. 2004. [Google Scholar]
  68. Partama, I.G.Y.; Kanno, A.; Ueda, M.; Akamatsu, Y.; Inui, R.; Sekine, M.; Yamamoto, K.; Imai, T.; Higuchi, T. Removal of Water-Surface Reflection Effects with a Temporal Minimum Filter for UAV-Based Shallow-Water Photogrammetry. Earth Surface Processes and Landforms 2018, 43, 2673–2682. [Google Scholar] [CrossRef]
  69. Le, C.; Li, Y.; Zha, Y.; Sun, D.; Huang, C.; Zhang, H. Remote Estimation of Chlorophyll a in Optically Complex Waters Based on Optical Classification. Remote Sensing of Environment 2011, 115, 725–737. [Google Scholar] [CrossRef]
  70. Cui, T.W.; Zhang, J.; Wang, K.; Wei, J.W.; Mu, B.; Ma, Y.; Zhu, J.H.; Liu, R.J.; Chen, X.Y. Remote Sensing of Chlorophyll a Concentration in Turbid Coastal Waters Based on a Global Optical Water Classification System. ISPRS Journal of Photogrammetry and Remote Sensing 2020, 163, 187–201. [Google Scholar] [CrossRef]
  71. Jackson, T.; Sathyendranath, S.; Mélin, F. An Improved Optical Classification Scheme for the Ocean Colour Essential Climate Variable and Its Applications. Remote Sensing of Environment 2017, 203, 152–161. [Google Scholar] [CrossRef]
  72. Yu, X.; Lee, Z.; Shen, F.; Wang, M.; Wei, J.; Jiang, L.; Shang, Z. An Empirical Algorithm to Seamlessly Retrieve the Concentration of Suspended Particulate Matter from Water Color across Ocean to Turbid River Mouths. Remote Sensing of Environment 2019, 235, 111491. [Google Scholar] [CrossRef]
  73. Zhang, F.; Li, J.; Shen, Q.; Zhang, B.; Tian, L.; Ye, H.; Wang, S.; Lu, Z. A Soft-Classification-Based Chlorophyll-a Estimation Method Using MERIS Data in the Highly Turbid and Eutrophic Taihu Lake. International Journal of Applied Earth Observation and Geoinformation 2019, 74, 138–149. [Google Scholar] [CrossRef]
  74. Doña, C.; Morant, D.; Picazo, A.; Rochera, C.; Sánchez, J.M.; Camacho, A. Estimation of Water Coverage in Permanent and Temporary Shallow Lakes and Wetlands by Combining Remote Sensing Techniques and Genetic Programming: Application to the Mediterranean Basin of the Iberian Peninsula. Remote Sensing 2021, 13, 652. [Google Scholar] [CrossRef]
  75. McEliece, R.; Hinz, S.; Guarini, J.-M.; Coston-Guarini, J. Evaluation of Nearshore and Offshore Water Quality Assessment Using UAV Multispectral Imagery. Remote Sensing 2020, 12, 2258. [Google Scholar] [CrossRef]
  76. Su, T.-C. A Study of a Matching Pixel by Pixel (MPP) Algorithm to Establish an Empirical Model of Water Quality Mapping, as Based on Unmanned Aerial Vehicle (UAV) Images. International Journal of Applied Earth Observation and Geoinformation 2017, 58, 213–224. [Google Scholar] [CrossRef]
  77. Zeng, C.; Richardson, M.; King, D.J. The Impacts of Environmental Variables on Water Reflectance Measured Using a Lightweight Unmanned Aerial Vehicle (UAV)-Based Spectrometer System. ISPRS Journal of Photogrammetry and Remote Sensing 2017, 130, 217–230. [Google Scholar] [CrossRef]
  78. Becker, R.H.; Sayers, M.; Dehm, D.; Shuchman, R.; Quintero, K.; Bosse, K.; Sawtell, R. Unmanned Aerial System Based Spectroradiometer for Monitoring Harmful Algal Blooms: A New Paradigm in Water Quality Monitoring. Journal of Great Lakes Research 2019, 45, 444–453. [Google Scholar] [CrossRef]
  79. Choo, Y.; Kang, G.; Kim, D.; Lee, S. A Study on the Evaluation of Water-Bloom Using Image Processing. Environ Sci Pollut Res 2018, 25, 36775–36780. [Google Scholar] [CrossRef] [PubMed]
  80. Olivetti, D.; Roig, H.; Martinez, J.-M.; Borges, H.; Ferreira, A.; Casari, R.; Salles, L.; Malta, E. Low-Cost Unmanned Aerial Multispectral Imagery for Siltation Monitoring in Reservoirs. Remote Sensing 2020, 12, 1855. [Google Scholar] [CrossRef]
  81. Doughty, C.L.; Ambrose, R.; Okin, G.; Cavanaugh, K. Characterizing Spatial Variability in Coastal Wetland Biomass across Multiple Scales Using UAV and Satellite Imagery. Remote Sensing in Ecology and Conservation 2021. [Google Scholar] [CrossRef]
Preprints 80381 g001
Figure 1. Study area with indicated study sites. n = 43
Figure 1. Study area with indicated study sites. n = 43
Preprints 80381 g001
Preprints 80381 g002
Figure 2. Processing workflow for selected UAV images. Arrows represent image processing and the dashed line arrows represent metadata.
Figure 2. Processing workflow for selected UAV images. Arrows represent image processing and the dashed line arrows represent metadata.
Preprints 80381 g002
Preprints 80381 g004
Figure 4. (a) RGB image before the sun-glint correction (b). Deglinted RGB image using Hedley’s method. (c) Threshold removed sun-glint area from NDWI image, recalculated to SD. (d) The boat removed from NDWI using Yen’s threshold and recalculated to SD.
Figure 4. (a) RGB image before the sun-glint correction (b). Deglinted RGB image using Hedley’s method. (c) Threshold removed sun-glint area from NDWI image, recalculated to SD. (d) The boat removed from NDWI using Yen’s threshold and recalculated to SD.
Preprints 80381 g004
Preprints 80381 g006
Figure 6. Mean r, RMSD and bias between corrected reflectances with each method and in situ measured reflectances, across all available wavelengths.
Figure 6. Mean r, RMSD and bias between corrected reflectances with each method and in situ measured reflectances, across all available wavelengths.
Preprints 80381 g006
Preprints 80381 g007
Figure 7. Scatterplots with linear regression line and 95 % degrees of freedom (green) of modeled and in situ SD. Best-fit line (1:1) (black) for each image dataset, using different pre-processing methods. (a) no correction, (b) threshold removed glint, (c) Goodman, (d) Hedley, (e) Joyce, (f) Lyzenga (n = 39).
Figure 7. Scatterplots with linear regression line and 95 % degrees of freedom (green) of modeled and in situ SD. Best-fit line (1:1) (black) for each image dataset, using different pre-processing methods. (a) no correction, (b) threshold removed glint, (c) Goodman, (d) Hedley, (e) Joyce, (f) Lyzenga (n = 39).
Preprints 80381 g007
Preprints 80381 g008
Figure 8. GAMs fitted smooth lines (blue lines) between the difference in SD (modeled vs. in situ) and graphs for each independent environmental factor: (a) Solar zenith (b) CDOM (c) Chl-a and (d) interaction plot of solar zenith angle, CDOM and predicted difference plane between in situ and modeled SD (m). The red dashed lines show a 95 % confidence interval for fitted lines. .
Figure 8. GAMs fitted smooth lines (blue lines) between the difference in SD (modeled vs. in situ) and graphs for each independent environmental factor: (a) Solar zenith (b) CDOM (c) Chl-a and (d) interaction plot of solar zenith angle, CDOM and predicted difference plane between in situ and modeled SD (m). The red dashed lines show a 95 % confidence interval for fitted lines. .
Preprints 80381 g008
Preprints 80381 g009
Figure 9. Figure 9. An SD transect of 45 mosaiced images in Kašučių Lake 2021 09 20. The dot with colour represents in situ SD of 4.3 m.
Figure 9. Figure 9. An SD transect of 45 mosaiced images in Kašučių Lake 2021 09 20. The dot with colour represents in situ SD of 4.3 m.
Preprints 80381 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated