Abstract
In order to solve the problem of high-resolution ISAR imaging under the condition of finite pulses, an improved smoothed L0 norm (SL0) sparse signal reconstruction ISAR imaging algorithm is proposed. Firstly, the ISAR imaging is transformed into the optimization problem of minimum L0 norm. Secondly, a single-loop structure is used instead of two loop layers in SL0 algorithm which increases the searching density of variable parameter to ensure the recovery accuracy. Finally, the compared step is added to ensure the optimization solution along the steepest descent gradient direction. The experimental results show that the proposed algorithm has better imaging effect.
1. Introduction
Inverse synthetic-aperture radar (ISAR) has the characteristics of all-weather, all-time, long-range, and high-resolution which is widely used in military and civil fields, such as to detect, locate, track, and estimate the parameters of the target [1–4]. The range and azimuth resolution of ISAR imaging are related to radar system bandwidth and imaging accumulation angle, respectively. The range resolution is usually improved by increasing the bandwidth of the transmitted signal, and the azimuth resolution is improved by increasing the observation angle range. In a long observation time, the scattering characteristics of noncooperative moving targets change greatly, and there may be a large Doppler time-varying, which is not conducive to imaging and makes it difficult to improve the azimuth resolution. So implementing imaging in a short time duration is meaningful.
In recent years, compressive sensing (CS) has become very popular in signal processing [5–8]. It is a new signal processing theory proposed by Donoho which is developed from the theory of signal sparse decomposition and approximation. The theory shows that as long as the signal is compressible or sparse in a certain transform domain, the original signal can be reconstructed with a small number of sampled signals. It is a technique proposed to improve signal separation ability using a prior sparse property information of the signal.
In radar imaging, the targets are generally composed of strong scattering points, which are sparsely distributed in the imaging region. So radar signal processing based on compressive sensing theory has been widely studied [9–12], such as SAR imaging [13], ISAR imaging [14, 15], and MIMO radar imaging [16]. The CS theory is not sensitive to data loss, which can effectively improve the imaging problem of radar imaging system in the case of data loss and improve the signal processing ability. The Bayesian sparse signal reconstruction ISAR imaging algorithm based on hierarchical prior model correlation prior is proposed in [17], which can obtain better focusing effect in noisy environment. However, the algorithm is computational complexity. By using the sparsity of ISAR target in azimuth direction, the target image can be obtained by solving the L1 norm minimization in [18]. Combination of local sparsity constraint and nonlocal total variation is discussed in [19]. A two-dimensional structure mode-coupled Bayesian (PC-SBL) ISAR imaging algorithm is proposed in [20], which can realize high-resolution ISAR imaging, but the algorithm has a large amount of computation, which is not conducive to real-time imaging.
Mohimani et al. proposed the smoothed L0 norm sparse signal recovery algorithm by using the continuous Gaussian function sequence with parameters as smoothing function to approximate the minimum L0 norm [21]. The steepest descent method and gradient projection principle by a double-layer loop is used to obtain the optimal solution. The outer layer establishes a sequence from large to small, and the inner layer uses the steepest descent method to obtain the approximate solution. The reconstruction speed of the algorithm is 2-3 times faster than basis pursuit algorithm when it has the same accuracy.
In order to achieve high-resolution ISAR imaging under the condition of finite pulses, an ISAR imaging algorithm with single-loop structure compared to SL0 sparse signal recovery algorithm is proposed. The imaging problem is transformed to a minimum L0 norm optimization problem to achieve high-resolution ISAR imaging. The steepest descent method should reduce the cost function at every step, but it is not really along the descending direction in the actual solving process. Therefore, in the above algorithm, the step of checking whether to descend is added in each iteration. If not, the midpoint of the previous point and the current point is taken for correction to ensure that the search direction is along the steepest descent direction. By using the improved algorithm, ISAR imaging is more intensive with limited pulse numbers. Real data ISAR images obtained using the proposed method is competitive to the several popular methods.
The paper is organized as follows. “ISAR Imaging Model” introduces the ISAR sparse imaging model. In “ISAR Imaging Algorithm,” the proposed reconstruction algorithm is introduced in detail. Simulation and real data ISAR imaging results are presented in “Experimental Results.” Finally, “Conclusion” provides the conclusion.
2. ISAR Imaging Model
Assuming that the target is located in far field, the radar transmits linear frequency modulization signal as where is the center frequency, is the slope, and is the pulse period. Then, the complex echo signal of scattering point is where is time delay and is signal amplitude. After range compression, the signal can be expressed as
After the envelope compensation, the target rotates around the coordinate axis at a constant angular velocity . Assuming the initial position of the scattering point is , the time delay is mainly determined by . After the Taylor expansion, is expressed as
Then, , where is light speed. The echo signal of range unite is where is the wavelength. After discretization, the above formula can be written as where , , and is the sampling interval. The received signal is written as a matrix; then where is the vector composed of scattering coefficient and is sparse matrix.
To solve , the following sparse optimization strategy is used: where is a small constant which is related to the noise variance. The ISAR imaging problem is transformed into the reconstruction problem of formula (8).
3. ISAR Imaging Algorithm
To solve , the smoothed L0 norm algorithm is proposed in [21]. The Gauss sequence is used to approach the L0 norm. The minimum L0 norm can be obtained by gradually reducing the value of control parameter . A two-layer method was proposed to solve the sparse signal recovery problem in the algorithm. In order to further improve the reconstruction effect of sparse signal, the one-order negative exponential function sequence is proposed as smoothed function to approach L0 norm in the paper. When , approximates L1 norm. When , approximates L0 norm. Therefore, the sparse solution can be searched with high probability at the beginning of iteration.
In the double-loop structure of SL0 algorithm, the solution of is only as the initial value of , so it is not necessary to solve the exact solution of in the inner loop. Therefore, the sparse signal reconstruction algorithm in this paper cancels the inner loop step and adopts the single loop instead of the double-loop structure in SL0 algorithm [52], the interval between control parameters and is reduced, only one gradient descent method is used to solve the minimum value of the cost function, and then the solution of is given as the initial value of . Due to the solution obtained by the steepest descent method is not necessarily in the feasible solution region, it is necessary to project the solution obtained by each iteration into the feasible solution region using the projection method. The proposed algorithm is called the improved single-loop smoothed L0 norm (ISSL0) algorithm. The ISSL0 algorithm can be expressed as follows:
Initialization is as follows: (1)Let be equal to the minimum L2 norm solution of , obtained by (2)Choose a suitable decreasing sequence for ,
for (i)Let and (ii)Minimize the function on the feasible set (iii)Let be the gradient of (iv)For , let (v)If , project back into the feasible set:(vi)Compare step
While update
If (vii).(3)Final solution is .
The steepest descent method should reduce the cost function at every step, but it is not the descending direction in the actual solution process. Therefore, in the above algorithm, the step of checking whether to descend is added in the process of each iteration. If not, the midpoint of the previous point and the current point are taken for correction to ensure that the search direction is along the steepest descent direction. By using the single loop and compared step, the algorithm ensures the reconstruction accuracy and the computation amount cannot increase. The choice of step size factor in the steepest descent method is important. For a large step size, it may not converge, but for a very small step size, the computation efficiency is low. The step size should decrease when the searching point approaches the minimum solution. For search step factor, because the search point is far away from the minimum value point at the initial value, a larger step size is selected. When the search point is gradually close to the minimum value, the search step size should be gradually reduced. Regulatory factors are used to adjust the step; let .
4. Experimental Results
4.1. Simulation 1: One-Dimensional Sparse Signal Recovery
In order to verify the reconstruction performance of the algorithm, the paper tests the algorithm through the MATLAB processing platform. The signal model with noise is , where is constructed by selecting its components from and is the sparse signal, the nonzero elements of sparse signal consist of random signal . For the SL0 method, the numbers of outer loop and inner loop are 20 and 10, respectively. For the ISSL0 algorithm, the loop number is 200. For the ISSL0 algorithm, the loop number is 200. The experiment was implemented 100 times. The MSE is defined as . Changing the sparsity ratio of the signal, the curve of reconstruction time, correction rate, and MSE for different algorithms of sparse signal are shown in Figures 1–3. As can be seen from Figure 1, the ISSL0 algorithm takes more time than the OMP [22] algorithm and SL0 algorithm and shorter than the Laplace [23] algorithm. From Figures 2 and 3, we can see that the performances of ISSL0 algorithm are competitive with other algorithms.
4.2. Simulation 2: ISAR Imaging Using Real Data
A set of real data of the Yak-42 plane is used to demonstrate the performance of the proposed ISAR imaging algorithm. The related parameter descriptions of the radar data are listed as follows: the carrier frequency is 10 GHz with signal bandwidth of 400 MHz, and the range resolution is 0.375 m. The center carrier frequency is 5.52 GHz, and the pulse repetition frequency is 50 Hz. 256 echo pulses are selected as experimental data. 16, 32, and 64 pulses are implemented. For the SL0 method, the numbers of outer loop and inner loop are 20 and 10, respectively. For the ISSL0 algorithm, the loop number is 200. The simulation results are compared visually and quantitatively to those images obtained by OMP, Laplace, SL0, and ISSL0 algorithms in Figures 4–6. It is noticeable that more amounts of pulses generally lead to better imagery results. When the number of pulses is 64, the four sparse signal reconstruction algorithms are better. For the SL0 algorithm, when 16 and 32 pulses are used, some information of scattered points is missing. When the number of pulses is small, a large number of false scattering points appear in the OMP and Laplace algorithms. Regardless of the number of pulses, the ISSL0 algorithm keeps the aircraft in a good shape without false scattered points, and the geometric representation of the target is clear.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
5. Conclusion
One improved sparse signal recovery ISAR imaging algorithm is proposed to ISAR imaging. A single-loop structure is proposed to solve the optimal solution. The revised step is added to ensure the searching direction is decrease compared to the traditional smoothed L0 norm recovery algorithm. The experiment results verify that the proposed algorithm can improve ISAR imaging quality.
Data Availability
All data sources are reliable and the corresponding author can be contacted.
Conflicts of Interest
The author declares no conflicts of interest.
Acknowledgments
This work was partially supported by the National Science Foundation of China (no. 61471191) and the Natural Science Research Project of Jiangsu Colleges and Universities (19KJB510043).