Abstract

A modified space-time adaptive processing (STAP) estimator is described in this paper. The estimator combines the incremental multiparameter (IMP) algorithm and the existing beam-space preprocessing techniques yielding a computationally cheap algorithm for the superresolution of multiple signals. It is a potential technique for the remote sensing of the ocean currents from the broadened first-order Bragg sea echo spectrum of shipborne high-frequency surface wave radar (HFSWR). Some simulation results and real-data analysis are shown to validate the proposed algorithm.

1. Introduction

The measurement of the near-surface currents is a very difficult task by using conventional methods, especially under some harsh sea conditions. Many advanced marine measurement instruments, such as drifting buoys [1] and acoustic current meters [2], have been used to collect sea state information. However, it would be very expensive to collect and interpret data from these devices for the sparse spatial sampling provided. Therefore, it is virtually impossible to form the current maps over a widespread ocean surface timely and accurately with these conventional meters.

In recent years, HFSWR has already become a powerful remote-sensing tool which receives increasing attention from oceanographers and research groups for its good ability to determine the large-scale sea state under all weather conditions [3, 4]. The measurement principle of HFSWR mainly depends on the Bragg resonant scattering theory and Doppler frequency effect theory in the sea echo spectrum of HFSWR. In the absence of ocean currents, the first-order Bragg lines would appear symmetrically above and below the zero Doppler frequency, which are caused by the ocean waves with precisely one-half wavelength of the radar moving towards and away from the radar. In practical application, some displacements will happen in the first-order Bragg lines since the near-surface currents always exist on the ocean [5].

As an extension of shore-based HFSWR, shipborne HFSWR not only inherits all the advantages of shore-based HFSWR but also shows some outstanding features, such as flexibility and mobility. However, some new problems emerge as the radar on board is a moving ship. One of the worst problems is that the first-order Bragg lines have been broadened into two pass bands (when the speed of ship is slow) or one low pass band (when the ship is sailing at a high speed) in the first-order Bragg sea echo spectrum of shipborne HFSWR [6]. For these cases, it will be hard to determine the ocean currents from the broadened sea echo spectrum. Furthermore, the moving ship yields different Doppler shifts to the different azimuth sea echoes. Thus, there exists a certain space-time coupling relation in the received sea echo spectrum of shipborne HFSWR.

The novel space-time IMP algorithm was proposed by Clarke and Spence [7] which was based on one-dimensional IMP algorithm and modified to detect and estimate multiple signals from the conventional beamwidth and (or) Doppler resolution bin. Although space-time IMP estimator effectively improves the robustness of the detection and estimation of multiple signals, the computational load of two-dimensional search process is too heavy for real-time application. Lately, Chadwick [8] used the eigendecomposition method instead of the full-search process to reduce the heavy computational burden in his proposed polarisation-sensitive IMP algorithm. However, for surface wave radar, this method becomes invalid since the returns of HFSWR are mostly vertical polarization sensitive components.

Many researchers, such as Shaw and Wilkins [9], used beam-space preprocessing technique to reduce the computational load and improve robust performance of high-resolution DOA (direction-of-arrival) estimation algorithms. Recently, Hassanien et al. [10, 11], proposed a new concept about beam-space preprocessing algorithm which was able to suppress the out-of-sector interferences accompanied with the updated array data. This technique was proven to be more robust than the aforementioned beam-space methods.

In this paper, we combine the space-time IMP algorithm and the adaptive beam-space preprocessing technique yielding a computationally cheap space-time adaptive estimator for detection, estimation and super-resolution of multiple signals. The proposed algorithm is validated by simulation results as well as experimental examples.

This paper is organized as follows. First we introduce the signal model. In the following section, the proposed algorithm and some simulation results are presented. The measurement of the near-surface currents of the ocean by shipborne HFSWR and some real-data analysis are presented in Section 4. The final part is the study conclusion.

2. Signal Model

Consider a uniform linear array (ULA) with omnidirectional antennas and the antenna spacing . If there are sources that have been received from a far field with different relative delays and attenuations. The received data is then given by [7] as follows: where denotes the data received from the th sensor at th sampling time, denotes the transpose operation, is the array manifold matrix, is the steering vector points to the direction , denotes wavelength of the radar, is a () diagonal matrix containing the signal magnitudes, is a () matrix comprising the normalized source waveforms, is the normalized source waveform, is the frequency of th source, denotes the Hermitian transpose operation, and is the () matrix comprising the zero mean and variance Gaussian noise.

