Next Article in Journal
Selecting Erosion- and Deposition-Dominated Zones in the Jezero Delta Using a Water Flow Model for Targeting Future In Situ Mars Surface Missions
Previous Article in Journal
Assessment of Habitat Quality in Arid Regions Incorporating Remote Sensing Data and Field Experiments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multiple Targets ISAR Imaging Method with Removal of Micro-Motion Connection Based on Joint Constraints

by
Hongxu Li
1,
Qinglang Guo
2,
Zihan Xu
1,
Xinfei Jin
1,
Fulin Su
1,* and
Xiaodi Li
1
1
School of Electronics and Information Engineering, Harbin Institute of Technology, Harbin 150006, China
2
CETC Academy of Electronics and Information Technology Group Co., Ltd., China Academic of Electronics and Information Technology, Beijing 100070, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(19), 3647; https://doi.org/10.3390/rs16193647
Submission received: 30 July 2024 / Revised: 25 September 2024 / Accepted: 26 September 2024 / Published: 29 September 2024

Abstract

:
Combining multiple data sources, Digital Earth is an integrated observation platform based on air–space–ground–sea monitoring systems. Among these data sources, the Inverse Synthetic Aperture Radar (ISAR) is a crucial observation method. ISAR is typically utilized to monitor both military and civilian ships due to its all-day and all-weather superiority. However, in complex scenarios, multiple targets may exist within the same radar antenna beam, resulting in severe defocusing due to different motion conditions. Therefore, this paper proposes a multiple-target ISAR imaging method with the removal of micro-motion connections based on the integration of joint constraints. The fully motion-compensated targets exhibit low rank and local similarity in the high-resolution range profile (HRRP) domain, while the micro-motion components possess sparsity. Additionally, targets display sparsity in the image domain. Inspired by this, we formulate a novel optimization by promoting the low-rank, the Laplacian, and the sparsity constraints of targets and the sparsity constraints of the micro-motion components. This optimization problem is solved by the linearized alternative direction method with adaptive penalty (LADMAP). Furthermore, the different motions of various targets degrade their inherent characteristics. Therefore, we integrate motion compensation transformation into the optimization, accordingly achieving the separation of rigid bodies and the micro-motion components of different targets. Experiments based on simulated data demonstrate the effectiveness of the proposed method.

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 X O Y 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
S 1 t r , t m = i = 1 Z I i rect t r 2 R i t m / c T P exp j 2 π f c t r 2 R i t m c exp j π γ t r 2 R i t m c 2
where t r and t m represent the fast time and the slow time, respectively. Z denotes the total number of scatterers. I i represents the amplitude of the i-th scatterer. f c , γ , and T p represent the basic frequency, the chirp rate, and the pulse width of the signal, respectively. R i 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
S 2 t r , t m = i = 1 Z I i rect t r 2 R i t m / c T P exp j 4 π γ c t r 2 R ref c Δ R i t m exp j 4 π f c c Δ R i t m exp j 4 π γ c 2 Δ R i 2 t m
where R r e f denotes the reference range of the dechirp signal. Δ R i is the difference between the instantaneous distance and the reference range of the i-th scatterer. Δ R i is written as R i in the following formula as different scatterers share the same reference range. Applying the fast Fourier transform along the t r direction to the (2), the HRRPs are written as
S 3 f r , t m = i = 1 Z I i sinc T P f r + 2 γ c R i t m · exp j 4 π f c c R i t m
where f r 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 R i and its translational and rotational components are written as
R i ( t m ) = R s ( t m ) + R l ( t m )
R s t m = x i sin θ t m + y i cos θ t m
R l t m = R 0 + v t m + 1 2 a t m 2 + + 1 n ! a n t m n
where R s and R l represent the rotational motion and the translational motion, respectively. x i and y i denote the abscissa and the ordinate of the i-th scatterer in the target coordinate. θ m is the rotational angle of the target. R 0 represents the initial distance between the target center and the radar. v, a, and a n 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
R s t m = x i θ t m + y i
R l t m = R 0 + v t m + 1 2 a t m 2
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, R l 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 X O Y coordinate system, which serves as the motion coordinate system for the micro-motion components themselves. Additionally, the origin of the X O Y coordinate system is set on the rigid body scatterer, meaning that this coordinate system rotates along with the rigid body.
Suppose that the scatterer K ( x k , y k ) on the rotating antenna rotates around its own rotational center with a rotational angle of θ and a rotational shaft of r k . The initial angle of the micro-motion components is θ k . The instantaneous distance between the scatterer and the radar is given as
R K t m = x K sin θ t m + y K cos θ t m + r K cos θ t m + θ k
According to the aforementioned approximation, we obtain
R K t m = x K θ t m + y K + r K cos θ t m + θ k
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 f r direction yields
s ( t r , t m ) = I exp ( j 4 π ( γ t r + f c ) c ( x θ ( t m ) + y ) )
which can be discretized as
s n , p = I exp j 4 π γ Δ t r c y n · exp j 4 π f c Δ θ c x p
where n and p are the index of the range cell and the pulse, respectively. Suppose that the size of the HRRP matrix is N × P , (12) can be written as
s n , p = I a T b
where a = [ a 1 , , a N ] and b = [ b 1 , , b P ] with a i = exp ( j 4 π γ y i Δ t r / c ) and b i = exp ( j 4 π f c x i Δ t r / c ) , respectively. For a target composed of U scattering points, its echo can be expressed as
s n , p = u U I u a u T b u
Based on the properties of rank, the maximum rank of this HRRP matrix does not exceed U . 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.

