Next Article in Journal
Novel Deep-Water Tidal Meter for Offshore Aquaculture Infrastructures
Previous Article in Journal
Virtual Reality for Shoulder Rehabilitation: Accuracy Evaluation of Oculus Quest 2
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tracking Multiple Targets Using Bearing-Only Measurements in Underwater Noisy Environments

Electronic and Electrical Department, Sungkyunkwan University, Suwon 03063, Korea
Sensors 2022, 22(15), 5512; https://doi.org/10.3390/s22155512
Submission received: 27 June 2022 / Revised: 17 July 2022 / Accepted: 22 July 2022 / Published: 24 July 2022
(This article belongs to the Section Intelligent Sensors)

Abstract

:
This article handles tracking multiple targets using bearing-only measurements in underwater noisy environments. For tracking multiple targets in underwater noisy environments, the Gaussian Mixture Probability Hypothesis Density (GM-PHD) filter provides good performance with its low computational load. Bearing-only measurements are passive and do not provide position information of a target. Note that the nonlinearity of the bearing-only measurements can be handled by Extended Kalman Filters (EKF) when applying the GM-PHD filter. However, range uncertainty of the target is large for bearing-only measurements. Thus, a single EKF leads to poor performance when it is applied in the GM-PHD. In this article, every bearing measurement gives birth to multiple target samples, which are distributed considering the feasible range of the passive sensor. Thereafter, every target sample is updated utilizing the measurement update step of the EKF. In this way, we run multiple EKFs associated to multiple target samples, instead of running a single EKF. To the best of our knowledge, our article is novel in tracking multiple targets in noisy environments, using the observer with bearing-only measurements. The effectiveness of the proposed GM-PHD is verified utilizing MATLAB simulations.

1. Introduction

Multiple Target Tracking (MTT) is applied as multiple targets are tracked with a camera [1] or sonar sensors [2,3]. The goal of MTT is to conjecture the target information from a series of measurement sets.
In this article, we consider a single observer, which has a passive sonar (receiver) array for measuring the bearing of sound emitted from a target. Array signal processing algorithms, such as the MUltiple SIgnal Classification (MUSIC) algorithm, have been used to estimate the bearing of sound generated from a target [4,5,6,7,8].
For underwater MTT, passive sonars are desirable, since they work in a stealthy manner and consume little energy. However, background noise (interference) may generate bearing measurements that do not originate from a true target. In addition, the birth and death of the target result in uncertain data association between target and bearing measurement. Therefore, data association in noisy environments is required for tackling underwater multiple target tracking problems.
Bearing-only measurements can be utilized for tracking a moving target [9,10,11], such as a target with a constant velocity [12]. Processing bearing-only measurements is a nonlinear estimation problem; thus, it requires various nonlinear filtering techniques for its solution [13,14,15,16,17].
For processing bearing-only measurements, one can apply the Extended Kalman filter (EKF), which is the nonlinear version of the Kalman filter linearization about an estimate of the current mean and covariance [11,14,18,19,20]. In the EKF, the state transition and observation models do not have to be linear functions of the state but may instead be differentiable functions. The EKF has both low computational complexity and reasonable estimation accuracy; therefore, it has been widely used for tracking based on bearing-only measurements [11,14,19,20]. In our paper, we also use the EKF for processing bearing-only measurements.
For handling an observer tracking a target using bearing-only measurements, the Range-Parameterized Extended Kalman Filter (RPEKF) was developed [19,20]. The RPEKF is a Gaussian-sum filter with multiple EKFs, each initialized at an estimated target range. In this way, the RPEKF reduces the initial range estimation error. The RPEKF assumes that the maximum sensing range of the observer is known in advance. Inspired by the RPEKF [19,20], our paper also assumes that the maximum sensing range of an observer is known in advance.
The RPEKF is suitable for tracking a single target using bearing-only measurements, and the RPEKF does not consider noisy environments. Distinct from the RPEKF, our paper considers tracking multiple targets using bearing-only measurements in noisy environments.
For tracking multiple targets in noisy environments, the Multiple Hypotheses Tracking (MHT) [21] and Joint Probability Data Association (JPDA) [22,23] filters were developed. However, the computational burden of these methods increases as the number of targets or interference points increases.
Recently, the Random Finite Set (RFS) was developed for avoiding explicit associations between measurement and target [3,24,25]. The Probability Hypothesis Density (PHD) filter [3,24,25] was invented using the RFS. The PHD is computationally effective, since it avoids data association between target and measurement. The authors of [25] invented the Gaussian Mixture Probability Hypothesis Density (GM-PHD) filter, which is time-efficient for solving the data association problem. For tracking multiple targets that maneuver in noisy environments, ref. [26] developed a multiple model GM-PHD filter. As far as we know, no paper on MTT has considered tracking multiple targets using bearing-only measurements.
In our article, the observer measures the bearing of sound emitted from a target in real time. Note that the nonlinearity of the bearing-only measurements can be handled by the EKF [18] when applying the GM-PHD filter. However, range uncertainty of the target is large for bearing-only measurements. This results in poor performance when a single EKF is utilized in the GM-PHD.
Therefore, in this article, every bearing measurement gives birth to multiple target samples, which are distributed considering the feasible range of the passive sensor. Note that the feasible range of the passive sensor is available since we assume that the maximum sensing range of an observer is known in advance. Thereafter, every target sample is updated using the measurement update step of the EKF. In this way, we run multiple EKFs associated to multiple target samples instead of running a single EKF.
As far as we know, this paper is novel in tracking multiple targets in noisy environments, as the observer uses bearing-only measurements. The effectiveness of the proposed GM-PHD is verified utilizing MATLAB simulations.
The remainder of this article is organized as follows. Section 2 introduces the bearing-only target tracking as the preliminary information of this article. Section 3 addresses the GM-PHD filter for bearing-only target tracking. Section 4 discusses MATLAB simulations for demonstrating the performance of the proposed filter. Conclusions are addressed in Section 5.