The () covariance matrix of the received data is given by where is the statistical expectation operator, is the () source covariance matrix, and is the identity matrix.

3. Modified STAP Estimator

3.1. Space-Time IMP

Space-time IMP is a two-dimensional maximum likelihood method which uses a set of space-time calibration response vectors to match with the received data. Thus, the primary objective of space-time IMP is to maximize the “signal plus noise” to “expected noise” power ratio (SNR). If the maximum output power is over the threshold, then a target is detected and the corresponding space-time calibration response vector is recorded. In order to reduce the sidelobe leakage of the detected targets and improve the detection and estimation of the potential signals in the residual data, the detected targets are removed from the original data through an orthogonal subspace projection matrix before each iterative stage.

According to the definition of SNR in space-time IMP, we have the following [7]: where is the () space-time calibration response matrix, is the vectorization operation, and is an () orthogonal projection matrix, which is given by where denotes a sum of space-time calibration response vectors corresponding to the detected signals; that is, projects the received data into a subspace orthogonal to the detected signals.

3.2. Adaptive Beam-Space Preprocessing

The adaptive beam-space preprocessing technique was first mentioned in [10] which used the updated data for adaptive suppression of out-of-sector interferences. This technique had been shown to be more robust than the aforementioned beam-space methods.

The primary objective of the data-adaptive beam-space preprocessing is to solve the optimal beam-space matrix design problem through minimizing the output power of the transformed data, which can be expressed as follows [12]: where is the trace of a matrix, is the () beam-space matrix, is the beam-space dimension, is the vector 2-norm, denotes all out-of-sector directions which are divided into angular grids, and are the angles corresponding to the in-of-sector and out-of-sector directions, and is the stopband attenuation parameter which should meet the requirement [12]

After the beam-space transformation, the array steering vector matrix and manifold matrix have already transformed into

Then, the () covariance matrix in beam-space should be rewritten as

Obviously, the dimension of the matrix is lower than that of . This fact is exploited in all beam-space-based methods to reduce computational load compared with the element-space algorithms [10].

3.3. Space-Time Adaptive Beam-Space IMP

In this section, we show how the conventional space-time estimator combines with the adaptive beam-space preprocessing technique to present a computationally cheap space-time adaptive beam-space IMP estimator.

Following the discussion above, the two-dimensional discriminants shown in (5) should be modified as where and has been reduced to a matrix in beam-space domain, which is given by where denotes a sum of space-time calibration response vectors corresponding to the detected signals in beam-space domain.

3.4. Threshold Setting

How to select an appropriate threshold to terminate the iterative process in IMP algorithm is very important. Theoretically, when all “significant peaks” have been detected and cleared out from the received data, there is only a completely flat plane in the residual scan [13]. However, it is impossible to accurately estimate the noise statistics from the limited received data. Furthermore, the definition of “significant peak” in IMP algorithm has not been clearly reported.

In the paper, we use a double-threshold setting method to ensure that the iterative process is halted timely. First, we check two successive scans before the next iteration. If the difference between the two scans is comparable to that of the “expected noise” level, that is, no “significant peak” appears during the last scan, then the iterative process should be halted. Besides, if the difference between two successive estimations has reached the preset threshold, which suggests that the iterations are estimating the same target, then the iterative process should be halted as well.

3.5. Simulation Results

Several simulation results are shown in this section to test the performance of the modified algorithm through comparing it with several conventional algorithms.

During the simulations, we assume that the radar works at , which contains an ULA with omnidirectional sensors and the elements are spaced one-half wavelength apart. The half power beamwidth is approximately . The number of snapshots and the beam-space dimension are chosen for our simulations. The adaptive beam-space matrix has been solved using the cvx optimization MATLAB toolbox. Since the minimum value in (9) is , we take the parameter for (8). Furthermore, the two simulation targets and are also selected in the simulations.

To define a successful experiment, we use the criterion mentioned in [14] if where and   () are, respectively, the estimated and truth values, then the two signals are successfully resolved.

