Next Article in Journal
Enhancing Identification of Meteorological and Biological Targets Using the Depolarization Ratio for Weather Radar: A Case Study of FAW Outbreak in Rwanda
Previous Article in Journal
Active Bidirectional Self-Training Network for Cross-Domain Segmentation in Remote-Sensing Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Adaptive Radar Target Detection Method Based on Alternate Estimation in Power Heterogeneous Clutter

by
Daipeng Xiao
,
Weijian Liu
*,
Hui Chen
,
Hao Li
and
Binbin Li
Wuhan Electronic Information Institute, Wuhan 430019, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(13), 2508; https://doi.org/10.3390/rs16132508
Submission received: 20 May 2024 / Revised: 24 June 2024 / Accepted: 5 July 2024 / Published: 8 July 2024
(This article belongs to the Section Engineering Remote Sensing)

Abstract

:
Multichannel radars generally need to utilize a certain amount of training samples to estimate the covariance matrix of clutter for target detection. Due to factors such as severe terrain fluctuations and complex electromagnetic environments, the training samples usually have different statistical characteristics from the data to be detected. One of the most common scenarios is that all data have the same clutter covariance matrix structure, while different data have different power mismatches, called power heterogeneous characteristics. For detection problems in the power heterogeneous clutter environments, we propose detectors based on alternate estimation, using the generalized likelihood ratio test (GLRT) criterion, Rao criterion, Wald criterion, Gradient criterion, and Durbin criterion. Monte Carlo simulation experiments and real data indicate that the detector based on the Rao criterion has the highest probability of detection (PD). Furthermore, when signal mismatch occurs, the detector based on the GLRT criterion has the best selectivity, while the detector based on the Durbin criterion has the most robust detection performance.

1. Introduction

A radar system is an electronic device that employs electromagnetic waves for target detection, which is extensively used in various aspects of life, such as military early warning, wireless communication, navigation, meteorology, and so on. One of the principal roles of radar is to detect targets [1,2,3,4]. Traditional radar target detection can be categorized into two steps: the first step is the filtering process, which involves suppressing clutter to enhance the signal-to-clutter ratio (SCR); the second step is to process the filtered data with a constant false alarm rate (CFAR). This method is complex and costly. In the past few years, adaptive detection has gradually come into people’s sight. Its core idea is to jointly use the data to be detected and training samples to construct an adaptive detector, which integrates filtering and CFAR processing [5]. Compared with the step-by-step detection method, adaptive detection technology not only simplifies the process but also improves detection performance.
It is important to mention that the clutter covariance matrix is typically unknown and requires estimation using training samples [6]. The traditional assumption model holds that the training samples and the data to be detected have the same statistical characteristics, following a complex Gaussian distribution [7], but this assumption often does not hold in complex electromagnetic environments. Due to factors such as clutter fluctuations in real scenarios, environments often show non-homogeneous characteristics, especially non-Gaussian environments. However, training samples should have some similarities to the data to be detected [8]; otherwise, they cannot provide useful information. In 1999, Kraut et al. put forward the generalized likelihood ratio test (GLRT) detector in a partially homogeneous environment, known as the adaptive coherence estimator [9]. In a partially homogeneous environment, the training samples share the same covariance matrix, but an unknown power mismatch exists between the training samples and the data to be detected. In a complicated and changeable real environment, this assumption is also too idealistic.
As the resolution of radar enhances, the quantity of scatterers within each resolution cell diminishes, leading to different levels of local roughness of clutter within different range resolution cells and causing larger amplitude fluctuations. In such cases, clutters exhibit non-Gaussian characteristics [10,11,12]. The distribution of amplitudes for non-Gaussian clutter exhibits a heavier tail, rendering the Gaussian model insufficient in accurately describing such clutter. The compound-Gaussian model [13,14,15] is an extensively accepted model for modeling non-Gaussian clutter, which can be characterized as the product of mutually independent slow-changing texture components and fast-changing speckle components, where the former is a non-negative random variable and the latter is a complex Gaussian random vector. Non-Gaussian clutter can be expressed as various compound-Gaussian models, depending on the specific texture component models. When the texture component follows a gamma distribution, the compound-Gaussian clutter is K distribution clutter. It can describe non-homogeneous clutter well at high resolution and small grazing angles. When the texture component follows an inverse gamma distribution, compound-Gaussian clutter is Pareto distribution clutter. It is more suitable for describing clutter with heavy tails. When the texture component follows an inverse Gaussian distribution, compound-Gaussian clutter is CG-IG distribution clutter. It has a better fitting capability for high-resolution sea clutter. Among them, the most common clutter is K distribution clutter [16,17,18]. However, not all types of real clutter can be modeled as a compound-Gaussian model. Additionally, some models that fit real clutter well do not conform to the compound-Gaussian model, such as Weibull distribution [19] and Log-Normal distribution [20]. Under the compound-Gaussian model, it is often necessary to assign a certain prior statistical distribution to the texture component [21], but real clutter environments may not necessarily correspond to the given distribution. In this case, the detector, designed according to the predetermined statistical distribution of the texture component, is bound to result in performance loss.
To solve the aforementioned problem, it is assumed in [22] that there also exists unknown power mismatches between the training samples, and the amount of mismatches are unknown non-random and non-negative quantities. This model is an important complement to the compound-Gaussian model, which can also be transformed from the compound-Gaussian model. Based on the GLRT criterion, together with the alternate estimation algorithm, Coluccia et al. proposed an adaptive detector in [22], which has the CFAR property and offers better detection performance compared with the well-known adaptive normalized matched filter (ANMF).
It is important to highlight that the detection problem for only rank-one point targets is considered in [22], and the detector design is solely based on the GLRT criterion. However, in practice, the signal might be categorized as a subspace signal, which is an extension of rank-one signals and refers to the signal located in the known subspace with unknown coordinates. Meanwhile, for target detection problems in real scenarios, there is often no optimal detector due to the presence of too many unknown parameters. In addition to the GLRT criterion, other useful detector design criteria also include the Rao criterion, Wald criterion, Gradient criterion, and Durbin criterion, among others. Therefore, in this paper, we extend the signal model in [22], extending the rank-one model to the subspace signal model, and design detectors based on the GLRT criterion, Rao criterion, Wald criterion, Gradient criterion, and Durbin criterion.
In addition, signal mismatch was not considered in [22]. In practical applications, factors such as array element errors, multipath effects, and rapid maneuvers of targets often lead to signal mismatch [23,24,25], meaning that the actual steering vector of the target does not match the value assumed by the radar system. There is a type of detector whose detection performance sharply degrades with signal mismatch increases, known as mismatch-sensitive detectors, which are usually applied in fields such as target tracking [26]. Conversely, another type of detector can still detect targets with a high probability when signal mismatch occurs, known as robust detectors, which are usually applied in fields such as searching models [27]. Therefore, in the performance evaluation, we also consider the performance of proposed detectors when signal mismatch occurs.
In summary, this paper investigates the detection problem for subspace signals in a power heterogeneous clutter environment. Utilizing the GLRT criterion, Rao criterion, Wald criterion, Gradient criterion, and Durbin criterion, in conjunction with the iterative concept of alternating estimation, we propose a variety of detectors, named the alternating estimation GLRT (AE-GLRT) detector, alternating estimation Rao (AE-Rao) detector, alternating estimation Wald (AE-Wald) detector, alternating estimation Gradient (AE-Gradient) detector, and alternating estimation Durbin (AE-Durbin) detector, respectively. Monte Carlo simulation experiments and IPIX radar data both demonstrate that the AE-Rao detector has the highest probability of detection (PD) in the absence of signal mismatch. When signal mismatch occurs, the AE-GLRT exhibits the best selectivity, while the AE-Durbin offers the most robust performance.
Notation: Vectors are represented by bold lowercase letters, whereas matrices are represented by bold uppercase letters. The symbols ( · ) T , ( · ) * , ( · ) H , ( · ) 1 , ( · ) , ( · ) , Re { · } , and | · | express transpose, conjugate, conjugate transpose, inverse, generalized inverse, partial derivative, the real part of a parameter, and the matrix determinant, respectively. ξ ^ denotes the estimation of ξ , E · indicates statistical expectation, and C N × M represents the N × M dimensional complex field.

