Next Article in Journal
Airborne S-Band SAR for Forest Biophysical Retrieval in Temperate Mixed Forests of the UK
Previous Article in Journal
Monitoring Urban Areas with Sentinel-2A Data: Application to the Update of the Copernicus High Resolution Layer Imperviousness Degree
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Successive Cancellation Method to Retrieve Sea Wave Components from Spatio-Temporal Remote Sensing Image Sequences

1
College of Automation, Harbin Engineering University, No. 145 Nantong Street, Harbin 150001, China
2
Department of Electrical and Computer Engineering, McMaster University, 1280 Main Street West, Hamilton, ON L8S 4K1, Canada
*
Author to whom correspondence should be addressed.
Submission received: 18 March 2016 / Revised: 1 July 2016 / Accepted: 11 July 2016 / Published: 20 July 2016

Abstract

:
In this paper, we consider retrieving individual wave components in a multi-directional sea wave model. To solve this problem, a currently and commonly used method is three-dimensional discrete Fourier transform (3D DFT) on the radar image sequence. However, the uniform frequency and the uniform wavenumber in a wavenumber frequency domain can not always strictly satisfy the dispersion relation, and the spectral leakage in both temporal and spatial domains exists due to the limited analysis area selected from an image sequence. As a result, the DFT method incurs undesirable error performance in retrieving directional wave components. By deeply investigating the data structure of the multi-directional sea wave model, we obtain a new and decomposable matrix representation for processing the wave components. Then, a novel successive cancellation method is proposed to efficiently and effectively extract individual wave components, whose frequency and wavenumber rigorously satisfy the liner dispersion relation. Thus, it avoids spectral leakage in the spatial domain. The algorithm is evaluated by using linear synthetic wave image sequences. The validity of the proposed novel algorithm is verified by comparing the retrieved parameters of amplitude, phase, and direction of the individual wave components with the simulated parameters as well as those obtained by using the 3D DFT method. In addition, the reconstructed sea field using the retrieved wave components is also compared with the simulated remote sensing images as well as those attained using the inverse 3D DFT method. All the simulation results demonstrate that our proposed algorithm is more effective and has better performance for retrieving individual wave components from the spatio-temporal remote sensing image sequences than the 3D DFT method.

Graphical Abstract

1. Introduction

The ocean wave is one of the most common ocean undulation phenomena and is critically important for ocean research and ocean exploitation [1,2,3,4,5,6,7]. Many offshore operations like the tensioning of a tanker, cargo transfer, and helicopter landing demand a deterministic prediction of the sea surface field in real-time [8,9,10,11,12,13,14,15]. In [13], the accurate real-time prediction of the sea surface field at the position of an operation platform is completely dependent on the retrieved sea wave components from a measured area in the vicinity of the platform [9,10,11,12,13,14,15,16]. In our paper, the wave component denotes the individual sine or cosine waves, since the sea surface can be seen as the superposition of numerous sine or cosine wave. The extracted sea wave components could be directly used as an input to initialize a wave propagation model, and then, to predict a real-time platform motion [13,14]. In addition, the sea wave components provide a low-cost approach to investigating the spatial and temporal development of 3D ocean surface fields [15] and the properties of individual wave components [17]. Hence, the accuracy of the retrieved sea wave components is one of the most important aspects in the real-time prediction of the sea surface field.
Thus far, most research on retrieving the sea wave components focuses on X-band marine radar for its economy and popularity, and the linear wave model because of its computational efficiency [9,10,11,13,14,15,16]. The X-band marine radar is commonly used to monitor moving vessels and to derive the statistical sea state parameters such as significant wave height, peak period, wave direction, bathymetry and ocean surface currents from coastal areas or moving platforms [17,18,19,20,21,22,23,24,25]. Now, with recent advances in computer technologies and marine radar techniques, it has the ability to retrieve the sea wave components from an X-band marine radar image sequence [16,19].
The method of retrieving the sea wave components, which can be fed into the wave propagation model, was first proposed in [11,12]. The 1D time series of the wave components were derived from the retrieved sea surface field. In [16,19], by using a 3D discrete Fourier transform (DFT) on the acquired X-band marine radar image sequence, a passband filter was used to select wave-related components in a wavenumber frequency domain. Since the X-band marine radar image is mainly modulated by hydrodynamic modulation, tilt modulation and shadowing, the sea clutter radar image intensities do not stand for the ocean surface field [1,2,3,4,5,6,7,9,15]. An empirical modulation transfer function (MTF) is needed to calibrate the amplitude spectrum [21,22,23,24,25]. Subsequently, an inverse 3D DFT method was used to reconstruct the sea surface field, which was used to retrieve the significant wave height. The sea surface field can also be obtained by summing all the individual wave components selected from the passband filter [11,12]. Compared with the method of summing all the individual wave components, the inverse 3D DFT method has the advantage of high accuracy and reconstructs the sea surface field quickly. In [3,17], the empirical methods to retrieve the ocean surface field in spatial and temporal domains were proposed. The radar image sequence was also decomposed into a wavenumber frequency domain with the 3D DFT method. In [3,17], the tilt modulation was considered as a major effect on retrieving the sea surface field. In [16,19], shadowing that has a significant influence on the radar image was taken into account at grazing incidence.
The individual wave components could be directly determined by the wave amplitude spectra after the filtering process and applying the modulation transfer function on the amplitude and phase spectra in [3,16,17,19]. In order to provide an effective prediction, the input requirement of the wave propagation model is that the frequency and wavenumber of the wave components should exactly match the dispersion relation. However, neither the 1D time series of wave components, which is retrieved from the sea surface field, nor the individual wave components, which can be directly retrieved from the frequency wavenumber spectrum after a 3D DFT on the radar image sequence, follow the requirement. Since the 3D DFT was taken on the radar image sequence in the spatial and temporal domain, the uniform frequency and the uniform wavenumber do not strictly satisfy the dispersion relation because of the limited analysis area selected from the image sequence, even though perfect linear synthetic wave data was used [13,14]. In addition, the spectral leakage existed in both the spatial and temporal domains, since the analysis area selected was finite [13,26], which may become a problem for estimating the wave surface field in the position of a platform. Furthermore, the empirical MTF requires an external equipment as the reference which increases the expense. Using the fixed MTF to different amplitude spectrum may also produce errors for retrieving wave components. Therefore, the resulting errors were introduced into the retrieved wave components, which could deteriorate the prediction of the wave surface field and the motion estimation of the offshore platform.
In order to make the retrieved wave components strictly satisfy the dispersion relation, various filtering methods have been put forward to manipulate the wavenumber and the frequency of the spectrum after taking the 3D DFT on the radar image sequence [13,16]. Nevertheless, the accuracy of the reconstructed wave field and the retrieved wave components were not optimal, since a bandpass filter was used and the extracted wave components were susceptible to the bandpass filter. To satisfy the requirement of the wave propagation model, a 2D DFT is applied to retrieve the directional wave components in [27]. The synthetic radar images, which were synthesized with linear waves while taking tilt modulation and shadowing, were used to verify its effectiveness. For this method, the significant wave height should be considered as prior information. In [10], a dynamic averaging and evolution scenario that integrates the inversion of the acquired radar images was used to provide a sea surface field prediction. A 4D Var assimilation method was proposed to predict the sea surface field in [28]. In [15], based on a multi-directional sea wave model, a method using the 2D DFT on the acquired radar image to predict the sea surface field at a desired location and time was proposed.
Our primary goal in this paper is to propose a new successive cancellation algorithm, which is an alternative to the 3D DFT method and fully based on theory, to efficiently and effectively extract the individual wave components without any reference measurements, whose frequency and wavenumber rigorously satisfy the linear dispersion relation and which can be directly used to initialize the wave propagation model.
This paper is organized as follows: Section 2 presents the linear dispersion relation and the sea wave model. The detailed derivation and the successive cancellation algorithm to recover the wave components based on the sea wave model are given in Section 3. In Section 4, the validity of the proposed method is investigated by computer simulation. Finally, the conclusion is summarized in Section 5.
Notation: Scalars, vectors and matrices are expressed by regular, bold lowercase letters and bold uppercase letters, respectively. A ( p , q ) stands for the ( p , q ) element of a matrix, A . The superscript ( · ) H , ( · ) T and ( · ) * indicate the Hermitian transpose, transpose and complex conjugate, respectively. F [ · ] denotes the DFT operator of the expression in brackets. | | · | | F stands for the Frobenius norm (F-norm or norm-two). Re [ · ] represents the real part of a complex number. C denotes the complex domain. C ¯ is a scalar and H ¯ is a matrix.

