Next Article in Journal
Depth Map Upsampling via Multi-Modal Generative Adversarial Network
Previous Article in Journal
Dynamic pH Sensor with Embedded Calibration Scheme by Advanced CMOS FinFET Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Maritime Moving Target Detection and Tracking Method for Airborne Forward-looking Scanning Radar

School of Communication and Information Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
*
Authors to whom correspondence should be addressed.
Sensors 2019, 19(7), 1586; https://doi.org/10.3390/s19071586
Submission received: 29 January 2019 / Revised: 29 March 2019 / Accepted: 29 March 2019 / Published: 2 April 2019
(This article belongs to the Section Remote Sensors)

Abstract

:
Maritime moving target detection and tracking through particle filter based track-before-detect (PF-TBD) has significant practical value for airborne forward-looking scanning radar. However, villainous weather and surging of ocean waves make it extremely difficult to accurately obtain a statistical model for sea clutter. As the likelihood ratio calculation in PF-TBD is dependent on the distribution of the clutter, the performance of traditional distribution-based PF-TBD seriously declines. To resolve these difficulties, this paper proposes a new target detection and tracking method, named spectral-residual-binary-entropy-based PF-TBD (SRBE-PF-TBD), which is independent from the prior knowledge of sea clutter. In the proposed method, the likelihood ratio calculation is implemented by first extracting the spectral residual of the input image to obtain the saliency map, and then constructing likelihood ratio through a binarization processing and information entropy calculation. Simulation results show that the proposed method had superior performance of maritime moving target detection and tracking.

1. Introduction

The rapid development of global trade has greatly promoted the demand for marine transport and tourism, which also increases marine accidents [1,2,3]. Therefore, marine search and rescue has become an important research issue and receives a lot of attention in recent years [4,5]. The primary task of marine search and rescue is to discover the target under complex sea condition quickly and accurately, especially small targets, such as lifeboats carrying survivors. The airborne forward-looking scanning radar is a main detection means for marine search and rescue because of its all-weather and all-time capacity [6,7,8,9]. Through the sequential scanning of radar, multi-frame data can be quickly acquired and the maritime target in the observation area can be monitored in real time.
In general, fine weather condition and calm sea surface benefit target detection and tracking because of the high signal clutter ratio (SCR) of the radar echo. In this situation, the target can be effectively detected by traditional constant false alarm rate (CFAR) technologies [10,11,12], and then the trajectory can be estimated through the tracking filter [13,14]. However, villainous weather and surging of ocean waves will lead to heavy-tailed clutter and fluctuation of target scattering. In addition to the weak scattering of the target, the SCR will become quite low. In this case, traditional detection and tracking method will generate massive false alarms, making it hard to detect the true target [15].
Different from the traditional detection and tracking method, the TBD method utilizes the motion characteristics of the target and processes consecutive scans jointly to realize the fine detection and tracking of low SCR targets. Therefore, the TBD method has drawn extensive attention [16]. TBD includes batch methods, such as dynamic programming [17] and Hough transform [18], and recursive methods based on the Bayesian approach, such as particle filter (PF) [19]. Because PF can effectively deal with nonlinear non-Gaussian problems, it adapts to the target detection and tracking under sea scene. The earliest use of the PF-TBD method in radar system is presented in [20]. Subsequently, Rutten et al. proposed a more efficient implementation method of PF-TBD [21], and analyzed the characteristics of the method in detail [22]. Recently, PF-TBD is applied in different radar systems, such as over-the-horizon radar (OTHR) [23], passive radars [24], SAR [25] and multi-station system [26]. However, these studies consider the Rayleigh clutter distribution, which is only applicable to low resolution radars or low sea conditions. When the high-resolution radar observes the sea surface at a small grazing angle, the clutter distribution presents a significant non-Gaussian characteristic [27]. In [28], the PF-TBD method based on K-distribution adapts to the non-Gaussian characteristic of the sea clutter amplitude, and the calculation of likelihood ratio is deduced in detail. However, due to the complexity of the K-distribution, the analytical solution of the likelihood ratio calculation cannot be obtained, and this method must be resolved by a time-consuming numerical integration operation. This is a great challenge for real-time demand in applications.
Moreover, the distribution-based methods need to assume the prior distribution of clutter. However, the complex ocean dynamics mechanism gives the sea surface non-stationary characteristic [29], which means it is difficult to accurately model the sea surface with a hypothetical distribution model. To reduce the dependence on the prior knowledge of clutter, the saliency detection approach is proposed based on visual saliency principle [30,31,32]. A fast saliency detection approach proposed in [33] obtains the saliency map by analyzing the log-spectrum of an input image and extracting the spectral residuals of the image in the frequency domain, and then realize target detection.
Inspired by the saliency detection approach, this paper proposes a new target detection and tracking method, named spectral-residual-binary-entropy-based PF-TBD (SRBE-PF-TBD), which is used to detect and track maritime moving target for airborne forward-looking scanning radar. In the proposed method, the likelihood ratio calculation is implemented by first extracting the spectral residual of the input image to obtain the saliency map, and then constructing likelihood ratio through a binarization processing and information entropy calculation.
The proposed SRBE-PF-TBD method does not require prior knowledge of the clutter, enhancing the adaptability and robustness. The effectiveness and performance of the proposed method are verified by simulation experiments based on the simulated data, which are obtained by a Swerling I target and the K-distributed clutter.
This paper is arranged as follows. Section 2 introduces the target dynamic model and forward-looking scanning radar measurement model. Section 3 shows the proposed SRBE-PF-TBD method. Section 4 presents the simulation results and performance analysis. Section 5 gives the conclusions.