2. Preliminary Information

We introduce the matrix notations utilized in this paper. Let C [ i ] define the i-th element in a column vector C . Let I n define the identity matrix with n rows and n columns. Let 0 n define the zero matrix with n rows and n columns. d i a g ( a , b , c , ) defines a diagonal matrix having a , b , c , as its diagonal elements in this order.

Bearing-Only Target Tracking

Let x k define the target state at sampling-stamp k. The target’s location and velocity at sampling-stamp k are defined as
x k = [ x k t , y k t , x ˙ k t , y ˙ k t ] .
Suppose that the maximum speed of a target is known a priori. Let v m a x define the target’s maximum speed.
Let T define the sampling interval. Let [ a x , a y ] indicate the target’s acceleration, which is not known in advance. Considering a target with Constant Velocity (CV) motion, its process model is
x k + 1 = F x k + μ k ,
where the state transition matrix F is
F = 1 T 0 1 I 2 .
In (2), μ k defines the process noise with mean 0 and variance Q , where
Q = Γ d i a g ( σ a x 2 , σ a y 2 ) Γ .
Here, d i a g ( σ a x 2 , σ a y 2 ) is a diagonal matrix having σ a x 2 , σ a y 2 as its diagonal elements in this order. Here, σ a presents the standard deviation for a. In addition, (4) utilizes
Γ = T 2 2 I 2 T I 2 .
Let [ x k o , y k o ] define the 2D position of the observer at sampling-stamp k. In addition, [ x ˙ k o , y ˙ k o ] defines the 2D velocity of the observer at sampling-stamp k. In bearing-only target tracking, observer maneuvering is required for target localization [27,28]. This implies that the observer changes its velocity [ x ˙ k o , y ˙ k o ] in order to obtain the observability on the target position.
Let r m a x define the maximum sensing range of the observer’s sonar. In addition, the sonar’s minimum sensing range is r m i n . Inspired by the RPEKF [19,20], we assume that both r m i n and r m a x are known in advance.
Once the relative distance between the target and the observer is within the interval [ 0 , r m a x ] , the bearing of the target can be measured by the observer. The bearing measurement at sampling-stamp k is
b ( x k ) = a t a n 2 ( y k t y k o , x k t x k o ) + e k .
Here, a t a n 2 ( y , x ) is the phase (angle) of a complex number x + j y . In addition, e k defines the measurement noise having a Gaussian distribution with mean 0 and standard deviation R .
Ignoring the noise term in (6), we obtain
x k t = ( x k t x k o , y k t y k o ) cos ( b ( x k ) ) + x k o , y k t = ( x k t x k o , y k t y k o ) sin ( b ( x k ) ) + y k o .

3. GM-PHD Filter

