Next Article in Journal
Using GIS, Remote Sensing, and Machine Learning to Highlight the Correlation between the Land-Use/Land-Cover Changes and Flash-Flood Potential
Previous Article in Journal
Remote Sensing Retrieval of Total Phosphorus in the Pearl River Channels Based on the GF-1 Remote Sensing Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neural Network Reflectance Prediction Model for Both Open Ocean and Coastal Waters

1
Joint Center for Earth Systems Technology, Department of Physics, University of Maryland, Baltimore County, 1000 Hilltop Circle Baltimore, Baltimore, MD 21250, USA
2
SSAI, NASA Goddard Space Flight Center, Code 616, Greenbelt, MD 20771, USA
3
MS 475 NASA Langley Research Center, Hampton, VA 23681-2199, USA
4
NASA Goddard Space Flight Center, Code 616, Greenbelt, MD 20771, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(9), 1421; https://doi.org/10.3390/rs12091421
Submission received: 1 April 2020 / Revised: 23 April 2020 / Accepted: 28 April 2020 / Published: 30 April 2020

Abstract

:
Remote sensing of global ocean color is a valuable tool for understanding the ecology and biogeochemistry of the worlds oceans, and provides critical input to our knowledge of the global carbon cycle and the impacts of climate change. Ocean polarized reflectance contains information about the constituents of the upper ocean euphotic zone, such as colored dissolved organic matter (CDOM), sediments, phytoplankton, and pollutants. In order to retrieve the information on these constituents, remote sensing algorithms typically rely on radiative transfer models to interpret water color or remote-sensing reflectance; however, this can be resource-prohibitive for operational use due to the extensive CPU time involved in radiative transfer solutions. In this work, we report a fast model based on machine learning techniques, called Neural Network Reflectance Prediction Model (NNRPM), which can be used to predict ocean bidirectional polarized reflectance given inherent optical properties of ocean waters. This supervised model is trained using a large volume of data derived from radiative transfer simulations for coupled atmosphere and ocean systems using the successive order of scattering technique (SOS-CAOS). The performance of the model is validated against another large independent test dataset generated from SOS-CAOS. The model is able to predict both polarized and unpolarized reflectances with an absolute error (AE) less than 0.004 for 99% of test cases. We have also shown that the degree of linear polarization (DoLP) for unpolarized incident light can be predicted with an AE less than 0.002 for 99% of test cases. In general, the simulation time of SOS-CAOS depends on optical depth, and required accuracy. When comparing the average speeds of the NNRPM against the SOS-CAOS model for the same parameters, we see that the NNRPM is able to predict the Ocean BRDF 6000 times faster than SOS-CAOS. Both ultraviolet and visible wavelengths are included in the model to help differentiate between dissolved organic material and chlorophyll in the study of the open ocean and the coastal zone. The incorporation of this model into the retrieval algorithm will make the retrieval process more efficient, and thus applicable for operational use with global satellite observations.

1. Introduction

The study of global ocean color measurements can help in understanding the ocean ecology and biogeochemistry which can subsequently help better understand the carbon cycle and its impact on climate change. Ocean polarized reflectances obtained from satellite sensors contain a lot of information about the constituents like colored dissolved organic matter (CDOM), sediments, phytoplankton, and pollutants in water. In order to retrieve this information, ocean color remote sensing algorithms retrieve remote sensing reflectances (Rse,sr−1) from the sensor measurement at the top of the atmosphere (TOA) through atmospheric correction, which removes the contribution of atmospheric and ocean surface effects from the total signal [1,2,3].
Once remote sensing reflectance is obtained, empirical and model based approaches [4,5,6,7,8,9,10,11] are developed to retrieve the inherent optical properties (IOPs). One such empirical approach is the band ratio method which statistically relates Rrs at two or more wavelengths to a particular water constituent. The algorithm quantifies water constituents as a function of blue and green light emanating from the ocean surface [4,12]. Another empirical approach is based on a line-height color index that is the difference between Rse in the green band and a reference formed linearly between remote sensing reflectance in the blue and red band [5]. The second approach performs better in areas with adjacency effects, such as sun glint or cloud edges. Both empirical approaches are successful for open ocean where the phytoplankton and detrital are the main constituents of water, co-varying with each other [13]. However, they often perform less successfully for coastal ocean where the optical constituents like CDOM and suspended sediments do not co-vary with water primary productivity.
The model based approach uses bio-optical models, most of which are empirically derived, to describe the relationship between in-water constituents and inherent optical properties (IOP’s). Radiative transfer models are then used to relate the reflectance spectra to the desired in-water constituents. Semi-analytical models fall within the category of model based algorithms. There are several semi-analytical models developed with different assumptions and limitations [6]. They essentially follow different empirical relationships of spectral shape functions or eigenvectors of the constituent absorption and scattering components. Thereafter, the optimized magnitudes or eigenvalues representing the relative contributions of each constituent is obtained by minimizing the difference between modelled and measured reflectance values (Rse) [14,15]. Those algorithms cannot predict the polarization of ocean reflectance and are not valid in the ultraviolet wavelength range.
Both empirical and model based methods depend on a reliable atmospheric correction algorithm. However, due to the presence of absorbing aerosols and complex optical properties of coastal waters, it is difficult to achieve successful atmospheric correction over coastal waters or when absorbing aerosols are involved [16,17]. To overcome the limitations of atmospheric correction in coastal and inland waters involving absorbing aerosols, the community proposes to use multi-angle multi-wavelength, polarized measurements to characterize aerosol and hydrosol properties. The polarized signal is measured by several multiangular polarimeters like POLDER [18], the Research Scanning Polarimeter (RSP) [19], the Airborne Multiangle SpectroPolarimetric Imager (AirMSPI) [20], the Spectropolarimeter for Planetary EXploration (SPEX) [21] and the HyperAngular Rainbow Polarimeter (HARP) [22]. it is worth noting that NASA’s Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission plans to carry SPEXone and HARP2 as well as the Ocean Color Instrument (OCI), which will acquire a plethora of unprecedented co-located datasets of polarimeter and ocean color measurements [23].
In order to take advantage of the polarimeter measurements, joint retrieval algorithms are developed to simultaneously retrieve information of aerosols and hydrosols [24,25,26]. The bio-optical models used by the aforementioned joint retrieval algorithms are for open ocean and thus are ill-suited for coastal waters. Recently, the Multi-Angle Polarimetric Ocean coLor (MAPOL) retrieval algorithm [27] has been designed specifically for the optically complex coastal waters and aerosols for improved atmospheric correction. MAPOL has been demonstrated over both open ocean and coastal waters, which shows promising applications in achieving atmospheric correction over complex scenes [27]. The joint retrieval algorithms, including MAPOL, need to call radiative transfer models iteratively in order to minimize the difference between measurements and model prediction, which is computationally expensive and difficult to achieve operational production. It is critical to develop a fast ocean reflectance model that can be used in the joint retrieval algorithms to achieve operational retrieval of aerosol and ocean color information.
Recently, Fan et al. [28] have created a neural network to model the ocean body reflection matrix based on 3 parameters. In this work, we present an ocean reflectance model based on machine learning techniques, hereafter referred to as Neural Network Reflectance Prediction Model (NNRPM), which can be used to replace the radiative transfer simulation of ocean waters and expedite the retrieval algorithm significantly. This model incorporates additional parameters that would enable it to be better used for coastal waters. Moreover, it can also be used in semi-analytical models that retrieve IOPs from water leaving radiance or equivalent remote sensing reflectance. The basic advantages of neural networks are exploited, namely their accuracy performance properties, their robustness to noise, and their rapidity of execution.
This paper is organized as follows: Section 2 describes the methodology and definitions used in this work; Section 3 shows the result of this study; Section 4 summarizes the conclusion reached in this study.