2. Models

In this section, the target dynamic model and the airborne forward-looking scanning radar measurement model are given, which are used in this paper.

2.1. Target Dynamic Model

Considering the maritime target moving in a two-dimensional plane, the state vector of the target in frame k is represented as x k = x k , x ˙ k , y k , y ˙ k T , where x k , y k is the position of the target and x ˙ k , y ˙ k denotes the velocity along x-axis and y-axis, respectively. Assuming a constant-velocity process model is used, the evolution of the target state can be written as
x k + 1 = F x k + w k
where
F = 1 τ 0 0 0 1 0 0 0 0 1 τ 0 0 0 1
τ is the frame sampling period, and w k represents the process noise, which obeys Gaussian distribution with mean zero and covariance Q :
Q = q × τ 3 / 3 τ 2 / 2 0 0 τ 2 / 2 τ 0 0 0 0 τ 3 / 3 τ 2 / 2 0 0 τ 2 / 2 τ
where q is the process noise coefficient [34].
To indicate the existence or non-existence of target in the observation area, the random variable ε k is introduced as the target existence state variable and ε k 0 , 1 , where ε k = 1 means there is a target in the observation area in frame k and vice versa. The above process can be modeled as a two-state Markov chain, as shown in Figure 1.
In Figure 1, P b i r t h is the probability of target birth and P d e a t h is the probability of target death, which are respectively equal to P ε k = 1 | ε k 1 = 0 and P ε k = 0 | ε k 1 = 1 . Thus, the two state Markov chain can be defined as
Π = 1 P b i r t h P b i r t h P d e a t h 1 P d e a t h

2.2. Airborne Forward-Looking Scanning Radar Measurement Model

The geometry illustration for airborne forward-looking scanning radar is shown in Figure 2. For the common overlapping area of multi-frame images, the observation area Ω k of frame k can be represented as
Ω k = σ x 1 , y 1 k σ x 1 , y 2 k σ x 1 , y N y k σ x 2 , y 1 k σ x 2 , y 2 k σ x 2 , y N y k σ x N x , y 1 k σ x N x , y 2 k σ x N x , y N y k
where σ x m , y n k m = 1 , 2 , , N x , n = 1 , 2 , , N y is the scattering coefficient of the point x m , y n , 0 in the observation area, and N x and N y are the numbers of resolution cells along x-axis and y-axis, respectively. Assume the airplane is moving along y-axis at speed v, the waveform of transmitted signal is s η , where η is the fast time, the scanning speed of radar beam is ω and the antenna pattern is w a t . At slow time t, the coordinate of radar is set as 0 , v t , H , and then the echo signal of the observation area in forward-looking scanning radar can be expressed as
u k η , t = m , n σ x m , y n k w a t s η 2 R x m , y n t c
where
R x m , y n t = R x m , y n 2 + v t 2 2 R x m , y n v t cos θ x m , y n cos φ
shows the instantaneous slant range of the point x m , y n , 0 , while R x m , y n is the initial slant range, θ x m , y n is the initial azimuth angle and φ is the grazing angle. In practice, due to the narrow azimuth beam, fast scanning speed and small grazing angle [7], the slant range can be approximated as
R x m , y n t R x m , y n v t cos φ
For the convenience of spatial data analysis, the time variable t and η can be represented by the spatial variables R and θ based on t = θ θ 0 θ θ 0 ω ω and η = 2 R 2 R c c , where θ 0 is initial scanning angle of the antenna. After the processing of matched filtering and range immigration correcting along the range dimension, the echo can be modeled as a convolution of the scattering coefficients and the two-dimensional system function [6], which is
u k R , θ = m , n σ ^ x m , y n k δ R R x m , y n , θ θ x m , y n h R , θ
where
σ ^ x m , y n k = σ x m , y n k · exp j 4 π λ R x m , y n t
and δ R R x m , y n , θ θ x m , y n is the Dirac function at R x m , y n , θ x m , y n . h R , θ is the two-dimensional system function, which can be written as
h R , θ = w a θ sinc 2 B c R
where B is the bandwidth of the transmitted signal.
In azimuth, the resolution of radar image can be significantly enhanced by super-resolution processing, contributing to an improved two-dimensional system function [35].
After coordinate transferring of Equation (9), the final image in spatial domain can be represented as
z k x , y = m , n σ ^ x m , y n k δ x x m , y y n h x m , y n
where h x m , y n is the two-dimensional system function in spatial domain. The imaging procedure of the forward-looking scanning radar is given in Figure 2.
Assuming the size of resolution cell in the observation area is Δ x × Δ y , the measurement data in frame k can be expressed as
z k = z k m , n : m = 1 , 2 , , N x , n = 1 , 2 , , N y
The hypothesis of the existence or non-existence of a target in resolution cell m , n can be defined as:
z k m , n = I k exp j ϕ k · h x m , y n + ς k m , n ε k = 1 ς k m , n ε k = 0
where I k is the amplitude of target, ϕ k is the random phase of target, and ς k m , n is the amplitude distribution of sea clutter, which is related with the sea condition and parameters of the surroundings.

3. The Proposed Method

