Next Article in Journal
Evaluation of the Moso Bamboo Age Determination Based on Laser Echo Intensity
Next Article in Special Issue
Bubble Plume Target Detection Method of Multibeam Water Column Images Based on Bags of Visual Word Features
Previous Article in Journal
Function-Based Troposphere Tomography Technique for Optimal Downscaling of Precipitation
Previous Article in Special Issue
3D Sparse SAR Image Reconstruction Based on Cauchy Penalty and Convex Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unambiguous ISAR Imaging Method for Complex Maneuvering Group Targets

1
Air and Missile Defense College, Air Force Engineering University, Xi’an 710051, China
2
Science College, Armed Police Engineering University, Xi’an 710051, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(11), 2554; https://doi.org/10.3390/rs14112554
Submission received: 21 April 2022 / Revised: 18 May 2022 / Accepted: 25 May 2022 / Published: 26 May 2022
(This article belongs to the Special Issue Radar and Sonar Imaging and Processing Ⅲ)

Abstract

:
In inverse synthetic-aperture radar (ISAR) imaging, it is essential to deal with the Doppler ambiguity of group targets with complex maneuvers in order to avoid the bias of target position towards the actual value. Simultaneously, migration through resolution cell (MTRC) under the Doppler ambiguity is unable to be compensated for as a preprocessing. Traditional ISAR imaging methods for maneuvering targets, however, are undesirable to handle the severe deformation and defocusing in the imaging results induced by the Doppler ambiguity and MTRC. In this paper, we propose a novel and effective ISAR imaging method to improve the imaging quality by removing the Doppler ambiguity and compensating for the MTRC. Specifically, we first model the echo as a multi-component cubic phase signal (m-CPS) and design a high-order instantaneous autocorrelation function–generalized scaled Fourier transform (HIAF–GSCFT) to process the echo. This is to estimate the rotational parameters without MTRC compensation. Then, a maximum weighted contrast algorithm is used to remove the Doppler ambiguity, followed by reconstructing the echo. Compared with the existing method, the proposed method can accurately estimate the rotational parameters under the existing MTRCs and achieves a high-quality ISAR image for group targets, with complex maneuvers without Doppler ambiguity. Experiments of simulated and measured datasets validate its effectiveness and robustness for single target and group targets.

1. Introduction

In traditional ISAR imaging algorithms for maneuvering targets, the difficulty of achieving high-quality imaging results for complex maneuvering group targets is mainly caused by Doppler ambiguity because the Doppler frequency of partial scattering centers exceeds the limit of radar pulse repetition frequency (PRF). For group targets, the Doppler ambiguity generally exists among marginal subtargets, thus yielding the following two problems:
  • It results in the difficulty of correcting migration through resolution cell (MTRC) [1] for accurate target rotation parameters estimation. Without the Doppler ambiguity, MTRC can be compensated for by methods such as Keystone [1], and only the Doppler diffusion caused by the higher-order phase is considered. In general, the echo of maneuvering targets can be seen as multi-component linear frequency modulation (LFM) signals in the slow-time domain [2]. ISAR imaging methods for LFM signals can be categorized into two types: time–frequency analysis methods, such as Wigner–Ville distribution [2,3,4], Radon–Wigner transform [5,6], and LV’s distribution (LVD) [7,8], and parameter estimation methods, such as matched Fourier transform (MFT) [9,10] and Chirp–Fourier transform [11,12]. With improvement in the maneuverability of the air targets, the echo can be further modeled as m-CPS [13,14]. Most ISAR imaging methods for m-CPS are based on parameter estimation, including cubic phase function (CPF) [15,16,17,18], high-order ambiguity function (HAF) [19,20,21,22,23], and discrete polynomial phase transformation (DPT) [24]. In recent years, some scholars proposed novel methods based on other well-performing time–frequency distributions [25,26,27]. With the Doppler ambiguity, however, MTRCs can no longer be effectively compensated for by the aforementioned methods. This will further cause difficulty in estimating the echo parameters to compensate rotational motion of targets. Therefore, it is of critical value to design an effective parameter estimation algorithm without corrected MTRC.
  • Scattering centers with Doppler ambiguity will appear in an incorrect position [28,29], which highlights the necessity of removing the Doppler ambiguity. To our best knowledge, however, few studies focus on dealing with this issue. Dr. Huang and Dr. Zhang [30] presented a hypothesis that Doppler ambiguity removal can be carried out by sparse reconstruction. However, their proposed method has limited effectiveness for the m-CPS model.
In this study, we propose a novel ISAR imaging method to obtain an unambiguous ISAR imaging result of complex maneuvering group targets. In this method, the echo parameters are estimated to compensate for the rotational motion of targets by high-order instantaneous autocorrelation function–generalized scaled Fourier transform (HIAF–GSCFT), and a maximum weighted contrast (MWC) based algorithm is used to remove the Doppler ambiguity. The echo is then reconstructed to achieve unambiguous ISAR imaging. Experimental results prove that our method has better performance for complex maneuvering targets with MTRCs and Doppler ambiguity than traditional methods and is robust with different SNRs.

2. Materials and Methods

2.1. Signal Model