2. Sea Wave Model

In this section, we briefly review the concepts of the linear wave theory and the sea wave model.
It is known that under the assumption of the linear wave theory, oceanic surface gravity waves are dispersive with the dispersion relation, given by [1,21,22]
w = g k tanh ( k h )
where w is an angular frequency, k represents the modulus of wavenumber, g is the gravity acceleration and h denotes a water depth. For gravity waves in infinitely deep water, Equation (1) can be simplified as
w = g k
The ocean surface waves can be seen as the superposition of a number of wave components with different frequencies, amplitudes, wave directions and phases. Based on a linear wave model [15,16,25], a wave field, H ¯ ( x , y , t ) , at spatial coordinates x, y, and a temporal coordinate t is given by
H ¯ ( x , y , t ) = n = N / 2 N / 2 m = 0 M 1 A ( w n , θ m ) · cos ( w n t x k n cos θ m y k n sin θ m ε m n )
where N is the frequency number, M is the number of wave component directions, A ( w n , θ m ) = 2 S ( w n , θ m ) w θ stands for a discrete amplitude, S ( w n , θ m ) is a directional wave spectrum, w n denotes a wave frequency, θ m is a wave component direction, k n is the modulus of wavenumber and ε m n is a uniformly distributed initial phase on an interval [ π , π ) . Since the wave images are composed of discrete data, we divide the wave frequency and wavenumber into N + 1 parts. The wave angular frequency can be written as
w n = n w
where w = w is angular frequency resolution. For presentation clarity, Equation (3) can be recasted as
H ¯ ( x , y , t ) = 1 2 [ H ( x , y , t ) + H * ( x , y , t ) ] = 1 2 n = N / 2 N / 2 C n ( x , y ) · e j n w t + 1 2 n = N / 2 N / 2 C n * ( x , y ) · e j n w t
where
H ( x , y , t ) = n = N / 2 N / 2 m = 0 M 1 A ( w n , θ m ) · [ cos ( n w t x k n cos θ m y k n sin θ m ε m n ) + j sin ( n w t x k n cos θ m y k n sin θ m ε m n ) ] = n = N / 2 N / 2 m = 0 M 1 A ( w n , θ m ) · e j ( x k n cos θ m + y k n sin θ m + ε m n ) · e j n w t
and
C n ( x , y ) = m = 0 M 1 A ( w n , θ m ) · e j ( x k n cos θ m + y k n sin θ m + ε m n )

3. Retrieving the Individual Wave Components

The sea surface field consists of a linear superposition of several individual wave components. If we can obtain each linear wave component whose angular frequency and wavenumber strictly match the dispersion relation, then the retrieved linear wave components can be directly utilized as an input to initialize the wave propagation model, and then to predict the sea surface field at the position of an operation platform and estimate a platform motion. Meanwhile, the sea surface field, which has free selection of temporal and spatial resolution, can be reconstructed by summing all the individual wave components. Our main purpose in this section is to propose a novel method to attain all the wave components by using the sea wave model.

3.1. 1D DFT on the Radar Image Sequence

When the spatial analysis area of a radar image is selected, the sea surface field can be sampled and written as a 2D matrix H ¯ t for time t. Suppose that M and N represent the sampling numbers in the x and y directions of the radar image, respectively. Then, the wave field can be represented as H ¯ t ( x p , y q ) , where p { 0 , 1 , , M 1 } and q { 0 , 1 , , N } . Thus, Equation (5) can be represented in a more compact matrix form as
H ¯ t = 1 2 ( H t + H t * ) = 1 2 n = N / 2 N / 2 C n e j n w t + n = N / 2 N / 2 C n * e j n w t = 1 2 n = N / 2 N / 2 C ¯ n e j n w t
where C ¯ n = C n + C n * , the ( p , q ) elements of H ¯ t and C ¯ n are H ¯ t ( x p , y q ) and C n ( x p , y q ) + C n * ( x p , y q ) , respectively. Now, we discretize the time t with a sample interval Δ t , which coincides with the radar rotation time. The frequency resolution can be represented as w = 2 π ( N + 1 ) Δ t and then, the radar image sequence H ¯ is achieved, where t s = 2 π s ( N + 1 ) w and s { N / 2 , N / 2 + 1 , , N / 2 } . After that, we take DFT at each location of the sea surface field so that
F [ H ¯ ] = DFT 1 2 n = N / 2 N / 2 C ¯ n e j n w t s = 1 2 s = N / 2 N / 2 n = N / 2 N / 2 C ¯ n e j n w 2 π s ( N + 1 ) w · e j 2 π n s N + 1 = N + 1 2 C ¯ n ,
from which an M by N Fourier coefficient matrix C ¯ n can be attained. In this case, the sea surface field is expressed as Fourier series.

3.2. Fourier Coefficient Matrix Decomposition

