Next Article in Journal
Recognition of Activities of Daily Living Based on Environmental Analyses Using Audio Fingerprinting Techniques: A Systematic Review
Next Article in Special Issue
Hyperspectral Image Classification for Land Cover Based on an Improved Interval Type-II Fuzzy C-Means Approach
Previous Article in Journal
Analysis of the Sensing Properties of a Highly Stable and Reproducible Ozone Gas Sensor Based on Amorphous In-Ga-Zn-O Thin Film
Previous Article in Special Issue
Comparison of Random Forest, k-Nearest Neighbor, and Support Vector Machine Classifiers for Land Cover Classification Using Sentinel-2 Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Proposal of a Method to Determine the Correlation between Total Suspended Solids and Dissolved Organic Matter in Water Bodies from Spectral Imaging and Artificial Neural Networks

by
Maurício R. Veronez
1,2,3,*,
Lucas S. Kupssinskü
3,
Tainá T. Guimarães
4,
Emilie C. Koste
1,
Juarez M. Da Silva
3,
Laís V. De Souza
2,
William F. M. Oliverio
3,
Rogélio S. Jardim
3,
Ismael É. Koch
3,
Jonas G. De Souza
3,
Luiz Gonzaga, Jr.
1,3,
Frederico F. Mauad
4,
Leonardo C. Inocencio
1 and
Fabiane Bordin
1
1
Advanced Visualization & Geoinformatics Lab—VizLab, Unisinos University, São Leopoldo 93022-750, Brazil
2
Graduate Programme in Geology, Unisinos University, São Leopoldo 93022-750, Brazil
3
Graduate Programme in Applied Computing, Unisinos University, São Leopoldo 93022-750, Brazil
4
Graduate Programme in Environmental Engineering Sciences, São Carlos Engineering School, University of São Paulo, São Carlos 13566-590, Brazil
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(1), 159; https://doi.org/10.3390/s18010159
Submission received: 16 November 2017 / Revised: 3 January 2018 / Accepted: 5 January 2018 / Published: 9 January 2018
(This article belongs to the Special Issue Analysis of Multispectral and Hyperspectral Data)

Abstract

:
Water quality monitoring through remote sensing with UAVs is best conducted using multispectral sensors; however, these sensors are expensive. We aimed to predict multispectral bands from a low-cost sensor (R, G, B bands) using artificial neural networks (ANN). We studied a lake located on the campus of Unisinos University, Brazil, using a low-cost sensor mounted on a UAV. Simultaneously, we collected water samples during the UAV flight to determine total suspended solids (TSS) and dissolved organic matter (DOM). We correlated the three bands predicted with TSS and DOM. The results show that the ANN validation process predicted the three bands of the multispectral sensor using the three bands of the low-cost sensor with a low average error of 19%. The correlations with TSS and DOM resulted in R2 values of greater than 0.60, consistent with literature values.

1. Introduction

Water quality monitoring involves data and water sample collection in the field and subsequent laboratory analysis. The effectiveness of the monitoring efforts depends on several factors, such as the frequency of sampling and the spatial distribution of parameters considered in the analysis [1,2,3]. Conventional water collection and analysis techniques are often costly and time-consuming, and may render water quality monitoring projects unfeasible. This is exacerbated by the locations of many water bodies, which may be in places that are difficult or dangerous to access, making on-site monitoring impossible. Therefore, there is a need to develop reliable and cost-effective spatial techniques for monitoring water quality that can be easily deployed [4].
The assessment of water quality based on remote sensing techniques is not new [5,6]. According to [7], several types of satellite-borne sensors have been used in the last three decades for monitoring aquatic environments. However most of these studies are based on assessing the quality of marine and coastal waters [5,6,8,9].
Challenges to advancing sustainable water quality management practices include the requirement for rapid and accurate monitoring and assessment of water quality in rivers, reservoirs, lakes, and estuaries, as well as knowledge of their spatiotemporal distribution [10,11,12]. Using appropriate technologies to effectively monitor these natural resources is, therefore, necessary. Remote sensing could be a complementary or alternative tool for monitoring water resources.
Remote sensing methods measure the changes in the spectral behavior of optically active components (OACs) in a water body, based on the variations in their concentration. OACs can interact with electromagnetic radiation (EMR), influencing the spectral behavior of water with their own spectral signatures, and altering its transparency. Such changes can be identified by sensor systems [12]. These substances can be identified through remote sensing techniques because their presence in a water body results in different EMR absorption and backscattering patterns, which are characteristic of each component. Suspended sediments and dissolved organic material are among the main parameters that influence the absorption and scattering of EMR in a water body [13].
Although satellite remote sensing technologies are effective in estimating the quality of water in large bodies, they still have the following limitations: (a) atmospheric interference in the images, affecting data interpretation and leading to a loss of confidence in the estimated water quality parameters; (b) low spatial resolution, making the detection of water pollution in small areas difficult; and (c) inefficient image review cycles that do not meet the need for rapid monitoring of these environments [11].
Sensors mounted on unmanned aerial vehicles (UAVs) have received increasing attention because they can provide high spatial resolution images, allow for the constant monitoring of the environment [14], and access the places that would otherwise be difficult or dangerous to visit [15]. Spectral sensors can be mounted on UAVs to provide images in different channels of the spectrum, which can be processed by various methods and algorithms, generating point clouds and surface models for various applications [16].
UAVs have been used to monitor the concentration of chlorophyll-a [17,18,19] and aquatic vegetation [20,21], spatialize suspended sediments [1], conduct environmental inspections [22], study flood events [23], conduct morphological and hydraulic characterizations [24], and monitor water pollution [25]. However, few studies have focused on determining the spectral behavior of water bodies using UAVs to estimate the concentration of multiple OACs. This is primarily due to the high costs of multi/hyperspectral sensors, which are often valued at USD 10,000 and USD 40,000.
Low-cost compact sensors are currently available on the market (~USD 600), with a modified filter sensitive to near-infrared radiation. This filter does not allow the passage of red-wavelength radiation (approximately between 600 and 700 nm), but it captures all other frequencies, including infrared, generating a charge-coupled device (CCD) with red channels that respond to the infrared band. These sensors have commonly been used in precision agriculture.
Thus, (a) to monitor water quality, we need multi/hyperspectral sensors that are sensitive to the electromagnetic spectrum between the visible and near infrared (400 to 800 nm) regions; (b) multi/hyperspectral sensors that can be used in UAVs are very expensive; and (c) low-cost sensors are not sensitive to the bands of the electromagnetic spectrum that are important for water quality monitoring.
These issues motivated us to propose an alternative method, based on an expert system for monitoring water quality using UAVs. We established the following hypothesis: it is possible to predict bands in the sensitive spectrum to monitor water quality using a low-cost sensor sensitive to the blue (B), green (G), and red (R) spectral bands. The proposed expert system is an artificial neural network (ANN).

