An Improved Clutter Suppression Method for Weather Radars Using Multiple Pulse Repetition Time Technique
This paper describes the implementation of an improved clutter suppression method for the multiple pulse repetition time (PRT) technique based on simulated radar data. The suppression method is constructed using maximum likelihood methodology in time domain and is called parametric time domain method (PTDM). The procedure relies on the assumption that precipitation and clutter signal spectra follow a Gaussian functional form. The multiple interleaved pulse repetition frequencies (PRFs) that are used in this work are set to four PRFs (952, 833, 667, and 513 Hz). Based on radar simulation, it is shown that the new method can provide accurate retrieval of Doppler velocity even in the case of strong clutter contamination. The obtained velocity is nearly unbiased for all the range of Nyquist velocity interval. Also, the performance of the method is illustrated on simulated radar data for plan position indicator (PPI) scan. Compared with staggered 2-PRT transmission schemes with PTDM, the proposed method presents better estimation accuracy under certain clutter situations.
The velocity and range ambiguity  is still one of the major limitations for Doppler weather radar observation. Given a radar wavelength , the relationship between the maximum unambiguous range and maximum unambiguous velocity can be presented as , where is the speed of light. With a uniform pulse repetition frequency (PRF) observation, can be shown as . That means large maximum unambiguous velocity can be obtained by increasing the PRF. Along with the increase of maximum unambiguous velocity comes a decrease in maximum range , and vice versa.
Various techniques have been proposed to overcome this problem [2–5], like dual-pulse repetition time (PRT) and staggered PRT [3, 6]. The dual-PRT and staggered dual-PRT techniques utilize two PRTs in an pulses transmission sequence. The dual-PRT method transmits two blocks of uniform PRT pulses with each block at a different PRT . The optimal velocity estimates can be recovered from the joint analysis of the two Doppler measurements. This method is suitable for the low antenna rotation rate, because the radial velocity remains the same during the pulse intervals. The staggered PRT method is not restricted by this condition. The PRT alternates between two values . The interval of the maximum unambiguous velocity was well extended and successively adopted for the Next Generation Weather Radar (NEXRAD) network . In order to improve the dealiasing success rate and obtain larger Nyquist velocity, the multiple PRT (-PRT) technique which utilizes multiple interleaved PRTs was proposed. Tabary et al.  presented a triple-PRT scheme implemented on the C-band French Radar Network. This technique works well when spectral width is smaller than 3 m/s and SNR is high. A simultaneous multiple pulse repetition frequency (SMPRT) code  was introduced and evaluated. This technique can provide enough information to produce a high-resolution measured spectrum for each range gate. Tahanout et al.  proposed an optimal 9-PRT scheme for the suitable shape reproduction of the power spectrum of the radar signal.
These techniques which used nonuniform time series make filtering of clutter from the radar signal more complicated. Standard clutter filters cannot be applied directly to the nonuniform sampling scheme. Banjanin and Zrnic  developed a scheme which consists of two filters that operate sequentially; the overall filter is time-varying with periodically changing coefficients. Cho and Chornoboy  introduced a finite impulse response time-varying filter that is applied to -PRT transmission. Siggia and Passarelli  proposed a Gaussian model adaptive processing (GMAP) algorithm that filters clutter and recovers weather spectrum after notching the spectrum around zero Doppler. Nguyen et al.  and Moisseev et al.  presented a parametric time domain method (PTDM) for staggered PRT observations, which can accurately estimate the spectral moments of precipitation echoes even with strong clutter.
In this work, we extend PTDM to -PRT observations. This paper is organized as follows. In Section 2, staggered 2-PRT and -PRT observation techniques are described. The PTDM for clutter suppression and spectral moments estimation is then presented in Section 3. Based on radar data simulations, Section 4 focuses on performance analysis of the proposed method. The error analyses of the proposed method and the staggered 2-PRT technique with PTDM for different noise and clutter scenarios are presented. In order to test the performance of the proposed method, different sample schemes techniques are used for plan position indicator (PPI) observations for comparison. In Section 5, discussions and conclusions are provided.
2. -PRT Techniques
2.1. Uniform PRT
Assuming that is the PRT of the traditional uniform PRT technique, the autocorrelation calculated from return signal time series at lag is . The maximum unambiguous velocity is . The estimated velocity can be deduced from  aswhere is argument function.
2.2. Staggered 2-PRT
The staggered PRT method applies a pulse transmission sequence that changes intervals between and . For , the maximum unambiguous velocity will be extended as
The autocorrelations calculated from return signal time series at lags and are and . The Doppler velocities and can be deduced from and , respectively, with (1).
Two approaches have been introduced to deal with the velocity aliasing. One algorithm uses the ratio of autocorrelations , given by
If the argument of exceeds the range of , the estimated velocity will be ambiguous. The other technique utilizes the maximum unambiguous velocity to revise the aliasing velocity. Suppose and are the maximum unambiguous velocity for uniform PRTs and . The dealiased velocity can be shown aswhere and are the correct Nyquist interval number and and are estimate errors.
Mostly, -PRT technique  can be regarded as an extension of the staggered 2-PRT technique. This technique chooses a transmission sequence which has several blocks, and every block is formed of different PRTs. Its transmitting scheme is . Similar to the staggered 2-PRT method, the maximum unambiguous velocity can be shown as
The relationship between the measured true velocity and the estimated velocities is given by
Finding the appropriate integers which minimize the error is the way for the measured true velocity retrieval.
3. Parametric Time Domain Method (PTDM)
Bringi and Chandrasekar  have shown that the real and imaginary parts of radar signals follow zero mean normal distribution. Suppose is samples of received radar signals in one radar resolution volume with sample time . Hence, the complex vector can be expressed as the sum of two real Gaussian vectors representing in-phase and quadrature components, respectively, like
The multivariate density function of  can be represented aswhere is the covariance matrix and is the sample covariance matrix. is the trace function. is the determinant of .
Assuming that Doppler spectra of clutter and precipitation obey Gaussian distribution, the observed spectrum of Doppler velocity can be expressed as the summation of independent spectra coming from precipitation , clutter , and white noise , aswhere and are the precipitation signal and clutter power, and are the precipitation and clutter spectrum width, respectively, is the noise power, is mean velocity, and is sample time. Then, the covariance function can be obtained by fast Fourier transform (FFT) aswhere is the vector of unknown parameters, . is the temporal variable.
For -PRT transmission, the element in the covariance matrix can be written aswhereHere, is integral function. is the remainder of a division.
Using the multivariate density function (8), the log-likelihood estimation of is as follows:The parameter vector can be estimated by minimizing .
4. Analysis of the Performance
4.1. Simulation of Radar Signal
The main task of this work is to discuss the accuracy of estimated velocity under different clutter conditions. So, different values of mean velocity, signal-to-noise ratio (SNR), and clutter-to-signal ratio (CSR) are chosen to generate the observed Doppler spectra. We select a periodic scheme with 4 different PRTs: 1050, 1200, 1500, and 1950 . Their greatest common factor is . The ratio of PRTs is 7/8/10/13. The corresponding PRFs can be written as 952, 833, 667, and 513 Hz. The parameters of the simulated X-band radar signal are shown in Table 1.
4.2. PTDM Implementation
Using the log-likelihood function (13) to retrieve the precipitation and clutter spectral moments, the convergence of a local minimum will lead to wrong estimations in some cases. In order to reduce the estimate errors, we should make sure that the results of minimum (13) come from the convergence of global minimum. According to Moisseev et al. , this problem can be solved by selecting the appropriate seed value, especially for mean velocity seed value selection which has the biggest effect on the convergence of the optimization procedure. Similar to Moisseev et al. , the search range of unknown parameters is given in Table 2. is the total signal power, and represents the system noise floor. Using the Nelder-Mead simplex method [17, 18], the estimated values of unknown parameters are retrieved. Figure 1 shows a scatter diagram of estimated velocities for 4-PRT technique using PTDM, when CSR = 20 dB. It can be revealed that the estimated velocities are accurately retrieved using this method. The mean velocities almost coincide with true velocities. The standard deviations and biases are presented in Figure 2. The results show that the standard deviations are less than 0.5 m/s and the biases are near zero without increase.
The noise and clutter environment is complex and variable in radar detection. The simulation can be carried out through two cases. The first one describes the conditions near the radar. The ground clutter echoes have various intensities. The other one considers the situations far away from the radar where ground clutter can be neglected. The simulated data are collected in the Nyquist velocity interval.
In the first type of simulation, the CSR changes from −20 dB to 40 dB. SNR is fixed to +40 dB during the computation. Figure 3 shows the curves of standard deviation and bias of velocity estimates versus CSR, when applying 2-PRT and -PRT schemes. The clutter suppression method used in both transmission schemes is PTDM. The biases increase with CSR at a similar rate for the two PRT schemes. The absolute values of biases do not exceed 0.5 m/s. For the range from −20 dB to 20 dB, the standard deviation are smaller than 1 m/s and obviously the 4-PRT schemes provide better estimation accuracy. The estimation error increases with CSR, especially from 20 dB to 30 dB. When CSR > 30 dB, the values tend to be almost steady near 1.8 m/s.
The second simulation is implemented for 0 dB ≤ SNR ≤ 40 dB and CSR = −20 dB. In Figure 4, the biases and standard deviations are shown for the PTDM velocity estimates. In the 15 dB, 40 dB interval, the measured velocities for both PRT schemes are sufficiently accurate. Then, the bias and standard deviation increase with the decrease of SNR and reach the maximum near SNR = 0 dB. For the range from 5 dB to 40 dB, the measured velocities are acceptable for the biases close to zero and the standard deviations are less than 0.6 m/s. The performance of 4-PRT schemes has better precision especially when the noise is strong (from 0 dB to 15 dB).
In order to test the effectiveness of the proposed method, the and data for plan position indicator (PPI) scan were simulated. Choosing a range-azimuth space of 30 km × 60°, the resolution is 150 m and 1°. The radial velocity of each cell in the background is set to one constant. The fan-shape range from the center to 12 km is defined as ground clutter with zero Doppler velocity and spectral width m/s. At the range from 8.5 km to 25 km, a circle region is set to simulate the wind shear in which the velocity field is unstable and changes rapidly in velocity magnitude and direction. The maximum velocity in wind shear field is 10 m/s and spectral width is equal to 4 m/s. From Figure 4, we have already known that the velocities estimates are desirable over a large range of SNR. This simulation is focused on comparing the observation results between different techniques under different CSR scenarios. SNR is fixed to 40 dB. Figures 5(a)–5(c) correspond to a situation with CSR = 10 dB, Figures 5(d)–5(f) to CSR = 20 dB, and Figures 5(g)–5(i) to CSR = 30 dB.
Figures 5(a), 5(d), and 5(g) illustrate the estimates of radial velocity sampled at single PRT 1050 using the pulse pair processing (PPP) method. The corresponding maximum unambiguous velocity is 7.14 m/s. The velocity aliasing can be seen obviously in the region of wind shear field. Because of the ground clutter, the estimation of velocity is highly influenced.
Figures 5(b), 5(e), 5(h), 5(c), 5(f), and 5(i) show the obtained velocity estimates for the 2-PRT (1050 and 1200 ) and 4-PRT (1050, 1200, 1500, and 1950 ) schemes with PTDM, respectively. The maximum unambiguous velocity of both transmission schemes is 50 m/s. The problem of velocity aliasing is well solved. Compared with Figures 5(b) and 5(e), the clutter is filtered out better in Figures 5(c) and 5(f). For the case CSR = 30 dB, the clutter filtering results are almost the same.
In this paper, a technique applying the parametric time domain method to the -PRT transmission schemes is proposed. Based on the assumption that the Doppler spectra of clutter and precipitation follow Gaussian shape, simulation studies are carried out using a case of -PRT schemes.
It is shown that the proposed method can effectively estimate Doppler velocity. Two cases are simulated to analyze the standard deviation and bias of measured velocities, which can help to compare the proposed method with 2-PRT transmission schemes with PTDM. The first one describes a situation where SNR is fixed to 40 dB and CSR varies. The biases are nearly the same for both of the methods. As long as CSR is lower than 20 dB, the standard deviation of measured velocity using the proposed method are almost three times less than the 2-PRT transmission schemes with PTDM. The second type of simulation is performed for CSR = −20 dB and 0 dB ≤ SNR ≤ 40 dB. The estimation error of measured velocities using the 4-PRT schemes is much better than that given by the 2-PRT schemes.
The PPI scan of measured velocity is built on simulated radar data with SNR = 40 dB and different CSRs. The results confirm that the proposed method not only solves Doppler velocity aliasing but also performs better in clutter suppression than 2-PRT schemes with PTDM especially when CSR < 30 dB.
The future works will focus on finding the optimal -PRT schemes which minimize the estimation error and its applications to real air observation. The method to improve the computing efficiency is also the next emphasis of our research.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
The authors wish to thank Dr. Wei Cheng for his continuous encouragement and support. The authors have benefited from in-house reports and information provided to them by the members of the Graduate Innovation Experimental Center.
R. J. Doviak and D. S. Zrnic, Doppler Radar and Weather Observations, Academic Press, 1993.
R. J. Doviak and D. Sirmans, “Doppler radar with polarization diversity,” Journal of the Atmospheric Sciences, vol. 30, no. 4, pp. 737–738, 1973.View at: Google Scholar
M. Sachidananda, D. S. Zrnic, and R. J. Doviak, “Signal design and processing techniques for WSR-88D ambiguity resolution. Part-1,” Tech. Rep., National Severe Storms Laboratory, 1997.View at: Google Scholar
M. Sachidananda, D. S. Zrnic, and R. J. Doviak, “Signal design and processing techniques for WSR-88D ambiguity resolution. Part-2,” Tech. Rep., National Severe Storms Laboratory, 1998.View at: Google Scholar
J. Pirttilä, M. S. Lehtinen, A. Huuskonen, and M. Markkanen, “A proposed solution to the range-Doppler dilemma of weather radar measurements by using the SMPRF codes, practical results, and a comparison with operational measurements,” Journal of Applied Meteorology, vol. 44, no. 9, pp. 1375–1390, 2005.View at: Publisher Site | Google Scholar
A. D. Siggia and R. E. Passarelli, “Gaussian model adaptive processing (GMAP) for improved ground clutter cancellation and moment calculation,” in Proceedings of the 3r European Conference on Radar Meteorology (ERAD '04), pp. 67–73, Visby, Sweden, 2004.View at: Google Scholar
C. M. Nguyen, D. N. Moisseev, and V. Chandrasekar, “A time domain clutter filter for staggered PRT and dual-PRF measurements,” in Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS '07), pp. 3325–3328, Barcelona, Spain, June 2007.View at: Publisher Site | Google Scholar
V. N. Bringi and V. Chandrasekar, Polarimetric Doppler Weather Radar: Principles and Applications, Cambridge University Press, Cambridge, UK, 2001.View at: Publisher Site