In this section, the proposed SRBE-PF-TBD method is specified. The PF-TBD used in the proposed method is based on Bayesian recursion, which is derived in detail in [21]. In the PF-TBD, the particles are divided into two parts: the continuing particles, which represent the continuing presence of the target, and the birth particles, which represent the newly appearing target. By using the likelihood ratio to update the mixing weights of the two part particles, these particles can be mixed into a large set, and then the target state is estimated by the large set of particles.
Therefore, the likelihood ratio calculation plays an import role in the PF-TBD. As described above, the traditional distribution-based likelihood ratio calculation is infeasible in the situation of non-stationary sea surface. The proposed method implements likelihood ratio calculation from the perspective of image information, which is independent from the prior knowledge of clutter. Based on information theory, the information of an image can be decomposed into two parts: the prior knowledge and the innovation [33]. The prior knowledge part represents the statistical invariant properties of the environment, which can be considered as the redundant information and should be suppressed for target detection. The innovation part is the novel part, which contains information about target. Inspired by this rationale, the proposed method uses spectral residual calculation to remove the redundant information, and then the saliency map can be obtained. Finally, the likelihood ratio expression is constructed through a binary processing and information entropy calculation, which is named as SRBE likelihood ratio.
Based on the above, this section first presents the calculation of spectral residual binary entropy (SRBE) likelihood ratio, and then shows the implementation of the proposed SRBE-PF-TBD method.

3.1. SRBE Likelihood Ratio Calculation

Given the airborne forward-looking scanning radar image z k m , n , m = 1 , 2 , , N x , n = 1 , 2 , , N y of frame k. The Fourier spectrum of the image can be calculated by
Z k f m , f n = F z k ( m , n )
where F [ · ] means the operation of two-dimensional Fourier transform, and the corresponding amplitude spectrum and phase spectrum can be denoted as
A k f m , f n = A Z k f m , f n
P k f m , f n = P Z k f m , f n
where A · means getting the amplitude of the input, while P · obtains the phase. Then, the log spectrum of the image amplitude can be expressed as
G k f m , f n = ln A k f m , f n
Thus, the spectral residual can be obtained by
R k f m , f n = G k f m , f n χ r f m , f n G k f m , f n
where χ r f m , f n is a local mean filter, and defined by a r × r matrix:
χ r f m , f n = 1 r 2 1 1 1 1 1 1 1 1 1
By performing two-dimensional inverse Fourier transform and Gaussian filtering, the saliency map in spatial domain can be expressed as
S k m , n = g m , n F 1 exp R k f m , f n + P k f m , f n 2
where F 1 · means the operation of two-dimensional inverse Fourier transform. The Gaussian filer g m , n is defined by
g m , n = 1 2 π γ 2 exp m 2 + n 2 2 γ 2
where γ is the filter parameter.
Then, the normalized saliency map can be obtained by
S ¯ k m , n = S k m , n max S k m , n
By a given binarization threshold T B , the normalized saliency map is binarized and the corresponding binary image B k is generated, which indicates the potential target area. Then, the information entropy of B k is calculated by
H k = p 0 k · log 2 p 0 k p 1 k · log 2 p 1 k
where p 0 and p 1 , respectively, indicate the probabilities of the pixel value 0 and 1 in the binary image B k .
The final SRBE likelihood ratio can be written as
L S R B E z k | ε k = p 1 k p 0 k · exp S ¯ k m , n H k
For the target with state x k , the coordinate of the corresponding cell can be calculated by
m x k = x k Δ x n x k = y k Δ y
where · denotes the operation of rounding down. The SRBE likelihood ratio corresponding to x k can be calculated by
L S R B E z k | x k , ε k = p 1 k p 0 k · exp S ¯ k m x k , n x k H k
The flowchart of SRBE likelihood ratio calculation is shown in Figure 3.
As shown in Figure 3, the proposed SRBE method has two key operations: one is obtaining the saliency map of the radar image by the spectral residual method, and the other is calculating the binary entropy of the normalized saliency map. Because the maritime target is generally composed of metal material and has dihedral and trihedral structures, the scattering of the target is normally higher than the scattering of the sea clutter, meaning that the target is more salient than the sea clutter in the radar image. Therefore, the saliency map can effectively indicate the potential target regions. Through the information theory, the binary entropy of the saliency map can be used to measure the cluster-degree of the saliency map. The larger the entropy value is, the smaller the cluster-degree is, indicating that the saliency map is more likely caused by the sea clutter scattering. Contrarily, the lower entropy value represents the larger cluster-degree, which means the higher probability of target existence. Therefore, the binary entropy of the saliency map can effectively indicate the confidence of target existence. Based on the above analysis, the SRBE likelihood is constructed as Equation (27) to update the particles, which utilizes the information of the potential target regions and the confidence of target existence in the radar image.

3.2. Implementation of SRBE-PF-TBD