2. Problem Formulation

According to the compound-Gaussian model, vector c C N × 1 represents the clutter in the cell under test (CUT), which can be denoted as
c = τ x
where τ represents the slowly changing texture component, which is generally a non-negative random variable and x represents the rapidly changing speckle component, following a circular symmetric complex Gaussian distribution with conditional covariance matrix E x x H | τ = τ R , where R C N × N represents the unknown clutter covariance matrix. The training samples can be expressed as c l = τ l x l , l = 1 , 2 , , L , whose conditional covariance matrix is E x l x l H | τ l = τ l R , l = 1 , 2 , , L .
Assume that the radar system possesses N channels. If a signal exists, it can be expressed as N × 1 vector s . For subspace signals, assuming s is located in the subspace H , it can be represented as s = H ξ . Therefore, in the power heterogeneous clutter environment, the detection problem can be represented by the following binary hypothesis test
H 0 : z C N N ( 0 , R ) z l C N N ( 0 , γ l R ) , l = 1 , 2 , , L H 1 : z C N N ( H ξ , R ) z l C N N ( 0 , γ l R ) , l = 1 , 2 , , L
where z C N × 1 is the data to be detected, z l C N × 1 represents the training samples, H C N × p is a known column-full-rank signal matrix, ξ C p × 1 denotes an unknown signal coordinate vector, R C N × N represents the unknown clutter covariance matrix, and γ 1 , γ 2 , , γ l > 0 are unknown parameters that denote the different power mismatch amounts between each training sample and the data to be detected.
As mentioned above, the detection model in (2) can be transformed from the compound-Gaussian model. Specifically, dividing all the clutter data by τ yields the clutter in the CUT as c = x with conditional covariance matrix R , where the value of τ can be estimated by [24]. Meanwhile, the training samples can be rewritten as c l = c l / τ . According to the definition of c l above, we set γ l = τ l / τ , l = 1 , 2 , , L , so the training samples can be written as c l = γ l x l , l = 1 , 2 , , L with conditional covariance matrix being γ l R , l = 1 , 2 , , L . In addition, this paper regards γ l , l = 1 , 2 , , L as definite unknown parameters, which can increase the applicability of this model; that is, it is not limited to a particular compound-Gaussian clutter with a predetermined texture distribution.

3. Detector Design

The joint probability density functions (PDFs) of the data to be detected and the training samples under hypotheses H 1 and H 0 can be expressed as
f 1 ( z , Z ; ξ , R , γ ) = 1 π N ( L + 1 ) l = 1 L γ l N | R | L + 1 exp z H ξ H R 1 z H ξ + l = 1 L z l H R 1 z l γ l
and
f 0 ( z , Z ; R , γ ) = 1 π N ( L + 1 ) l = 1 L γ l N | R | L + 1 exp z H R 1 z + l = 1 L z l H R 1 z l γ l
respectively, where Z = z 1 , z 2 , z L denotes the matrix of training samples and γ = γ 1 , γ 2 , , γ L T represents the column vector of power mismatches between the training samples and the data to be detected. The following subsections will design detectors according to the GLRT criterion, Rao criterion, Wald criterion, Gradient criterion, and Durbin criterion.

3.1. AE-GLRT Detector