Just as we have seen from Equation (6), by taking the 1D DFT on the radar image sequence H ¯ , we obtain C ¯ n . The ( p , q ) element of C ¯ n can be explicitly rewritten as
C ¯ n ( p , q ) = C n ( p , q ) + C n * ( p , q ) = m = 0 M 1 e j x p k n cos θ m · d m , n · e j y q k n sin θ m + m = 0 M 1 e j x p k n cos θ m · d m , n * · e j y q k n sin θ m ,
where d m , n = A ( w n , θ m ) · e j ε m n , which contains the amplitude and the phase information, are called modulated parameters in this paper. Based on the expression of Equation (7), it would be important to notice that the Fourier coefficient matrix C ¯ n can be decomposed into two matrices terms C n and C n * and each term can be subsequently decomposed as the product of matrices, which is represented as
C ¯ n = U n D n V n + U n * D n * V n *
where C n = U n D n V n and C n * = U n * D n * V n * . D n is the modulated parameters matrix that is equal to diag d 0 , n d 1 , n d M 1 , n , U n = u 0 ( n ) u 1 ( n ) u M 1 ( n ) C M × M and V n = v 0 ( n ) v 1 ( n ) v M 1 ( n ) T C M × N . Vector u m ( n ) and v m ( n ) are given by u m ( n ) = e j x 0 k n cos θ m e j x 1 k n cos θ m e j x M 1 k n cos θ m T C M × 1 and v m ( n ) = e j y 0 k n sin θ m e j y 1 k n sin θ m e j y N 1 k n sin θ m C 1 × N , respectively. Therefore, the sea surface elevation with frequency w n denoted by Equation (8) can be rewritten as
C ¯ n = u 0 ( n ) d 0 , n v 0 ( n ) + u 1 ( n ) d 1 , n v 1 ( n ) + + u M 1 ( n ) d M 1 , n v M 1 ( n ) + u 0 ( n ) * d 0 , n * v 0 ( n ) * + u 1 ( n ) * d 1 , n * v 1 ( n ) * + + u M 1 ( n ) * d M 1 , n * v M 1 ( n ) *
Here, it is observed clearly that each term of Equation (9) only involves one wave component direction θ m , which corresponds to two terms u m ( n ) d m , n v m ( n ) and u m ( n ) * d m , n * v m ( n ) * . For the fixed frequency w n , the sea surface elevation C ¯ n and the wavenumber k n are given. Thus, Equation (9) is a function of modulated parameters d m , n and d m , n * and wave component direction θ m . Suppose that the two terms containing θ m are moved to the left side of Equation (9) and then the right side is the terms involving θ m , m { 0 , 1 , , M 1 } and m m . For the left side, we take an F-norm on it, and, then, by solving a optimization problem, we take the derivative with respect to modulated parameters d m , n and d m , n * . Thus, a relation between modulated parameters and direction θ m can be achieved. Subsequently, substituting the relation into F-norm expression, it becomes a function that only involves θ m . In the previous section, we have discretized the wave component direction into the M direction. Therefore, we can find a direction θ m that minimizes the F-norm of θ m expression by searching the M discrete component directions. After that, the value of modulated parameters d m , n and d m , n * can be obtained through the above achieved relation and the chose direction θ m , and, then, the amplitude A ( w n , θ m ) and the phase ε m n of the sea wave component are determined.
Up to now, the left side is known to us. Then, we continue to shift the terms containing direction θ m + 1 to the left side. Similarly, all the directions, amplitudes and initial phases of the wave components can be achieved by repeating the above process. The detailed derivation and analysis of the successive cancellation algorithm are given in Appendix A. The complete successive cancellation algorithm for extracting all the wave components is now depicted below.

3.3. Successive Cancellation Algorithm

By investigating the data structure uncovered by the matrix decomposition in Equation (9) and the above analysis, it motivates us to propose the following successive cancellation algorithm, which first optimizes the direction parameters θ m by iteratively moving the terms containing the direction parameters θ m and then estimates the corresponding modulated parameters d m , n .
  • Given R 0 ( n ) = C ¯ n , initially retrieve modulated parameters d 0 , n , d 0 , n and direction parameter θ 0 by solving the following optimization problem:
    min R 0 ( n ) u 0 ( n ) d 0 , n v 0 ( n ) + u 0 ( n ) * d 0 , n * v 0 ( n ) * F
    where R m ( n ) = C ¯ n u m ( n ) d m , n v m ( n ) u m ( n ) * d m , n * v m ( n ) * (more detail is in Appendix A).
  • Suppose that all previous extractions of d m , n , d m , n and θ m for m = 0 , 1 , , k are completely correct. Now, successively retrieve the modulated parameters d k + 1 , n , d k + 1 , n , and direction parameter θ k + 1 by solving the following optimization problem:
    min R k ( n ) u k + 1 ( n ) d k + 1 , n v k + 1 ( n ) + u k + 1 ( n ) * d k + 1 , n * v k + 1 ( n ) * F
    for k = 1 , 2 , , M 2 , where R k ( n ) = C ¯ n m = 0 k u m ( n ) d m , n v m ( n ) + u m ( n ) * d m , n * v m ( n ) * In addition, the solution in Equation (10) can be efficiently obtained by the following two steps:
    (a)
    The direction parameters θ k + 1 are extracted by solving the optimization problem below:
    min θ k + 1 n R k ( n ) u k + 1 ( n ) d k + 1 , n v k + 1 ( n ) + u k + 1 ( n ) * d k + 1 , n * v k + 1 ( n ) * F
    (b)
    The modulated parameters d k + 1 , n and d k + 1 , n are determined by d ¯ k + 1 ( n ) = A k + 1 ( n ) 1 b k + 1 ( n )
    where
    d ¯ k + 1 ( n ) = d k + 1 , n d k + 1 , n * T A k + 1 ( n ) = Tr u k + 1 ( n ) v k + 1 ( n ) v k + 1 ( n ) H u k + 1 ( n ) H Tr u k + 1 ( n ) * v k + 1 ( n ) * v k + 1 ( n ) H u k + 1 ( n ) H , Tr u k + 1 ( n ) v k + 1 ( n ) v k + 1 ( n ) T u k + 1 ( n ) T Tr u k + 1 ( n ) * v k + 1 ( n ) * v k + 1 ( n ) T u k + 1 ( n ) T
    and
    b k + 1 ( n ) = u k + 1 ( n ) H C ¯ n v k + 1 ( n ) H u k + 1 ( n ) T C ¯ n v k + 1 ( n ) T
Thus, the parameters of amplitude A ( w n , θ m ) , initial phase ε m , n and direction θ m of the individual wave components of the sea wave model can be obtained by the above successive cancellation algorithm. In the algorithm, for each index m, it traverses all the M wave component directions to search for the suitable discrete direction value θ m , and then it iterates M times to attain all the wave component directions. Hence, the overall complexity of the algorithm is O ( M 2 ) .

4. Results and Analysis of Experiment

In this section, we carry out radar image simulations at different sea state conditions to examine the validity of our proposed method for retrieving the individual wave components. In order to simplify the problem, the perfect linear synthetic wave data is used. The performance of our method recovering the amplitude and the phase parameters of the individual wave components is presented. The sea surface field is reconstructed by summing all the retrieved wave components achieved. Moreover, the performance of the 3D DFT method retrieving the sea wave components and the sea surface field reconstructed by using the 3D inverse DFT method are presented as a reference. The reconstruction accuracy of the sea surface field of our successive cancellation algorithm is shown by comparing with that of the other method. The experiment and analysis results are shown as follows.