2. Materials and Methods

The method that we are proposing can be structured according to the following steps: overflight with the UAV and processing of the acquired images; LANDSAT 8 (OLI) image processing; collection of water samples and laboratory analysis; training and validation of the ANN; generation of NDVI and NDWI images by the ANN; and establishment of mathematical models to correlate the NDVI and NDWI indices with TSS and DOM. The flowchart of the proposed method is depicted in Figure 1 and detailed in the following items.

2.1. Field Site

The lake at Unisinos University, South Brazil, was selected as the study site (Figure 2). The lake is a small artificial water body and receives rainwater runoff of the campus, as it is located at the lowest altitude of the university’s terrain. Due to this, several inorganic and organic compounds reach this aquatic environment through the rainwater run-off into the drainage area. These compounds can be found in the lake in the form of suspended solids or dissolved organic matter. It is approximately 320 m and 170 m at its longest and widest points, respectively, with a surface area of approximately 25,000 m2 and a depth of 4 m at the center. This lake was chosen because the equipment and resources available at the university could be easily accessed for this study.

2.2. Data Acquisition

Field sampling was performed in March 2017 and consisted of a UAV flight over the lake and an in situ collection of water samples at 21 sampling points, all occurring on the same day. Figure 3 shows the positions of the 21 sampling points.
The samplings were carried out in the superficial layer of water. The laboratory analysis for determination of total suspended solids (TSS) and dissolved organic matter (DOM) were made by the gravimetric method using the solids series (suspended and dissolved, fixed and volatile). DOM was obtained through the dissolved and volatile fraction of the samples. We sampled 5 L of water at each sampling point (a sample simple for each point), and the TSS and DOM analyses were performed in triplicate [26].
The UAV used for image acquisition was a hexacopter (Figure 4), on which a Canon RGB sensor was mounted. Figure 4 presents the remote sensing system used in the acquisition of spectral images.
The UAV was flown at an altitude of 150 m with an average speed of 5 m/s, generating 46 images obtained with the Canon sensor. To georeference the images, we established a network with six ground control points (GCPs), evenly distributed throughout the study area, and obtained their geodetic coordinates by the real-time kinematic (RTK) method using the Global Navigation Satellite System (GNSS). We adopted SIRGAS 2000 (Geocentric Reference System for the Americas) as the referencing system in Universal Transverse Mercator (UTM) projection, zone 22S.
Processing of the UAV data was a simple, three-step workflow. In the first step, we created a controlled mosaic and determined the accuracy of the model after applying GCPs. In the second step, we created a point cloud from the pixels of the images using the structure from motion technique [27]. In the third step, we transformed the point cloud into a digital surface model (DSM) used as the basis for the orthorectification process. Thus, in terms of image quality, we generated an orthophoto with a pixel size of 0.04 × 0.04 m2.

2.3. Proposal of an Artificial Neural Network to Predict Spectral Bands of Images Acquired by a UAV

