Next Article in Journal
A Self-Supervised Learning Approach for Extracting China Physical Urban Boundaries Based on Multi-Source Data
Next Article in Special Issue
Enhancing Spatial Variability Representation of Radar Nowcasting with Generative Adversarial Networks
Previous Article in Journal
Sensitivity of Grassland Coverage to Climate across Environmental Gradients on the Qinghai-Tibet Plateau
Previous Article in Special Issue
Technical Evaluation of Precipitation Forecast by Blending Weather Radar Based on New Spatial Test Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Optimization Strategy of Sidelobe Suppression for Pulse Compression Weather Radar

1
School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China
2
The Key Laboratory of Electronic and Information Technology in Satellite Navigation, Ministry of Education, Beijing Institute of Technology, Beijing 100081, China
3
The Aerospace Newsky Technology Co., Ltd., Wuxi 214000, China
4
Beijing Institute of Technology Chongqing Innovation Center, Chongqing 401120, China
5
Chongqing Key Laboratory of Novel Civilian Radar, Chongqing 401120, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(12), 3188; https://doi.org/10.3390/rs15123188
Submission received: 25 May 2023 / Revised: 15 June 2023 / Accepted: 18 June 2023 / Published: 19 June 2023
(This article belongs to the Special Issue Processing and Application of Weather Radar Data)

Abstract

:
The solid-state transmitters are widely adopted for weather radars, where pulse compression is operated to provide the required sensitivity and range resolution. Therefore, effective sidelobe suppression strategies must be employed, especially for weather observation. Currently, many methods can suppress the sidelobe to a very low level in the case of point targets or uniformly distributed targets. However, in strong convection weather process, the weather echo amplitude lies in a wide dynamic range and the main lobe of weak target is prone to being contaminated by the sidelobe of strong target, causing the degradation of weather fundamental data estimation, even generating artifacts and affecting the quantitative precipitation evaluation. In this paper, we propose a novel strategy which is the further processing of a general pulse compression radar to mitigate the effects of sidelobes. The proposed method is called the predominant component extraction (PCE), in which the re-weighting processing is operated after pulse compression, and then the echo of each bin is optimized and its energy will approach the real targets in each bin. It can improve the estimation of weak signals or even eliminate the artifact at the edge of the scene. Numerical simulation experiments and real-data verifications are implemented to validate the feasibility and superiority. It is noted that the proposed method has no requirement on the transmitted waveform and can be realized only by adding a step after pulse compression in the actual system.

Graphical Abstract

1. Introduction