4.1. Radar Image Sequence Simulated with Single Wave

For the primary test, the radar image sequence is simulated with a monochromatic wave. The wave parameters with the significant wave height A = 2 m, wave length λ = 160.2 m ( k n = 0.0392 rad/m), wave direction θ = 90 , wave period T = 10.13 s ( w n = 0.6203 rad/s) and initial phase ε = 72 are used. In this study, a 3D 128 × 127 × 127 image sequence is simulated. The spatial resolution of the image is Δ x = Δ y = 7.5 m. For the time domain of the wave image sequence, the time resolution is Δ t = 2.47 s. After 1D DFT on the radar image sequence at each position, the Fourier coefficient matrix is the input of our novel iterative algorithm. Based on our method, all the parameters of the wave component including the component direction, the amplitude, initial phase, the wave angular frequency and the wavenumber can be almost exactly recovered. However, it should mentioned that, in our experiment, there is a 180 degree direction ambiguity, since the minimized n = N / 2 , n 0 N / 2 R m ( n ) F 2 always corresponds to two different wave directions. The reason might be that the F-norm used in our method possibly causes this ambiguity. In practice, the dominant wave direction can be easily known beforehand, and the directions of most wave components are close to it. Thus, we can choose the one that is relatively close to the dominant wave direction. It is validated that the novel successive cancellation algorithm has the capability to estimate the single wave component from the radar image sequence.

4.2. The Simulation of Radar Image Sequence with Random Wave Field

Real waves in nature are always random and irregular. Our acquired radar image sequence which does not satisfy the input requirement has only 32 continuous images. Hence, we have to use the simulated radar image sequence to certify the effectiveness of our method. The simulation of the sea field is commonly based on the Joint North Sea Wave Project (JoNSWaP) spectrum S ( f ) [16,25]. Since S ( f ) is in 1D form and the multi-directional characteristic of ocean waves, a Mitsuyasu-type directional spreading function G ( f , θ ) is used:
S ( f , θ ) = S ( f ) G ( f , θ )
G ( f , θ ) = 1 π 2 2 s 1 Γ 2 ( 2 s + 1 ) Γ ( 2 s + 1 ) c o s 2 s ( ( θ θ p ) / 2 )
where S ( f , θ ) is the directional wave spectrum, θ p is the dominant wave direction, and Γ is the gamma function. The parameter s leads to the angular distribution and varies with the frequency. The wave conditions are initialized with the significant wave height H s = 3.5 m, the dominant wave direction θ p = 104 and the peak wave period T p = 10.99 s. The angular frequency and the direction of the wave components are uniformly discretized with w = 0.03 rad/s and θ = 11 , respectively. Then, the simulated amplitude of the wave components can be obtained by the directional wave spectrum, which is shown in Figure 1a. The temporal resolution of the wave image sequence is about Δ t = 1.65 s. The image sequence contains 127 continuous images, where N = 126 . This means that about 4000 individual wave components are considered. Thus, the 3D 128 × 127 × 127 radar image sequence is generated by summing all the linear wave components as Equation (3). The spatial resolution is Δ x = Δ y = 7.5 m. The hydrodynamic and tilt modulation have a minor impact on the imaging mechanism at grazing incidence, while the shadowing modulation plays an important role [16,19]. Hence, we only take shadowing into account in our simulation. Here, we assume that an antenna is 20 m above the mean sea level and the location of the selected region from the radar platform is 300 m. Since the shadowing mask is time variant and the real radar image is non-stationary, it is unreasonable to retrieve parameters by taking 1D DFT on the image sequence without any calibration in subsequent processing. Here, we use ten fixed shadowing masks to produce ten radar image sequences based on the simulated sea surface field sequence, and then we average the retrieved parameters from the ten radar image sequences in our simulation. In our future research, we will develop a time frequency analysis method to directly retrieve parameters from the marine radar image sequence. The simulated radar image is shown in Figure 1b at t = 50 s.
The 3D DFT method and our method are both used to retrieve the amplitude and the phase from the radar image sequence. Figure 2 describes the retrieved amplitude of the wave components with these two different methods. Since the energy points do not strictly follow the dispersion relation in Equation (2) after a 3D DFT on the radar image, a passband filter is used to extract the amplitude and the phase of the wave components. The estimated amplitude using the 3D DFT method is shown in Figure 2a. Figure 2b presents the estimated results of the amplitude using our algorithm.
To verify the retrieved accuracy of the amplitude, the error between the simulated amplitude and the retrieved amplitude is presented in Figure 3. Figure 3a shows the error between the simulated amplitude and the amplitude retrieved by the 3D DFT method. Figure 3b shows the error between the simulated amplitude and the amplitude retrieved by our method. From Figure 2 and Figure 3, it can be observed that the retrieved amplitude using the 3D DFT has a larger error than that of our method. Therefore, we can make the conclusion that our method has a relatively better performance than the 3D DFT method in terms of amplitude estimation.
Similarly, the phase of the individual wave components can be also retrieved with our method and the 3D DFT method. Figure 4 presents the phase error of the retrieved individual wave components using these two methods. In order to better describe the phase error, Figure 4 only presents the phase error of the dominant wave components, which means that the phase of individual wave component having a small amplitude is not taken into account. The difference between the simulated initial phase and the retrieved phase based on the 3D DFT method is given in Figure 4a, while the difference between the simulated initial phase and the recovered phase based on our method is given in Figure 4b. From Figure 4, we can observe that the retrieved phase using 3D DFT method incurs larger errors than that of our successive cancellation method. Hence, we can say that our method also shows a relatively better performance than the 3D-DFT method in terms of phase estimation.
Despite the fact that we have compared the retrieved amplitude and phase between the 3D DFT method and our method, it is not reasonable to assess the performance only by these two parameters. On the other hand, the sea surface field based on the 3D DFT method can be reconstructed by the inverse 3D DFT on the wave spectrum or summing all the retrieved individual wave components, and the inverse 3D DFT method is more accurate than the other method [11], since the superposition of the inaccurate wave components that are extracted by a traditional dispersion relation bandpass filter [2,13] would deteriorate the reconstructed wave field. Hence, we can evaluate the performance of the retrieved wave components by the reconstructed field. To verify the effectiveness of our newly proposed algorithm, we only need to compare the sea surface field reconstructed by summing all individual wave components retrieved from our successive cancellation algorithm with that reconstructed by the inverse 3D DFT method on the wave spectrum after a bandpass filtering. The simulated wave field and the reconstructed wave fields described by these two methods at t = 50 s are given in Figure 5. Figure 5a is simulated wave field at t = 50 s using the amplitude in Figure 1a. The red color denotes wave crest and the blue color represents wave trough. Figure 5b shows the reconstructed sea surface field by the inverse 3D DFT method. The sea surface field reconstructed by summing all the individual wave components retrieved based on our successive cancellation method is presented in Figure 5c. It is observed clearly that the reconstructed sea wave fields using the two methods are all similar to the simulated sea wave field in Figure 5a. In addition, the reconstruction error of the sea surface field of the two methods is shown in Figure 6.
In order to further verify the validity of the proposed algorithm, 500 realizations with different sea states and directional resolution θ = 11 are used. The average reconstruction error is expressed by
ε r r ( x , y , t ) = 1 N R i = 1 N R | η i ( x , y , t ) η ^ i ( x , y , t ) | 2
where N R = 500 is the number of realizations, η i ( x , y , t ) and η ^ i ( x , y , t ) represent the synthetic wave field and the recovered wave field, respectively and ε r r ( x , y , t ) denotes the average reconstruction error. Taking the synthetic wave image sequences as the ground truth, the average reconstruction error with the two methods are shown in Figure 7 at three different times for t = 0 s, t = 50 s and t = 100 s. Figure 7a–c show the average reconstruction error between the synthetic sea field and the recovered wave field with the 3D DFT method. Figure 7d–f present the average reconstruction error between the synthetic sea field and the recovered wave field with our method. It can be observed that the reconstructed sea surface field using the retrieved individual wave components always shows a minor error compared with that of the inverse 3D DFT method at different times. It also implies that the successive cancellation method has a better performance in terms of exactly recovering sea wave components whose amplitudes and initial phases are almost successfully extracted using our method.
Since the 3D DFT method assigns energy to the components whose wavenumber k and angular frequency w do not necessarily satisfy the linear dispersion relation and the spectral leakage exists in both temporal and spatial domains, this results in the fact that the amplitude and initial phase extracted by the 3D DFT method have large errors [3,13,14,16,17,19,26]. However, our new successive cancellation algorithm presents smaller error in extracting the wave components, since we attain k from w through the linear dispersion relation directly, which necessarily assures that our method has the capability to accurately extract wave components whose frequency and wavenumber rigorously satisfy the liner dispersion relation from the radar image sequence. In addition, in our method, the spectral leakage in the spatial domain is avoided, since we only take 1D DFT on the temporal domain and apply the successive cancellation algorithm in the spatial domain, which is alternative to the 2D DFT method.