Neural networks were used in the following two modeling stages: ANN training and validation. Our working hypothesis is that it is possible to structure an ANN capable of predicting the spectral bands used in water quality monitoring in dams/lakes using the three bands of the Canon sensor (blue (B), green (G), and red (R)) as the ANN input. In order to structure the output of the network we used the bands b4 (0.64–0.67 μm), b5 (0.85–0.88 μm), and b9 (1.36–1.38 μm) of the LANDSAT 8 (OLI) satellite. The images were acquired from the United States National Aeronautics and Space Administration. We performed atmospheric correction of the images using the method described in [28,29]. We chose to predict bands b4, b5, and b9 because they are the most appropriate spectral bands to study water quality in dams and lakes [13,19].
As the images were georeferenced with the GNSS system, we standardized the pixel sizes of the images from the Canon sensor (0.04 m) based on those of the three bands of the satellite LANDSAT 8 (OLI) (30 m). Each pixel of the ANN input image corresponded to that of the output image. Thus, the B, G, and R Canon sensor images were the ANN input, and the b4, b5, and b9 LANDSAT 8 (OLI) images were the output. For ANN training, 80% of the pixels were used and, for validation, 20% were used, with randomly-defined data [30].
The ANN was trained using Adaptive Moment Estimation (Adam), which is an optimization algorithm that can be used instead of the classical stochastic gradient descent procedure to update network weights iteratively. Adam is simple to implement, computationally efficient, requires few memory requirements, has a good convergence rate, and is appropriate for machine learning problems [31,32].
During this training phase, as in most studies with neural networks [32], the definition of the best ANN topology was achieved through numerous experiments, in which the numbers of intermediate layers and neurons per layer, function of activation, and number of training cycles were tested. The activation function is important because it controls the way each individual neuron is activated and propagates the information through to the next layer. We used the sigmoid activation function because it increases the effectiveness of the backpropagation algorithm; thus, it is suitable for use in multilayer networks [33]. To ensure that the ANN achieved its results with the lowest possible amount of training effort, we normalized the intensity of the pixels to 0–1 [34].
After determining the topology that provided the lowest mean square error during supervised training, we validated the model with 20% of the data so that we can compare the known levels of intensity of the LANDSAT 8 (OLI) bands with those provided by the ANN.

2.4. Data Analysis

To evaluate the consistency of the results, we studied the correlation between the intensities of the known pixels from the ANN for each LANDSAT 8 (OLI) band. We calculated the coefficient of determination (R²) at a significance level of 5%.
Once the neural network was validated, we re-introduced the three-band Canon sensor images with their original resolution of 0.04 m2. By introducing the Canon sensor’s images into the ANN with their original resolution, we obtained the b4, b5, and b9 LANDSAT 8 (OLI) 8 images with pixels of 0.04 m2 as outputs. We georeferenced these images and generated the NDVI (Normalized Difference Vegetation Index) (Equation (1)) and NDWI (Normalized Difference Water Index) (Equation (2)):
N D V I =   b 5 b 4 b 5 + b 4
N D W I =   G b 5 G + b 5
NDVI was initially developed to evaluate vegetation biomass with higher reflectance in the near-infrared spectral range rather than in the red band. Thus, NDVI values vary between −1 and +1, with negative values corresponding to the presence of water, values close to zero representing exposed soil, and positive values indicating vegetated surfaces (ranging from not dense to dense, represented by an increase in the value of the index). NDWI was developed to measure the spectral characteristics of water using the green band in place of the red band, allowing water to be identified in terrestrial and vegetated surfaces. The NDWI values also vary from −1 to +1, with positive values indicating water bodies and negative values indicating vegetation and soil [35].
We plotted the sample points at which we obtained the total suspended solids (TSS) and dissolved organic matter (DOM) concentrations on the georeferenced NDVI and NDWI images. Four sample points were considered outliers because they were outside the range of two standard deviations (above or below) of the mean of the set [19]. They were excluded from the set, and the remaining seventeen data points for the analysis of correlations. We determined the correlations between TSS and DOM and the NDVI and NDWI values obtained from the images generated by the ANN.

3. Results and Discussion

In this section, we present the results of this research in three sub-sections: the results of laboratory analysis, the results of the ANN processing, and the correlations between TSS and DOM and the NDVI and NDWI values obtained from the images generated by the ANN.

3.1. Results of Laboratory Analysis

The laboratory analysis generated satisfactory results to determine the concentrations of TSS and DOM in water samples. Table 1 summarizes these results.

3.2. Results of Artificial Neural Network Processing

During the training process, using 80% of the pixels of LANDSAT (OLI) 8 and Canon sensor images, we tested different network structures, and varying numbers of intermediate layers, neurons per layer, and training cycles. Such a practical experimentation is common in studies that use neural networks [32,34] and the best topology identified was 3-4-6-3, with 100 training cycles at a computational cost of 288 s, considering the computer used had the following configuration: processor—Intel® Core ™ i7-3770 CPU @ 3.40 GHz ×8; memory—12 GB DDR3 1333 MHz RAM.
Once the best topology for the ANN was defined, we validated the ANN using 20% of the data from the pixels of LANDSAT (OLI) 8 and the Canon sensor images and compared the intensity of the known pixels (bands b4, b5, and b9) with those generated by the ANN. We adopted 80% of the data for training and 20% for validation, according to [30]. We calculated the mean square error and conducted a linear regression between the known pixel intensity values and those predicted by the ANN for the three output bands of the network. It is important to emphasize that we also evaluated nonlinear regression models. However, when evaluating the result of the adjustment by the least square method, the linear model was more effective than nonlinear models, presenting a higher value for the determination coefficient R2 and a lower value for the mean square error (MSR). The results are presented in Table 2.
We used the Canon sensor images with their original resolutions (pixel of 0.04 m2) and introduced the three bands (B, G, R) in the trained and validated ANN, obtaining the three bands, b4, b5, and b9, with a computational cost of 30.5 s that generated NDVI and NDWI images (Figure 5).

3.3. Correlations