Pulse compression technology can meet the requirements of high-range resolution and long-range detection at the same time, which has been widely used in weather radar systems with the application of solid-state transmitter.
The commonly used pulse compression waveforms include linear frequency modulation (LFM), nonlinear frequency modulation (NLFM), phase coded signal [1,2,3], etc. The LFM waveform which is the most widely used has high-range sidelobe (−13.26 dB) after matched filtering. A window function is often utilized to suppress range sidelobes, resulting in the expansion of the main lobe and the loss of the signal to noise ratio (SNR) [4,5]. However, the reflectivity of precipitations ranges from about −10 to 75 dBZ [6,7] and the radar sensitivity is one of the most critical factors for distributed precipitations in weather observations. Therefore, in order to reduce the SNR loss of LFM signal after windowing, the NLFM waveform was developed decades ago. The concept of NLFM waveform was put forward by Fowle and Brandon in 1959 [8,9]. The energy spectrum (i.e., the square of the spectrum) of NLFM can be designed as a window function, namely, the effect of windowing can be transferred to signal modulation, and then the matching filter can be performed directly to achieve low sidelobe performance, in order to avoid the SMR loss caused by windowing processing [7]. The peak sidelobe ratio (PSLR) can be reduced to below −40 dB after matched filtering, which can achieve the same results of window function weighting processing [10]. At the same time, NLFM has better detection rate characteristics and more accurate range detection performance than LFM [11].
However, some researches have shown that the adjacent reflectivity in extreme weather targets varies dramatically and the gradient often ranges from 30 to 40 dBZ within 1 km, sometimes even reaching 55 dBZ/km [12,13,14]. In order to realize the detection of weather targets with high gradient-reflectivity phenomenon, a very effective sidelobe suppression strategy must be adopted to avoid the artifact caused by range sidelobes [15], while the sidelobe performances of NLFM waveform are not enough to meet the accuracy of the quantitative detection of distributed scatterers, and the sidelobe needs to be further suppressed.
To reach the ultra-low range sidelobes, a mismatch filter can be used at the cost of losing SNR [5,16,17]. Argenti et al. [12] designed transmit waveforms and receive filters using the quadratic nonlinear optimization method by minimizing both PSLR and integrated sidelobe ratio (ISLR) of the waveform at the receiver output, and the PSLR and ISLR can reach −80 dB and −70 dB, respectively with the loss of resolution degradation. Beauchamp et al. [18] discussed the optimal design of pulse compression waveform/filter pairs for use with near-nadir spaceborne radar in low Earth orbit for the observation of clouds and precipitation. It was demonstrated that the LFM waveforms provide superior performance over NLFM waveforms for the application subject to unmitigated Doppler shifts and the PSLR and ISLR could reach −56 dB and −34 dB, respectively using the minimum integrated sidelobe (ISL) mismatch filter. Kurdzo et al. [15] designed NLFM waveform using a genetic algorithm that took into account individual system characteristics and performance measures in order to design a low SNR loss (high power efficiency) waveform for use with weather radar utilizing pulse compression. In addition, the waveform was implemented in the X-band transportable solid-state dual-polarized weather radar system (PX-1000), and the PSLR and ISLR in the actual system can reach −52 dB and −37 dB, respectively. For weather radar, the received signal of one radar resolution volume (RRV) is the sum of scattered signals from the ensemble of particles in this RRV. Therefore, the velocities of the particles in RVV continuously distribute in an interval and the Doppler spectrum or power spectral density (PSD) generally has Gaussian shape [19]. The spectral moments, i.e., reflectivity, mean radial velocity, and spectrum width (ZVW) are called fundamental weather parameters and are defined as the top three order moments of PSD. Bharadwaj et al. [20] used frequency diversity waveform and minimum ISL mismatch filter for pulse compression. The simulation results show that the estimation errors of fundamental weather parameters, i.e., reflectivity, mean radial velocity, and spectrum width (ZVW) for weak targets will increase when the reflectivity changes greatly, which had been verified using CASA (Center for Collaborative Adaptive Sensing of the Atmosphere) X-band dual-polarized radar.
The sidelobe of the point target from the above-mentioned pulse compression filters can be compressed to a very low level, while it can only be achieved for point or uniformly distributed targets, such as layered precipitation. In the process of strong convection weather phenomenon, the reflectivity will appear as a large gradient, and the sidelobe energy in the strong target range bins will be superimposed into the weak target range bins. When the reflectivity gradient reaches a certain range (such as 40 dBZ/km), the sidelobe energy of the strong target is equal to the main lobe energy of the weak target which will be contaminated, causing the unacceptable estimation errors of ZVW, even generating false targets and affecting the estimation of precipitation [15,21].
Let us assume that the probing of the meteorological object is conducted with frequency modulated waveform (such as LFM or NLFM) at a fixed antenna orientation. After pulse compression of the return signal by the optimal filter (OF), such as matched filter or mismatched filter, some realization of a random echo-signal is received. This signal is a mixture of the amplitudes of the main and sidelobes of the response of the OF, resulting in the increase in estimation error of ZVW. To isolate the reaction of the OF only by the main lobe, in this paper, an artificial realization of the echo-signal is initially modeled. This realization consists only of rectangular pulses with amplitudes obtained by smoothing the amplitudes of the real realization of the echo-signal, which is called basic waveform transfer function (B-WTF) to describe the main lobe signals. Then, considering the parameters of the waveform and OF, the frequency-modulated pulse with amplitudes obtained by smoothing the amplitudes of the real realization of the echo-signal is modeled to describe the sidelobe signals, which is called pulse-compression waveform transfer function (PC-WTF) and is supplemented by the main lobe signals. Thereafter, the cost function of each range bin can be constructed by B-WTF and PC-WTF for the calculations of optimal extraction matrix (OEM) for re-weighting of original real realization. The optimized echoes are used to calculate the relevant data quality index according to some criterions. If the criterions are not satisfied, the iteration procedure is repeated with slightly modified initial conditions until the matching criterion is met. Through the above-mentioned procedure, the predominant signal in each range bin can be extracted from the mixed signal and the sidelobes are cleared.
This method is called the predominant component extraction (PCE), which is operated after pulse compression for a normal weather radar with LFM or NLFM waveform to improve sidelobe suppression. After PCE processing, the energy of optimized echoes will be close to the main lobe energy of the actual targets, achieving the sidelobe suppression of the weather targets with large gradient reflectivity.
When the reflectivity gradient is large, the proposed method can obtain the weak target signal and reduce the estimation errors of ZVW. This optimization effect is particularly clear for the large gradient reflectivity scene in the typhoon or other strong convective weather phenomenon. In addition, this method can significantly improve the data quality and eliminate the artifact at the edge of the scene, which contributes to obtaining the accurate estimation of precipitation and other parameters. It is worth mentioning that the proposed method is operated after pulse compression, and it has no requirement on the radar transmitting waveform, which can be realized by adding a step after pulse compression in the actual radar system.
This article is organized as follows. The optimization model is presented in Section 2. Section 3 describes the calculation of the OEM and the iterative optimization process of the PCE. Some results of simulations and verifications based on the real measured data are shown in Section 4. The discussion is drawn in Section 5. Finally, Section 6 provides the conclusion and research perspectives.

2. Problem Statement and Modeling