2. Methods

Our goal is to build a polarized reflectance model for ocean waters. The ocean water model is the same as the one used in the MAPOL algorithm [27,29]. In order to model the optical properties of coastal waters the ocean water is assumed to consist of four components: pure sea water, phytoplankton particles, non-algae particles (NAP), and CDOM. The absorption coefficient a w ( λ ) and scattering coefficient b w ( λ ) for pure sea water are from the experimental data [30,31,32]. The backscattering fraction for water is 0.5 [33]. We denote the absorption coefficient for phytoplankton as a ph ( λ ) , the total absorption coefficient for both CDOM and NAP as a dg ( λ ) , the total backscattering coefficient and total backscattering fraction for both phytoplankton and NAP as b bp ( λ ) and B p ( λ ) . The scattering coefficient of CDOM is assumed to be negligible. The bio-optical model can be written as [29]:
a ph ( λ ) = A ph ( λ ) [ Chla ] E ph ( λ ) , a dg ( λ ) = a dg ( 440 ) exp [ S dg ( λ 440 ) ] , b bp ( λ ) = b bp ( 660 ) ( λ / 660 ) S bp , B p ( λ ) = B p ( 660 ) ( λ / 660 ) S Bp ,
where wavelength ( λ ) is in nm; a ph ( λ ) is determined by [Chla], Aph, and Eph; [Chla] is in units of mg/m3. Values of the coefficients Aph, and Eph are obtained from Bricaud et al. [34]. Sdg is the exponential spectral slope for adg; Sbp and SBp are the polynomial spectral slopes for bbp(λ) and Bp(λ) respectively. The total absorption coefficient is a(λ) = aw(λ) + aph(λ) + adg(λ), and the total backscattering coefficient is bb(λ) = 0.5bw(λ) + bbp(λ). The inherent optical properties (IOPs), aph(λ), adg(λ), bbp(λ), aw(λ), and bw(λ) are used to describe the ocean color spectrum in ocean bio-optical models [11,15,35]. Generalised IOP (GIOP) model is also used in the earth science community [35,36].
The Fournier-Forand (FF) [37] phase function is used to represent the phase functions of both phytoplankton and NAP, which is in agreement with various in-situ volume scattering function (VSF) measurements [38]. The scattering coefficient of the phytoplankton and NAP particles can be derived from b p = b bp / B p . The determination of the FF phase function [39] is done by backscattering fraction Bbp [39]. The total phase function for the coastal water system is obtained by mixing the FF phase function with the pure seawater phase function weighted by the scattering coefficients of the constituents. The normalized Mueller matrix is from measurements from average ocean waters [40,41,42]. The ocean water is assumed to be vertically homogeneous with a depth of 200 m. The minimum and maximum values of the parameters [29] considered are given in Table 1.
Here, we are focusing on the coastal waters so that the particle size is primarily large. In addition to the aforementioned parameters, viewing zenith angles range from 0 to 87.5° with an interval of 2.5°. The viewing azimuth angles are from 0 to 180° with an interval of 5°. The incident zenith angles are taken from 0 to 85°. The wavelengths used are 385, 410, 440, 470, 555, 670, 760, 863.5, and 870 nm. The aforementioned wavelengths include the ultraviolet, visible and near infrared (VNIR) part of the electromagnetic spectrum. The wavelengths 410, 470, 555, 670, and 863.5 nm are used as part of research scanning polarimeter (RSP) and the rest are adopted based on other existing and new ocean color sensors. The above mentioned bio-optical parameters, wavelengths, incident zenith angles (θi), viewing zenith and azimuth angles were provided as inputs to the radiative transfer code based on successive order of scattering for coupled atmosphere and ocean systems (SOS-CAOS) [43,44]. There is no atmosphere in our radiative transfer simulation as we are focused on the ocean reflectance model. Glint from the ocean surface can be predicted analytically, which is excluded in this simulation. The output of the radiative transfer simulations were considered as the true or desired value while training the neural network.
The objective of the neural network is to predict the non-dimensional effective Mueller matrix [45], M as a function of the wavelength λ , viewing zenith angle θ v , viewing azimuth angle ϕ v , incident zenith angle θ i , and incident azimuth angle ϕ i :
M ( λ , θ v , θ i , ϕ v ϕ i ) = M 11 M 12 M 13 M 14 M 21 M 22 M 23 M 24 M 31 M 32 M 33 M 34 M 41 M 42 M 43 M 44 .
We denote the radiance column vector as L = ( I , Q , U , V ) T , where ( I , Q , U , V ) are the Stokes parameters representing the polarization state of a light beam and the superscript T stands for transpose.
Once M is obtained, the upwelling radiance vector L ( θ v , ϕ v ) can be calculated from incident downwelling radiance vector L ( θ i , ϕ i ) . For an arbitrarily polarized incident light, the reflected radiance can be calculated from the following relation:
L ( θ v , ϕ v ) = 1 π M ( θ v , θ i , ϕ v ϕ i ) L ( θ i , ϕ i ) c o s θ i d c o s θ i d ϕ i ,
where we have implicitly assumed the wavelength dependency of both M and L. Hereafter, we neglect the fourth column and row of M, due to their small contribution for oceanic environment [46]. The 3 × 3 Muller matrix M can be obtained by considering the following incident polarizations: ( I 1 , Q 1 , U 1 ) as complete unpolarization (1,0,0), ( I 2 , Q 2 , U 2 ) as horizontal polarization (1,1,0) and ( I 3 , Q 3 , U 3 ) as 4 5 polarization (1,0,1) to get:
M = π F 0 c o s θ i I w 1 I w 2 I w 1 I w 3 I w 1 Q w 1 Q w 2 Q w 1 Q w 3 Q w 1 U w 1 U w 2 U w 1 U w 3 U w 1 ,
where ( I w 1 , Q w 1 , U w 1 ) , ( I w 2 , Q w 2 , U w 2 ) , ( I w 3 , Q w 3 , U w 3 ) represents remote sensing reflectance just above the ocean water surface corresponding to the three polarized incident beams. F 0 is the incident irradiance. The elements of the matrix in Equation (4), are obtained by solving SOS-CAOS [43,44]. As part of the neural network training, all the elements (except M11) of the matrix in Equation (4) are normalized by M11.
We have created a supervised, model based, batch learning, neural network algorithm with the architecture of a multilayer perceptron (MLP) using MATLAB’s Neural Network package [47]. An MLP is a feedforward and backpropagated artificial neural network [48]. This neural network approximates the function relating the various inputs of the network to the output [49] and would be used to predict ocean bidirectional reflectance distribution function in this work [50]. The network was trained with 10 inputs (Chla, a dg , b bp , B p , S Bp , S bp , S dg , θ i , θ v and ϕ v ) and the output reflectance for every wavelength being studied.
Figure 1 shows the structure of the neural network used in this work. We arrived at the above architecture by using the general guidelines in [51] as the base and then modifying it as per the requirements of our dataset. The input layer has 10 inputs x 1 to x 10 . The hidden layer consists of 4 layers each having 40 nodes. The output layer has one node. Each node has a summation operator and an activation function within it except the input nodes. The output of every node in a layer is multiplied by a weight and summed with the output of other weighted nodes in that layer and a bias to be fed as input to the activation function. The output of the activation function becomes the input of a node in the following layer. we have used tangent sigmoid as our activation function:
f = tansig ( x ) = 2 ( 1 + exp ( 2 * x ) ) 1 ,
where x denotes the sum of weighted nodes of the layer and the bias.
For each training instance this neural network algorithm first makes a prediction on the forward pass, measures the error, then goes through each layer in reverse to measure the error contribution from each previous connection (on the reverse pass), and finally tweaks each of the connection weights to reduce the error as part of Levenberg-Marquardt optimization algorithm [52]. So, essentially, we start off with random weight values and then improve it gradually one step at a time such that each step reduces the error until the algorithm converges to a minimum. The Levenberg-Marquardt optimization algorithm can be represented as:
w k + 1 = w k J T J + μ I 1 J T e ,
where w k is the weight/bias of previous iteration and wk+1 is the weight/bias of next iteration. J is the Jacobian matrix and JT is its transpose. J contains first derivatives of the network errors with respect to the weights and biases, and e is a vector of network errors. I is the identity matrix. The parameter μ determines direction as well as the length of the particular iteration step. It is a value that speeds up or slows down how quickly an algorithm learns. Technically, it determines the size of step an algorithm takes when moving towards a global minimum. Levenberg-Marquardt algorithm combines the Gauss-Newton and gradient descent optimizer methods. For larger values of μ, this algorithm will use gradient descent as the primary optimizer and for very small values of μ, Gauss-Newton would act as the primary optimizer. Here we have used μ = 10−4.
After each iteration of training, the performance of neural network was monitored on the test dataset by calculating the mean square error. In order to avoid over-fitting, the training was stopped when the mean square error continued to increase for 6 consecutive iterations on the test dataset [51]. The model equation can be represented as:
y = b 5 [ 1 ] + Σ k = 1 N h w 5 [ 1 , k ] . f b 4 [ k ] + Σ k = 1 N h w 4 [ k , k ] . f { b 3 [ k ] + Σ k = 1 N h w 3 [ k , k ] . f [ b 2 [ k ] + Σ k = 1 N h w 2 [ k , k ] . f ( b 1 [ k ] + Σ l = 1 N i w 1 [ k , l ] . x l ) ] } ,
where y is the output reflectance, predicted by the network. Ni is the number of input parameters, Nh is number of neurons in each hidden layer, w 1 [ k , l ] , w 2 [ k , k ] , w 3 [ k , k ] , w 4 [ k , k ] , and w 5 [ 1 , k ] are the weights of the four hidden layers and the output layer respectively, and b 1 [ k ] , b 2 [ k ] , b 3 [ k ] , b 4 [ k ] , and b 5 [ 1 ] are the biases of the four hidden layers and the output layer respectively. f is the tangent sigmoid activation function given in Equation (5). Figure 2 shows the high level flow diagram of the neural network.
Feature sequences for the various input parameters were generated using a quasi-random number sequence, the Halton sequence [53]. First introduced by Halton [54], it is constructed using a deterministic method and is of low discrepancy. Training and test datasets corresponding to the above sequences were generated by SOS-CAOS [43,44]. The neural network was trained and validated using the training and test datasets.
The error measurement for the element M11 is computed using the relation below:
Percentage   Error   ( PE ) = 100 × | Ture   Value Predicted   Value | Ture   Value .
In order to avoid the numerical difficulty of dividing by zero, the error measurement for all other normalized elements (Mij/M11) is computed using the relation below:
Absolute   Error   ( AE ) = | Ture   Value Predicted   Value | .
The above errors are used to calculate the empirical cumulative distribution function, which gives the probability of the number of samples within a certain (error) value, which is used to illustrate the performance of the neural network. If { X 1 , , X n } is an independently and identically distributed sample from an unknown distribution function, F ( x ) = P ( X x ) , then the empirical cumulative distribution function ( F n ^ ( x ) ⇒ ECDF) is defined as [55]:
ECDF = 1 n i = 1 n 1 { X i x } ,
where,
1 { X i x } = 1 , if X i x 0 , otherwise .
Another metric used here to evaluate the performance of NNRPM is the Pearson correlation coefficient (R) [56]:
R = i = 1 n ( x i x ¯ ) ( y i y ¯ ) i = 1 n ( x i x ¯ ) 2 i = 1 n ( y i y ¯ ) 2 ,
where n is the sample size, xi and yi are sample points, x ¯ and y ¯ are sample means of two variables x and y respectively.