The results of the correlation analysis for NDVI and NDWI with the concentrations of TSS and DOM are presented in Table 3. The best results for the coefficient of determination (R2), root mean square error (RMSE), and standard deviation (SD) were from the polynomial regressions.
The results obtained and presented in Table 3 are consistent with those reported in previous studies. In the literature, R2 values of the correlations between NDVI/NDWI and TSS were around 0.60 [1,21]. In our study, R² values were 0.63 (NDVI) and 0.77 (NDWI), respectively. R2 values are higher in our study because our method allows the prediction of spectral bands with the same geometric resolution (pixel with centimetric size) as the resolution of bands R, G, and B used in the ANN input. For the correlations of NDVI/NDWI with DOM, our results are promising, but we did not find any studies with which to compare our results.
The consistency of the values shown in the Table 3 is evidenced when the scatter plot between the known and simulated values of TSS and DOM is established (Figure 6).
NDVI and NDWI are commonly used for images with near-infrared bands of approximately 800 to 900 nm [36]. Previous studies [37] mapped the spectral behavior of turbid water and water rich in suspended solids, and the spectral signature for TSS had reflectance peaks close to 550, 650, and 700 nm. The same patterns were found by other researchers [19]. However, water bodies in those studies have OACs with different characteristics, which make it possible to evaluate other compositions of spectral bands to determine NDVI and NDWI. For example, [38] proposed a method to obtain NDWI spectral bands ranging from 0.86 μm to 1.24 μm. In our lake, the TSS and DOM average concentrations are smaller than those in previous studies. Based on the work of [38] we modified Equations (1) and (2) and determined NDVI and NDWI using Equations (3) and (4):
N D V I =   b 9 b 4 b 9 + b 4
N D W I =   G b 9 G + b 9
In Equations (3) and (4), we replaced band b5 (0.85–0.88 μm) with band b9 (1.36–1.38 μm) and extended the spectral range of the infrared to show the presence of TSS and DOM because mean concentrations of these elements were small (Figure 7).
Based on the NDVI and NDWI images shown in Figure 7, the results of correlation of these indices with TSS and DOM are presented in Table 4.
Replacing b5 by b9 in NDVI and NDWI equations slightly increased the R² of the correlation of TSS and DOM with NDVI, and slightly decreased RMSE and SD (Table 4). The reason for the small increase is that band 9 provided a greater contrast of NDVI, indicating vegetative areas and the water body in Figure 7A. R2 of the correlation of TSS and DOM with NDWI showed a slight decrease. In Figure 7B, if we use the spectral band b9 instead of the b5 band, the contrast of the NDWI image worsened (compared to Figure 5B). The consistency of the results in Table 4 is also verified when the scatter plot between the known and simulated values of TSS and DOM is established (Figure 8).
We synthesized all these results in Figure 9. The correlation of TSS with NDVI (Figure 9A) obtained with bands b5 and b9 show similar values for R2 (line), with values greater than 0.6, and for RMSE (bar chart), with values of approximately 1.7 mg/L. The correlation of TSS with NDWI showed a significant improvement in relation to the NDVI, with an increase of R² (around 0.77) and a decrease in RMSE (lower than 1.4 mg/L), also with similar results for bands b5 and b9. The results related to DOM (Figure 9B) followed the same trend as those related to TSS. We observed stronger correlations with NDWI (R2 of 0.60) and a significant decrease in RMSE compared to the correlations with NDVI.
The improvement in the results of the correlation with NDWI is in agreement with the results of Mohanty et al. [39] that used Indian remote sensing satellite images (IRS-IB) and the obtained R2 values between TSS and NDWI were around 0.6. Our results were higher than those of [39] because our method allows generating multispectral images (centimeter pixel size) from an expert system based on an ANN. The pixel size of multispectral images predicted by the ANN will always be equal to the pixel size of the images from a low-cost sensor (spectral bands R, G, B) because these are the bands used as ANN input data. The coefficients of correlation of TSS with NDVI in our study are also higher than those in Masocha et al. [40]. Our results for DOM were consistent, but we did not find any studies with which to compare our results.
Even with the good results found in this article, the proposed method has limitations that are highlighted below:
  • The water sampling must be performed in periods under good light conditions. This condition is important to obtain high-quality UAV images and spectral bands of satellites without clouds;
  • It can be difficult to obtain orthophotos by UAV for large dams/lakes. Generally, water bodies are quite homogeneous, making it difficult to determine the homologous points between the images during the processing. The homologous points are fundamental for creating a point cloud from the pixels of the images using the structure from motion technique and to transform the point cloud into a digital surface model (DSM) used as the basis for the orthorectification process [27];
  • Our approach did not consider the seasonal variation of TSS and DOM. The homogeneity of the physical and chemical conditions of the water is a reflection of the temporal conditions of the bodies of water. Thus, our approach can be improved by incorporating in the neural network input the climate variable characteristics of the study site (for example, wind velocity, rainfall, temperature, etc.).

4. Conclusions