We utilize the GM-PHD filter in [25]. The multitarget state space is defined as F ( X ) . At sampling-stamp k, target sample set X k and bearing measurement set Z k are defined as
X k = { x k , 1 , x k , 2 , , x k , N ( k ) } F ( X ) .
Z k = { z k , 1 , z k , 2 , , z k , M ( k ) } .
Here, N ( k ) defines the number of target samples and M ( k ) defines the number of bearing measurements at sampling-stamp k.
Every x k X k defines the target pose (position and velocity) addressed in (8). For a given state X k 1 at sampling-stamp k 1 , every x k 1 X k 1 survives at sampling-stamp k with probability P S ( x k 1 ) . Therefore, the survived target is represented using a RFS S k | k 1 ( x k 1 ) .
At sampling-stamp k, a new target appears by spontaneous birth or by spawning from an existing target at sampling-stamp k 1 . In addition, spontaneous birth sets  Γ k and spawned target sets B k | k 1 ( x k 1 ) are modeled as RFS at sampling-stamp k. Birth and spawn of a new target from bearing measurements are introduced in Section 3.
At sampling-stamp k, X k is
X k = ζ X k 1 S k | k 1 ( ζ ) ( ζ X k 1 B k | k 1 ( ζ ) ) Γ k .
In addition, at sampling-stamp k, a target is detected by bearing sensors with probability P D ( x k ) . Every target x k X k generates its associated bearing measurement b ( x k ) at sampling-stamp k. See (6) for the equation of b ( x k ) . In addition, the bearing sensor measures an interference point at sampling-stamp k. An interference point is modeled as a RFS K k . Therefore, the RFS measurement set Z k is
Z k = K k ( x X k b ( x ) ) .
Let p k ( X k | Z 1 : k ) define the posterior density, f k | k 1 define the transition density, and g k ( Z k | X ) define the likelihood of every target. The posterior is calculated as
p k | k 1 ( X k | Z 1 : k 1 ) = f k | k 1 ( X k | X ) p k 1 ( X | Z 1 : k 1 ) μ ( d X ) .
p k ( X k | Z 1 : k ) = g k ( Z k | X k ) p k | k 1 ( X k | Z 1 : k 1 ) g k ( Z k | X ) p k | k 1 ( X | Z 1 : k 1 ) μ ( d X ) .
Here, μ is an appropriate reference measure on F ( X ) [29].
It is assumed that an interference is independent of target-originated measurements. In addition, it is assumed that an interference point and predicted RFS are distributed using a Poisson distribution.
The survival and detection probabilities do not depend on the state vector; hence, we utilize P S , k and P D , k instead of P S ( x k ) and P D ( x k ) , respectively.
Let N ( x ; m , P ) denote the Gaussian density with mean m and covariance P. The birth intensity function of a newborn target at sampling-stamp k is
γ k ( x ) = i = 1 J γ , k w γ , k ( i ) N ( x ; m γ , k ( i ) , P γ , k ( i ) )
where J γ , k represents the number of newborn Gaussian components, and w γ , k ( i ) , m γ , k ( i ) , and P γ , k ( i ) are the weight, mean, and covariance, respectively, of the i-th newborn Gaussian component.
Let v k | k define the posterior density intensity at sampling-stamp k. The posterior intensity propagates utilizing
v k | k 1 ( x ) = P S , k f k | k 1 ( x | ζ ) v k 1 ( ζ ) d ζ + γ k ( x ) .
v k | k ( x ) = ( 1 P D , k ) v k | k 1 ( x ) + z Z k P D , k g k ( z | x ) v k | k 1 ( x ) λ k c k ( z ) + P D , k g k ( z | ζ ) v k | k 1 ( ζ ) d ζ .
Here, λ k is the average number of interference points at sampling-stamp k. In addition, c k ( z ) defines the probability distribution of every interference point.
In the proposed GM-PHD, the prediction step is linear and the measurement update step is nonlinear. Thus, we only have to linearize the measurement update step. In our paper, the EKF is used for the measurement update step, since the EKF has both low computational complexity and reasonable estimation accuracy [11,14,19,20].
We perform the measurement update, once we obtain the newest measurement set Z k at sampling-stamp k. Since the measurement model in (6) is not linear, we can utilize the EKF in [18] for the measurement update step.
Under the GM-PHD [25,30], (15) is replaced by
v k | k 1 ( x ) = i = 1 J k | k 1 w k | k 1 ( i ) N ( x ; m k | k 1 ( i ) , P k | k 1 ( i ) ) = v S , k | k 1 ( x ) + γ k ( x )
where
v S , k | k 1 ( x ) = i = 1 J k 1 P S , k w k 1 ( i ) N ( x ; m k | k 1 ( i ) , P k | k 1 ( i ) ) .
Here, J k 1 is the number of components of the intensity at sampling-stamp k 1 . Moreover, we utilize
m k | k 1 ( i ) = F m k 1 | k 1 ( i ) .
In addition, we utilize
P k | k 1 ( i ) = F P k 1 | k 1 ( i ) F + Q .
Further, (16) is replaced by
v k | k ( x ) = ( 1 P D , k ) v k | k 1 ( x ) + i = 1 J k | k 1 w k | k ( i ) ( z ) N ( x ; m k | k ( i ) , P k | k ( i ) ) .
In (21), we use
w k | k ( i ) ( z ) = P D , k g k ( z | m k | k 1 ( i ) ) w k | k 1 ( i ) λ k c k ( z ) + P D , k l = 1 J k | k 1 g k ( z | m k | k 1 ( l ) ) w k | k 1 ( l ) .
Here, one uses
g k ( z | m k | k 1 ( i ) ) = N ( z ; b ( m k | k 1 ( i ) ) , S k ( i ) ) ,
where
S k ( i ) = H P k | k 1 ( i ) H + R .
In (21), we use
m k | k ( i ) = m k | k 1 ( i ) + K ( z k , g b ( m k | k 1 ( i ) ) ) .
Here, K is
K = P k | k 1 ( i ) H ( S k ( i ) ) 1 .
In addition, H = b ( x ) x where x = m k | k 1 ( i ) .
In (21), we further use
P k | k ( i ) = ( I K H ) P k | k 1 ( i ) ,
where I is the identity matrix. The detailed derivation of (21) is addressed in [25,30].
To reduce the computational load, Gaussian components must be pruned and merged [30]. Then, the posterior intensity is represented as
v k | k ( x ) = i = 1 J k w k ( i ) N ( x ; m k | k ( i ) , P k | k ( i ) ) .
Here, w k ( i ) is the weight of Gaussian distribution and J k is the number of components of the intensity at sampling-stamp k.