3. Results

In Figure 3, we present a plot between error percentage and Halton sequences. Here, error percentage represents the percentage of test cases with absolute error greater than 0.004. At 8 K Halton sequences, the error percentage reached its minimum value and then varies little after that. This established the following relationship with error percentage and training Halton combinations (N):
Error % 1 N .
In Figure 4, we have shown the performance of the neural network NNRPM by the linear regression plot between the true value calculated by SOS-CAOS and the predicted value predicted by NNRPM for elements M 22 / M 11 and M 32 / M 11 at wavelength 555 nm, for 1000 Halton sequence independently generated for validation.
Figure 5 and Figure 6 show the polar plots for the true value (SOS-CAOS), predicted value (NNRPM), and percentage error for the first element of the first column (M11) for two different sets of input parameters for wavelengths 385 and 555 nm. Figure 7 and Figure 8 show the polar plots for the true value (SOS-CAOS), predicted value (NNRPM), and absolute error for the reduced first element of the second column (M12/M11) for the same cases as shown in Figure 5 and Figure 6. For each polar plot, the zenith point is at the center, and the viewing zenith angle of 67.5° is at the rim of the polar plot. The viewing azimuthal angles are shown around the plots.
In Figure 9, we have shown ECDF versus PE at wavelengths 385 and 555 nm, which shows that ECDF approaches 1 asymptotically as PE increases and becomes very close to 1 when PE is larger than 5%.
In Table 2, Table 3 and Table 4, the performance of the neural networks corresponding to each element of the Mueller matrix from Equation (4), for every wavelength being studied, has been shown. For evaluating the prediction performance of the degree of linear polarization (DoLP = Q 2 + U 2 / I ), we tested the AE for the expression M 21 2 + M 31 2 / M 11 , which is equal to DoLP when the incident light is unpolarized. The results are shown in Table 2. For the first element M11, the NNRPM is able to predict values with percentage error (PE) less than 5% for 99% of the test cases. For the rest of the normalized elements, the NNRPM is able to predict values with absolute error (AE) less than 0.004 for 99% of the test cases. DoLP for unpolarized incident light can be predicted with an absolute error (AE) less than 0.002 for 99% of the test cases.
The CPU time taken by the SOS-CAOS RT code for the prediction of ocean BRDF is 3 min for a case simulated with the following parameters:
  • Order of scattering: 40
  • Gaussian quadrature for ocean: 80
  • Optical depth of the system (for wavelength 555 nm): 6986.19
  • Scattering coefficient of phytoplankton and NAP particles, bp: 34.468 m−1
  • System specifications: 2.8 GHz Intel Core i7 and memory of 16 GB 1600 MHz DDR3