The results of this study demonstrate the applicability of a methodological model that uses an expert ANN system to monitor water quality by UAVs. We corroborate our hypothesis because we were able to predict the bands of a multispectral sensor through the insertion of bands from a low-cost sensor using an ANN validation process.
The results indicate that our proposed method for predicting multispectral bands based on an expert system and a low-cost sensor mounted on a UAV may be a useful tool for monitoring water quality for sites with low concentrations of TSS and DOM. Traditional methods of water quality monitoring in dams/lakes are time-consuming, costly, and often impractical in large bodies of water. Thus, the present study provides a better alternative for monitoring water quality using multispectral images.
For the water body we studied, the best spectral range to monitor TSS and DOM concentrations was 530–1380 nm. We established the correlations of TSS and DOM with NDVI and NDWI obtained from the spectral bands predicted by the ANN (b4: 0.64–0.67 μm), (b5: 0.85–0.88 μm), and (b9: 1.36–1.38 μm). The correlation coefficient for TSS was higher than that found for DOM and the R2 values were the highest for the correlations with NDWI. For correlations between TSS and the indices, R2 in this study was higher than that reported in the literature. Regarding the correlations with DOM concentrations, our results were consistent, but we could not find any studies to make comparisons.
The method described in this article can be applied to monitor other water bodies. With this method, it is easy to mount low-cost sensors on any unmanned aerial vehicle, making it possible to monitor water quality in dams and lakes.

Acknowledgments

M.R.V. thanks the Brazilian Council for Scientific and Technological Development (CNPq) for the research grant (process 309399/2014-9). L.V.S. thanks PROSUP/CAPES for the financial support of the PhD scholarship. L.S.K., J.M.S., and I.É.K. thank PROSUP/CAPES for the financial support of the MSc scholarship. T.T.G. thanks DS/CAPES for the financial support of the MSc scholarship. J.M.S. thanks CNPq for the financial support of the MSc scholarship.

Author Contributions