By incorporating the SRBE likelihood ratio into the PF-TBD, the specific implementation of the proposed SRBE-PF-TBD is described as follows.
  • Calculation of birth density
    A set of N b birth particles is drawn from the proposal density, whose positions are uniformly distributed within some highest measurements [22]. The sampled birth particles for frame k can be represented as
    x ( b ) i k q x ( b ) k | ε k = 1 , ε k 1 = 0 , z k , i 1 , 2 , , N b
    where q · is the proposal density, and the symbol “∼” means the operation of drawing samples from the proposal density.
    Based on the likelihood ratio proposed in Equation (27), the non-normalized weights of these birth particles can be calculated by
    ω b i k = L S R B E z k | x b i k , ε k = 1 p x b i k | ε k = 1 , ε k 1 = 0 N b q x b i k | ε k = 1 , ε k 1 = 0 , z k
    where the prior probability density p x b i k | ε k = 1 , ε k 1 = 0 is a uniform distribution. Then, the normalized weights can be obtained by
    ω ˜ b i k = ω b i k i = 1 N b ω b i k
  • Calculation of continuing density
    Since the auxiliary particle filter (APF) utilizes the information of measurement in the current frame, the knowledge of current observation is incorporated into the proposal mechanism, and particles are not moved blindly in the state space, making the current target state estimated with high reliability [36]. Considering these benefits, this paper adopts APF for continuing particle filtering.
    A set of N c continuing particles is drawn using the dynamic equation described in Equation (1), which can be expressed as
    μ ( c ) i k p x ( c ) k | x ( c ) k 1 , ε k = 1 , ε k 1 = 1 , i 1 , 2 , , N c
    where p · represents the transition density from frame k 1 to frame k.
    The non-normalized weights of these particles are calculated with the current measurements as below.
    b ( c ) i k = 1 N c L S R B E z k | μ ( c ) i k , ε k = 1
    The normalized weights of continuing particles can be calculated by
    b ˜ c i k = b c i k i = 1 N c b c i k
    By using weights b ˜ c i k , the particles x c i k 1 i = 1 N c in frame k 1 are resampled to obtain x c i j k 1 j = 1 N c , where i j j = 1 N c denotes the index of particle x c i k 1 in frame k 1 .
    Then, particles for frame k by Equation (1) and the resampled continuing particles in frame k 1 are again drawn, which can be expressed as
    x ( c ) j k p x ( c ) k | x ( c ) i j k 1 , ε k = 1 , ε k 1 = 1 , j 1 , 2 , , N c
    The non-normalized weights of newly resampled particles can be calculated by
    ω ( c ) j k = L S R B E z k | x ( c ) j k , ε k = 1 L S R B E z k | μ ( c ) i j k , ε k = 1
    The normalized weights can be obtained by
    ω ˜ c j k = ω c j k j = 1 N c ω c j k
  • Calculation of mixing probabilities
    The mixing probability of birth density is calculated by
    M b = P b i r t h 1 p ε k 1 = 1 | z 1 : k 1 i = 1 N b ω b i k
    Similarly, the mixing probability of continuing density is calculated by
    M c = 1 P d e a t h p ε k 1 = 1 | z 1 : k 1 i = 1 N c ω c i k
    Then, the normalized mixing probabilities can be expressed by
    M ˜ b = M b M b + M c
    M ˜ c = M c M b + M c
  • Calculation of scaled particle weights
    The particle weights are scaled according to the mixing probabilities, and the scaled particle weights of birth density and continuing density are obtained:
    ω ^ b i k = M ˜ b ω ˜ b i k
    ω ^ c i k = M ˜ c ω ˜ c i k
    Then, the particles of birth density and continuing density are combined into one large particle set:
    x j i k , ω ^ j i k | i 1 , . . . , N j , j b , c
    The particles of the large set are resampled to obtain N c particles with uniform weights:
    x ^ c i k , 1 1 N c N c | i 1 , . . . , N c
  • Calculation of probability of existence
    The probability of existence in frame k can be expressed as
    p ^ ε k = 1 | z 1 : k = M b + M c M b + M c + P d e a t h p ε k 1 = 1 | z 1 : k 1 + 1 P b i r t h 1 p ε k 1 = 1 | z 1 : k 1
  • Estimation of target state
    The detection threshold of existence is assumed as P t h . When p ^ ε k = 1 | z 1 : k > P t h , the hypothesis of target existence is accepted and the target state will be estimated by
    x ^ k = 1 N c i = 1 N c x ^ c i k

4. Experimental Results

To verify the proposed method, simulation experiments of maritime moving target detection and tracking based on simulated data were conducted.
The simulated data of airborne forward-looking scanning radar consisted of 30 frames, whose frame sampling period was τ = 1 s , and per frame contained N x × N y = 60 × 60 bins. The target appeared in frame K b = 6 , disappeared in frame K d = 21 , and moved according to the dynamic equation (Equation (1)), while the process noise coefficient was q = 0 . 001 . The initial state of the target was set as x 0 = 12 . 2 1 . 2 15 . 1 0 . 8 T .
There were 4000 particles, where the numbers used for birth density and continuing density were N b = 2000 and N c = 2000 , respectively. The initial particle velocities were sampled uniformly from interval 0 , 2 , and the particle positions were uniformly distributed within the 200 highest measurements. The probabilities of birth and death are defined as P b i r t h = P d e a t h = 0 . 05 , while the detection threshold was set as P t h = 0 . 5 . The binarization threshold T B was set to 0 . 5 . All simulation results were averaged over 100 Monte Carlo trials.

4.1. Simulated Data

In practice, since airborne forward-looking scanning radar is used for long-range detection with small grazing angle, it is reasonable to model the sea clutter amplitude as K-distribution [29], which can be expressed as
p ς k = 2 β Γ α β ς k 2 α K α 1 β ς k
where α is the shape parameter of the distribution, while β represents the scale parameter, and K α 1 · is the modified second-kind Bessel function of order α 1 .
Due to the movement of the ocean, the RCS of the maritime moving target varies between radar scans. Therefore, Swerling type I model was used to model the maritime moving target in our simulation [37]. The amplitude of the target can be written as
p I k = 2 I k σ ¯ exp I k 2 σ ¯
where σ ¯ denotes the mean squared target amplitude.
Thus, the SCR can be defined as
SCR = 10 log 10 σ ¯ P c l u t t e r
where P c l u t t e r means the mean power of the clutter in frame k and can be expressed by
P c l u t t e r = 4 α β 2