The corresponding time taken by NNRPM is 0.0296 seconds which is mostly stable for an arbitrary case.

4. Discussion

In order to determine the optimum number of training cases required, we experimented with various sets of Halton sequences. We found that accuracy increases as the number of Halton combinations used in the training increases until it reaches 8000. From Figure 3, we can see that error initially decreases as the training data is increased and reaches a minimum for a certain volume of training data. This result is similar to the error’s inverse relationship with number of training instances mentioned in Gallant et al. [57]. Figure 4 shows a linear relationship between the true value represented by SOS-CAOS and the predicted value represented by NNRPM. The value of R-squared being close to 1 indicates clear correlation between the true value (SOS-CAOS) and the predicted value (NNRPM). The slope (m) being close to 1 and bias (y) being close to 0 further confirms the agreement between NNRPM with SOS-CAOS.
The Bidirectional Reflectance Distribution Function (BRDF) is the quantity which geometrically characterizes the reflecting properties of a point on a surface [58]. The angular distribution of the water leaving signals are illustrated through the polar plots in Figure 5, Figure 6, Figure 7 and Figure 8. The SOS-CAOS and NNRPM polar plots provide a pictorial representation of how close the neural network predicted value is to the actual SOS-CAOS value for the specific scenario under consideration. In Figure 5 and Figure 6, the accuracy of the model for the element M11 is illustrated by the Percentage error polar plots which has a maximum value around 3%. For test cases with the smaller values of adg, the values of reflectance will be higher for both wavelengths, and vice versa. In Figure 7 and Figure 8, the accuracy of the model for the element M12/M11 is illustrated by the Absolute error polar plots which has a maximum value around 0.0009. We can see that the model is able to predict both positive and negative results equally well.
Figure 9 is the visual representation of how quickly the cumulative distribution function increases to 1. It shows the probability of an unknown test dataset having a percentage error below any chosen threshold. Here we can say that the probability of getting a percentage error of less than 5% across the entire test dataset is close to 1. This further confirms the accuracy of the model on an unseen test dataset.
The simulation time of SOS-CAOS depends on the order of scattering and the number of Gaussian quadratures used for angular integration, which is in turn determined by the required accuracy. Moreover, it also depends upon the total optical depth of the system. Using the parameters mentioned in the results section, we can see that the incorporation of NNRPM would improve the speed of retrieval algorithms by a factor of 6000. For different ocean reflectance cases, the CPU time varies for the SOS-CAOS solution. However the order of magnitude is about the same. Thus the comparison provides a general guideline on the improvement in CPU time.
The future incorporation of this model into the retrieval algorithms will make the retrieval process more efficient, and thus applicable for operational use with global satellite observations. Future work could also examine the necessity of adding more IOPs to these models to simulate the behavior of different water bodies.

5. Conclusions

Ocean polarized reflectances obtained from satellite sensors contains information about the constituents of the upper ocean euphotic zone, such as colored dissolved organic matter (CDOM), sediments, phytoplankton types, and pollutants. Traditional ocean color remote sensing algorithms retrieve the constituent information by first performing atmospheric correction, then deriving IOPs from the resultant remote sensing reflectances. However, atmospheric correction algorithms are still a challenging task for scenes involving coastal waters and absorbing aerosols. In order to overcome this difficulty, joint retrieval algorithms of aerosols and ocean color have been designed to use information in polarized measurements. Joint retrieval algorithms need to solve the radiative transfer equation in atmospheric and ocean coupled system repeatedly, which is resource prohibitive, making them impractical for operational production of satellite data.
In order to resolve this problem, we have built a forward model using machine learning techniques, called Neural Network Reflectance Prediction Model (NNRPM), which can be used to predict ocean bidirectional reflectance, given inherent optical properties of ocean waters. The accuracy of the model is verified with the test cases generated from radiative transfer code based on the successive order of scattering technique (SOS-CAOS). The NNRPM model is able to predict the reflectance values for the first element with percentage error less than 5% for 99% of the test cases. For the rest of the normalized elements, the NNRPM is able to predict reflectance values with an absolute error less than 0.004 for 99% of the test cases. For unpolarized incident light, the model is able to predict DoLP with an AE less than 0.002 for 99% of test cases. When comparing the average speeds of the NNRPM against the SOS-CAOS model, we see that the NNRPM is able to predict the Ocean BRDF 6000 times faster than SOS-CAOS. The application of this neural network model in inverse algorithms would expedite the retrieval process and make them more efficient for use in global satellite observations. It can also be used in traditional semi-analytical models of retrieving ocean water IOPs, since our bio-optical model is more general than most existing models.

Author Contributions

Conceptualization, L.M., P.-W.Z. and Y.H.; methodology, L.M., P.-W.Z. and M.G.; software, L.M.; validation, L.M.; formal analysis, L.M.; investigation, L.M. and P.-W.Z.; resources, P.-W.Z.; data curation, L.M.; writing—original draft preparation, L.M.; writing—review and editing, L.M., P.-W.Z., M.G., Y.H., B.A.F. and P.J.W.; visualization, L.M. and P.-W.Z.; supervision, P.-W.Z.; project administration, P.-W.Z.; funding acquisition, L.M. and P.-W.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by NASA Headquarters under the NASA Earth and Space Science Fellowship Program - Grant 80NSSC17K0366.

Acknowledgments