4.3. The Experiment with High Directional Resolution

When the more wave directions are taken in the directional wave spectrum, the angular frequency and the wave component direction are uniformly discretized with w = 0.03 rad/s and θ = 5.5 , respectively. The simulated amplitude and the recovered amplitude are given in Figure 8. Figure 8a is the simulated amplitude based on the JoNSWaP spectrum. The amplitude recovered with our method is shown in Figure 8b. With the increase of the directional resolution, the experiment shows that it is unable to resolve all the directions, since the directions are too close to be resolved using our method. It can be seen that not all the amplitudes can be recovered and the recovered amplitude presents an obvious error. In practice, there is an large error in the retrieved phase, which is not presented here.
In spite of the fact that our method is unable to recover every wave component, the reconstructed sea field using the obtained inaccurate wave components is also presented. The simulated sea field and reconstructed sea surface field with the two methods at t = 100 s are shown in Figure 9. Figure 9a is the simulated sea field at t = 100 s. The recovered sea field with the inverse 3D DFT method is shown in Figure 9b. Figure 9c is the reconstructed sea field by summing all the attained inaccurate wave components based on our method. The result is interesting, since the reconstructed sea surface field based on our method is in accordance with the simulated sea field in Figure 9a as well as that based on the inverse 3D DFT method in Figure 9b. Moreover, the reconstruction error of the two methods is shown in Figure 10. The experimental result is that the reconstruction error of our proposed method is smaller than that of the inverse 3D DFT method, which is the same as the above in Figure 8.
Therefore, when more wave directions are taken in the directional wave spectrum, the experiment shows that our method fails to resolve all the directions of the individual wave components, since the directions are too close to be resolved using the successive cancellation method. However, the reconstructed sea surface field is still accurate. Based on the experimental result, we conjecture that each retrieved inaccurate wave component might be the superposition of a few of wave components simulated in some certain range of angles.

4.4. Results Analysis

The accuracy of our method is investigated by analyzing the extracted wave components and the sea wave field that is reconstructed by the sum of all the achieved wave components. Under the condition of low directional resolution, the experiments above show that our proposed successive cancellation algorithm is capable of accurately extracting the parameters of the wave component including the amplitude and phase from the simulated sea wave image sequences, which has less error than that of the 3D DFT method. With high directional resolution, our method can still accurately reconstruct the sea wave field. However, it should be mentioned clearly that these results are based on the perfect data. The accuracy and effectiveness of our algorithm are verified with a rather generic and idealized test. Therefore, more experiments are needed to verify its accuracy in various sea conditions.

5. Conclusions

In this paper, we have proposed a novel successive cancellation algorithm to extract wave components by revealing the decomposable data structure of the sea wave model. The detailed analysis and derivation are presented. Compared with the traditional 3D DFT method, our proposed method only takes DFT in the temporal domain, thus avoiding the spectral leakage in the spatial domain. The computer simulations have demonstrated that it is capable of accurately extracting wave components from the radar image sequence, and it is more accurate than the currently and commonly used 3D DFT method. However, these results are based on the perfect data. The accuracy and effectiveness of our algorithm are verified with a rather generic and idealized test. More experiments are needed to verify its accuracy in various sea conditions. Future research will focus on more realistic radar images.
Our method takes much time to reconstruct the sea field by summing all the wave components, but the angular frequency w and wavenumber k of the extracted wave component strictly satisfy the linear dispersion relation and can be directly fed into the wave propagation model. Hence, our method can be used to investigate wave propagation in the future. In addition, another approach to solve the ambiguity is to use norm-one optimization or some other signal processing methods, which is under our future investigation.

Acknowledgments

The work was supported by the National Natural Science Foundation of China (No. 11405035), the Special Funding Project of Marine Nonprofit Industry Research (No. 201405022-1) the China Scholarship Council (CSC) and the Natural Sciences and Engineering Research Council of Canada (NSERC). The authors would like to thank the anonymous reviewers for their very valuable suggestions and comments that have improved the quality of work.

Author Contributions

Yanbo Wei and Jian-Kang Zhang conceived and designed the experiments; Yanbo Wei performed the experiments; Yanbo Wei and Zhizhong Lu analyzed the data; Jian-Kang Zhang and Zhizhong Lu contributed reagents/materials/analysis tools; Yanbo Wei and Jian-Kang Zhang wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