4.2. Simulation Results

4.2.1. Feasibility Analysis

The feasibility of the proposed method was demonstrated when SCR was set to 10 dB . Figure 4 shows the radar images of Frames 1, 5, 10, 15, 20 and 25. The target exists in Frames 10, 15 and 20, and the positions of the target are marked by red rectangles in corresponding images. As shown in Figure 4, the target was difficult to distinguish from the clutter background.
Figure 5 illustrates the distribution of particles calculated by the proposed SRBE-PF-TBD method, where the red triangles show the true positions of the target. When there was no target in Frames 1 and 5, particles were distributed according to the saliency map of the corresponding frame. In Frames 10, 15 and 20, where the target existed, the particles clustered around the target. When the target did not exist in Frame 25, the particles were distributed according to the saliency map again.
The true target trajectory and the estimated trajectory are shown in Figure 6. In the first two frames where the target existed, the target positions were not well estimated. However, with the accumulation of frames, the SRBE-PF-TBD method tracked the target effectively, which verified the effectiveness of SRBE-PF-TBD method.

4.2.2. Performance Evaluation

To analyze the performance of the proposed method, this section first gives the definition of the evaluation criteria [28,38,39], and then prsents the performance evaluation of the proposed method.
  • Root Mean Squared Error (RMSE) of Estimated Position
    The RMSE is calculated by
    RMSE k = 1 M i = 1 M x ^ i k x i k 2 + y ^ i k y i k 2 k = K b , . . . , K d 1
    where M is the number of Monte Carlo experiments, and x ^ i k , y ^ i k and x i k , y i k denote the estimated position and true position of the target in the ith Monte Carlo experiment, respectively.
  • Average RMSE
    The average RMSE is calculated by
    RMSE ¯ = 1 K Δ k = K b K d 1 RMSE ¯ k
    where K Δ = K d K b .
  • Detection Probability
    P d = 1 M i = 1 M K n u m p ^ ε k = 1 | z 1 : k > P t h K Δ k = K b , . . . , K d 1
    where K n u m p ^ ε k = 1 | z 1 : k > P t h is the number of frames that satisfy the decision p ^ ε k = 1 | z 1 : k > P t h in the ith Monte Carlo experiment.
  • False-alarm Probability
    P f a = 1 M i = 1 M K n u m p ^ ε k = 1 | z 1 : k > P t h K K Δ k = 1 , . . . , K b 1 K d , . . . , K
  • Earliest Effective Detection Frame
    F e = 1 M i = 1 M F i
    where F i represents the frame number that the target is first detected in the ith Monte Carlo experiment.
Figure 7 shows the probability density functions of the K-distribution corresponding to five different sets of parameters, while Figure 8 demonstrates the performance of the proposed method under these different K-distribution parameters. Figure 8a–d, respectively, shows the detection probabilities, false-alarm probabilities, earliest effective detection frame and average RMSE of the proposed method with a fixed detection threshold under different SCRs. The results in Figure 8a,c,d show that, with the increase of SCR, the detection probability increased, while the earliest effective detection frame and the average RMSE declined. In the case of the same SCR, the proposed method performed better with larger shape parameters. As for the false-alarm probabilities in Figure 8b, the false-alarm probabilities basically remained consistent under different SCRs when the distribution parameters of the K-distribution were fixed. Meanwhile, the false alarm probability decreased when the shape parameter α became larger. These are reasonable because, when the shape parameter of the K-distribution increased, the heavy-tailed effect of K-distribution became weaker. One extreme case was that, when the shape parameter tended to infinity, the K-distribution degenerated into a Rayleigh distribution.
Then, the proposed method was compared with other two methods: the K-based method [28] and the ESIR method [21]. Figure 9, Figure 10 and Figure 11, respectively, show the probabilities of existence and RMSEs calculated by the three methods in cases of SCR = 8 dB, 12 dB and 16 dB. The shape parameter of K-distribution was set as α = 2 and the scale parameter aws set as β = 2 2 , thus P c l u t t e r = 1 .
In the results of probability of existence, the frames of true target existence are represented by the black dash line. From the results of probabilities of existence, the K-based method and the proposed method had better target existence estimation with higher SCR. In the case of the same SCR, the proposed method performed better than the K-based method. The ESIR method had the best estimation of target existence when there was a target. However, in the case of target non-existence, the ESIR method still gave a high probability of existence, representing a large false alarm. This is because the ESIR method is based on Rayleigh distribution, which will generate large false alarm when the clutter has the heavy-tailed feature.
In the results of RMSEs, the RMSE bound is illustrated by the black dash line. The results show that the RMSEs of these methods decreased with the increase of SCR, which means the estimated positions of the three methods gradually converged to the true position of the target with the accumulation of target existence time. In comparison, the estimation accuracy of the proposed method was superior to the other two methods.
Figure 12 shows the detection probabilities of the three methods under different SCRs in the case the false-alarm probability was equal to 10 2 . In Figure 12a, the K-distribution parameters were set as α = 2 and β = 2 2 . In this situation, the K-based method performed best when the SCR was below 8 dB, while the proposed method became superior to the other two methods with the increasing SCR. Meanwhile, it is worth noting that the detection probabilities of the three method were all smaller than 0.1 when the SCR was lower than 8 dB, indicating the inability to achieve reliable target detection in practical applications. In Figure 12b, the K-distribution parameters were α = 10 and β = 2 2 , which means a smoother K-distributed clutter compared with Figure 12a. While the SCR was lower than 5 dB in Figure 12b, the detection probability of the ESIR method was slightly higher than the other two methods, which was attributed to the approximate Rayleigh distributed clutter in this case. When the SCR was beyond 5 dB, the proposed method had better detection performance than the other methods. Based on these results, it can be concluded that the proposed method can achieve better detection performance under more extensive clutter distribution parameters compared with the other two methods.