Birth or Spawn of New Target Samples Based on Bearing Measurements

Recall that at every sampling-stamp k, a new target sample is generated by spontaneous birth or by spawning from an existing target sample at sampling-stamp k 1 . Spontaneous birth sets Γ k and spawned target sets B k | k 1 ( x k 1 ) are modeled as RFS at sampling-stamp k.
In this article, the observer measures the bearing of sound emitted from a target in real time. Bearing-only measurements are passive and do not provide position information of a target. Therefore, in this article, every bearing measurement gives birth to multiple target samples, which are distributed considering the feasible range of the passive sensor.
We update the mean m and covariance P of the newly born targets in Γ k as follows. Suppose that we generate B target samples using z k , g ( g { 1 , 2 , , M ( k ) } ), the g-th bearing measurement at sampling-stamp k. Let t b ( b { 1 , 2 , , B } ) define the b-th target sample generated on the bearing line at sampling-stamp k.
The feasible range of the true target is in the interval [ r m i n , r m a x ] . This interval is divided into B subintervals [ r m i n + ( b 1 ) s R , r m i n + b s R ] , where s R = r m a x r m i n B ( b { 1 , 2 , , B } ).
The range for the center of every subinterval is r b = r m i n + ( b 0.5 ) s R ( b { 1 , 2 , , B } ). In addition, the velocity of every newly born target sample is zero. Under (7), the mean m of a newly born target sample t b ( b { 1 , 2 , , B } ) is
m ( t b ) = r b cos ( z k , g ) + x k o r b sin ( z k , g ) + y k o 0 0
where r b = r m i n + ( b 0.5 ) s R .
We next derive the covariance P of t b . Under (6), the variance of the bearing error is R 2 . In addition, since every subinterval length is s R , the variance of the range error is s R 2 4 .
The error covariance related to the 2D position of t b ( b { 1 , 2 , , B } ) is
P s ( t b ) = D d i a g ( s R 2 4 , R 2 ) D
where D is
D = m ( t b ) [ 1 ] r b , m ( t b ) [ 1 ] z k , g m ( t b ) [ 2 ] r b , m ( t b ) [ 2 ] z k , g .
Under (30), we further obtain
D = cos ( z k , g ) , r b sin ( z k , g ) sin ( z k , g ) , r b cos ( z k , g ) .
The covariance P of t b ( b { 1 , 2 , , B } ) is derived as
P ( t b ) = P s ( t b ) , 0 2 0 2 , d i a g ( v m a x 2 , v m a x 2 ) .
In the case where the relative distance between a newly born target sample, e.g., t b ( b { 1 , 2 , , B } ), and an existing target sample, e.g., x k 1 , is less than a certain threshold α > 0 , we say that the born target sample t b is spawned from the existing target sample x k 1 . This implies that we have t b B k | k 1 ( x k 1 ) .
We update the mean m and covariance P of the newly spawned target sample t b B k | k 1 ( x k 1 ) as follows. We derive the covariance P of t b using (33). In addition, we use
m ( t b ) = x k 1 [ 1 ] x k 1 [ 2 ] 0 0 .

4. MATLAB Simulations

We address MATLAB simulation results for verifying the effectiveness of the proposed GM-PHD filter. Target bearings are measured at every d t = 1 s. In addition, the simulation running time is 500 s.
Target detection probability is 0.98. The target’s maximum speed is v m a x = 30 m/s. Considering the sonar sensing range, we set r m a x = 10,000 m. In addition, the sonar’s minimum sensing range is set to r m i n = 500 m. The standard deviation of bearing measurement noise is 1 degree.
In bearing-only target tracking, observer maneuvering is required for improving the tracking convergence [27,28]. Since the observer uses bearing-only measurements, the tracking convergence depends on the bearing rate of a target with respect to the observer. In other words, the geometry of a target with respect to the observer determines the tracking convergence of the target.
In the simulations, the observer maneuvers as follows. The observer moves with a constant speed v o b s = 10 m/s, and it starts from the origin. Before 125 s, the observer moves with course 0 degree measured clockwise from the east direction. From 125 to 250 s, the observer moves with course 135 degrees measured clockwise from the east direction. From 250 to 375 s, the observer moves with course 0 degree measured clockwise from the east direction. After 375 s, the observer moves with course 90 degrees measured clockwise from the east direction.
See Figure 1 for the simulation scenario. In this figure, the observer’s path is marked with yellow colors. Every target moves with a constant velocity. In Figure 1, we present the direction of target motion using symbols. A circle symbol indicates the start point of a target motion and a diamond symbol indicates the end point of a target motion.
Since the sonar’s maximum sensing range is 10 km, we deploy a target initially so that it can be detected by the observer. The position of each target at every second is marked with distinct colors. Target 1 (red) starts from [7500, 7500] and moves with velocity [3, −5] in m/s. Target 2 (blue) starts from [−2500, −2500] and moves with velocity [0, 5] in m/s. Targets 1 (red) and 2 (blue) are generated at sampling-stamp 0, and target 3 (green) is spawned from target 1 at 30 s. Target 3 (green) moves with velocity [−5, −2] in m/s.
For clear presentation of the observer’s path, Figure 1 simulates the case where there is no interference point. In other words, Figure 1 shows the simulation result without interference points. In this figure, the estimated target position is marked with magenta triangles. Before the observer maneuver at 125 s, true target positions are not tracked by target estimations (magenta triangles). However, as time goes on, true target positions are tracked by target estimations.
Considering the simulation of the scenario in Figure 1, Figure 2 shows the simulation result with interference points. Recall that λ k (16) is the average number of interference points at sampling-stamp k. We set λ k = 1 , which implies that, on average, one interference point is generated at each sampling-stamp. In Figure 2, the estimated target position is marked with magenta triangles. Due to interference points, many false targets are generated surrounding the observer. Before the observer maneuver at 125 s, true target positions are not tracked by target estimations (magenta triangles). However, as time goes on, true target positions are tracked by target estimations.
Considering the simulation of Figure 2, Figure 3 shows the bearing (in radians) change with respect to time (in seconds). In this figure, true target bearing measurements are marked with circles. In addition, interference points are marked with crosses. The estimated target bearing is marked with magenta circles. Before the observer maneuver at 125 s, true target bearings are not tracked by the estimated target bearing (magenta circles). However, as time goes on, true target bearings are tracked by estimated target bearings.