The GLRT design criterion can be formulated as [28]
max ξ , R , γ f 1 ( z , Z ; ξ , R , γ ) max R , γ f 0 ( z , Z ; R , γ ) H 1 H 0 σ
where σ denotes the detection threshold set based on the predesigned PFA. Under hypothesis H 1 , taking the partial derivative of (3) with respect to R to be zero, we can obtain the maximum likelihood estimation (MLE) of R for a given γ l as
R ^ 1 = 1 L + 1 z H ξ z H ξ H + l = 1 L z l z l H γ l
Utilizing the same method, we can derive the MLE of R under hypothesis H 0 as
R ^ 0 = 1 L + 1 z z H + l = 1 L z l z l H γ l
Substituting (6) and (7) into (3) and (4), respectively, we can obtain
l 1 ( z , Z ; ξ , R ^ 1 , γ ) = L + 1 e π N ( L + 1 ) l = 1 L γ l N ( z H ξ ) ( z H ξ ) H + l = 1 L z l z l H γ l ( L + 1 )
and
l 0 ( z , Z ; R ^ 0 , γ ) = L + 1 e π N ( L + 1 ) l = 1 L γ l N z z H + l = 1 L z l z l N γ l ( L + 1 )
Directly solving for the MLE of ξ and γ l , l = 1 , 2 , , L does not yield a closed-form solution. Therefore, we use the alternate estimation algorithm to address this problem. For ease of expression, we specify the following symbols: γ ^ l 1 , ( 0 ) and γ ^ l 1 , ( t max ) represent the initial value and final estimation of γ l under hypothesis H 1 , respectively. ξ ^ ( t max ) denotes the final estimation of ξ under hypothesis H 1 ; γ ^ l 0 , ( 0 ) and γ ^ l 0 , ( t max ) represent the initial value and final estimation of γ l under hypothesis H 0 , respectively.
The estimation process is as follows: First, we let the initial value of γ l under hypothesis H 1 be γ ^ l 1 , ( 0 ) , l = 1 , 2 , , L [22], and ξ ^ ( 1 ) can be obtained by (A4); then, we obtain γ ^ l 1 , ( 1 ) , l = 1 , 2 , , L based on (A8). The process is repeated in a loop until the maximum number of iterations t max is reached. In this paper, we set it to 15, whose validity will be shown in the fourth section. Algorithm 1 details the steps of the alternate estimation iterative algorithm.
Algorithm 1 The algorithm of alternate estimation.
Input:  γ ^ l 0 , ( 0 ) , l = 1 , 2 , , L , γ ^ l 1 , ( 0 ) , l = 1 , 2 , , L and t max
Output:  γ ^ l 0 , ( t max ) , l = 1 , 2 , , L , γ ^ l 1 , ( t max ) , l = 1 , 2 , , L and ξ ^ ( t max )
 1:
Let the initial value of γ l under hypothesis H 1 be γ ^ l 1 , ( 0 ) = z l H z z H z l , l = 1 , 2 , , L ;
 2:
Substituting γ ^ l 1 , ( 0 ) , l = 1 , 2 , , L into (A4) yields the value of ξ ^ ( 1 ) ;
 3:
Substituting ξ ^ ( 1 ) into (A8) results in γ ^ 1 1 , ( 1 ) , γ ^ 2 1 , ( 1 ) , , γ ^ L 1 , ( 1 ) ;
 4:
Substituting γ ^ 1 1 , ( 1 ) , γ ^ 2 1 , ( 1 ) , , γ ^ L 1 , ( 1 ) into (A4) results in ξ ^ ( 2 ) ;
 5:
Substituting ξ ^ ( 2 ) into (A8) results in γ ^ 1 1 , ( 2 ) , γ ^ 2 1 , ( 2 ) , , γ ^ L 1 , ( 2 ) . Terminating the process when the number of iterations reaches t max . Finally, we can obtain γ ^ 1 1 , ( t max ) , γ ^ 2 1 , ( t max ) , , γ ^ L 1 , ( t max ) , and ξ ^ ( t max ) ;
 6:
Setting ξ ^ = 0 under hypothesis H 0 , and then repeating the above steps to obtain γ ^ 1 0 , ( t max ) , γ ^ 2 0 , ( t max ) , , γ ^ L 0 , ( t max ) .
When obtaining the estimations of ξ and γ l , l = 1 , 2 , , L from Appendix A, we can substitute (6), (7), (A4), (A8), and (A10) into (5), resulting in the AE-GLRT as
t AE - GLRT = l = 1 L γ ^ l 0 , ( t max ) l = 1 L γ ^ l 1 , ( t max ) N L + 1 z z H + l = 1 L z l z l H γ ^ l 0 , ( t max ) z H ξ ^ ( t max ) z H ξ ^ ( t max ) H + l = 1 L z l z l H γ ^ l 1 , ( t max )

3.2. AE-Rao Detector

The Rao design criterion can be expressed as [29]
t Rao = ln f 1 ( z , Z ) Θ r Θ = Θ ^ 0 T I 1 Θ ^ 0 Θ r Θ r ln f 1 ( z , Z ) Θ r * Θ = Θ ^ 0
where I 1 Θ ^ 0 Θ r Θ r is the value at Θ = Θ ^ 0 of the Θ r , Θ r block of the inverse of the fisher information matrix (FIM) I Θ . It is known by the matrix partitioning lemma
I 1 ( Θ ) Θ r , Θ r = I Θ r , Θ r ( Θ ) I Θ r , Θ s ( Θ ) I Θ s , Θ s 1 ( Θ ) I Θ s , Θ r ( Θ ) 1
For the complex parameter Θ , its FIM is defined as
I Θ = E ln f ( z ; Θ ) Θ * ln f ( z ; Θ ) Θ T
which can be divided into
I ( Θ ) = I Θ r , Θ r ( Θ ) I Θ r , Θ s ( Θ ) I Θ s , Θ r ( Θ ) I Θ s , Θ s ( Θ )
with
Θ = Θ r T , Θ s T T
being a ( N 2 + L + p ) × 1 dimensional parameter vector. Θ r = ξ denotes the relevant paramater and Θ s = γ T , ve c T ( R ) T represents the irrelevant parameter. Θ ^ 0 = Θ ^ r 0 T , Θ ^ s 0 T T indicates the MLE of parameter Θ under hypothesis H 0 with Θ ^ r 0 = 0 p × 1 and Θ ^ s 0 = γ ^ 0 T , ve c T ( R ^ ) T . Θ ^ 1 = Θ ^ r 1 T , Θ ^ s 1 T T indicates the MLE of parameter Θ under hypothesis H 1 with Θ ^ r 1 = ξ ^ and Θ ^ s 1 = γ ^ 1 T , ve c T ( R ^ ) T . According to (3), the partial derivatives of ln f 1 z , Z with respect to ξ and ξ * are
ln f 1 z , Z Θ r = ln f 1 z , Z ξ = z H ξ H R 1 H T
and
ln f 1 z , Z Θ r * = ln f 1 z , Z ξ * = H H R 1 z H ξ
respectively. According to the definition of (14), substituting (16) and (17) into (13) results in
I Θ r , Θ r ( Θ ) = E ln f ( z , Z ) Θ r * ln f ( z , Z ) Θ r T = H H R 1 E z H ξ z H ξ H R 1 H = H H R 1 H
I Θ r , Θ s ( Θ ) can be written as
I Θ r , Θ s ( Θ ) = I ξ , γ ( Θ ) I ξ , vec ( R ) ( Θ )
We can derive the following two equations according to (17)
I ξ , γ ( Θ ) = E 2 ln f 1 ( z , Z ) ξ * γ T = E γ T H H R 1 z H ξ = 0 p × L
and
I ξ , vec ( R ) ( Θ ) = E 2 ln f 1 ( z , Z ) ξ * ve c T ( R ) = E ve c T ( R ) H H R 1 z H ξ = 0 p × N 2
Substituting (18), (19), (20), and (21) into (12) results in
I 1 ( Θ ) Θ r Θ r = H H R 1 H 1
Substituting (16), (17) and (22) into (11) yields the AE-Rao for the given R as
t AE - Rao R = z H R 1 H ( H H R 1 H ) 1 H H R 1 z
Substituting the final alternate estimation of γ 0 into (7) results in the MLE of R under hypothesis H 0 as
R ^ 0 = 1 L + 1 l = 1 L z l z l H γ ^ l 0 , ( t max ) + z z H
where the value of γ ^ l 0 , ( t max ) can be obtained by (A9) and (A10). For convenience, we define Σ 0 = l = 1 L z l z l H γ ^ l 0 , ( t max ) + z z H . Substituting (24) into (23) and ignoring the constant terms results in the AE-Rao as
t AE - Rao = z H Σ 0 1 H H H Σ 0 1 H 1 H H Σ 0 1 z