We move the terms containing direction θ m on the right side of Equation (9) to the left side and let matrix R m ( n ) = C ¯ n u m ( n ) d m , n v m ( n ) u m ( n ) * d m , n * v m ( n ) * . Then, the square of the F-norm of R m ( n ) is expressed as
R m ( n ) F 2 = Tr C ¯ n u m ( n ) d m , n v m ( n ) u m ( n ) * d m , n * v m ( n ) * C ¯ n u m ( n ) d m , n v m ( n ) u m ( n ) * d m , n * v m ( n ) * H = Tr ( C ¯ n C ¯ n H ) 2 Re · Tr C ¯ n u m ( n ) d m , n v m ( n ) + u m ( n ) * d m , n * v m ( n ) * H + Tr u m ( n ) d m , n v m ( n ) + u m ( n ) * d m , n * v m ( n ) * u m ( n ) d m , n v m ( n ) + u m ( n ) * d m , n * v m ( n ) * H = Tr ( C ¯ n C ¯ n H ) 2 Re d m , n * Tr C ¯ n v m ( n ) H u m ( n ) H + d m , n Tr C ¯ n v m ( n ) T u m ( n ) T + d m , n d m , n * Tr u m ( n ) v m ( n ) v m ( n ) H u m ( n ) H + d m , n * d m , n Tr u m ( n ) * v m ( n ) * v m ( n ) T u m ( n ) T + d m , n d m , n Tr u m ( n ) v m ( n ) v m ( n ) T u m ( n ) T + d m , n * d m , n * Tr u m ( n ) * v m ( n ) * v m ( n ) H u m ( n ) H = Tr ( C ¯ n C ¯ n H ) 2 Re d m , n * u m ( n ) H C ¯ n v m ( n ) H + d m , n u m ( n ) T C ¯ n v m ( n ) T + d m , n d m , n q = 0 N e j 2 y q k n sin θ m p = 0 M 1 e j 2 x p k n cos θ m + d m , n * d m , n * q = 0 N e j 2 y q k n sin θ m p = 0 M 1 e j 2 x p k n cos θ m + d m , n d m , n * M ( N + 1 ) + d m , n * d m , n M ( N + 1 ) .
We let d ¯ m ( n ) = [ d m , n d m , n * ] T , which is the modulated parameter vector. Then, we have
R m ( n ) F 2 = d ¯ m ( n ) H A m ( n ) d ¯ m ( n ) 2 Re d ¯ m ( n ) H b m ( n ) + Tr ( C ¯ n C ¯ n H ) ,
where A m ( n ) is a 2 by 2 Hermitian matrix given by
A m ( n ) = M ( N + 1 ) q = 0 N e j 2 y q k n sin θ m p = 0 M 1 e j 2 x p k n cos θ m q = 0 N e j 2 y q k n sin θ m p = 0 M 1 e j 2 x p k n cos θ m M ( N + 1 ) ,
and b m ( n ) is a 2 by 1 vector defined by
b m ( n ) = u m ( n ) H C ¯ n v m ( n ) H u m ( n ) T C ¯ n v m ( n ) T .
The R m ( n ) F 2 in Equation (A1) is a function of modulated parameter d ¯ m ( n ) and wave component direction θ m . If we take the first order derivative of Equation (A1) with respect to d ¯ m ( n ) H and set it to equal to zero, we can obtain
A m ( n ) d ¯ m ( n ) b m ( n ) = 0 .
From Appendix B, we do not need to consider the case n = 0 , since it is highly impractical and we know that A m ( n ) is always invertible for n 0 . Therefore, when n 0 , we have
d ¯ m ( n ) = A m ( n ) 1 b m ( n ) .
Substituting Equation (A3) into Equation (A1) yields
R m ( n ) F 2 = b m ( n ) H A m ( n ) H b m ( n ) 2 Re b m ( n ) H A m ( n ) H b m ( n ) + Tr ( C ˜ n C ˜ n H ) = b m ( n ) T A m ( n ) T b m ( n ) * + Tr ( C ˜ n C ˜ n H ) .
In Equation (A4), the F-norm R m ( n ) F 2 only involves the direction parameter θ m . In the previous section, we have discretized the wave components into M directions. For each frequency plane n, we could find an appropriate direction θ m that minimizes Equation (A4) through the M directions. Thus, we would totally have N directions. Theoretically, min R m ( n ) F 2 for each n should have the same direction θ m for all the N frequency planes except n = 0 . In order to decrease the estimation error of the wave component direction, we search for a suitable direction θ m that minimizes n = N / 2 , n 0 N / 2 R m ( n ) F 2 by searching all the M directions. Thus far, the wave component direction θ m is chosen and d ¯ m ( n ) can be determined according to Equation (A3) after deriving the direction θ m . Subsequently, for n 0 , we can attain the amplitude A ( w n , θ m ) , which is equal to the modulus of d m , n , and the phase ε m , n , which is equal to the argument a r g ( d m , n ) .
Now, the terms involving θ m are canceled and the matrix R m ( n ) is known. Thus, we shift the terms embodying θ m + 1 and update R m ( n ) . Then, we attain R m + 1 ( n ) = R m ( n ) u m + 1 ( n ) d m + 1 , n v m + 1 ( n ) u m + 1 ( n ) * d m + 1 , n * v m + 1 ( n ) * . By repeating the above procedure, we recover all the parameters of Equation (3) including the amplitude A ( w n , θ m ) , the wave direction θ m and the phase ε m n . Now, the complete successive cancellation algorithm for extracting all the wave components is depicted.

Appendix B

Assume x p = p Δ x and y q = q Δ y , where Δ x and Δ y are the spatial resolution along x- and y-axis direction, respectively. In practice, we usually take Δ y = Δ x = Δ . It is observed clearly that each sum polynomial in matrix A m ( n ) is a sum of geometric series. Suppose that all the common ratio of geometric series e j 2 Δ k n sin θ m , e j 2 Δ k n sin θ m , e j 2 Δ k n cos θ m and e j 2 Δ k n cos θ m are not equal to 1, it means that sin θ m , cos θ m and wavenumber k n are not equal to 0. Then, matrix A m ( n ) can be rewritten as
A m ( n ) = M ( N + 1 ) 1 e j 2 ( N + 1 ) Δ k n sin θ m 1 e j 2 Δ k n sin θ m · 1 e j 2 M Δ k n cos θ m 1 e j 2 Δ k n cos θ m 1 e j 2 ( N + 1 ) Δ k n sin θ m 1 e j 2 Δ k n sin θ m · 1 e j 2 M Δ k n cos θ m 1 e j 2 Δ k n cos θ m M ( N + 1 ) .
The determinant of A m ( n ) is that
| A m ( n ) | = M 2 ( N + 1 ) 2 1 e j 2 ( N + 1 ) Δ k n sin θ m 1 e j 2 Δ k n sin θ m · 1 e j 2 ( N + 1 ) Δ k n sin θ m 1 e j 2 Δ k n sin θ m · 1 e j 2 M Δ k n cos θ m 1 e j 2 Δ k n cos θ m · 1 e j 2 M Δ k n cos θ m 1 e j 2 Δ k n cos θ m = M 2 ( N + 1 ) 2 1 cos ( 2 M Δ k n cos θ m ) 1 cos ( 2 Δ k n cos θ m ) · 1 cos ( 2 ( N + 1 ) Δ k n sin θ m ) 1 cos ( 2 Δ k n sin θ m ) = M 2 ( N + 1 ) 2 sin 2 ( M Δ k n cos θ m ) sin 2 ( Δ k n cos θ m ) · sin 2 ( ( N + 1 ) Δ k n sin θ m ) sin 2 ( Δ k n sin θ m ) > 0 .
Since sin θ m and cos θ m can not be equal to 0 simultaneously, it is observed clearly that | A m ( n ) | is greater than 0 when k n 0 . Hence, A m ( n ) is always invertible for n 0 . However, the determinant of | A m ( n ) | = 0 when k n = 0 . According to Equations (2) and (4), we attain
k n = w n 2 g = n 2 w 2 g = k n 2 .
Then, we have n = 0 and w n = 0 when k n = 0 . In practice, the case that the frequency and wavenumber of wave components are equal to 0 does not exist. Thus, we need not to consider the case for n = 0 .