Considering one range-direction signal, the normalized echo matrix of unit amplitude transmitted modulation waveform, such as LFM or NLFM after pulse compression is defined as the pulse-compression waveform transfer function (PC-WTF) s M F , and the normalized echo matrix with unit amplitude rectangular-pulse transmitted waveform is defined as the basic waveform transfer function (B-WTF) s S P , which are expressed as follows:
s M F = [ s M F , 1 s M F , 2 s M F , k s M F , K ] = [ s M F , 1 ( 1 ) , s M F , 1 ( 2 ) , , s M F , 1 ( p ) , , s M F , 1 ( P ) s M F , 2 ( 1 ) , s M F , 2 ( 2 ) , , s M F , 2 ( p ) , , s M F , 2 ( P ) s M F , k ( 1 ) , s M F , k ( 2 ) , , s M F , k ( p ) , , s M F , k ( P ) s M F , K ( 1 ) , s M F , K ( 2 ) , , s M F , K ( p ) , , s M F , K ( P ) ]
s S P = [ s S P , 1 s S P , 2 s S P , k s S P , K ] = [ s S P , 1 ( 1 ) , s S P , 1 ( 2 ) , , s S P , 1 ( p ) , , s S P , 1 ( P ) s S P , 2 ( 1 ) , s S P , 2 ( 2 ) , , s S P , 2 ( p ) , , s S P , 2 ( P ) s S P , k ( 1 ) , s S P , k ( 2 ) , , s S P , k ( p ) , , s S P , k ( P ) s S P , K ( 1 ) , s S P , K ( 2 ) , , s S P , K ( p ) , , s S P , K ( P ) ]
where each row of (1) and (2) represents the range direction sample sequence and P is the number of range bins of echoes. Each column represents the normalized amplitude of each target applied to the corresponding range bin, and K is the number of targets. In contrast to point targets, the entire radar beam is usually filled with weather targets; therefore, the number of range bins can be treated as equivalent to the number of targets, i.e., P = K . In addition, s M F , k is the PC-WTF of the kth target and s S P , k is the B-WTF of the kth target as follows:
s M F , k = [ s M F , k ( 1 ) , s M F , k ( 2 ) , , s M F , k ( p ) , , s M F , k ( P ) ]
s S P , k = [ s S P , k ( 1 ) , s S P , k ( 2 ) , , s S P , k ( p ) , , s S P , k ( P ) ]
Due to the absence of sidelobes of rectangular-pulse waveform and the presence of sidelobes of pulse compression waveform, s M F , k ( p ) and s S P , k ( p ) which are the members of (3) and (4), respectively, can be expressed as follows:
s M F , k ( p ) = { a 0 p = k b ( p ) p k
s S P , k ( p ) = { a 0 p = k 0 p k
where a 0 is the main lobe value after pulse compression of kth target’s echo, k is the range bin position where the main lobe is located, and b ( p ) is the sidelobe value at other range bin locations of kth target’s echo.
Therefore, the values of a 0 and b ( p ) are only related to the waveform parameters, the radar system parameters, and the method of pulse compression, which are independent of the target characteristics. In particular, as long as the waveform parameters, the radar system parameters, and the method of pulse compression are determined, s M F and s S P will be determined. Considering the PC-WTF as an example, its schematic diagram is shown in Figure 1:
Let us assume that A = [ A 1 A k A K ] is the backscattering coefficient matrix, and A k is the complex backscattering coefficient of the kth target, where the amplitude is related to the scattering intensity of the target and the phase includes the information of doppler velocity of the target. The received echo is the superposition from all targets as shown in Figure 2, which can be expressed as (7):
P 0 = A s M F = [ k = 1 K A k s M F , k ( 1 ) , k = 1 K A k s M F , k ( 2 ) , , k = 1 K A k s M F , k ( P ) ]
It can be seen from (7) that the energy of each range bin is the superposition of the main lobe of this bin and the sidelobes from other bins. When the energy of the extra added sidelobes is comparable to that of the original main lobe, it will affect the estimations of ZVW. Therefore, it is considered to re-optimize the signal of each range bin and its energy will approach that of the actual targets in each bin; therefore, we call it the predominant component extraction (PCE) method.
Considering the operation of the weighting processing of the signals in each range bin, we can obtain the weighted signal as follows:
P e s t = A ( W o p t s M F )
where W o p t is the OEM with K × P dimension as follows:
W o p t = [ w o p t , 1 ( 1 ) w o p t , 2 ( 1 ) w o p t , p ( 1 ) w o p t , P ( 1 ) w o p t , 1 ( 2 ) w o p t , 2 ( 2 ) w o p t , p ( 2 ) w o p t , P ( 2 ) w o p t , 1 ( k ) w o p t , 2 ( k ) w o p t , p ( k ) w o p t , P ( k ) w o p t , 1 ( K ) w o p t , 2 ( K ) w o p t , p ( K ) w o p t , P ( K ) ]
where w o p t , p = [ w o p t , p ( 1 ) w o p t , p ( K ) ] T is the optimization coefficient of pth range bin and is Hadamard product. The signal is weighed using w o p t , p in each range bin in order that the energy will approach that of the actual targets. Next, an optimization problem is established by B-WTF to obtain W o p t .
Since the echoes of rectangular pulse signal only have the main lobe and there is no sidelobe to be superimposed in the other range bins, the calculated signal energies are optimal in theory. Therefore, we can use rectangular waveforms to model the expected output echo:
P exp = A s S P = [ k = 1 K A k s S P , k ( 1 ) , k = 1 K A k s S P , k ( 2 ) , , k = 1 K A k s S P , k ( P ) ] = [ A 1 s S P , 1 ( 1 ) , A 2 s S P , 2 ( 2 ) , , A P s S P , P ( P ) ]
In order to obtain P exp , W o p t needs to meet the following condition:
min W o p t p ( W o p t ) s . t . p ( W o p t ) = P e s t P exp 2 2
To date, the optimization model has been established.
It should be noted that, in this paper, we assume that the echoes are entirely from meteorological targets and the patterns of targets are determined by the backscattering coefficient matrix and waveform transfer function. However, the PCE method can be operated regardless of how the target is distributed. We do not need to know the value of backscattering coefficient matrix in advance when modeling PC-WTF and B-WTF. As long as the waveform and radar system parameters are determined, PC-WTF and B-WTF can be modeled, which can be realized in the actual processing process. After modeling PC-WTF and B-WTF, what we need to solve is how to optimize (7) to make it closer to (10), which needs to be realized by the PCE method proposed in this paper. The block diagram of the PCE algorithm is presented in Figure 3.

3. Optimization Process

3.1. The Solution of OEM

After obtaining the above-mentioned optimization model, W o p t can be calculated by the derivative of p ( W o p t ) and setting the derivative to zero. Next, we will introduce the calculation process in detail.
Expanding (11), we can achieve:
min W o p t p ( W o p t ) s . t . p ( W o p t ) = A ( W o p t s M F ) A s S P 2 2
For distributed meteorological targets, P = K can be assumed and A is a matrix that does not change in one coherent processing interval (CPI). At this time, a unit matrix I can be introduced to facilitate the derivative calculation of Hadamard product. Then, the above-mentioned optimization problem can be written as:
min W o p t p ( W o p t ) s . t . p ( W o p t ) = ( W o p t I s M F ) s S P 2 2
Additionally, W o p t can be calculated from d [ p ( W o p t ) ] d W o p t = 0 , as shown in (14), where the detailed calculation process is in Appendix A.
W o p t = ( I s M F s M F I T ) 1 ( I s M F s S P )

3.2. Iteration Process

After obtaining W o p t , the echoes can be optimized through formula (8). However, the initial matrix A needs to be known. Thereafter, the data quality can be improved through iteration. If the actual received echo is S r , we can express it as:
S r = A s M F
A can be calculated as:
A = [ ( s M F T ) 1 S r T ] T
Then, the estimated echo by the PCE method can be expressed as:
S e s t = A ( W o p t s M F )
Thereafter, the OEM of (17) is recalculated and optimized iteratively until the stop criterion is met. Next, we set the appropriate stop criterion to make the data quality meet the requirements.
Since the proposed method in this paper is not aimed at the point target sidelobe, the performance of PSLR and ISLR cannot be compared. Therefore, other indicators that can reflect the data quality need to be set as the stop criterion.
  • Data error iteration condition:
The root mean square error (RMSE) reflects the degree of data deviation from the ground truth. In the quantitative analysis of simulation, the ground truth is given; therefore, the stop criterions can be set according to the RMSE. The RMSE of one range direction can be calculated as:
R M S E E = 1 P p = 1 P | E e s t ( p ) E g t ( p ) | 2
where E e s t represents the ZVW (i.e., reflectivity Z e s t , radial velocity V e s t , and spectral width W e s t ) calculated from the optimized echoes and E g t is the ground truth of ZVW.
Since the energy optimization effect of each range bin is our concern in this paper, we can set the RMSE of the reflectivity less than a certain threshold ε 1 as the stop criterion, as follows:
R M S E Z e s t < ε 1
  • Data fluctuation iteration condition:
In real-data verification, the ground truth cannot be obtained, while the moving average can be used as the ground truth to calculate SD. Therefore, the SD can reflect the fluctuation of radar data and we call it data fluctuation, which can be used to set the stop criterion. The standard deviation (SD) of ZVW in the pth range bin can be calculated by the data on n range bins before and n range bins after, and the first n range bins and the last n range bins do not participate in the calculation of SD [22]:
S D p = 1 2 n + 1 j = 0 2 n + 1 ( E p + j n E ¯ p ) 2 , 0 < n < P 2 , n < p < P n E ¯ p = 1 2 n + 1 j = 0 2 n + 1 E p + j n , 0 < n < P 2 , n < p < P n
where E represents the ZVW (i.e., reflectivity Z e s t , radial velocity V e s t , and spectral width W e s t ) calculated from the optimized echoes. E ¯ p is the moving average calculated from 2 n + 1 range bins.
Setting the percentage of data with the SD less than 1 in the total data greater than ε 2 as the stop criterion, we can make the data fluctuation meet the requirements. This percentage is defined as the data fluctuation qualification rate AR, as follows:
A R = sum ( S D p < 1 ) P × 100 % > ε 2 , 0 < ε 2 < 100 %
In summary, Algorithm 1 can be summarized as follows.
Algorithm 1: The PCE algorithm
  • STEP 1: Model the PC-WTF and B-WTF according to the radar system and waveform parameters.
  • STEP 2: Calculate W o p t according to (14).
  • STEP 3: Calculate the initial matrix A according to (16).
  • STEP 4: Calculate the first optimization echo S e s t according to (17) and estimate the spectral moments, i.e., Z e s t , V e s t , and W e s t .
  • STEP 5: In the quantitative analysis of simulation stage, calculate the RMSE and AR of Z e s t and set the iteration thresholds ε 1 and ε 2 Judge whether the stop criterion of (19) and (21) are met. In the real-data verification stage, calculate the AR of Z e s t and set the iteration threshold ε 2 Judge whether the stop criterion of (21) is met.
  • STEP 6: Repeat steps 3 to 5 until the stop criterion is satisfied.

4. Results

Before the distributed target simulation, first, the point target scene is simulated and the performances of PSLR and ISLR of four different pulse compression methods are compared, i.e., LFM waveform with windowing, LFM waveform with mismatch filtering, NLFM waveform with windowing, and NLFM waveform with mismatch filtering, wherein the NLFM signal is constructed according to [23] and the mismatch filter is constructed according to the method in [24]. Then, the input ZVW is used to simulate one-dimensional echoes of distributed targets in the range direction. The pulse compression method with the best performance is used for processing. The pulse pair processing method (PPP) is used for spectral moment estimation [21]. This processing set can explain the influence of sidelobe superposition of distributed meteorological targets. Thereafter, the PCE method is used and the RMSE and AR are calculated through (18) and (21) to quantitatively analyze the advantages of the proposed method. Finally, the real in-phase and quadrature-phase (I/Q) data of an actual supercell from a ground-based weather radar and global precipitation measurement (GPM) precipitation data are obtained, which are applied to the proposed method for verification.

4.1. Verifications Based on Numerical Simulation Experiments

4.1.1. Sidelobe Superposition Effects

First, the echo of point target is simulated, and the parameters are shown in Table 1. The processing results using different pulse compression methods are in Figure 4, where the window function is Hamming window function.
The performances of the four pulse compression methods are shown in Table 2, where “W” refers to the windowing process. It can be indicated that the windowing process will cause serious SNR loss and the loss of SNR can be mitigated by using NLFM waveform. In addition, the performances of the combination of NLFM waveform and mismatch filter are best, which can suppress PSLR and ISLR below −60 dB and −35 dB, respectively.
Next, the one-dimensional echoes of distributed targets with different reflectivity gradients are simulated. The simulation parameters are shown in Table 3, where the SNR of the input signal is 60 dB. The target number is equal to 350, which is equivalent to the total range bin number. The range revolution is 37.5 m and the detection distance is 13.125 km. In addition, the weak targets account for 7.5 km and the strong targets account for 3.6 km. The transmit waveform is NLFM in [23] and the mismatch filter in [24] is used for pulse compression. The input Z with different gradients is shown in Figure 5. In addition, the reflectivity, mean velocity, and spectral width can be obtained by the PPP method in Figure 6, Figure 7 and Figure 8, respectively.
Although the performance of the pulse compression method is good enough, when the reflectivity gradient is large (such as 40 dBZ/km), the energy of the sidelobe of the strong target is equal to or even greater than the energy of the weak target, resulting in the estimation errors of the ZVW in the weak target area.

4.1.2. Optimization Process

In order to solve the above-mentioned problems, the PCE algorithm is used for each range bin. The comparison results are shown in Figure 9, Figure 10 and Figure 11. The RMSE and AR of ZVW before and after optimization can be calculated through (18) and (21), respectively. The results are shown in Table 4.
It can be indicated that the RMSE and AR can be improved by the proposed method. Moreover, we calculate the RMSE and AR changing with the gradient of reflectivity in the range of 10~40 dBZ/km. The results are shown in Figure 12 and Figure 13. When the reflectivity gradient is large (≥30 dBZ/km), the Z and V estimations are seriously deviated from the ground truth due to the fact that the RMSEs are large. At this time, the data are invalid and the calculated AR is not referential. The proposed method can greatly reduce the RMSE, in order that the invalid data can become valid. When the reflectivity gradient is not large (<30 dBZ/km), the RMSE is acceptable and the data are valid. At this time, we focus on the results of AR. The ARs of optimized data are more acceptable and the data quality is effectively improved.
In addition, we compared the SNR loss caused by the proposed method. In contrast to point targets, the distributed targets cannot directly compare the energy of the main lobe. Therefore, the SNR loss is defined as the mean difference SNR before and after optimization:
S N R l o s s = 1 P p = 1 P [ S N R b e f o r e ( p ) S N R a f t e r ( p ) ]
The SNR loss changing with reflectivity gradient is shown in Figure 14. It can be indicated that with the increase in reflectivity gradient, SNR loss will become serious. Therefore, the proposed method is to improve the data quality at the cost of loss of SNR.

4.2. Verifications Based on Real Data

4.2.1. Verifications by Ground-Based Weather Radar Data

The radar is located at WRCP station (32.75°N, 119.35°E) with an altitude of 32.3 m and works in the conical scanning mode with the azimuth angle of 0 to 360 degrees. The operating parameters are listed in Table 5.
Supercells, as one of the important mesoscale weather systems, can form severe convective weather, such as heavy precipitation, thunderstorm wind, hail, and tornado. In the radar map, the supercell appears as a tightly organized image of high reflectivity, which may have a hook echo. The in-phase and quadrature(I/Q) data originated from a supercell appearing in 22:52 UTC on 17 July 2020. The echo data contained 902 CPI, 32 pulses per CPI, and 652 range gates per pulse. The transmit waveform is LFM and the matched filtering and Hamming window are used.
The ZVWs of the supercell are shown in Figure 15 and the results after optimization by PCE are shown in Figure 16. Due to the large scale of the supercell, the contrast result is not clear in the complete image. Since the proposed method improved the reflectivity (i.e., echo power) significantly, we enlarged the local reflectivity for comparison as shown in Figure 17. We can see that the artifacts at the scene edge are basically removed after optimization and the target edge is clearer due to the fact that the sidelobe is suppressed.
The ARs of the ZVW are calculated in Table 6. It can be indicated that the percentages of ZVW with the SD less than 1 in total data are increased after optimization, which indicated the effective improvement of the data quality.

4.2.2. Verifications by Quantitative Precipitation Estimation (QPE) Data

In order to verify the application value of the proposed method, we use the Z–I relationship to estimate the precipitation in positions A and B, in which Z is the reflectivity calculated by the unoptimized data and optimized data, respectively, and I is the precipitation intensity in mm/h. The same Z–I relationship is used for both, which is Z = 70 I 1 . 38 fitted from a supercell in [25].
The ground truth of precipitation is the global precipitation measurement (GPM) from Goddard Earth Sciences Data and Information Services Center (GES) DISK) [26] with a time resolution of 0.5 h and a spatial resolution of 0.1°× 0.1°.
In positions A and B, we select 3 × 3 pixels with an adjacent distance of 10 km. For each pixel, the reflectivity of the range bins within the four surrounding spatial resolutions (0.1° × 0.1°) is extracted to calculate the precipitations which are averaged successively, and the precipitation estimates at the pixel can be obtained. The selected region is shown in Figure 18. The precipitations at nine pixel points in the two scenes are calculated using the optimized and unoptimized reflectivity and the results are compared as shown in Figure 19. Then, the differences between the estimation and ground truth for optimized and unoptimized data are calculated to obtain the estimation errors as shown in Figure 20.
It is indicated that the precipitations calculated by the optimized reflectivity are more accurate and the RMSE is effectively decreased for the convective region (Position A) that has large gradient reflectivity, which indicates that the proposed method can improve the quality of precipitation estimation. For the stratified precipitation with low rainfall (Position B), the improvement effect is not clear. Therefore, the proposed method is more suitable for strong convection weather targets.
To date, the application value of the proposed method in this paper is illustrated by the real data of ground-based weather radar and QPE data. This method can eliminate artifacts at the edge of the scene and improve the data quality of ZVW and precipitation.

5. Discussion

In the numerical simulation experiments, the performances of the combination of NLFM waveform and mismatch filter are best, which can suppress PSLR and ISLR below −60 dB and −35 dB, respectively. Then, this pulse compression method is used in the distributed targets with different reflectivity gradients, while the energy of the sidelobe of the strong target is equal to or even greater than the energy of the weak target when the reflectivity gradient is large, resulting in the estimation errors of the ZVW in the weak target area, as shown in Figure 6, Figure 7 and Figure 8. It is indicated that the point target performance (such as PSLR and ISLR) of the pulse compression method is good enough, while it cannot meet the requirement of the high-gradient reflectivity distributed targets. At this time, if we operate the PCE algorithm after pulse compression, the weak targets can be reconstructed and the ZVW can be calculated more accurately, as shown in Figure 12 and Figure 13. In addition, the RMSE and AR results in Figure 12 and Figure 13 can verify the better performance of the PCE algorithm.
Furthermore, in real-data verification, the ground-based weather radar data and QPE data are used to validate the feasibility and superiority of the proposed method. First, the PCE algorithm can eliminate the artifact in the scene and the contour edge of reflectivity is smoother, as shown in Figure 16. In addition, the percentages of ZVW with the SD less than 1 in total data are increased after optimization, as shown in Table 6. These results indicated the effective improvement of the data quality after PCE algorithm. Finally, the precipitations calculated by the optimized reflectivity are more accurate, especially for the convective region as shown in Figure 19 and Figure 20, indicating the advancement of the proposed method in the precipitation estimation.

6. Conclusions

In this paper, a novel sidelobe suppression strategy based on the extraction and iteration of weather radar called PCE is proposed. The cost function is constructed by modeling the transfer function of each range bin for the calculations of OPE. Through PCE processing, the energy of optimized echoes will be close to the main lobe energy of the actual targets in this range bin, achieving the sidelobe suppression of the weather targets with large gradient reflectivity. The proposed method is operated after pulse compression and it has no requirement on the radar transmitting waveform, which can be realized by adding a step after pulse compression in the actual radar system.
It is indicated from numerical simulation experiments that when the reflectivity gradient is large, the proposed method can greatly reduce the estimation errors for weak target, which is at the cost of loss of SNR. When the reflectivity gradient is not large, the data qualities are also effectively improved. The real-data verifications indicate that the proposed method can eliminate artifacts at the edge of the scene and improve the data quality of ZVW and precipitation.
Furthermore, although we only used the measured radar data of supercells to verify the proposed method in this paper, the method can also be applied to other strong convective weather phenomenon such as eyewall of a hurricane and a convective cell. These will be studied in the follow-up work.

Author Contributions

Conceptualization, J.H. and X.D.; methodology, J.H. and X.D.; software, J.H.; validation, J.H. and X.D.; formal analysis, J.H. and X.D.; investigation, J.H. and W.T.; resources, J.H. and C.H.; data curation, J.H. and K.F.; writing—original draft preparation, J.H.; writing—review and editing, J.H. and X.D.; visualization, J.H. and J.L.; supervision, X.D. and C.H.; project administration, W.T.; funding acquisition, X.D. and C.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “the Special Fund for Research on National Major Research Instruments (NSFC Grant Nos. 61827901)” and “the National Science Fund for Distinguished Young Scholars (Grant Nos. 62225104)” and “the National Natural Science Foundation of China (Grant Nos. 61971039, 61960206009, 61971037)” and “the Natural Science Foundation of Chongqing (Grant No. cstc2020jcyj-msxmX0621)” and “Distinguished Young Scholars of Chongqing (Grant No. cstc2020jcyj-jqX0008)” and “National Ten-thousand Talents Program ‘Young top talent’ (Grant No. W03070007)”.

Data Availability Statement

The data are unavailable due to privacy.

Acknowledgments

The authors would like to thank Aerospace Newsky Technology Co., Ltd. for providing real ground-based weather radar data.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The optimal extraction matrix can be calculated as follows.
According to (13) p ( W o p t ) can be expanded as:
p ( W o p t ) = ( W o p t I s M F ) s S P 2 2 = [ ( W o p t I s M F ) s S P ] [ ( W o p t I s M F ) s S P ] T = ( W o p t I s M F ) 2 2 2 ( W o p t I s M F ) s S P s S P 2 2
Considering the derivative and setting it to zero, the optimal W o p t can be obtained.
Next, we consider the derivative of each term in:
S = W o p t I s M F f = W o p t I s M F 2 2 = S : S g = 2 ( W o p t I s M F ) s S P = 2 S s S P
where “:” represents the matrix inner product.
  • Calculation of d S :
According to the property of the derivative of Hadamard product, d S can be calculated as:
d S = ( d W o p t I ) s M F + W o p t I d s M F
Since s M F is independent of W o p t , we can achieve:
d S = ( d W o p t ) I s M F
  • Calculation of d f d W o p t :
According to the properties of matrix inner product and matrix differential [27]:
A : B = t r ( A T B )
d f = t r [ ( d f d X ) T d X ]
Then, d f can be calculated from (A4) and (A5):
d f = d S : S + S : d S = 2 ( W o p t I s M F ) : [ ( d W o p t ) I s M F ] = 2 ( W o p t I s M F s M F I T ) : ( d W o p t ) = t r [ 2 ( W o p t I s M F s M F I T ) T ( d W o p t ) ]
Comparing (A7) with (A6), we can achieve:
d f d W o p t = 2 W o p t I s M F s M F I T
  • Calculation of d g d W o p t :
Similarly, d g d W o p t can be calculated as:
d g d W o p t = 2 ( d S d W o p t ) s S P + 2 S d s S P d W o p t
Since s S P is independent of W o p t , we can achieve:
d g d W o p t = 2 ( d ( W o p t I s M F ) d W o p t ) s S P = 2 ( d ( W o p t I ) d W o p t s M F + W o p t I d s M F d W o p t ) s S P = 2 I s M F s S P
In addition, since s S P is independent of W o p t , the derivative of s S P 2 2 is zero. Therefore, in combination with (A8) and (A10), the derivative of (A1) can be calculated as:
d [ p ( W o p t ) ] d W o p t = d f d W o p t d g d W o p t = 2 W o p t I s M F s M F I 2 I s M F s S P
Additionally, W o p t can be calculated from d [ p ( W o p t ) ] d W o p t = 0 and can be expressed as:
W o p t = ( I s M F s M F I T ) 1 ( I s M F s S P )

References

  1. Skolnik, M.I. Radar Handbook, 3rd ed.; McGraw-Hill: New York, NY, USA, 2008. [Google Scholar]
  2. Richards, M.A. Fundamentals of Radar Signal Processing; McGraw-Hill Processional: New York, NY, USA, 2005. [Google Scholar]
  3. Cook, E.; Bernfeld, M. Radar Signals—An Introduction to Theory and Applications; Artech House: Norwood, MA, USA, 1993. [Google Scholar]
  4. Griffiths, H.D.; Vinagre, L. Design of Low-Sidelobe Pulse Compression Waveforms. Electron. Lett. 1994, 30, 1004–1005. [Google Scholar] [CrossRef]
  5. Cilliers, J.E.; Smit, J.C. Pulse Compression Sidelobe Reduction by Minimization of Lp-Norms. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 1238–1247. [Google Scholar] [CrossRef]
  6. Bringi, V.N.; Chandrasekar, V. Polarimetric Doppler Weather Radar: Principles and Applications; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  7. Pang, C.; Hoogeboom, P.; Le Chevalier, F.; Russchenberg, H.W.; Dong, J.; Wang, T.; Wang, X. A Pulse Compression Waveform for Weather Radars with Solid-State Transmitters. IEEE Geosci. Remote Sens. Lett. 2015, 12, 2026–2030. [Google Scholar] [CrossRef]
  8. Brandon, P.S. The Design of a Nonlinear Pulse Compression System to Give a Low Loss High Resolution Radar Performance. Marconi Rev. 1973, 36, 1–45. [Google Scholar]
  9. Fowle, E. The Design of FM Pulse Compression Signals. IEEE Trans. Inform. Theory 1964, 10, 61–67. [Google Scholar] [CrossRef]
  10. Patton, L.K. On the Satisfaction of Modulus and Ambiguity Function Constraints in Radar Waveform Optimization for Detection. Ph.D. Thesis, Wright State University, Dayton, OH, USA, 2009. [Google Scholar]
  11. Luszczyk, M.; Labudzinski, A. Side Lobe Level Reduction for Complex Radar Signals with Small Base. In Proceedings of the 2012 13th International Radar Symposium, Warsaw, Poland, 23–25 May 2012; pp. 146–149. [Google Scholar]
  12. Argenti, F.; Facheris, L. Radar Pulse Compression Methods Based on Nonlinear and Quadratic Optimization. IEEE Trans. Geosci. Remote Sens. 2021, 59, 3904–3916. [Google Scholar] [CrossRef]
  13. Hu, J.; Dong, X.; Hu, C. Research on Nonlinear Wind Field Retrieval Method Based on Modified VAD Analysis. J. Signal Process. 2021, 37, 284–291. [Google Scholar] [CrossRef]
  14. Wang, X.; Zhai, W.; Greco, M.; Gini, F. Cognitive Sparse Beamformer Design in Dynamic Environment via Regularized Switching Network. IEEE Trans. Aerosp. Electron. Syst. 2022, 59, 1816–1833. [Google Scholar] [CrossRef]
  15. Kurdzo, J.M.; Cheong, B.L.; Palmer, R.D.; Zhang, G.; Meier, J.B. A Pulse Compression Waveform for Improved-Sensitivity Weather Radar Observations. J. Atmos. Ocean. Technol. 2014, 31, 2713–2731. [Google Scholar] [CrossRef] [Green Version]
  16. Ackroyd, M.H.; Ghani, F. Optimum Mismatched Filters for Sidelobe Suppression. IEEE Trans. Aerosp. Electron. Syst. 1973, AES-9, 214–218. [Google Scholar] [CrossRef]
  17. Cohen, M.N.; Baden, J.M. Pulse Compression Coding Study; Georgia Tech Research Institute: Atlanta, GA, USA, 1983. [Google Scholar]
  18. Beauchamp, R.M.; Tanelli, S.; Peral, E.; Chandrasekar, V. Pulse Compression Waveform and Filter Optimization for Spaceborne Cloud and Precipitation Radar. IEEE Trans. Geosci. Remote Sens. 2017, 55, 915–931. [Google Scholar] [CrossRef]
  19. Cavallaro, S. Statistical Properties of Polarimetric Weather Radar Returns for Nonuniformly Filled Beams. IEEE Geosci. Remote Sens. Lett. 2017, 14, 1584–1588. [Google Scholar] [CrossRef]
  20. Bharadwaj, N.; Chandrasekar, V. Wideband Waveform Design Principles for Solid-State Weather Radars. J. Atmos. Ocean. Technol. 2012, 29, 14–31. [Google Scholar] [CrossRef] [Green Version]
  21. Dong, X.; Hu, J.; Hu, C.; Chen, Z.; Li, Y. Rapid Identification and Spectral Moment Estimation of Non-Gaussian Weather Radar Signal. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–12. [Google Scholar] [CrossRef]
  22. Li, S.; Yang, M.; Li, L. Evaluation on Data Quality of X-Band Dual Polarization Weather Radar Based on Standard Deviation Analysis. J. Arid Meteorol. 2019, 37, 467–476. [Google Scholar]
  23. Levanon, N.; Mozeson, E. Radar Signals; Wiley: Hoboken, NJ, USA, 2004. [Google Scholar]
  24. Sun, Y.; Liu, Q.; Cai, J.; Long, T. A Novel Weighted Mismatched Filter for Reducing Range Sidelobes. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 1450–1460. [Google Scholar] [CrossRef]
  25. He, K.; Fan, Q.; Li, K. Z-R Relation with Its Application to Typhoon Precipitation in Zhoushan. J. Appl. Meteorol. 2007, 18, 573–576. [Google Scholar]
  26. Huffman, G.J.; Stocker, E.F.; Bolvin, D.T.; Nelkin, E.J.; Tan, J. GPM IMERG Final Precipitation L3 Half Hourly 0.1 Degree x 0.1 Degree V06 2019; Goddard Earth Sciences Data and Information Services Center (GES DISC): Greenbelt, MD, USA, 2019. [Google Scholar]
  27. Zhang, X. Matrix Analysis and Applications; Tsinghua University Press: Beijing, China, 2013. [Google Scholar]
Figure 1. The schematic diagram of PC-WTF.
Figure 1. The schematic diagram of PC-WTF.
Remotesensing 15 03188 g001
Figure 2. The schematic diagram of received echoes.
Figure 2. The schematic diagram of received echoes.
Remotesensing 15 03188 g002
Figure 3. The block diagram of the PCE algorithm.
Figure 3. The block diagram of the PCE algorithm.
Remotesensing 15 03188 g003
Figure 4. The pulse compression results of point target.
Figure 4. The pulse compression results of point target.
Remotesensing 15 03188 g004
Figure 5. The input Z with different gradients. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Figure 5. The input Z with different gradients. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Remotesensing 15 03188 g005
Figure 6. The reflectivity calculated from unoptimized echoes. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Figure 6. The reflectivity calculated from unoptimized echoes. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Remotesensing 15 03188 g006
Figure 7. The mean velocity calculated from unoptimized echoes. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Figure 7. The mean velocity calculated from unoptimized echoes. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Remotesensing 15 03188 g007
Figure 8. The spectral width calculated from unoptimized echoes. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Figure 8. The spectral width calculated from unoptimized echoes. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Remotesensing 15 03188 g008
Figure 9. The comparison results of reflectivity. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Figure 9. The comparison results of reflectivity. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Remotesensing 15 03188 g009
Figure 10. The comparison results of mean velocity. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Figure 10. The comparison results of mean velocity. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Remotesensing 15 03188 g010
Figure 11. The comparison results of spectral width. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Figure 11. The comparison results of spectral width. (a) Gradient of Z is 15 dBZ/km; (b) gradient of Z is 25 dBZ/km; (c) gradient of Z is 40 dBZ/km.
Remotesensing 15 03188 g011
Figure 12. The RMSE of ZVW changing with the gradient of reflectivity. (a) Reflectivity; (b) mean velocity; (c) spectral width.
Figure 12. The RMSE of ZVW changing with the gradient of reflectivity. (a) Reflectivity; (b) mean velocity; (c) spectral width.
Remotesensing 15 03188 g012
Figure 13. The AR of ZVW changing with the gradient of reflectivity. (a) Reflectivity; (b) mean velocity; (c) spectral width.
Figure 13. The AR of ZVW changing with the gradient of reflectivity. (a) Reflectivity; (b) mean velocity; (c) spectral width.
Remotesensing 15 03188 g013
Figure 14. The SNR loss changing with reflectivity gradient.
Figure 14. The SNR loss changing with reflectivity gradient.
Remotesensing 15 03188 g014
Figure 15. The ZVW of the supercell. ((a) Reflectivity; (b) velocity; (c) spectral width).
Figure 15. The ZVW of the supercell. ((a) Reflectivity; (b) velocity; (c) spectral width).
Remotesensing 15 03188 g015
Figure 16. The results after optimization by PCE. ((a) Reflectivity; (b) velocity; (c) spectral width).
Figure 16. The results after optimization by PCE. ((a) Reflectivity; (b) velocity; (c) spectral width).
Remotesensing 15 03188 g016
Figure 17. Local reflectivity before and after optimization. ((a) Not optimized; (b) optimized).
Figure 17. Local reflectivity before and after optimization. ((a) Not optimized; (b) optimized).
Remotesensing 15 03188 g017
Figure 18. Selected region for QPE.
Figure 18. Selected region for QPE.
Remotesensing 15 03188 g018
Figure 19. The comparison results of precipitations. (a) Position A; (b) position B.
Figure 19. The comparison results of precipitations. (a) Position A; (b) position B.
Remotesensing 15 03188 g019
Figure 20. The estimation errors of precipitations. (a) Position A; (b) position B.
Figure 20. The estimation errors of precipitations. (a) Position A; (b) position B.
Remotesensing 15 03188 g020
Table 1. The simulation parameters of point target.
Table 1. The simulation parameters of point target.
ParametersValues
Frequency (GHz)13.6
Band width (MHz)4
Pulse width (μs)64
Sample frequency (MHz)16
Target distance (km)15
Table 2. The performances of the four pulse compression methods.
Table 2. The performances of the four pulse compression methods.
LFM
(MF + W)
LFM
(WMMF)
NLFM
(MF + W)
NLFM
(WMMF)
PSLR (dB)−53.34−55.31−60.97−67.89
ISLR (dB)−25.66−27.84−28.94−37.23
SNR loss (dB)6.121.843.160.14
Table 3. The simulation parameters of one-dimensional echoes of distributed targets with different reflectivity gradients.
Table 3. The simulation parameters of one-dimensional echoes of distributed targets with different reflectivity gradients.
ParametersValuesParametersValues
Frequency (GHz)13.6Pulse width (μs)64
Band width (MHz)4Sample frequency (MHz)16
Accumulation pulse number 64 PRF (Hz)1860/1395
Target number350Gradient of Z (dBZ/km)15/25/40
Target spectral width (m/s)Randi [1, 5]Target velocity (m/s)Int [−7, 7]
Table 4. The RMSE and AR of ZVW (not optimized/optimized).
Table 4. The RMSE and AR of ZVW (not optimized/optimized).
Gradient   of   Z
(ΔZ, dBZ/km)
152540
Z (dBZ)RMSE (dBZ)2.65/2.323.30/2.3410.42/4.23
AR (%)24.57/50.0132.37/45.6628.90/28.61
V (m/s)RMSE (m/s)2.02/1.142.93/1.204.13/1.61
AR (%)39.02/85.8428.61/82.9536.99/49.42
W (m/s)RMSE (m/s)1.92/0.721.99/0.692.01/1.56
AR (%)36.13/50.8745.66/48.8439.60/59.54
Table 5. The operating parameters of ground-based weather radar.
Table 5. The operating parameters of ground-based weather radar.
ParametersValueParametersValue
Frequency (GHz)5.5Peak Power (W)1710
Antenna gain (dB, T/R)34/38Beam width (°, AZ/EL)1/3
Noise figure (dB)3Pulse number32
Pulse width (μs)100Bandwidth (MHz)1
Dual PRF (Hz)900/1200Elevation angle (°)2
Table 6. The AR of the ZVW.
Table 6. The AR of the ZVW.
ZVW
The AR before optimization (%)62.6864.9076.39
The AR after optimization (%)70.6670.2084.54
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hu, J.; Dong, X.; Tian, W.; Hu, C.; Feng, K.; Lu, J. A Novel Optimization Strategy of Sidelobe Suppression for Pulse Compression Weather Radar. Remote Sens. 2023, 15, 3188. https://doi.org/10.3390/rs15123188

AMA Style

Hu J, Dong X, Tian W, Hu C, Feng K, Lu J. A Novel Optimization Strategy of Sidelobe Suppression for Pulse Compression Weather Radar. Remote Sensing. 2023; 15(12):3188. https://doi.org/10.3390/rs15123188

Chicago/Turabian Style

Hu, Jiaqi, Xichao Dong, Weiming Tian, Cheng Hu, Kai Feng, and Jun Lu. 2023. "A Novel Optimization Strategy of Sidelobe Suppression for Pulse Compression Weather Radar" Remote Sensing 15, no. 12: 3188. https://doi.org/10.3390/rs15123188

APA Style

Hu, J., Dong, X., Tian, W., Hu, C., Feng, K., & Lu, J. (2023). A Novel Optimization Strategy of Sidelobe Suppression for Pulse Compression Weather Radar. Remote Sensing, 15(12), 3188. https://doi.org/10.3390/rs15123188

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