Emilie C. Koste, Tainá T. Guimarães, Mauricio R. Veronez, Leonardo C. Inocencio, Luiz Gonzaga Jr., and Fabiane Bordin were responsible for collecting and processing the spectral data and the images obtained from the UAV. Lucas S. Kupssinskü, Juarez M. da Silva, Laís V. de Souza, William F. M. Oliverio, Rogélio S. Jardim, Ismael É. Koch, and Jonas G. de Souza computationally implemented the artificial neural network and performed the training and validation steps. Tainá T. Guimarães, Emilie C. Koste, and Frederico F. Mauad were responsible for water collection at the predetermined sampling points and for the laboratory analysis. Mauricio R. Veronez, Emilie C. Koste, and Tainá T. Guimarães wrote the paper. Mauricio R. Veronez and Frederico F. Mauad reviewed the paper. All the authors have read and approved the paper final version.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Roig, H.L.; Ferreira, A.M.R.; Menezes, P.H.B.J.; Marotta, G.S. Uso de câmeras de baixo custo acopladas a veículos aéreos leves no estudo do aporte de sedimentos no Lago Paranoá. In Proceedings of the Anais XVI Simpósio Brasileiro de Sensoriamento Remoto—SBSR, Foz do Iguaçu, Brazil, 13–18 April 2013; INPE: Foz do Iguaçu, Brazil, 2013; pp. 9332–9339. (In Portuguese). [Google Scholar]
  2. George, D.G. The airborne remote sensing of phytoplankton chlorophyll in the lakes and tarns of the English Lake District. Int. J. Remote Sens. 1997, 18, 1961–1975. [Google Scholar] [CrossRef]
  3. Dekker, A.G.; Vos, R.J.; Peters, S.W.M. Analytical algorithms for lake water TSM estimation for retrospective analyses of TM and SPOT sensor data. Int. J. Remote Sens. 2002, 23, 15–35. [Google Scholar] [CrossRef]
  4. Wu, J.-L.; Ho, C.-R.; Huang, C.-C.; Srivastav, A.L.; Tzeng, J.-H.; Lin, Y.-T. Hyperspectral Sensing for Turbid Water Quality Monitoring in Freshwater Rivers: Empirical Relationship between Reflectance and Turbidity and Total Solids. Sensors 2014, 14, 22670–22688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Han, L.; Rundquist, D.C. Spectral characterization of suspended sediments generated from two texture classes of clay soil. Int. J. Remote Sens. 1996, 17, 643–649. [Google Scholar] [CrossRef]
  6. Gitelson, A.; Garbuzov, G.; Szilagyi, F.; Mittenzwey, K.-H.; Karnieli, A.; Kaiser, A. Quantitative remote sensing methods for real-time monitoring of inland waters quality. Int. J. Remote Sens. 1993, 14, 1269–1295. [Google Scholar] [CrossRef]
  7. Santini, F.; Alberotanza, L.; Cavalli, R.M.; Pignatti, S. A two-step optimization procedure for assessing water constituent concentrations by hyperspectral remote sensing techniques: An application to the highly turbid Venice lagoon waters. Remote Sens. Environ. 2010, 114, 887–898. [Google Scholar] [CrossRef]
  8. Thiemann, S.; Kaufmann, H. Determination of Chlorophyll Content and Trophic State of Lakes Using Field Spectrometer and IRS-1C Satellite Data in the Mecklenburg Lake District, Germany. Remote Sens. Environ. 2000, 73, 227–235. [Google Scholar] [CrossRef]
  9. Doxaran, D.; Froidefond, J.-M.; Lavender, S.; Castaing, P. Spectral signature of highly turbid waters: Application with SPOT data to quantify suspended particulate matter concentrations. Remote Sens. Environ. 2002, 81, 149–161. [Google Scholar] [CrossRef]
  10. Teodoro, A.C.; Veloso-Gomes, F.; Gonçalves, H. Statistical Techniques for Correlating Total Suspended Matter Concentration with Seawater Reflectance Using Multispectral Satellite Data. J. Coast. Res. 2008, 24, 40–49. [Google Scholar] [CrossRef]
  11. Liu, Z.; Wu, J.; Yang, H. Developing unmanned airship onboard multispectral imagery system for quick-response to drinking water pollution. In Proceedings of the International Society for Optics and Photonics, Yichang, China, 30 October–1 November 2009; Volume 7494, p. 74940L. [Google Scholar]
  12. Ferreira, M.S.; Ennes, R.; Galo, M.L.B.T. Estudo do comportamento espectral das águas da planície de inundação do Alto Rio Paraná baseada em técnicas de análise de dados espectrais. In Proceedings of the 3rd Simpósio Brasileiro de Ciências Geodésicas e Tecnologias da Geoinformação, Recife, Brazil, 27–30 July 2010; pp. 1–9. (In Portuguese). [Google Scholar]
  13. Campbell, G.; Phinn, S.R.; Dekker, A.G.; Brando, V.E. Remote sensing of water quality in an Australian tropical freshwater impoundment using matrix inversion and MERIS images. Remote Sens. Environ. 2011, 115, 2402–2414. [Google Scholar] [CrossRef] [Green Version]
  14. Gago, J.; Douthe, C.; Coopman, R.E.; Gallego, P.P.; Ribas-Carbo, M.; Flexas, J.; Escalona, J.; Medrano, H. UAVs challenge to assess water stress for sustainable agriculture. Agric. Water Manag. 2015, 153, 9–19. [Google Scholar] [CrossRef]
  15. Hardin, P.J.; Hardin, T.J. Small-Scale Remotely Piloted Vehicles in Environmental Research. Geogr. Compass 2010, 4, 1297–1311. [Google Scholar] [CrossRef]
  16. Berni, J.A.J.; Zarco-Tejada, P.J.; Suarez, L.; Fereres, E. Thermal and Narrowband Multispectral Remote Sensing for Vegetation Monitoring from an Unmanned Aerial Vehicle. IEEE Trans. Geosci. Remote Sens. 2009, 47, 722–738. [Google Scholar] [CrossRef] [Green Version]
  17. Cândido, A.K.A.A.; Filho, A.C.P.; Haupenthal, M.R.; da Silva, N.M.; de Sousa Correa, J.; Ribeiro, M.L. Water Quality and Chlorophyll Measurement through Vegetation Indices Generated from Orbital and Suborbital Images. Water. Air Soil Pollut. 2016, 227, 224. [Google Scholar] [CrossRef]
  18. Su, T.-C.; Chou, H.-T. Application of Multispectral Sensors Carried on Unmanned Aerial Vehicle (UAV) to Trophic State Mapping of Small Reservoirs: A Case Study of Tain-Pu Reservoir in Kinmen, Taiwan. Remote Sens. 2015, 7, 10078–10097. [Google Scholar] [CrossRef]
  19. Guimarães, T.T.; Veronez, M.R.; Koste, E.C.; Gonzaga, L., Jr.; Bordin, F.; Inocencio, L.C.; Larocca, A.P.C.; de Oliveira, M.Z.; Vitti, D.C.; Mauad, F.F. An Alternative Method of Spatial Autocorrelation for Chlorophyll Detection in Water Bodies Using Remote Sensing. Sustainability 2017, 9, 416. [Google Scholar] [CrossRef]
  20. Flynn, K.F.; Chapra, S.C. Remote Sensing of Submerged Aquatic Vegetation in a Shallow Non-Turbid River Using an Unmanned Aerial Vehicle. Remote Sens. 2014, 6, 12815–12836. [Google Scholar] [CrossRef]
  21. Luo, J.; Li, X.; Ma, R.; Li, F.; Duan, H.; Hu, W.; Qin, B.; Huang, W. Applying remote sensing techniques to monitoring seasonal and interannual changes of aquatic vegetation in Taihu Lake, China. Ecol. Indic. 2016, 60, 503–513. [Google Scholar] [CrossRef]
  22. Silva, W.F.; Silva, L.S.; Malta, É.A.; de Oliveira Gondim, R.; Warren, M.S. Avaliação de uso de Veículo Aéreo Não Tripulado—VANT em atividades de fiscalização da Agência Nacional de Águas. In Proceedings of the Anais XVII Simpósio Brasileiro de Sensoriamento Remoto—SBSR, João Pessoa, Brazil, 25–29 April 2015; INPE: João Pessoa, Brazil, 2015; pp. 1791–1798. (In Portuguese). [Google Scholar]
  23. Tauro, F.; Olivieri, G.; Petroselli, A.; Porfiri, M.; Grimaldi, S. Flow monitoring with a camera: A case study on a flood event in the Tiber River. Environ. Monit. Assess. 2016, 188, 118. [Google Scholar] [CrossRef] [PubMed]
  24. Tamminga, A.; Hugenholtz, C.; Eaton, B.; Lapointe, M. Hyperspatial Remote Sensing of Channel Reach Morphology and Hydraulic Fish Habitat Using an Unmanned Aerial Vehicle (UAV): A First Assessment in the Context of River Research and Management. River Res. Appl. 2015, 31, 379–391. [Google Scholar] [CrossRef]
  25. Zang, W.; Lin, J.; Wang, Y.; Tao, H. Investigating small-scale water pollution with UAV Remote Sensing Technology. In Proceedings of the World Automation Congress, Puerto Vallarta, Mexico, 24–28 June 2012; pp. 1–4. [Google Scholar]
  26. American Public Health Association (APHA). Standard Methods for Examination of Water and Wastewater; American Public Health Association: Washington, DC, USA, 1995. [Google Scholar]
  27. Tonkin, T.N.; Midgley, N.G. Ground-Control Networks for Image Based Surface Reconstruction: An Investigation of Optimum Survey Designs Using UAV Derived Imagery and Structure-from-Motion Photogrammetry. Remote Sens. 2016, 8, 786. [Google Scholar] [CrossRef]
  28. Gusso, A.; Cafruni, C.; Bordin, F.; Veronez, M.R.; Lenz, L.; Crija, S. Multi-Temporal Patterns of Urban Heat Island as Response to Economic Growth Management. Sustainability 2015, 7, 3129–3145. [Google Scholar] [CrossRef]
  29. Poursanidis, D.; Chrysoulakis, N.; Mitraka, Z. Landsat 8 vs. Landsat 5: A comparison based on urban and peri-urban land cover mapping. Int. J. Appl. Earth Obs. Geoinf. 2015, 35, 259–269. [Google Scholar] [CrossRef]
  30. Haykin, S.S. Neural Networks and Learning Machines, 3rd ed.; Pearson: Upper Saddle River, NJ, USA, 2009; p. 936. [Google Scholar]
  31. Kingma, D.P.; Ba, J. Adam: A Method for Stochastic Optimization. arXiv, 2014; arXiv:14126980. [Google Scholar]
  32. Brownlee, J. Deep Learning with Python: Understand Your Data, Create Accurate Models an Work Projects End-to-End; v1.4; Jason Brownlee: Melbourne, Australia, 2016; p. 170. [Google Scholar]
  33. Chen, Y.; Bastani, F. ANN with two-dendrite neurons and its weight initialization. In Proceedings of the IJCNN International Joint Conference on Neural Networks, Baltimore, MD, USA, 7–11 June 1992; Volume 3, pp. 139–146. [Google Scholar]
  34. Sola, J.; Sevilla, J. Importance of input data normalization for the application of neural networks to complex industrial problems. IEEE Trans. Nucl. Sci. 1997, 44, 1464–1468. [Google Scholar] [CrossRef]
  35. Mcfeeters, S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  36. Ke, Y.; Im, J.; Lee, J.; Gong, H.; Ryu, Y. Characteristics of Landsat derived NDVI by comparison with multiple satellite sensors and in-situ observations. Remote Sens. Environ. 2015, 164, 298–313. [Google Scholar] [CrossRef]
  37. Sváb, E.; Tyler, A.N.; Preston, T.; Présing, M.; Balogh, K.V. Characterizing the spectral reflectance of algae in lake waters with high suspended sediment concentrations. Int. J. Remote Sens. 2005, 26, 919–928. [Google Scholar] [CrossRef]
  38. Gao, B. NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 1996, 58, 257–266. [Google Scholar] [CrossRef]
  39. Mohanty, P.K.; Pal, S.R.; Mishra, P.K. Monitoring Ecological Conditions of a Coastal Lagoon using IRS Data: A Case Study in Chilka, East Coast of India. J. Coast. Res. 2001, 34, 459–469. [Google Scholar]
  40. Masocha, M.; Murwira, A.; Magadza, C.H.D.; Hirji, R.; Dube, T. Remote sensing of surface water quality in relation to catchment condition in Zimbabwe. Phys. Chem. Earth Parts 2017. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the proposed method.