Generally, the movement of a target consists of translation and rotation. Since in this research, we mainly discuss rotational motion compensation, in this section, we directly analyze the echo after translational motion compensation. In this case, the turntable model of ISAR imaging is shown in Figure 1.
We set the transmitted signal as a linear frequency modulation signal
s T ( t ^ , t m ) = rect ( t ^ / T p ) e j 2 π f c t e j π μ t ^ 2
where T p is the pulse width, t ^ is the fast time, f c is the carrier frequency, μ is the chirp rate, t = t ^ + t m is the total time, and t m is the slow time, with an integer multiplying the pulse repetition period.
After the stretch processing [31], the echo of the scattering center p is
s p ( f r , t m ) = A p rect ( f r B ) e j 4 π c ( f c + f r ) Δ R p
where B is the bandwidth, and f r is the fast time–frequency defined by
f r = μ ( t ^ 2 R r e f / c )
where R r e f is the reference range, and c is the light speed.
Given that the target performs complex maneuvers, the rotation angle of the scattering center p is
θ ( t m ) = ω t m + 1 2 a t m 2 + 1 6 b t m 3
where ω is the angular velocity of rotation, a is the angular acceleration, and b is the second-order angular acceleration.
Based on planar wavefront approximation that is reasonable for ISAR imaging, Δ R p in Equation (2) can be expressed as
Δ R p ( t m ) = x p cos ( ω t m + 1 2 a t m 2 + 1 6 b t m 3 ) + y p sin ( ω t m + 1 2 a t m 2 + 1 6 b t m 3 )
where y p and x p represent the Y–X coordinates of p, respectively, in the coordinate system in Figure 1.
Here, Equation (5) can be approximated by Taylor series expansion as
Δ R p ( t m ) = x p + y p ω t m + 1 2 ( y p a x p ω 2 ) t m 2 + 1 6 ( y p b 3 x p a ω ω 3 ) t m 3
For the simplicity of notation, Equation (6) is written as
Δ R p ( t m ) = ϕ 1 , p + ϕ 2 , p t m + 1 2 ϕ 3 , p t m 2 + 1 6 ϕ 4 , p t m 3
where ϕ 1 , p = x p , ϕ 2 , p = y p ω , ϕ 3 , p = y p a x p ω 2 , and ϕ 4 , p = y p b 3 x p a ω ω 3 .
Considering the influence of Doppler ambiguity, the echo signal of complex maneuvering group targets containing N scattering centers is
S R ( f r , t m ) = p = 1 N s p ( f r , t m ) = p = 1 N A p rect ( f r B ) e j 4 π c ( f c + f r ) ( ϕ 1 , p + ( ϕ 2 , p + d a n p λ 2 P R F ) t m + 1 2 ϕ 3 , p t m 2 + 1 6 ϕ 4 , p t m 3 )
where d a n p [ 2 , 1 , 0 , 1 , 2 ... ] is the number of Doppler ambiguity [30], and λ is the wavelength.
It is worth noting that the Doppler ambiguity shifts the unique correspondence between the azimuth information of the scattering centers and its Doppler frequency. In other words, scattering centers with the same Doppler frequency may also have different MTRC slopes, as shown in Figure 2.
In Figure 2a, the position of the scattering center with MTRC is x = 18   m and y = 40   m . In contrast, the position of the scattering center in Figure 2b is x = 18   m and y = 60   m . Although with the same Doppler frequency, the MTRC slopes of these two scattering centers are completely different, making the traditional MTRC correction method fail to achieve an ideal result.

2.2. Parameters Estimation Based on HIAF–GSCFT

