Sea Clutter Covariance Matrix Estimation and Its Application to Whitening Filter

Article information

J. Electromagn. Eng. Sci. 2021;21(2):134-142
Publication date (electronic) : 2021 April 30
doi : https://doi.org/10.26866/jees.2021.21.2.134
1KOMPSAT Development Division, Korea Aerospace Industries Ltd., Sacheon, Korea
2Department of Electronic Convergence Engineering, Kwangwoon University, Seoul, Korea
3Department of Electronics & Software, Renault Samsung Motors, Yongin, Korea
4Radar Research Institute, Hanwha Systems, Yongin, Korea
5The 3rd R&D Institute, Agency for Defense Development, Daejeon, Korea
*Corresponding Author: Hoongee Yang (e-mail: hgyang@kw.ac.kr)
Received 2020 May 4; Revised 2020 August 21; Accepted 2020 December 15.

Abstract

The accurate estimation of clutter covariance matrix (CCM) is essential in designing a radar detector/filter to suppress sea clutter. This estimation might not be easily accomplished because of the scarcity of valid training vectors adjacent to the range cell under test (CUT). We propose a new CCM estimation algorithm that is derived by modeling time-series clutter returns into a clutter Doppler spectrum in the frequency domain and exploiting mutual independence among spectral components. To justify its excellence over the conventional sample covariance matrix (SCM) algorithm, we design two filters—a maximum signal-to-interference-plus-noise ratio (SINR)-based filter and a whitening filter—that use the estimated CCMs and compare their performance in a numerically simulated sea clutter scenario. Comparisons are made by showing the eigenvector spectra of the estimated CCMs and the frequency responses and outputs of the filters. Moreover, SINRs at the target Doppler bin are examined and compared with a theoretical, analytically derived SINR.

I. Introduction

Detecting a low velocity, small target in sea clutter, especially with low-grazing angle radar, is an important subject in maritime reconnaissance. The existence of sea clutter has been acknowledged as one of the major obstacles that hinders successful target detection. An ideal detector can inherently estimate clutter statistics, such as the clutter covariance matrix (CCM), the clutter power of cell under test (CUT), and range cells adjacent to the CUT, on the fly and apply those parameters to its signal processing component to maintain a constant false alarm rate (CFAR) [1, 2].

It has been said that the stationarity of clutter returns is as short as 100 ms at most [3]. Moreover, because of the heterogeneity of the sea environment, there are not enough training vectors that share the same CCM structure as a primary data vector [4]. These will make it difficult to estimate the CCM and consequently suppress clutter returns.

Many studies have been conducted to investigate the statistical properties of sea clutter, either by analyzing measurement data or by theoretical analysis. The statistical properties include the modeling of Doppler spectra of clutter returns, the spatial and temporal correlation of speckle and texture, clutter amplitude probability density function related with shape parameter and scale parameter, etc. [5, 6]. It turns out that a compound-Gaussian model has been found to be the best fit for sea clutter, especially at a low-grazing angle and high-resolution radar [7].

A detection scheme that works by suppressing clutter returns has been developed through the exploitation of the statistical properties of sea clutter. The generalized likelihood ratio test (GLRT) detector, referred to as a normalized matched filter (NMF), was derived in Gaussian noise, assuming the CCM structure was known [8]. The adaptive GLRT detector structure using the estimated CCM was derived and demonstrated to ensure the CFAR property with respect to the CCM [1].

There have been many studies proposing subspace-based clutter suppression techniques. Mathematical tools such as orthogonal projection (OP), singular value decomposition (SVD), principle eigenvector, etc. were applied to the constructed clutter subspace [9, 10].