Figure 1. Flowchart of the proposed method.
Sensors 18 00159 g001
Figure 2. Study area.
Figure 2. Study area.
Sensors 18 00159 g002
Figure 3. Positions of the sampling points.
Figure 3. Positions of the sampling points.
Sensors 18 00159 g003
Figure 4. Hexacopter used for lake mapping.
Figure 4. Hexacopter used for lake mapping.
Sensors 18 00159 g004
Figure 5. Images generated by the compositions of the spectral bands obtained by ANN: (A) NDVI image; and (B) NDWI image.
Figure 5. Images generated by the compositions of the spectral bands obtained by ANN: (A) NDVI image; and (B) NDWI image.
Sensors 18 00159 g005
Figure 6. Scatter plot of SST and DOM known and simulated considering the b4 and b5 bands in obtaining the NDVI and NDWI indices.
Figure 6. Scatter plot of SST and DOM known and simulated considering the b4 and b5 bands in obtaining the NDVI and NDWI indices.
Sensors 18 00159 g006
Figure 7. Images generated by the compositions of the spectral band b9 instead of the band b5 obtained by ANN: (A) NDVI image; and (B) NDWI image.
Figure 7. Images generated by the compositions of the spectral band b9 instead of the band b5 obtained by ANN: (A) NDVI image; and (B) NDWI image.
Sensors 18 00159 g007
Figure 8. Scatter plot of SST and DOM known and simulated considering the b4 and b9 bands in obtaining the NDVI and NDWI indices.
Figure 8. Scatter plot of SST and DOM known and simulated considering the b4 and b9 bands in obtaining the NDVI and NDWI indices.
Sensors 18 00159 g008
Figure 9. Comparison between the correlation results.
Figure 9. Comparison between the correlation results.
Sensors 18 00159 g009
Table 1. Results of TSS and DOM analysis.
Table 1. Results of TSS and DOM analysis.
ParameterAverage (mg/L)Minimum/Location (mg/L)Maximum/Location (mg/L)Standard Deviation (mg/L)
TSS13.659.33 (P02)20 (P15)3.07
DOM38.054.67 (P19)175 (P12)39.76
Table 2. Linear regression between known LANDSAT 8 (OLI) pixel intensity and the prediction from the proposed ANN model.
Table 2. Linear regression between known LANDSAT 8 (OLI) pixel intensity and the prediction from the proposed ANN model.
Pixel IntensityLinear Equation (1)R2 (2)MSR (3)
b4y = 0.6271x + 31.21250.82106.80%
b5y = 0.8051x + 20.07810.840118.15%
b9y = 0.7097x + 25.45210.812010.37%
(1) Linear regression conducted between known pixel intensity and ANN-predicted pixel intensity. (2) Determination coefficient. (3) Mean square error calculated by dividing the sum of the squared differences between known pixel intensities and ANN-predicted values per mean pixel intensity value.
Table 3. Polynomial regression results.
Table 3. Polynomial regression results.
CorrelationPolynomial EquationR2 (1)RMSE (2)SD (2)
TSS × NDVITSS = 26.7 × NDVI2 + 30.9 × NDVI + 20.00.631.692.50
TSS × NDWITSS = 82.2 × NDWI2 − 131.5 × NDWI + 63.50.771.342.36
DOM × NDVIDOM = 208.573 × NDVI3 + 229.6 × NDVI2 + 51.3 × NDVI + 27.00.525.483.87
DOM × NDWIDOM = −1929.1 × NDWI3 + 3940.8 × NDWI² − 2592.4 × NDWI + 576.70.607.355.32
(1) Correlations between the water quality parameters (TSS and DOM) and NDVI and NDVI indices obtained with the spectral bands b4 and b5 predicted by ANN. (2) Root mean square error (RMSE) and standard deviation (SD) are expressed as mg/L.
Table 4. Polynomial regression results.
Table 4. Polynomial regression results.
CorrelationPolynomial EquationR2 (1)RMSE (2)SD (2)
TSS × NDVITSS = 45.4 × NDVI2 + 43.1 × NDVI + 20.90.651.652.33
TSS × NDWITSS = 68.7 × NDWI2 − 111.2 × NDWI + 56.10.761.332.54
DOM × NDVIDOM = 244.9 × NDVI3 + 186.2 × NDVI2 + 7.0 × NDVI + 21.80.545.034.47
DOM × NDWIDOM = −2119.5 × NDWI3 + 4559.1 × NDWI2 − 2760.4 × NDWI + 603.60.594.235.28
(1) Correlations between the water quality parameters (TSS and DOM) and NDVI and NDVI indices obtained with the spectral bands b4 and b9 predicted by ANN. (2) Root mean square error (RMSE) and standard deviation (SD) are expressed as mg/L.