Figures 1 and 2 illustrate the probability of source resolution and their root mean square error (RMSE) versus SNR in Doppler domain, respectively. The conventional space-time IMP, space-time beam-space IMP which combines space-time IMP with discrete fourier transform (DFT) matrix beam-space processing technique [15], and 64 points and 256 points FFT results are used for the comparison in the figures. We find that the beam-space-based algorithms show better resolution and smaller RMSE than the other algorithms in resolving the two simulation targets. Thus, it is reasonable to conclude that the beam-space-based methods require less observation time but maintain high Doppler accuracy compared with the conventional algorithms. By the way, all the simulation results shown in this section have averaged over 1000 independent Monte Carlo experiments.

Figures 3 and 4 are the probability of source resolution and their RMSE versus SNR in azimuth domain, respectively. According to these figures, the beam-space-based algorithms show better performance than the other algorithms. Comparing these methods, the space-time adaptive beam-space IMP algorithm shows the highest robust and lowest RMSE in resolving the two simulation targets. Thus, it is reasonable to conclude that the proposed algorithm requires smaller antenna array and lower SNR threshold for detection and estimation of multiple signals when compared with the conventional DOA algorithms.

4. Shipborne HFSWR

4.1. Space-Time Coupling Relation

In [6], Xie et al. had proven that the first-order Bragg spectrums were broadened along the azimuthal directions from the real-data analysis. The authors concluded that there was a space-time coupling relation existed in the first-order Bragg sea echo spectrum of shipborne HFSWR.

Assuming that both the transmitting and the receiving antennas of HFSWR are mounted on a ship which is moving in the positive direction of the x-axis at a constant speed , as shown in Figure 5.

In the absence of ocean current, the space-time coupling relation in the first-order Bragg sea echo spectrum of shipborne HFSWR can be expressed as follows [16]: where is the azimuth direction, are the first-order Bragg frequencies, the positive and negative signs are, respectively, the Bragg waves moving towards and away from the radar.

4.2. Current Measurement

In the presence of ocean currents, the first-order Bragg lines in (16) are shifted from the theoretical positions. The displacements are proportional to the radial current velocities. Thus, (16) should be rewritten as [17]

As shown in (17), the first-order Bragg lines are related to the azimuth directions as well as the speed of ship and currents. Therefore, the first-order Bragg peaks are not only displaced, but also broadened into two pass bands (for slow ship speed case) in the first-order Bragg sea echo spectrum of shipborne HFSWR.

4.3. Real-Data Analysis

The real data used in this paper was received from the shipborne HFSWR experiments conducted on the Yellow Sea of China on September 8, 1998 [6]. Figures 6 and 7 show the space-time coupling relation in the sixth range bin of the real-data file 1128 (containing 7 channels 256 samples 32 range bins and the ship speed was about 4.8 m/s), which is processed through the DFT plus DBF (Digital Beamforming) cascade processing and the proposed algorithm. As shown in the figures, the first-order Bragg lines are broadened along the azimuth directions, which tally well with the theoretical lines in (16). The displacements may be caused by ocean currents or interferences.

Table 1 is an example of the DOA and Doppler estimations. The modified STAP estimator is used to process the above-mentioned real data within the section between the azimuth and . Since there was no information about the sea state recorded during the experiment, we here assume that the detected targets near the theoretical first-order Bragg lines are considered as ocean currents and their displacements are proportional to their radial velocities. Based on this assumption, five currents have been detected and estimated. All of them are very close to the theoretical first-order Bragg frequencies and their corresponding radial velocities are calculated in the table.

Table 2 is another example for better understanding the robustness of the proposed algorithm, where we add a simulation target to the real data used previously, as shown in Figure 8. The added simulation target is detected correctly and the estimations of DOA and Doppler frequency are within and of the true signal position.

5. Conclusion

In this paper, we have introduced a modified space-time adaptive processing estimator that can be used for the detection, estimation, and superresolution of multiple signals. The method combines the conventional IMP method and the existing adaptive beam-space preprocessing techniques yielding a computationally cheap algorithm for estimating the near-surface currents of the ocean from the broadening of the first-order Bragg sea echo spectrum of shipborne HFSWR. The proposed algorithm is validated by simulation results as well as experimental examples.

Acknowledgment

This work is supported by the State Key Program of National Natural Science of China (Grant no. 61132005).