In adaptive detectors, the estimation of the CCM has been of paramount importance in determining clutter suppression performance. The CCM in the NMF was replaced by the estimate of the sample covariance matrix (SCM) or by the normalized sample covariance matrix (NSCM) to guarantee the CFAR property with respect to the CCM or to the statistics of the texture or clutter power [4, 11]. In [12], an eigen-analysis-based CCM estimation algorithm called the regularized estimation method was proposed, and it was accessed through the adaptive NMF. Through effectively setting regularization parameters, this method tackled the problem of the scarcity of training vectors. This problem has also been dealt with in the context of the space-time adaptive processing (STAP) algorithm. Knowledge-aided STAP algorithms were proposed in this context, utilizing a priori data for the clutter map rather than many training vectors [13, 14]. To adapt to a more heterogeneous environment, a single data set (SDS) STAP algorithm that does not use any training vectors was proposed [15].

In this paper, we propose a new CCM estimation algorithm that does not require many training vectors and compare it with the conventional SCM-based algorithm, paying attention to the performance variation, which depends on the number of training vectors.

This paper is organized as follows. In Section II-1, we present the modeling of the sea CCM using discrete Fourier transform (DFT). In Section II-2, we propose a new CCM estimation algorithm. To compare it with the conventional SCM, two filters—a maximum signal-to-interference-plus-noise ratio (SINR)-based filter and a whitening filter—using the CCM estimates are introduced and analyzed in Section III-1 and Section III-2. In Section IV, through numerical simulations, we review the performance of the proposed algorithm, showing the two filters’ outputs. Finally, our conclusions are stated in Section V.

II. CCM Modeling and Estimation

1. Sea CCM Modeling

We consider a shipborne radar scenario where a radar makes N consecutive pulse transmissions and corresponding receptions during period T into a radar scene. To focus on sea CCM modeling, the target signature in the radar scene is neglected in this section.

Let cj[n] denote a clutter return of the jth range cell for the nth pulse. For N consecutive transmissions, we obtain the clutter return vector of the jth range cell, {cj[n]}(n = 0, 1, ···, N − 1). Let {Cj[k]}(k=-N2,,N2-1) denote the DFT of {cj[n]}, where k is related to Doppler frequency fk=kNT. According to [5], Cj[k] is modelled as the product of two random variables: texture and speckle. Using the Wiener–Khinchin theorem, the DFT of the autocorrelation function of cj[n] corresponds to the clutter power spectral density (PSD), which is approximately E[|Cj(k)|2]. According to the DFT theory, cj[n] can be expressed as

(1) cj[n]=k=-N2N2-1Cj[k]exp [j(2πN)kn]

Eq. (1) can be rewritten in vector notation as

(2) cj=k=-N2N2-1γj,kNuk

where cj ∈ ℂN×1 denotes cj ≡ [cj[0], …, cj[N−1]]T(superscript T denotes transpose) and for notational simplicity, rj,k is used in place of Cj[k]. Meanwhile, uk in Eq. (2) denotes

(3) uk=1N[ej(2πN)k(0),ej(2πN)k(n),ej(2πN)k(N-1)]T

It is noted that uk ∈ ℂN×1 is, in fact, a Doppler steering vector of fk=kNT with intensity |uk| = 1. The overall interference vector vj can be expressed as

(4) vj=cj+nj

where nj is a complex additive white Gaussian noise vector. Using Eq. (2) and the statistical independence between cj and nj, the theoretical covariance matrix for vj, i.e., Rj ∈ ℂN×N, is expressed as

(5) Rj=E[vjvjH]k=-N2N2-1m=-N2N2-1E[γj,kγj,m*](Nuk)(Num)H+E[njnjH].

When manipulating Eq. (5), we use the orthogonality ukum for km and uknj. According to [16], a collection of random variables—namely, γj,ks—can be assumed to be mutually independent, having zero mean. That is,

(6) E [γj,kγj,m*]=δk,mE[|γj,k|2]

where δk,m is the delta function defined as δk,m = 1 for k = m and δk,m = 0 otherwise.

Using Eq. (6), Eq. (5) is further simplified as

(7) RjNk=-N2N2-1ηj,k2ukukH+σ2I

where ηj,k2=E[|γj,k|2] is the clutter PSD at fk=kNT, σ2 is noise variance, and I is the identity matrix.

Rj in Eq. (7) can be rewritten as the form of the eigen-decomposition, that is,

(8) RjUΛjUH

