1. Introduction
Inverse Synthetic Aperture Radar (ISAR) is widely utilized to generate high-resolution images of military or civilian targets under all-day, all-weather conditions [
1,
2,
3]. The pulse-Doppler (PD) radar typically transmits linear frequency-modulated (LFM) signals with a large bandwidth to achieve high range resolution. The high azimuth resolution is obtained through the synthetic aperture formed by the relative motion between the target and the radar. Since the imaging targets are generally non-cooperative, such relative motion may introduce undesirable range migration. Therefore, many methods have been proposed in the past few decades to compensate for range migration, thereby achieving excellent performance [
4,
5].
However, in the modern situation (e.g., ships moving in a formation), multiple targets may appear within the same radar antenna beam, impairing the mono-target assumption. Additionally, micro-motion components are widespread among various targets. Induced by the rotating antennas or propellers, the so-called micro-Doppler effect may cause interference stripes in the imaging result, significantly degrading the image quality. Under the multi-target assumption, different targets may be connected by these interference stripes, making it difficult to achieve focused images. Therefore, separating and imaging multiple targets under the micro-Doppler connection is a great challenge.
The existing multi-target ISAR imaging methods typically fall into three categories: time-frequency (TF)-based [
6,
7,
8,
9,
10,
11,
12,
13,
14], segmentation-based [
15,
16,
17,
18,
19], and parameter-based methods [
20,
21,
22,
23,
24].
TF-based methods first extract Doppler variation along the cross-range direction. Then, different targets’ TF curves are filtered in the TF domain according to the estimated variation. Based on the assumption that different targets possess their own chirp rates, Li et al. [
6] first transformed the signal into the fractional Fourier transform domain. Then, the CLEAN algorithm was utilized to separate different targets. Zhang et al. [
7] integrated the Radon transform with the Reassigned Smoothed Pseudo Wigner–Ville Distribution (RSPWVD), achieving significantly improved separation. Xiao et al. [
10] first employed the Short-Time Fourier transform (STFT) on the overlapped echo. Then, different TF curves were filtered using a mask suitable for different chirp rates. As a result, they applied the inverse STFT on these curves to reconstruct the finally separated echoes. Nevertheless, the separation efficiency usually depends on the estimation accuracy. Moreover, many TF methods, such as Cohen’s classes, suffer from a serious cross term, severely impairing their effectiveness.
Segmentation-based methods first separate different targets in the image or the high-resolution range profiles (HRRPs) domain. For the image domain, separated images are transformed back into the HRRPs domain, where motion compensation and imaging are further performed. Bai et al. [
15] first simultaneously compensated for the range migration to form a “bulk” image of targets. Then, different targets were segmented through clustering methods and normalized cuts. Finally, these separated images were transformed back to the HRRPs domain to achieve the refined compensation, accordingly obtaining the well-focused targets. For the HRRPs domain, different targets are first separated in the HRRPs domain to meet the single-target assumption. Following that, the conventional imaging algorithms are then applied to the separated HRRPs. Kong et al. [
16] employed the modified envelope correlation method and the Radon transform to align the HRRPs. Then, they separated different targets and further employed the Discrete chirp-Fourier transform (DCFT) to achieve final imaging. However, overlapping HRRPs and images limit the application of the segmentation-based and the HRRPs-based algorithms.
As for parameter-based methods, they typically regard the translational motion of each target as a polynomial. Then, coefficient parameters are searched through various iterative methods. Liu et al. [
20] utilized the particle swarm optimization (PSO) to search coefficients, consequently compensating for the range curvature. Then, a modified CLEAN technique was applied to extract different targets. Nevertheless, the parametric methods are usually sensitive with respect to the searching step as they are all based on the parameter space.
The aforementioned methods achieve outstanding results within their specific assumptions. However, the micro-Doppler effect is not considered in their assumptions. The micro-Doppler effect, caused by micro-motion components, usually creates interference stripes along the cross-range dimension in imaging results. Over the past few decades, numerous algorithms have been developed to eliminate this interference [
25,
26]. Recently, micro-Doppler removal methods based on Robust Principal Component Analysis (RPCA) have received increasing attention [
27]. Zhou et al. [
28] utilized RPCA to decompose the target into rigid body and micro-motion components. Zhang et al. [
29] extended RPCA based on the adjacent similarity among HRRPs, accordingly achieving excellent separation.
Inspired by these methods, this paper proposes a multiple-target ISAR imaging method with the removal of micro-motion connections based on the integration of joint constraints. First, micro-motion components are roughly filtered out in the image domain, and then transformed back into the HRRPs domain. Then, we apply the Radon transform to the micro-motion components to roughly estimate the target’s range migration, further forming the coarse range migration transformation. Following that, we construct a constrained cost function consisting of low-rank and sparsity. The cost function is solved using the linearized alternating direction method with adaptive penalty (LADMAP) [
30] to ensure convergence, further achieving rough separation of different targets and micro-motion components. The coarse separated targets are then used to estimate a more accurate range migration transformation through a correlation-based method. Subsequently, we construct a tri-constrained cost function consisting of low-rank, sparsity, and local similarity. Finally, different targets and their micro-motion components are entirely separated. Additionally, we analyze the periodicity of the micro-motion components to ensure complete separation.
The main contributions of this paper are summarized as follows:
We formulate a constrained optimization, consisting of the low-rank and the sparsity. Such optimization is solved by the LAMAP method to achieve the coarse separation. Consequently, the range migration is accurately estimated by a correlation-based method.
We establish a triple-constrained problem especially for the multi-target scenario with micro-Doppler connections. Specifically, we formulate a constrained optimization, consisting of the low-rank, the sparsity, and the Laplacian constraint. Such optimization is solved by the LAMAP method to achieve complete separation.
We modified the mono-target assumption and extended it to a multi-target scenario with micro-Doppler connections. Under this framework, the micro-Doppler effect is eliminated, and the targets are separated. In the proposed framework, different norms correspond to different signal characteristics, which effectively demonstrates the applicability of our approach.
The remainder of this article is organized as follows:
Section 2 models the signal model of multiple targets under the micro-Doppler scenarios.
Section 3 analyzes the signal characteristics.
Section 4 gives a detailed description of the proposed procedure.
Section 5 evaluates the performances of the proposed method through the simulated data. Finally, the conclusions are given in
Section 6.
2. Imaging Model
In this section, we will briefly introduce the multi-target ISAR signal model with micro-Doppler connections.
The imaging geometry model of multiple targets is presented in
Figure 1. The PD radar is located at the origin point
O of the
coordinate system. Blue lines represent the radar antenna beam. Moving in various directions, different targets exist within the same antenna beam, resulting in multiple range migrations along the line of sight (LOS) of the radar. Collectively referred to as translational motion, these range migrations may result in severe defocusing in imaging results. Additionally, the target’s orientation rotates by a certain angle around its rotational center during the imaging interval. This rotation generates a synthetic aperture, leading to high resolution along the cross-range dimension of the imaging results.
In ISAR imaging, a target is generally regarded as a cluster of scatterers. Therefore, the signal reflected from a target is considered as the summation of the signals reflected by each scatterer. The PD radar transmits linear frequency-modulated (LFM) signals with a large bandwidth. Hence, the signal received by the radar receiver can be written as
where
and
represent the fast time and the slow time, respectively.
Z denotes the total number of scatterers.
represents the amplitude of the
i-th scatterer.
,
, and
represent the basic frequency, the chirp rate, and the pulse width of the signal, respectively.
denotes the instantaneous distance between the
i-th scatterer and the radar. Typically, dechirp processing is employed on the reflected signal to reduce the bandwidth. Therefore, the echo after dechirp processing can be derived as
where
denotes the reference range of the dechirp signal.
is the difference between the instantaneous distance and the reference range of the
i-th scatterer.
is written as
in the following formula as different scatterers share the same reference range. Applying the fast Fourier transform along the
direction to the (
2), the HRRPs are written as
where
represents the range frequency of the target. As previously mentioned, the motion of a target consists of the translational motion and the rotational motion. The overall motion
and its translational and rotational components are written as
where
and
represent the rotational motion and the translational motion, respectively.
and
denote the abscissa and the ordinate of the
i-th scatterer in the target coordinate.
is the rotational angle of the target.
represents the initial distance between the target center and the radar.
v,
a, and
represent the target’s velocity, acceleration, and
n-th order acceleration, respectively. For most non-cooperative targets, the imaging interval is relatively short, leading to small rotation angles. Therefore, we obtain
The rotational motion represented by (
7) leads to high resolution along the cross-range dimension, while the translational motion represented by (
8) results in severe defocusing. In addition, this paper focuses on the separation of multiple targets with different motion characteristics within a smaller range. Under the multi-target assumption,
belongs to its own coordinate systems. Therefore, conventional range alignment methods fail to achieve unified motion compensation for multiple targets.
In addition to rigid body parts, most non-cooperative targets contain micro-motion components that rotate at high speeds (e.g., engines in aircraft, propellers, and antennas on ships). The motion of these micro-motion components is usually superimposed on the rotation of the rigid body. The Doppler shifts caused by the micro-motion components are also superimposed on the echoes of the rigid body. This Doppler shift induced by micro-motion components is known as the micro-Doppler effect. As illustrated in
Figure 2, we place the rotational center of the target’s rigid body at the origin of the XOY coordinate system. The rotational center of the micro-motion components is set at the origin of the
coordinate system, which serves as the motion coordinate system for the micro-motion components themselves. Additionally, the origin of the
coordinate system is set on the rigid body scatterer, meaning that this coordinate system rotates along with the rigid body.
Suppose that the scatterer
on the rotating antenna rotates around its own rotational center with a rotational angle of
and a rotational shaft of
. The initial angle of the micro-motion components is
. The instantaneous distance between the scatterer and the radar is given as
According to the aforementioned approximation, we obtain
It is seen that the cosine term exhibits multiple periods within the imaging interval as the variation of is much greater than that of . These periods result in interference stripes in the imaging results. In the next section, we will describe the procedure of simultaneously separating targets and micro-Doppler interferences in detail.
3. Analysis of the Signal Characteristics
In this section, we will analyze the unique characteristics of the micro-motion components and the rigid body signal. First, the low-rank property of the rigid body in the HRRPs domain is theoretically determined by the number of the scattering points. Assuming that the translational motion of the target has been fully compensated, performing a fast Fourier transform (FFT) on a single scatterer in (
3) along the
direction yields
which can be discretized as
where
n and
p are the index of the range cell and the pulse, respectively. Suppose that the size of the HRRP matrix is
, (
12) can be written as
where
and
with
and
, respectively. For a target composed of
scattering points, its echo can be expressed as
Based on the properties of rank, the maximum rank of this HRRP matrix does not exceed . Thus, an echo with well-compensated motion can be considered low rank. However, undesirable motion disrupts this characteristic. Therefore, we leverage the inherent low-rank nature of the echo to separate targets with different motions.
For the sparsity of rigid body signals in the image domain and micro-motion components in the HRRPs domain, the sparsity of a rigid body signal in the image domain is theoretically determined by the number of scatterers. In contrast, micro-motion components typically occupy several range cells and form stripes across the entire Doppler domain in the image domain, resulting in lower sparsity. Conversely, the signal of the micro-motion components in the HRRPs domain occupies only a few range cells, demonstrating strong sparsity.
We conducted an experiment to verify such characteristics.
Figure 3 shows the range profiles and the imaging results of the rigid body and the micro-motion components.
Figure 3a,b present the range profiles and the imaging results of the ridge body, while
Figure 3c,d show those for the micro-motion components. In the HRRPs domain, the rate of non-zero pixels for the rigid body and micro-motion components is 0.22 and 0.02, respectively, confirming the sparsity of the micro-motion components in the HRRPs domain. Moreover, the absolute value of the rate of non-zero pixels may be dependent on the selected range after pulse compression, while the relative value depends on the size ratio between the rigid body and the micro-motion components. Therefore, micro-motion components exhibit a much stronger sparsity in the HRRPs domain compared to the rigid body. In the image domain, the rate of non-zero pixels for the rigid body and the micro-motion components is 0.01 and 0.06, respectively, verifying the sparsity of the rigid body in the image domain.
Additionally, rigid body signals exhibit strong similarity in adjacent pulses, with this similarity weakening as the distance increases. In contrast, micro-motion components do not possess this local similarity.
Figure 4 presents the autocorrelation of the rigid body and micro-motion components in the HRRP domain.
Figure 4a,c show the correlation coefficients for the rigid body and micro-motion components, respectively. To present this more clearly, we extracted the correlation coefficients of the 256th pulse and illustrated them in
Figure 4b,d. It can be observed that the HRRP matrix of the rigid body exhibits strong similarity in a few adjacent pulses, with the similarity diminishing as the pulses get farther apart. The HRRP matrix of the micro-motion components does not show this property, further validating the local similarity of the rigid body in the HRRP domain. In the next section, we will utilize these properties to construct our framework about the multiple targets separation with removal of the micro-motion connection.
5. Simulated Results and Analysis
In this section, we will validate the effectiveness of our algorithm through simulated experiments. First, we will introduce the simulation model of the radar and the targets. The radar is assumed to transmit LFM signals with a central frequency of
and a bandwidth of
, respectively. The pulse repetition frequency and pulse width are
and
, as shown in
Table 1, respectively. For the targets, the scatterer model and antenna model are illustrated in
Figure 6. The red scatterers represent the rigid body of the target, while the blue scatterers represent the micro-motion components of the target. Due to the generally small size of the antenna, the rotational radius of the micro-motion components is set as 2 m.
The motion parameters of the rigid body and the micro-motion components are presented in
Table 2. In our assumption, two targets have range overlap along the line of sight (LOS), further causing overlap in the HRRPs domain. Additionally, we add 0 dB complex Gaussian noise to the range profiles without loss of generality.
Figure 7 presents the range profiles with micro-Doppler interference and the multi-target imaging results. It can be seen that under such conditions, the range profiles of different targets severely overlap. Additionally, different targets are connected by interference stripes caused by the micro-Doppler effect in the image domain. As mentioned in
Section 4, we first apply a fixed-threshold filter in the image domain to obtain imaging results of the micro-motion components. Then, these imaging results are transformed back into the HRRPs domain to allow the Radon transform to estimate the motion parameters of targets.
Figure 8 shows the micro-motion components and their Radon transform results.
Figure 8a displays the filtered micro-motion components in the HRRPs domain. Since the translational motion of the micro-motion components is consistent with that of the rigid body, the motion parameters of the rigid body are obtained by estimating the motion parameters of the micro-motion components. The Radon transform is naturally utilized to estimate the motion of different targets.
Figure 8b shows the result of applying the Radon transform to
Figure 8a. The two peaks circled in red indicate the motion parameters of these two targets. Based on these motion parameters, we generate range migration transformation and input them into the algorithm, as shown in Algorithm 1.
The results of the preliminary separation are presented in
Figure 9.
Figure 9a,b show the range profiles of different targets, while
Figure 9c,d present imaging results of these targets. For clarity, we have magnified several range cells in
Figure 9e,f. It can be observed that the range profiles of different targets still exhibit slight range migration, resulting in defocusing in the imaging results. Additionally, signals from the micro-motion components remain in the imaging results. Therefore, further range alignment and target separation are required.
The MCRA algorithm is used to accurately estimate the range migration of the targets to generate range migration transformation.
Figure 10 shows the range profiles of different targets after range alignment and their enlarged views.
Figure 10a,d display the range profiles and the enlarged view of target1, while
Figure 10b,d display that of target2. Based on these range migration corrections, we completely separate different targets and the micro-Doppler by solving the triple-constraint cost function. The range profiles of different targets and their imaging results are shown in
Figure 11.
Figure 11a,b display the range profiles and the imaging result of target1, while
Figure 11a,b display those of target2. It is clearly visible that different targets are thoroughly separated, and the micro-Doppler interference is suppressed. The image contrast and the image entropy of different stages are given in
Table 3. It is seen that our proposed method achieves the entire separation step-by-step.
Additionally, we obtain the rotation period of the micro-motion components by performing an auto-correlation transformation, further validating the effectiveness of the separation.
Figure 12 shows the aligned micro-motion components in the HRRPs domain and their auto-correlated results.
Figure 12a,b indicate the aligned range profiles and the auto-correlated result of target1, while
Figure 12c,d show those of target2. It is seen that the differences of the peaks in
Figure 12b,d are about 198 and 159. The pulse repetition interval of the radar is
. Therefore, the rotational velocities of these two components are
and
. The theoretical rotational velocities are about
and
. The estimated periods are close to the theoretical periods, confirming the effectiveness of our separation procedure.
To thoroughly evaluate the effectiveness of the proposed framework, we selected three comparison methods: a TF-based method, a segmentation-based method, and a parameter-based method. Moreover, to further discuss the effectiveness of our algorithm, we analyze scenarios where both the rigid body and the micro-motion components of different targets overlap in the image domain.
Figure 13 displays the separation results using these methods, alongside the outcomes of our proposed method.
Figure 13a shows the separation result using the TF-based method, which is notably affected by cross-term interference and target overlap, leading to reduced separation performance.
Figure 13b presents the separation result from the segmentation-based method, where performance is significantly impacted by the broadening and overlapping of targets in the Doppler domain.
Figure 13c demonstrates the parameter-based method’s separation results, revealing that interference from different targets lowers the parameter estimation accuracy and thus hampers separation effectiveness.
Figure 13d illustrates the separation result using the proposed method, indicating the effectiveness of the proposed method. Furthermore,
Table 4 summarizes the separation performance of the different methods. Our proposed method achieves the lowest image entropy and highest image contrast, indicating its superior ability in target separation.
To evaluate the performance of our algorithm under varying signal-to-noise ratios (SNRs),
Figure 14 illustrates the separation results at −5 dB, 0 dB, 5 dB, and 10 dB SNR levels. The algorithm achieves satisfactory results within the 0–10 dB range. However, its effectiveness significantly diminishes when the SNR drops below 0 dB.
To demonstrate the impact of the initial value of
on the separation performance,
Figure 15 shows the curve related to the separation contrast and the initial value of
. The initial value of the
matrix varies from 0 to 10. It is seen that the initial value of
has a minimal influence on the accuracy of the algorithm, demonstrating the robustness of our framework. Additionally, The variation in the value of the cost function is presented in
Figure 16, which demonstrates that our algorithm exhibits reasonable convergence.