4. Multiple Targets Separation Method with Removal of Micro-Motion Connection Based on Dual Optimization

In this section, we will explain our proposed multiple targets method with removal of micro-motion connection based on dual optimization in detail. First, we roughly filter out the micro-motion components in the image domain and transform them back to the HRRPs domain. Then, we estimate the motion parameters of these filtered micro-motion components using the Radon transform. These parameters are then utilized to construct a coarse motion compensation function. Following that, we integrate the coarse range migration transformation with the dual-constrained optimization model to form a cost function. This cost function is solved by utilizing LADMAP to coarsely separate targets and micro-motion components. Subsequently, we employ a correlation-based range alignment algorithm to measure the motion parameters of the coarsely separated targets. These parameters are utilized to construct a more precious range migration transformation. This more accurate range migration transformation is then incorporated into a triple-constrained optimization model. Similarly, this cost function is solved using LADMAP to ensure convergence. Finally, different targets and micro-motion components are completely separated. Figure 5 illustrates the detailed procedure of our proposed algorithm. We will elaborate on each step of this process in the following subsections.

4.1. Constructing Coarse Range Migration Transformation

First, we transform the multi-target echo with micro-Doppler effects into the image domain and use a constant threshold to roughly filter out the micro-motion components. Then, the micro-motion components are transformed back into the HRRPs domain. Since the micro-motion components and the targets share the same motion, the Radon transform is allowed to roughly extract the motion parameters of the targets. Based on these motion parameters, a coarse range migration transformation is constructed.

4.2. Coarse Separating through Double-Constraint Optimization

Motion-compensated targets exhibit low-rank characteristics in the HRRPs domain and sparsity in the image domain. Micro-Doppler components also exhibit sparsity as they occupy fewer range cells. Since the rotational speed of the micro-motion components is generally much higher than that of the rigid body, the low-rank characteristics of the micro-motion components in the HRRPs domain are weaker than those of the rigid body. Based on these inherent properties, we construct a dual-constrained cost function as follows:
d = 1 D rank L d + d = 1 D S d 0 + S M D 0 s . t . H = d = 1 D R d L d + S M D L d = IFFT S d
where L d and S d denote the rigid body in the HRRPs domain and in the image domain with the target index of d. S M D represents the micro-motion components in the HRRPs domain. H denotes the original multi-target echo with micro-Doppler effects. R d is the range migration transformation of the d-th target. Since (15) is a non-convex problem, rank L d and S d 0 are generally relaxed to L d and S d 1 . This leads to
d = 1 D L d + d = 1 D S d 1 + S M D 1 s . t . H = d = 1 D R d L d + S M D L d = IFFT S d
whose augmented Lagrange function is accordingly written as
L L d , S d , S M D , Y d , Y E = d = 1 D L d + d = 1 D S d 1 + S M D 1 + d = 1 D Y d , L d IFFT S d + d = 1 D μ d 2 L d IFFT S d F 2 + Y E , H R 1 L 1 R 2 L 2 S M D + μ E 2 H R 1 L 1 R 2 L 2 S M D F 2
where · , · 1 , and · F denote the nuclear norm, l 1 -norm, and Frobenius norm, respectively. < · > denotes the inner product operator. Y d and Y E represent the Lagrange multiplier matrices of the d-th target and the overall summation equation, respectively. Suppose that two targets emerge in the same radar antenna beam (i.e., D = 2 ), the augmented Lagrange function is decomposed by the LADMAP into the following sub-problems as
L d k + 1 = min L d L L 1 , L 2 , Y k S d k + 1 = min S d L S 1 , S 2 , Y k S M D k + 1 = min S M D L L 1 , L 2 , Y k Y 1 k + 1 = Y 1 k + μ 1 L 1 k IFFT S 1 Y 2 k + 1 = Y 2 k + μ 2 L 2 k IFFT S 2 Y E k + 1 = Y E k + μ E H R 1 L 1 R 2 L 2 S M D
where k denotes the iteration index. Through alternately solving the sub-problems in (18), we make (17) solvable. Furthermore, the LADMAP algorithm decomposes and linearizes the complex optimization problem into subproblems, ensuring the convergence of the solution process. To update L 1 , we fix all variables except L 1 , and subsequently derive the augmented Lagrange function as
L L 1 = L 1 + Y 1 , L 1 IFFT S 1 + μ 1 2 L 1 IFFT S 1 F 2 + Y H , H R 1 L 1 R 2 L 2 S M D + μ E 2 H R 1 L 1 R 2 L 2 S M D F 2
Based on the LAMDAP, we approximate (19) at L 1 = L 1 ( k ) as
L L 1 L 1 + L 1 q L 1 k , L 1 L 1 k + η 2 L 1 L 1 k F 2 = L 1 + η 2 L 1 L 1 k + L 1 q L 1 k / η F 2
where L 1 q L 1 k is written as
q L 1 = Y 1 + μ 1 L 1 IFFT S 1 R 1 Y E μ E R 1 H R 1 L 1 R 2 L 2 S M D
Utilizing the singular value shrinkage threshold (SVT) operator, the function has a closed-form solution as
L 1 k + 1 = 1 / η L 1 k 1 / η q L 1
where y x denote shrinking the singular matrix of x utilizing the soft threshold y [31]. This operator can be specifically written as
y x = U ζ y σ V H
where σ denote the singular matrix, i.e., x = U σ V H . U and V denote the unitary matrix about x . ζ y σ denotes the soft threshold operator, which can be written as
ζ y σ = sgn σ max σ y , 0
Since the process of updating L 1 is the same as L 2 , we omit it here and derive the augmented Lagrange function of S 1 as
L S 1 = S 1 1 + Y 1 , L 1 IFFT S 1 + μ 1 2 L 1 IFFT S 1 F 2
Similarly, the gradient of a part of (25) is written as
q S 1 = FFT Y 1 μ 1 FFT L 1 IFFT S 1
which leads (25) to be approximated as
L S 1 S 1 + S 1 q S 1 k , S 1 S 1 k + η 2 S 1 S 1 k F 2 = S 1 + η 2 S 1 S 1 k + S 1 q S 1 k / η F 2
Utilizing (24), the close-form solution of (27) is given as
S 1 k + 1 = 1 / η S 1 k 1 / η q S 1
As previously mentioned, we omit the processing of S 1 .
As for the range profiles of the micro-motion components S M D , the augmented Lagrange function is written as
L S M D = S M D 1 + Y E , H R 1 L 1 R 2 L 2 S M D + μ E 2 H R 1 L 1 R 2 L 2 S M D F 2
The gradient is also derived as
q S M D = Y E μ E H R 1 L 1 R 2 L 2 S M D
The close-form solution is derived as
S M D k + 1 = 1 / η S M D k 1 / η q S M D
Finally, the Lagrange multipliers Y 1 , Y 2 , and Y E are updated as
Y 1 k + 1 = Y 1 k + μ 1 L 1 IFFT S 1
Y 2 k + 1 = Y 2 k + μ 2 L 2 IFFT S 2
Y E k + 1 = Y E k + μ E H R 1 L 1 R 2 L 2 S M D
To illustrate more clearly, such a separation model is summarized in Algorithm 1. ε is set as 1.05 to accelerate the convergence. Additionally, the employed LADMAP algorithm iteratively refines the Lagrange multipliers and adaptively adjusts during computation. Therefore, we initialize the Lagrange multiplier matrix Y 1 , Y 2 , and Y E as the zero matrix in Algorithm 1.
Algorithm 1 Double-constraint-based coarse separation
Input: 
Range profiles H , Coarse motion compensation transformation of target1 R 1 · , Coarse motion compensation transformation of target2 R 2 · .
Output: 
HRRPs of target1 L 1 , HRRPs of target2 L 2 , HRRPs of micro-motion components S M D , Imaging result of target1 S 1 , Imaging result of target2 S 2 .
Initialization: 
L 1 0 = L 2 0 = S 1 ( 0 ) = S 2 ( 0 ) = S 1 ( 0 ) = Y E ( 0 ) = Y 1 ( 0 ) = Y 2 ( 0 ) = 0 N × M , μ 1 = μ 2 = μ E = 10 4 , ε = 1.05 .
 1:
while not converged do
 2:
    The update of the rigid body in the HRRPs domain.
 3:
     q L 1 = Y 1 + μ 1 L 1 IFFT S 1 R 1 Y E
             μ E R 1 H R 1 L 1 R 2 L 2 S M D
 4:
     L 1 k + 1 = 1 / η L 1 k 1 / η q L 1
 5:
     q L 2 = Y 2 + μ 2 L 2 IFFT S 2 R 2 Y E
             μ E R 2 H R 1 L 1 R 2 L 2 S M D
 6:
     L 2 k + 1 = 1 / η L 2 k 1 / η q L 2
 7:
    The update of the rigid body in the image domain.
 8:
     q S 1 = FFT Y 1 μ 1 FFT L 1 IFFT S 1
 9:
     S 1 k + 1 = 1 / η S 1 k 1 / η q S 1
10:
     q S 2 = FFT Y 2 μ 2 FFT L 2 IFFT S 2
11:
     S 2 k + 1 = 1 / η S 2 k 1 / η q S 2
12:
    The update of the micro-motion components in the HRRPs domain.
13:
     q S M D = Y E μ E H R 1 L 1 R 2 L 2 S M D
14:
     S M D k + 1 = 1 / η S M D k 1 / η q S M D
15:
    The update of the the Lagrange multiplier matrices.
16:
     Y 1 k + 1 = Y 1 k + μ 1 L 1 IFFT S 1
17:
     Y 2 k + 1 = Y 2 k + μ 1 L 2 IFFT S 2
18:
     Y E k + 1 = Y E k + μ E H R 1 L 1 R 2 L 2 S M D
19:
     μ 1 = ε μ 1 , μ 2 = ε μ 2 , μ E = ε μ E
20:
     k = k + 1 ;
21:
end while

4.3. Constructing Accurate Range Migration Transformation

