Next Article in Journal
Improved Landsat Operational Land Imager (OLI) Cloud and Shadow Detection with the Learning Attention Network Algorithm (LANA)
Previous Article in Journal
High Spatiotemporal Estimation of Reservoir Evaporation Water Loss by Integrating Remote-Sensing Data and the Generalized Complementary Relationship
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Three-Dimensional Signal Source Localization with Angle-Only Measurements in Passive Sensor Networks

1
National Key Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
2
China Academy of Launch Vehicle Technology, Beijing 100076, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(8), 1319; https://doi.org/10.3390/rs16081319
Submission received: 27 February 2024 / Revised: 6 April 2024 / Accepted: 7 April 2024 / Published: 9 April 2024
(This article belongs to the Topic Multi-Sensor Integrated Navigation Systems)

Abstract

:
Some passive sensors can provide only relative angles of a signal source. To obtain the signal source location, multiple passive sensors can be constructed into a passive sensor network through communication links. This paper investigates the source localization problem with angle-only measurements in three-dimensional space. First, we present an intersection localization method, which estimates the target position by minimizing the sum of distances between lines formed by angle-only measurements. It has the same target position estimate as the widely used least-squares (LS) method, but with a lower computational cost. Furthermore, considering the differences in measurement accuracy of sensors, the weighted least-squares (WLS) algorithm can achieve better localization performance than the LS method. Unfortunately, since the coefficient matrix and the noise vector are correlated, the WLS method is biased. The bias-compensation WLS (BCWLS) method is also presented in this paper to reduce the bias by estimating the correlation between the coefficient matrix and the pseudolinear noise vector. To evaluate the performance of the presented algorithms, numerical simulations are conducted, indicating that the superiority of the intersection localization method in computational cost and the superiority of the BCWLS method in localization accuracy.

1. Introduction

Signal source localization has applications in many fields, such as surveillance, guidance, and tracking [1,2,3]. In terms of whether to actively transmit electromagnetic signals, sensors for target localization can be divided into active sensors and passive sensors. Active sensors can achieve single-station positioning by transmitting signals to measure the distance to the target [4]. Passive sensors do not emit electromagnetic signals; therefore, they have a strong concealment [5]. Generally speaking, the types of passive sensors’ measurements include time delay (TD), Doppler shift (DS), angle of arrival (AOA), angle rate, or their combinations [6,7]. Among them, signal source localization with angle-only measurements has attracted considerable attention for many years, and it has been applied in many fields, including radar, sonar, navigation, and communications [8,9,10,11]. The signal source position cannot be obtained by the angle-only measurement in one snapshot with one passive sensor. Therefore, in order to estimate the signal source position, a widely used approach is to connect passive sensors distributed at different positions through communication links, and then design appropriate positioning algorithms to estimate the source position.
In the three-dimensional (3D) space, each passive sensor can measure the relative azimuth angle and elevation angle of the signal source, and output them as angle-only measurements. Therefore, it is obvious that each angle-only measurement can form a line with the sensor position as the endpoint. In the absence of noise, the signal source position are on a straight line, and the lines formed by the angle-only measurements of different sensors from the same signal source intersect at the signal source position. However, in practice, due to the presence of measurement error, including the angle measurement error and the sensor self-positioning error, even the lines formed by angle measurements from the same target are unlikely to intersect at the same point, but they intersect in a small volume. The small volume is determined by the level of the measurement noise. In [11], following this concept, a test statistic is constructed using the distance between the lines formed by angle-only measurements to complete the data association of angle-only measurements from different sensors.
Although the concept of target localization may be simple, it is challenging to obtain a robust and efficient estimate of the target position. Due to the nonlinearity of the signal source position and the angle-only measurement, signal source localization based on angle-only measurements is a nonlinear estimation problem [12,13]. To solve this problem, a pseudolinear equation regarding the signal source position is constructed by linearizing the angle measurement equation [6,8,10,14]. The least-squares (LS) method is an efficient and intuitive method to solve the pseudolinear overdetermined equation [15]. It is widely used because of its low complexity and stability. However, in the constructed pseudolinear overdetermined equation, the coefficient matrix is also contaminated by the measurement noise, which is not considered in the LS method [16,17]. The total least-squares (TLS) method not only considers the error of the data vector but also takes into account the perturbation of the coefficient matrix [18,19]. It is a natural generalization of the LS method [14,20]. In addition, the signal source position can also be estimated based on the spatial relationship between the source position and the lines formed by angle-only measurements [11,21,22,23].
The LS and TLS methods do not take into account the difference in the measurement noise levels of different sensors and they assume that the weights for different measurements are the same. Furthermore, to obtain a more accurate target position estimate, the weighted least-squares algorithm has been presented by considering the difference in measurement noise levels between different sensors [8,24]. The WLS algorithm is widely used due to its stability and accuracy [8,15,25,26]. Although the WLS method provides a more accurate target position estimate compared to the LS method, the WLS method is biased due to the presence of measurement noise not only in the pseudolinear noise vector but also in the coefficient matrix, especially if the angle measurement noise level is high [27]. To reduce the bias of the target position estimate, a series of methods have been investigated [28,29]. In 2D space, the bias-compensation operation is studied, which can reduce the bias to a certain extent [28]. In addition, from the perspective of reducing the noise of the coefficient matrix, the instrumental variables constructed in [29] can partly overcome the bias. Therefore, if the target is stationary or the observations of the sensors are synchronized, the above series of localization algorithms can be used to locate the target. However, if the target is moving and the sensor observations are not synchronized, it may be inappropriate to ignore the target movement.
To estimate the target motion state under asynchronous sensor observations, the gross LS method and the linear LS method are presented in [30]. In addition, a series of filtering algorithms have also been applied to estimate the target motion state with angle-only measurements [31]. Because of the nonlinearity of angle measurements and the target position, it is necessary to use a nonlinear tracking algorithm to estimate the target state. First, the standard Kalman filter algorithm has been used to process the pseudolinear equations constructed by the angle-only measurements, called pseudolinear Kalman filter (PLKF) [32,33]. In addition, many nonlinear filtering algorithms have also been used for target tracking with angle-only measurements. The extended Kalman filter (EKF) approximately obtains a linear observation equation through a first-order Taylor expansion. However, due to the neglect of high-order terms, it faces a filtering divergence problem when the measurement noise level is high [34]. Sigma-point Kalman filtering algorithms assume that approximating a distribution is easier than approximating a nonlinear function itself. They approximate the distribution of the target state by selecting a set of sigma points and use a Kalman filter to update the target state, such as the cubature Kalman filter (CKF) [35] or the unscented Kalman filter (UKF) [36,37]. In addition, the particle filter is also used for target tracking with angle-only measurement [38].
In this paper, we investigate the source localization problem with angle-only measurements from the passive sensor network connected by communication links in one snapshot. Firstly, we formulate an algorithm named intersection localization method, which is obtained by minimizing the distances between lines formed by the angle-only measurements. The LS method is also studied in this paper, which estimates the target position by solving the equation constructed using the angle measurements. By comparing the results of the presented intersection localization method and the LS method, we prove that the target position estimates of the intersection localization method and the LS method are the same. However, compared with the LS method, the coefficient matrix size of the proposed intersection localization method is smaller. Therefore, it requires smaller multiplication and addition operations. In theory, the computational cost of the intersection localization method is lower. Furthermore, the TLS method is studied, which takes into account not only errors in the data vector but also errors in the coefficient matrix. In practice, different sensors may have different levels of measurement noise. The WLS method is determined by considering the difference in measurement noise of different sensors. It is found that the WLS method is biased due to the presence of noise not only in the pseudolinear noise vector but also in the coefficient matrix. Therefore, we further formulate a bias-compensation WLS (BCWLS) algorithm by compensating the positioning bias of the WLS algorithm. To analyze the localization performance, the Cramér–Rao lower bound (CRLB) of the target position estimators based on angle-only measurements is also derived in this paper.
In numerical simulation, we analyze the impacts of angle measurement noise and sensor self-positioning noise on the target positioning accuracy of several algorithms. Numerical results show that the performance of the BCWLS method is closer to the CRLB performance among the above algorithms, proving the performance improvement of the BCWLS algorithm. In addition, it is verified that the intersection localization method and the LS method have the same target position estimate. Furthermore, we compare the running time of several algorithms under the same conditions, and the results verify the inference that the intersection localization algorithm requires a lower computational cost than the LS algorithm. Simulation results validate the performance of the proposed methods and demonstrate the improvement compared with the previous algorithms. In summary, the main contributions of this article can be highlighted as follows:
(1)
We present an intersection localization method with angle-only measurements in a passive sensor network. The presented intersection localization method has the same target position solution as the widely used and computationally efficient LS method but has a lower computational cost.
(2)
We present a bias computation WLS estimator for target localization using angle-only measurements. The presented method can compensate for the positioning bias of the WLS method, thereby improving the accuracy of the target position estimate.
(3)
We derive the CRLB for estimating the target position under the condition of sensor self-positioning error. Numerical simulations evaluate the superiority of the intersection localization method in terms of computational cost and the superiority of the BCWLS method in terms of localization accuracy.
The rest of this article is sectioned as follows. In Section 2, we construct a target localization scenario with angle-only measurements in a passive sensor network and provide a series of target localization methods. In Section 2.2, we present an intersection localization method. In Section 2.3, Section 2.4 and Section 2.5, the LS, TLS, and WLS methods are introduced, respectively. In Section 2.6, we analyze the solution of the WLS method and further present a bias-compensation WLS method. Section 2.7 gives the CRLB of the target position estimators with angle-only measurements in passive sensor network. Section 3 examines the performance of the above methods via numerical simulations. Section 4 summarizes the results of the simulation results and provides the advantages and disadvantages of the two presented algorithms. Finally, Section 5 concludes this article.