Due to Doppler ambiguity, the echo parameters estimation needs to be carried out when the MTRC cannot be compensated for. Thus, HIAF–GSCFT is designed to estimate the echo parameters. Inspired by the autocorrelation function [32], we define HIAF as
R ( f r , t m ) = S R ( f r , t m + ς 1 + ς 2 ) S R * ( f r , t m ς 1 + ς 2 ) S R * ( f r , t m + ς 1 ς 2 ) S R ( f r , t m ς 1 ς 2 )
where ς 1 and ς 2 are two delay constants. It is worth noting that f r and t m should be interpolated for the sake of the same number of sampling points. Combined with Equation (8), Equation (9) can be written as
R ( f r , t m ) = p = 1 N A p 4 e j 4 π c ( f c + f r ) ( 4 ϕ 3 , p ς 1 ς 2 + 4 ϕ 4 , p ς 1 ς 2 t m ) + R c r o s s ( f r , t m )
where R c r o s s ( f r , t m ) represents the cross-terms created by HIAF. In Equation (10), we ignore rect ( f r / B ) because in HIAF and phase separation, rect ( f r / B ) is invariant.
To ensure that GSCFT can effectively estimate the parameters, the phases only related to t m need to be eliminated. We define this process as phase separation (PS), which can be expressed as
PS [ R ( f r , t m ) ] = F T t r [ | IFT f r [ R ( f r , t m ) ] | ]
where t r is the corresponding time variable of f r after inverse Fourier transform, F T t r represents the Fourier transform in t r , IFT f r represents the inverse Fourier transform in f r , and * represents the conjugate operation of the result after IFT.
After PS, the echo phase is
PS [ R ( f r , t m ) ] = p = 1 N A p 4 e j 4 π f r c ( 4 ϕ 3 , p ς 1 ς 2 + 4 ϕ 4 , p ς 1 ς 2 t m ) + R P S c r o s s ( f r , t m )
where R P S c r o s s ( f r , t m ) is the cross-terms processed by PS.
According to [15] and [33], GSCFT is defined as
T ( f [ ξ Δ m e Υ m f ] , Δ m ) = Υ m f [ g ( Δ m ) e j 2 π ψ Δ m e Υ m f ] e j 2 π ( ξ Δ m e Υ m f f [ ξ Δ m e Υ m f ] ) d ( Δ m e Υ m f ) = g ( Δ m ) δ ( f [ ξ Δ m e Υ m f ] ψ ξ )
where ξ represents the zoom factor to avoid spectrum aliasing, f [ ξ Δ m e Υ m f ] is the scaled frequency domain with respect to Υ m f , e and f are constant power, g ( Δ m ) is a function of Δ m , and ψ is an empirical parameter.
According to Equation (12), the symbol in Equation (13) can be replaced by
Δ m = f r , Υ m = t m , e = f = 1
g ( Δ m ) e j 2 π ψ Δ m e Υ m f = PS [ R ( f r , t m ) ]
To this end, the result of GSCFT is
T ( f r , f [ ξ f r t m ] ) = p = 1 N A p 4 δ ( f [ ξ f r t m ] + 8 ϕ 4 , p ς 1 ς 2 c ξ ) e j 16 π c ϕ 3 , p ς 1 ς 2 f r + T c r o s s ( f r , f [ ξ f r t m ] )
where T c r o s s ( f r , f [ ξ f r t m ] ) is the cross-terms processed by GSCFT.
Clearly, the first term of Equation (16) on the right side achieves energy focusing in f [ ξ f r t m ] domain, and the positions of the peaks are only related to ϕ 4 , p . After inverse Fourier transform in f r domain, T ( f r , f [ ξ f r t m ] ) can be expressed as
IFT f r [ T ( f r , f [ ξ t r t m ] ) ] = p = 1 N A p 4 δ ( f [ ξ f r t m ] + 8 ϕ 4 , p ς 1 ς 2 c ξ ) δ ( t r 8 ϕ 3 , p ς 1 ς 2 c ) + T F T c r o s s ( t r , f [ ξ f r t m ] )
In Equation (17), the echo of each individual scattering center has formed a peak on the t r f [ ξ f r t m ] domain. Simultaneously, the cross-terms cannot achieve energy accumulation through HIAF–GSCFT (see a detailed analysis of the cross-term in Appendix A). Therefore, the peak detection to estimate ϕ 3 , p and ϕ 4 , p is defined by
{ t r m , f [ ξ f r t m ] m } = arg max | IFT f r [ T ( t r , f [ ξ f r t m ] ) ] |
ϕ ^ 3 , p m = c t r m 8 ς 1 ς 2 , ϕ ^ 4 , p m = c ξ f [ ξ f r t m ] m 8 ς 1 ς 2
where t r m and f [ ξ f r t m ] m represent the coordinates of the maximum value, respectively.
Moreover, GSCFT, ϕ ^ 3 , p m , and ϕ ^ 4 , p m can also be used to estimate ϕ 1 , p m and ϕ 2 , p m . We first compensate for the higher-order phase of the dominant scattering center and perform PS by
Φ ( f r , t m ) = e j 4 π c ( f r + f c ) ( 1 2 ϕ ^ 3 , p m t m 2 + 1 6 ϕ ^ 4 , p m t m 3 )
S R P S ( f r , t m ) = PS [ S R ( f r , t m ) Φ ( f r , t m ) ] = A p m e 4 π c ( ϕ 1 , p m f r + ϕ 2 , p m f r t m ) + S o t h e r s ( f r , t m )
where S o t h e r s ( f r , t m ) represents the echo components of nondominant scattering centers.
After GSCFT and Fourier transform, S R P S ( f r , t m ) becomes
T R ( t r , f [ ξ f r t m ] ) = IFT f r [ GSCFT [ S R ( f r , t m ) ] ] = A p m δ ( t r 2 ϕ 1 , p m c ) δ ( f [ ξ f r t m ] + 2 ϕ 2 , p m c ξ ) + T o t h e r s ( t r , f [ ξ f r t m ] )
where T o t h e r s ( t r , f [ ξ f r t m ] ) is S o t h e r s ( f r , t m ) processed by GSCFT and Fourier transform (detailed analysis is provided in Appendix A).
Similarly, peak detection is used to obtain ϕ ^ 1 , p m and ϕ ^ 2 , p m by
{ t r m , f [ ξ f r t m ] m } = arg max | T R ( t r , f [ ξ f r t m ] ) |
ϕ ^ 1 , p m = c t r m 2 , ϕ ^ 2 , p m = c ξ f [ ξ f r t m ] m 2
In order to ensure the accuracy of parameter estimation, we use a grid-searching manner for f [ ξ f r t m ] in [ P R F / 2 , P R F / 2 ] .

2.3. Doppler Ambiguity Removal Based on MWC

Given the Doppler ambiguity, the relationship between ϕ ^ 2 , p m and the actual value is
ϕ 2 , p m = ϕ ^ 2 , p m + d a n p m λ 2 P R F
where d a n p m is the Doppler ambiguity number of the dominant scattering center.
MWC is thus proposed to determine d a n p m . We first construct Equation (26) to multiply with the echo and transform the result into the t r f d domain
f 4 ( f r , t m , d a n ) = e j 4 π c f r ( ( ϕ ^ 2 , p m + d a n λ 2 P R F ) t m + 1 2 ϕ ^ 3 , p m t m 2 + 1 6 ϕ ^ 4 , p m t m 3 )
I d a ( t r , f d , d a n ) = | FFT 2 D [ S R ( f r , t m ) f 4 ( f r , t m , d a n ) ] |
where d a n [ 2 , 1 , 0 , 1 , 2 ... ] is used to determine d a n p m .
In Equation (27), when d a n = d a n p m , the MTRC will be completely compensated for. In this case, the image quality of I d a ( t r , f d , d a n ) is better than d a n d a n p m .
To prevent the interference caused by the nondominant scattering centers, we define the weighted contrast as
WC ( d a n ) = ( k ( t r , f d ) Ι d a ( t r , f d , d a n ) μ ) d t r d f d N r N d μ 2
Ι d a ( t r , f d , d a n ) = I d a ( t r + 2 ϕ ^ 1 , p m c , f d + 2 ϕ ^ 2 , p m λ , d a n )
μ = 1 N r N d k ( t r , f d ) I d a ( t r , f d , d a n ) d t r d f d
where k ( t r , f d ) is the weighting coefficient defined by
k ( t r , f d ) = 1 2 π σ ( t r ) e f d 2 2 σ ( t r ) 2
where σ ( t r ) is proportional to the standard deviation of the data in each range cell.
Now, the dominant scattering center is moved to the middle of I d a ( t r , f d , d a n ) and has a larger weight than the other scattering centers. Thus, d a n p m can be determined by the maximum of WC ( d a n ) by
{ d a ^ n p m } = max { WC ( d a n ) }
and the estimated Doppler frequency is
ϕ ^ 2 , p m = ϕ ^ 2 , p m + d a ^ n p m λ 2 P R F