References

  1. Young, I.; Rosenthal, W.; Ziemer, F. A three dimensional analysis of marine radar images for the determination of ocean wave directionality and surface currents. J. Geophys. Res. 1985, 90, 142–149. [Google Scholar] [CrossRef]
  2. Gangeskar, R. Ocean current estimated from X-band radar sea surface images. IEEE Trans. Geosci. Remote Sens. 2002, 40, 783–792. [Google Scholar] [CrossRef]
  3. Dankert, H.; Horstmann, J.; Rosenthal, W. Wind- and wave-field measurements using marine X-band radar-image sequences. J. Ocean. Eng. 2005, 30, 534–542. [Google Scholar] [CrossRef]
  4. Huang, W.; Gill, E.; An, J. Iterative least-squares-based wave measurement using X-band nautical radar. IET Radar Sonar Navig. 2014, 8, 853–863. [Google Scholar] [CrossRef]
  5. Serafino, F.; Lugni, C.; Soldovie, F. A novel strategy for the surface current determination from marine X-band radar data. IEEE Geosci. Remote Sens. Lett. 2010, 7, 231–235. [Google Scholar] [CrossRef]
  6. Shao, W.; Li, X.; Sun, J. Ocean wave parameters retrieval from TerraSAR-X images validated against buoy measurements and model results. Remote Sens. 2015, 7, 12815–12828. [Google Scholar] [CrossRef]
  7. Carratelli, E.P.; Chapron, B.; Dentale, F.; Reale, F. Simulating the Influence of Wave Whitecaps on SAR Images; European Space Agency: Paris, France, 2008. [Google Scholar]
  8. Carratelli, E.P.; Dentale, F.; Reale, F. Reconstruction of SAR Wave Image Effects through Pseudo Random Simulation; European Space Agency: Montreux, Switzerland, 2007. [Google Scholar]
  9. Dannenberg, J.; Naaijen, P.; Boom, H.V.D.; Reichert, K. The on board and motion estimator OWME. In Proceedings of the Twentieth (2010) International Offshore and Polar Engineering Conference, Beijing, China, 20–25 June 2010; pp. 424–431.
  10. Wijaya, A.P.; Naaijen, P.; Andonowati; Groesen, E. Reconstruction and future prediction of the sea surface from radar observations. Ocean Eng. 2015, 106, 261–270. [Google Scholar] [CrossRef]
  11. Reichert, K.; Dannenberg, J.; Boom, H.V.D. X-Band radar derived sea surface elevation maps as input to ship motion forecasting. In Proceedings of the IEEE OCEANS, Sydney, Australia, 24–27 May 2010.
  12. Hessner, K.; Reichert, K. Sea surface elevation maps obtained with a nautical X-Band radar—Examples from WaMoS II stations. In Proceedings of the 10th International Workshop on Wave Hindcasting and Forecasting and Coastal Hazard Symposium, Turtle Bay, HI, USA, 11–16 November 2007.
  13. Naaijen, P.; Blondel-Couprie, E. Reconstruction and Prediction of Short-Crested Seas Based on the Application of a 3D-DFT on Synthetic Waves: Part 1-Reconstruction. In Proceedings of the ASME 2012 31st International Conference on Ocean, Offshore and Arctic Engineering, Rio de Janeiro, Brazil, 1–6 July 2012.
  14. Blondel-Couprie, E.; Naaijen, P. Reconstruction and Prediction of Short-Crested Seas Based on the Application of a 3D-DFT on Synthetic Waves: Part 2-Prediction. In Proceedings of the OMAE2012, Rio de Janeiro, Brazil, 1–6 July 2012.
  15. Belmont, M.K.; Christmas, J.; Dannenberg, J.; Hilmer, T.; Duncan, J.; Duncan, J.M. An examination of the feasibility of linear deterministic sea wave prediction in multidirectional seas using wave profiling radar: Theory, simulation, and sea trials. J. Atmos. Ocean. Technol. 2014, 31, 1601–1614. [Google Scholar] [CrossRef]
  16. Nieto-Borge, J.C.; Rodriguez, G.; Hessner, K.; Gonzalez, P.I. Inversion of marine radar images for surface wave analysis. J. Atmos. Ocean. Technol. 2004, 21, 1291–1300. [Google Scholar] [CrossRef]
  17. Dankert, H.; Rosenthal, W. Ocean surface determination from X-band radar-image sequences. J. Geophys. Res. 2004, 109. [Google Scholar] [CrossRef]
  18. Shen, C.; Huang, W.; Gill, E.W.; Carrasco, R.; Horstmann, J. An algorithm for surface current retrieval from X-band marine radar images. Remote Sens. 2015, 7, 7753–7767. [Google Scholar] [CrossRef]
  19. Nieto-Borge, J.C.; Hessner, K.; Jarabo-Amores, P.; MataMoya, D. Signal-to-noise ratio analysis to estimate ocean wave heights from X-band marine radar image time series. IET Radar Sonar Navig. 2008, 2, 35–41. [Google Scholar] [CrossRef]
  20. An, J.; Huang, W.; Gill, E.W. A Self-Adaptive Wavelet-Based Algorithm for Wave Measurement Using Nautical Radar. IEEE Trans. Geosci. Remote Sens. 2015, 53, 567–577. [Google Scholar]
  21. Senet, C.; Seemann, J.; Ziemer, F. The near-surface current velocity determined from image sequences of the sea surface. IEEE Trans. Geosci. Remote Sens. 2001, 39, 492–505. [Google Scholar] [CrossRef]
  22. Gangeskar, R. An algorithm for estimation of wave height from shadowing in X-band radar sea surface images. IEEE Trans. Geosci. Remote Sens. 2014, 52, 3373–3381. [Google Scholar] [CrossRef]
  23. Nieto-Borge, J.C.; Guedes Soares, C. Analysis of directional wave fields using X-band navigation radar. Coast. Eng. 2000, 40, 375–391. [Google Scholar] [CrossRef]
  24. Ludeno, G.; Reale, F.; Dentale, F.; Carratelli, E.P.; Natale, A.; Soldovieri, F.; Serafino, F. An X-band radar system for bathymetry and wave field analysis in a harbour area. Sensors 2015, 15, 1691–1707. [Google Scholar] [CrossRef] [PubMed]
  25. Chuang, L.Z.-H.; Wu, L.-C. Study of wave group velocity estimation from inhomogeneous sea-surface image sequences by spatiotemporal continuous wavelet transform. IEEE J. Ocean. Eng. 2014, 39, 444–457. [Google Scholar] [CrossRef]
  26. Serafino, F.; Lugni, C.; Nieto-Borge, J.C.; Soldovieri, F. A simple strategy to mitigate the aliasing effect in X-band marine radar data: Numerical results for a 2D case. Sensors 2011, 11, 1009–1027. [Google Scholar] [CrossRef] [PubMed]
  27. Naaijen, P.; Wijaya, A.P. Phase resolved wave prediction from synthetic radar images. In Proceedings of the 33rd International Conference on Ocean, Offshore and Arctic Engineering, San Francisco, CA, USA, 8–13 June 2014.
  28. Aragh, S.; Nwogu, O. Variation assimilating of synthetic radar data into a pseudo-spectral wave model. J. Coast. Res. 2008, 52, 235–244. [Google Scholar] [CrossRef]