where U=[u1u2uN] is an N×N matrix whose columns are the eigenvectors of Rj, and Λj = diag{λj,1 λj,2λj,N} is a diagonal matrix in which λj,k = j,k2 + σ2 is an eigenvalue corresponding to eigenvector uk ∈ ℂN×1.

2. Proposed Algorithm for CCM Estimation

Let j denote an estimate for Rj obtained by the SCM, which used to be a primary choice for CCM estimation. j is then expressed as

(9) R˜j=1Llj,l=j-L2j+L2vlvlH

where L is the number of valid training vectors adjacent to the CUT and vl is the interference vector at the lth range cell.

Because of a finite number of training vectors, j cannot be ideally eigen-decomposed, as expressed in Eq. (8).

Let zl be the DFT of vl. zl can then be written as

(10) zl=1NDHvl.

In Eq. (10), D=1N[d-N2,d-N2+1,,dN2-1]T is an N×N Fourier transformation matrix where

(11) dl=[1,ei2πNl,,ei2πNl(N-1)].

Let Cj=[γj,-N2γj,-N2+1γj,N2-1]T, i.e., the DFT of cj. zl is then zl=1N(Cl+nl), where nl=1NDHnl.

Using Eq. (10), j in Eq. (9) can be rewritten as

(12) R˜j=D(NLlzlzlH)DH.

If there were an infinite number of training vectors available, NLlzlzlH in Eq. (12) would be approximately expressed as

(13) NLlzlzlHN[E[|γj,1|2]+σ2E[γj,1γj,N*]E[γj,Nγj,1*]E[|γj,N|2]+σ2]

The off-diagonal terms of the matrix in Eq. (13) should theoretically be zero because of the mutual independence of γj,ks, as mentioned in Eq. (6). Although NLlzlzlH cannot be a diagonal matrix if it is calculated with a finite number of training vectors, we could manipulate it into a diagonal matrix by using the Hadamard product with identity matrix I. Thus, the proposed CCM estimate, hereafter denoted as j, becomes

(14) R^j=D(NLlzlzlHI)DH

where ∘ denotes the Hadamard product.

In summary, the proposed CCM estimation procedures are described as follows: i) obtaining a collection of training vector vls from the range bins adjacent to the CUT, ii) obtaining zls by taking the DFT of vls, iii) calculating (NLlzlzlH)I, and finally, iv) calculating D{(NLlzlzlH)I}DH.

III. Filter Construction Using Estimated CCM

1. Max SINR Filter

We now consider a scenario of the presence of a moving target in the target scene. Since we are concerned with a case with a short observation interval, the target is assumed to have a constant Doppler frequency.

Let the clutter power only reside in [−p1, p2] in the k domain. λj,k in (8) is then expressed as