2.4. ISAR Imaging Method Based on HIAF–GSCFT and Doppler Ambiguity Removal

In this section, a HIAF–GSCFT- and MWC-based ISAR imaging method is introduced for complex maneuvering group targets with the following procedure:
Step 1
Obtain S R ( f r , t m ) by processing the original echo signal using stretch processing;
Step 2
Estimate ϕ 3 , p m and ϕ 4 , p m using HIAF–GSCFT;
Step 3
Estimate ϕ 1 , p m and ϕ 2 , p m with PS, GSCFT, inverse Fourier transform, and the compensation function constructed by the estimated echo parameters;
Step 4
Remove the Doppler ambiguity by using the MWC and determine the actual Doppler frequency of the dominant scattering center;
Step 5
Estimate the amplitude of the dominant scattering center by the least-squares method as
A ^ p m = S R ( f r , t m ) Θ * ( f r , t m ) | Θ ( f r , t m ) | 2
Θ ( f r , t m ) = e j 4 π c ( f c + f r ) ( ϕ ^ 1 , p m + ϕ ^ 2 , p m t m + 1 2 ϕ ^ 3 , p m t m 2 + 1 6 ϕ ^ 4 , p m t m 3 )
Step 6
Delete the echo of the dominant scattering center from S R ( f r , t m ) as
S r m ( f r , t m ) = A ^ p e j 4 π c ( f c + f r ) ( ϕ ^ 1 , p m + ϕ ^ 2 , p m t m + 1 2 ϕ ^ 3 , p m t m 2 + 1 6 ϕ ^ 4 , p m t m 3 )
S R ( f r , t m ) = S R ( f r , t m ) S r m ( f r , t m )
Step 7
Repeat steps (2)–(6) till the residual energy of the echo reach the energy threshold;
Step 8
Reconstruct the echo by
S r e c ( f r , t m ) = p = 1 N A ^ p m e j 4 π c ϕ ^ 1 , p m f r e j 4 π λ ϕ ^ 2 , p m t m
Eventually, we interpolate the slow time grid and perform a 2D Fourier transform on S r e c ( f r , t m ) to obtain an ISAR imaging result as
S r e c ( t r , f d ) = p = 1 N A ^ p m δ ( t r + 2 ϕ ^ 1 , p m c ) δ ( f d + 2 ϕ ^ 2 , p m λ )
The flowchart of the proposed method is shown in Figure 3.

3. Results

In this section, we verify the proposed method from three aspects: the performance of parameter estimation, the capability of single target imaging, and the effectiveness of imaging for complex maneuvering group targets.

3.1. Validation of the Parameters Estimation Method Based on HIAF–GSCFT

As an m-CPS has two components, its signal model is similar to (8), and the rotational parameters are set as
Component 1: ϕ 1 , 1 = 20 , ϕ 2 , 1 = 30 , ϕ 3 , 1 = 15 , ϕ 4 , 1 = 25 ;
Component 2: ϕ 1 , 2 = 25 , ϕ 2 , 2 = 15 , ϕ 3 , 2 = 20 , ϕ 4 , 2 = 25 .
After HIAF–GSCFT for the m-CPS, the result is shown in Figure 4.
In Figure 4, CR and QCR represent the chirp rate and the quadratic chirp rate, respectively.
The two components achieve two individual energy accumulation. According to Figure 4, we estimate the parameters as ϕ ^ 3 , 1 = 14.91 , ϕ ^ 4 , 1 = 24.94 , ϕ ^ 3 , 2 = 19.88 , and ϕ ^ 4 , 2 = 24.94 , corresponding to relative estimation errors as 0.6%, 0.24%, 0.6%, and 0.24%, respectively. This is caused by discrete sampling intervals. Considering the small scale, however, these errors can be ignored.
Using the estimated parameters to compensate for the high-order phases of component 1 and component 2, respectively, and perform GSCFT, the results under two compensation cases are shown in Figure 5.
In Figure 5a, only component 1 achieves effective energy focusing because its high-order phases are compensated for correctly. Thus, we have ϕ ^ 1 , 1 = 19.88 and ϕ ^ 2 , 1 = 29.96 . Similarly, the frequency and the initial phase of component 2 can be estimated from Figure 5b as ϕ ^ 1 , 2 = 14.91 and ϕ ^ 2 , 2 = 24.94 respectively. The relative errors of these estimated parameters are 0.6%, 0.24%, 0.6%, and 0.24%, respectively.
From these results, it is fundamental that our proposed HIAF–GSCFT-based parameter estimation method has a strong capability of achieving accurate parameter estimation of m-CPS, with high-order coupling between f r and t m .

3.2. Validation of the Capability of Single Target Imaging