2. Materials and Methods

2.1. Signal Model of Angle-Only Measurements

We consider a 3D target localization scenario, which includes M targets and a passive sensor network. This passive sensor network consists of N passive sensors that can only provide the angle information of the target, such as passive radar, photoelectric sensor, and infrared sensor. In 3D space, angle-only measurements consist of azimuth and elevation angles. The true positions of the targets need to be estimated by the angle-only measurements. The position of each passive sensor is measured by its mounted positioning device, such as inertial navigation systems and the Global Positioning System (GPS). At time t, the true position of the nth sensor is denoted as p n o ( t ) = [ x n , s o ( t ) , y n , s o ( t ) , z n , s o ( t ) ] T , n = 1 , 2 , , N . The elements x n , s o ( t ) , y n , s o ( t ) and z n , s o ( t ) , respectively, represent the coordinates of the nth sensor on the x-axis, y-axis, and z-axis in the common coordinate system at time t. The notation [ · ] T denotes the transpose operation. The true position of the mth target at time t is denoted as g m o ( t ) = [ x m , g o ( t ) , y m , g o ( t ) , z m , g o ( t ) ] T , m = 1 , 2 , , M . Figure 1 shows the 3D target localization scenario.
Assume that for the nth sensor, there are a total of N n observations. Each observation can receive radiation signals from multiple signal sources. The time instant of each observation can be defined as t k , n , k = 1 , , N n . For each radiation signal, the passive sensor can measure its AOA, termed angle-only measurement. At time instant t k , n , the measurement of the nth sensor position can be expressed as
p k , n ( t k , n ) = p n o ( t k , n ) + Δ p n ( t k , n ) = [ x n , s ( t k , n ) , y n , s ( t k , n ) , z n , s ( t k , n ) ] T
where the vector Δ p n ( t k , n ) is the sensor self-positioning noise. Without loss of generality, it is assumed to follow a zero-mean Gaussian distribution with the covariance matrix R s , n ( t k , n ) = E ( Δ p n ( t k , n ) Δ p n T ( t k , n ) ) , namely Δ p n ( t k , n ) N ( 0 , R s , n ( t k , n ) ) . The notation E ( · ) is the expectation operation. At time instant t k , n , the real position of the mth target can be denoted as
g m o ( t k , n ) = [ x m , g o ( t k , n ) , y m , g o ( t k , n ) , z m , g o ( t k , n ) ] T .
Assuming that the nth sensor receives a total of L k , n signals corresponding to L k , n angle-only measurements during the kth observation. Therefore, the angle-only measurements of the kth observation of the nth sensor can be indexed by a triple ( l , k , n ) , l = 1 , , L k , n . It should be noted that due to the presence of false alarms and missed detections, there may be a situation where L k , n is not equal to the number of targets M. In this way, we can establish a one-to-one index for all angle-only measurements of all sensors.
The set of angle-only measurements of the nth sensor at instant t k , n can be represented as L k , n , and
L k , n = { j | j = ( l , k , n ) } .
The cardinality of L k , n satisfies L k , n = | L k , n | , where the notation | · | represents the cardinality of the set. Denote the set L n and L , which represent the set of measurements by the nth sensor and the set of measurements by N sensors, respectively. The sets L n and L are
L n = k = 1 N n L k , n , L = n = 1 N L n
where the notation ∪ represents the union operation. Therefore, we define N s as the total number of the angle-only measurements of N sensors,
N s = | L | = n = 1 N k = 1 N n L k , n .
It is necessary to associate the measurements of N sensors before target localization. Each angle-only measurement can only come from one of M targets or a false alarm. We define a set M = { 0 , 1 , , M } as the associated results’ index set, where the index 0 represents the false alarm, and the index 1 , , M represent target 1 , , M , respectively. Each angle-only measurement corresponds only to one member of the set M . In our previous work [11], we constructed a test statistic using the minimum distance between the lines formed by the angle-only measurements of different sensors and then achieved the data association of the angle-only measurements. Then, we can obtain a mapping ψ : L M , which represents that any element in L corresponds to a unique element in M . Therefore, after data association, the index set L can be divided into M + 1 disjoint set, denoted as B 0 , B 1 , , B M , where B 0 represents the index set of angle-only measurements from the false alarm, and B m represents the index set of angle-only measurements from the mth target, where B m is
B m = { j | ψ ( j ) = m , j L }
where the sets B i and B j satisfy B i B j = , i , j M , i j . The notation ∩ represents the intersection operation of sets. In this paper, we study the target localization problem with angle-only measurements from multiple passive sensors. Therefore, for simplicity, we assume that the data association of angle-only measurements has been completed, and the mapping ψ is known. For the mth target, its angle-only measurements index set is B m . We assume that | B m | = L m . For the angle-only measurement index n B m , assuming that the sensor position measurement of its corresponding sensor can be expressed as
p n = p n o + Δ p n
where the vector p n represents the sensor position measurement, the vector p n o represents the true sensor position, and the vector Δ p n represents the sensor self-positioning noise assumed to be a zero-mean Gaussian distribution with covariance matrix R p , n .
The positions of the sensors with respect to the L m angle-only measurements can be expressed as
p = p o + Δ p
where p = [ p 1 T , p 2 T , , p L m T ] T , p o = [ p 1 o T , p 2 o T , , p L m o T ] T , and Δ p = [ Δ p 1 o T , Δ p 2 o T , , Δ p L m o T ] T . The covariance matrix of Δ p is R p = blkdiag ( R p , 1 , R p , 2 , , R p , L m ) , where the notation blkdiag ( A 1 , A 2 , , A N ) represents the operation of constructing a block diagonal matrix with the matrices A 1 , A 2 , , A N .
The true value of the azimuth angle θ n o and elevation angle φ n o corresponding to the nth angle-only measurement can be expressed as [34]
θ n o = f θ ( g m o , p n o ) = tan 1 y m , g o y n , s o x m , g o x n , s o
φ n o = f φ ( g m o , p n o ) = tan 1 z m , g o z n , s o ( x m , g o x n , s o ) 2 + ( y m , g o y n , s o ) 2
where g m o = [ x m , g o , y m , g o , z m , g o ] T represents the target position, θ n o [ π , π ] , φ n o [ π 2 , π 2 ] , and tan 1 is the 4-quadrant arctangent function.
In practice, the angle-only measurements θ n and φ n can be expressed as
θ n = θ n o + Δ θ n = f θ ( g m o , p n ) + Δ θ n
φ n = φ n o + Δ θ n = f φ ( g m o , p n ) + Δ φ n
where the elements Δ θ n and Δ φ n represent the angle measurement noise, θ n [ π , π ] , and φ n [ π 2 , π 2 ] . Without loss of generality, we assume that Δ θ n and Δ φ n are zero-mean Gaussian noises with variances σ Δ θ n 2 and σ Δ φ n 2 , namely, Δ θ n N ( 0 , σ Δ θ n 2 ) and Δ φ n N ( 0 , σ Δ φ n 2 ) , respectively.
From (11) and (12), due to the presence of sensor positioning error, (11) and (12) can be expressed as
θ n = θ n o + e θ , n + Δ θ n = θ n o + δ θ , n
φ n = φ n o + e φ , n + Δ φ n = φ n o + δ φ , n
where the elements e θ , n and e φ , n denote the measurement noise of the azimuth and elevation angles, respectively, caused by the sensor self-positioning error. They can be expressed as
e θ , n = f θ ( g m o , p n ) p n T Δ p n , e φ , n = f φ ( g m o , p n ) p n T Δ p n .
Therefore, the elements δ θ , n and δ φ , n represent the total measurement error of the azimuth angle and elevation angle, respectively. Assuming that the sensor self-positioning noise and the angle measurement noise are independent of each other, the variances of δ θ , n and δ φ , n are
σ θ , n 2 = σ Δ θ n 2 + b θ , n T R p , n b θ , n
σ φ , n 2 = σ Δ φ n 2 + b φ , n T R p , n b φ , n
where
b θ , n = f θ ( g m o , p n ) p n , b φ , n = f φ ( g m o , p n ) p n
The nth angle-only measurement vector is θ n = [ θ n , φ n ] T , and its corresponding measurement error vector is δ θ , n = [ δ θ , n , δ φ , n ] T . The covariance matrix of δ θ , n is R θ , n = diag ( σ θ , n 2 , σ φ , n 2 ) , where the notation diag ( a 1 , a 2 , , a N ) represents the operation of constructing a diagonal matrix using a 1 , a 2 , , a N .
The angle-only measurement vector of the mth target is
θ = θ o + δ θ
where θ = [ θ 1 T , θ 2 T , , θ L m T ] T , θ o = [ θ 1 o T , θ 2 o T , , θ L m o T ] T , and δ θ = [ δ θ , 1 T , δ θ , 2 T , , δ θ , L m T ] T . The covariance matrix of δ θ is R θ = blkdiag ( R θ , 1 , R θ , 2 , , R θ , L m ) .

2.2. Intersection Localization Method

It is obvious that each angle-only measurement of a passive sensor forms a line in 3D space with the sensor position as the endpoint. Under the condition that this measurement is not a false alarm, and there is no measurement noise, the target is on that line. With multiple angle-only measurements from different sensors, the target can be located. The line formed by the nth measurement can be expressed as
L n : x n = p n + α n e n , α n R
where the element α n is a distance parameter indicating the distance to p n , the vector e n = [ e n , x , e n , y , e n , z ] T denotes the normalized direction vector, and
e n = [ cos ( θ n ) cos ( φ n ) , sin ( θ n ) cos ( φ n ) , sin ( φ n ) ] T
where e n = 1 , and the notation · over a vector denotes the 2 -norm.
In [11], the minimum distance between two lines is used to solve the distance parameters of the two lines. In this part, we extend the algorithm such that the distance parameters α = [ α 1 , , α L m ] T are solved as a whole.
For simplicity, assume that the elements in B m can be expressed by B m = { 1 , , L m } . The square of the Euclidean distance between two points on two lines L i , L j is
d i , j = p j p i + ( e i , e j ) [ α i , α j ] T 2 = P I i , j 1 + E I i , j α 2 = α T I i , j E T E I i , j α + 2 α T I i , j E T P I i , j 1 + 1 T I i , j P T P I i , j 1
where P = [ p 1 , , p L m ] R 3 × L m , E = [ e 1 , , e L m ] R 3 × L m , the vector 1 R L m × 1 denotes an all-one vector of length L m , and the elements of matrix I i , j R L m × L m can be expressed as
[ I i , j ] k , m = 1 k = m = j 1 k = m = i 0 else .
It should be noted that the distance d i , j is actually the square of the Euclidean distance, instead of the Euclidean distance itself. To extend the minimum total distance between all the lines, instead of two lines, we can formulate the optimization problem as
min α d = 1 i < j L m d i , j
where only one ( i , j ) pair is taken into account in the summation operation since d i , j = d j , i .
It can be proved that the total distance can be expressed by
d = 1 i < j L m d i , j = α T K α + 2 α T K ¯ 1 + 1 T K ˜ 1
where
K = 1 i < j L m I i , j E T E I i , j = K T
K ¯ = 1 i < j L m I i , j E T P I i , j
K ˜ = 1 i < j L m I i , j P T P I i , j = K ˜ T .
It is evident that K = K T , and K ˜ = K ˜ T . For any matrix G R L m × L m , we have
[ I i , j G I i , j ] k , m = G i , i k = m = i G i , j k = i , m = j , or ,   k = j , m = i G j , j k = m = j 0 else .
With this equation, we can verify that
1 i < j L m I i , j G I i , j = L Diag ( G ) G
where the notation Diag ( · ) corresponds to a square matrix representing a diagonal matrix with the same diagonal elements as the input matrix.
Note that all the diagonal elements of E E T are 1 since e T e = 1 ; then,
K = 1 i < j L m I i , j E T E I i , j = L m I E T E
K ¯ = L m Diag ( E T P ) E T P
K ˜ = L m Diag ( P T P ) P T P .
To solve for the parameter α that minimizes d, we take the derivative of (25) with respect to α ,
d d d α = 2 K α + 2 K ¯ 1 .
Setting (34) to zero and solving the equation, we obtain a solution
d d d α = 0 α opt = K 1 K ¯ 1
Putting α opt into (25), we obtain the minimal distance d min as
d min = 1 T ( K ˜ K ˜ K 1 K ˜ ) 1 .
In particular, if
K ˜ = K ˜ K 1 K ˜
then
d min = 0 .
In practice, the size of K is L m × L m , and thus the computational cost to calculate K 1 is huge. According to the Sherman–Woodbury–Morrison equation,
( A + B C ) 1 = A 1 A 1 B ( I + C A 1 B ) 1 C A 1
From K in (31), we let A L m I , B E T / L m , and C E ; then, the matrix inversion operation K 1 can be simplified as
K 1 = 1 L m I + 1 L m 2 E T ( I E E T / L m ) 1 E
where the only matrix inverse operation is over a 3 × 3 matrix I E E T / L m .
Consequently,
K 1 K ¯ = Diag ( E T P ) 1 L m E T P + 1 L m 2 E T ( I E E T / L m ) 1 E ( L m Diag ( E T P ) E T P ) = Diag ( E T P ) 1 L m E T P + 1 L m E T ( I E E T / L m ) 1 E Diag ( E T P ) + 1 L m E T ( I E E T / L m ) 1 ( I E E T / L m I ) P = Diag ( E T P ) + 1 L m E T ( I E E T / L m ) 1 ( E Diag ( E T P ) P ) .
The optimal weight can be written as
α opt = K 1 K ¯ 1 = diagvec ( E T P ) 1 L m E T ( I E E T / L m ) 1 ( E diagvec ( E T P ) P 1 )
where the notation diagvec ( · ) with a square matrix input represents a vector formed by the diagonal elements of the input matrix.
For the ith line, the point corresponding to α i is
x ( t i ) = p i + α i e i
which, however, may be different for different lines. Therefore, we take the mean as the final estimate of the target position g m o , namely,
g ^ m = 1 L m i = 1 L m x ( t i ) = 1 L m ( P 1 + E α opt ) = 1 L m ( P 1 E K 1 K ¯ 1 ) = ( L m I E E T ) 1 ( P 1 E diagvec ( E T P ) ) = ( L m I E E T ) 1 i = 1 L m ( I e i e i T ) p i .

2.3. Least-Squares Method

The LS method is a very classic and widely used method to solve linear overdetermined equations [15]. In this section, we use the LS method to solve the pseudolinear equations constructed by the angle measurement equations to estimate the target position. Firstly, let us derive the pseudolinear equations about a target position.
In (9), taking the tangent on both sides of the equation and replacing tan θ n o with sin θ n o / cos θ n o , it can be rewritten as [8]
u θ o T p n o = u θ o T g m o
where u θ o = [ sin θ n o , cos θ n o , 0 ] T .
Replace p n o and θ n o with their corresponding measurement values p n and θ n , respectively. If the angle measurement noise is small, we have approximations of sin δ θ , n δ θ , n and cos δ θ , n 1 , then
sin θ n = sin ( θ n o + δ θ , n ) sin θ n o + δ θ , n cos θ n o
cos θ n = cos ( θ n o + δ θ , n ) cos θ n o δ θ , n sin θ n o .
Taking (46) and (47) into (45), the pseudolinear equation obtained from the azimuth angle can be expressed as
z θ , n = A θ , n g m o + η θ , n
where z θ , n = u θ , n T p n , η θ , n d n o cos φ n o δ θ , n + u θ , n o T Δ p n , A θ , n = u θ , n T , u θ , n = [ sin θ n , cos θ n , 0 ] T , d n o = p n o g m o . In (48), we have used the geometric relationship ( x g , m o x n , s o ) cos θ n o + ( y g , m o y n , s o ) sin θ n o = d n o cos φ n o .
Similarly, Equation (10) can be rewritten as [8]
u φ , n o T p n o = u φ , n o T g m o
where u φ , n o = [ cos θ n o sin φ n o , sin θ n o sin φ n o , cos φ n o ] T . Similarly, from (46) and (47), we also have the following approximation about the elevation angle
sin φ n = sin ( φ n o + δ φ , n ) sin φ n o + δ φ , n cos φ n o
cos φ n = cos ( φ n o + δ φ , n ) cos φ n o δ φ , n sin φ n o .
Putting (46), (47), (50), (51), and the noisy sensor position p n into (49), the pseudolinear equation obtained by the elevation angle can be expressed as
z φ , n = A φ , n g m o + η φ , n
where z φ , n = u φ , n T p n , A φ , n = u φ , n T , and u φ , n = [ cos θ n sin φ n , sin θ n sin φ n , cos φ n ] T , η φ , n d n o δ φ , n + u φ , n o T Δ p n . Combining (48) and (52) in matrix form, we obtain a pseudolinear equation,
z n = A n g m o + η n
where z n = [ z θ , n , z φ , n ] T , and A n = [ u θ , n , u φ , n ] T R 2 × 3 . The vector η n denotes the pseudolinear noise vector,
η n = [ η θ , n , η φ , n ] T D n o δ θ , n + A n o Δ p n
where D n o = diag ( d n o cos φ n o , d n o ) and A n o = [ u θ , n o , u φ , n o ] T . The covariance matrix of η n is
R η , n = E [ η n , η n T ] = D n o R θ , n D n o T + A n o R p , n A n o T .
Combining (55) formed by L m angle-only measurements, we have
z = A g m o + η
where the coefficient matrix A and the data vector z can be expressed as
z = [ z 1 T , z 2 T , , z L m T ] T
A = [ A 1 T , A 2 T , , A L m T ] T .
The pseudolinear noise vector η is
η = [ η 1 T , η 2 T , , η L m T ] T D o δ θ + B o Δ p
where D o = blkdiag ( D 1 o , D 2 o , , D L m o ) and B o = blkdiag ( A 1 o , A 2 o , , A L m o ) . The covariance matrix of vector η can be expressed as
R η = blkdiag ( R η , 1 , R η , 2 , , R η , L m ) .
With L m observations, there are a total of 2 L m equations about the target position g m o , and there are three unknown parameters in g m o . Therefore, if L m 2 , g m o has a solution, its LS solution can be obtained by minimizing η , which can be expressed as [15]
g ^ m , LS = ( A T A ) 1 A T z .
By rewriting the target position g m o as ( A T A ) 1 A T A g m o , the bias of the LS estimate g ^ m , LS can be expressed as
E [ Δ g m , L S ] = E [ g m o g ^ m , LS ] = E [ ( A T A ) 1 A T η ] .
We can see that if the measurement noise is negligible and the pseudolinear noise vector satisfies E [ η ] = 0 , the LS estimator can be seen as approximately unbiased.
It can be proved that the following equations exist
I e n e n T = A n T A n
A n T z = A n T A n p n = ( I e n e n T ) p n .
Therefore, in (44) and (61), we have the following equations
( L m I E E T ) 1 = ( A T A )
i = 1 L m ( I e i e i T ) p i = A T z .
Therefore, the target position estimate g ^ m of the intersection localization method and the LS estimate g ^ m , L S are equal. Due to the fact that the size of A in the LS method is 2 L m × 3 , and the size of E in the intersection localization method is 3 × L m , their computational costs are different. In theory, the intersection localization method requires fewer multiplication and addition operations than the LS method.

2.4. Total Least-Squares Method

It should be noted that the construction of matrix A and vector z is based on the angle measurements and the sensor position measurements. Therefore, measurement noise not only exists in vector z but also in matrix A , and Equation (56) can actually be expressed as [18]
( A o + Δ A ) g m o = z o + Δ z
where A = A o + Δ A , and z = z o + Δ z . The matrix A o and the vector z o represent A and z with the noise values replaced by the true values, respectively.
By transferring items, (67) can be expressed as
( A o + Δ A | z o + Δ z ) g m o 1 = 0 .
The TLS solution of the target position can be obtained by minimizing the Frobenius norm of the matrix [ Δ A , Δ z ] . Therefore, the TLS estimate g ^ m , TLS can be obtained by solving the following constrained optimization problem, as in [18],
argmin g ^ m , TLS , Δ A TLS , Δ z TLS Δ A Δ z F s . t . ( A Δ A ) g m o = z Δ z
where the notation · F with a matrix entry represents the Frobenius norm of the matrix. The TLS estimate g ^ m , TLS of the target position g m o can be obtained by singular value decomposition (SVD) of the augmented matrix [ A , z ] [19].
Denote C = [ A , z ] R 2 L m × 4 , and the SVD of the matrix C is
C = U Σ V T
where the matrix Σ = diag ( σ 1 , , σ 4 ) . σ 1 , , σ 4 are the singular values of C and satisfy σ 1 σ 2 σ 4 . The matrix V can be partitioned as
V = V 11 v 12 v 21 v 22
where V 11 R 3 × 3 , v 12 R 3 × 1 , v 21 R 1 × 3 , and v 22 R 1 × 1 .
It has been proved in [19] that the TLS solution of (69) exists if and only if the element v 22 is non-singular, i.e., v 22 0 . In this case, the TLS estimate can be expressed as
g ^ m , TLS = v 12 v 22 1 .

2.5. Weighted Least-Squares Method

The LS method applies the same weight to different measurements. However, in practice, the accuracy of different measurements may be different and may be known a priori. We assume that the distributions of angle-only measurement noise and sensor self-positioning noise are known a priori. Therefore, the target localization accuracy can be improved by giving proper weights to different measurements, which is the WLS method [8,14,15].
In (56), by considering the covariance of η , the cost function of the WLS formulation is
J = ( z A g m o ) T W ( z A g m o )
where W is the weighting matrix,
W = R η 1 = blkdiag ( W 1 , W 2 , , W L m )
where W n = R η , n 1 , n = 1 , 2 , , L m .
The partial derivative of the cost function J for the target position g m o is
J g m o = 2 A T W z + 2 A T W A g m o .
Let J / g m o = 0 ; then, the WLS solution of the target position is [8]
g ^ m , WLS = ( A T W A ) 1 A T W z .
Rewriting the target position as ( A T W A ) 1 A T W A g m o , we obtain the estimate error of g ^ m , WLS as
Δ g m , WLS = g ^ m , WLS g m o = ( A T W A ) 1 A T W η .
The covariance matrix of g ^ m , WLS is [29]
R g , WLS = E [ Δ g m , WLS Δ g m , WLS T ] ( A o T W A o ) 1 A o T W E [ η η T ] W T A o ( A o T W A o ) 1 = ( A o T W A o ) 1
where we have used E [ η η T ] = W 1 .
It should be noted that the weighting matrix W requires the unknown true position of the target via D o and B o . To solve this problem, we obtain an initial target position estimate by using the LS method. The weighting matrix W can be obtained by this initial solution; then, the WLS estimate can be derived by (76).

2.6. Bias-Compensation WLS Method

This section firstly analyzes the bias issue of the WLS method introduced in Section 2.5. Then, a bias-compensated WLS method is studied by estimating the bias of the WLS estimate g ^ m , WLS .
The expectation of Δ g m , WLS in (77) can be expressed as
δ g , WLS = E [ Δ g m , WLS ] = E [ ( A T W A ) 1 A T W η ] ( A T W A ) 1 E [ A T W η ]
Since both the matrix A and the pseudolinear noise vector η contain measurement noise, the matrix A and the noise vector η are correlated. Therefore, E [ A T W η ] 0 , so that the WLS estimator is biased.
Since the matrix W is a block diagonal matrix, the expectation E [ A T W η ] can be rewritten as
E [ A T W η ] = n = 1 L m E [ A n T W n η n ] n = 1 L m m n .
Because the weighting matrix W n is the inverse of the covariance matrix R n , it is a symmetric matrix that can be expressed as
W n = a n b n b n c n
Putting (81) into (80), the vector m n can be written as
m n = g 1 , n + g 2 , n
where
g 1 , n = a n a 1 , n b θ , n T + b n ( a 2 , n b θ , n T + a 3 , n b φ , n T ) R p , n u θ , n o d n o cos φ n o σ Δ θ n 2 ( a n a 1 , n + b n a 2 , n )
g 2 , n = b n a 1 , n b θ , n T + c n ( a 2 , n b θ , n T + a 3 , n b φ , n T ) R p , n u φ , n o + c n d n o a 3 , n σ Δ φ n 2
where
a 1 , n = [ cos θ n o , sin θ n o , 0 ] T
a 2 , n = [ sin θ n o sin φ n o , cos θ n o sin φ n o , 0 ] T
a 3 , n = [ cos θ n o cos φ n o , sin θ n o cos φ n o , sin φ n o ] T
Therefore, according to (77), the bias-compensation WLS estimate of the target position is
g ^ m , W L S b c = g ^ m , W L S δ g , W L S .
It should be noted that the calculation of δ g , W L S depends on the true azimuth angle θ n o , elevation angle φ n o , and distance d n o . To ensure accuracy, we can use the WLS solution g ^ m , W L S of the target position to estimate δ g , W L S and then derive the final solution by (88).

2.7. Cramér–Rao Lower Bound

It is well-known that the CRLB establishes a lower bound on the performance of estimators [8,14,24]. Therefore, it is widely used to calculate the best estimation accuracy in theory and evaluate the performance of estimators. For the scenario in this article, with the presence of the sensor self-positioning noise, the unknown parameters include the target position g m and the positions of the sensors p o . Therefore, under the assumption that the angle-only measurements’ noise and the sensor self-positioning noise are independent and all subject to the zero-mean Gaussian distribution, we would like to derive the CRLB of the parameter vector Θ = [ g m o T , p o T ] T R ( 3 L m + 3 ) × 1 . The measurement vector consists of L m angle-only measurements θ and corresponding position measurements p of sensors, denoted by m = [ θ T , p T ] T R 5 L m × 1 . Under the given parameters Θ , the logarithm of the joint probability density function of m is
l n p ( m | Θ ) = l n p ( θ | Θ ) + l n p ( p | Θ ) = κ 1 2 ( θ θ o ) T R θ 1 ( θ θ o ) 1 2 ( p p o ) T R p 1 ( p p o )
where κ is a constant.
According to (89), the Fisher information matrix (FIM) of Θ is
FIM ( Θ ) = E l n 2 p ( m | Θ ) Θ Θ T
The CRLB matrix of an estimate of Θ is the inverse of FIM ( Θ ) :
CRLB ( Θ ) = FIM 1 ( Θ ) = θ o g m o T R θ 1 θ o g m o θ o g m o T R θ 1 θ o p o θ o p o T R θ 1 θ o g m o θ o p o T R θ 1 θ o p o + R p 1 1 = X Z Z T Y 1
where
θ o g m o = blkdiag θ 1 o g m o , θ 2 o g m o , , θ L m o g m o = blkdiag ( B 1 , B 2 , , B L m ) R 3 L m × 2 L m
θ o p o = blkdiag θ 1 o p 1 o , θ 2 o p 2 o , , θ L m o p L m o = blkdiag ( B 1 , B 2 , , B L m ) R 3 L m × 2 L m
where
B n = θ n o g m o = θ n o g m o , φ n o g m o
θ n o g m o = y n , s o y m , g o ( x m , g o x n , s o ) 2 + ( y m , g o y n , s o ) 2 , x m , g o x n , s o ( x m , g o x n , s o ) 2 + ( y m , g o y n , s o ) 2 , 0 T
φ n o g m o = cos θ n o sin φ n o d n o , sin θ n o sin φ n o d n o , sin φ n o d n o T
According to the block matrix inversion formula, we can obtain the CRLB of target position g m as
CRLB ( g m o ) = ( X Z Y 1 Z T ) 1
According to the matrix inversion lemma, (97) can be rewritten as
CRLB ( g m o ) = X 1 + X 1 Z ( Y Z T X 1 Z ) 1 Z T X 1
where X 1 represents the CRLB of the target position estimate if there is no self-positioning error of the sensors, and X 1 Z ( Y Z T X 1 Z ) 1 Z T X 1 is the increase in the target positioning error caused by the sensor self-positioning error.

3. Results

In this section, we compare the target localization performance of the above estimators in different scenarios and different noise levels. We first considered a scenario that included four sensors and one target, and the sensors could only provide the angle information of the target. For simplicity and without loss of generality, false alarms and missed detections were not considered. As shown in Table 1, we established the positions of the sensors and the target.
We assumed that the sensor self-positioning noise levels of the four sensors were the same, and their covariance matrices was denoted as R s , n = diag ( σ s 2 , σ s 2 , σ s 2 ) , n { 1 , 2 , 3 , 4 } , where σ s represents the sensor self-positioning noise level. It was also assumed that the angle measurement noises of each sensor, including azimuth and elevation angles, followed a zero-mean Gaussian distribution and were independent of each other. The variances of angle measurements were denoted as σ Δ θ n 2 = σ Δ φ n 2 = σ θ 2 , n { 1 , 2 , 3 , 4 } , where σ θ represents the angle measurement noise level. All the following simulations were performed on a laptop with an Intel core i7-12700H and 32-GB RAM. The software version was MATLAB R2020a. In the simulations, except for the TLS method using the built-in singular value decomposition (SVD) function in MATLAB, no other specialized toolboxes were used.

3.1. Statistical Metrics

The bias norm (BNorm) and root-mean-square error (RMSE) were used to evaluate the target localization accuracy of the estimators. BNorm reflects the bias performance of the estimators, and a smaller BNorm indicates a smaller localization bias. BNorm and RMSE can be calculated by L Monte Carlo simulations as
BNorm = 1 L l = 1 L ( g ^ ( l ) g o )
RMSE = 1 L l = 1 L g ^ ( l ) g o 2 1 / 2
where g ^ ( l ) represents the estimated value of the real target position g o at the lth Monte Carlo simulation. At different error conditions, we set the number of Monte Carlo simulations to L = 20 , 000 . The square root of the trace of the CRLB matrix of (98) was referred to as root CRLB and was used as the theoretical boundary for the RMSE of the target position estimators.

3.2. The Impacts of Sensor Self-Positioning Noise and Angle Measurement Noise

To evaluate the influence of the measurement noise level on the target localization accuracy, we compared the RMSE and BNorm performance of the above algorithms under different error conditions. First, we evaluated the impact of the angle measurement noise level on the target localization accuracy of each algorithm. We fixed the sensor self-positioning noise level σ s at 10 m, and the angle measurement noise level varied from 0.1° to 6.1° at intervals of 0.5°. Figure 2 illustrates the impact of angle measurement noise level on the target localization performance of the above target position estimators.
Form Figure 2, it appears that the proposed intersection localization algorithm and the LS algorithm had the same BNorm and RMSE, verifying the inference that the target position solutions of the two algorithms were equal. The BNorm of the TLS and WLS algorithms was very close, slightly better than that of LS algorithm. In addition, the BCWLS algorithm had the smallest BNorm among several algorithms, which showed the effectiveness of the proposed algorithm. We can see that the BNorm of all the methods was less than 5 m if the angle measurement noise level was less than 1.1°. However, with the increase in σ θ , the BNorm of the other methods, other than the BCWLS algorithm, increased rapidly. When the σ θ was 4.1°, the BNorm of the WLS and TLS algorithms was close to 40 m, the BNorm of the LS algorithm and intersection localization algorithm was close to 60 m, while the BNorm of BCWLS algorithm was still less than 5 m. Form Figure 2b, it can be seen that the proposed BCWLS algorithm had the best RMSE among the above algorithms, and its RMSE was much closer to the root CRLB.
Next, we evaluated the impact of sensor self-positioning accuracy on target localization performance. We fixed the angle measurement noise level σ θ at 2°, and the sensor self-positioning noise level σ s varied from 5 m to 30 m, at intervals of 5 m. Figure 3 illustrates the impact of the sensor self-positioning noise level on the target localization performance of the above target position estimators. Figure 3a gives the BNorm results of the estimators. It shows that the BCWLS algorithm had the best BNorm performance among the considered algorithms. The BNorm results of the TLS and WLS methods were very close to each other and slightly better than the BNorm of the LS method and intersection localization method. In addition, we can see that the BNorm of the above various algorithms remained stable as the sensor self-positioning noise level σ s increased. The reason is that the bias of the target position estimators is caused by the nonlinearity of the angle-only measurement and the target position. The impacts of the sensor self-positioning noise level on the RMSE of several estimators are shown in Figure 3b. It is obvious that the BCWLS method had the best RMSE results among several algorithms. The RMSE of the TLS algorithm, LS algorithm, and intersection localization method was very close and larger than that of the WLS method. In addition, the LS algorithm and intersection localization algorithm had the same RMSE and BNorm results, which further verified the inference that both algorithms had the same target position estimates.
To illustrate the advantages of the BCWLS algorithm in terms of bias, we directly show the target position estimates of the WLS and BCWLS algorithms. Figure 4 shows the 100 target position estimates of the two methods under the conditions σ θ = 4 ° and σ s = 10 m. The shapes ◯ and □ represent the target position estimates obtained by the WLS and BCWLS methods, respectively. The shape Δ represents the true position of the target. The shape ✩ is the mean of the 100 target position estimates obtained by the WLS or BCWLS method. It should be noted that the measurements used for the 100 target position estimates by these two methods were the same. We can see that the mean of the target position estimates by the WLS method deviates from the true position of target. In contrast, the mean of the target position estimates by BCWLS method is closer to the true value of the target position. Thus, the BCWLS method generally has a better localization performance than the WLS method.

3.3. The Impact of the Number of Sensors

It is well known that target localization performance depends to a large extent on the number of sensors and the geometry of the sensors and the target. Therefore, we changed the number of sensors to three, removed the fourth sensor in Table 1, and kept the positions of the remaining sensors and the target unchanged.
Figure 5 illustrates the target localization performance of the considered algorithms with σ s fixed at 10 m and σ θ varying from 0.1° to 6.1° in the case of three sensors. Form Figure 2a and Figure 5a, it is obvious that after the removal of sensor # 4 , the BNorm of several algorithms changed under different angle measurement noise levels. Under the condition of σ θ = 4.1 ° , when the number of sensors was four, the BNorm of the LS algorithm and intersection localization algorithm was about 45 m, and the BNorm of the BCWLS algorithm was about 5 m. However, after removing sensor # 4 , the BNorm of the intersection localization algorithm, LS, TLS, and WLS methods was very close, only about 30 m, and the BNorm of BCWLS was also about 5 m. Therefore, according to the results, increasing the number of sensors may not necessarily improve the BNorm performance of target position estimators, and the geometry between the sensors and the target needs to be considered. Form Figure 2b and Figure 5b, it is obvious that after removing sensor # 4 , the RMSE of BCWLS algorithm was also the smallest among the considered algorithms. In addition, after removing sensor # 4 , under the different conditions of angle measurement noise level, the root CRLB of the target position estimates remained almost unchanged compared to that without removal. For example, in the case of σ θ = 4.1 ° , before removing sensor # 4 , the root CRLB of target position estimates was 243.2 m, and after removing sensor # 4 , the root CRLB was 252.1 m, with only a difference of about 8.9 m.
Figure 6 shows the target localization performance of the considered algorithms at a fixed σ θ of 2° and σ s varying from 5 m to 30 m in the case of three sensors. Form Figure 6a and Figure 3a, we can see that the BNorm of the BCWLS algorithm was the smallest among the tested algorithms. In addition, by comparing Figure 6a and Figure 3a, we can see that after removing sensor # 4 , the BNorm of the intersection localization method and LS method decreased approximately from 14 m to 7 m, and the BNorm of the TLS and WLS methods decreased approximately from 10 m to 7 m. Therefore, the results further demonstrate that increasing the number of sensors may not necessarily improve the BNorm performance of target position estimates, and it is necessary to consider the geometry of the sensors and the target. Form Figure 6b and Figure 3b, we can see that the RMSE of these methods was very close before and after removing sensor # 4 . The RMSE of the BCWLS method was also the smallest among the considered algorithms in different cases of sensor self-positioning noise levels. The results verify the superiority of the positioning accuracy of the BCWLS method among the above methods.

3.4. The Computational Cost

According to the analysis in Section 2.3, as the size of the coefficient matrices E and A of the intersection localization method and LS method is 3 × L m and 2 L m × 3 , respectively, the intersection localization method theoretically involves fewer multiplication and addition operations and has a lower computational cost. In order to study the computational cost of the above algorithms, we recorded the required computational time of these algorithms under the same conditions.
Table 2 shows the total time of various algorithms across 20,000 Monte Carlo runs, where the positions of the sensors and the target were set as in Table 1. It illustrates that the time of the intersection localization method was the shortest. Compared with the LS algorithm, the intersection localization algorithm reduced the total time by about 178.16 ms in 20,000 simulations, with a reduction ratio of approximately 36 % . Therefore, the target position estimates of the intersection localization algorithm and the LS algorithm were the same, but the computational cost of the intersection localization algorithm was smaller. Since the target position estimates were obtained by an SVD of the augmented matrix C , the computational cost of the TLS algorithm was higher than that of the LS algorithm. In addition, the time of the BCWLS algorithm increased approximately 385.78 ms compared to that of the WLS algorithm, because the BCWLS algorithm has a bias-compensation operation compared to the WLS algorithm, resulting in a higher computational cost.

4. Discussion

In this paper, we first presented an intersection localization method to estimate the target position based on angle-only measurements. A theoretical analysis showed that this method had the same target position solution as the LS method, but due to the smaller size of the coefficient matrix, the computational cost was lower. The simulation results verified that the intersection localization method and the LS method had the same target position solution and a lower computational cost under the same conditions. Since the size of the coefficient matrices for both the intersection localization and the LS methods is proportional to the number of angle-only measurements, the more angle-only measurements there are, the more computational cost of the intersection localization decreases compared to the LS method.
In addition, through a theoretical analysis, we found that the WLS target localization method with angle-only measurements was biased due to the correlation between the coefficient matrix and the noise vector. Therefore, we further presented a BCWLS method by estimating the bias of the WLS method. This method could reduce the localization bias to a certain extent. The simulation results showed that in different scenarios, the BCWLS method had a smaller BNorm and RMSE compared to the intersection localization method, LS method, TLS method, and WLS method. In addition, the results showed that the BCWLS method could approximately achieve the CRLB. However, due to the addition of a bias-compensation step in the BCWLS method, the computational cost was higher than that of the WLS method.
In order to investigate the impact of the number of sensors on the target positioning accuracy, we removed one sensor from the original sensor network in a simulation and compared the positioning accuracy of various algorithms before and after removal. The results showed that the RMSE and BNorm of various algorithms under different error conditions before and after removal were very close, and the BNorm of various algorithms after removal was even smaller. Therefore, when increasing the number of the sensors to improve the target positioning accuracy, it should be noted that the geometry of the sensors and the target plays an important role. If the placement of new sensors is improper, it is likely to have the opposite effect.

5. Conclusions

This paper studied the source localization problem with angle-only measurements in a passive sensor network. We first presented an intersection localization method that was obtained by minimizing the distances between lines formed by angle-only measurements. Starting from the angle measurement formula, we studied the LS algorithm by solving the target position equations constructed from angle measurements. Comparing the closed-form solutions of the intersection localization algorithm and the LS algorithm, we proved that the two algorithms had the same target position solution. However, since the coefficient matrix of the intersection localization method was smaller, its computational cost was lower than that of the LS method. Furthermore, we studied the TLS method, which takes into account not only the errors in the data vector but also the errors in the coefficient matrix. In contrast, the LS method only considers the error in the data vector. The intersection localization method, LS, and TLS methods do not take into account the difference in measurements errors of different sensors. In practice, the measurement noise level of each sensor may be different. We studied the WLS method by considering the difference in measurement accuracy of sensors. Since the coefficient matrix and the pseudolinear noise vector of the WLS method are correlated, the WLS method is biased, especially at high measurement noise level. To reduce the bias, we presented a BCWLS method by estimating the correlation between the coefficient matrix and pseudolinear noise vector. The BCWLS method had a higher accuracy of the target position estimate than the WLS method. On the other hand, due to the addition of bias-compensation steps, the computational cost of the BCWLS method also increased. Finally, we derived the CRLB of the target localization based on angle-only measurements to evaluate the positioning performance of the above algorithms.
The numerical simulations showed that the intersection localization algorithm and LS algorithm had the same localization results, verifying the theoretical derivation. We also analyzed the impacts of the measurement noise level on the target localization performance of the above target position estimators, including angle measurement noise and sensor self-positioning noise. The target localization performance advantages of the BCWLS method were verified by a numerical simulation. Furthermore, the running time of the considered methods was compared, and the results showed that the intersection localization method had the lowest computational cost.
In this paper, we assumed that the self-positioning noise and the angle measurement noise both followed a zero-mean Gaussian distribution. However, in practice, this assumption may not be valid. In particular, if the passive sensor is mounted on a motion platform, the sensor self-positioning noise may not be zero-mean. Therefore, in the future, we will study the target localization based on angle-only measurements with a non-Gaussian noise background.

Author Contributions

Conceptualization, L.W. and S.Z.; data curation, S.Z.; writing—original draft, L.W., M.G. and P.Z.; writing—review and editing, L.W., J.Y. and X.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work is partially supported by the National Natural Science Foundation of China (62371384).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Olson, E.; Leonard, J.J.; Teller, S. Robust Range-Only Beacon Localization. IEEE J. Ocean. Eng. 2006, 31, 949–958. [Google Scholar] [CrossRef]
  2. Hamdollahzadeh, M.; Amiri, R.; Behnia, F. Moving Target Localization in Bistatic Forward Scatter Radars: Performance Study and Efficient Estimators. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 1582–1594. [Google Scholar] [CrossRef]
  3. Sadeghi, M.; Behnia, F.; Amiri, R.; Farina, A. Target Localization Geometry Gain in Distributed MIMO Radar. IEEE Trans. Signal Process. 2021, 69, 1642–1652. [Google Scholar] [CrossRef]
  4. Famoriji, O.J.; Shongwe, T. Source Localization of EM Waves in the Near-Field of Spherical Antenna Array in the Presence of Unknown Mutual Coupling. Wirel. Commun. Mob. Comput. 2021, 2021, 3237219. [Google Scholar] [CrossRef]
  5. Wang, C.L.; Yang, R.J. Research on Airborne Infrared Passive Location Method Based on Orthogonal Multi-station Triangulation. Laser Infrared 2007, 11, 1184–1187+1191. [Google Scholar]
  6. Kazemi, S.A.R.; Amiri, R.; Behnia, F. Efficient Closed-Form Solution for 3-D Hybrid Localization in Multistatic Radars. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 3886–3895. [Google Scholar] [CrossRef]
  7. Yang, J.; Liu, C.; Huang, J.; Hu, D.; Zhao, C. Overcoming Unknown Measurement Noise Powers in Multistatic Target Localization: A Cyclic Minimization and Joint Estimation Algorithm. Radioengineering 2023, 32, 415–424. [Google Scholar] [CrossRef]
  8. Wang, Y.; Ho, K.C. An Asymptotically Efficient Estimator in Closed-Form for 3-D AOA Localization Using a Sensor Network. IEEE Trans. Wirel. Commun. 2015, 14, 6524–6535. [Google Scholar] [CrossRef]
  9. Bai, G.; Liu, J.; Song, Y.; Zuo, Y. Two-UAV Intersection Localization System Based on the Airborne Optoelectronic Platform. Sensors 2017, 17, 98. [Google Scholar] [CrossRef]
  10. Peng, S.; Zhao, Q.; Ma, Y.; Jiang, J. Research on the technology of cooperative dual-station position based on passive radar system. In Proceedings of the 2020 3rd International Conference on Unmanned Systems (ICUS), Harbin, China, 27–28 November 2020; pp. 611–614. [Google Scholar] [CrossRef]
  11. Chen, Y.; Wang, L.; Zhou, S.; Chen, R. Signal Source Positioning Based on Angle-Only Measurements in Passive Sensor Networks. Sensors 2022, 22, 1554. [Google Scholar] [CrossRef]
  12. Doğançay, K. 3D Pseudolinear Target Motion Analysis From Angle Measurements. IEEE Trans. Signal Process. 2015, 63, 1570–1580. [Google Scholar] [CrossRef]
  13. Nguyen, N.H.; Doğançay, K.; Kuruoglu, E.E. An Iteratively Reweighted Instrumental-Variable Estimator for Robust 3-D AOA Localization in Impulsive Noise. IEEE Trans. Signal Process. 2019, 67, 4795–4808. [Google Scholar] [CrossRef]
  14. Jia, T.; Wang, H.; Shen, X.; Jiang, Z.; He, K. Target localization based on structured total least squares with hybrid TDOA-AOA measurements. Signal Process. 2018, 143, 211–221. [Google Scholar] [CrossRef]
  15. Wells, D.E.; Krakiwsky, E.J. The Method of Least Squares; Department of Surveying Engineering, University of New Brunswick Canada: Fredericton, NB, Canada, 1971; Volume 18. [Google Scholar]
  16. Doğançay, K. On the bias of linear least squares algorithms for passive target localization. Signal Process. 2004, 84, 475–486. [Google Scholar] [CrossRef]
  17. Doğançay, K. Bearings-only target localization using total least squares. Signal Process. 2005, 85, 1695–1710. [Google Scholar] [CrossRef]
  18. Golub, G.H.; van Loan, C.F. An Analysis of the Total Least Squares Problem. SIAM J. Numer. Anal. 1980, 17, 883–893. [Google Scholar] [CrossRef]
  19. Markovsky, I.; Van Huffel, S. Overview of total least-squares methods. Signal Process. 2007, 87, 2283–2302. [Google Scholar] [CrossRef]
  20. Dogancay, K. Relationship Between Geometric Translations and TLS Estimation Bias in Bearings-Only Target Localization. IEEE Trans. Signal Process. 2008, 56, 1005–1017. [Google Scholar] [CrossRef]
  21. Wu, W.; Jiang, J.; Fan, X.; Zhou, Z. Performance analysis of passive location by two airborne platforms with angle-only measurements in WGS-84. Infrared Laser Eng. 2015, 44, 654–661. [Google Scholar]
  22. Frew, E.W. Sensitivity of cooperative target geolocalization to orbit coordination. J. Guid. Control. Dyn. 2008, 31, 1028–1040. [Google Scholar] [CrossRef]
  23. Zhang, Y.; Li, Y.; Qi, G.; Sheng, A. Cooperative Target Localization and Tracking with Incomplete Measurements. Int. J. Distrib. Sens. Netw. 2014, 2014. [Google Scholar] [CrossRef]
  24. Pang, F.; Doğançay, K.; Nguyen, N.H.; Zhang, Q. AOA Pseudolinear Target Motion Analysis in the Presence of Sensor Location Errors. IEEE Trans. Signal Process. 2020, 68, 3385–3399. [Google Scholar] [CrossRef]
  25. Ghaderpour, E.; Antonielli, B.; Bozzano, F.; Scarascia Mugnozza, G.; Mazzanti, P. A fast and robust method for detecting trend turning points in InSAR displacement time series. Comput. Geosci. 2024, 185, 105546. [Google Scholar] [CrossRef]
  26. Hu, J.; Li, Z.; Ding, X.; Zhu, J.; Sun, Q. Spatial–temporal surface deformation of Los Angeles over 2003–2007 from weighted least squares DInSAR. Int. J. Appl. Earth Obs. Geoinf. 2013, 21, 484–492. [Google Scholar] [CrossRef]
  27. Sun, Y.; Ho, K.C.; Wan, Q. Eigenspace Solution for AOA Localization in Modified Polar Representation. IEEE Trans. Signal Process. 2020, 68, 2256–2271. [Google Scholar] [CrossRef]
  28. Shao, H.J.; Zhang, X.P.; Wang, Z. Efficient Closed-Form Algorithms for AOA Based Self-Localization of Sensor Nodes Using Auxiliary Variables. IEEE Trans. Signal Process. 2014, 62, 2580–2594. [Google Scholar] [CrossRef]
  29. Doğançay, K. Passive emitter localization using weighted instrumental variables. Signal Process. 2004, 84, 487–497. [Google Scholar] [CrossRef]
  30. Zhou, S.; Wang, L.; Liu, R.; Chen, Y.; Peng, X.; Xie, X.; Yang, J.; Gao, S.; Shao, X. Signal Source Localization with Long-Term Observations in Distributed Angle-Only Sensors. Sensors 2022, 22, 9655. [Google Scholar] [CrossRef] [PubMed]
  31. Wang, L.; Zhou, S.; Chen, Y.; Xie, X. Bias compensation Kalman filter for 3D angle-only measurements target traking. In Proceedings of the International Conference on Radar Systems (RADAR 2022), Edinburgh, UK, 24–27 October 2022; Volume 2022, pp. 546–551. [Google Scholar] [CrossRef]
  32. Rao, S.K. Pseudo-linear estimator for bearings-only passive target tracking. IEEE Proc. Radar Sonar Navig. 2001, 148, 16–22. [Google Scholar] [CrossRef]
  33. He, S.; Wang, J.; Lin, D. Three-Dimensional Bias-Compensation Pseudomeasurement Kalman Filter for Bearing-Only Measurement. J. Guid. Control. Dyn. 2018, 41, 2678–2686. [Google Scholar] [CrossRef]
  34. Mallick, M.; Morelande, M.; Mihaylova, L.; Arulampalam, S.; Yan, Y. Angle-Only Filtering in Three Dimensions; Wiley: Hoboken, NJ, USA, 2016; pp. 3–42. [Google Scholar]
  35. Arasaratnam, I.; Haykin, S. Cubature Kalman Filters. IEEE Trans. Autom. Control. 2009, 54, 1254–1269. [Google Scholar] [CrossRef]
  36. Julier, S.; Uhlmann, J. Unscented filtering and nonlinear estimation. Proc. IEEE 2004, 92, 401–422. [Google Scholar] [CrossRef]
  37. Mao, W.L. GPS Interference Mitigation Using Derivative-free Kalman Filter-based RNN. Radioengineering 2016, 25, 518–526. [Google Scholar] [CrossRef]
  38. Yu, J.Y.; Coates, M.J.; Rabbat, M.G.; Blouin, S. A Distributed Particle Filter for Bearings-Only Tracking on Spherical Surfaces. IEEE Signal Process. Lett. 2016, 23, 326–330. [Google Scholar] [CrossRef]
Figure 1. Three-dimensional target localization scenario.
Figure 1. Three-dimensional target localization scenario.
Remotesensing 16 01319 g001
Figure 2. The BNorm and RMSE of the target position estimators with σ s fixed at 10 m and σ θ varying from 0.1° to 6.1°.
Figure 2. The BNorm and RMSE of the target position estimators with σ s fixed at 10 m and σ θ varying from 0.1° to 6.1°.
Remotesensing 16 01319 g002
Figure 3. The BNorm and RMSE of the target position estimators with σ θ fixed at 2° and σ s varying from 5 m to 30 m.
Figure 3. The BNorm and RMSE of the target position estimators with σ θ fixed at 2° and σ s varying from 5 m to 30 m.
Remotesensing 16 01319 g003
Figure 4. The target position estimates of the WLS and BCWLS algorithms on the X-O-Y plane and X-O-Z plane. ◯ represents the WLS estimates of the target position. □ represents the BCWLS estimates of the target position. Δ represents the true position of the target. ✩ represents the mean of target position estimates.
Figure 4. The target position estimates of the WLS and BCWLS algorithms on the X-O-Y plane and X-O-Z plane. ◯ represents the WLS estimates of the target position. □ represents the BCWLS estimates of the target position. Δ represents the true position of the target. ✩ represents the mean of target position estimates.
Remotesensing 16 01319 g004
Figure 5. The BNorm and RMSE of the target position estimators with σ s fixed at 10 m and σ θ varying from 0.1° to 6.1° in the case of 3 sensors.
Figure 5. The BNorm and RMSE of the target position estimators with σ s fixed at 10 m and σ θ varying from 0.1° to 6.1° in the case of 3 sensors.
Remotesensing 16 01319 g005
Figure 6. The BNorm and RMSE of the target position estimators with σ θ fixed at 2° and σ s varying from 5 m to 30 m in the case of 3 sensors.
Figure 6. The BNorm and RMSE of the target position estimators with σ θ fixed at 2° and σ s varying from 5 m to 30 m in the case of 3 sensors.
Remotesensing 16 01319 g006
Table 1. Positions of sensors and target.
Table 1. Positions of sensors and target.
Position (m)
Sensor #1 p 1 o = [ 0 , 0 , 0 ] T
Sensor #2 p 2 o = [ 3000 , 3000 , 0 ] T
Sensor #3 p 3 o = [ 1000 , 1000 , 0 ] T
Sensor #4 p 4 o = [ 1000 , 1000 , 0 ] T
Target #1 g 1 o = [ 1000 , 2000 , 1000 ] T
Table 2. The total computation time of the intersection localization method, LS, TLS, WLS, and BCWLS algorithms across 20,000 Monte Carlo runs, absolute and relative.
Table 2. The total computation time of the intersection localization method, LS, TLS, WLS, and BCWLS algorithms across 20,000 Monte Carlo runs, absolute and relative.
MethodIntersectionLSTLSWLSBCWLS
Absolute time (ms)312.05490.21675.601130.041515.82
Relative time0.6411.382.313.09
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

Wang, L.; Zhou, S.; Gong, M.; Zhao, P.; Yang, J.; Sui, X. Three-Dimensional Signal Source Localization with Angle-Only Measurements in Passive Sensor Networks. Remote Sens. 2024, 16, 1319. https://doi.org/10.3390/rs16081319

AMA Style

Wang L, Zhou S, Gong M, Zhao P, Yang J, Sui X. Three-Dimensional Signal Source Localization with Angle-Only Measurements in Passive Sensor Networks. Remote Sensing. 2024; 16(8):1319. https://doi.org/10.3390/rs16081319

Chicago/Turabian Style

Wang, Linhai, Shenghua Zhou, Min Gong, Pengfei Zhao, Jian Yang, and Xin Sui. 2024. "Three-Dimensional Signal Source Localization with Angle-Only Measurements in Passive Sensor Networks" Remote Sensing 16, no. 8: 1319. https://doi.org/10.3390/rs16081319

APA Style

Wang, L., Zhou, S., Gong, M., Zhao, P., Yang, J., & Sui, X. (2024). Three-Dimensional Signal Source Localization with Angle-Only Measurements in Passive Sensor Networks. Remote Sensing, 16(8), 1319. https://doi.org/10.3390/rs16081319

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