3.3. AE-Wald Detector

The Wald design criterion can be expressed as [29]
t Wald = Θ ^ r 1 Θ r 0 H I 1 Θ ^ 1 Θ r , Θ r 1 Θ ^ r 1 Θ r 0
Nulling the derivative of (3) with respect to ξ yields the MLE of ξ for the given R as
ξ ^ = H H R 1 H 1 H H R 1 z
Plugging (22) and (27) into (26) results in the AE-Wald for the given R as
t AE - Wal d R = z H R 1 H ( H H R 1 H ) 1 H H R 1 z
Substituting the final alternate estimations of γ 1 and ξ into (6) results in the MLE of R under hypothesis H 1 as
R ^ 1 = 1 L + 1 l = 1 L z l z l H γ ^ l 1 , ( t max ) + ( z H ξ ^ ( t max ) ) ( z H ξ ^ ( t max ) ) H
where the values of γ ^ l 1 , ( t max ) and ξ ^ ( t max ) can be obtained by iterating (A4) and (A8) to the t max . For convenience, we define Σ 1 = l = 1 L z l z l H γ ^ l 1 , ( t max ) + ( z H ξ ^ ( t max ) ) ( z H ξ ^ ( t max ) ) H . Substituting (29) into (28) and ignoring the constant terms lead to the AE-Wald as
t AE - Wald = z H Σ 1 1 H H H Σ 1 1 H 1 H H Σ 1 1 z

3.4. AE-Gradient Detector

The Gradient design criterion can be written as [30]
t Gradient = Re ln f 1 ( z , Z ) Θ r T Θ = Θ ^ 0 Θ ^ r 1 Θ r 0
Substituting (16) and (27) into (31) results in the AE-Gradient for the given R as
t AE - Gradient R = z H R 1 H ( H H R 1 H ) 1 H H R 1 z
The MLE of R under hypotheses H 1 and H 0 are (29) and (24), respectively. Substituting (30) and (24) into (32), simplifying with the definition of Σ 0 and Σ 1 , and ignoring the constant terms yield the AE-Gradient as
t AE - Gradient = z H Σ 0 1 H H H Σ 1 1 H 1 H H Σ 1 1 z

3.5. AE-Durbin Detector

The Durbin design criterion can be formulated as [30]
t Durbin = Θ ^ r , 01 Θ r 0 H I Θ ^ 0 Θ r , Θ r I 1 Θ ^ 0 Θ r , Θ r I Θ ^ 0 Θ r , Θ r Θ ^ r , 01 Θ r 0
where Θ ^ r , 01 is the MLE of Θ r under hypothesis H 1 when Θ s is set to its MLE under hypothesis H 0 . Substituting (18), (22) and (27) into (34) results in the Durbin detector as
t AE - Durbin = z H R ^ 1 1 H H H R ^ 1 1 H 1 H H R ^ 0 1 H H H R ^ 1 1 H 1 H H R ^ 1 1 z
In the same way, substituting (29) and (24) into (35), simplifying with the definition of Σ 0 and Σ 1 , and ignoring the constant terms yield the AE-Durbin as
t AE - Durbin = z H Σ 1 1 H H H Σ 1 1 H 1 z H Σ 0 1 H H H Σ 1 1 H 1 H H Σ 1 1 z

4. Performance Evaluation

We evaluate the proposed detectors with simulated and measured data in this section. To measure the relative variation of the clutter covariance matrix R across consecutive iterations, we define the relative variation function
Δ R = R ^ ( t ) F R ^ ( t 1 ) F R ^ ( t 1 ) F
where · F represents the Frobenius norm of a matrix and R ^ ( t ) denotes the estimation of R after t iterations.
Figure 1 illustrates the relative change in the estimation of the covariance matrix R over successive iterations under hypotheses H 0 and H 1 . Obviously, as the number of iterations increases, the degree of relative change in R gradually decreases. When the number of iterations exceeds 15, the degree of relative change gradually slows down. Considering the time complexity of the algorithm, the number of iterations for the power mismatches γ l and the coordinate vector ξ is set to be 15.

4.1. Simulation Experiment

First, we set the number of system channels to 8, the number of training samples to 16, the dimension of signal subspace to 3, and the expected PFA to 10 3 . A total of 100 / PFA Monte Carlo simulations are run to derive the detection threshold, and 10 4 independent repeated tests are used to calculate the PD.
In the simulation experiment, the speckle component follows a zero-mean complex Gaussian random process, and its covariance matrix is [31]
R j , k = δ j k , 1 j , k N
where δ denotes the one-lag correlation coefficient, set to 0.95.
In this paper, the SCR is defined as [32]
SC R out = ξ H H H R 1 H ξ
In the experiment, we set ξ = α ξ 0 , where H and ξ 0 are randomly generated, but fixed after generation. α is the set regulatory factors to meet the given SCR as
α = SCR out ξ 0 H H H R 1 H ξ 0
In addition, the detection performance when the signal mismatch occurs is also discussed in this paper. For subspace signals, signal mismatch refers to the fact that the actual signal s 0 is not entirely located in the assumed subspace H . For the purpose of quantifying the degree of signal mismatch, we employ the generalized cosine square [33]
co s 2 ϕ = s 0 H R 1 H ( H H R 1 H ) 1 H H R 1 s 0 s 0 H R 1 s 0
We consider K distribution clutter, whose texture component follows a gamma distribution with the PDF as [18]
p τ ( τ ) = 1 Γ v v u v τ v 1 exp v u τ
where v is the shape parameter, representing the degree of the clutter trailing, and u is the scale parameter, representing the size of clutter power. Note that the smaller the shape parameter, the heavier the clutter trailing degree and the stronger the non-Gaussian characteristics. In the experiment, we set v to 1. Under the condition of a given texture component τ , the condition PDF of the clutter amplitude is
f ( r | τ ) = 2 r τ exp r 2 τ
From the full probability formula, the PDF of the clutter amplitude is
f ( r ) = 0 p ( τ ) f r | τ d τ
Substituting (42) and (43) into (44) results in
f ( r ) = 0 p τ f r | τ d τ = 0 v v Γ v u v τ v 1 2 r τ exp v u τ r 2 τ d τ = 2 2 v 1 Γ v u 2 r u v K v 1 2 r u
where K v · represents the modified Bessel function of the second kind and v-th order. The cumulative distribution function (CDF) of the K distribution can be obtained by integrating (45) as
F ( r ) = 1 1 2 v 1 Γ ( v ) 2 r u v K v 2 r u
In Figure 2, we use simulated data to compare the variation in PDs of the proposed detectors with respect to SCR; we can find that AE-GLRT, AE-Rao, AE-Wald, and AE-Gradient have similar PDs, but the performance of AE-Rao is slightly better than the other three proposed detectors. Specifically, when the SCR = 27 dB , AE-Rao can detect almost exactly for the given PFA. However, the PD of AE-Durbin is relatively low, when the SCR = 27 dB, PD is close to 0.91; when the SCR = 35 dB, PD is close to 1.
Figure 3 depicts the variation in PDs when signal mismatch occurs at SCR = 35 dB. Signal mismatch refers to the case that the actual signal s 0 is not entirely located within the nominal signal subspace spanned by the columns of H . A quantity to measure the amount of signal mismatch is given in (41). Note that the case cos 2 ϕ = 1 indicates that the actual signal completely lies in the signal subspace assumed by the radar system. In contrast, cos 2 ϕ = 0 means that the actual signal is orthogonal to the signal subspace assumed by the radar system. We use Figure 3 to measure the robustness or selectivity of the proposed detectors to mismatched signals. The robustness of a detector refers to the ability of the detector to detect targets with a high probability when signal mismatch occurs. The selectivity of a detector refers to the ability to decrease the PD of a detector when signal mismatch occurs. Figure 3 indicates that AE-Rao and AE-GLRT have better selectivity, while AE-Durbin is the most robust.
Figure 4 demonstrates the PDs at different SNRs with v as a parameter. It is evident that as the shape parameter v decreases, indicating stronger non-Gaussian characteristics of the clutter, there is an enhancement in the PDs of all detectors. Moreover, as indicated in Figure 4a,b, the PD of AE-Rao remains the highest.

4.2. Real Data Experiment

We utilize a set of sea clutter measured data with a range resolution of 60 m from the IPIX database in 1998, named “19980204_223753_ANTSTEP.mat”, including 60,000 pulses and 28 range cells. We selected the 13th range cell as the CUT, divided the data within it into 200 equally spaced segments, and used the histogram density estimation method to plot the empirical PDF and CDF curves of the real data. To verify the rationality of the selected measured data, we used the fractional moment estimation method (MoFM) [34,35,36] for the real data to calculate the parameters of three common compound-Gaussian distributions, namely, the shape parameters and scale parameters of K distribution, Pareto and CG−IG distributions, and plotted the amplitude PDF and CDF curves of the three compound-Gaussian distributions. Specifically, to reduce the error caused by abnormal samples, for K-distributed clutter, we use the first moment and the half-order moment of the clutter samples for parameter estimation, which are denoted by m ^ 1 and m ^ 0.5 , with the expressions as follows
m ^ 0.5 = Γ 1.25 Γ v + 0.25 Γ v 2 u 1 2 m ^ 1 = Γ 1.5 Γ v + 0.5 Γ v 2 u
where (47) serves as an implicit function of v. It is worth mentioning that the gamma function is defined as Γ ( z ) = 0 + t z 1 e t d t . So we can eliminate u, unfold, and simplify (47), learning that v has a unique solution. Next, we use the lookup table method or employ the f s o l v e · function in MATLAB to find the solution. As shown in Figure 5 and Figure 6, the K distribution provides a better fit to the real data.
To further measure the fitting degree of the three compound-Gaussian distributions mentioned above to the real data, Kolmogorov-Smirnov (KS) distance is commonly used as a metric. In addition, we plot the fitting results of the clutter amplitude distribution of real data in Figure 7. Here, the vertical axis represents the 1 CDF domain, which is equivalent to the meaning of the PFA [37]. As the PFA we set is 10 2 , we narrow the scope of the KS test to around 10 2 , and the corresponding expression of the KS distance can be written as
K S ( v , u ; v ^ , u ^ ) = max r ( 2.5 , 4.5 ) F ( r ; v , u ) F ( r ; v ^ , u ^ )
where F ( r ; v , u ) represents the empirical CDF obtained from the real data and F ( r ; v ^ , u ^ ) represents the fitting CDF obtained through parameter estimation based on the selected distribution model. Table 1 indicates that the KS distance for the K distribution fit is 0.0011 for the PFA we set. As we know, a smaller KS distance indicates a better fit between the empirical CDF of clutter amplitude and the selected clutter model. So, we can draw a conclusion that the K distribution still offers a good fit. It is worth noting that using real data with good fitting results from the K distribution for experiments aims to correspond better with simulation experiments. In this paper, the power mismatch of the training samples obtained through iterative estimation does not need to consider the specific distribution of clutter.
We set PFA = 10 2 due to the limited availability of pulses in the real data, and selected the 13 L / 2 , 12 and 14 , 13 + L / 2 range cells near the CUT as training samples. The other parameters are consistent with the simulation experiment.
In Figure 8, we graph the PDs of the proposed detectors across various SCRs in the absence of signal mismatch. The observation reveals that the results from the real data are basically consistent with the simulation experiment. The detection performances of AE-GLRT, AE-Rao, AE-Wald, and AE-Gradient are comparable and significantly better than AE-Durbin. Specifically, AE-Rao still has the highest PD, and the performances of AE-Wald and AE-Gradient slightly surpass that of AE-GLRT.
In Figure 9, utilizing the measured data from the IPIX radar, we compare and analyze the variation in the PDs of the proposed detectors in relation to the degree of signal mismatch. It can be observed that the result of the measured data is also basically consistent with the simulation experiment, indicating that AE-Durbin has the strongest robustness; when the signal mismatch angle is approximately 53 degrees, AE-Durbin can still detect the target with a PD of 0.9, while both AE-GLRT and AE-Rao have good selectivity, but AE-GLRT has stronger selectivity, indicating a better ability to suppress mismatched signals.
Meanwhile, we also conduct experiments with “19980204_224024_ANTSTEP.mat” real data, which is well fitted by the Pareto distribution, and the results are basically consistent with the above. We can reasonably infer that the detectors proposed in this paper have a certain universality in power heterogeneous clutter environments.

5. Conclusions

In this article, we investigated the adaptive detection problem of subspace signals within power heterogeneous clutter environments. Based on the GLRT, Rao, Wald, Gradient, and Durbin criteria, we used an alternate estimation iterative algorithm to design five adaptive detectors: AE-GLRT, AE-Rao, AE-Wald, AE-Gradient, and AE-Durbin. The simulation experiments and measured data illustrated that AE-Rao had a relatively higher PD when there was no signal mismatch. However, in the case of signal mismatch, AE-GLRT had the best selectivity, whereas AE-Durbin had the most robust detection performance.

Author Contributions

Conceptualization, H.C.; methodology, D.X.; software, H.L.; validation, B.L.; formal analysis, H.C.; investigation, W.L.; writing—original draft preparation, D.X.; writing—review and editing, W.L.; visualization, H.L.; supervision, W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under grant 62071482.

Data Availability Statement

The data utilized in this study are available from the first or corresponding authors upon a reasonable request.

Acknowledgments

The authors would like to express appreciation for the valuable comments of the associate editors and anonymous reviewers.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

In order to find the final value of the alternate estimation of ξ under hypothesis H 1 , we set
f ( ξ ) = ( z H ξ ) ( z H ξ ) H + A ( t )
where A ( t ) = l = 1 L z l z l H γ ^ l 1 , ( t ) and the superscript t indicates the t-th estimation of the matrix A .
By the matrix determinant lemma, for matrices A , B , C , and D of compatible dimensions, it holds that
A + B D H = A I + D H A 1 B
so (A1) can be simplified to
f ( ξ ) = A ( t ) 1 + z H ξ H A ( t ) 1 z H ξ
Calculating the partial derivative of (A3) with respect to ξ and setting it to zero result in
ξ ^ ( t + 1 ) = H H A ( t ) 1 H 1 H H A ( t ) 1 z
where the superscript t + 1 indicates the ( t + 1 ) -th estimation of ξ . Next, solve for γ ^ h , h = 1 , 2 , , L . From (8), we define
f ( γ h ) = γ h N ( z H ξ ^ ( t + 1 ) ) ( z H ξ ^ ( t + 1 ) ) H + l = 1 h 1 z l z l H γ ^ l 1 , ( t + 1 ) + l = h + 1 L z l z l H γ ^ l 1 , ( t ) + z h z h H γ h L + 1 = γ h N B h 1 , ( t + 1 ) + z h z h H γ h L + 1 = γ h N B h 1 , ( t + 1 ) L + 1 1 + z h H B h 1 , ( t + 1 ) 1 z h γ h L + 1
where
B h 1 , t + 1 = ( z H ξ ^ ( t + 1 ) ) ( z H ξ ^ ( t + 1 ) ) H + l = 1 h 1 z l z l H γ ^ l 1 , ( t + 1 ) + l = h + 1 L z l z l H γ ^ l 1 , ( t )
The above simplification used (A2). Upon analyzing (A5), it can be seen that when γ h 0 , f ( γ h ) ; when γ h , f ( γ h ) . Therefore, a certain value within 0 , exists that minimizes f ( γ h ) . Taking the partial derivative of (A5) with respect to γ h results in
f ( γ h ) γ h = N γ h N 1 B h 1 , ( t + 1 ) L + 1 1 + 1 γ h z h H B h 1 , ( t + 1 ) 1 z h L + 1 + γ h N B h 1 , ( t + 1 ) L + 1 L + 1 1 + 1 γ h z h H B h 1 , ( t + 1 ) 1 z h L × 1 γ h 2 z h H B h 1 , ( t + 1 ) 1 z h
We can further obtain that
γ ^ h 1 , ( t + 1 ) = L + 1 N N z h H B h 1 , ( t + 1 ) 1 z h
where h = 1 , 2 , , L . To date, we can give an initial value of γ ^ l 1 , ( 0 ) , l = 1 , 2 , , L and loop from (A4) and (A8) to the termination condition, and we can obtain the final estimation of ξ ^ ( t max ) and γ ^ l 1 , ( t max ) , l = 1 , 2 , , L under hypothesis H 1 .
Obviously, we can set ξ = 0 p × 1 under hypothesis H 0 . Then, we also need to calculate the estimations of γ l , l = 1 , 2 , , L under hypothesis H 0 , which is similar to the derivation under hypothesis H 1 , so we directly give the following expressions
B h 0 , t + 1 = z z H + l = 1 h 1 z l z l H γ ^ l 0 , ( t + 1 ) + l = h + 1 L z l z l H γ ^ l 0 , ( t )
and
γ ^ h 0 , ( t + 1 ) = L + 1 N N z h H B h 0 , ( t + 1 ) 1 z h

References

  1. Zhang, C.; Xu, Y.; Chen, W.; Chen, B.; Gao, C.; Liu, H. IfCMD: A Novel Method for Radar Target Detection under Complex Clutter Backgrounds. Remote Sens. 2024, 16, 2199. [Google Scholar] [CrossRef]
  2. Xue, J.; Fan, Z.; Xu, S.; Liu, J. Persymmetric Adaptive Radar Target Detection in CG-LN Sea Clutter Using Complex Parameter Suboptimum Tests. IEEE Trans. Geosci. Remote Sens. 2023, 61, 5111211. [Google Scholar] [CrossRef]
  3. Liu, W.; Han, H.; Liu, J.; Li, H.; Li, K.; Wang, Y.L. Multichannel radar adaptive signal detection in interference and structure nonhomogeneity. Sci. China Inf. Sci. 2017, 60, 112302. [Google Scholar] [CrossRef]
  4. Xue, J.; Yan, J.; Pan, M.; Xu, S. Knowledge-aided adaptive Gradient test for radar targets in correlated compound Gaussian sea clutter with lognormal texture. IEEE Geosci. Remote Sens. Lett. 2023, 20, 3509105. [Google Scholar] [CrossRef]
  5. Liu, W.; Liu, J.; Hao, C.; Gao, Y.; Wang, Y.L. Multichannel adaptive signal detection: Basic theory and literature review. Sci. China Inf. Sci. 2022, 65, 121301. [Google Scholar] [CrossRef]
  6. Kelly, E.J. An Adaptive Detection Algorithm. IEEE Trans. Aerosp. Electron. Syst. 1986, 22, 115–127. [Google Scholar] [CrossRef]
  7. 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]
  8. Liu, W.; Zhang, Z.; Liu, J.; Shang, Z.; Wang, Y.L. Detection of a rank-one signal with limited training data. Signal Process. 2021, 186, 108120. [Google Scholar] [CrossRef]
  9. Kraut, S.; Scharf, L.L. The CFAR adaptive subspace detector is a scale-invariant GLRT. IEEE Trans. Signal Process. 1999, 47, 2538–2541. [Google Scholar] [CrossRef]
  10. Xu, S.; Shi, X.; Xue, J.; Shui, P. Adaptive subspace detection of range-spread target in compound Gaussian clutter with inverse Gaussian texture. Digit. Signal Process. 2018, 81, 79–89. [Google Scholar] [CrossRef]
  11. Xue, J.; Xu, S.; Shui, P. Near-optimum coherent CFAR detection of radar targets in compound-Gaussian clutter with inverse Gaussian texture. Signal Process. 2020, 166, 107236. [Google Scholar] [CrossRef]
  12. Xue, J.; Xu, S.; Shui, P. Knowledge-based adaptive detection of radar targets in generalized Pareto clutter. Signal Process. 2018, 143, 106–111. [Google Scholar] [CrossRef]
  13. Balleri, A.; Nehorai, A.; Wang, J. Maximum likelihood estimation for compound-Gaussian clutter with inverse gamma texture. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 775–779. [Google Scholar] [CrossRef]
  14. Ollila, E.; Tyler, D.E.; Koivunen, V.; Poor, H.V. Compound–Gaussian clutter modeling with an inverse Gaussian texture distribution. IEEE Signal Process. Lett. 2012, 19, 876–879. [Google Scholar] [CrossRef]
  15. Shang, X.; Song, H. Radar detection based on compound-Gaussian model with inverse gamma texture. IET Radar Sonar Navig. 2011, 5, 315–321. [Google Scholar] [CrossRef]
  16. Jiang, Q.; Wu, Y.; Liu, W.; Zheng, D.; Jian, T.; Gong, P. Subspace-based distributed target detection in compound-Gaussian clutter. Digit. Signal Process. 2023, 140, 104141. [Google Scholar] [CrossRef]
  17. Dong, Y. Optimal coherent radar detection in a K-distributed clutter environment. IET Radar Sonar Navig. 2012, 6, 283–292. [Google Scholar] [CrossRef]
  18. Shui, P.l.; Liu, M.; Xu, S.w. Shape-parameter-dependent coherent radar target detection in K-distributed clutter. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 451–465. [Google Scholar] [CrossRef]
  19. Chen, Z.; Chen, A.; Liu, W.; Ma, X. CFAR Detection in Nonhomogeneous Weibull Sea Clutter for Skywave OTHR. IEEE Geosci. Remote Sens. Lett. 2023, 20, 3508005. [Google Scholar] [CrossRef]
  20. Trunk, G.; George, S. Detection of targets in non-Gaussian sea clutter. IEEE Trans. Aerosp. Electron. Syst. 1970, 5, 620–628. [Google Scholar] [CrossRef]
  21. Xu, S.W.; Xue, J.; Shui, P.L. Adaptive detection of range-spread targets in compound Gaussian clutter with the square root of inverse Gaussian texture. Digit. Signal Process. 2016, 56, 132–139. [Google Scholar] [CrossRef]
  22. Coluccia, A.; Orlando, D.; Ricci, G. A GLRT-like CFAR detector for heterogeneous environments. Signal Process. 2022, 194, 108401. [Google Scholar] [CrossRef]
  23. Liu, W.; Liu, J.; Huang, L.; Yang, Z.; Yang, H.; Wang, Y.L. Distributed target detectors with capabilities of mismatched subspace signal rejection. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 888–900. [Google Scholar] [CrossRef]
  24. Tang, P.; Wang, Y.L.; Liu, W.; Du, Q.; Wu, C.; Chen, W. A tunable detector for distributed target detection in the situation of signal mismatch. IEEE Signal Process. Lett. 2019, 27, 151–155. [Google Scholar] [CrossRef]
  25. Liu, W.; Liu, J.; Gao, Y.; Wang, G.; Wang, Y.L. Multichannel signal detection in interference and noise when signal mismatch happens. Signal Process. 2020, 166, 107268. [Google Scholar] [CrossRef]
  26. Coluccia, A.; Ricci, G. A tunable W-ABORT-like detector with improved detection vs. rejection capabilities trade-off. IEEE Signal Process. Lett. 2014, 22, 713–717. [Google Scholar] [CrossRef]
  27. Conte, E.; De Maio, A.; Ricci, G. GLRT-based adaptive detection algorithms for range-spread targets. IEEE Trans. Signal Process. 2001, 49, 1336–1348. [Google Scholar] [CrossRef]
  28. Liu, B.; Chen, B.; Michels, J.H. A GLRT for multichannel radar detection in the presence of both compound Gaussian clutter and additive white Gaussian noise. Digit. Signal Process. 2005, 15, 437–454. [Google Scholar] [CrossRef]
  29. Sun, M.; Liu, W.; Liu, J.; Hao, C. Rao and wald tests for target detection in coherent interference. IEEE Trans. Aerosp. Electron. Syst. 2021, 58, 1906–1921. [Google Scholar] [CrossRef]
  30. Besson, O. Adaptive detection of Gaussian rank-one signals using adaptively whitened data and Rao, gradient and Durbin tests. IEEE Signal Processing Letters 2023, 30, 399–402. [Google Scholar] [CrossRef]
  31. Gini, F.; Greco, M. Covariance matrix estimation for CFAR detection in correlated heavy tailed clutter. Signal Process. 2002, 82, 1847–1859. [Google Scholar] [CrossRef]
  32. Liu, W.; Wang, Y.L.; Liu, J.; Huang, L.; Hao, C. Performance analysis of adaptive detectors for point targets in subspace interference and Gaussian noise. IEEE Trans. Aerosp. Electron. Syst. 2017, 54, 429–441. [Google Scholar] [CrossRef]
  33. Yufeng, C.; Yongliang, W.; Weijian, L.; Qinglei, D.; Xichuan, Z.; Xuhui, L. A tunable adaptive detector for distributed targets when signal mismatch occurs. J. Syst. Eng. Electron. 2023, 34, 873–878. [Google Scholar]
  34. Xue, J.; Xu, S. Parameters estimation based on moments and Nelder-Mead algorithm for compound-Gaussian clutter with inverse Gaussian texture. In Proceedings of the 2016 IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC), Hong Kong, China, 5–8 August 2016; pp. 1–5. [Google Scholar]
  35. Iskander, D.R.; Zoubir, A.M. Estimation of the parameters of the K-distribution using higher order and fractional moments [radar clutter]. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 1453–1457. [Google Scholar] [CrossRef]
  36. Gini, F.; Greco, M. Texture modelling, estimation and validation using measured sea clutter data. IEE Proc.-Radar Sonar Navig. 2002, 149, 115–124. [Google Scholar] [CrossRef]
  37. Chan, H. Radar sea-clutter at low grazing angles. IEE Proc. F (Radar Signal Process.) 1990, 137, 102–112. [Google Scholar] [CrossRef]
Figure 1. The relative variation in the estimation of the covariance matrix across successive iterations.
Figure 1. The relative variation in the estimation of the covariance matrix across successive iterations.
Remotesensing 16 02508 g001
Figure 2. Probability of detection (PD) versus signal-to-clutter ratio (SCR) without signal mismatch with simulated data.
Figure 2. Probability of detection (PD) versus signal-to-clutter ratio (SCR) without signal mismatch with simulated data.
Remotesensing 16 02508 g002
Figure 3. PD versus cos 2 ϕ with signal mismatch with simulated data. SCR = 35 dB .
Figure 3. PD versus cos 2 ϕ with signal mismatch with simulated data. SCR = 35 dB .
Remotesensing 16 02508 g003
Figure 4. PD versus SCR without signal mismatch with simulated data for different v.
Figure 4. PD versus SCR without signal mismatch with simulated data for different v.
Remotesensing 16 02508 g004
Figure 5. Compound-Gaussian distribution PDF fitting to empirical PDF of real data.
Figure 5. Compound-Gaussian distribution PDF fitting to empirical PDF of real data.
Remotesensing 16 02508 g005
Figure 6. Compound-Gaussian distribution CDF fitting to empirical CDF of real data.
Figure 6. Compound-Gaussian distribution CDF fitting to empirical CDF of real data.
Remotesensing 16 02508 g006
Figure 7. Compound-Gaussian distribution modeling results of real data.
Figure 7. Compound-Gaussian distribution modeling results of real data.
Remotesensing 16 02508 g007
Figure 8. PD versus SCR without signal mismatch with IPIX radar data.
Figure 8. PD versus SCR without signal mismatch with IPIX radar data.
Remotesensing 16 02508 g008
Figure 9. PD versus cos 2 ϕ with signal mismatch with IPIX radar data. SCR = 25 dB .
Figure 9. PD versus cos 2 ϕ with signal mismatch with IPIX radar data. SCR = 25 dB .
Remotesensing 16 02508 g009
Table 1. Kolmogorov-Smirnov (KS) distances for each distribution fit.
Table 1. Kolmogorov-Smirnov (KS) distances for each distribution fit.
DistributionKS Distance
K distribution0.0011
Pareto distribution0.0066
CG-IG distribution0.0013
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

Xiao, D.; Liu, W.; Chen, H.; Li, H.; Li, B. An Adaptive Radar Target Detection Method Based on Alternate Estimation in Power Heterogeneous Clutter. Remote Sens. 2024, 16, 2508. https://doi.org/10.3390/rs16132508

AMA Style

Xiao D, Liu W, Chen H, Li H, Li B. An Adaptive Radar Target Detection Method Based on Alternate Estimation in Power Heterogeneous Clutter. Remote Sensing. 2024; 16(13):2508. https://doi.org/10.3390/rs16132508

Chicago/Turabian Style

Xiao, Daipeng, Weijian Liu, Hui Chen, Hao Li, and Binbin Li. 2024. "An Adaptive Radar Target Detection Method Based on Alternate Estimation in Power Heterogeneous Clutter" Remote Sensing 16, no. 13: 2508. https://doi.org/10.3390/rs16132508

APA Style

Xiao, D., Liu, W., Chen, H., Li, H., & Li, B. (2024). An Adaptive Radar Target Detection Method Based on Alternate Estimation in Power Heterogeneous Clutter. Remote Sensing, 16(13), 2508. https://doi.org/10.3390/rs16132508

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