We adopted Yak-42 measured data to validate the imaging ability for a single target of the proposed method. The outline of Yak-42 is shown in Figure 6.
Parameters of the radar are 5.52 GHz carrier frequency, 100 Hz pulse repetition frequency, and 400 MHz bandwidth. We added white Gaussian noise (SNR = 5 dB and 0 dB) and MTRC to the data and compared our proposed method with range-Doppler (RD) [1], HAF–ICPF [21], and integrated parametric cubic phase function-reversing Wigner–Ville distribution (IPCPF–RWVD) [15], where RD is a classical ISAR imaging algorithm, HAF–ICPF is an effective imaging method for complex maneuvering target, and IPCPF–RWVD is a robust imaging method for a complex maneuvering target. The results of the above four methods are in Figure 7 and Figure 8.
Specifically, the RD method was severely defocused, especially on the body and wings of the aircraft. The HAF–ICPF method slightly improved the imaging quality, but the scattering centers of the body were still defocused. The IPCPF–RWVD method had better performance than HAF–ICPF because the IPCPF–RWVD method had a lower degree of nonlinearity. However, the result of the IPCPF–RWVD method was still defocused because MTRC affects the accuracy of parameter estimation. In contrast, our method effectively overcame the negative influence of MTRC in parameter estimation and obtained clear imaging results, as shown in Figure 7d and Figure 8d. In summary, this experiment proved the strong ability for high-quality ISAR imaging for a single target when MTRC is uncalibrated.

3.3. Validation of Imaging for Complex Maneuvering Group Targets

Next, we used synthetic data to validate the effectiveness of the proposed method for complex maneuvering group targets. As shown in Figure 9a, the simulated aircraft consists of 330 scattering centers, with a length of 18.67 m and a wingspan of 24 m. The amplitude of the scattering centers obeys the Rayleigh distribution, with σ = 0.5 . The group targets were composed of five simulated aircraft, and the spatial distribution is given in Figure 9b.
Parameters of radar and group targets are shown in Table 1 and Table 2, respectively.
For SNR = 5 dB, the ISAR imaging result obtained with the proposed method is shown in Figure 10f. In comparison, the results by RD, RD–Keystone, SPWVD, and HAF–ICPF, as well as the results reported in [30], are provided in Figure 10a–e, respectively. SPWVD is an effective imaging method for the LFM signal model, and the method in [30] is an imaging method for group targets.
When using the RD method, the Doppler ambiguity resulted in targets 1 and 5 appearing in the wrong positions, and the results were completely defocused. The result of RD–Keystone was even worse than RD because Keystone is invalid with MTRC under the Doppler ambiguity. SPWVD and HAF–ICPF methods improved the imaging results at different but limited levels of target defocusing; however, they did not have any effects on removing the Doppler ambiguity. The method in [30] was unable to remove the Doppler ambiguity of every scattering center because [30] is not ideal when applied to the m-CPS model. In contrast, our proposed method completely eliminated the Doppler ambiguity and defocus, as shown in Figure 10f. In summary, the proposed method effectively achieved the best quality imaging result for complex maneuvering group targets among all of the compared methods.
Finally, we validated the imaging quality of the proposed method under different SNRs. We adopted entropy and contrast to measure the quality of imaging results. After 50 Monte Carlo experiments, the SNR versus entropy curves and SNR versus contrast curves of the above six methods were drawn, which are shown in Figure 11.
We found that the imaging quality of RD and RD–Keystone were invalid. This is because the two methods were unable to compensate for the influence caused by the nonuniform rotation. SPWVD and HAF had better performance than the above two methods but were still not ideal because of the existence of Doppler ambiguity and MTRC. The performance of the method in [30] was the closest to that of the proposed method. However, the method was unable to correct the Doppler ambiguity of complex maneuvering group targets. In comparison, the imaging quality of our proposed method is significantly superior to all the other methods with different SNRs. This experiment effectively revealed the robustness of our method against different noise levels.

4. Conclusions

In this research, we presented a novel and effective ISAR imaging method for complex maneuvering group targets. In this method, the echo parameters were estimated via HIAF–GSCFT, and the MWC algorithm was used to remove the Doppler ambiguity. The proposed method effectively estimated rotational parameters in the presence of MTRC and improved the accuracy of Doppler ambiguity removal. Moreover, our method could be applied for complex maneuvering group targets. Experimental results demonstrated the effectiveness and robustness of this method with different SNRs. In future research, we plan to further improve the proposed method, by reducing the nonlinearity of the algorithm to improve the algorithm’s noise robustness and reduce the computational complexity.

Author Contributions

D.H. provided research ideas, F.L. completed simulation experiments and thesis writing, X.G. refined and polished the manuscript, and C.F. performed the final proofreading of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Postdoctoral Science Foundation of China under grant 2019M661508, Shaanxi Provincial Fund Youth Project of China under grant 2019JQ-497 and Aviation Science Fund of China under grant 201920096001 (Funder: Darong Huang). The APC was funded by 2019JQ-497.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the request of the funder.

Acknowledgments