The hardware used in the computational studies is part of the UMBC High Performance Computing Facility (HPCF). The facility is supported by the U.S. National Science Foundation through the MRI program (grant nos. CNS-0821258 and CNS-1228778) and the SCREMS program (grant no. DMS-0821311), with additional substantial support from the University of Maryland, Baltimore County (UMBC).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Antoine, D.; Morel, A. A multiple scattering algorithm for atmospheric correction of remotely sensed ocean colour (MERIS instrument): Principle and implementation for atmospheres carrying various aerosols including absorbing ones. Int. J. Remote Sens. 1999, 20, 1875–1916. [Google Scholar] [CrossRef]
  2. Fukushima, H.; Higurashi, A.; Mitomi, Y.; Nakajima, T.; Noguchi, T.; Tanaka, T.; Toratani, M. Correction of atmospheric effect on ADEOS/OCTS ocean color data: Algorithm description and evaluation of its performance. J. Oceanogr. 1998, 54, 417–430. [Google Scholar] [CrossRef]
  3. Gordon, H.R.; Wang, M. Retrieval of water-leaving radiance and aerosol optical thickness over the oceans with SeaWiFS: A preliminary algorithm. Appl. Opt. 1994, 33, 443–452. [Google Scholar] [CrossRef] [PubMed]
  4. Dierssen, H.M. Perspectives on empirical approaches for ocean color remote sensing of chlorophyll in a changing climate. Proc. Natl. Acad. Sci. USA 2010, 107, 17073–17078. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Hu, C.; Lee, Z.; Franz, B. Chlorophyll aalgorithms for oligotrophic oceans: A novel approach based on three band reflectance difference. J. Geophys. Res. Oceans 2012, 117, 1–25. [Google Scholar] [CrossRef] [Green Version]
  6. Shanmugam, P.; Ahn, Y.H.; Ryu, J.H.; Sundarabalan, B. An evaluation of inversion models for retrieval of inherent optical properties from ocean color in coastal and open sea waters around Korea. J. Oceanogr. 2010, 66, 815–830. [Google Scholar] [CrossRef]
  7. Boss, E.; Roesler, C. Over constrained linear matrix inversion with statistical selection. In Remote Sensing of Inherent Optical Properties: Fundamentals, Tests of Algorithms, and Applications; Volume 5 of Reports of the International Ocean Colour Coordinating Group; Lee, Z.P., Ed.; IOCCG: Dartmouth, NS, Canada, 2006; pp. 57–62. [Google Scholar]
  8. Garver, S.A.; Siegel, D.A. Inherent optical property inversion of ocean color spectra and its biogeochemical interpretation: 1. Time series from the Sargasso Sea. J. Geophys. Res. Oceans 1997, 102, 18607–18625. [Google Scholar] [CrossRef]
  9. Maritorena, S.; Siegel, D.A.; Peterson, A.R. Optimization of a semianalytical ocean color model for global-scale applications. Appl. Opt. 2002, 41, 2705–2714. [Google Scholar] [CrossRef]
  10. Gordon, H.R.; Brown, O.B.; Evans, R.H.; Brown, J.W.; Smith, R.C.; Baker, K.S.; Clark, D.K. A semianalytic radiance model of ocean color. J. Geophys. Res. Atmos. 1988, 93, 10909–10924. [Google Scholar] [CrossRef]
  11. IOCCG. Remote Sensing of Inherent Optical Properties: Fundamentals, Tests of Algorithms, and Applications. In Reports of the International Ocean Colour Coordinating Group; Lee, Z.P., Ed.; IOCCG: Dartmouth, NS, Canada, 2006; Volume 5. [Google Scholar]
  12. Ligi, M.; Kutser, T.; Kallio, K.; Attila, J.; Koponen, S.; Paavel, B.; Soomets, T.; Reinart, A. Testing the performance of empirical remote sensing algorithms in the Baltic Sea waters with modelled and in situ reflectance data. Oceanologia 2017, 59, 57–68. [Google Scholar] [CrossRef] [Green Version]
  13. Morel, A.; Gentili, B. A simple band ratio technique to quantify the colored dissolved and detrital organic material from ocean color remotely sensed data. Remote Sens. Environ. 2009, 113, 998–1011. [Google Scholar] [CrossRef]
  14. Franz, B.A.; Werdell, P.J. A generalized framework for modeling of inherent optical properties in ocean remote sensing applications. Proc. Ocean Opt. Anchorage Alaska 2010, 27, 1–3. [Google Scholar]
  15. Werdell, P.J.; McKinna, L.I.; Boss, E.; Ackleson, S.G.; Craig, S.E.; Gregg, W.W.; Lee, Z.; Maritorena, S.; Roesler, C.S.; Rousseaux, C.S.; et al. An overview of approaches and challenges for retrieving marine inherent optical properties from ocean color remote sensing. Prog. Oceanogr. 2018, 160, 186–212. [Google Scholar] [CrossRef] [PubMed]
  16. IOCCG. Remote Sensing of Ocean Colour in Coastal, and Other Optically-Complex, Waters. In Reports of the International Ocean Colour Coordinating Group; Sathyendranath, S., Ed.; IOCCG: Dartmouth, NS, Canada, 2000; Volume 3. [Google Scholar]
  17. IOCCG. Atmospheric Correction for Remotely-Sensed Ocean-Colour Products. In Reports of the International Ocean Colour Coordinating Group; Wang, M., Ed.; IOCCG: Dartmouth, NS, Canada, 2010; Volume 10. [Google Scholar]
  18. Deschamps, P.Y.; Bréon, F.M.; Leroy, M.; Podaire, A.; Bricaud, A.; Buriez, J.C.; Seze, G. The POLDER mission: Instrument characteristics and scientific objectives. IEEE Trans. Geosci. Remote Sens. 1994, 32, 598–615. [Google Scholar] [CrossRef]
  19. Cairns, B.; Russell, E.E.; Travis, L.D. Research scanning polarimeter: Calibration and ground-based measurements. In Polarization: Measurement, Analysis, and Remote Sensing II. Int. Soc. Opt. Photonics 1999, 3754, 186–196. [Google Scholar]
  20. Diner, D.J.; Xu, F.; Garay, M.J.; Martonchik, J.V.; Rheingans, B.E.; Geier, S.; Davis, A.; Hancock, B.R.; Jovanovic, V.M.; Bull, M.A.; et al. the Airborne Multiangle SpectroPolarimetric Imager (AirMSPI): A new tool for aerosol and cloud remote sensing. Atmos. Meas. Tech. 2013, 6, 2007–2025. [Google Scholar] [CrossRef] [Green Version]
  21. Snik, F.; Rietjens, J.H.; Van Harten, G.; Stam, D.M.; Keller, C.U.; Smit, J.M.; Laan, E.C.; Verlaan, A.L.; Ter Horst, R.; Navarro, R.; et al. SPEX: The spectropolarimeter for planetary exploration. Space Telescopes and Instrumentation 2010: Optical, Infrared, and Millimeter Wave. Int. Soc. Opt. Photonics 2010, 7731. [Google Scholar] [CrossRef] [Green Version]
  22. Martins, J.V.; Fernandez-Borda, R.; McBride, B.; Remer, L.; Barbosa, H.M. The Harp Hype Ran Gular Imaging Polarimeter and the Need for Small Satellite Payloads with High Science Payoff for Earth Science Remote Sensing. In Proceedings of the IGARSS 2018— IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 22–27 July 2018; pp. 6304–6307. [Google Scholar]
  23. Werdell, P.J.; Behrenfeld, M.J.; Bontempi, P.S.; Boss, E.; Cairns, B.; Davis, G.T.; Franz, B.A.; Gliese, U.B.; Gorman, E.T.; Hasekamp, O.; et al. the Plankton, Aerosol, Cloud, ocean Ecosystem mission: Status, science, advances. Bull. Am. Meteorol. Soc. 2019, 100, 1775–1794. [Google Scholar] [CrossRef]
  24. Chowdhary, J.; Cairns, B.; Mishchenko, M.I.; Hobbs, P.V.; Cota, G.F.; Redemann, J.; Rutledge, K.; Holben, B.N.; Russell, E. Retrieval of aerosol scattering and absorption properties from photopolarimetric observations over the ocean during the CLAMS experiment. J. Atmos. Sci. 2005, 62, 1093–1117. [Google Scholar] [CrossRef]
  25. Hasekamp, O.P.; Litvinov, P.; Butz, A. Aerosol properties over the ocean from PARASOL multiangle photopolarimetric measurements. J. Geophys. Res. Atmos. 2011, 116, 2156–2202. [Google Scholar] [CrossRef] [Green Version]
  26. Xu, F.; Dubovik, O.; Zhai, P.W.; Diner, D.J.; Kalashnikova, O.V.; Seidel, F.C.; Litvinov, P.; Bovchaliuk, A.; Garay, M.J.; van Harten, G.; et al. Joint retrieval of aerosol and water-leaving radiance from multispectral, multiangular and polarimetric measurements over ocean. Atmos. Meas. Tech. 2016, 9, 2877–2907. [Google Scholar] [CrossRef] [Green Version]
  27. Gao, M.; Zhai, P.W.; Franz, B.; Hu, Y.; Knobelspiesse, K.; Werdell, P.J.; Ibrahim, A.; Cairns, B.; Chase, A. Inversion of multi-angular polarimetric measurements over open and coastal ocean waters: A joint retrieval algorithm for aerosol and water leaving radiance properties. Atmos. Meas. Tech. 2019, 12, 3921–3941. [Google Scholar] [CrossRef] [Green Version]
  28. Fan, C.; Fu, G.; Di Noia, A.; Smit, M.; Rietjens, H.H.J.; Ferrare, A.R.; Burton, S.; Li, Z.; Hasekamp, P.O. Use of a Neural Network-Based Ocean Body Radiative Transfer Model for Aerosol Retrievals from Multi-Angle Polarimetric Measurements. Remote Sens. 2019, 11, 2877. [Google Scholar] [CrossRef] [Green Version]
  29. Gao, M.; Zhai, P.W.; Franz, B.; Hu, Y.; Knobelspiesse, K.; Werdell, P.J.; Ibrahim, A.; Xu, F.; Cairns, B. Retrieval of aerosol properties and water-leaving reflectance from multi-angular polarimetric measurements over coastal waters. Opt. Express 2018, 26, 8968–8989. [Google Scholar] [CrossRef] [PubMed]
  30. Kou, L.; Labrie, D.; Chylek, P. Refractive indices of water and ice in the 0.65-to 2.5-μm spectral range. Appl. Opt. 1993, 32, 3531–3540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Pope, R.M.; Fry, E.S. Absorption spectrum (380–700 nm) of pure water. II. Integrating cavity measurements. Appl. Opt. 1997, 36, 8710–8723. [Google Scholar] [CrossRef]
  32. Zhang, X.; Hu, L. Scattering by pure seawater at high salinity. Opt. Express 2009, 17, 12685–12691. [Google Scholar] [CrossRef]
  33. Morel, A. Optical properties of pure water and pure sea water. Opt. Asp. Oceanogr. 1974, 1, 1–24. [Google Scholar]
  34. Bricaud, A.; Morel, A.; Babin, M.; Allali, K.; Claustre, H. Variations of light absorption by suspended particles with chlorophyll a concentration in oceanic (case 1) waters: Analysis and implications for bio-optical models. J. Geophys. Res. Oceans 1998, 103, 31033–31044. [Google Scholar] [CrossRef]
  35. Werdell, P.J.; Franz, B.A.; Bailey, S.W.; Feldman, G.C.; Claustre, H. Generalized ocean color inversion model for retrieving marine inherent optical properties. Appl. Opt. 2013, 52, 2019–2037. [Google Scholar] [CrossRef]
  36. Werdell, P.J.; McKinna, L.I.W. Sensitivity of Inherent Optical Properties From Ocean Reflectance Inversion Models to Satellite Instrument Wavelength Suites. Front. Earth Sci. 2019, 7, 54. [Google Scholar] [CrossRef] [PubMed]
  37. Fournier, G.R.; Forand, J.L. Analytic phase function for ocean water. In Ocean Optics XII; International Society for Optics and Photonics: Bergen, Norway, 1994; pp. 194–201. [Google Scholar]
  38. Sullivan, J.M.; Twardowski, M.S. Angular shape of the oceanic particulate volume scattering function in the backward direction. Appl. Opt. 2009, 48, 6811–6819. [Google Scholar] [CrossRef] [PubMed]
  39. Mobley, C.D.; Sundman, L.K. Phase function effects on oceanic light fields. Appl. Opt. 2002, 41, 1035–1050. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Voss, K.J.; Fry, E.S. Measurement of the Mueller matrix for ocean water. Appl. Opt. 1984, 23, 4427–4439. [Google Scholar] [CrossRef]
  41. Kokhanovsky, A.A. Parameterization of the Mueller matrix of oceanic waters. J. Geophys. Res. Oceans 2003, 108, 3175. [Google Scholar] [CrossRef]
  42. Zhai, P.W.; Knobelspiesse, K.; Ibrahim, A.; Franz, B.A. Water-leaving contribution to polarized radiation field over ocean. Opt. Express 2017, 25, A689–A708. [Google Scholar] [CrossRef]
  43. Zhai, P.W.; Hu, Y.; Trepte, C.R.; Lucker, P.L. A vector radiative transfer model for coupled atmosphere and ocean systems based on successive order of scattering method. Opt. Express 2009, 17, 2057–2079. [Google Scholar] [CrossRef]
  44. Zhai, P.W.; Yongxiang, H.; Chowdhary, J.; Trepte, C.R.; Lucker, P.L.; Josset, D.B. A vector radiative transfer model for coupled atmosphere and ocean systems with a rough interface. J. Quant. Spectrosc. Radiat. Transf. 2010, 111, 1025–1040. [Google Scholar] [CrossRef]
  45. Lawless, R.; Xie, Y.; Yang, P.; Kattawar, G.W.; Laszlo, I. Polarization and effective Mueller matrix for multiple scattering of light by nonspherical ice crystals. Opt. Express 2006, 14, 6381–6393. [Google Scholar] [CrossRef]
  46. Kattawar, G.W.; Adams, C.N. Stokes vector calculations of the submarine light field in an atmosphere-ocean with scattering according to a Rayleigh phase matrix: Effect of interface refractive index on radiance and polarization. Limnol. Oceanogr. 1989, 34, 1453–1472. [Google Scholar] [CrossRef]
  47. MathWorks. Deep Learning Toolbox; The MathWorks Inc.: Natick, MA, USA, 2019. [Google Scholar]
  48. Smith, J.A. LAI inversion using a back-propagation neural network trained with a multiple scattering model. IEEE Trans. Geosci. Remote Sens. 1993, 31, 1102–1106. [Google Scholar] [CrossRef] [Green Version]
  49. Smith, J.A. Multilayer feedforward networks are universal approximators. Neural Netw. 1989, 2, 359–366. [Google Scholar]
  50. Abuelgasim, A.A.; Gopal, S.; Strahler, A.H. Forward and inverse modelling of canopy directional reflectance using a neural network. Int. J. Remote Sens. 1998, 19, 453–471. [Google Scholar] [CrossRef]
  51. Gron, A. Hands-On Machine Learning with Scikit-Learn and TensorFlow: Concepts, Tools, and Techniques to Build Intelligent Systems, 1st ed.; O’Reilly Media, Inc.: Newton, MA, USA, 2017. [Google Scholar]
  52. Marquardt, D.W. An algorithm for least-squares estimation of nonlinear parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  53. Loyola, R.; Diego, G.; Pedergnana, M.; Gimeno, S.G. Smart sampling and incremental function learning for very large high dimensional data. Neural Netw. 2016, 78, 75–87. [Google Scholar] [CrossRef] [Green Version]
  54. Halton, J.H. On the efficiency of certain quasi-random sequences of points in evaluating multi-dimensional integrals. Numer. Math. 1960, 2, 84–90. [Google Scholar] [CrossRef]
  55. Pernot, P.; Savin, A. Probabilistic performance estimators for computational chemistry methods: The empirical cumulative distribution function of absolute errors. J. Chem. Phys. 2018, 148, 241707. [Google Scholar] [CrossRef] [Green Version]
  56. James, G.; Witten, D.; Hastie, T.; Tibshirani, R. An Introduction to Statistical Learning: With Applications in R; Springer Publishing Company: Berlin, Germany, 2014. [Google Scholar]
  57. Gallant, A.R.; White, H. There exists a neural network that does not make avoidable mistakes. In Proceedings of the International Conference on Neural Networks, San Diego, CA, USA, 24–27 July 1988. [Google Scholar]
  58. Nicodemus, F.E.; Richmond, J.C.; Hsia, J.J.; Ginsberg, I.W.; Limperis, T. Geometrical Considerations and Nomenclature for Reflectance; NBS Monograph; National Bureau of Standards: Washington, DC, USA, 1992. [Google Scholar]