4.1. Monte-Carlo (MC) Simulations

A total of 100 Monte-Carlo (MC) simulations are utilized for verifying the outperformance of the proposed filter. Let M C = 100 denote the number of total MC simulations. We set λ k = 1 , which implies that, on average, one interference point is generated at each sampling-stamp. We generate B = 50 target samples using every bearing measurement at every sampling-stamp. It takes 14 s to run a single MC simulation.
Tracking accuracy is checked with the cardinality (number of targets) and the Optimal Subpattern Assignment (OSPA) metric [31,32]. We derive the cardinality as the mean of the estimated number of targets under M C = 100 MC simulations. We also derive the mean of OSPA using M C = 100 MC simulations.
Recall that J k in (28) is the number of components of the intensity. The cardinality at sampling-stamp k is defined as
c a r d k = j = 1 M C J k M C .
The true cardinality at sampling-stamp k indicates the true number of targets detected at sampling-stamp k.
We briefly address the Optimal Subpattern Assignment (OSPA) metric in [31,32]. Let d ( c ) ( x , y ) = m i n ( c , x y ) denote the distance between x and y . Here, c denotes a cut-off parameter in [31,32]. In our simulations, we use c = 2 10 4 .
Let M ( k ) = { m 1 ( k ) , m 2 ( k ) , , m J k ( k ) } denote the set of target position estimates at sample-step k. Recall that J k denotes the number of components of the intensity. Suppose that there are N ( k ) targets at sample-step k. Let Tgt ( k ) = { tgt 1 ( k ) , tgt 2 ( k ) , , tgt N ( k ) ( k ) } denote the true target positions at sample-step k. The OSPA distance between the target estimates M ( k ) and the true target positions Tgt ( k ) is defined as
o s p a d ( k ) = 1 N ( k ) ( m i n q P N ( k ) i = 1 J k d ( c ) ( m i ( k ) , tgt q ( i ) ( k ) ) + c ( N ( k ) J k ) ) .
Here, P N ( k ) denotes the set of permutations on { 1 , 2 , , N ( k ) } . The OSPA metric can be computed efficiently by using the Hungarian method for optimal point assignment [31,32].
The mean of OSPA at sample-step k is
o s p a k = j = 1 M C o s p a d ( k ) M C .
See Figure 4 for the OSPA o s p a k and the cardinality c a r d k of the proposed GM-PHD. Before 30 s pass, there are two targets (Targets 1 and 2). Target 3 is spawned from target 1 at 30 s. In bearing-only target tracking, observer maneuver is required for target localization [27,28]. The observer maneuvers at 125 s. See that the OSPA decreases gradually after the observer maneuver. In Figure 4, T r u e plots the number of bearing measurements associated with true targets. In addition, E s t i m a t e is the mean of the estimated number of targets using 100 MC simulations. As time goes on, OSPA gradually approaches to zero and E s t i m a t e converges to T r u e in cardinality. In other words, the cardinality error converges to almost zero as time goes on.

4.1.1. Comparison with the Case Where a Single EKF Is Used in the GM-PHD

Note that the nonlinearity of the bearing-only measurements can be handled by the EKF when applying the GM-PHD filter. However, range uncertainty of the target is large for bearing-only measurements ( r m a x = 10 km in MATLAB simulations). This results in poor performance when a single EKF is used in the GM-PHD. Setting B = 1 corresponds to the case where we apply a single EKF in the GM-PHD. As we increase B, we distribute more target samples considering the maximum sensing range of the bearing sensor. Consider the scenario in Figure 2. We next check the effect of changing B, the number of target samples generated using every bearing measurement at every sampling-stamp.
Considering the case where B = 1 is utilized, Figure 5 plots the OSPA and the cardinality of the proposed GM-PHD. Compared with the case where B = 50 is used (Figure 4), the estimation accuracy decreases considerably. This shows the effectiveness of the proposed sample distribution approach used in our paper. It takes 6 s to run a single MC simulation, as we utilize B = 1 . Compared with the case where B = 50 is used (Figure 4), the computational load decreases.

