Next Article in Journal
Availability of Automatic Identification System (AIS) Based on Spectral Analysis of Mean Time to Repair (MTTR) Determined from Dynamic Data Age
Next Article in Special Issue
SEAN: A Simple and Efficient Attention Network for Aircraft Detection in SAR Images
Previous Article in Journal
Towards Single-Component and Dual-Component Radar Emitter Signal Intra-Pulse Modulation Classification Based on Convolutional Neural Network and Transformer
Previous Article in Special Issue
A Multi-Pulse Cross Ambiguity Function for the Wideband TDOA and FDOA to Locate an Emitter Passively
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Micro-Doppler Curves Extraction of Space Target Based on Modified Synchro-Reassigning Transform and Ridge Segment Linking

School of Automation, Central South University, Changsha 410083, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(15), 3691; https://doi.org/10.3390/rs14153691
Submission received: 11 June 2022 / Revised: 24 July 2022 / Accepted: 28 July 2022 / Published: 2 August 2022
(This article belongs to the Special Issue Microwave Remote Sensing for Object Detection)

Abstract

:
The micro-movement feature is recognized as one of the practical features of space target recognition in academic circles. The separation of the micro-Doppler curve of the scattering center is the key to feature extraction and parameter estimation, which depends on the time–frequency analysis method. The existing techniques have low separation accuracy and adaptability when there are overlap and noise in the time–frequency domain. This paper proposes a micro-Doppler feature extraction algorithm of a space target based on the modified synchro-reassigning transform (MSRT) and ridge segment linking. The MSRT can eliminate repeated assignment problems, has more accurate micro-Doppler frequency estimates than the synchro-reassigning transform, and has lower computational complexity than second-order synchronous compression and synchronous extraction transforms. The re-linking of the ridge realizes the correct connection of the micro-Doppler curves of each scattering center. The simulation data and the electromagnetic calculation data verify the method’s effectiveness.

1. Introduction

Micro-motion reflects the electromagnetic scattering, geometric structure, and motion characteristics of a space target [1], so the m-D feature is widely used in the recognition task of space targets. Separating the m-D curves from the scattering center of the target is the basis of parameter estimation and feature extraction.
Micro-motion features can be obtained from the radar cross-section (RCS) and the joint time–frequency (JTF) distribution, high-resolution range profiles (HRRPs), and the range-instantaneous Doppler (RID) image [2]. Restricted by real-time and feasibility requirements [3], most existing micro-motion studies use narrow-band radar. In this case, the separation of the m-D curves is roughly divided into three categories. The first is dictionary learning [4], which involves high-dimensional parameter search and has high computational complexity. The second is the decomposition class, such as empirical mode decomposition (EMD) [5,6], ensemble empirical mode decomposition (EEMD) [7], chirplet decomposition [8], and independent component analysis (ICA) [9]. However, the m-D signals of each scattering center have severe overlap in the time–frequency domain (TFD), which reduces the effectiveness of this type of algorithm. The third is the time–frequency image domain method, which extracts the m-D curves from the m-D spectrum obtained by the time–frequency analysis technology [10]. Therefore, its performance depends on a high-precision time–frequency representation (TFR) and a robust ridge detection and association algorithm.
Classical time–frequency analysis (TFA) methods are mainly divided into linear TFA methods and nonlinear TFA methods. Linear TFA methods, such as the short-time Fourier transform (STFT) [11] and S-transform (ST) [12], are limited by Heisenberg’s uncertainty principle, and the time–frequency resolution is low. Reference [13] proposed a multi-order short-time fractional Fourier transform (STFRFT) time–frequency algorithm using the time–frequency characteristics of each micro-Doppler component signal, which has better performance in micro-Doppler time–frequency analysis. However, there is a high computational complexity of order search. The nonlinear TFA methods have serious cross-term interference, such as Wigner–Ville distribution (WVD) [14]. Apart from this, the parametric TFA [15,16] and demodulation TFA [17] aim to improve the energy concentration by demodulating time-varying signals according to an extended parametric mathematical model, the computational complexity is high, and Heisenberg’s uncertainty principle must limit their results. Recently, an attempt has been made to enhance the TFR through a deep neural network model named TFA-Net, which achieved better results [18]. However, the signal cannot be reconstructed, and there are still many limitations in the practical application of the model. Some advanced post-processing methods are effective measures to solve the low time–frequency resolution, such as the synchrosqueezing transform (SST) [19], second-order SST (SST2) [20], synchroextracting transform (SET) [21], second-order SET (SET2), and their higher-order forms [22,23]. These TFAs cannot extract the TF coefficients on the instantaneous frequency (IF) trajectories when dealing with signals with strongly nonlinear IFs or relatively close IFs. The synchro-reassigning transform (SRT) [24] can reassign these TF coefficients to the IF trajectory by the derivative of the amplitude spectrum and the three-step selection rule, which eliminates the repeated assignment problem (RAP). However, the estimated IF is a first-order approximation to the signal at each transient, so the SRT cannot adapt to strongly time-varying signals.
M-D curve separation can use the birth–death of IF (BDIF) [25] estimation strategy on the TFD or the local peak detection and component linking (LPDCL) algorithm [26]. In the presence of overlapping components in the TFD, it may lead to tracking errors and interleaving between components. The Viterbi algorithm can separate the signals sequentially according to the energy under a low signal-to-noise ratio (SNR) [10]. However, it only considers the time–frequency amplitude and the absolute frequency change, and there may be the problem of incorrect signal correlation in the cross-region, as well as the computational complexity is high. In some literature, multi-component signals are regarded as target traces, and then, these traces are tracked by a modified Kalman filter [27], particle filter [28], high-order multi-frame track-before-detect [29], etc. Reference [30] proposed a simple ridge reorganization algorithm considering the relationship between the components. However, due to the selection of the thresholds, the algorithm has low adaptability. Parametric algorithms, such as the Wigner–Hough transform [31] and Viterbi–Hough joint algorithm [32], cannot be used due to unparameterized or higher-order parameters of the signal. In recent years, micro-Doppler feature extraction under complex conditions such as the occlusion effect [33] or discontinuous observation [34] has also been considered.
We propose an m-D curves’ separation algorithm of a space target based on the modified synchro-reassigning transform (MSRT) and ridge segment linking to solve the above problems. This paper mainly has the following two contributions:
  • This paper proposes the MSRT. We deduce the objective function W ( x ) of the MSRT, establish a two-step rearrangement rule to approximate the second-order local instantaneous frequency of the signal, and solve the problem that the SRT cannot be applied to strong time-varying signals.
  • We propose a novel ridge segment linking strategy. We make full use of the relationship between the components and the ridge information of the intersecting interval to realize the practical and robust association of the ridges of each component and solve the modal mismatch problem of the extracted ridges.
The remainder of this paper is organized as follows. Section 2 establishes the radar echo model of the cone-shaped target. Section 3 introduces the theory and implementation of the proposed MSRT. In Section 4, the details of the ridge segment linking algorithm are depicted. Section 5 verifies the effectiveness of the proposed method using simulation and electromagnetic calculation data. Finally, some conclusions are drawn in Section 6.

2. Radar Echo Model of Space Target

We built the micro-motion model of the space cone target shown in Figure 1. The initial distance between the radar and center of mass O is R 0 . Taking O as the coordinate origin, the symmetry axis is the z axis to establish a local coordinate system x , y , z , which remains stationary with the target. A reference coordinate system X , Y , Z is established with the re-entry direction as the positive Z axis, translating with the target. Its attitude remains stationary with the radar. During the flight of the target, while spinning around the O z axis, it will make a conical rotation around the O Z axis. The angle between the O z and O Z axis is the precession angle θ n when θ n remains stable during the movement; it is called precession, and when it swings within a specific range, it is called nutation.
The angle between the projection of the line of sight of the radar on the X O Y plane and the X axis is β n . Assume that the O z axis is on the Y O Z plane at the initial moment. For the position vector r b of any scattering point in the body coordinate at t, its position vector r R in the radar coordinate system can be expressed as
r R = R 0 n + R c t R n t R s t r b ,
where n = sin θ n cos β n , sin θ n sin β n , cos θ n is the unit direction vector of the line of sight for the radar in the reference coordinate system and R s t , R n t , and R c ( t ) are the transformation matrix of the spin, the initial rotation of the nutation angle, and the cone spin, respectively.
Assume that the angular velocity of the spin is ω s and the angular velocity of the cone spin is ω c . The average of the cone spin angle is θ c ; the angular velocity of the wobble for the symmetry axis is ω w ; the wobble amplitude is θ w ; the initial wobble phase is ψ w . According to Rodrigues’ rotation formula, we can obtain
R s ( t ) = I + e s sin ω s t + e s · e s 1 cos ω s t ,
R n ( t ) = I + e n sin θ t + e n · e n 1 cos θ t ,
R c ( t ) = I + e c sin ω c t + e c · e c 1 cos ω c t ,
where θ t and the oblique symmetric matrix are, respectively,
θ t = θ c + θ w cos ω w t + ψ w ,
e s = 0 1 0 1 0 0 0 0 0 ,
e n = 0 0 0 0 0 1 0 1 0 ,
e c = 0 1 0 1 0 0 0 0 0 .
A flat-bottomed cone-type target generally has three scattering centers, the ideal scattering center A at the top of the cone and the sliding scattering centers B and C at the edge of the radar incident plane and the bottom of the cone. The distances from the center of mass to the top and the bottom are h 1 and h 2 , and the direction vector of the spin axis in the reference coordinate system n s = R c t R n t R s t 0 , 0 , 1 T . Based on the geometric relationship, the distance vectors from the three scattering centers to the radar are calculated by
r A = R 0 n + h 1 n s ,
r B = R 0 n h 2 n s + r n × n s × n s ,
r c = R 0 n h 2 n s r n × n s × n s .
Due to the occlusion effect, some scattering does not contribute to the radar echo under a certain angle of view. The occlusion of the scattering center depends on the angle of view β and the semi-cone angle γ . The specific correspondence is shown in Table 1, where Y (or N) indicate that it is (or is not) occluded. In addition, the change rule of the angle of view is
β t = a cos n T n s .
Finally, the narrowband radar echo of the cone target is
S t = I N t σ I exp j 4 π f 0 r I C ,
where f 0 represents the carrier frequency of the transmitted signal, C is the speed of light, I N t represents the set of the scattering that is not occluded at time t, and · is the modulo operation.

3. The Proposed MSRT

The echo signal defined in the previous section behaves as a non-stationary multi-component signal, the modes of which will be reconstructed in this section through the MSRT. Firstly, we introduce the basic theory of the SRT. Secondly, the relationship between the SRT and the instantaneous frequency estimation operator is analyzed, and the basic idea of the second-order SRT is explained in detail. Finally, the implementation of the MSRT is given.

3.1. The Basic Theory of the SRT

For a multicomponent signal x t ,
x t = k = 1 K A k t exp j 2 π ϕ k t ,
where A k t and ϕ k ( t ) represent the instantaneous amplitude (IA) and instantaneous frequency (IF) of the k-th mode, respectively. Its STFT is defined as
V x g ( ω , t ) = x ( τ ) g ( τ t ) e j ω τ d τ ,
where τ and ω are the integral variable with respect to time and angular frequency, the symbol “∗” indicates the complex conjugate, and g ( t ) L 2 ( R ) is a real Gaussian window function.
By extracting the local maximum of amplitude spectrum P = | V x g ( ω , t ) | , we can obtain the approximate ideal time–frequency representation; P = | V x g ( ω , t ) | satisfies
P ω = δ V x g ( ω , t ) δ ω = 0 2 P ω 2 = δ 2 V x g ( ω , t ) δ ω 2 < 0 .
Then, the SRT is defined as
S R T η , t = V x g ( η , t ) s a t i s f y i n g ( 16 ) 0 n o t s a t i s f y i n g ( 16 ) .
and the mode x k t can be reconstructed by
x k ( t ) 1 g ^ 0 η ϕ k ( t ) < d S R T ( η , t ) d η ,
where g ^ is the Fourier transform of the window function g and ϕ k ( t ) is the instantaneous frequency of the k-th component at time t.
Next, we analyze the essence of the SRT. In [21], the two-dimensional IF can be obtained by
ω ˜ x ( η , t ) = t V x g ( η , t ) 2 i π V x g ( η , t ) = η 1 2 i π V x g ( η , t ) V x g ( η , t ) ,
where η is the integral variable with respect to frequency and V x g ( η , t ) is the STFT of the signal x ( t ) when the window function is the derivative of g ( t ) to t. Then, the IF estimation operator is defined as
ω ^ x ( η , t ) = ω ˜ x ( η , t ) = η 1 2 i π V x g ( η , t ) V x g ( η , t ) .
When the window function is g t = e π t 2 σ 2 , then the first derivative of g t to time g t = 2 π σ 2 t g t , ω ^ x ( η , t ) has a more concise expression.
ω ^ x = η + 1 σ 2 V x t g V x g .
Next, we analyze the characteristics of the zeros of ω ^ x η . ω ^ x η can be rewritten as
ω ^ x η = 1 σ 2 V x t g V x g = η V x g 2 i π σ 2 V x g = 1 4 π σ 2 η V x g 2 V x g 2 .
This shows that the points η 0 , t corresponding to the zeros satisfy η V x g t , η 0 2 = 0 and ω ^ x can estimate the local maximum of V x g t , η 0 2 . We can easily obtain η V x g t , η 0 = 0 ; it can be found that the essence of the SRT is to find the local maximum of V x g t , η 0 , that is the zeros of ω ^ x η . While the SRT completely eliminates the RAP through the three-step rule, it also has the limitation that the signal is required to be a pure harmonic signal.

3.2. MSRT

In order to overcome the limitation that the first-order IF estimator ω ^ x ( η , t ) is only suitable for pure harmonic signals, Reference [19] introduces the complex time delay operator t ˜ f ( t , η ) and the complex frequency modulation operator q ˜ f , which are
t ˜ x ( t , η ) = t η V x g ( t , η ) 2 i π V x g ( t , η ) = t + V x t g ( t , η ) V x g ( t , η ) ,
and
q ˜ f = η ω ˜ x η t ˜ x = 1 2 i π V x g V x t g V x g V x t g V x g 2 V x g V x t 2 g V x t g 2 ,
where q ˜ f η , t is the estimate of ϕ η , t and V x g , V x t g , and V x t 2 g are the STFTs of x ( t ) when the window functions are g , t g , and t 2 g , respectively.
Then, the second-order complex IF ω ˜ x [ 2 ] is defined by
ω ˜ x [ 2 ] = ω ˜ x + q ˜ x × t t ˜ x if η t ˜ x 0 ω ˜ x otherwise ,
and the second-order IF estimator ω ^ x [ 2 ] = ω ˜ x [ 2 ] .
ω ^ x [ 2 ] = ω ^ x + 1 2 i π V x g 2 + V x g V x t g V x g V x t g V x g V x t 2 g V x t g 2 V x t g V x g = ω ^ x + 1 2 i π V x g V x t g V x g V x t 2 g V x t g 2 1 σ 2 V x t g V x g = ω ^ x + 1 2 i π V x g V x t g V x g V x t 2 g V x t g 2 ω ^ x η = η + 1 2 π V x g V x t g V x g V x t 2 g V x t g 2
ω ^ x [ 2 ] η is equivalent to
ω ^ x [ 2 ] η = 1 2 π V x g V x t g V x g V x t 2 g V x t g 2 = 1 2 π V x g 2 V x t g V x t 2 g V x t g 2 V x g V x t g V x g V x t 2 g V x t g 2 2 = 1 2 π V x g 2 V x t g η V x t g V x t g 2 V x g η V x g V x g V x t 2 g V x t g 2 2 = 1 2 π V x g 2 η V x t g 2 V x t g 2 η V x g 2 V x g V x t 2 g V x t g 2 2
Assume W η , t = V x g 2 η V x t g 2 V x t g 2 η V x g 2 ; its discrete form is W m , n . Then, the point η 0 , t corresponding to the zeros of ω ^ f [ 2 ] η satisfies W η 0 , t = 0 . Based on the idea of the SRT, we can define a similar two-step rule to find the zeros of W η , t ,
W ( m , n ) < ξ W ( m 1 , n ) W ( m , n ) W ( m + 1 , n ) W ( m , n )
and
W ( m 1 , n ) > W ( m , n ) W ( m , n ) > W ( m + 1 , n ) ,
where ξ is the calculation error; we set ξ = 5 in our research.
The MSRT is defined as
M S R T m , n = V x g ( m , n ) s a t i s f y i n g ( 28 ) , a n d ( 29 ) 0 n o t s a t i s f y i n g ( 28 ) , a n d ( 29 ) ,

3.3. Implementation of the MSRT

The MSRT needs to calculate W η , t , where η V x t g 2 and η V x g 2 can be obtained by
η V x t g 2 = 2 V x t g η V x t g + 2 V x t g η V x t g = 2 V x t g η V x t g + 2 V x t g η V x t g ,
η V x g 2 = 2 V x g η V x g + 2 V x g η V x g ,
where η V x g = 2 i π V x t g , η V x t g = 2 i π V x t 2 g , so we need to calculate the STFTs V x g , V x t g , V x t 2 g when the window functions are g , t g , and t 2 g . We employ the standard fast discrete Fourier transform (FFT) for efficient computation. Consider a time domain discrete signal x n , n = 0 , , N 1 and a time domain discrete window function g m , m = 0 , , L 1 ; x [ n ] is the uniform sampling of the continuous signal x [ t ] at time t n = t 0 + n T , and T is the sampling period, then
V x g [ · , n ] = F F T x n g V x t g [ · , n ] = F F T x n t g V x t 2 g [ · , n ] = F F T x n t 2 g ,
where ⊙ represents elementwise multiplication.
Then, the MSRT can be calculated by Algorithm 1, and the computational complexity of the MSRT is about O 3 L N log 2 ( N ) .
Algorithm 1 Fast calculation of the discrete MSRT.
Input: x ( n )
Output: M S R T
  1:
Initialization: ξ , T x ( m , n ) 0
  2:
Calculate V x g , V x t g , and V x t 2 g based on (33)
  3:
η V x g 2 i π V x t g , η V x t g 2 i π V x t 2 g
  4:
Calculate η V x t g 2 and η V x g 2 based on (31) and (32)
  5:
W V x g 2 η V x t g 2 V x t g 2 η V x g 2
  6:
for m = 2 to L 1  do
  7:
   for  n = 1 to N do
  8:
     if (28) and (29) are satisfied then
  9:
         M S R T ( m , n ) = V x g ( m , n )
10:
     end if
11:
   end for
12:
end for

4. Ridge Segment Linking and Mode Reconstruction

M-D curves of the space target are severely overlapped in the time–frequency domain. A popular multi-ridge detection algorithm proposed in [19] can only be used for non-overlapping signals. The extracted ridges may not precisely match the modalities, but the ridges represent the global TF characteristics of the multi-component signal. Therefore, we can re-link the extracted time–frequency ridge line segments to achieve micro-Doppler curve separation and modal reconstruction of each scattering center of the space target. The process includes three steps: ridge detection, ridge segment linking, and mode reconstruction.

4.1. Ridge Detection

When the number of modes k is known, the ridges can be obtained by minimizing the following function E.
E ( ϕ ) = k = 1 K + T F R t , ϕ k ( t ) 2 d t + + λ · ϕ k ( t ) 2 + β · ϕ k ( t ) 2 d t ,
where k = 1 K t , ϕ k ( t ) is an estimate of the IF trajectory in the TF plane and λ and β are the two regularization parameters. In this paper, the algorithm in [19] is used to extract the micro-Doppler curve, and several random initializations are used to improve the robustness. When a component is present, its surrounding energy is set to zero and stops when the energy is below a certain threshold.

4.2. Ridge Segment Linking

First, determine the intersection interval t s , t e according to the distance between the ridges, and merge the adjacent intervals. Considering the distortion of the IF of each mode in the intersection interval, we discard the ridges in the intersection interval. Then, relink the ridge segment according to the rate of change of the ridges. The rate of change k l m , k r n , and k l r m n at t s , t e , and t s + t e / 2 can be approximated by
k l m = I F l m t s I F l m t s Δ t / Δ t k r n = I F r n t e + Δ t I F r n t e / Δ t k l r m n = I F m n t s + t e 2 I F m n t s + t e 2 Δ t / Δ t ,
where m , n = 1 , , I . I is the number of modes; Δ t represents the sampling time interval; I F l m and I F r n represent the m-th ridge on the left and the n-th ridge on the right, respectively; I F m n can be obtained by cubic spline interpolation using I F m and I F n . In practical applications, the average change rate with ridge segments within a specific time range increases the robustness. The matching degree matrix C m of different ridge segments is established to relink the ridges before and after the intersecting interval.
C m ( m , n ) = max k l m k r n , k r n k l r m n , k l r m n k l m max k l m , k r n , k l r m n .
Then, we can determine the correct connection pair ( p , q ) belonging to the same mode by minimizing the matching matrix C m .
p , q = arg min ( m , n ) 1 , 2 , , I C m m , n .
This process is clearly shown in Figure 2.

4.3. Mode Reconstruction

The MSRT is an estimate of the IF of a second-order local approximation of the signal. Without loss of generality, the signal is locally approximated at each transient with a Gaussian-modulated linear chirp signal x t with a Gaussian signal in amplitude and a second-order polynomial in phase.
x t = A x e t t 0 2 / 2 s 2 e j 2 π a + b t + 1 2 c t 2 ,
where A x represents the amplitude, t 0 represents the time center, and s is the standard deviation of the time distribution of the Gaussian window function.
According to δ V x g / δ t and δ V x g / δ t , we can deduce that
V x g = q V x t g q t + p j ω V x g ,
V x g = q V x t g q t + p j w V x g .
where p = t 0 / s 2 + j b and q = 1 / s 2 j c . g and g represent the first and second derivatives of the window function g ( t ) with respect to time t. V x g , V x t g , and V x t g , respectively, represent the short-time Fourier transform of x t using the window function g t , t g t and t g t . Combining (39) and (40), we can obtain
q t + p = V x g V x t g V x g V x t g V x g V x t g V x g V x t g + j w q = V x g V x g V x g V x g V x g V x t g V x g V x t g ,
Furthermore, the signal x ( t ) can be reconstructed by the M S R T result.
x ( t ) = M S R T ω , t s 2 π 4 q + 1 s 2 e x p 1 2 Re q t + p q + 1 s 2 2

5. Simulation and Verification

The proposed algorithm obtains the TFR of the signal through the MSRT and then realizes the practical separation of multi-component echo signals through ridge extraction and ridge segment linking. In this section, the effectiveness and robustness of the algorithm are verified by the simulation data in Section 2, bat echo, and electromagnetic calculation data. The simulation and verification were performed on AMD Ryzen 7 4800H which manufactured by TSMC, Hsinchu, Taiwan and Matlab 2022a.

5.1. Detection Process and Parameter Setting

Considering real-time and feasibility requirements, single-frequency continuous wave radar is used to observe the space target. Figure 3a shows the detection scene. After the X-band radar has selected the target through the wave gate, the proposed MSRT method is used to obtain the m-D spectrum. After ridge extraction and ridge segment connection, the m-D curve of each scatter is obtained, which is used to extract micro-motion features for target identification such as precession frequency, precession angle, etc. The overall processing flow of the radar signal is shown in Figure 3b.
We considered a flat-bottomed cone-shaped target with a rotationally symmetric structure, as shown in Figure 4a. Its height is 1.6 m, and the radius of the bottom surface is 0.4 m. The target has three scattering centers, and the radar echo model was established according to Section 2. The carrier frequency of the radar transmit signal is 10 GHz; the bandwidth was not considered. The simulation assumes that the translation has been fully compensated by some method [35]. For the complete observation of the m-D frequency and the Doppler frequency brought by the residual translation, the sampling rate was selected to be 1000 Hz. The time–frequency window length is only obtained by the optimization method at the beginning of the simulation. When the radar carrier frequency and accumulation time are determined, the variation range of the optimal window length is small. The specific simulation parameters are shown in Table 2.

5.2. Verification on Simulation Data

According to the theoretical analysis, when the horizon angle is 60°, the scattering center C is always in the occlusion state, and the signal-to-noise ratio S N R = 15 dB, where the definition of SNR is shown in (43). The simulation time was 1 s. To ensure the accuracy of the comparison results, the source codes of the comparison TFA were refactored in the same way. The STFT result of the target echo in Figure 4b is consistent with the theoretical analysis.
SNR = 20 log 10 x 2 x n 2 ,
where 2 stands for the l 2 -norm, x represents the radar echo, and x n represents the noise signal.
Figure 5a–f show the m-D spectrum of the radar echo when the target precesses, which are the results of the SST, SST2, SET, SET2, SRT, and MSRT, respectively. It can be seen that the energy concentration of these post-processing methods is significantly better than that of the STFT. Compared with the SST, SET, and SST2, SRT and SRT2 can obtain a higher energy concentration; this is because the three-step rearrangement rule of the SRT can accurately identifytime–frequency ridges, and its concentration is comparable to that of the SET2.
To quantitatively evaluate the performance of the proposed MSRT, we compared the Rényi entropy and running time of different time–frequency representations. The Rényi entropy is widely used in the evaluation of time–frequency energy concentration. The smaller the entropy value is, the higher the time–frequency concentration is. The third-order Rényi entropy can be calculated by (44). Table 3 shows the Rényi entropy and running time of different TFA methods. The three-step rule of the SRT can accurately locate the local maximum represented by the amplitude spectrum and has a particular ability to suppress noise, so the Rényi entropy of the SRT is the smallest. The MSRT is slightly higher than the SRT. The reason is that the constraint of the two-step rule is weaker than that of the SRT, and the MSRT can estimate the second-order IF. For strong time-varying signals, the IF is not entirely concentrated in the local maximum of the amplitude spectrum. In addition, the MSRT can obtain similar estimation accuracy of the IF to that of the SET2 and SST2. The MSRT only needs to calculate the STFTs of three window functions. At the same time, its time complexity is much lower than that of the SST2 and SET2.
R = 1 2 log 2 R 2 | S ( t , ω ) | 3 d ω d t R 2 | S ( t , ω ) | d ω d t ,
where S ( t , ω ) represents the results of the TFA methods such as the STFT and SET.
Figure 6a,b show the m-D curve extraction results and relinking results, respectively. The ridge extraction algorithm cannot be applied to overlapping multi-component signals, and the extracted ridges have correlation errors in the intersection interval. Still, the global features of the m-D curve can be accurately extracted. The ridge segment connection can effectively correlate the Doppler curves of each scattering center, which will be helpful for the subsequent estimation of the fretting parameters. Figure 6c,d show the difference between the reconstructed signals of the scattering centers A and B and the original theoretical signals, and the signal-to-noise ratio of the reconstructed signals can reach more than 14 dB.
For the convenience of calculation, the time–frequency representation of the scattering center A was used to test the instantaneous frequency estimation accuracy of the SRT and MSRT. We used the root-mean-squared error between the estimated IF and the theoretical IF to evaluate the IF estimation accuracy and the output SNR to evaluate the signal reconstruction ability. According to Figure 7 and Table 4, it can be seen that the IF estimated by the SRT is entirely consistent with the local maximum of the amplitude spectrum. Near the maximum value of the iF, the MSRT and SRT have similar errors. At other times, the IF estimation based on the MSRT has higher accuracy. Since the second-order approximation of the signal is considered during modal reconstruction, the output SNR is far higher than the SRT.
Whether the m-D curves can be correctly correlated will seriously affect the micro-motion feature extraction process of space targets. Next, two-component and three-component signals were used to compare and verify the separation algorithm of Viterbi [10], the RPRG [30], and the proposed ridge segment linking. According to the analysis of the results shown in Figure 8a–c, the Viterbi algorithm only considers the absolute change of the amplitude and frequency of the spectrum. When the energy of each component is similar, the algorithm cannot effectively track the current Doppler trajectory at the intersection of the different modes. The RPRG is effective when appropriate parameters are selected. It can effectively correlate each mode in Figure 8b. However, when the same parameters are applied to the three-component signal in Figure 8e, the correlation error occurs, and the number of signal components is required for preliminary information. The proposed ridge segment linking makes full use of the information of the intersection interval and has a more robust constraint matching matrix to make ridge segment linking more stable.
Finally, we verified the noise robustness of the MSRT. We added different levels of Gaussian white noise to the radar echo signal, and the SNR ranged from 0 dB to 20 dB. Figure 9a shows the Rényi entropy calculated by different TFA methods. It can be seen that under different SNRs, the Rényi entropy of the MSRT is only higher than that of the SRT. This is because the MSRT rearranges the TF energy through the two-step rule, which can eliminate the RAP, effectively avoiding many noise ridges and having higher noise robustness. In Figure 9b, the output signal-to-noise ratios of the SST2-, SET2-, and MSRT-reconstructed signals are much higher than other methods. Under different input signal-to-noise ratio conditions, the output SNR of the MSRT reconstruction results has less fluctuation. With the decrease of the SNR, the Rényi entropy and the output SNR decrease. This is because the noise has three influences on the entire process. First, the noise is randomly distributed in the whole T-F plane, which makes the energy concentration of the TFR result lower. Second, the noise causes a deviation between the ridge and the IF of the signal, resulting in a decrease in the output SNR of the reconstructed signal. Third, high-energy noise may lead to the failure of the ridge segment connection.

5.3. Application to Bat Echo

We used the bat echo [20] collected by Rice University to verify the adaptability of the MSRT proposed in this paper to the actual signal. The sampling rate of the bat echo used is 140 kHz, with a total of 400 sampling points. The waveform and spectrum are shown in Figure 10. It is not easy to understand how bats use their echoes to identify and locate targets.
Figure 11a shows the STFT of the bat echo. The nonlinear behaviors of bat echolocation precisely in the time–frequency domain are evident. The echo contains four modes, and the duration and frequency of each mode are easily available from the TFR. Figure 11b–e show the zooms of the STFT/MSRT/SST2 and SET2, respectively; the STFT has the lowest concentration, and the other three post-processing methods are significantly better than the STFT. The positioning accuracy of the MSRT for the instantaneous frequency of each mode is comparable to that of the SET2, but it introduces less noise.
In Figure 12, the IFs of each mode are estimated from the MSRT using the ridge extraction algorithm (see Figure 12a). Each mode was reconstructed from its TF coefficients. The reconstructed signals of each mode were superimposed to obtain the reconstructed echoes with tiny reconstruction errors, indicating that the MSRT has a good signal reconstruction capability, i.e., reversibility.

5.4. Application to Electromagnetic Calculation Data

When the cone-shaped target was nutated, we verified the proposed algorithm with electromagnetic calculation data. The cone rotation frequency was 1 Hz, and the simulation time was 5 s. The other simulation parameters are shown in Table 2. First, CADFEKO was used to establish the cone-shaped target shown in Figure 13a, and the static echo at full attitude was calculated. According to the Section 2, the field of view angle change law when the target is nutated (see in Figure 13b) was obtained. Then, the static data at the corresponding posture were extracted from the full-pose static echo to form the dynamic echo. In order to avoid the influence of the occlusion effect on the scattering intensity, the amplitude of the echo signal was normalized, that is S n t = S t / S t . The short-time Fourier transform and MSRT were obtained as shown in Figure 13c,d, respectively, and it can be seen that the C scattering center is in an occluded state, which is consistent with the theoretical analysis. The MSRT has a higher time–frequency energy concentration. Figure 13e,f show that the proposed algorithm can extract the complete micro-Doppler curves corresponding to A and B and achieve the correct correlation.
For more complex objects, the radar echo behaves as a multi-component signal. The proposed algorithm obtains the TFR of the echo through the MSRT, and the ridge segment linking is used to associate and reorganize the ridges of each mode. As long as the various modes of the target echo signal do not completely overlap in the time–frequency domain, that is the requirements of separability are met [17], the algorithm is still applicable.

6. Conclusions

Aiming at the problems of the low precision of separation curves and easy association errors of overlapping m-D curves in the existing methods, we proposed an m-D curve separation algorithm with a space target based on the MSRT and ridge segment linking. Based on the idea of the SRT, the second-order frequency estimation operator is introduced to correct the SRT, which assists in rearranging the time–frequency coefficients of strongly time-varying signals. The MSRT can obtain a relatively high time–frequency concentration, strong noise robustness, and low computational complexity. Then, the m-D curve of the target scattering center is extracted from the MSRT, and the position of the cross-interval is adaptively determined according to the distance of each scattering center at each moment. Then, the m-D curve segments are correlated and reorganized by ridge segment linking. The simulation and experimental results verify the effectiveness of the algorithm. When the signal is completely submerged in noise, the denoising method of the radar signal or the more robust TFA method can improve the effectiveness of the proposed algorithm.

Author Contributions

Conceptualization, D.Y. and X.W.; methodology, X.W.; software, J.L. and X.W.; validation, Z.P.; writing—original draft preparation, X.W.; writing—review and editing, D.Y. and J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key Research and Development Program of China grant number 2021YFC3090402, Defense Advance Research Funder grant number 514010402, and National Natural Science Foundation of China grant number 62171475.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, V.; Li, F.; Ho, S.S.; Wechsler, H. Micro-Doppler effect in radar: Phenomenon, model, and simulation study. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 2–21. [Google Scholar] [CrossRef]
  2. Tian, X.; Bai, X.; Xue, R.; Qin, R.; Zhou, F. Fusion Recognition of Space Targets with Micro-Motion. IEEE Trans. Aerosp. Electron. Syst. 2022. [Google Scholar] [CrossRef]
  3. Chen, J.; Zhang, J.; Jin, Y.; Yu, H.; Liang, B.; Yang, D.G. Real-Time Processing of Spaceborne SAR Data With Nonlinear Trajectory Based on Variable PRF. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5205212. [Google Scholar] [CrossRef]
  4. Li, G.; Varshney, P.K. Micro-Doppler parameter estimation via parametric sparse representation and pruned orthogonal matching pursuit. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 4937–4948. [Google Scholar] [CrossRef]
  5. Zhao, Y.; Su, Y. The extraction of micro-Doppler signal with EMD algorithm for radar-based small UAVs’ detection. IEEE Trans. Instrum. Meas. 2019, 69, 929–940. [Google Scholar] [CrossRef]
  6. Wang, Z.; Luo, Y.; Li, K.; Yuan, H.; Zhang, Q. Micro-Doppler Parameters Extraction of Precession Cone-Shaped Targets Based on Rotating Antenna. Remote Sens. 2022, 14, 2549. [Google Scholar] [CrossRef]
  7. Zhu, N.; Hu, J.; Xu, S.; Wu, W.; Zhang, Y.; Chen, Z. Micro-Motion Parameter Extraction for Ballistic Missile with Wideband Radar Using Improved Ensemble EMD Method. Remote Sens. 2021, 13, 3545. [Google Scholar] [CrossRef]
  8. Zhang, W.; Xie, J.; Li, H. Micro-Doppler separation method in ISAR imaging based on short-time chirplet decomposition. J. Appl. Remote Sens. 2021, 15, 036515. [Google Scholar] [CrossRef]
  9. Choi, I.; Kang, K.; Kim, K.; Park, S. Use of ICA to Separate Micro-Doppler Signatures in ISAR Images of Aircraft That Has Fast-Rotating Parts. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 234–246. [Google Scholar] [CrossRef]
  10. Djurović, I.; Stanković, L. An algorithm for the Wigner distribution based instantaneous frequency estimation in a high noise environment. Signal Process. 2004, 84, 631–643. [Google Scholar] [CrossRef]
  11. Djurović, I.; Popović-Bugarin, V.; Simeunović, M. The STFT-based estimator of micro-Doppler parameters. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1273–1283. [Google Scholar] [CrossRef]
  12. Li, B.; Zhang, Z.; Zhu, X. Adaptive S-Transform with Chirp-Modulated Window and Its Synchroextracting Transform. Circuits Syst. Signal Process. 2021, 40, 5654–5681. [Google Scholar] [CrossRef]
  13. Pang, C.; Han, Y.; Hou, H.; Liu, S.; Zhang, N. Micro-Doppler Signal Time-Frequency Algorithm Based on STFRFT. Sensors 2016, 16, 1559. [Google Scholar] [CrossRef] [PubMed]
  14. Boashash, B.; Black, P. An efficient real-time implementation of the Wigner-Ville distribution. IEEE Trans. Acoust. Speech Signal Process. 1987, 35, 1611–1618. [Google Scholar] [CrossRef]
  15. Yang, Y.; Peng, Z.; Zhang, W.; Meng, G. Parameterised time–frequency analysis methods and their engineering applications: A review of recent advances. Mech. Syst. Signal Process. 2019, 119, 182–221. [Google Scholar] [CrossRef]
  16. Ai, X.; Xu, Z.; Wu, Q.; Liu, X.; Xiao, S. Parametric Representation and Application of Micro-Doppler Characteristics for Cone-Shaped Space Targets. IEEE Sens. J. 2019, 19, 11839–11849. [Google Scholar] [CrossRef]
  17. Meignen, S.; Pham, D.H.; McLaughlin, S. On demodulation, ridge detection, and synchrosqueezing for multicomponent signals. IEEE Trans. Signal Process. 2017, 65, 2093–2103. [Google Scholar] [CrossRef] [Green Version]
  18. Pan, P.; Zhang, Y.; Deng, Z.; Fan, S.; Huang, X. TFA-Net: A Deep Learning-Based Time-Frequency Analysis Tool. IEEE Trans. Neural Netw. Learn. Syst. 2022. [Google Scholar] [CrossRef]
  19. Meignen, S.; Singh, N. Analysis of Reassignment Operators Used in Synchrosqueezing Transforms: With an Application to Instantaneous Frequency Estimation. IEEE Trans. Signal Process. 2021, 70, 216–227. [Google Scholar] [CrossRef]
  20. Oberlin, T.; Meignen, S.; Perrier, V. Second-Order Synchrosqueezing Transform or Invertible Reassignment? Towards Ideal Time-Frequency Representations. IEEE Trans. Signal Process. 2015, 63, 1335–1344. [Google Scholar] [CrossRef] [Green Version]
  21. Yu, G.; Yu, M.; Xu, C. Synchroextracting transform. IEEE Trans. Ind. Electron. 2017, 64, 8042–8054. [Google Scholar] [CrossRef]
  22. Bao, W.; Li, F.; Tu, X.; Hu, Y.; He, Z. Second-order synchroextracting transform with application to fault diagnosis. IEEE Trans. Instrum. Meas. 2020, 70, 3509409. [Google Scholar] [CrossRef]
  23. Pham, D.H.; Meignen, S. High-Order Synchrosqueezing Transform for Multicomponent Signals Analysis—With an Application to Gravitational-Wave Signal. IEEE Trans. Signal Process. 2017, 65, 3168–3178. [Google Scholar] [CrossRef] [Green Version]
  24. Li, M.; Wang, T.; Kong, Y.; Chu, F. Synchro-reassigning transform for instantaneous frequency estimation and signal reconstruction. IEEE Trans. Ind. Electron. 2021, 69, 7263–7274. [Google Scholar] [CrossRef]
  25. McAulay, R.; Quatieri, T. Speech analysis/synthesis based on a sinusoidal representation. IEEE Trans. Acoust. Speech Signal Process. 1986, 34, 744–754. [Google Scholar] [CrossRef] [Green Version]
  26. Rankine, L.; Mesbah, M.; Boashash, B. IF estimation for multicomponent signals using image processing techniques in the time–frequency domain. Signal Process. 2007, 87, 1234–1250. [Google Scholar] [CrossRef] [Green Version]
  27. Ren, K.; Du, L.; Lu, X.; Zhuo, Z.; Li, L. Instantaneous Frequency Estimation Based on Modified Kalman Filter for Cone-Shaped Target. Remote Sens. 2020, 12, 2766. [Google Scholar] [CrossRef]
  28. Hong, L.; Wang, X.; Liu, S. Micro-Doppler curves extraction based on high-order particle filter track-before detect. IEEE Geosci. Remote Sens. Lett. 2019, 16, 1550–1554. [Google Scholar] [CrossRef]
  29. Zhu, N.; Xu, S.; Li, C.; Hu, J.; Fan, X.; Wu, W.; Chen, Z. An Improved Phase-Derived Range Method Based on High-Order Multi-Frame Track-Before-Detect for Warhead Detection. Remote Sens. 2022, 14, 29. [Google Scholar] [CrossRef]
  30. Chen, S.; Dong, X.; Xing, G.; Peng, Z.; Zhang, W.; Meng, G. Separation of overlapped non-stationary signals by ridge path regrouping and intrinsic chirp component decomposition. IEEE Sens. J. 2017, 17, 5994–6005. [Google Scholar] [CrossRef]
  31. Barbarossa, S. Analysis of multicomponent LFM signals by a combined Wigner-Hough transform. IEEE Trans. Signal Process. 1995, 43, 1511–1515. [Google Scholar] [CrossRef]
  32. Ding, Y.; Liu, R.; She, Y.; Jin, B.; Peng, Y. Micro-Doppler Trajectory Estimation of Human Movers by Viterbi–Hough Joint Algorithm. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5113111. [Google Scholar] [CrossRef]
  33. Zhou, Y.; Chen, Z.; Zhang, L.; Xiao, J. Micro-Doppler Curves Extraction and Parameters Estimation for Cone-Shaped Target with Occlusion Effect. IEEE Sens. J. 2018, 18, 2892–2902. [Google Scholar] [CrossRef]
  34. Shi, Y.; Jiu, B.; Liu, H. Optimization-Based Discontinuous Observation Strategy for Micro-Doppler Signature Extraction of Space Cone Targets. IEEE Access 2019, 7, 58915–58929. [Google Scholar] [CrossRef]
  35. Chen, J.; Xing, M.; Yu, H.; Liang, B.; Peng, J.; Sun, G.C. Motion Compensation/Autofocus in Airborne Synthetic Aperture Radar: A Review. IEEE Geosci. Remote Sens. Mag. 2022, 10, 185–206. [Google Scholar] [CrossRef]
Figure 1. Micro-motion model of cone target with three scattering centers A, B and C.
Figure 1. Micro-motion model of cone target with three scattering centers A, B and C.
Remotesensing 14 03691 g001
Figure 2. Ridge-segment linking process. (a) Determine the intersection interval. (b) Calculate the rate of change of the ridges. (c) Ridge-segment linking.
Figure 2. Ridge-segment linking process. (a) Determine the intersection interval. (b) Calculate the rate of change of the ridges. (c) Ridge-segment linking.
Remotesensing 14 03691 g002
Figure 3. Detection scene and signal processing flow. (a) Detection scene of space target with three scattering centers A, B and C. (b) Processing flow.
Figure 3. Detection scene and signal processing flow. (a) Detection scene of space target with three scattering centers A, B and C. (b) Processing flow.
Remotesensing 14 03691 g003
Figure 4. Target model and TFD. (a) A flat-bottomed cone-shaped target with a rotationally symmetric structure. A, B and C are the three scattering centers of the cone-shaped target. (b) STFT result.
Figure 4. Target model and TFD. (a) A flat-bottomed cone-shaped target with a rotationally symmetric structure. A, B and C are the three scattering centers of the cone-shaped target. (b) STFT result.
Remotesensing 14 03691 g004
Figure 5. Micro-Doppler spectrum through different time–frequency analysis methods. (a) SST. (b) SET. (c) SRT. (d) SST2. (e) SET2. (f) MSRT.
Figure 5. Micro-Doppler spectrum through different time–frequency analysis methods. (a) SST. (b) SET. (c) SRT. (d) SST2. (e) SET2. (f) MSRT.
Remotesensing 14 03691 g005
Figure 6. Ridge detection and mode reconstruction results. (a) Ridge detection result. (b) Ridge segment linking. (c) Mode reconstruction result of scatter A. (d) Mode reconstruction result of scatter B.
Figure 6. Ridge detection and mode reconstruction results. (a) Ridge detection result. (b) Ridge segment linking. (c) Mode reconstruction result of scatter A. (d) Mode reconstruction result of scatter B.
Remotesensing 14 03691 g006
Figure 7. IF estimates and errors. (a) The instantaneous frequency estimate of the scatter A by different TFAs. (b) IF estimation error.
Figure 7. IF estimates and errors. (a) The instantaneous frequency estimate of the scatter A by different TFAs. (b) IF estimation error.
Remotesensing 14 03691 g007
Figure 8. Ridge correlation results by different methods. (ac) Result of ridge segment linking, the RPRG, and the Viterbi algorithm for a two-component signal. (df) Result of ridge segment linking, the RPRG, and the Viterbi algorithm for a three-component signal.
Figure 8. Ridge correlation results by different methods. (ac) Result of ridge segment linking, the RPRG, and the Viterbi algorithm for a two-component signal. (df) Result of ridge segment linking, the RPRG, and the Viterbi algorithm for a three-component signal.
Remotesensing 14 03691 g008
Figure 9. Noise robustness verification results of the MSRT. (a) Rényi entropy under different SNR conditions. (b) Output SNR of reconstruction results under different SNR conditions.
Figure 9. Noise robustness verification results of the MSRT. (a) Rényi entropy under different SNR conditions. (b) Output SNR of reconstruction results under different SNR conditions.
Remotesensing 14 03691 g009
Figure 10. The base information of the bat signal. (a) Waveform. (b) Spectrum.
Figure 10. The base information of the bat signal. (a) Waveform. (b) Spectrum.
Remotesensing 14 03691 g010
Figure 11. Results of different TFAs. (a) STFT result. (b) Zoom of the STFT result. (c) Zoom of the MSRT result. (d) Zoom of the SST2 result. (e) Zoom of the SET2 result.
Figure 11. Results of different TFAs. (a) STFT result. (b) Zoom of the STFT result. (c) Zoom of the MSRT result. (d) Zoom of the SST2 result. (e) Zoom of the SET2 result.
Remotesensing 14 03691 g011
Figure 12. The estimated IF and mode reconstruction of bat echo. (a) The estimated IF trajectories based on the MSRT. (b) Mode reconstruction results and error.
Figure 12. The estimated IF and mode reconstruction of bat echo. (a) The estimated IF trajectories based on the MSRT. (b) Mode reconstruction results and error.
Remotesensing 14 03691 g012
Figure 13. Results of the proposed algorithm on electromagnetic computational data. (a) Cone-shaped target model with three scattering centers. (b) The field of view angle change law. (c) Result by the STFT. (d) Result by the MSRT. (e) Ridge extraction. (f) Ridge segment linking.
Figure 13. Results of the proposed algorithm on electromagnetic computational data. (a) Cone-shaped target model with three scattering centers. (b) The field of view angle change law. (c) Result by the STFT. (d) Result by the MSRT. (e) Ridge extraction. (f) Ridge segment linking.
Remotesensing 14 03691 g013aRemotesensing 14 03691 g013b
Table 1. Occlusion of scattering centers under different semi-cone angles γ and the angle of view β .
Table 1. Occlusion of scattering centers under different semi-cone angles γ and the angle of view β .
β t 0 β t < γ γ β t < π / 2 π / 2 β t < π γ π γ β t < π
ANNNY
BNNNN
CNYNN
Table 2. Simulation parameter settings.
Table 2. Simulation parameter settings.
ParameterValue
Carrier frequency f 0 10 GHz
Sampling rate f s 1000 Hz
Sampling time t1 s
Average of viewing angle θ n 60
Azimuth angle β n 90
Angular velocity of the cone ω c 5 π rad/s
Cone angle θ c 8
Angular velocity of the wobble ω w 2 π rad/s
Amplitude of the wobble θ w 2
Semi-cone angle γ 14
Table 3. Comparison with different time–frequency analysis methods.
Table 3. Comparison with different time–frequency analysis methods.
TFRSTFTSSTSST2SETSET2SRTMSRT
Rényi19.256917.166414.746814.539113.145912.559413.0054
Time (s)1.6451.18642.85001.03152.28980.38940.9914
RMSE2.15613.15881.59842.15611.55462.15781.5581
Table 4. RMSE between the estimated IF and the theoretical IF and the output SNR of the reconstructed signal.
Table 4. RMSE between the estimated IF and the theoretical IF and the output SNR of the reconstructed signal.
TFRSTFTSSTSST2SETSET2SRTMSRT
RMSE2.15613.15881.59842.15611.55462.15781.5581
Out SNR (dB)4.80724.794314.13883.378314.32974.807214.2471
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, D.; Wang, X.; Li, J.; Peng, Z. Micro-Doppler Curves Extraction of Space Target Based on Modified Synchro-Reassigning Transform and Ridge Segment Linking. Remote Sens. 2022, 14, 3691. https://doi.org/10.3390/rs14153691

AMA Style

Yang D, Wang X, Li J, Peng Z. Micro-Doppler Curves Extraction of Space Target Based on Modified Synchro-Reassigning Transform and Ridge Segment Linking. Remote Sensing. 2022; 14(15):3691. https://doi.org/10.3390/rs14153691

Chicago/Turabian Style

Yang, Degui, Xing Wang, Jin Li, and Zhenghong Peng. 2022. "Micro-Doppler Curves Extraction of Space Target Based on Modified Synchro-Reassigning Transform and Ridge Segment Linking" Remote Sensing 14, no. 15: 3691. https://doi.org/10.3390/rs14153691

APA Style

Yang, D., Wang, X., Li, J., & Peng, Z. (2022). Micro-Doppler Curves Extraction of Space Target Based on Modified Synchro-Reassigning Transform and Ridge Segment Linking. Remote Sensing, 14(15), 3691. https://doi.org/10.3390/rs14153691

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