ABSTRACT
In this study, the effect of parameter β in Kaiser Window on the SNR values of MST radar echoes is investigated. Use is made of six sets of multibeam observations of the lower atmosphere made by the Indian Mesosphere-Stratosphere-Troposphere (MST) radar. The effects of data weighting with the variation in the parameter β of the Kaiser window applied to the inphase and quadrature components of the radar echo samples prior to Fourier transformation are presented. It is observed that the increase of β, increases the SNR in the higher bins and a good improvement is reported. At the same time in the lower bins a slight decrease in SNR is observed. From these observations, it is concluded that the Kaiser window can be used with β greater than or equal to 6, to taper the data for spectral analysis in place of a rectangular window. The results also show that, the improvement of SNR of noisy data due to the effect of side lobe reduction and demands for the design of optimal windows.
PDF Abstract XML References Citation
How to cite this article
DOI: 10.3923/ajsr.2008.203.212
URL: https://scialert.net/abstract/?doi=ajsr.2008.203.212
INTRODUCTION
Harmonic analysis with the Discrete Fourier Transform (DFT) plays a central role in radar signal processing. The significance of using data weighting windows with the DFT (Harris, 1978; Kay, 1988; Marple, 1987) plays an important role in resolving the frequency components of the signal buried under the noise. Since the use of an in-appropriate window can lead to corruption of the principal spectral parameters it is instructive to consider the criteria by which the choice of data weighting window to be used is made (Hooper, 1999). This paper presents the effect of β in Kaiser window on the SNR of radar returns and proposed an optimum value of β with which data may be weighted using Kaiser window.
The Data Weighting Windows
Windowing
It is well known (Harris, 1978; Kay, 1988; Marple, 1987) that the application of FFT to a finite length data gives rise to leakage and picket fence effects. Weighting the data with suitable windows can reduce these effects. However the use of the data windows other than the rectangular window affects the bias, variance and frequency resolution of the spectral estimates (Kay, 1988; Marple, 1987). In general, the variance of the estimate increases with the use of a window. An estimate is to be consistent if the bias and the variance both tend to zero as the number of observations is increased. Thus, the problem associated with the spectral estimation of a finite length data by the FFT techniques is the problem of establishing efficient data windows or data smoothing schemes.
Data windows are used to weight time series of the in-phase and quadrature components of the radar return samples prior to applying the DFT. The observed Doppler spectra therefore represent convolutions of the Fourier transforms of the original signals with those of the data weighting windows, projected onto the discrete (angular) frequencies (Harris, 1978).
Spectral Leakage
For signal frequencies, observed through the rectangular window, which do not correspond exactly to one of the sampling frequencies, the pattern is shifted such that non-zero values are projected onto all sampling frequencies. This phenomenon of spreading signal power from the nominal frequency across the entire width of the observed spectrum is known as spectral leakage (Harris, 1978).
The radar returns considered to be composed of a quasi-monotonic (atmospheric) signal superimposed on a background of white noise. As might be expected, since the signal does not correspond exactly to one of the sampling frequencies, the forms of the signal portions of the spectra follow those of the envelopes of the side lobe maxima. Spectral leakage from the signal therefore exceeds noise level, evaluated by the method (Hildebrand and Sekhon, 1974) and a corresponding underestimate of signal-to-noise ratio.
The Kaiser Window
It is a relatively simple approximation of the prolate spheroidal functions. For discrete time the Kaiser window is expressed as (Alan and Schafer, 1988; Kay, 1988):
![]() | (1) |
Where:
N | = | The window length, controls the main lobe width |
β | = | Controls the amplitude for the side lobes |
IO(x) | = | The modified zeroth-order Bessel function |
In contrast to the other windows, the Kaiser window has two parameters: The length of the sequence N and a shape parameter β. As the length of the window is fixed to 512 data points here, the shape parameter β can be varied. As the parameter β increases the side lobe level of the frequency response decreases (Alan and Schafer, 1988). In this paper we have investigated the effect of β on the SNR of MST radar data. The results are presented in Fig. 1a-f and 2a-f.
Kaiser-Window Applied to Atmospheric Radar Signals
Wind profile detection of a MST Radar signal meant the measurement of Dopplers of the signal due to scattering of the atmospheric elements. Atmospheric Radar signal meant the signal received by the Radar due to the back scattering property of the atmospheric layers, stratified or turbulent. The back-scattered signal from the atmospheric layers is very small in terms of power with which it was emitted. The received back-scattered signals otherwise called as Radar returns are associated with Gaussian noise. The noise dominates the signal as the distance between the Radar and the target increases and this leads to a decrease in Signal to Noise ratio. This makes the detection of the signal difficult. Doppler profile information is obtained from the power spectrum using Fast Fourier Transform. Frequency characteristics of the back-scattered signals of the Radar s are analyzed with power spectrum, which specifies the spectral characteristics of a signal in frequency domain.
![]() | |
Fig. 1: | Average SNR for the MST Radar data collected from NARL, Gadanki, India on 8th July 2002 |
The specifications of the data are shown in Table 1. The SNR analysis is performed on MST Radar data corresponds to the lower stratosphere obtained from the NARL, Gadanki, India, on 8th July 2002. The Radar was operated in Zenith X, Zenith Y, North, South, West and East with an angle of 10° from the vertical, direction. The data obtained from the six directions are used to carry on the analysis. The complete implementation of the scheme using C ++ and Matlab, to study the effect of β on the SNR of the radar returns can be put as follows.
• | Compute the Kaiser window with β = 0 using Eq. 1. |
• | Taper the radar data with the Kaiser window parameters specified in step (a). |
• | Perform the Fourier analysis of the above tapered data (Appendix A) |
• | Compute the SNR using the procedure (Anandan, 1987; Hooper, 1999) (Appendix B) |
• | Compute the average of SNR minimas (Mean of - dB SNR Values) |
• | Compute the average of SNR maximas (Mean of + dB SNR Values) |
• | Up date the value of β and repeat the steps b-f. |
Table 1: | Specifications of the MST Radar, India data on which the analysis is performed |
![]() | |
E10y = East West polarization with off-zenith angle of 100; W10y = East West polarization with off-zenith angle of 100; N10x = North South polarization with off-zenith angle of 100; S10x = North South polarization with off-zenith angle of 100 |
RESULTS AND DISCUSSION
The SNR computation (Hildebrand and Sekhon, 1974; Anandan, 2000, 1987) (Appendix B) for the six sets of Radar data is carried on and presented in Fig. 1a-f and 2a-f. From the Fig. 1a-f, it is inferred that when β = 0, which corresponds to the rectangular window, the average SNR in the negative side is minimum and increases with the increase in β steadily. The SNR increases slowly when β is varied from 4 to 6, There is no appreciable change in average SNR value for further increase in β.
On the other hand in all the six sets of data, the average positive SNR has not suffered any major change when β is increased. More over, a slight decrease in average positive SNR is observed when β is increased beyond 6.
From the Fig. 2a-f it is also seen that for the lowermost 50 bins the average SNR in the negative and positive sides are not improved appreciably. Moreover there is slight and marginal decrease in both SNRs is observed. For the middle 50 bins and the uppermost 50 bins the increase in average negative SNR values is almost 5-6 dB when β is around 6. Further improvement is also seen when β is increased beyond 6. This result is important since the back-scattered signal from the middle and uppermost bins is very week and improvement in SNR is highly desirable in spectral estimation.
Noting the above observations, it is concluded that the Kaiser window can be used with β greater than or equal to 6 to taper the data for spectral analysis in place of a rectangular window. The results also suggest that, the effect of side lobe reduction in the improvement of SNR of noisy data and demands for the design of optimal windows.
![]() | |
Fig. 2: | Average SNR for Lower, middle and upper bins data collected from NARL, Gadanki, India on 8 th July 2002 |
ACKNOWLEDGMENTS
The first author is thankful to Dr. S. Narayana Rao, Director, NARL, Gadanki, India for permitting to use the facilities at the Radar center. He is thankful to Dr. A.S.R. Reddy, Professor, LBR College of Engineering, Mylavaram, India for his support and help. He is also thankful to Dr. I. Gopal Reddy, Director and Dr. O. Mahesh, Principal, Priyadarshini College of Engineering, Nellore, India for their continued support and encouragement to carry out the present study.
APPENDIX A
The Algorithm presented here to compute the power spectrum estimation is based on the procedure described in the references (Anandan, 2000, 1987).
Power Spectrum Estimation
Step 1
Read the raw data (complex time samples) X(n)
Step 2
Normalize the Pre-Processed data
Step 3
Taper the Normalized data with Kaiser window described in Eq. 1
(d(n) = x(n)*w (n) | (A1) |
Step 4
Compute the Normalized FFT of d (n)
![]() | (A2) |
Step 5
Compute the Power Spectrum Pi
Pi = | D(i)|2 = 0,
N-1 | (A3) |
Step 6
Perform Incoherent Integration (i.e averaging of the power spectrum number of times).
![]() | (A4) |
Where, M is the number of spectra integrated. In this work M = 1 is taken.
Step 7
Power spectrum cleaning is performed to remove Clutter/DC value or zero mean velocities using the basic operation as:
![]() | (A5) |
N/2 corresponds to zero frequency
APPENDIX B
Noise Level Estimation
The Algorithm presented here for Noise level estimation is based on the procedure described in the references (Anandan, 1987, 2000).
Step 1
Reorder the spectrum {Pi, i = 0,... N-1} in ascending order and form the sequence {Ai, i = 0,... N-1} and Ai < Aj for i < j
Step 2
Compute
![]() | (B1) |
![]() | (B2) |
and if Qn > 0,
![]() |
Where, M = 1, the number of spectra that were averaged
Step 3
Compute the Noise Level.
![]() |
Step 4
Reorder the spectrum to its correct index of frequency in the following manner.
Spectral index |
0 | Ambiguous frequency |
1 | Maximum negative frequency |
. | |
. | |
N/2 | Zero frequency |
. | |
. | |
N-1 | Maximum positive frequency |
Step 5
Subtract noise level L from spectrum to form
Step 6
• | Find l the index of the peak value in the spectrum, |
![]() |
• | Find m the minimum index of positive point corresponding to the valley point of the detected peak Doppler spectrum. |
![]() |
• | Find n the maximum index of positive point corresponding to the valley point of the detected peak Doppler spectrum. |
![]() |
Step 7
Compute the zeroth moment or the total power in the Doppler spectrum.
![]() | (B5) |
Step 8
Compute the Signal to Noise Ratio as
![]() | (B6) |
Where: |
N | = | The number of FFT points. |
Step 9
Repeat the Procedures described in Appendix-A and Appendix-B for all the bins of MST radar data.
Step 10
Compute the Mean of -dB SNR Values.
Step 11
Compute the Mean of +dB SNR Values.
Step 12
Update β value and repeat the steps in Appendix A and B.
Step 13
Plot the Mean of -dB SNR Values and Mean of +dB SNR Values with respect to β.
REFERENCES
- Harris, F.J., 1978. On the use of windows for harmonic analysis with the discrete fourier transform. Proc. IEEE, 66: 51-83.
CrossRef - Hildebrand, P.H. and R.S. Sekhon, 1974. Objective determination of the noise level in Doppler spectra. J. Applied Meteorol., 13: 808-811.
CrossRefDirect Link - Hooper, D.A., 1999. Signal and noise level estimation for narrow spectral width returns observed by the Indian MST radar. Radio Sci., 34: 859-870.
CrossRefDirect Link