4.1.2. The Effect of Changing B

Considering the case where B = 10 is utilized, Figure 6 plots the OSPA o s p a k and the cardinality c a r d k of the proposed GM-PHD. Compared with the case where B = 1 is used (Figure 5), the estimation accuracy increases considerably. It takes 6 s to run a single MC simulation. Compared with the case where B = 50 is used, the computational load decreases.
Considering the case where B = 20 is utilized, Figure 7 plots the OSPA and the cardinality of the proposed GM-PHD. Compared with the case where B = 1 is used (Figure 5), the estimation accuracy increases considerably. It takes 7 s to run a single MC simulation. Compared with the case where B = 50 is used, the computational load decreases.
Table 1 shows the computational load (running time for one MC simulation) as B varies. Recall that setting B = 1 corresponds to the case where we apply a single EKF in the GM-PHD. Table 1 shows that as B increases, the computational load increases in general. However, considering both the computational load and the estimation accuracy (Figure 6), setting B = 10 is desirable.

5. Conclusions

This article handles tracking multiple targets using bearing-only measurements in underwater noisy environments. We apply GM-PHD filters to solve our MTT problem based on bearing-only measurements. Since bearing-only measurements do not provide position information of a target, we make every bearing measurement give birth to B target samples, which are distributed considering the feasible range of the passive sensor. The effectiveness of the proposed GM-PHD is verified utilizing MATLAB simulations.
In the future, we will extend the proposed tracking filter so that one can handle the case where multiple observers are used for tracking multiple targets in noisy environments under bearing-only measurements. As we use multiple observers, false targets may appear as bearing lines of multiple observers’ intersects. Thus, we require filters for removing false targets effectively.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Luo, W.; Xing, J.; Milan, A.; Zhang, X.; Liu, W.; Kim, T.K. Multiple object tracking: A literature review. Artif. Intell. 2021, 293, 103448. [Google Scholar] [CrossRef]
  2. Ozer, E.; Akar, A.O.; Hocaoglu, A.K. Passive sonar multiple target tracking with different resampling algorithms. In Proceedings of the 2018 26th Signal Processing and Communications Applications Conference (SIU), Izmir, Turkey, 2–5 May 2018; pp. 1–4. [Google Scholar]
  3. Sheng, X.; Chen, Y.; Guo, L.; Yin, J.; Han, X. Multitarget Tracking Algorithm Using Multiple GMPHD Filter Data Fusion for Sonar Networks. Sensors 2018, 18, 3193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Mohanna, M.; Rabeh, M.L.; Zieur, E.M.; Hekala, S. Optimization of MUSIC algorithm for angle of arrival estimation in wireless communications. NRIAG J. Astron. Geophys. 2013, 2, 116–124. [Google Scholar] [CrossRef] [Green Version]
  5. Qian, C.; Huang, L.; So, H. Computationally efficient ESPRIT algorithm for direction-of-arrival estimation based on Nyström method. Signal Process. 2014, 94, 74–80. [Google Scholar] [CrossRef]
  6. Gupta, P.; Kar, S. MUSIC and improved MUSIC algorithm to estimate direction of arrival. In Proceedings of the 2015 International Conference on Communications and Signal Processing (ICCSP), Melmaruvathur, India, 2–4 April 2015; pp. 0757–0761. [Google Scholar]
  7. Kim, J. Direction of Arrival Estimation Using Four Isotropic Receivers. IEEE Instrum. Meas. Mag. 2021, 24, 77–81. [Google Scholar] [CrossRef]
  8. Wen, F.; Javed, U.; Yang, Y.; He, D.; Zhang, Y. Improved subspace direction-of-arrival estimation in unknown nonuniform noise fields. In Proceedings of the 2016 Fourth International Conference on Ubiquitous Positioning, Indoor Navigation and Location Based Services (UPINLBS), Shanghai, China, 2–4 November 2016; pp. 230–233. [Google Scholar]
  9. Kim, J. Obstacle Information Aided Target Tracking Algorithms for Angle-Only Tracking of a Highly Maneuverable Target in Three Dimensions. IET Radar Sonar Navig. 2019, 13, 1074–1080. [Google Scholar] [CrossRef]
  10. Kim, J. Observer manoeuvre control to track multiple targets considering Doppler-bearing measurements in threat environments. IET Radar Sonar Navig. 2019, 13, 2158–2165. [Google Scholar] [CrossRef]
  11. Kim, J. Maneuvering target tracking of underwater autonomous vehicles based on bearing-only measurements assisted by inequality constraints. Ocean Eng. 2019, 189, 106404. [Google Scholar] [CrossRef]
  12. Nardone, S.; Lindgren, A.; Gong, K. Fundamental properties and performance of conventional bearings-only target motion analysis. IEEE Trans. Autom. Control 1984, 29, 775–787. [Google Scholar] [CrossRef]
  13. Clark, J.; Vinter, R.; Yaqoob, M. The shifted Rayleigh filter for bearings only tracking. In Proceedings of the 2005 7th International Conference on Information Fusion, Philadelphia, PA, USA, 25–28 July 2005; Volume 1, p. 8. [Google Scholar]
  14. Kim, J.; Suh, T.; Ryu, J. Bearings-only target motion analysis of a highly manoeuvring target. IET Radar Sonar Navig. 2017, 11, 1011–1019. [Google Scholar] [CrossRef]
  15. Jiang, H.; Cai, Y. Bearings-only tracking with a Gaussian-sum based ensemble Kalman filter. In Proceedings of the 2017 29th Chinese Control And Decision Conference (CCDC), Chongqing, China, 28–30 May 2017; pp. 4823–4828. [Google Scholar]
  16. Li, L.; Xie, W.; Liu, Z. Bearings-only maneuvering target tracking based on truncated quadrature Kalman filtering. Int. J. Electron. Commun. 2015, 69, 281–289. [Google Scholar]
  17. Ristic, B.; Houssineau, J.; Arulampalam, S. Robust target motion analysis using the possibility particle filter. IET Radar Sonar Navig. 2019, 13, 18–22. [Google Scholar] [CrossRef]
  18. Ristic, B.; Arulampalam, S.; Gordon, N. Beyond the Kalman Filter: Particle Filters for Tracking Applications; Artech House: Norwood, MA, USA, 2004. [Google Scholar]
  19. Peach, N. Bearings-only tracking using a set of range-parameterised extended Kalman filters. IEE Proc. Control Theory Appl. 1995, 142, 73–80. [Google Scholar] [CrossRef]
  20. Karlsson, R.; Gustafsson, F. Recursive Bayesian estimation: Bearing-only applications. IEE Proc. Radar Sonar Navig. 2005, 152, 305–313. [Google Scholar] [CrossRef] [Green Version]
  21. Blackman, S. Multiple hypothesis tracking for multiple target tracking. IEEE Aerosp. Electron. Syst. Mag. 2004, 19, 5–18. [Google Scholar] [CrossRef]
  22. He, S.; Shin, H.S.; Tsourdos, A. Joint Probabilistic Data Association Filter with Unknown Detection Probability and Clutter Rate. Sensors 2018, 18, 269. [Google Scholar] [CrossRef] [Green Version]
  23. Musicki, D.; Kaune, R.; Koch, W. Mobile Emitter Geolocation and Tracking Using TDOA and FDOA Measurements. IEEE Trans. Signal Process. 2009, 58, 1863–1874. [Google Scholar] [CrossRef]
  24. Vo, B.N.; Singh, S.; Doucet, A. Sequential Monte Carlo methods for multitarget filtering with random finite sets. IEEE Trans. Aerosp. Electron. Syst. 2005, 41, 1224–1245. [Google Scholar]
  25. Vo, B.N.; Ma, W.K. The Gaussian Mixture Probability Hypothesis Density Filter. IEEE Trans. Signal Process. 2006, 54, 4091–4104. [Google Scholar] [CrossRef]
  26. Wang, X.; Han, C. An improved multiple model GM-PHD filter for maneuvering target tracking. Chin. J. Aeronaut. 2013, 26, 179–185. [Google Scholar] [CrossRef] [Green Version]
  27. Oshman, Y.; Davidson, P. Optimization of observer trajectories for bearings-only target localization. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 892–902. [Google Scholar] [CrossRef]
  28. Baek, S.S.; Kwon, H.; Yoder, J.A.; Pack, D. Optimal Path Planning of a Target-Following Fixed-Wing UAV Using Sequential Decision Processes. In Proceedings of the Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Tokyo, Japan, 3–7 November 2013; pp. 2959–2962.
  29. Mahler, R. Multitarget Bayes filtering via first-order multitarget moments. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 1152–1178. [Google Scholar] [CrossRef]
  30. Gao, Y.; Jiang, D.; Zhang, C.; Guo, S. A Labeled GM-PHD Filter for Explicitly Tracking Multiple Targets. Sensors 2021, 21, 3932. [Google Scholar] [CrossRef] [PubMed]
  31. Schuhmacher, D.; Vo, B.T.; Vo, B.N. A Consistent Metric for Performance Evaluation of Multi-Object Filters. IEEE Trans. Signal Process. 2008, 56, 3447–3457. [Google Scholar] [CrossRef] [Green Version]
  32. Nguyen, H.V.; Rezatofighi, H.; Vo, B.N.; Ranasinghe, D.C. Distributed Multi-Object Tracking Under Limited Field of View Sensors. IEEE Trans. Signal Process. 2021, 69, 5329–5344. [Google Scholar] [CrossRef]
Figure 1. For clear presentation of the observer’s path, this figure simulates the case where there is no interference point. The observer’s path is marked with yellow colors. Target 1 (red) starts from [7500, 7500] and moves with velocity [3, −5] in m/s. Target 2 (blue) starts from [−2500, −2500] and moves with velocity [0, 5] in m/s. Targets 1 (red) and 2 (blue) are generated at sampling-stamp 0, and target 3 (green) is spawned from target 1 at 30 s. Target 3 (green) moves with velocity [−5, −2] in m/s. We present the direction of target motion using symbols. A circle symbol indicates the start point of a target motion, and a diamond symbol indicates the end point of a target motion. The estimated target position is marked with magenta triangles.
Figure 1. For clear presentation of the observer’s path, this figure simulates the case where there is no interference point. The observer’s path is marked with yellow colors. Target 1 (red) starts from [7500, 7500] and moves with velocity [3, −5] in m/s. Target 2 (blue) starts from [−2500, −2500] and moves with velocity [0, 5] in m/s. Targets 1 (red) and 2 (blue) are generated at sampling-stamp 0, and target 3 (green) is spawned from target 1 at 30 s. Target 3 (green) moves with velocity [−5, −2] in m/s. We present the direction of target motion using symbols. A circle symbol indicates the start point of a target motion, and a diamond symbol indicates the end point of a target motion. The estimated target position is marked with magenta triangles.
Sensors 22 05512 g001
Figure 2. This figure shows the simulation result with interference points. We set λ k = 1 , which implies that, on average, one interference point is generated at each sampling-stamp. The estimated target position is marked with magenta triangles. Before the observer maneuver at 125 s, true target positions are not tracked by target estimations (magenta triangles). However, as time goes on, true target positions are tracked by target estimations.
Figure 2. This figure shows the simulation result with interference points. We set λ k = 1 , which implies that, on average, one interference point is generated at each sampling-stamp. The estimated target position is marked with magenta triangles. Before the observer maneuver at 125 s, true target positions are not tracked by target estimations (magenta triangles). However, as time goes on, true target positions are tracked by target estimations.
Sensors 22 05512 g002
Figure 3. The bearing (in radians) change with respect to time (in seconds). In this figure, true target bearings are marked with circles. In addition, interference points are marked with crosses. The estimated target bearing is marked with magenta circles. Before the observer maneuver at 125 s, true target bearings are not tracked by the estimated target bearing (magenta circles). However, as time goes on, true target bearings are tracked by estimated target bearings.
Figure 3. The bearing (in radians) change with respect to time (in seconds). In this figure, true target bearings are marked with circles. In addition, interference points are marked with crosses. The estimated target bearing is marked with magenta circles. Before the observer maneuver at 125 s, true target bearings are not tracked by the estimated target bearing (magenta circles). However, as time goes on, true target bearings are tracked by estimated target bearings.
Sensors 22 05512 g003
Figure 4. The OSPA and the cardinality of the proposed GM-PHD in the case where B = 50 is utilized. As time goes on, OSPA gradually approaches to zero and E s t i m a t e converges to T r u e in cardinality.
Figure 4. The OSPA and the cardinality of the proposed GM-PHD in the case where B = 50 is utilized. As time goes on, OSPA gradually approaches to zero and E s t i m a t e converges to T r u e in cardinality.
Sensors 22 05512 g004
Figure 5. The OSPA and the cardinality of the proposed GM-PHD in the case where B = 1 is utilized. Setting B = 1 corresponds to the case where we apply a single EKF in the GM-PHD. Compared with the case where B = 50 is used (Figure 4), the estimation accuracy decreases considerably.
Figure 5. The OSPA and the cardinality of the proposed GM-PHD in the case where B = 1 is utilized. Setting B = 1 corresponds to the case where we apply a single EKF in the GM-PHD. Compared with the case where B = 50 is used (Figure 4), the estimation accuracy decreases considerably.
Sensors 22 05512 g005
Figure 6. The OSPA and the cardinality of the proposed GM-PHD in the case where B = 10 is utilized.
Figure 6. The OSPA and the cardinality of the proposed GM-PHD in the case where B = 10 is utilized.
Sensors 22 05512 g006
Figure 7. The OSPA and the cardinality of the proposed GM-PHD in the case where B = 20 is utilized.
Figure 7. The OSPA and the cardinality of the proposed GM-PHD in the case where B = 20 is utilized.
Sensors 22 05512 g007
Table 1. Computational load analysis.
Table 1. Computational load analysis.
B OneMCtime
16 s
106 s
207 s
5014 s
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, J. Tracking Multiple Targets Using Bearing-Only Measurements in Underwater Noisy Environments. Sensors 2022, 22, 5512. https://doi.org/10.3390/s22155512

AMA Style

Kim J. Tracking Multiple Targets Using Bearing-Only Measurements in Underwater Noisy Environments. Sensors. 2022; 22(15):5512. https://doi.org/10.3390/s22155512

Chicago/Turabian Style

Kim, Jonghoek. 2022. "Tracking Multiple Targets Using Bearing-Only Measurements in Underwater Noisy Environments" Sensors 22, no. 15: 5512. https://doi.org/10.3390/s22155512

APA Style

Kim, J. (2022). Tracking Multiple Targets Using Bearing-Only Measurements in Underwater Noisy Environments. Sensors, 22(15), 5512. https://doi.org/10.3390/s22155512

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