Figure 1. Structure of the neural network.
Figure 1. Structure of the neural network.
Remotesensing 12 01421 g001
Figure 2. High level flow diagram for NNRPM.
Figure 2. High level flow diagram for NNRPM.
Remotesensing 12 01421 g002
Figure 3. Error% versus Halton sequences for reduced M22/M11 element at 410 nm.
Figure 3. Error% versus Halton sequences for reduced M22/M11 element at 410 nm.
Remotesensing 12 01421 g003
Figure 4. Linear regression between the true reflectance values generated by SOS-CAOS and those predicted by NNRPM for elements (a) M22/M11 and (b) M32/M11 at 555 nm for 1000 Halton sequences.
Figure 4. Linear regression between the true reflectance values generated by SOS-CAOS and those predicted by NNRPM for elements (a) M22/M11 and (b) M32/M11 at 555 nm for 1000 Halton sequences.
Remotesensing 12 01421 g004
Figure 5. Polar plots of SOS-CAOS, NNRPM and PE for M11 at wavelength 385 nm. Input parameters for Top Row: θi = 55.91°, Chla = 20.43 mg/m3, adg(440) = 7.36 m−1, bbp = 0.097 m−1, Bp = 0.015, SBp = = −0.157 nm−1, Sbp = 0.258 nm−1,Sdg = 0.019 nm−1. Input parameters for Bottom Row: θi = 43.02°, Chla = 24.77 mg/m3, adg(440) = 0.004 m−1, bbp = 0.099 m−1, Bp = 0.046, SBp = −0.106 nm−1, Sbp = 0.023 nm−1, Sdg = 0.014 nm−1.
Figure 5. Polar plots of SOS-CAOS, NNRPM and PE for M11 at wavelength 385 nm. Input parameters for Top Row: θi = 55.91°, Chla = 20.43 mg/m3, adg(440) = 7.36 m−1, bbp = 0.097 m−1, Bp = 0.015, SBp = = −0.157 nm−1, Sbp = 0.258 nm−1,Sdg = 0.019 nm−1. Input parameters for Bottom Row: θi = 43.02°, Chla = 24.77 mg/m3, adg(440) = 0.004 m−1, bbp = 0.099 m−1, Bp = 0.046, SBp = −0.106 nm−1, Sbp = 0.023 nm−1, Sdg = 0.014 nm−1.
Remotesensing 12 01421 g005
Figure 6. Polar plots of SOS-CAOS, NNRPM and PE for M11 at wavelength 555 nm. Input parameters for Top Row: θi = 55.91°, Chla = 20.43 mg/m3, adg(440) = 0.258 m−1, bbp = 0.088 m−1, Bp = 0.016, SBp = −0.157 nm−1, Sbp = 0.258 nm−1, Sdg = 0.019 nm−1. Input parameters for Bottom Row: θi = 43.02°, Chla = 24.77 mg/m3, adg(440) = 0.0004 m−1, bbp = 0.098 m−1, Bp = 0.048, SBp = −0.106 nm−1, Sbp = 0.023 nm−1, Sdg = 0.014 nm−1.
Figure 6. Polar plots of SOS-CAOS, NNRPM and PE for M11 at wavelength 555 nm. Input parameters for Top Row: θi = 55.91°, Chla = 20.43 mg/m3, adg(440) = 0.258 m−1, bbp = 0.088 m−1, Bp = 0.016, SBp = −0.157 nm−1, Sbp = 0.258 nm−1, Sdg = 0.019 nm−1. Input parameters for Bottom Row: θi = 43.02°, Chla = 24.77 mg/m3, adg(440) = 0.0004 m−1, bbp = 0.098 m−1, Bp = 0.048, SBp = −0.106 nm−1, Sbp = 0.023 nm−1, Sdg = 0.014 nm−1.
Remotesensing 12 01421 g006
Figure 7. Polar plots of SOS-CAOS, NNRPM and AE for M12/M11 at wavelength 385 nm. Input parameters for Top Row: θi = 55.91°, Chla = 20.43 mg/m3, adg(440) = 7.36 m−1, bbp = 0.097 m−1, Bp = 0.015, SBp = −0.157 nm−1, Sbp = 0.258 nm−1, Sdg = 0.019 nm−1. Input parameters for Bottom Row: θi = 43.02°, Chla = 24.77 mg/m3, adg(440) = 0.004 m−1, bbp = 0.099 m−1, Bp = 0.046, SBp = −0.106 nm−1, Sbp = 0.023 nm−1, Sdg = 0.014 nm−1.
Figure 7. Polar plots of SOS-CAOS, NNRPM and AE for M12/M11 at wavelength 385 nm. Input parameters for Top Row: θi = 55.91°, Chla = 20.43 mg/m3, adg(440) = 7.36 m−1, bbp = 0.097 m−1, Bp = 0.015, SBp = −0.157 nm−1, Sbp = 0.258 nm−1, Sdg = 0.019 nm−1. Input parameters for Bottom Row: θi = 43.02°, Chla = 24.77 mg/m3, adg(440) = 0.004 m−1, bbp = 0.099 m−1, Bp = 0.046, SBp = −0.106 nm−1, Sbp = 0.023 nm−1, Sdg = 0.014 nm−1.
Remotesensing 12 01421 g007
Figure 8. Polar plots of SOS-CAOS, NNRPM and AE for M12/M11 at wavelength 555 nm. Input parameters for Top Row: θi = 55.91°, Chla = 20.43 mg/m3, adg(440) = 0.258 m−1, bbp = 0.088 m−1, Bp = 0.016, SBp = −0.157 nm−1, Sbp = 0.258 nm−1, Sdg = 0.019 nm−1. Input parameters for Bottom Row: θi = 43.02°, Chla = 24.77 mg/m3, adg(440) = 0.0004 m−1, bbp = 0.098 m−1, Bp = 0.048, SBp = −0.106 nm−1, Sbp = 0.023 nm−1, Sdg = 0.014 nm−1.
Figure 8. Polar plots of SOS-CAOS, NNRPM and AE for M12/M11 at wavelength 555 nm. Input parameters for Top Row: θi = 55.91°, Chla = 20.43 mg/m3, adg(440) = 0.258 m−1, bbp = 0.088 m−1, Bp = 0.016, SBp = −0.157 nm−1, Sbp = 0.258 nm−1, Sdg = 0.019 nm−1. Input parameters for Bottom Row: θi = 43.02°, Chla = 24.77 mg/m3, adg(440) = 0.0004 m−1, bbp = 0.098 m−1, Bp = 0.048, SBp = −0.106 nm−1, Sbp = 0.023 nm−1, Sdg = 0.014 nm−1.
Remotesensing 12 01421 g008
Figure 9. ECDF vs PE for M11 at wavelengths (a) 385 and (b) 555 nm.
Figure 9. ECDF vs PE for M11 at wavelengths (a) 385 and (b) 555 nm.
Remotesensing 12 01421 g009
Table 1. Minimum and maximum values of parameters in the bio-optical models mentioned in Equation (1).
Table 1. Minimum and maximum values of parameters in the bio-optical models mentioned in Equation (1).
ParameterChlaadg(440)bbp(660)Bp(660)SdgSbpSBp
Unitmg/m3m−1m−1 nm−1nm−1nm−1
Min∼0.0∼0.0∼0.0∼0.00.010.0−0.2
Max30.02.50.10.050.020.50.2
Table 2. Percentage of cases for element M11 with Percentage error (PE) less than 5% and percentage of cases for elements M21/M11 and M31/M11 with Absolute error (AE) less than 0.004, and Degree Of Linear Polarization.
Table 2. Percentage of cases for element M11 with Percentage error (PE) less than 5% and percentage of cases for elements M21/M11 and M31/M11 with Absolute error (AE) less than 0.004, and Degree Of Linear Polarization.
M11M21/M11M31/M11DoLP
Wavelength (nm)PE < 5%AE < 0.004AE < 0.004AE < 0.002
38599.6499.7899.8399.23
41099.2599.6899.8298.90
44099.3999.7599.7398.97
47099.2499.7799.7098.98
55599.2999.7899.7198.79
67099.2999.7599.8099.02
863.599.5799.9399.9299.71
87099.5599.9499.9299.72
Table 3. Percentage of cases for elements M12/M11, M22/M11, and M32/M11 with Absolute error (AE) less than 0.004.
Table 3. Percentage of cases for elements M12/M11, M22/M11, and M32/M11 with Absolute error (AE) less than 0.004.
M12/M11M22/M11M32/M11
Wavelength (nm)AE < 0.004AE < 0.004AE < 0.004
38599.6799.4399.56
41099.3899.6499.40
44099.5099.3099.37
47099.2699.1399.08
55598.6899.5499.24
67099.5399.7099.667
863.599.9799.6899.79
87099.9299.7199.72
Table 4. Percentage of cases for elements M13/M11, M23/M11, and M33/M11 with Absolute error (AE) of elements less than 0.004.
Table 4. Percentage of cases for elements M13/M11, M23/M11, and M33/M11 with Absolute error (AE) of elements less than 0.004.
M13/M11M23/M11M33/M11
Wavelength (nm)AE < 0.004AE < 0.004AE < 0.004
38598.8999.3499.42
41098.8399.2399.39
44098.7899.0599.29
47098.0098.8699.17
55598.2098.3199.32
67098.5699.5299.63
863.599.9599.7399.74
87099.9299.6999.74