4.2.3. Computational Efficiency

Computational efficiency is a key indicator for real-time processing. To compare the computational efficiency of the three methods, all experiments were done in MATLAB R2018a, while the computer CPU was Intel (R) Core (TM) i5-7500 3.40 GHz, with 16 GB RAM.
Each frame contained 60 × 60 bins, the number of Monte Carlo experiments was 300, and the average running time of single frame was calculated by
t ¯ f = 1 300 i = 1 300 t f
where t f is the running time of single frame.
The computational consumption results are given in Table 1, which shows the average running time of the three methods with different numbers of particles. The results demonstrate that the proposed method and the ESIR method were obviously superior to the K-based method in term of running time. The K-based method calculated the likelihood ratio by a numerical solution of integral operation, causing a huge computational consumption. The computational efficiency of ESIR method was slightly higher than the proposed method, but, as shown in previous experimental results, the ESIR method only worked under the Rayleigh distribution clutter and caused high false alarm under the non-Gaussian clutter background.
All these experimental results demonstrate that the proposed method had fine detection and tracking performance and high computational efficiency.

5. Conclusions

In this paper, an effective maritime moving target detection and tracking method motivated by the human vision system is presented for airborne forward-looking scanning radar. The key idea of the proposed method is to use the saliency information of radar image to construct the SRBE likelihood ratio, which is independent from the prior knowledge of sea clutter. The proposed method effectively deals with the difficulty in target detection and tracking caused by non-Gaussian and non-stationary sea clutter. The performance of the proposed method was verified by Monte Carlo experiments using K-distribution clutter and Swerling type I target model. In the simulation experiments, the proposed method was evaluated under different K-distribution clutter parameters and compared with other two methods. The experimental results demonstrate that the proposed method had superior performance of target detection and tracking. Besides, the proposed method was computationally efficient and suitable for real-time application.
There are several continuing areas of research for improving the proposed method. The proposed method has been verified by simulated data, and future research will be to apply the proposed method to measured data. In addition, the extension of single-target situation to multi-target situation is the field to be researched.

Author Contributions