(15) λj,k={Nηj,k2+σ2k=-p1,-p1+1,p2σ2elsewhere

In this section, we intend to design an adaptive filter that maximizes the output SINR, which hereafter will be referred to simply as the max SINR filter. In this case, for the desired Doppler steering vector sr with unity intensity, optimum weight wj,r becomes [17]

(16) wj,r=Rj-1sr=(k=1Nλj,k-1ukukH)sr

For kr, sr is orthogonal to uk. In this case, wj,r is expressed as

(17) wj,r={1Nηj,r2+σ2urr=-p1,,p21σ2urelsewhere

Let the target Doppler frequency be fq, i.e., k = q. In this case, the received radar return from the jth range cell, xj, can be expressed as

(18) xj=ζjNuq+k=-p1k=p2γj,kNuk+nj

where ζj is the target reflectivity. ζj is assumed constant in the following derivation because the change in target reflectivity for short processing intervals might be negligible. To examine the SINR of the filter output at the target bin, that is, the qth Doppler bin, we let fq be located somewhere within the clutter PSD, i.e., −p1qp2. The max SINR filter output at fr, namely, yj,r, then becomes

(19) yj,r=wj,rHxj={Nζjβr+Nγj,r+urHnjβrr=qNγj,r+urHnjβrrq,-p1rp2urHnjσ2elsewhere

where βr = j,r2 + σ2. The output SINR at the qth Doppler bin becomes

(20) SINR=|Nζjβq|2E{|Nγj,q+uqHnjβq|2}=N|ζj|2Nηj,q2+σ2

In deriving Eq. (20), we use the mutual independence between γj,q and uqHnj. Initially, the filter input SINR at the qth Doppler bin is |ζj|2E[|γj,q|2]=|ζj|2ηj,q2. As expected, when the noise power is negligible, the filter output SINR is approximately the same as the filter input SINR. This result is quite obvious because the mathematical expression of the clutter spectral component at the qth Doppler bin is the same as that of the target return at the qth Doppler bin, except for being scaled by an unknown complex amplitude γj,q. Hence, the filter weight wj,r cannot make differences in the SINRs before and after filtering. Moreover, when noise power is negligible, N does not contribute to SINR improvement, as will be shown in Section IV.

Now, we examine the change in clutter power before and after filtering at the target-free bins. Let the rth Doppler bin be one out of the target-free region, that is, rq and −p1rp2. Employing some manipulations, the clutter power at the rth Doppler bin, which is the mean square value of yj,r, E[|yj,r|2], is expressed as 1Nηj,r2+σ2, whereas the clutter power at the rth Doppler bin before filtering is ηj,r2. This result implies that the bigger ηj,r2 and N are, the larger the clutter suppression ratio becomes.

2. Whitening Filter

In this section, we introduce the whitening filter, whose weight is dictated by wj,r=Rj-12sr(as opposed to the max SINR filter, which uses wj,r=Rj-1sr) and present analytical comparisons of the two filters in terms of the output SINR at the target bin and the interference suppression at the target-free bins.

Following the same procedures leading to Eq. (19) in Section III-1, the whitening filter output at fr, yj,r, is

(21) yj,r=wj,rHxj={Nζjβr+Nγj,r+urHnjβrr=qNγj,r+urHnjβrrq,-p1rp2urHnjσe.w.

where βr = j,r2 + σ2, which is defined the same way it was in Eq. (19). As expected, at the target-free bins where rq, the mean square value of yj,r in (21) becomes approximately 1, which implies that the interference PSD is whitened. The SINR at the qth target Doppler bin is approximately N|ζj|2Nηj,q2+σ2, which is the same as that of the max SINR filter.

Let us compare two filter outputs in terms of the interference suppression at the target-free bins. First, for the clutter-plus-noise bins in which rq and −p1rp2, the E[|yj,r|2]s of the max SINR filter and the whitening filter are approximately 1Nηj,r2 and 1, respectively. Since the initial clutter power at the rth Doppler bin is ηj,r2, their interference suppression ratios are approximately 1Nηj,r4 and 1ηj,r2

, respectively. Thus, depending on whether j,r2 is greater than 1, one filter outperforms the other.

We next examine two filter outputs at the noise-only Doppler bins that do not include clutter components. Employing some calculations, the output noise power E[|yj,r|2]s of the max SINR filter and the whitening filter are approximately 1σ2 and 1, respectively. To compare suppression performance in the noise-only region, we compare the output noise power after equating the target return power at two filter outputs. This can be done by multiplying 1βr by the whitening filter output yj,r in Eq. (21). The calculated noise power ratio is Nηj,r2+σ2σ2, which is always greater than 1. Thus, the whitening filter always suppresses more noise than the max SINR, which will be evident from the simulation results in Section IV.

To conclude, the whitening filter is generally more advantageous than the max SINR filter in that it outperforms noise suppression in the noise-only region while providing comparable results in terms of SINR at the target bin and in the extent of clutter suppression in the clutter plus noise region. In addition, it provides interference whitening in the target-free region, which may simplify the CFAR algorithm that will be applied to the filter output.

IV. Simulations

Before delving into the simulation details, we will briefly describe how the sea clutter returns were simulated. We essentially applied the procedure proposed in [5, 18] and simulated two kinds of clutter returns, VV polarized and HH polarized data, for 64 consecutive pulse transmissions. To simulate VV polarized data, the following steps are taken. Since the clutter spectrum in VV polarization is known to resemble a Gaussian function, we tried to assign the mean Doppler frequency and the spectrum width of the Gaussian function that characterize the function shape. To do this, we assumed that the wind speed is 12.5 m/s, operating frequency is 3 GHz, antenna look angle relative to the wind direction is 0°, relative texture value is 1, and platform-induced Doppler frequency is 0 Hz. With the parameter values assumed, the mean Doppler frequency of the Gaussian function becomes roughly 50 Hz. Its spectrum width was made arbitrarily large such that a moving target signature resides within the clutter spectrum. If a radar platform has motion, the platform-induced Doppler frequency is no longer 0 Hz, and a shift in the mean Doppler frequency is created. We will deal with this nonzero Doppler scenario in a later part of this section. For the chosen function shape, we next obtained the autocorrelation function by taking the inverse DFT (IDFT) of the Gaussian function. The sample values of the autocorrelation function were used as the finite impulse response (FIR) filter weights. As the input sequence to the FIR filter, we used a sequence of complex samples whose in-phase and quadrature components are random variables with zero mean and a variance of one-half. The output sequence of the FIR filter corresponds to the time-series clutter/speckle return data in VV polarization. We used the same texture value for the 64 clutter speckle components, as the interval is relatively short. The dashed line in Fig. 1 shows the initially chosen Gaussian-shaped Doppler spectrum. Following the afore-mentioned procedures, we generated 25 sets of the time-series speckle data cj[n] (n = 0, ···, 63), which is the output of the FIR filter. The solid line in Fig. 1 shows an ensemble average of |Cj[k]|2, which is the clutter PSD.

Fig. 1

Initial function and speckle PSD in VV polarization.

To simulate HH polarized return data, we first determine a Lorentzian function characterized by a small number of mean Doppler frequencies. We arbitrarily chose two mean Doppler frequencies, −150 Hz and 200 Hz. We next obtained the speckle Doppler function by convolving the Lorentzian function with an arbitrarily chosen Gaussian function. Since HH polarization returns are related to scatterers moving faster than Bragg scatterers, the typical function shape of the Doppler spectrum cannot be specified in such a way as to associate it with sea state and radar parameters. Thus, we do not enumerate the associated parameter values. For a given Doppler function, the rest of the procedure to simulate the time-series speckle return data is the same as before. The dashed line and the solid line in Fig. 2, respectively, show the initially chosen Doppler function and the PSD of the simulated HH polarized speckle data.

Fig. 2

Initial function and speckle PSD in HH polarization.

The first scenario concerns a short interval case of 64 radar returns in VV and HH polarization. A shipborne radar with a single transceiver operating at 3 GHz sends a sequence of pulses at 1 kHz of pulse repetition frequency (PRF). Within the target scene, there exists a moving target whose arbitrarily chosen Doppler frequency is 109 Hz. Considering the short interval, we assumed the target reflectivity was constant.

To calculate j using the SCM shown in Eq. (9) with the VV polarized data, we used 25 training vectors adjacent to the CUT. We next decomposed j into eigenvectors and eigenvalues. For all eigenvectors, we took the DFT and obtained their spectra. Fig. 3(a) shows the spectra arranged in such a manner that the spectrum of an eigenvector associated with a larger eigenvalue lay to the left of the x axis. We next applied the proposed algorithm to obtain j. Fig. 3(b) shows the plots of the spectra of the eigenvectors in j, which are arranged as in Fig. 3(a). We note that from the spectra shown in Fig. 3(b), the eigenvectors look to have a form of a Doppler steering vector of a single frequency, which is particularly apparent in the eigenvectors associated with larger eigenvalues. In order to determine the results that could be achieved if there were a sufficient number of training vectors, we arbitrarily generated 250 training vectors, calculated the j using the SCM, and finally obtained the eigenvector spectra shown in Fig. 3(c). More eigenvectors with larger eigenvalues appear to have the form of Doppler steering vectors of a single frequency. Fig. 3(d) is the plot of the eigen vector spectra obtained with the HH polarized data by applying the proposed method. Two dominant Doppler frequencies can be observed in the spectra of the eigenvectors associated with the large eigenvalues.

Fig. 3

Spectra of eigenvectors. (a) CCM estimate by SCM (25 training vectors; VV). (b) CCM estimate by proposed algorithm (VV). (c) CCM estimate by SCM (250 training vectors; VV). (d) CCM estimate by proposed algorithm (HH).

We next compared the frequency responses of the four filters constructed with CCM estimates by the VV polarized data. For convenience, three of these filters—whose weights are j1sr, j1sr, and R^j-12sr—obtained with 25 training vectors are called filter A, filter B, and filter C, respectively. The filter with weight j1sr using 250 training vectors is called filter D. Fig. 4(a)–(d) show the frequency responses of filters A–D, respectively. The frequency responses of filter B and filter C appeared smoother around the clutter-dominant region than those of filter A, which implies that filter A might not effectively suppress clutter components. As expected, filter D, which was obtained with many training vectors, exhibits better performance than filter A.

Fig. 4

(a) Frequency response of filter A. (b) Frequency response of filter B. (c) Frequency response of filter C. (d) Frequency response of filter D.

We next added a target signature to the generated interference, VV polarized data plus additive noise, and examined the filter outputs. The target return power was set to be small enough that the input signal-to-clutter ratio (SCR) and input clutter-to-noise ratio (CNR) were −27 dB and 20 dB, respectively. As shown in Fig. 5(a), the target peak could not be observed at the target Doppler bin of filter A. The target peak could be observed at the right Doppler frequency in Fig. 5(b), (c), and (d), which are the outputs of filter B, filter C, and filter D, respectively. The target peak shown in Fig. 5(c) looked more dominant because the interference in the target-free region was whitened.

Fig. 5

Filter outputs in VV polarization. (a) Output of filter A. (b) Output of filter B. (c) Output of filter C. (d) Output of filter D.

Fig. 6(a)–(d) show the outputs of filter A, filter B, filter C, and filter D, respectively, for the HH polarized data scenario. We may easily observe the same results as we did in the VV polarization case.

Fig. 6

Filter outputs in HH polarization. (a) Output of filter A. (b) Output of filter B. (c) Output of filter C. (d) Output of filter D.

As mentioned earlier, the platform-induced Doppler frequency tends to shift the mean Doppler frequency of the clutter spectrum [18]. To reflect this aspect, we examined the filter outputs, moving around the moving target Doppler frequency while the clutter spectrum remained fixed. From a radar detector standpoint, this is basically equivalent to shifting the clutter mean Doppler frequency. Fig. 7(a) and (b) denote three SINRs obtained from the VV and HH polarized data, respectively.

Fig. 7

Three SINRs in (a) VV polarization data and (b) HH polarization data.

The three SINR curves shown in Fig. 7 are the theoretical SINR curve and two other SINR curves obtained with the CCM estimates by the proposed and the conventional/SCM method, respectively. Calculation of SINRs was made at arbitrarily selected different Doppler frequencies in either the clutter plus noise region or the noise-only region of the clutter spectrum. As the theoretical SINR was obtained by Eq. (20), that is, by using a true CCM, it must put a ceiling on the SINR. For both the VV and the HH polarized data, the SINRs by the proposed method are closer to the theoretical SINRs than are the SINRs by the conventional method. This is quite obvious in the clutter plus noise region (see Figs. 1 and 2), which implies that the filter designed with the CCM from the proposed method effectively suppresses clutter. However, no noticeable difference among the three SINRs in the noise-only region can be observed, as the filter performance has no effect in this region.

V. Conclusion

In this paper, we proposed a CCM estimation algorithm in a sea clutter environment. To develop the algorithm, we modified the time-series clutter returns into the frequency domain clutter spectrum and exploited the mutual independence among the complex-valued Doppler spectral components. To compare the proposed algorithm with the SCM algorithm, two filters—the max SINR filter and the whitening filter—that require the CCM were used. We first analytically derived the SINR at the target bin and the clutter suppression ratio in clutter and noise-only regions. Next, we numerically simulated sea clutter, applied the conventional and proposed estimation algorithms, and examined the outputs of the filters designed with the estimated CCMs. Comparisons were made by showing the eignvector spectra of the estimated CCMs and the frequency responses and outputs of the filters. Moreover, SINRs at target Doppler bins were examined with a theoretical, analytically derived SINR.

The filter using the proposed CCM estimate, regardless of the filter type, always showed better performance in suppressing clutter than the filter using the SCM estimate. Moreover, the whitening filter was more advantageous, especially in the noise-only region.

Acknowledgments

This work was supported by a Research Grant from Kwangwoon University in 2020 and by the Agency for Defense Development (No. UD170020FD).

References

1. Kelly EJ. An adaptive detection algorithm. IEEE Transactions on Aerospace and Electronic Systems AES-22(2):115–127. 1986;
2. Conte E, Longo M, Lops M, Ullo SL. Radar detection of signals with unknown parameters in K-distributed clutter. IEEE Proceedings F (Radar and Signal Processing) 138(2):131–138. 1991;
3. Shui PL, Shi YL. Subband ANMF detection of moving targets in sea clutter. IEEE Transactions on Aerospace and Electronic Systems 48(4):3578–3593. 2012;
4. Greco M, Stinco P, Gini F. Impact of sea clutter nonstationarity on disturbance covariance matrix estimation and CFAR detector performance. IEEE Transactions on Aerospace and Electronic Systems 46(3):1502–1513. 2010;
5. Watts S. Modeling and simulation of coherent sea clutter. IEEE Transactions on Aerospace and Electronic Systems 48(4):3303–3317. 2012;
6. Greco M, Gini F, Rangaswamy M. Statistical analysis of measured polarimetric clutter data at different range resolutions. IEEE Proceedings–Radar, Sonar and Navigation 153(6):473–481. 2006;
7. Conte E, Bisceglie MD, Galdi C, Ricci G. A procedure for measuring the coherence length of the sea texture. IEEE Transactions on Instrumentation and Measurement 46(4):836–841. 1997;
8. Conte E, Lops M, Ricci G. Asymptotically optimum radar detection in compound-Gaussian clutter. IEEE Transactions on Aerospace and Electronic Systems 31(2):617–625. 1995;
9. Poon MWY, Khan RH, Le-Ngoc S. A singular value decomposition (SVD) based method for suppressing ocean clutter in high frequency radar. IEEE Transactions on Signal Processing 41(3):1421–1425. 1993;
10. Sira SP, Cochran D, Papandreou-Suppappola A, Morrell D, Moran W, Howard S. A subspace-based approach to sea clutter suppression for improved target detection. In : Proceedings of 2006 40th Asilomar Conference on Signals, Systems and Computers. Pacific Grove, CA; 2006; p. 752–756.
11. Conte E, Lops M, Ricci G. Adaptive detection schemes in compound-Gaussian clutter. IEEE Transactions on Aerospace and Electronic Systems 34(4):1058–1069. 1998;
12. Kammoun A, Couillet R, Pascal F. Optimal design of the adaptive normalized matched filter detector using regularized Tyler estimators. IEEE Transactions on Aerospace and Electronic Systems 54(2):755–769. 2018;
13. Melvin W, Wicks M, Antonik P, Salama Y, Li P, Schuman H. Knowledge-based space-time adaptive processing for airborne early warning radar. IEEE Aerospace and Electronic Systems Magazine 13(4):37–42. 1998;
14. Jeon H, Chung Y, Chung W, Kim J, Yang H. Clutter covariance matrix estimation using weight vectors in knowledge-aided STAP. Electronics Letters 53(8):560–562. 2017;
15. Chen ZY, Zhou Y, Zhang LR, Li Q, Gu YB. A robust single data set-STAP algorithm. In : Proceedings of IET International Radar Conference. Hangzhou, China; 2015;
16. Theodoridis S, Chellappa R. Academic Press Library in Signal Processing: Communications and Radar Signal Processing New York, NY: Academic Press; 2014.
17. Guerci JR. Space-Time Adaptive Processing for Radar Norwood, MA: Artech House; 2003.
18. Kemkemian S, Lupinski L, Corretja V, Cottron R, Watts S. Performance assessment of multi-channel radars using simulated sea clutter. In : Proceedings of 2015 IEEE Radar Conference (RadarCon); 2015; p. 1015–1020.

Biography

Sanghyun Choi received B.S. and M.S. degrees in electric wave engineering from Kwangwoon University, Seoul, Korea, in 2017 and 2019, respectively. He currently works at Korea Aerospace Industries Ltd., Sacheon, Korea, where he is engaging in research on and development of satellite electrical systems. His research interests include the generation and suppression of sea clutter and moving target imaging in a MIMO radar environment.

Hoongee Yang received a B.S. degree in electronic engineering from Yonsei University, Seoul, Korea, in 1985 and M.S. and Ph.D. degrees in electrical and computer engineering from the State University of New York at Buffalo, Amherst, NY, in 1987 and 1992, respectively. Since 1993, he has been with the Department of Electronic Convergence Engineering, Kwangwoon University, Seoul, Korea. His current research interests are radar signal processing, array signal processing, and SAR/ISAR imaging.

Jimin Song received B.S. and M.S. degrees in electric wave engineering from Kwangwoon University, Seoul, Korea, in 2017 and 2019, respectively. He is currently working in the Electronics & Software Development Department of Renault Samsung Motors, Yongin, Korea, where he is being engaged in the development of a control logic for safety & seat systems. His current interests are radar imaging techniques and sensor fusion using radar and vision data.

Hyeonmu Jeon received a B.S. degree in electronic engineering and a Ph.D. degree in radio science and engineering from Kwangwoon University, Seoul, Korea, in 2017. Since 2017, he has been with the Radar Research Institute at Hanwha Systems, where he is currently a Senior Engineer. His main research interests are array signal processing, space time adaptive processing, and synthetic aperture radar.

Jongmann Kim received a B.S. degree in electronic and electrical engineering from Kyungpook National University, Daegu, Korea, in 2002 and an M.S. degree in electronic and electrical engineering from Pohang University of Science and Technology (POSTECH), Pohang, Korea, in 2004. Since 2005, he has been with the Agency for Defense Development (ADD), Daejeon, Korea, where he is currently a Principle Researcher. His current research interests include radar signals and data processing.

Yongseek Chung received his B.S., M.S., and Ph.D. degrees in electrical engineering from Seoul National University, Seoul, Korea, in 1989, 1991, and 2000, respectively. From 1991 to 1996, he was with the Living System Laboratory at LG Electronics. From 1998 to 2000, he was a Teaching Assistant in electrical engineering at Seoul National University. From 2001 to 2002, he was with Syracuse University, Syracuse, NY. From 2003 to 2005, he has been a faculty member with the Department of Communication Engineering, Myongji University, Kyunggi, Korea. He is currently a professor in the Department of Electronic Convergence Engineering, Kwangwoon University, Seoul, Korea. His current interests are radar signal processing and microwave imaging techniques.

Article information Continued

Fig. 1

Initial function and speckle PSD in VV polarization.

Fig. 2

Initial function and speckle PSD in HH polarization.

Fig. 3

Spectra of eigenvectors. (a) CCM estimate by SCM (25 training vectors; VV). (b) CCM estimate by proposed algorithm (VV). (c) CCM estimate by SCM (250 training vectors; VV). (d) CCM estimate by proposed algorithm (HH).

Fig. 4

(a) Frequency response of filter A. (b) Frequency response of filter B. (c) Frequency response of filter C. (d) Frequency response of filter D.

Fig. 5

Filter outputs in VV polarization. (a) Output of filter A. (b) Output of filter B. (c) Output of filter C. (d) Output of filter D.

Fig. 6

Filter outputs in HH polarization. (a) Output of filter A. (b) Output of filter B. (c) Output of filter C. (d) Output of filter D.

Fig. 7

Three SINRs in (a) VV polarization data and (b) HH polarization data.