Figure 1. The simulated amplitude of the wave field and the generated radar image. (a) the simulated amplitude A ( w n , θ m ) of the wave components; (b) the simulated radar image at t = 50 s using the amplitude in Figure 1a.
Figure 1. The simulated amplitude of the wave field and the generated radar image. (a) the simulated amplitude A ( w n , θ m ) of the wave components; (b) the simulated radar image at t = 50 s using the amplitude in Figure 1a.
Remotesensing 08 00607 g001
Figure 2. The retrieved amplitude of the wave components. (a) the retrieved amplitude based on the 3D discrete Fourier transform (DFT) method; (b) the retrieved amplitude based on the novel method.
Figure 2. The retrieved amplitude of the wave components. (a) the retrieved amplitude based on the 3D discrete Fourier transform (DFT) method; (b) the retrieved amplitude based on the novel method.
Remotesensing 08 00607 g002
Figure 3. The amplitude error of the retrieved wave components. (a) the amplitude error based on the 3D DFT method; (b) the amplitude error based on the novel method.
Figure 3. The amplitude error of the retrieved wave components. (a) the amplitude error based on the 3D DFT method; (b) the amplitude error based on the novel method.
Remotesensing 08 00607 g003
Figure 4. The phase error of the dominant wave components. (a) the phase error based on the 3D DFT method; (b) the phase error based on the novel method.
Figure 4. The phase error of the dominant wave components. (a) the phase error based on the 3D DFT method; (b) the phase error based on the novel method.
Remotesensing 08 00607 g004
Figure 5. The sea surface field. (a) the simulated sea surface field at t = 50 s using the amplitude in Figure 1a; (b) the reconstructed sea surface field at t = 50 s based on the inverse 3D DFT method; (c) the reconstructed sea surface field at t = 50 s by summing all the retrieved wave components.
Figure 5. The sea surface field. (a) the simulated sea surface field at t = 50 s using the amplitude in Figure 1a; (b) the reconstructed sea surface field at t = 50 s based on the inverse 3D DFT method; (c) the reconstructed sea surface field at t = 50 s by summing all the retrieved wave components.
Remotesensing 08 00607 g005
Figure 6. The reconstruction error of the sea surface field. (a) the reconstruction error based on the inverse 3D DFT method; (b) the reconstruction error by summing all the retrieved wave components.
Figure 6. The reconstruction error of the sea surface field. (a) the reconstruction error based on the inverse 3D DFT method; (b) the reconstruction error by summing all the retrieved wave components.
Remotesensing 08 00607 g006
Figure 7. The reconstruction error of the sea surface field. (ac) show the average reconstruction error of the reconstructed sea surface field based on the inverse 3D DFT method at t = 0 s, t = 50 s and t = 100 s, respectively; (df) denote the average reconstruction error of the reconstructed sea field by summing all the retrieved wave components at t = 0 s, t = 50 s and t = 100 s, respectively.
Figure 7. The reconstruction error of the sea surface field. (ac) show the average reconstruction error of the reconstructed sea surface field based on the inverse 3D DFT method at t = 0 s, t = 50 s and t = 100 s, respectively; (df) denote the average reconstruction error of the reconstructed sea field by summing all the retrieved wave components at t = 0 s, t = 50 s and t = 100 s, respectively.
Remotesensing 08 00607 g007
Figure 8. The amplitude of the simulated wave components and the retrieved wave components. (a) the amplitude of the simulated wave components based on the directional wave spectrum; (b) the amplitude of the retrieved wave components based on the novel method with high directional resolution.
Figure 8. The amplitude of the simulated wave components and the retrieved wave components. (a) the amplitude of the simulated wave components based on the directional wave spectrum; (b) the amplitude of the retrieved wave components based on the novel method with high directional resolution.
Remotesensing 08 00607 g008
Figure 9. The sea surface field. (a) the simulated sea surface field; (b) the reconstructed sea surface field based on the 3D DFT method; (c) the reconstructed sea surface field by summing all the retrieved wave components.
Figure 9. The sea surface field. (a) the simulated sea surface field; (b) the reconstructed sea surface field based on the 3D DFT method; (c) the reconstructed sea surface field by summing all the retrieved wave components.
Remotesensing 08 00607 g009
Figure 10. The reconstruction error of the sea surface field with high directional resolution. (a) the reconstruction error based on the inverse 3D DFT method; (b) the reconstructed error by summing all the retrieved wave components.
Figure 10. The reconstruction error of the sea surface field with high directional resolution. (a) the reconstruction error based on the inverse 3D DFT method; (b) the reconstructed error by summing all the retrieved wave components.
Remotesensing 08 00607 g010

Share and Cite

MDPI and ACS Style

Wei, Y.; Zhang, J.-K.; Lu, Z. A Novel Successive Cancellation Method to Retrieve Sea Wave Components from Spatio-Temporal Remote Sensing Image Sequences. Remote Sens. 2016, 8, 607. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8070607

AMA Style

Wei Y, Zhang J-K, Lu Z. A Novel Successive Cancellation Method to Retrieve Sea Wave Components from Spatio-Temporal Remote Sensing Image Sequences. Remote Sensing. 2016; 8(7):607. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8070607

Chicago/Turabian Style

Wei, Yanbo, Jian-Kang Zhang, and Zhizhong Lu. 2016. "A Novel Successive Cancellation Method to Retrieve Sea Wave Components from Spatio-Temporal Remote Sensing Image Sequences" Remote Sensing 8, no. 7: 607. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8070607

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