This research was supported by the National Postdoctoral Science Foundation of China under Grant 2019M661508, Shaanxi Provincial Fund Youth Project of China under Grant 2019JQ-497, and Aviation Science Fund of China under Grant 201920096001.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this section, we conduct a detailed analysis of the changes in the cross-terms and nondominant echo components during the HIAF–GSCFT process. In Equation (17), we give three different forms of the components below.
R c r o s s 1 ( f r , t m ) = A p A q A r A s e j 4 π c ( f r + f c ) ( ( ϑ 2 , p q ϑ 2 , r s ) + ς 2 ( ϑ 3 , p q + ϑ 3 , r s ) + ς 2 2 ( ϑ 4 , p q ϑ 4 , r s ) 2 ) t m e j 4 π c ( f r + f c ) ( ( ϑ 3 , p q ϑ 3 , r s ) 2 + ς 2 ( ϑ 4 , p q + ϑ 4 , r s ) 2 ) t m 2 . e j 4 π c ( f r + f c ) ( ( ϑ 4 , p q ϑ 4 , r s ) 6 ) t m 3
R c r o s s 2 ( f r , t m ) = A p 2 A r A s e j 4 π c ( f r + f c ) ( ϑ 2 , r s + ς 2 ϑ 3 , r s + ς 2 2 ϑ 4 , p q 2 + 2 ς 1 ς 2 ϕ 4 , p 2 ς 1 ϕ 3 , p ) t m . e j 4 π c ( f r + f c ) ( ϑ 3 , r s + ς 2 ϑ 4 , r s 2 ς 2 ϕ 4 , p ) t m 2 e j 4 π c ( f r + f c ) ( ϑ 4 , r s 6 ) t m 3
R c r o s s 3 ( f r , t m ) = A p 2 A q 2 e j 4 π c ( f r + f c ) ( 2 ς 1 ( ϕ 3 , p ϕ 3 , q ) + 2 ς 1 ς 2 ( ϕ 4 , p + ϕ 4 , q ) ) t m e j 4 π c ( f r + f c ) ς 1 ( ϕ 4 , p ϕ 4 , q ) t m 2
where p, q, r, and s are four different scattering centers on the target. When they are used as the subscripts of a parameter, it represents that the parameter belongs to the echo of the corresponding scattering center. ϑ 2 , p q , ϑ 3 , p q and ϑ 4 , p q are written in Equation (A4), respectively.
{ ϑ 2 , p q = ϕ 2 , p ϕ 2 , q + ς 1 ( ϕ 3 , p + ϕ 3 , q ) + ς 1 2 2 ( ϕ 4 , p ϕ 4 , q ) ϑ 3 , p q = ϕ 3 , p ϕ 3 , q + ς 1 ( ϕ 4 , p + ϕ 4 , q ) 2 ϑ 4 , p q = ϕ 4 , p ϕ 4 , q 6
In Equation (21), s o t h e r s ( f r , t m ) is given by
s o t h e r s ( f r , t m ) = A q e j 4 π c ( f r + f c ) ( ϕ 1 , q + ϕ 2 , q t m ) e j 4 π c ( f r + f c ) ( 1 2 ( ϕ 3 , q ϕ ^ 3 , p ) t m 2 + 1 6 ( ϕ 4 , q ϕ ^ 4 , p ) t m 3 )
For two different scattering centers p and q in the imaging scene, their position relationship should at least satisfy x p x q 0 or y p y q 0 . Combining Equation (6), the parameters of p and q satisfy at least ϕ 1 , p ϕ 1 , q 0 or ϕ 2 , p ϕ 2 , q 0 and are subject to ϕ 3 , p ϕ 3 , q 0 and ϕ 4 , p ϕ 4 , q 0 . Furthermore, we can obtain ϑ 2 , p q 0 , ϑ 3 , p q 0 , ϑ 4 , p q 0 while ϑ 2 , p q ϑ 3 , p q ϑ 4 , p q . After PS processing, the signal models of the cross terms and nondominant echo components can be summarized into by
{ Ω 1 ( f r , t m ) = e j 4 π f r c ( η + κ t m + ρ t m 2 + χ t m 3 ) Ω 2 ( f r , t m ) = e j 4 π f r c ( η + κ t m + ρ t m 2 )
According to Equation (13), the prerequisite for using GSCFT to achieve energy focusing is that g must be related to Δ m only. In Equation (A6), however, f r is coupled with t m 2 and t m 3 . Moreover, the high-order phase of t m prevents the signal from forming a peak through inverse Fourier transform. In summary, the cross-term and the nondominant echo components cannot achieve energy focusing through GSCFT and Fourier transform, indicating that they have no effects on the parameter estimation.

References

  1. Xing, M.; Wu, R.; Lan, J.; Bao, Z. Migration Through Resolution Cell Compensation in ISAR Imaging. IEEE Geosci. Remote Sens. Lett. 2004, 1, 141–144. [Google Scholar] [CrossRef]
  2. Berizzi, F.; Mese, E.; Diani, M.; Martorella, M. High-resolution ISAR imaging of maneuvering targets by means of the range instantaneous Doppler technique: Modeling and performance analysis. IEEE Trans. Image Process. 2001, 10, 1880–1890. [Google Scholar] [CrossRef] [PubMed]
  3. Barbarossa, S. Analysis of multicomponent LFM signals by a combined Wigner-Hough transform. IEEE Trans. Signal Process. 1995, 43, 1511–1515. [Google Scholar] [CrossRef]
  4. Wang, Y.; Jiang, Y. ISAR Imaging of Maneuvering Target Based on the L-Class of Fourth-Order Complex-Lag PWVD. IEEE Trans. Geosci. Remote Sens. 2010, 48, 1518–1527. [Google Scholar] [CrossRef]
  5. Wood, J.; Barry, D. Radon transformation of time-frequency distributions for analysis of multicomponent signals. IEEE Trans. Signal Process. 1994, 42, 3166–3177. [Google Scholar] [CrossRef]
  6. Li, W.-C.; Wang, X.-S.; Wang, G.-Y. Scaled Radon-Wigner Transform Imaging and Scaling of Maneuvering Target. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 2043–2051. [Google Scholar] [CrossRef]
  7. Yang, T.L.; Yang, L.; BI, G. ISAR cross-range scaling algorithm based on LVD. In Proceedings of the 2016 39th International Conference on Telecommunications and Signal Processing (TSP), Vienna, Austria, 27–29 June 2016; pp. 643–664. [Google Scholar]
  8. Li, X.L.; Cui, G.L.; Kong, L.J. Fast Non-Searching Method for Maneuvering Target Detection and Motion Parameters Esti-mation. IEEE Trans. Signal Process. 2016, 64, 2232–2244. [Google Scholar] [CrossRef]
  9. Guo, B.; Sun, H.; Yin, W.; Zeng, H.; Chen, G.; Deng, L. ISAR Imaging Algorithm of Maneuvering Target Based on Matching Flourier Transform. In Proceedings of the 2018 IEEE 4th International Conference on Computer and Communications (ICCC), Chengdu, China, 7–10 December 2018; pp. 1710–1714. [Google Scholar] [CrossRef]
  10. Wang, Y.; Zhang, Q.; Zhao, B. ISAR imaging of target with complex motion based on novel approach for the parameters estimation of multi-component cubic phase signal. Multidimens. Syst. Signal Process. 2018, 29, 1285–1307. [Google Scholar] [CrossRef]
  11. Wang, B.; Xu, S.; Wu, W.; Hu, P.; Chen, Z. Adaptive ISAR Imaging of Maneuvering Targets Based on a Modified Fourier Transform. Sensors 2018, 18, 1370. [Google Scholar] [CrossRef] [Green Version]
  12. Zuo, L.; Wang, B. ISAR Imaging of Non-Uniform Rotating Targets Based on Optimized Matching Fourier Transform. IEEE Access 2020, 8, 64324–64330. [Google Scholar] [CrossRef]
  13. O’Shea, P.J. A Fast Algorithm for Estimating the Parameters of a Quadratic FM Signal. IEEE Trans. Signal Process. 2004, 52, 385–393. [Google Scholar] [CrossRef] [Green Version]
  14. O’Shea, P. Improving Polynomial Phase Parameter Estimation by Using Nonuniformly Spaced Signal Sample Methods. IEEE Trans. Signal Process. 2012, 60, 3405–3414. [Google Scholar] [CrossRef]
  15. Lv, Q.; Su, T.; He, X. An ISAR Imaging Algorithm for Nonuniformly Rotating Targets with Low SNR Based on Modified Bilinear Parameter Estimation of Cubic Phase Signal. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 3108–3124. [Google Scholar] [CrossRef]
  16. Zheng, J.; Liu, H.; Liu, Z.; Liu, Q.H. ISAR Imaging of Ship Targets Based on an Integrated Cubic Phase Bilinear Autocorrelation Function. Sensors 2017, 17, 498. [Google Scholar] [CrossRef] [Green Version]
  17. Zhu, L. Quadratic Frequency Modulation Signals Parameter Estimation Based on Product High Order Ambiguity Func-tion-Modified Integrated Cubic Phase Function. Information 2019, 10, 140. [Google Scholar] [CrossRef] [Green Version]
  18. Li, L.; Yan, L.; Li, D.; Liu, H.; Zhang, C. A Novel ISAR Imaging Method for Maneuvering Target Based on AM-QFM Model Under Low SNR Environment. IEEE Access 2019, 7, 140499–140512. [Google Scholar] [CrossRef]
  19. Wang, Y.; Kang, J.; Jiang, Y. ISAR Imaging of Maneuvering Target Based on the Local Polynomial Wigner Distribution and Integrated High-Order Ambiguity Function for Cubic Phase Signal Model. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 2971–2991. [Google Scholar] [CrossRef]
  20. Huang, P.; Liao, G.; Yang, Z.; Xia, X.-G.; Ma, J.; Zheng, J. Ground Maneuvering Target Imaging and High-Order Motion Parameter Estimation Based on Second-Order Keystone and Generalized Hough-HAF Transform. IEEE Trans. Geosci. Remote Sens. 2016, 55, 320–335. [Google Scholar] [CrossRef]
  21. Wang, Y. Inverse synthetic aperture radar imaging of manoeuvring target based on range-instantaneous-Doppler and range-instantaneous-chirp-rate algorithms. IET Radar Sonar Navig. 2013, 6, 921–928. [Google Scholar] [CrossRef]
  22. Jing, F.L.; Si, W.J.; Wang, Y. Parameter estimation of polynomial-phase signal using the hybrid LvHAF. In Proceedings of the 2017 Progress in Electromagnetics Research Symposium—Fall (PIERS—FALL), Singapore, 19–22 November 2017; pp. 417–426. [Google Scholar]
  23. Zhang, J.; Su, T.; Zheng, J.; He, X. Parameter estimation of CFM signals based on MICPF-HAF. Electron. Lett. 2018, 54, 456–458. [Google Scholar] [CrossRef]
  24. Wang, Y.; Huang, X.; Cao, R. Novel Approach for ISAR Cross-Range Scaling Based on the Multidelay Discrete Polynomial- Phase Transform Combined with Keystone Transform. IEEE Trans. Geosci. Remote Sens. 2019, 58, 1221–1231. [Google Scholar] [CrossRef]
  25. Jin, K.; Lai, T.; Wang, Y.; Xing, X.; Zhao, Y. Parameter estimation of quadratic frequency modulated signal based on three-dimensional scaled Fourier transform. IET Radar, Sonar Navig. 2019, 13, 1689–1696. [Google Scholar] [CrossRef]
  26. Qian, J.; Huang, S.; Wang, L.; Bi, G.; Yang, X. Super-Resolution ISAR Imaging for Maneuvering Target Based on Deep-Learning-Assisted Time–Frequency Analysis. IEEE Trans. Geosci. Remote Sens. 2021, 60, 1–14. [Google Scholar] [CrossRef]
  27. Huang, P.; Xia, X.-G.; Zhan, M.; Liu, X.; Liao, G.; Jiang, X. ISAR Imaging of a Maneuvering Target Based on Parameter Estimation of Multicomponent Cubic Phase Signals. IEEE Trans. Geosci. Remote Sens. 2021, 60, 1–18. [Google Scholar] [CrossRef]
  28. Pan, X.; Wang, W.; Feng, D.; Liu, Y.; Fu, Q.; Wang, G. On deception jamming for countering bistatic ISAR based on sub-Nyquist sampling. IET Radar Sonar Navig. 2014, 8, 173–179. [Google Scholar] [CrossRef]
  29. Pan, X.; Liu, J.; Chen, J.; Xie, Q.; Ai, X. Sub-Nyquist Sampling Jamming Against Chirp-ISAR with CS-D Range Compression. IEEE Sensors J. 2017, 18, 1140–1149. [Google Scholar] [CrossRef]
  30. Huang, D.; Zhang, L.; Xing, M.; Bao, Z. Doppler ambiguity removal and ISAR imaging of group targets with sparse decomposition. IET Radar Sonar Navig. 2016, 10, 1711–1719. [Google Scholar] [CrossRef]
  31. Caputi, W.J. Stretch: A Time-Transformation Technique. IEEE Trans. Aerosp. Electron. Syst. 1971, AES-7, 269–278. [Google Scholar] [CrossRef]
  32. Bao, Z.; Xing, M.D.; Wang, T. Radar Imaging Technology; Publishing House of Electronics Industry: Beijing, China, 2005; pp. 264–266. [Google Scholar]
  33. Zheng, J.; Su, T.; Zhu, W.; Zhang, L.; Liu, Z.; Liu, Q.H. ISAR Imaging of Nonuniformly Rotating Target Based on a Fast Parameter Estimation Algorithm of Cubic Phase Signal. IEEE Trans. Geosci. Remote Sens. 2015, 53, 4727–4740. [Google Scholar] [CrossRef]
Figure 1. Turntable model of ISAR imaging.
Figure 1. Turntable model of ISAR imaging.
Remotesensing 14 02554 g001
Figure 2. Comparison of range profiles: (a) range profile of unambiguous scattering centers; (b) range profile of ambiguous scattering centers.
Figure 2. Comparison of range profiles: (a) range profile of unambiguous scattering centers; (b) range profile of ambiguous scattering centers.
Remotesensing 14 02554 g002
Figure 3. Flowchart of the proposed method.
Figure 3. Flowchart of the proposed method.
Remotesensing 14 02554 g003
Figure 4. Result of HIAF–GSCFT.
Figure 4. Result of HIAF–GSCFT.
Remotesensing 14 02554 g004
Figure 5. Results of GSCFT: (a) using the estimated parameters of component 1 to compensate for the echo; (b) using the estimated parameters of component 2 to compensate for the echo.
Figure 5. Results of GSCFT: (a) using the estimated parameters of component 1 to compensate for the echo; (b) using the estimated parameters of component 2 to compensate for the echo.
Remotesensing 14 02554 g005
Figure 6. The outline of Yak-42.
Figure 6. The outline of Yak-42.
Remotesensing 14 02554 g006
Figure 7. ISAR imaging results for Yak-42 by different methods under SNR = 5 dB: (a) the result of RD method; (b) the result of HAF–ICPF method; (c) the result of IPCPF–RWVD method; (d) the result of the proposed method.
Figure 7. ISAR imaging results for Yak-42 by different methods under SNR = 5 dB: (a) the result of RD method; (b) the result of HAF–ICPF method; (c) the result of IPCPF–RWVD method; (d) the result of the proposed method.
Remotesensing 14 02554 g007aRemotesensing 14 02554 g007b
Figure 8. ISAR imaging results for Yak-42 by different methods under SNR = 0 dB: (a) the result of RD method; (b) the result of HAF–ICPF method; (c) the result of IPCPF–RWVD method; (d) the result of the proposed method.
Figure 8. ISAR imaging results for Yak-42 by different methods under SNR = 0 dB: (a) the result of RD method; (b) the result of HAF–ICPF method; (c) the result of IPCPF–RWVD method; (d) the result of the proposed method.
Remotesensing 14 02554 g008
Figure 9. Simulated group targets: (a) simulated aircraft; (b) spatial distribution of simulated aircraft.
Figure 9. Simulated group targets: (a) simulated aircraft; (b) spatial distribution of simulated aircraft.
Remotesensing 14 02554 g009
Figure 10. ISAR imaging results obtained from different methods: (a) RD; (b) RD–Keystone; (c) SPWVD; (d) HAF–ICPF; (e) the method of [30]; (f) proposed method.
Figure 10. ISAR imaging results obtained from different methods: (a) RD; (b) RD–Keystone; (c) SPWVD; (d) HAF–ICPF; (e) the method of [30]; (f) proposed method.
Remotesensing 14 02554 g010aRemotesensing 14 02554 g010b
Figure 11. SNR versus entropy curve and SNR versus contrast curve of the six methods: (a) SNR versus entropy curves; (b) SNR versus contrast curves [30].
Figure 11. SNR versus entropy curve and SNR versus contrast curve of the six methods: (a) SNR versus entropy curves; (b) SNR versus contrast curves [30].
Remotesensing 14 02554 g011
Table 1. Parameters of radar.
Table 1. Parameters of radar.
Frequency9.6 GHz
Bandwidth2 GHz
PRF200 Hz
Pulse width60 μs
Sampling rate13 MHz
Sampling number512
Table 2. Parameters of group targets.
Table 2. Parameters of group targets.
Target CodeTarget 1, 5Target 2, 3, 4
Velocity/m/s453.78471.24
A = Acceleration/m/s269.8174.17
Jerk/ m/s343.6347.99
Rotation angle/◦5.776.04
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, F.; Huang, D.; Guo, X.; Feng, C. Unambiguous ISAR Imaging Method for Complex Maneuvering Group Targets. Remote Sens. 2022, 14, 2554. https://doi.org/10.3390/rs14112554

AMA Style

Liu F, Huang D, Guo X, Feng C. Unambiguous ISAR Imaging Method for Complex Maneuvering Group Targets. Remote Sensing. 2022; 14(11):2554. https://doi.org/10.3390/rs14112554

Chicago/Turabian Style

Liu, Fengkai, Darong Huang, Xinrong Guo, and Cunqian Feng. 2022. "Unambiguous ISAR Imaging Method for Complex Maneuvering Group Targets" Remote Sensing 14, no. 11: 2554. https://doi.org/10.3390/rs14112554

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