By solving the coarse cost function, different targets and micro-motion components are preliminarily separated. Consequently, we employ a correlation-based range alignment method to estimate the range migration of different targets. Specifically, we use the maximum-correlation range alignment (MCRA) algorithm to estimate the precise range migration of different targets. These parameters are then used to generate more accurate r range migration transformation.

4.4. Accurate Separating through Triple-Constraint Optimization

As previously mentioned, micro-motion components exhibit high rotational speeds, significantly reducing their similarity between adjacent pulses. Typically, rigid bodies have lower rotational speeds, with a rotation angle of 2 degrees within the imaging interval. This characteristic motivates us to incorporate local similarity into the multi-target optimization. Specifically, we integrate Laplacian regularization into the cost function, forming a triple-constrained cost function as follows:  
d = 1 D rank L d + β 1 Tr L d BL d H + S d 0 + S M D 0 s . t . H = d = 1 D R d L d + S M D L d = IFFT S d
where β denotes the penalty parameter of the Laplacian regularization. B denotes the graph Laplacian matrix [32], which is given as
B = D W
where D and W represent the degree matrix and the weight matrix, respectively. These two matrices are further given as
D i i = j W i j
W i j = 1 i = j a i j < b 0 otherwise
where a and b are set as 0.8 and 3. Additionally, the cost function is similarly relaxed as
d = 1 D rank L d + β 1 Tr L d BL d H + S d 0 + S M D 1 s . t . H = d = 1 D R d L d + S M D L d = IFFT S d
Utilizing the LADMAP, we derive the augmented Lagrange function as
L L d , S d , S M D , Y d , Y H = d = 1 D L d + β 1 Tr L d BL d H + S d 1 + S M D 1 + d = 1 D Y d , L d IFFT S d + μ d 2 L d IFFT S d F 2 + Y E , H R 1 L 1 R 2 L 2 S M D + μ E 2 H R 1 L 1 R 2 L 2 S M D F 2
Also, the cost function is decomposed as
L d k + 1 = min L d L L 1 , L 2 , Y k S d k + 1 = min S d L S 1 , S 2 , Y k S M D k + 1 = min S M D L L 1 , L 2 , Y k Y 1 k + 1 = Y 1 k + μ 1 L 1 k IFFT S 1 Y 2 k + 1 = Y 2 k + μ 2 L 2 k IFFT S 2 Y E k + 1 = Y E k + μ E H R 1 L 1 R 2 L 2 S M D
Fixing all the variables except L 1 , we derive the augmented Lagrange function of L 1 as
L L 1 = L 1 + β 1 Tr L 1 BL 1 H + Y 1 , L 1 IFFT S 1 + μ 1 2 L 1 IFFT S 1 F 2 + Y E , H R 1 L 1 R 2 L 2 S M D + μ E 2 H R 1 L 1 R 2 L 2 S M D F 2
The gradient of a part of (42) at L 1 = L 1 ( k ) is written as
q L 1 = β L 1 k B + L 1 k B T + Y 1 + μ 1 L 1 IFFT S 1 R 1 Y 3 μ E R 1 H R 1 L 1 R 2 L 2 S M D
Therefore, L 1 is reasonably approximated as
L 1 k + 1 = 1 / η L 1 k 1 / η q L 1
Also, the strategy of updating L 2 is similar to L 1 . As for updating S 1 , S 2 , S M D , Y 1 , Y 2 , and Y E , the rule is the same as Section 4.2. In summary, we give the overall steps in Algorithm 2. It is noted that we perform 12 multiplications, 22 FFT, and 2 SVD for Algorithm 1. As for Algorithm 2, we perform 16 multiplications, 22 FFT, and 2 SVD in each iteration. Suppose that the size of the HRRP matrix is N × N , the computational complexity in a single iteration is about O N 3 + N 2 log N . Furthermore, the algorithm typically converges within 70 iterations, which ensures that it does not impose a significant computational burden.
Algorithm 2 Triple-constraint-based accurate separation
Input: 
Range profiles H , Coarse motion compensation transformation of target1 R 1 · , Coarse motion compensation transformation of target2 R 2 · .
Output: 
HRRPs of target1 L 1 , HRRPs of target2 L 2 , HRRPs of micro-motion components S M D , Imaging result of target1 S 1 , Imaging result of target2 S 2 .
Initialization: 
L 1 0 = L 2 0 = S 1 ( 0 ) = S 2 ( 0 ) = S 1 ( 0 ) = Y E ( 0 ) = Y 1 ( 0 ) = Y 2 ( 0 ) = 0 N × M , μ 1 = μ 2 = μ E = 10 4 , β = 10 8 , ε = 1.05 .
 1:
while not converged do
 2:
    The update of the rigid body in the HRRPs domain.
 3:
     q L 1 = β L 1 k B + L 1 k B T + Y 1 + μ 1 L 1 IFFT S 1
             R 1 Y E μ E R 1 H R 1 L 1 R 2 L 2 S M D
 4:
     L 1 k + 1 = 1 / η L 1 k 1 / η q L 1
 5:
     q L 2 = β L 2 k B + L 2 k B T + Y 2 + μ 2 L 2 IFFT S 2
             R 2 Y E μ E R 2 H R 1 L 1 R 2 L 2 S M D
 6:
     L 2 k + 1 = 1 / η L 2 k 1 / η q L 2
 7:
    The update of the rigid body in the image domain.
 8:
     q S 1 = FFT Y 1 μ 1 FFT L 1 IFFT S 1
 9:
     S 1 k + 1 = 1 / η S 1 k 1 / η q S 1
10:
     q S 2 = FFT Y 2 μ 2 FFT L 2 IFFT S 2
11:
     S 2 k + 1 = 1 / η S 2 k 1 / η q S 2
12:
    The update of the micro-motion components in the HRRPs domain.
13:
     q S M D = Y E μ E H R 1 L 1 R 2 L 2 S M D
14:
     S M D k + 1 = 1 / η S M D k 1 / η q S M D
15:
    The update of the the Lagrange multiplier matrices.
16:
     Y 1 k + 1 = Y 1 k + μ 1 L 1 IFFT S 1
17:
     Y 2 k + 1 = Y 2 k + μ 1 L 2 IFFT S 2
18:
     Y E k + 1 = Y E k + μ E H R 1 L 1 R 2 L 2 S M D
19:
     μ 1 = ε μ 1 , μ 2 = ε μ 2 , μ E = ε μ E
20:
     k = k + 1 ;
21:
end while

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 f 0 = 10 GHz and a bandwidth of B = 200 MHz , respectively. The pulse repetition frequency and pulse width are PRF = 800 Hz and T 0 = 100 μ s , 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 0.001 s . Therefore, the rotational velocities of these two components are 15.87 rad / s and 19.75 rad / s . The theoretical rotational velocities are about 15.71 rad / s and 18.84 rad / s . 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 Y d on the separation performance, Figure 15 shows the curve related to the separation contrast and the initial value of Y d . The initial value of the Y d matrix varies from 0 to 10. It is seen that the initial value of Y d 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.

6. Conclusions

This paper proposes a multi-target ISAR imaging method based on joint constraints with removal of the micro-motion connection. We first estimate the rough translational motion through filtered micro-motion components, reasonably transforming the micro-Doppler interference stripes into the motion information. Following that, we construct the range migration transformation and input this transformation into the double-constraint optimization, thereby roughly separating the targets and the micro-motion components. The separated targets are then utilized to measure the accurate motions and to form the accurate range migration transformation through the MCRA method. Subsequently, the triple-constraints optimization is employed to entirely separate the targets and the micro-motion components. The effectiveness of the proposed method is validated by the simulated experiments.

Author Contributions

Conceptualization, H.L. and F.S.; methodology, H.L.; software, H.L.; validation, H.L.; investigation, H.L.; resources, F.S.; data curation, H.L.; writing—original draft preparation, H.L. and X.L.; writing—review and editing, H.L.; visualization, H.L., X.J., Z.X. and Q.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank the anonymous reviewers and editors for their valuable comments, which were crucial in improving the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Anger, S.; Jirousek, M.; Dill, S.; Peichl, M. Inverse Synthetic Aperture Radar Imaging of Space Targets Using Wideband Pseudo-Noise Signals with Low Peak-to-Average Power Ratio. Remote Sens. 2024, 16, 1809. [Google Scholar] [CrossRef]
  2. Chen, V.; Martorella, M. Inverse Synthetic Aperture Radar Imaging: Principles, Algorithms and Applications; IET: London, UK, 2014. [Google Scholar]
  3. Long, B.; Tang, P.; Wang, F.; Jin, Y.Q. 3-D Reconstruction of Space Target Based on Silhouettes Fused ISAR-Optical Images. IEEE Trans. Geosci. Remote Sens. 2024, 62, 5213919. [Google Scholar] [CrossRef]
  4. Gao, Y.; Xing, M.; Li, Y.; Sun, W.; Zhang, Z. Joint translational motion compensation method for ISAR imagery under low SNR condition using dynamic image sharpness metric optimization. IEEE Trans. Geosci. Remote Sens. 2021, 60, 5108515. [Google Scholar] [CrossRef]
  5. Wang, J.; Zhang, L.; Du, L.; Yang, D.; Chen, B. Noise-robust motion compensation for aerial maneuvering target ISAR imaging by parametric minimum entropy optimization. IEEE Trans. Geosci. Remote Sens. 2019, 57, 4202–4217. [Google Scholar] [CrossRef]
  6. Li, Y.; Fu, Y.; Li, X.; Li, L.W. An ISAR imaging method for multiple moving targets based on fractional Fourier transformation. In Proceedings of the 2009 IEEE Radar Conference, Pasadena, CA, USA, 4–8 May 2009; IEEE: Washington, DC, USA, 2009; pp. 1–6. [Google Scholar]
  7. Zhang, D.; Chen, W.; Wang, D. ISAR imaging of multiple moving targets based on RSPWVD-Hough transform. In Proceedings of the 2008 Asia-Pacific Microwave Conference, Hong Kong, China, 16–20 December 2008; IEEE: Washington, DC, USA, 2008; pp. 1–4. [Google Scholar]
  8. Luo, X.; He, Q.; Lv, G.; Shang, C. TSAR imaging of multiple targets based on adaptive Gaussian chirplet decomposition. In Proceedings of the 2006 CIE International Conference on Radar, Shanghai, China, 16–19 October 2006; IEEE: Washington, DC, USA, 2006; pp. 1–4. [Google Scholar]
  9. Fan, L.; Yiming, P.; Shunji, H. Multi-target imaging processing algorithms of ISAR based on time-frequency analysis. In Proceedings of the 2006 CIE International Conference on Radar, Shanghai, China, 16–19 October 2006; IEEE: Washington, DC, USA, 2006; pp. 1–4. [Google Scholar]
  10. Xiao, D.; Su, F.; Wu, J. Multi-target ISAR imaging based on image segmentation and short-time Fourier transform. In Proceedings of the 2012 5th International Congress on Image and Signal Processing, Chongqing, China, 16–18 October 2012; IEEE: Washington, DC, USA, 2012; pp. 1832–1836. [Google Scholar]
  11. Wang, A.; Mao, Y.; Chen, Z. Imaging of multitargets with ISAR based on the time-frequency distribution. In Proceedings of the Proceedings of ICASSP’94. IEEE International Conference on Acoustics, Speech and Signal Processing, Adelaide, SA, Australia, 19–22 April 1994; IEEE: Washington, DC, USA, 1994; p. V-173. [Google Scholar]
  12. Zhang, Y.; Xu, N.; Li, N.; Guo, Z. A Multi-Domain Joint Novel Method for ISAR Imaging of Multi-Ship Targets. Remote Sens. 2023, 15, 4878. [Google Scholar] [CrossRef]
  13. Chen, J.; Xiao, H.; Fan, H.; Song, Z. ISAR imaging of multiple moving targets using signals separation. In Proceedings of the Proceedings 2013 International Conference on Mechatronic Sciences, Electric Engineering and Computer (MEC), Shenyang, China, 20–22 December 2013; IEEE: Washington, DC, USA, 2013; pp. 1156–1159. [Google Scholar]
  14. Zhao, J.; Zhang, Y.Q.; Wang, X.; Wang, S.; Shang, F. A novel method for ISAR imaging of multiple maneuvering targets. Prog. Electromagn. Res. M 2019, 81, 43–54. [Google Scholar] [CrossRef]
  15. Bai, X.; Zhou, F.; Xing, M.; Bao, Z. A novel method for imaging of group targets moving in a formation. IEEE Trans. Geosci. Remote Sens. 2011, 50, 221–231. [Google Scholar] [CrossRef]
  16. Kong, L.; Zhang, W.; Zhang, S.; Zhou, B. Radon transform and the modified envelope correlation method for ISAR imaging of multi-target. In Proceedings of the 2010 IEEE Radar Conference, Arlington, VA, USA, 10–14 May 2010; IEEE: Washington, DC, USA, 2010; pp. 637–641. [Google Scholar]
  17. Xiao, D.; Su, F.; Wu, J. A method of ISAR imaging for multiple targets. In Proceedings of the 2012 IEEE 11th International Conference on Signal Processing, Beijing, China, 21–25 October 2012; IEEE: Washington, DC, USA, 2012; Volume 3, pp. 2011–2015. [Google Scholar]
  18. Jun, S.; Xiaoling, Z.; Shuwei, H. Multi-target ISAR imaging method. In Proceedings of the 2005 IEEE International Geoscience and Remote Sensing Symposium, Seoul, Republic of Korea, 29–29 July 2005; IGARSS’05. IEEE: Washington, DC, USA, 2005; Volume 7, pp. 4745–4748. [Google Scholar]
  19. Chen, J.; Xiao, H.; Song, Z.; Fan, H. Simultaneous ISAR imaging of group targets flying in formation. Chin. J. Aeronaut. 2014, 27, 1554–1561. [Google Scholar] [CrossRef]
  20. Liu, L.; Zhou, F.; Tao, M.; Zhang, Z. A novel method for multi-targets ISAR imaging based on particle swarm optimization and modified CLEAN technique. IEEE Sens. J. 2015, 16, 97–108. [Google Scholar] [CrossRef]
  21. Liu, L.; Zhou, F.; Guo, Y.; Tao, M.; Sun, P.; Zhang, Z. A novel method for multi-targets ISAR imaging. In Proceedings of the 2015 IEEE 5th Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Singapore, 1–4 September 2015; IEEE: Singapore, 2015; pp. 314–319. [Google Scholar]
  22. Choi, G.; Park, S.; Kim, H.; Kim, K. ISAR imaging of multiple targets based on particle swarm optimization and hough transform. J. Electromagn. Waves Appl. 2009, 23, 1825–1834. [Google Scholar] [CrossRef]
  23. Zhao, J.; Zhang, M.; Wang, X.; Nie, D. Parameters estimation and ISAR imaging of multiple maneuvering targets based on an order reduction method for cubic chirps. J. Electromagn. Waves Appl. 2017, 31, 1658–1675. [Google Scholar] [CrossRef]
  24. Chen, J.; Xiao, H.; Song, Z.; Fan, H. Imaging targets moving in formation using parametric compensation. EURASIP J. Adv. Signal Process. 2014, 2014, 1–9. [Google Scholar] [CrossRef]
  25. Li, J.; Ling, H. Application of adaptive chirplet representation for ISAR feature extraction from targets with rotating parts. IEE Proc.-Radar Sonar Navig. 2003, 150, 284–291. [Google Scholar] [CrossRef]
  26. Yuan, B.; Chen, Z.; Xu, S. Micro-Doppler analysis and separation based on complex local mean decomposition for aircraft with fast-rotating parts in ISAR imaging. IEEE Trans. Geosci. Remote Sens. 2013, 52, 1285–1298. [Google Scholar] [CrossRef]
  27. Lu, L.; Chen, P.; Wu, L. A RPCA-based ISAR imaging method for micromotion targets. Sensors 2020, 20, 2989. [Google Scholar] [CrossRef] [PubMed]
  28. Zhou, W.; Yeh, C.; Jin, R.; Li, Z.; Song, S.; Yang, J. ISAR imaging of targets with rotating parts based on robust principal component analysis. IET Radar Sonar Navig. 2017, 11, 563–569. [Google Scholar] [CrossRef]
  29. Zhang, S.; Liu, Y.; Li, X.; Hu, D. Removal of micro-Doppler effect of ISAR image based on Laplacian regularized nonconvex low-rank representation. IEEE Trans. Image Process. 2021, 30, 6446–6458. [Google Scholar] [CrossRef] [PubMed]
  30. Lin, Z.; Liu, R.; Su, Z. Linearized alternating direction method with adaptive penalty for low-rank representation. Adv. Neural Inf. Process. Syst. 2011, 24. [Google Scholar]
  31. Cai, J.F.; Candès, E.J.; Shen, Z. A singular value thresholding algorithm for matrix completion. SIAM J. Optim. 2010, 20, 1956–1982. [Google Scholar] [CrossRef]
  32. Yin, M.; Gao, J.; Lin, Z. Laplacian regularized low-rank representation and its applications. IEEE Trans. Pattern Anal. Mach. Intell. 2015, 38, 504–517. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The imaging geometry of multiple targets.
Figure 1. The imaging geometry of multiple targets.
Remotesensing 16 03647 g001
Figure 2. The imaging geometry of the micro-motion components.
Figure 2. The imaging geometry of the micro-motion components.
Remotesensing 16 03647 g002
Figure 3. Range profiles and the imaging results of the ridge body and the micro-motion components. (a) Range profiles of the ridge body. (b) Imaging result of the ridge body. (c) Range profiles of the micro-motion components. (d) Imaging result of the micro-motion components.
Figure 3. Range profiles and the imaging results of the ridge body and the micro-motion components. (a) Range profiles of the ridge body. (b) Imaging result of the ridge body. (c) Range profiles of the micro-motion components. (d) Imaging result of the micro-motion components.
Remotesensing 16 03647 g003
Figure 4. Range profiles and the imaging results of the ridge body and the micro-motion components. (a) Correlation coefficients of the ridge body. (b) Correlation coefficients of 256th pulse of the ridge body. (c) Correlation coefficients of the micro-motion components. (d) Correlation coefficients of 256th pulse of the micro-motion components.
Figure 4. Range profiles and the imaging results of the ridge body and the micro-motion components. (a) Correlation coefficients of the ridge body. (b) Correlation coefficients of 256th pulse of the ridge body. (c) Correlation coefficients of the micro-motion components. (d) Correlation coefficients of 256th pulse of the micro-motion components.
Remotesensing 16 03647 g004
Figure 5. The procedure of the proposed multi-target separation method.
Figure 5. The procedure of the proposed multi-target separation method.
Remotesensing 16 03647 g005
Figure 6. The scatterer model and antenna model of targets.
Figure 6. The scatterer model and antenna model of targets.
Remotesensing 16 03647 g006
Figure 7. Range profiles and their imaging results. (a) Range profiles of targets. (b) Imaging result of targets.
Figure 7. Range profiles and their imaging results. (a) Range profiles of targets. (b) Imaging result of targets.
Remotesensing 16 03647 g007
Figure 8. The micro-motion components and their Radon transform results. (a) Filtered micro-motion components. (b) The Radon transform of the filtered components.
Figure 8. The micro-motion components and their Radon transform results. (a) Filtered micro-motion components. (b) The Radon transform of the filtered components.
Remotesensing 16 03647 g008
Figure 9. Results of the coarse separation. (a) Coarse range profiles of target1. (b) Coarse range profiles of target2. (c) Coarse imaging result of target1. (d) Coarse imaging result of target1. (e) Enlarged view of target1. (f) Enlarged view of target2.
Figure 9. Results of the coarse separation. (a) Coarse range profiles of target1. (b) Coarse range profiles of target2. (c) Coarse imaging result of target1. (d) Coarse imaging result of target1. (e) Enlarged view of target1. (f) Enlarged view of target2.
Remotesensing 16 03647 g009
Figure 10. Aligned range profiles and their enlarged views. (a) Aligned coarse range profiles of target1. (b) Enlarged view of coarse range profiles of target1. (c) Aligned coarse range profiles of target1. (d) Enlarged view of coarse range profiles of target2.
Figure 10. Aligned range profiles and their enlarged views. (a) Aligned coarse range profiles of target1. (b) Enlarged view of coarse range profiles of target1. (c) Aligned coarse range profiles of target1. (d) Enlarged view of coarse range profiles of target2.
Remotesensing 16 03647 g010
Figure 11. Results of accurate separation. (a) Accurate range profiles of target1. (b) Accurate range profiles of target2. (c) Accurate imaging result of target1. (d) Accurate imaging result of target2.
Figure 11. Results of accurate separation. (a) Accurate range profiles of target1. (b) Accurate range profiles of target2. (c) Accurate imaging result of target1. (d) Accurate imaging result of target2.
Remotesensing 16 03647 g011
Figure 12. Range profiles of micro-motion components of targets and their auto-correlation results. (a) Aligned range profiles of micro-motion components of target1. (b) The auto-correlation result of micro-motion components of target1. (c) Aligned range profiles of micro-motion components of target2. (d) The auto-correlation result of micro-motion components of target2.
Figure 12. Range profiles of micro-motion components of targets and their auto-correlation results. (a) Aligned range profiles of micro-motion components of target1. (b) The auto-correlation result of micro-motion components of target1. (c) Aligned range profiles of micro-motion components of target2. (d) The auto-correlation result of micro-motion components of target2.
Remotesensing 16 03647 g012
Figure 13. The separation results of different methods and the proposed method. (a) The separation result using the TF-based method. (b) The separation result using the segmentation-based method. (c) The separation result using the parameter-based method. (d) The separation result using the proposed method.
Figure 13. The separation results of different methods and the proposed method. (a) The separation result using the TF-based method. (b) The separation result using the segmentation-based method. (c) The separation result using the parameter-based method. (d) The separation result using the proposed method.
Remotesensing 16 03647 g013
Figure 14. The separation results under different SNR conditions. (a) The separation result under the −5 dB SNR condition. (b) The separation result under the 0 dB SNR condition. (c) The separation result under the 5 dB SNR condition. (d) The separation result using the proposed method.
Figure 14. The separation results under different SNR conditions. (a) The separation result under the −5 dB SNR condition. (b) The separation result under the 0 dB SNR condition. (c) The separation result under the 5 dB SNR condition. (d) The separation result using the proposed method.
Remotesensing 16 03647 g014
Figure 15. The curve related to the separation contrast and the initial value of Y d .
Figure 15. The curve related to the separation contrast and the initial value of Y d .
Remotesensing 16 03647 g015
Figure 16. The variation in the value of the cost function.
Figure 16. The variation in the value of the cost function.
Remotesensing 16 03647 g016
Table 1. System parameters of the simulated radar.
Table 1. System parameters of the simulated radar.
Motion ParametersValue
Central frequency f 0 10 GHz
Bandwidth B 200 MHz
Pulse repetition frequency PRF 800 Hz
Pulse width T 0 100 μ s
Table 2. Motion parameters of rigid body and micro-motion components.
Table 2. Motion parameters of rigid body and micro-motion components.
System ParametersTarget1Target2Antenna1Antenna2
X-direction position 3000 m 2900 m Same as target1Same as target2
Y-direction position 2000 m 2140 m Same as target1Same as target2
X-direction velocity 7 m / s 8 m / s Same as target1Same as target2
Y-direction velocity 0 m / s 0 m / s Same as target1Same as target2
Rotational speed 0.02 rad / s 0.02 rad / s 2 π rad / s 2.5 π rad / s
Table 3. The image contrast and the image entropy in different stages.
Table 3. The image contrast and the image entropy in different stages.
StagesTarget IndexImage ContrastImage Entropy
Original signalAll targets 12.9 4.9
Coarse seprationTarget1 15.8 3.3
Target2 23.1 9.4
Accurate seprationTarget1 43.4 3.3
Target2 41.8 9.4
Table 4. The image contrast and the image entropy of different methods.
Table 4. The image contrast and the image entropy of different methods.
Different MethodsTarget IndexImage ContrastImage Entropy
TF-based methodTarget1 21.8 8.2
Target2 21.5 10.1
Segmentation-based methodTarget1 14.8 6.4
Target2 14.3 9.7
Parameter-based methodTarget1 28.1 7.7
Target2 29.7 10.3
The proposed methodTarget1 43.4 3.3
Target2 41.8 9.4
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

Li, H.; Guo, Q.; Xu, Z.; Jin, X.; Su, F.; Li, X. A Multiple Targets ISAR Imaging Method with Removal of Micro-Motion Connection Based on Joint Constraints. Remote Sens. 2024, 16, 3647. https://doi.org/10.3390/rs16193647

AMA Style

Li H, Guo Q, Xu Z, Jin X, Su F, Li X. A Multiple Targets ISAR Imaging Method with Removal of Micro-Motion Connection Based on Joint Constraints. Remote Sensing. 2024; 16(19):3647. https://doi.org/10.3390/rs16193647

Chicago/Turabian Style

Li, Hongxu, Qinglang Guo, Zihan Xu, Xinfei Jin, Fulin Su, and Xiaodi Li. 2024. "A Multiple Targets ISAR Imaging Method with Removal of Micro-Motion Connection Based on Joint Constraints" Remote Sensing 16, no. 19: 3647. https://doi.org/10.3390/rs16193647

APA Style

Li, H., Guo, Q., Xu, Z., Jin, X., Su, F., & Li, X. (2024). A Multiple Targets ISAR Imaging Method with Removal of Micro-Motion Connection Based on Joint Constraints. Remote Sensing, 16(19), 3647. https://doi.org/10.3390/rs16193647

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