W.H. proposed the idea of the method; J.P. and W.H. conceived and designed the experiments; Q.Z. and W.H. performed the experiments; Y.H. and J.Y. analyzed the data; J.Y. contributed materials and analysis tools; and W.H. wrote the paper. All authors have approved the content of the submitted manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 61671117, in part by the Collaborative Innovation Center of Information Sensing and Understanding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jung, C.Y.; Yoo, S.L. Optimal Rescue Ship Locations Using Image Processing and Clustering. Symmetry 2019, 11, 32. [Google Scholar] [CrossRef]
  2. Surís-Regueiro, J.C.; Garza-Gil, M.D.; Varela-Lafuente, M.M. Marine economy: A proposal for its definition in the European Union. Mar. Policy 2013, 42, 111–124. [Google Scholar] [CrossRef]
  3. Dwyer, L. Emerging ocean industries: Implications for sustainable tourism development. Tour. Mar. Environ. 2018, 13, 25–40. [Google Scholar] [CrossRef]
  4. Rahmes, M.; Chester, D.; Hunt, J.; Chiasson, B. Optimizing cooperative cognitive search and rescue UAVs. Proc. SPIE 2018, 10643, 106430T. [Google Scholar]
  5. Varlamis, I.; Tserpes, K.; Sardianos, C. Detecting Search and Rescue missions from AIS data. In Proceedings of the 2018 IEEE 34th International Conference on Data Engineering Workshops (ICDEW), Paris, France, 16–20 April 2018; pp. 60–65. [Google Scholar]
  6. Tan, K.; Li, W.; Pei, J.; Huang, Y.; Yang, J. An I/Q-Channel Modeling Maximum Likelihood Super-Resolution Imaging Method for Forward-Looking Scanning Radar. IEEE Geosci. Remote Sens. Lett. 2018, 15, 863–867. [Google Scholar] [CrossRef]
  7. Zha, Y.; Huang, Y.; Sun, Z.; Wang, Y.; Yang, J. Bayesian deconvolution for angular super-resolution in forward-looking scanning radar. Sensors 2015, 15, 6924–6946. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Jakobsson, A.; Zhang, Y.; Huang, Y.; Yang, J. Wideband Sparse Reconstruction for Scanning Radar. IEEE Trans. Geosci. Remote Sens. 2018, 56, 1–14. [Google Scholar] [CrossRef]
  9. Pei, J.; Huang, Y.; Huo, W.; Zhang, Y.; Yang, J.; Yeo, T.S. SAR Automatic Target Recognition Based on Multiview Deep Learning Framework. IEEE Trans. Geosci. Remote Sens. 2018, 56, 2196–2210. [Google Scholar] [CrossRef]
  10. Kuttikkad, S.; Chellappa, R. Non-Gaussian CFAR Techniques for Target Detection in High Resolution SAR Images. In Proceedings of the 1st International Conference on Image Processing, Austin, TX, USA, 13–16 November 1994; pp. 910–914. [Google Scholar]
  11. Aalo, V.A.; Peppas, K.P.; Efthymoglou, G. Performance of CA-CFAR detectors in nonhomogeneous positive alpha-stable clutter. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 2027–2038. [Google Scholar] [CrossRef]
  12. El-Darymli, K.; McGuire, P.; Power, D.; Moloney, C.R. Target detection in synthetic aperture radar imagery: A state-of-the-art survey. J. Appl. Remote Sens. 2013, 7, 071598. [Google Scholar] [CrossRef]
  13. Blackman, S.; Popoli, R. Design and Analysis of Modern Tracking Systems (Artech House Radar Library); Artech House: Boston, MA, USA, 1999. [Google Scholar]
  14. Ramachandra, K. Kalman Filtering Techniques for Radar Tracking; Marcel Dekker: New York, NY, USA, 2000. [Google Scholar]
  15. Hadzagic, M.; Michalska, H.; Lefebvre, E. Track-before detect methods in tracking low-observable targets: A survey. Sens. Trans. Mag. 2005, 54, 374–380. [Google Scholar]
  16. Davey, S.J.; Rutten, M.G.; Gordon, N.J. Track-before-detect techniques. In Integrated Tracking, Classification, and Sensor Management; John Wiley & Sons: Hoboken, NJ, USA, 2013; pp. 311–362. [Google Scholar]
  17. Arnold, J.; Shaw, S.; Pasternack, H. Efficient target tracking using dynamic programming. IEEE Trans. Aerosp. Electron. Syst. 1993, 29, 44–56. [Google Scholar] [CrossRef]
  18. Moyer, L.R.; Spak, J.; Lamanna, P. A multi-dimensional Hough transform-based track-before-detect technique for detecting weak targets in strong clutter backgrounds. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 3062–3068. [Google Scholar] [CrossRef]
  19. Davey, S.J.; Rutten, M.G. A comparison of three algorithms for tracking dim targets. In Proceedings of the Information, Decision and Control, Adelaide, Australia, 12–14 February 2007; pp. 342–347. [Google Scholar]
  20. Boers, Y.; Driessen, J. Particle filter based detection for tracking. In Proceedings of the 2001 American Control Conference, Arlington, VA, USA, 25–27 June 2001; Volume 6, pp. 4393–4397. [Google Scholar]
  21. Rutten, M.G.; Gordon, N.J.; Maskell, S. Efficient particle-based track-before-detect in Rayleigh noise. In Proceedings of the International Conference on Information Fusion, Stockholm, Sweden, 28 June–1 July 2004; pp. 693–700. [Google Scholar]
  22. Rutten, M.G.; Ristic, B.; Gordon, N.J. A comparison of particle filters for recursive track-before-detect. In Proceedings of the 2005 8th International Conference on Information Fusion, Philadelphia, PA, USA, 25–28 July 2005; Volume 1, p. 7. [Google Scholar]
  23. Su, H.; Shui, P.; Liu, H.; Bao, Z. Particle filter based track-before-detect algorithm for over-the-horizon radar target detection and tracking. Chin. J. Electron. 2009, 18, 59–64. [Google Scholar]
  24. Lehmann, F. Recursive Bayesian filtering for multitarget track-before-detect in passive radars. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 2458–2480. [Google Scholar] [CrossRef]
  25. Gao, H.; Li, J. Detection and tracking of a moving target using SAR images with the particle filter-based track-before-detect algorithm. Sensors 2014, 14, 10829–10845. [Google Scholar] [CrossRef]
  26. Huang, D.; Xue, A.; Guo, Y. A particle filter track-before-detect algorithm for multi-radar system. Elektronika ir Elektrotechnika 2013, 19, 3–8. [Google Scholar] [CrossRef]
  27. Brekke, E.; Hallingstad, O.; Glattetre, J. Tracking small targets in heavy-tailed clutter using amplitude information. IEEE J. Ocean. Eng. 2010, 35, 314–329. [Google Scholar] [CrossRef]
  28. Jiang, H.; Yi, W.; Cui, G.; Kong, L.; Yang, X. Tracking targets in K clutter via particle filter. In Proceedings of the 2015 IEEE Radar Conference (RadarCon), Arlington, VA, USA, 10–15 May 2015; pp. 0350–0355. [Google Scholar]
  29. Ward, K.; Tough, R.; Watts, S. Sea Clutter: Scattering, the K Distribution and Radar Performance; Institution of Engineering and Technology: London, UK, 2013. [Google Scholar]
  30. Itti, L.; Koch, C.; Niebur, E. A model of saliency-based visual attention for rapid scene analysis. IEEE Trans. Pattern Anal. Mach. Intell. 1998, 20, 1254–1259. [Google Scholar] [CrossRef] [Green Version]
  31. Borji, A.; Itti, L. State-of-the-art in visual attention modeling. IEEE Trans. Pattern Anal. Mach. Intell. 2013, 35, 185–207. [Google Scholar] [CrossRef] [PubMed]
  32. Pei, J.; Huang, Y.; Huo, W.; Miao, Y.; Zhang, Y.; Yang, J. Synthetic Aperture Radar Processing Approach for Simultaneous Target Detection and Image Formation. Sensors 2018, 18, 3377. [Google Scholar] [CrossRef]
  33. Hou, X.; Zhang, L. Saliency detection: A spectral residual approach. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, 2007 (CVPR’07), Minneapolis, MN, USA, 17–22 June 2007; pp. 1–8. [Google Scholar]
  34. Piwowarski, R. Dynamic Waveform Design for Track-Before-Detect Algorithms in Radar. Master’s Thesis, Arizona State University, Phoenix, AZ, USA, 2011. [Google Scholar]
  35. Zhang, Y.; Zhang, Y.; Li, W.; Huang, Y.; Yang, J. Super-resolution surface mapping for scanning radar: Inverse filtering based on the fast iterative adaptive approach. IEEE Trans. Geosci. Remote Sens. 2018, 56, 127–144. [Google Scholar] [CrossRef]
  36. Pitt, M.K.; Shephard, N. Filtering via simulation: Auxiliary particle filters. Econ. Pap. 1999, 94, 590–599. [Google Scholar] [CrossRef]
  37. Richards; Mark, A. Fundamentals of Radar Signal Processing; Tata McGraw-Hill Education: New York, NY, USA, 2005. [Google Scholar]
  38. Davey, S.J.; Rutten, M.G.; Cheung, B. A Comparison of Detection Performance for Several Track-before-Detect Algorithms. Eur. J. Adv. Signal Process. 2007, 2008, 1–10. [Google Scholar] [CrossRef]
  39. Guo, S.; Cao, Z.; Feng, C.; Gao, S. Track before detect algorithm based on particle filter. In Proceedings of the 2011 Cross Strait Quad-Regional Radio Science and Wireless Technology Conference, Harbin, China, 26–30 July 2011; Volume 2, pp. 1590–1593. [Google Scholar] [CrossRef]