Share and Cite

MDPI and ACS Style

R. Veronez, M.; Kupssinskü, L.S.; T. Guimarães, T.; Koste, E.C.; Da Silva, J.M.; De Souza, L.V.; Oliverio, W.F.M.; Jardim, R.S.; Koch, I.É.; De Souza, J.G.; et al. Proposal of a Method to Determine the Correlation between Total Suspended Solids and Dissolved Organic Matter in Water Bodies from Spectral Imaging and Artificial Neural Networks. Sensors 2018, 18, 159. https://doi.org/10.3390/s18010159

AMA Style

R. Veronez M, Kupssinskü LS, T. Guimarães T, Koste EC, Da Silva JM, De Souza LV, Oliverio WFM, Jardim RS, Koch IÉ, De Souza JG, et al. Proposal of a Method to Determine the Correlation between Total Suspended Solids and Dissolved Organic Matter in Water Bodies from Spectral Imaging and Artificial Neural Networks. Sensors. 2018; 18(1):159. https://doi.org/10.3390/s18010159

Chicago/Turabian Style

R. Veronez, Maurício, Lucas S. Kupssinskü, Tainá T. Guimarães, Emilie C. Koste, Juarez M. Da Silva, Laís V. De Souza, William F. M. Oliverio, Rogélio S. Jardim, Ismael É. Koch, Jonas G. De Souza, and et al. 2018. "Proposal of a Method to Determine the Correlation between Total Suspended Solids and Dissolved Organic Matter in Water Bodies from Spectral Imaging and Artificial Neural Networks" Sensors 18, no. 1: 159. https://doi.org/10.3390/s18010159

APA Style

R. Veronez, M., Kupssinskü, L. S., T. Guimarães, T., Koste, E. C., Da Silva, J. M., De Souza, L. V., Oliverio, W. F. M., Jardim, R. S., Koch, I. É., De Souza, J. G., Gonzaga, L., Jr., Mauad, F. F., Inocencio, L. C., & Bordin, F. (2018). Proposal of a Method to Determine the Correlation between Total Suspended Solids and Dissolved Organic Matter in Water Bodies from Spectral Imaging and Artificial Neural Networks. Sensors, 18(1), 159. https://doi.org/10.3390/s18010159

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop