Next Article in Journal
High-Efficiency Forward Modeling of Gravitational Fields in Spherical Harmonic Domain with Application to Lunar Topography Correction
Next Article in Special Issue
Enhancing Digital Twins with Human Movement Data: A Comparative Study of Lidar-Based Tracking Methods
Previous Article in Journal
AMHFN: Aggregation Multi-Hierarchical Feature Network for Hyperspectral Image Classification
Previous Article in Special Issue
Mapping Building Heights at Large Scales Using Sentinel-1 Radar Imagery and Nighttime Light Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Doppler-Spread Space Target Detection Based on Overlapping Group Shrinkage and Order Statistics

School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(18), 3413; https://doi.org/10.3390/rs16183413
Submission received: 11 July 2024 / Revised: 11 September 2024 / Accepted: 11 September 2024 / Published: 13 September 2024
(This article belongs to the Special Issue Remote Sensing: 15th Anniversary)

Abstract

:
The Doppler-spread problem is commonly encountered in space target observation scenarios using ground-based radar when prolonged coherent integration techniques are utilized. Even when the translational motion is accurately compensated, the phase resulting from changes in the target observation attitude (TOA) still leads to extension of the target’s echo energy across multiple Doppler cells. In particular, as the TOA change undergoes multiple cycles within a coherent processing interval (CPI), the Doppler spectrum spreads into equidistant sparse line spectra, posing a substantial challenge for target detection. Aiming to address such problems, we propose a generalized likelihood ratio test based on overlapping group shrinkage denoising and order statistics (OGSos-GLRT) in this study. First, the Doppler domain signal is denoised according to its equidistant sparse characteristics, allowing for the recovery of Doppler cells where line spectra may be situated. Then, several of the largest Doppler cells are integrated into the GLRT for detection. An analytical expression for the false alarm probability of the proposed detector is also derived. Additionally, a modified OGSos-GLRT method is proposed to make decisions based on an increasing estimated number of line spectra (ENLS), thus increasing the robustness of OGSos-GLRT when the ENLS mismatches the actual value. Finally, Monte Carlo simulations confirm the effectiveness of the proposed detector, even at low signal-to-noise ratios (SNRs).

1. Introduction

Space surveillance plays a crucial role in military actions and space combat scenarios, holding a pivotal position in terms of space situational awareness (SSA) [1]. However, the detection of distant or small-sized space targets poses significant challenges for ground-based radar systems due to their weak signals, which often require long-time coherent integration to enhance the signal-to-noise ratio (SNR). Numerous studies based on coherent integration have focused on correcting the range migration (RM) and Doppler frequency migration (DFM) caused by translational motions to accumulate echo energy [2,3,4,5,6,7,8,9,10]. In reality, space debris has lost attitude and orbit control, exhibiting uncontrolled tumbling behavior [11]. The Delta 4 rocket body (NORAD ID: 40535) was observed exhibiting such tumbling by a long-baseline bistatic radar system using an L-band radio telescope network [12]. Spin, precession, and nutation are typical micro-motions in the midcourse phase of ballistic targets [13], all of which cause the target observation attitude (TOA) to vary relative to the radar line of sight. As the integration time increases, the impact of the TOA resulting from micro-motions or observation geometries on accumulation performance is non-negligible. Although targets may be unresolved in range cells in a narrowband radar system, irregular changes in the TOA could still induce the Doppler spread. On one hand, TOA changes become more prominent with a long coherent processing interval (CPI); on the other hand, the Doppler resolution improves with increased integration time [14]. Consequently, the target energy may disperse across multiple Doppler cells, necessitating secondary accumulation in the Doppler domain to enhance the detection capability. Assuming that translational motions are precisely compensated, the aforementioned issues arise when detecting Doppler-spread targets, for which theories related to high-resolution range-spread target detection are often adopted.
For spread target detection under the influence of Gaussian noise, if the target scattering characteristics are completely known, an ideal likelihood ratio detector can be constructed, enabling a theoretical analysis of the detection performance. However, the scattering information of targets is usually not fully available in practice, leading to unrealizable test statistics, in which case suboptimal detectors are sought. An energy detector [15], which is independent of prior information, has been used to perform the non-coherent integration of echoes within a certain distance range before detection, demonstrating effective performance in scenarios with evenly distributed energy. But, the energy detector also integrates noise cells that lack target echo signals. This is particularly problematic for targets with sparse energy distribution, where the improvement in signal-to-noise ratio (SNR) is reduced and the detection performance degrades. An M/N detector [16], which operates through binary comparisons, assuming that the number of points exceeding the first threshold within N data points is no less than M, generally exhibited superior detection capabilities in scenarios with sparsely distributed energy. Despite that, the detection performance of an M/N detector heavily depends on the second threshold selection, which will cause significant quantization errors when not appropriately selected. Dai et al. [17] re-arranged the data in the detection window according to their energy, and then integrated the range cells and carried out detection according to their order. This method also suffered from the setting of the detection threshold due to the complexity of the probability distribution of order statistics. In addition, the features of the cross-time–frequency distribution (TFD) of a pair of adjacent received signals were exploited in range-spread target detectors [18,19], which were shown to outperform traditional detectors. Consecutive high-resolution range profiles (HRRPs) have also been utilized for the detection of range-spread targets [20,21]. The detection methods related to the HRRPs leverage the stability of the range profile of wideband radar targets. However, when noise or clutter power is high, the SNR suffers significant loss, which results in the deterioration of detection performance.
In terms of adaptive detection, DeMaio [22] addressed the problem related to the polarimetric adaptive detection of range-spread targets with an unknown covariance matrix. Tang et al. [23] extensively investigated the adaptive detection of range-spread targets with noise covariance sharing block-diagonal structures, and proposed two categories of invariant detectors through adopting the invariance principle. Adaptive Rao and Wald tests were researched in the context of an unknown covariance matrix in [24]. Moreover, an adaptive double-threshold (DT) detection method was proposed in [25] that fully utilizes the features of HRRPs. However, this approach is computationally intensive and experiences performance degradation at low SNRs.
Furthermore, many detection algorithms have been developed for range-spread targets in scenes with various clutter. For a Gaussian clutter scene, a double-threshold constant false alarm rate (DT-CFAR) detector was proposed by exploiting the scatterer density and amplitude in [26]. Assuming that a set of secondary data is available, adaptive Rao and Wald tests were proposed in [27], resulting in a two-step detection strategy; according to this, an adaptive detector was also proposed by exploiting the persymmetry of the unknown covariance matrix in [28]. For a non-Gaussian clutter scene, Gerlach [29] proposed a scatterer density-dependent generalized likelihood ratio test (SDD-GLRT) using prior knowledge of the spatial scatterer density and the clutter covariance matrix. He et al. [30] proposed order statistics GLRT (OS-GLRT), which essentially assumes that the number of scatterers and secondary data is known. SDD-GLRT and OS-GLRT use prior knowledge of the spatial scatterer density, which is unattainable for targets in space surveillance that are always non-cooperative. Therefore, the detection performance will decrease when there is a mismatch in prior information. To avoid reliance on prior information, Ren et al. [31,32] proposed an adaptive scatterer estimation GLRT (ASCE-GLRT) by modeling the problem of scatterer estimation in terms of a sparse signal representation, and proposed an adaptive Doppler steering matrix estimation GLRT (ADSME-GLRT) using a sparse Bayesian optimization model, both with available secondary data. Both the methods perform well in detecting sparsely distributed signals but fail to fully exploit the structural information of the sparse signals.
In this work, we investigate the Doppler spread target detection of TOA changes over multiple periods within a CPI. Relevant targets include rapidly spinning space debris, precessional warheads, tumbling rocket bodies, and satellites. After compensating for translational motion and coherent integration, the echo energy is sparsely dispersed in equidistant Doppler cells, which poses significant challenges for target detection. To solve this problem, we first adopt a denoising algorithm. The overlapping group shrinkage (OGS) algorithm was initially proposed by Chen and Selesnick [33] to address signal denoising when large-amplitude coefficients form groups. Chen and Selesnick [34] further introduced non-convex regularization terms to make the total cost function convex, enhancing its performance under sparsity and refining the OGS algorithm. Unlike traditional sparse denoising methods, which treat each coefficient independently, OGS considers groups of coefficients. This grouping helps to preserve the relationships between signal components, which is crucial for accurately recovering the underlying signal structure. He et al. [35] integrated a periodicity-induced OGS (POGS) algorithm into the penalty function to mimic the periodicity of sparse groups for fault detection in rotating machinery. Li et al. [36] devised a non-convex penalty function based on an elastic net and the Lp norm, proposing a sparsity-enhanced POGS method for detecting weak periodic impulses. Moreover, although order-statistic-based detection methods pose challenges related to their computational complexity and sensitivity to parameter selection, they offer robust performance in noisy and cluttered environments with minimal prior information [17,30].
The main contributions of this article are summarized as follows:
(1)
An interval-adaptive OGS-based denoising algorithm is developed for Doppler-spread line spectra with equidistant sparsity, which performs well in restoring the Doppler cells occupied by line spectra at low SNRs;
(2)
A GLRT based on the order statistics of a denoised Doppler sequence is proposed for target detection, and an analytical expression for the false alarm probability with any combination in the index set is derived.
The remaining sections of this paper are organized as follows: Section 2 presents the signal model. In Section 3, the OGS denoising algorithm is introduced, the denoising process of line spectra with equidistant sparse constraints is deduced, and interval-adaptive steps (as well as the case where the TOA undergoes a non-integer cycle) are provided. Section 4 derives the generalized likelihood ratio detector guided by the order-denoised Doppler signals and provides an analytical expression for the false alarm probability. Section 5 comprehensively evaluates the proposed detector’s performance through Monte Carlo simulations and comparisons with several typical GLRT-based detectors. Finally, Section 6 summarizes the paper.

2. Signal Model

It is assumed that the target observation attitude (TOA) exhibits periodicity over a coherent processing interval (CPI) T I with period T σ satisfying T r T σ T I , where T r represents the pulse repetition time (PRT). For such targets, a common scattering characteristic model is the sinusoidal frequency-modulated signal [37]. Then, the temporal variation in the target’s complex scattering factor can be decomposed into several harmonic components. To develop a model suitable for more general periodic signals, this paper relaxes the constraints on the specific positions and amplitudes of the harmonics and retains only the constraint of equidistance [38]. Consequently, the discrete form of the complex scattering factor of the space target in the slow-time dimension can be expressed as
x n = l = 0 L 1 A l exp j 2 π l 0 + l n T r 2 π l 0 + l n T r T σ T σ + φ l , n = 1 , , N ,
where N is the number of pulses; L represents the harmonic number, which also corresponds to the number of line spectra; A l and φ l denote the amplitude and initial phase of the l-th harmonic, respectively; and l 0 signifies the starting position of the line spectra. Therefore, when the range and Doppler migration caused by translational motion are precisely compensated, the Doppler spectrum of the multi-pulse coherent integration exhibits sparse equidistant line spectra. The spectrum of the complex scattering factor during a CPI is
X f = T I T σ l = l 0 l 0 + L 1 X σ l T σ s i n c T I f ,
where ⊗ denotes a convolution operator; X σ f represents the spectrum corresponding to the complex scattering factor of a single TOA period; X f can be viewed as sampling X σ f at intervals of 1 1 T σ T σ with a resolution of 1 1 T I T I ; and its discrete form manifests as L equidistant line spectra with an interval of I = T I T I T σ T σ . For convenience, T I T I T σ T σ is first assumed to be an integer in this context, while the case where T I T I T σ T σ is a non-integer is discussed in detail in Section 3.4. Figure 1 depicts the normalized echo signals in the slow-time and Doppler dimensions at an SNR of −12 dB with eight line spectra. In this study, the SNR is defined as the ratio of the total energy of all line spectra to the average noise energy before coherent integration, denoted as
S N R = 10 log 1 σ 2 l = 0 L 1 A d l 2 ,
where σ 2 signifies the average noise energy in the slow-time domain.
The following research is carried out in the Doppler domain, in which the Doppler noisy observation vector is represented as
Y = X + W ,
where Y = y 1 , , y N T C N denotes the Doppler noisy observation vector, X = x 1 , , x N T C N represents the equidistant sparse line spectra in the Doppler domain, and W = w 1 , , w N T C N is the complex Gaussian white noise, the variance of which, denoted as σ w 2 , is assumed to be estimated from reference Doppler cells without target echoes.

3. Denoising Based on Overlapping Group Shrinkage

In this section, the fundamental theory of the overlapping group shrinkage (OGS) algorithm is briefly introduced. Subsequently, the denoising process for a complex sequence based on OGS constrained by equidistant sparsity and a solution to the convex problem are presented. Then, an adaptive interval adjustment strategy is proposed for the situation where the intervals between the line spectra are unknown. Finally, we discuss the case where the TOA within a CPI presents a non-integer number of cycles.

3.1. Overlapping Group Shrinkage

The OGS algorithm can preserve the correlation between signals, effectively inducing a specific sparse structure in them. Its objective function is typically formulated as a convex function and is solved using the majorize–minimization (MM) iterative algorithm, the core principle of which is successive upper-bound minimization. It involves devising a series of approximate optimization functions to bound the original function from above and then iteratively minimizing this sequence to converge toward the optimal solution of the original objective function. This problem can be described as estimating a sparse signal x from a noisy observation vector y as follows:
x = arg min x R N J 0 x = 1 2 y x 2 2 + λ n = 1 N ϕ k = 1 K x n + k 2 1 / 2 ; a ,
where J 0 represents the objective function, K is the size of a group, λ is a regularization parameter, and ϕ : R R + denotes the penalty function, with the parameter a inducing corresponding sparsity in the optimization problem, which is subject to the conditions listed in Appendix A.
Common penalty functions used to promote sparsity include the logarithmic penalty, arctangent penalty, and first-order rational penalty. When the parameter of the penalty function approaches zero, the aforementioned penalty functions approximate the absolute value function, making the absolute value function a type of penalty function that satisfies the specified conditions. With the same value of a, the arctangent penalty function exhibits the most concavity, with non-convex regularization being more effective in accurately recovering sparse signals than convex regularization. Therefore, this study proceeds with subsequent research based on the arctangent penalty, as presented in Equation (6). Detailed analyses focused on parameter selection for a and λ have been presented in [35,39].
ϕ x ; a = 2 3 a arctan 1 + 2 a x 3 π 6 ,

3.2. Denoising of Line Spectra with Equidistant Sparsity

To induce equidistant sparsity in the Doppler sequence, a 0 , 1 binary-weighted sequence is introduced into the OGS model:
x = arg min x R N J 1 x = 1 2 y x 2 2 + λ n = 1 N ϕ k = 1 K b k x n + k 2 1 / 2 ; a ,
b = 1 , 0 , , 0 I , 1 , 0 , , 0 I , , 1 , 0 , , 0 I K = M I
where b k s refer to elements of b , b is a binary-weighted sequence, and M is the number of intervals contained in one group. Following this, for the convenient direct application of the aforementioned reality-based theoretical derivations, the complex Doppler signal is represented as an extended signal stacked with its real and imaginary parts. Equation (4) can then be reformulated as
Y ¯ = X ¯ + W ¯ ,
where
Y ¯ = Δ Re Y T Im Y T T = y ¯ 1 , y ¯ 2 , , y ¯ 2 N T R 2 N ;
X ¯ and W ¯ have the same definitions as above. The extended signal X ¯ is a stack of the real and imaginary parts of the complex domain signal X ; as such, the first and second halves of X ¯ possess the same sparse structure. Therefore, we constrain the two parts of the extended signal with the same weighted sequence and and substitute Equation (9) into Equation (7), yielding
X ¯ = arg min X ¯ R 2 N J ( X ¯ ) = 1 2 Y ¯ X ¯ 2 2 + λ n = 1 N ϕ k = 0 K 1 b k x ¯ n + k 2 + x ¯ N + n + k 2 1 / 2 ; a .
J is strictly convex if 0 a 1 λ M , which is proven in detail in Appendix B. The detailed derivation for solving problem (11) using the MM method is as follows:
u i + 1 = arg min u H u , u i ,
where i denotes the iteration index and H : R 2 N × R 2 N R is a majorizer of J , which satisfies
H u , v J u H v , v = J v , u R N .
A majorizer φ : R × R R of the penalty function ϕ was given in [34]:
φ ( u , v ) = u 2 v 2 2 v ϕ v + ϕ v ,
Hence, we define H as
H X ¯ , X ¯ i = 1 2 Y ¯ X ¯ 2 2 + λ n = 1 N φ ξ n , ζ n .
In Equation (15), ξ n = k = 0 K 1 b k x ¯ n + k 2 + x ¯ N + n + k 2 1 / 2 and ζ n = k = 0 K 1 b k x ¯ n + k i 2 + x ¯ N + n + k i 2 1 / 2 represent the l 2 norms of the weighted group for the current and previous iterations, respectively. Therefore, when the binary weight sequence b is given, ζ n is only related to the previous iteration value X ¯ i . Substituting Equation (14) into Equation (15) and combining the terms unrelated to X ¯ yields
H X ¯ , X ¯ i = 1 2 Y ¯ X ¯ 2 2 + 1 2 λ n = 1 N ξ n 2 ζ n 2 ζ n ϕ ζ n + ϕ ζ n = 1 2 Y ¯ X ¯ 2 2 + 1 2 λ n = 1 N ϕ ζ n ζ n ξ n 2 + C X ¯ i .
Here, C X ¯ i represents a constant that is not dependent on X ¯ . Further transformations applied to Equation (16) yield
H X ¯ , X ¯ i = 1 2 Y ¯ X ¯ 2 2 + 1 2 λ n = 1 N x ¯ n 2 + x ¯ N + n 2 r n X ¯ i + C X ¯ i ,
where
r n X ¯ i = j = 0 K 1 b j ϕ ζ n j ζ n j .
Expanding Y ¯ and X ¯ separately and combining the terms unrelated to X ¯ , Equation (17) is rewritten as
H X ¯ , X ¯ i = n = 1 N 1 + λ r n X ¯ i 2 x ¯ n 2 + x ¯ N + n 2 y ¯ n x ¯ n y ¯ N + n x ¯ N + n + C X ¯ i , Y ¯ .
Similarly, C X ¯ i , Y ¯ is a constant only depending on X ¯ i and Y ¯ . Equation (19) provides an explicit minimizer. Finally, the MM iteration with respect to Equation (12) is given by
x ¯ n i + 1 = y ¯ n 1 + λ r n X ¯ i , 1 n N y ¯ n 1 + λ r n N X ¯ i ,    N + 1 n 2 N .

3.3. Interval-Adaptive Steps

In practice, the interval between line spectra is generally unknown. Therefore, we propose a traversal estimation method based on maximizing the sum of the largest l 2 norms of weighted groups, allowing us to update the interval in each optimization iteration. Specifically, the l 2 norms of the weighted groups are first sorted in ascending order:
ζ 1 ζ n ζ N ,
Then, the actual interval is subject to maximization of the sum of the h largest terms of ζ n ; that is,
I = arg max I I min , I max , I N n = N h + 1 N ζ n ,
where I min , I max is the prior information regarding the intervals between line spectra.
The complete interval-adaptive denoising algorithm based on OGS for line spectra with equidistant sparsity is summarized in Algorithm 1. Here, X denotes the denoised signal, which exhibits an equidistant sparsity in the Doppler domain. Moreover, Figure 2 illustrates a Doppler sequence with S N R = 20 dB, L = 8 , and I = 16 before and after denoising. The local zoom in Figure 2a reveals that only two main line spectra are visible before denoising, while Figure 2b shows the accurate recovery of all Doppler cells occupied by the line spectra, thus demonstrating the ability of the denoising algorithm at low SNRs. Therefore, from Figure 2c, it is easy to conclude that the interval-adaptive denoising algorithm rapidly converges to the actual interval.
Algorithm 1 OGS-based interval-adaptive denoising algorithm
Input: The noisy observation vector Y ; regularization parameter λ ;
   the penalty function ϕ ; prior information of the interval I min , I max ;
Pre-processing with Equation (10);
Initialization: X ¯ 0 = Y ¯ ;
Procedure:
For each i 0 , i t e r do
For each I = I min , I max , I N do
  Compute b with Equation (8);
   ζ n = k = 0 K 1 b k x ¯ n + k i 2 + x ¯ N + n + k i 2 1 / 2 ;
  Sort ζ n in the ascending order;
  Estimate interval utilizing Equation (22);
End for
 Update weight sequence with Estimated interval;
 Compute r n X ¯ i with Equation (18);
 Iterate x ¯ n i + 1 using Equation (20);
End for
Output: X = X ¯ 1 : N + j X ¯ N + 1 : 2 N .

3.4. Non-Integer Period Case

It is often not feasible to guarantee that a CPI contains integer cycles of TOA variations. In this subsection, the case of non-integer cycles of TOA variations within a CPI is considered.
Let I = T I I = T I T σ T σ , where I is a non-integer value. In this case, the line spectra are no longer completely equidistantly sparse in the Doppler domain, with the interval between the l-th and l + 1 -th line spectra being
I l = r o u n d l I r o u n d l 1 I , l = 1 , , L 1
that is, I l I , I . Here, r o u n d x denotes the integer closest to x; x and x represent the integer flooring and ceiling of x, respectively.
Nevertheless, the interval between the line spectra oscillates between the two values, resembling a pseudo-equidistant sparsity. According to Section 3.2, the binary weight sequence—which utilizes a unit-step sliding window to capture the sparse structure of the global line spectra—does not necessarily provide a perfect match with the line spectra. To enhance the sparsity of the line spectra, we increase the count of element 1 in the binary weight sequence as follows:
b = 1 , 1 , 0 , , 0 I , 1 , 1 , 0 , , 0 I , , 1 , 1 , 0 , , 0 I K = M I ,
thereby enabling better adaptation to subtle variations in the interval between the line spectra during the sliding window process. Figure 3 illustrates a Doppler sequence with S N R = 20 dB, L = 8 , and I = 17 , 17 , 16 , 17 , 16 , 16 , 16 before and after denoising. The graph shows that the OGS denoising algorithm can also effectively reconstruct the Doppler cells of the line spectra when the interval oscillates between the two adjacent integers at low SNRs.

4. OGSos-GLRT Detector

In this section, a generalized likelihood ratio test detector based on overlapping group shrinkage denoising and order statistics (OGSos-GLRT) is presented, as well as its schematic. Furthermore, the probability of a false alarm for a given index set of Doppler cells integrated in the detector is derived.

4.1. The Proposed Detector

The problem of detecting line spectra with equidistant sparsity in the Doppler domain can be formulated as a binary hypothesis:
H 0 : y n = w n H 1 : y n = x n + w n , n = 1 , . . . , N ,
The joint probability density functions (pdfs) under H 0 / H 1 are, respectively, denoted as
P Y H 0 Y H 0 = 1 π N σ w 2 N exp Y H Y σ w 2 ,
P Y H 1 Y X , H 1 = 1 π N σ w 2 N exp Y X H Y X σ w 2 ,
where the superscript H represents the complex conjugated transpose.
In order to avoid collapsing loss—which refers to processing in which additional noise samples are applied to the integrator—specific Doppler cells are integrated in the likelihood ratio detector. Based on the OGS denoising algorithm introduced in the previous section, the Doppler sequence is denoised to a series of line spectra, and the Doppler cells occupied by the most likely target harmonics are restored. Then, the initial Doppler cells corresponding to some of the largest denoised Doppler observations are selected for the likelihood ratio detector. The generalized likelihood ratio can be expressed as
Λ O G S o s G L R T = P Y H 1 Y X ^ , H 1 P Y H 0 Y H 0 = n G 1 π σ w 2 exp y n 2 σ w 2 n G 1 π σ w 2 exp y n x n 2 σ w 2 = n G y n 2 < > H 0 H 1 γ ,
where γ is the detection threshold and G is the index set of Doppler cells integrated in the detector.

4.2. False Alarm Probability

This subsection outlines the process of deducing the formula that correlates the false alarm probability with the detection threshold for OGSos-GLRT. The denoising algorithm discussed above diversifies the cases of the composition of set G in the Doppler cells. Here, the false alarm probability for any G set can be derived based on the order statistics. Initially, let z n = y n σ w σ w 2 2 2 , n = 1 , , N , where z n denotes the energy of the Doppler sequence. Then, z n follows an exponential distribution with a pdf:
P z z = 1 2 exp z 2 .
The energy order statistics of the noisy observations can be described by ordering z n s as
0 z 1 z n z N ,
and the z n no longer adheres to independent and identically distributed (IID) assumptions. However, δ z n = N n + 1 z n + 1 z n conforms to IID, having a pdf consistent with that of z n [40]. Let v n = δ z n δ z n c n c n , n = 1 , , N and c n = N n + 1 N n + 1 N n , S N n , S , where S denotes the indices of z n : n G in z n and N n , S represents the count of elements in S that are greater than or equal to n; that is, z n : n G in z n and z n : n S in z n belong to the same set of variables and correspond in a one-to-one manner. Then, we have n S z n = n = 1 N v n , and the test statistics of OGSos-GLRT can be rewritten as
Λ O G S o s G L R T = n = 1 N v n < > H 0 H 1 γ ,
where γ is the new detection threshold. As previously stated, v n s are independent but not identically distributed, with their pdfs and moment generation functions (MGFs) being pairs of Laplace transforms, respectively, denoted as [30]
P v n v n H 0 = c n 2 exp c n v n 2 ,
M v n s = E exp v n s = c n 2 s + c n .
Hence, the MGF of the test statistics of OGSos-GLRT can be denoted as the product of M v n s ; that is,
M Λ O G S o s G L R T s = n = 1 N M v n s = n = 1 N c n 2 s + c n .
Furthermore, the pdf of Λ O G S o s G L R T under H 0 is the inverse Laplacian transformation of its MGF, which can be expressed as
P Λ O G S o s G L R T τ | H 0 = L 1 n = 1 N c n c n 2 2 s + c n c n 2 2 = L 1 i = 1 k j = 1 m i a i , j s + c i c i 2 2 j = i = 1 k j = 1 m i a i , j τ j 1 j 1 ! exp c i τ 2 ,
where L 1 denotes the inverse Laplacian operator; 1 k N is the number of poles; m i represents the multiplicity of c i c i 2 2 , satisfying 1 m i N and m 1 + + m k = N . a i , j denotes the coefficients of the rational fraction decomposition, which are expressed as follows:
a i , j = n = 1 N c n 2 · 1 m i j ! · d m i j d s m i j i = 1 k 1 s + c i c i 2 2 m i · s + c i c i 2 2 m i s = c i c i 2 2 .
Overall, the false alarm probability of OGSos-GLRT with the detection threshold of γ is given by Equation (37):
P f a = P Λ O G S o s G L R T γ H 0 = γ i = 1 k j = 1 m i a i , j τ j 1 j 1 ! exp c i τ 2 d τ = i = 1 k j = 1 m i r = 0 j 1 exp c i 2 γ 2 c i r + 1 a i , j γ j 1 r j 1 r ! .
Finally, a schematic of the proposed detection algorithm is shown in Figure 4. Figure 5a illustrates the probabilities of a false alarm versus the detection threshold for random G , obtained through both the analytical expression and Monte Carlo simulations. The result confirms the validity of Equation (37). In addition, the CFAR characteristic of the detection method is also examined by Monte Carlo simulations under different noise power conditions. Figure 5b illustrates the curves of the false alarm probability versus the noise standard deviation with a preset normalized detection threshold, and demonstrates the CFAR characteristic of the algorithm.

5. Simulations

In this section, Monte Carlo simulations are carried out for MF-GLRT, OS-GLRT, and NLSD-GLRT in addition to OGSos-GLRT in order to comparatively evaluate the performance of the proposed detector. MF-GLRT represents a GLRT, with the actual Doppler cells occupied by the line spectra integrated in the test statistics; it is an ideal likelihood ratio detector, which is generally difficult to achieve in practice. Furthermore, OS-GLRT represents the energy order statistics GLRT [30], while NLSD-GLRT represents a GLRT that is not dependent on the line spectra characteristics [29]. The test statistics for the three mentioned detectors are defined as
Λ M F G L R T = n K z n ,
Λ O S G L R T = n = N L ^ + 1 N z n ,
Λ N L S D G L R T = n = 1 N z n .
where K is the index set of Doppler cells occupied by the line spectra.

5.1. Simulation Parameters

Except for cases in which specific parameters are investigated for their impact on detection performance and explicitly stated, the simulation parameters were set as shown in Table 1. In the table, the number of pulses, T r , and T σ were all set based on parameters commonly used in real-world space target surveillance scenes. As Equation (2), these parameters determine the interval of the line spectra, which is the ratio of the CPI to the TOA period. The number of line spectra is determined by the electrical size of the space target (the ratio of the target size to the signal wavelength), of which, the impact on the proposed algorithm is analyzed later. In radar systems, a higher false alarm rate may undermine the stability and efficiency of the entire system while a lower false alarm rate might result in the missed detection of actual targets. In practical applications, the false alarm probability is typically set between 10 6 and 10 3 . As a rule of thumb, the minimum number of Monte Carlo simulations required for radar signal processing is 10 / P f a , while a preferable number would be 100 / P f a . However, for a compromise between accuracy and computational cost, P f a and the number of Monte Carlo simulations are set to 10 4 and 10 5 , respectively. For space target modeling, several typical line spectra models were considered, as presented in Table 2. Model 1 represents uniformly distributed harmonic energy, Model 2 represents uniformly distributed harmonic energy with certain harmonics missing, Model 3 represents the target Doppler spectrum spread as a sinc-shaped envelope with equidistant sparsity, and Model 4 represents the target Doppler spectrum as containing multiple fluctuating envelopes with equidistant sparsity.

5.2. Effect of Regularization Parameter

For brevity in subsequent discussions, we refer to the case where the estimated number of line spectra (ENLS) matches the actual one as the match case (i.e., L = L ^ ); otherwise, it is referred to as the mismatch case. First, we evaluated the effect of the regularization parameter λ in the OGS denoising algorithm on the performance of OGSos-GLRT. Simulations in the match case were conducted with λ ranging from 45 to 65 in increments of 5, and other parameters set as shown in Table 1 were conducted. Figure 6 shows the probability of detection (Pd) curves of OGSos-GLRT versus the SNR. It can be observed that, within a certain range, the regularization parameter λ does not significantly affect the detection performance of OGSos-GLRT. This can be attributed to the denoising algorithm for equidistant sparse line spectra, which effectively restores the indices of the line spectra within a relatively large range of λ . Therefore, without loss of generality, λ = 55 was adopted in the following simulations.

5.3. Detection Performance for Target Models

Next, the impact of different line spectra models on the detection performance of OGSos-GLRT and OS-GLRT was examined, with this factor evidently not affecting MF-GLRT or NLSD-GLRT. Models 1–4 in Table 2 were subjected to simulations in the match case, using OGSos-GLRT and OS-GLRT with the parameter settings shown in Table 1. For comparison, the simulations also involved MF-GLRT and NLSD-GLRT for Model 1. Figure 7 compares the detection performance of the four methods on different models, which reveals that, across the various models, the detection performance of OGSos-GLRT remained relatively consistent, whereas OS-GLRT showed a certain degree of performance decline, particularly with Model 1. This indicates that OGSos-GLRT exhibits greater robustness than OS-GLRT across different line spectra models.

5.4. Effect of Collapsing Loss

Additionally, due to the presence of collapsing loss, the detection number of Doppler cells evidently affected the performance of OGSos-GLRT, OS-GLRT, and NLSD-GLRT, while MF-GLRT remained unaffected. Figure 8 shows the detection curves of the three algorithms in the match case for the four models, with Doppler cells N doubling sequentially from 320. The results indicate that, as the number of Doppler cells increased, the collapsing loss also increased, leading to a decline in the detection performance of OGSos-GLRT, OS-GLRT, and NLSD-GLRT. Notably, OGSos-GLRT and OS-GLRT mitigated the collapsing loss to varying degrees, with OGSos-GLRT demonstrating superior performance in alleviating the collapsing loss. In addition, by comparing the subfigures in Figure 8, the robustness of OGSos-GLRT to models is again demonstrated, which is consistent with Figure 7.
To quantify the collapsing loss, we define the additional SNR required relative to MF-GLRT as the collapsing loss for OGSos-GLRT, OS-GLRT, and NLSD-GLRT in the match case, using the given false alarm probability and detection probability. Figure 9 presents the collapsing losses for OGSos-GLRT, OS-GLRT, and NLSD-GLRT with Model 3 at P f a = 10 4 and P d = 95 % , of which, the transverse axis represents the number of Doppler cells and is normalized to multiples of 320. The figure shows that the collapsing losses of OGSos-GLRT, OS-GLRT, and NLSD-GLRT linearly increased with the number of Doppler cells multiplying, with these increases being 0.2 dB, 0.3 dB, and 1.4 dB, respectively.

5.5. Effect of TOA Period

Similarly, Figure 10 illustrates the Pds of OGSos-GLRT for Model 3 in the match case, with the number of Doppler cells consistent with Figure 8. However, we adjusted the TOA periods to 20, 10, 5, and 2.5 seconds in this simulation. As the Doppler cells increased and the TOA period decreased, the equidistant line spectra became more sparsely distributed in the Doppler domain. Figure 10 indicates that, as the sparsity of the line spectra was enhanced, the detection performance of OGSos-GLRT did not degrade due to the collapsing loss. However, the result for OS-GLRT remained consistent with that depicted in Figure 8, still suffering from the collapsing loss. This suggests that OGSos-GLRT not only exhibits better robustness with respect to targets having different TOA periods but also demonstrates superior detection performance for Doppler line spectra with higher sparsity.

5.6. Detection Performance for Numbers of Line Spectra

Furthermore, the detection performance of OGSos-GLRT, OS-GLRT, and MF-GLRT for Model 1 and with different numbers of line spectra in the match case was also investigated. Figure 11 presents the Pds versus the SNR of the above three detectors, with the numbers of line spectra being L = 4 , 8 , 16 , 32 . The results indicate that, as the number of line spectra increased, the detection performance of MF-GLRT, OGSos-GLRT, and OS-GLRT decreased by varying degrees. This can be attributed to the following: for a given SNR, a larger number of line spectra means that each line spectrum has less energy and that more noise energy is introduced. Two extreme cases are as follows: (1) L = 1 —in this case, the target Doppler spectrum does not spread and coherent integration suffers no loss, resulting in the highest detection probability; and (2) L = N —in this case, the target fully spreads its energy evenly across each Doppler cell, leading to the worst detection performance. Nonetheless, the detection performance of OGSos-GLRT surpassed that of OS-GLRT and approached that of MF-GLRT as the number of line spectra increased.

5.7. Detection Performance in Mismatch Cases

In practice, the number of line spectra of space targets is often unknown. Therefore, the ENLS is used instead of the actual number when utilizing OGSos-GLRT or OS-GLRT for target detection, which might cause mismatch loss. To investigate the detection performance of OGSos-GLRT in the mismatch case, simulations for Model 1 using the OGSos-GLRT and OS-GLRT algorithms with L ^ = 2 , 4 , 8 , 12 , 16 , 32 were carried out, and the Pds are shown in Figure 12. It can be observed that the detection performance of both OGSos-GLRT and OS-GLRT decreased in the mismatch case, especially when the ENLS was significantly smaller than the actual value. The subplot in Figure 12 quantitatively illustrates the mismatch loss of OGSos-GLRT and OS-GLRT for a given detection probability of 0.95 , indicating that the mismatch loss of OGSos-GLRT was about 2.5 dB when L ^ = 2 . Despite the superior robustness of OS-GLRT with respect to variation in the ENLS, the overall detection performance of OGSos-GLRT remained better than that of OS-GLRT.

5.8. Modified OGSos-GLRT

To enhance the robustness of OGSos-GLRT in the mismatch case, we iteratively cycled through the ENLS (from 1 to a predefined maximum value L ^ m ) during the detection process. Each cycle involved the generation of test statistics Λ L ^ and detection thresholds γ L ^ , followed by decision making. If Λ L ^ γ L ^ , then we decided H 1 and terminated the loop; otherwise, the next iteration began with L ^ + 1 . Iteration was continued until L ^ = L ^ m and, if Λ L ^ remained less than γ L ^ at this point, then we decided H 0 . The aforementioned detection algorithm is referred to as the modified OGSos-GLRT, the schematic of which is outlined in Figure 13. Figure 14 illustrates the Pds of OGSos-GLRT in the match case and the modified OGSos-GLRT for Model 1 with varying numbers of line spectra. The results indicate that the detection performance of the modified OGSos-GLRT closely aligned with that of OGSos-GLRT, with some slight advantages observed at low SNRs. This is attributed to the fact that, under a low SNR, the detection performance may deteriorate when Doppler cells containing line spectra are completely masked by noise or cells containing only noise are integrated into the test statistics. However, the modified OGSos-GLRT partially mitigates such loss through leveraging the minimal ENLS for detection.

6. Conclusions

In this study, we presented a robust generalized likelihood ratio test based on overlapping group shrinkage denoising and order statistics (OGSos-GLRT). This method addresses the problem related to the detection of a Doppler spectrum spreading into equidistant sparse line spectra resulting from periodic changes in TOA during the prolonged coherent integration of space targets. To elaborate, we initially utilize the orbital information of space targets to perform the translational compensation and coherent integration of target echoes, yielding corresponding Doppler sequences. Subsequently, leveraging the equidistant sparse characteristics of the Doppler sequences, we employ the OGS algorithm for denoising. Following denoising, the sequences undergo energy sorting, identifying some of the largest values that correspond to the Doppler cell indices before denoising, which are then integrated into the GLRT for the detection of targets. Furthermore, in order to enhance the robustness of the proposed algorithm in terms of the ENLS, we sequentially evaluate the test statistics against the detection threshold using an increasing ENLS until a decision is reached.
We evaluated the performance of the proposed algorithm through Monte Carlo simulations. First, the detection performance of OGSos-GLRT was found to be superior to that of OS-GLRT and NLSD-GLRT and approached that of MF-GLRT, especially for targets with a large number of line spectra. Second, OGSos-GLRT exhibited robustness with respect to the regularization parameter in the OGS denoising algorithm—an advantage that also applies to space targets with different line spectra models and TOA periods. Furthermore, compared with the OS-GLRT and NLSD-GLRT detectors, OGSos-GLRT was the most effective in reducing the collapsing loss due to the presence of Doppler cells that only contain noise. In addition, OGSos-GLRT performed well for line spectra with more sparsity. Despite showing a higher performance loss than OS-GLRT when the ENLS mismatches the actual value, the overall detection performance of OGSos-GLRT was still superior to that of OS-GLRT, and the mismatch loss was effectively addressed by the modified OGSos-GLRT.
This study considers the case of precise compensation for target translational effects. Even if there is a first-order residual translational phase in the target echoes, the line spectra only undergo a Doppler shift in the frequency domain. Notably, the OGS denoising algorithm employs a sliding window across the entire sequence and, therefore, its denoising performance is translation-invariant; that is, it is dependent on the positions of line spectra and does not affect the detection performance. However, when higher-order phase residues are present, the equidistant sparsity of the Doppler spectrum may be disrupted, which will lead to a decline in detection performance. Future research will focus on situations involving residual higher-order phases resulting from translational motions.

Author Contributions

Conceptualization, L.B. and T.F.; methodology, L.B. and T.F.; software, L.B.; validation, L.B., J.H. and S.Z.; formal analysis, L.B. and J.H.; investigation, L.B. and S.Z.; resources, L.B. and D.C.; data curation, L.B. and H.C.; writing—original draft preparation, L.B., H.C. and D.C.; writing—review and editing, T.F. and S.Z.; visualization, L.B. and J.H.; supervision, T.F. and D.C.; project administration, L.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
TOATarget observation attitude
CPICoherent processing interval
SNRSignal-to-noise ratio
GLRTGeneralized likelihood ratio test
ENLSEstimated number of line spectra
OGSOverlapping group shrinkage
SSASpace situational awareness
RMRange migration
DFMDoppler frequency migration
TFDTime–frequency distribution
HRRPHigh-resolution range profile
CFARConstant false alarm rate
SDD-GLRTScatterer density-dependent GLRT
OS-GLRTOrder statistics GLRT
ASCE-GLRTAdaptive scatterers estimation GLRT
ASCE-GLRTAdaptive Doppler steering matrix estimation GLRT
OGSos-GLRTOverlapping group shrinkage and order statistics GLRT
MF-GLRTMatch filter GLRT
NLSD-GLRTNon-line spectra-dependent GLRT

Appendix A

In this appendix, the following assumptions on the penalty function are made:
  • ϕ is continuous on R ;
  • ϕ is even;
  • ϕ has a unit slope at zero;
  • ϕ has a second derivative with a lower bound of a on R 0 ;
  • ϕ increases monotonically and is concave on the positive x-axis;
  • ϕ equals the absolute value function when a = 0 .

Appendix B

In this appendix, the derivation of the condition under which parameter a makes the objective function a strictly convex function J is given.
J ( X ¯ ) = 1 2 n = 1 N y ¯ n x ¯ n 2 + y ¯ N + n x ¯ N + n 2 + λ n = 1 N ϕ ξ n ; a ,
J ( X ¯ ) = 1 2 n = 1 N x ¯ n 2 + x ¯ N + n 2 + R X ¯ + λ n = 1 N ϕ ξ n ; a ,
where R X ¯ is linear with X ¯ . For a strictly convex function, the addition of a linear term preserves its strict convexity. Let the sum of the remaining two terms in J be Q ; that is,
Q ( X ¯ ) = 1 2 n = 1 N x ¯ n 2 + x ¯ N + n 2 + λ n = 1 N ϕ ξ n ; a .
According to k = 0 K 1 b k = M , Q can be rewritten as
Q ( X ¯ ) = 1 2 M n = 1 N k = 0 K 1 b k x ¯ n 2 + x ¯ N + n 2 + λ n = 1 N ϕ ξ n ; a = 1 2 M n = 1 N K + 1 k = 0 K 1 b k x ¯ n + k 2 + x ¯ N + n + k 2 + λ n = 1 N K + 1 ϕ ξ n ; a = n = 1 N K + 1 1 2 M k = 0 K 1 b k x ¯ n + k 2 + x ¯ N + n + k 2 + λ ϕ ξ n ; a = n = 1 N K + 1 W ( x ¯ n ) ,
where W ( x ¯ n ) = 1 2 M k = 0 K 1 b k x ¯ n + k 2 + x ¯ N + n + k 2 + λ ϕ ξ n ; a . Then, Q is strictly convex if W is strictly convex. We rewrite W as:
W ( x ¯ n ) = 1 2 M k k : b k = 1 b k x ¯ n + k 2 + x ¯ N + n + k 2 + λ ϕ k k : b k = 1 b k x ¯ n + k 2 + x ¯ N + n + k 2 1 / 2 ; a ,
and let ξ n = k k : b k = 1 b k x ¯ n + k 2 + x ¯ N + n + k 2 1 / 2 and γ = 1 M . Then, Equation (A5) can be transformed as follows:
F ξ = γ 2 ξ 2 + λ ϕ ξ ; a ,
with its second derivative
F ξ = γ + λ ϕ ξ ; a , ξ 0 .
When ξ 0 , the sufficient condition for F ξ > 0 is ϕ ξ ; a > γ γ λ λ ; when ξ = 0 , F 0 < F 0 + , as stated in Lemma A [34], F is strictly convex. In conclusion, the sufficient condition for F to be strictly convex is 0 < a 1 λ M . When a = 0 , ϕ x ; 0 = x is regarded as an extreme case. Ultimately, when 0 a 1 λ M , F is a strictly convex function; that is, J is strictly convex.

References

  1. Lin, B.; Wang, J.; Wang, H.; Zhong, L.; Yang, X.; Zhang, X. Small Space Target Detection Based on a Convolutional Neural Network and Guidance Information. Aerospace 2023, 10, 426. [Google Scholar] [CrossRef]
  2. Perry, R.; DiPietro, R.; Fante, R. SAR imaging of moving targets. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 188–200. [Google Scholar] [CrossRef]
  3. Kirkland, D. Imaging moving targets using the second-order keystone transform. IET Radar Sonar Navig. 2011, 5, 902–910. [Google Scholar] [CrossRef]
  4. Xu, J.; Yu, J.; Peng, Y.N.; Xia, X.G. Radon-Fourier Transform for Radar Target Detection, I: Generalized Doppler Filter Bank. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 1186–1202. [Google Scholar] [CrossRef]
  5. Kong, L.; Li, X.; Cui, G.; Yi, W.; Yang, Y. Coherent Integration Algorithm for a Maneuvering Target With High-Order Range Migration. IEEE Trans. Signal Process. 2015, 63, 4474–4486. [Google Scholar] [CrossRef]
  6. Li, X.; Cui, G.; Kong, L.; Yi, W. Fast Non-Searching Method for Maneuvering Target Detection and Motion Parameters Estimation. IEEE Trans. Signal Process. 2016, 64, 2232–2244. [Google Scholar] [CrossRef]
  7. Huang, P.; Liao, G.; Yang, Z.; Xia, X.G.; Ma, J.T.; Ma, J. Long-Time Coherent Integration for Weak Maneuvering Target Detection and High-Order Motion Parameter Estimation Based on Keystone Transform. IEEE Trans. Signal Process. 2016, 64, 4013–4026. [Google Scholar] [CrossRef]
  8. Li, X.; Cui, G.; Yi, W.; Kong, L. Sequence-Reversing Transform-Based Coherent Integration for High-Speed Target Detection. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1573–1580. [Google Scholar] [CrossRef]
  9. Huang, P.; Xia, X.G.; Liao, G.; Yang, Z.; Zhang, Y. Long-Time Coherent Integration Algorithm for Radar Maneuvering Weak Target With Acceleration Rate. IEEE Trans. Geosci. Remote Sens. 2019, 57, 3528–3542. [Google Scholar] [CrossRef]
  10. Zhang, J.; Ding, T.; Zhang, L. Longtime Coherent Integration Algorithm for High-Speed Maneuvering Target Detection Using Space-Based Bistatic Radar. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5100216. [Google Scholar] [CrossRef]
  11. Uno, K.; Matsuoka, T.; Uchida, A.; Yoshida, K. Structure from Motion-based Motion Estimation and 3D Reconstruction of Unknown Shaped Space Debris. arXiv 2024, arXiv:2408.01035. [Google Scholar]
  12. Serrano, A.; Kobsa, A.; Uysal, F.; Cerutti-Maori, D.; Ghio, S.; Kintz, A.; Morrison, R.L., Jr.; Welch, S.; van Dorp, P.; Hogan, G.; et al. Long baseline bistatic radar imaging of tumbling space objects for enhancing space domain awareness. IET Radar Sonar Navig. 2024, 18, 598–619. [Google Scholar] [CrossRef]
  13. Zhou, X.; Wang, Y.; Yeh, C.; Lu, X. Precession Parameter Estimation From Wideband Measurements for 3-D ISAR Imaging of Cone-Shaped Targets. IEEE Geosci. Remote Sens. Lett. 2022, 19, 4501305. [Google Scholar] [CrossRef]
  14. Lin, L.; Sun, G.; Cheng, Z.; He, Z. Long-Time Coherent Integration for Maneuvering Target Detection Based on ITRT-MRFT. IEEE Sens. J. 2020, 20, 3718–3731. [Google Scholar] [CrossRef]
  15. Van Der Spek, G.A. Detection of a Distributed Target. IEEE Trans. Aerosp. Electron. Syst. 1971, AES-7, 922–931. [Google Scholar] [CrossRef]
  16. Hughes, P. A High-Resolution Radar Detection Strategy. IEEE Trans. Aerosp. Electron. Syst. 1983, AES-19, 663–667. [Google Scholar] [CrossRef]
  17. Dai, F.; Liu, H.; Wu, S. Order-Statistic-Based Detector for Range Spread Target. J. Electron. Inf. Technol. 2009, 31, 2488–2492. [Google Scholar]
  18. Shui, P.L.; Liu, H.W.; Bao, Z. Range-Spread Target Detection Based on Cross Time-Frequency Distribution Features of Two Adjacent Received Signals. IEEE Trans. Signal Process. 2009, 57, 3733–3745. [Google Scholar] [CrossRef]
  19. Zuo, L.; Li, M.; Zhang, X.; Zhang, P.; Wu, Y. CFAR Detection of Range-Spread Targets Based on the Time-Frequency Decomposition Feature of Two Adjacent Returned Signals. IEEE Trans. Signal Process. 2013, 61, 6307–6319. [Google Scholar] [CrossRef]
  20. Shui, P.L.; Xu, S.W.; Liu, H.W. Range-Spread Target Detection using Consecutive HRRPs. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 647–665. [Google Scholar] [CrossRef]
  21. Yang, X.; Wen, G.; Ma, C.; Hui, B.; Ding, B.; Zhang, Y. CFAR Detection of Moving Range-Spread Target in White Gaussian Noise Using Waveform Contrast. IEEE Geosci. Remote Sens. Lett. 2016, 13, 282–286. [Google Scholar] [CrossRef]
  22. De Maio, A. Polarimetric adaptive detection of range-distributed targets. IEEE Trans. Signal Process. 2002, 50, 2152–2159. [Google Scholar] [CrossRef]
  23. Tang, M.; Rong, Y.; Zhou, J.; Li, X.R. Invariant Adaptive Detection of Range-Spread Targets Under Structured Noise Covariance. IEEE Trans. Signal Process. 2017, 65, 3048–3061. [Google Scholar] [CrossRef]
  24. Sun, S.; Liu, J.; Liu, W.; Jian, T. Robust detection of distributed targets based on Rao test and Wald test. Signal Process. 2021, 180, 107801. [Google Scholar] [CrossRef]
  25. Chen, X.; Gai, J.; Liang, Z.; Liu, Q.; Long, T. Adaptive Double Threshold Detection Method for Range-Spread Targets. IEEE Signal Process. Lett. 2022, 29, 254–258. [Google Scholar] [CrossRef]
  26. Gu, X.; Tao, J.; You, H. Double Threshold CFAR Detector of Range-spreadTarget and Its Performance Analysis. J. Electron. Inf. Technol. 2012, 34, 1318–1323. [Google Scholar] [CrossRef]
  27. Shuai, X.; Kong, L.; Yang, J. Adaptive detection for distributed targets in Gaussian noise with Rao and Wald tests. Sci. China Inf. Sci. 2012, 55, 1290–1300. [Google Scholar] [CrossRef]
  28. Liu, J.; Liu, W.; Tang, B.; Zheng, J.; Xu, S. Distributed Target Detection Exploiting Persymmetry in Gaussian Clutter. IEEE Trans. Signal Process. 2019, 67, 1022–1033. [Google Scholar] [CrossRef]
  29. Gerlach, K. Spatially distributed target detection in non-Gaussian clutter. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 926–934. [Google Scholar] [CrossRef]
  30. He, Y.; Jian, T.; Su, F.; Qu, C.; Gu, X. Novel Range-Spread Target Detectors in Non-Gaussian Clutter. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 1312–1328. [Google Scholar] [CrossRef]
  31. Ren, Z.; Yi, W.; Zhao, W.; Kong, L. Range-Spread Target Detection Based on Adaptive Scattering Centers Estimation. IEEE Trans. Geosci. Remote Sens. 2023, 61, 5100414. [Google Scholar] [CrossRef]
  32. Ren, Z.; Yi, W.; Kong, L.; Farina, A.; Orlando, D. Adaptive Range and Doppler Distributed Target Detection in Non-Gaussian Clutter. IEEE Trans. Signal Process. 2023, 71, 2376–2390. [Google Scholar] [CrossRef]
  33. Chen, P.Y.; Selesnick, I.W. Group-Sparse Signal Denoising: Non-Convex Regularization, Convex Optimization. IEEE Trans. Signal Process. 2014, 62, 3464–3478. [Google Scholar] [CrossRef]
  34. Chen, P.Y.; Selesnick, I.W. Translation-invariant shrinkage/thresholding of group sparse signals. Signal Process. 2014, 94, 476–489. [Google Scholar] [CrossRef]
  35. He, W.; Ding, Y.; Zi, Y.; Selesnick, I.W. Sparsity-based algorithm for detecting faults in rotating machines. Mech. Syst. Signal Process. 2016, 72–73, 46–64. [Google Scholar] [CrossRef]
  36. Li, Z.; Li, J.; Ding, W.; Cheng, X.; Meng, Z. A sparsity-enhanced periodic OGS model for weak feature extraction of rolling bearing faults. Mech. Syst. Signal Process. 2022, 169, 108733. [Google Scholar] [CrossRef]
  37. Wang, Q.; Xing, M.; Lu, G.; Bao, Z. Single Range Matching Filtering for Space Debris Radar Imaging. IEEE Geosci. Remote Sens. Lett. 2007, 4, 576–580. [Google Scholar] [CrossRef]
  38. Bell, M.; Grubbs, R. JEM modeling and measurement for radar target identification. IEEE Trans. Aerosp. Electron. Syst. 1993, 29, 73–87. [Google Scholar] [CrossRef]
  39. Zhao, Z.; Wu, S.; Qiao, B.; Wang, S.; Chen, X. Enhanced Sparse Period-Group Lasso for Bearing Fault Diagnosis. IEEE Trans. Ind. Electron. 2019, 66, 2143–2153. [Google Scholar] [CrossRef]
  40. David, H.A.; Nagaraja, H.N. Order Statistics, 3rd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2003. [Google Scholar]
Figure 1. Normalized echo signals with S N R = 12 dB, L = 8 : (a) slow−time domain. (b) Doppler domain.
Figure 1. Normalized echo signals with S N R = 12 dB, L = 8 : (a) slow−time domain. (b) Doppler domain.
Remotesensing 16 03413 g001
Figure 2. Doppler sequence with S N R = 12 dB, L = 8 , I = 16 : (a) before denoising. (b) After denoising. (c) Interval-updating process.
Figure 2. Doppler sequence with S N R = 12 dB, L = 8 , I = 16 : (a) before denoising. (b) After denoising. (c) Interval-updating process.
Remotesensing 16 03413 g002
Figure 3. Doppler sequence with S N R = 12 dB, L = 8 , I = 17 , 17 , 16 , 17 , 16 , 16 , 16 : (a) before denoising. (b) After denoising.
Figure 3. Doppler sequence with S N R = 12 dB, L = 8 , I = 17 , 17 , 16 , 17 , 16 , 16 , 16 : (a) before denoising. (b) After denoising.
Remotesensing 16 03413 g003
Figure 4. Schematic of the proposed detector.
Figure 4. Schematic of the proposed detector.
Remotesensing 16 03413 g004
Figure 5. Probability of false alarm: (a) comparison between the analytical expression and Monte Carlo simulations. (b) Sensitivity analysis with respect to noise power.
Figure 5. Probability of false alarm: (a) comparison between the analytical expression and Monte Carlo simulations. (b) Sensitivity analysis with respect to noise power.
Remotesensing 16 03413 g005
Figure 6. Pds versus SNR of OGSos−GLRT with Model 1 for λ = 45 , 50 , 55 , 60 , 65 .
Figure 6. Pds versus SNR of OGSos−GLRT with Model 1 for λ = 45 , 50 , 55 , 60 , 65 .
Remotesensing 16 03413 g006
Figure 7. Pds versus SNR of MF−GLRT, OGSos−GLRT, OS−GLRT, and NLSD−GLRT for Models 1−4.
Figure 7. Pds versus SNR of MF−GLRT, OGSos−GLRT, OS−GLRT, and NLSD−GLRT for Models 1−4.
Remotesensing 16 03413 g007
Figure 8. Pds versus SNR of OGSos−GLRT, OS−GLRT, and NLSD−GLRT with the 4 models for N = 320 , 640 , 1280 , 2560 .
Figure 8. Pds versus SNR of OGSos−GLRT, OS−GLRT, and NLSD−GLRT with the 4 models for N = 320 , 640 , 1280 , 2560 .
Remotesensing 16 03413 g008
Figure 9. Collapsing loss versus number of Doppler cells for OGSos−GLRT, OS−GLRT, and NLSD−GLRT with Model 3 at P f a = 10 4 and P d = 95 % .
Figure 9. Collapsing loss versus number of Doppler cells for OGSos−GLRT, OS−GLRT, and NLSD−GLRT with Model 3 at P f a = 10 4 and P d = 95 % .
Remotesensing 16 03413 g009
Figure 10. Pds versus SNR of OGSos−GLRT and OS−GLRT with Model 3 for T σ = 20 , 10 , 5 , 2.5 s.
Figure 10. Pds versus SNR of OGSos−GLRT and OS−GLRT with Model 3 for T σ = 20 , 10 , 5 , 2.5 s.
Remotesensing 16 03413 g010
Figure 11. Pds versus SNR of MF−GLRT, OGSos−GLRT, and OS−GLRT with Model 1 for L = 4 , 8 , 16 , 32 .
Figure 11. Pds versus SNR of MF−GLRT, OGSos−GLRT, and OS−GLRT with Model 1 for L = 4 , 8 , 16 , 32 .
Remotesensing 16 03413 g011
Figure 12. Pds versus SNR of OGSos−GLRT and OS−GLRT with Model 1 for L = 8 , L ^ = 2 , 4 , 8 , 12 , 16 , 32 ; mismatch loss of OGSos−GLRT and OS−GLRT versus L ^ with P d = 95 % .
Figure 12. Pds versus SNR of OGSos−GLRT and OS−GLRT with Model 1 for L = 8 , L ^ = 2 , 4 , 8 , 12 , 16 , 32 ; mismatch loss of OGSos−GLRT and OS−GLRT versus L ^ with P d = 95 % .
Remotesensing 16 03413 g012
Figure 13. Schematic of the modification to OGSos−GLRT.
Figure 13. Schematic of the modification to OGSos−GLRT.
Remotesensing 16 03413 g013
Figure 14. Pds versus SNR of OGSos−GLRT (in the match case) and modified OGSos−GLRT with Model 1 for L = 4 , 8 , 16 , 32 .
Figure 14. Pds versus SNR of OGSos−GLRT (in the match case) and modified OGSos−GLRT with Model 1 for L = 4 , 8 , 16 , 32 .
Remotesensing 16 03413 g014
Table 1. Parameters used in simulations.
Table 1. Parameters used in simulations.
ParametersValues
Regularization parameter λ 55
Number of pulses4000
Pulse repetition time T r 40 ms
TOA period T σ 10 s
Number of line spectra L8
Number of Doppler cells N640
Probability of false alarm P f a 10 4
Number of Monte Carlo simulations 10 5
Table 2. Typical models of Doppler-spread targets.
Table 2. Typical models of Doppler-spread targets.
ModelsDoppler Cells
1 2 3 4 5 6 7 8
111111111
210110011
30.30.60.9110.90.60.3
40.510.50.250.250.510.5
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

Bu, L.; Fu, T.; Chen, D.; Cao, H.; Zhang, S.; Han, J. Doppler-Spread Space Target Detection Based on Overlapping Group Shrinkage and Order Statistics. Remote Sens. 2024, 16, 3413. https://doi.org/10.3390/rs16183413

AMA Style

Bu L, Fu T, Chen D, Cao H, Zhang S, Han J. Doppler-Spread Space Target Detection Based on Overlapping Group Shrinkage and Order Statistics. Remote Sensing. 2024; 16(18):3413. https://doi.org/10.3390/rs16183413

Chicago/Turabian Style

Bu, Linsheng, Tuo Fu, Defeng Chen, Huawei Cao, Shuo Zhang, and Jialiang Han. 2024. "Doppler-Spread Space Target Detection Based on Overlapping Group Shrinkage and Order Statistics" Remote Sensing 16, no. 18: 3413. https://doi.org/10.3390/rs16183413

APA Style

Bu, L., Fu, T., Chen, D., Cao, H., Zhang, S., & Han, J. (2024). Doppler-Spread Space Target Detection Based on Overlapping Group Shrinkage and Order Statistics. Remote Sensing, 16(18), 3413. https://doi.org/10.3390/rs16183413

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