Figure 1. Diagram of state transition.
Figure 1. Diagram of state transition.
Sensors 19 01586 g001
Figure 2. Geometry illustration for forward-looking scanning radar and the flowchart of imaging processing.
Figure 2. Geometry illustration for forward-looking scanning radar and the flowchart of imaging processing.
Sensors 19 01586 g002
Figure 3. Flowchart of SRBE likelihood ratio calculation.
Figure 3. Flowchart of SRBE likelihood ratio calculation.
Sensors 19 01586 g003
Figure 4. Sampled images at different frames.
Figure 4. Sampled images at different frames.
Sensors 19 01586 g004
Figure 5. Distribution of particles corresponding with frames in Figure 4.
Figure 5. Distribution of particles corresponding with frames in Figure 4.
Sensors 19 01586 g005
Figure 6. The estimation results of the proposed method.
Figure 6. The estimation results of the proposed method.
Sensors 19 01586 g006
Figure 7. Probability density functions of K-distribution with different shape and scale parameters.
Figure 7. Probability density functions of K-distribution with different shape and scale parameters.
Sensors 19 01586 g007
Figure 8. Performances under different K-distribution parameters.
Figure 8. Performances under different K-distribution parameters.
Sensors 19 01586 g008
Figure 9. Probabilities of existence and Results of RMSE in case of SCR = 8 dB.
Figure 9. Probabilities of existence and Results of RMSE in case of SCR = 8 dB.
Sensors 19 01586 g009
Figure 10. Probabilities of existence and Results of RMSE in case of SCR = 12 dB.
Figure 10. Probabilities of existence and Results of RMSE in case of SCR = 12 dB.
Sensors 19 01586 g010
Figure 11. Probabilities of existence and Results of RMSE in case of SCR = 16 dB.
Figure 11. Probabilities of existence and Results of RMSE in case of SCR = 16 dB.
Sensors 19 01586 g011
Figure 12. Detection probabilities of the three methods in case of Pfa = 10−2.
Figure 12. Detection probabilities of the three methods in case of Pfa = 10−2.
Sensors 19 01586 g012
Table 1. Computational consumption.
Table 1. Computational consumption.
Number of ParticlesProposed MethodESIR MethodK-Based Method
40005.57 ms5.51 ms2.19 s
60008.79 ms7.20 ms3.27 s
800011.64 ms9.06 ms4.39 s
10,00015.23 ms11.46 ms5.4 s

Share and Cite

MDPI and ACS Style

Huo, W.; Pei, J.; Huang, Y.; Zhang, Q.; Yang, J. A New Maritime Moving Target Detection and Tracking Method for Airborne Forward-looking Scanning Radar. Sensors 2019, 19, 1586. https://doi.org/10.3390/s19071586

AMA Style

Huo W, Pei J, Huang Y, Zhang Q, Yang J. A New Maritime Moving Target Detection and Tracking Method for Airborne Forward-looking Scanning Radar. Sensors. 2019; 19(7):1586. https://doi.org/10.3390/s19071586

Chicago/Turabian Style

Huo, Weibo, Jifang Pei, Yulin Huang, Qian Zhang, and Jianyu Yang. 2019. "A New Maritime Moving Target Detection and Tracking Method for Airborne Forward-looking Scanning Radar" Sensors 19, no. 7: 1586. https://doi.org/10.3390/s19071586

APA Style

Huo, W., Pei, J., Huang, Y., Zhang, Q., & Yang, J. (2019). A New Maritime Moving Target Detection and Tracking Method for Airborne Forward-looking Scanning Radar. Sensors, 19(7), 1586. https://doi.org/10.3390/s19071586

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