Share and Cite

MDPI and ACS Style

Mukherjee, L.; Zhai, P.-W.; Gao, M.; Hu, Y.; A. Franz, B.; Werdell, P.J. Neural Network Reflectance Prediction Model for Both Open Ocean and Coastal Waters. Remote Sens. 2020, 12, 1421. https://doi.org/10.3390/rs12091421

AMA Style

Mukherjee L, Zhai P-W, Gao M, Hu Y, A. Franz B, Werdell PJ. Neural Network Reflectance Prediction Model for Both Open Ocean and Coastal Waters. Remote Sensing. 2020; 12(9):1421. https://doi.org/10.3390/rs12091421

Chicago/Turabian Style

Mukherjee, Lipi, Peng-Wang Zhai, Meng Gao, Yongxiang Hu, Bryan A. Franz, and P. Jeremy Werdell. 2020. "Neural Network Reflectance Prediction Model for Both Open Ocean and Coastal Waters" Remote Sensing 12, no. 9: 1421. https://doi.org/10.3390/rs12091421

APA Style

Mukherjee, L., Zhai, P. -W., Gao, M., Hu, Y., A. Franz, B., & Werdell, P. J. (2020). Neural Network Reflectance Prediction Model for Both Open Ocean and Coastal Waters. Remote Sensing, 12(9), 1421. https://doi.org/10.3390/rs12091421

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