Next Article in Journal
Integrating Metal-Oxide-Decorated CNT Networks with a CMOS Readout in a Gas Sensor
Next Article in Special Issue
Sensorial Systems Applied to Intelligent Spaces
Previous Article in Journal
Dual Super-Systolic Core for Real-Time Reconstructive Algorithms of High-Resolution Radar/SAR Imaging Systems
Previous Article in Special Issue
Ultrasonic Array for Obstacle Detection Based on CDMA with Kasami Codes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Error Estimation for the Linearized Auto-Localization Algorithm

Centro de Automática y Robótica (CAR), Consejo Superior de Investigaciones Científicas (CSIC)-UPM, Ctra. Campo Real km 0.2, La Poveda-Arganda del Rey, 28500, Madrid, Spain
*
Author to whom correspondence should be addressed.
Sensors 2012, 12(3), 2561-2581; https://doi.org/10.3390/s120302561
Submission received: 6 January 2012 / Revised: 17 February 2012 / Accepted: 20 February 2012 / Published: 24 February 2012
(This article belongs to the Special Issue Sensorial Systems Applied to Intelligent Spaces)

Abstract

: The Linearized Auto-Localization (LAL) algorithm estimates the position of beacon nodes in Local Positioning Systems (LPSs), using only the distance measurements to a mobile node whose position is also unknown. The LAL algorithm calculates the inter-beacon distances, used for the estimation of the beacons’ positions, from the linearized trilateration equations. In this paper we propose a method to estimate the propagation of the errors of the inter-beacon distances obtained with the LAL algorithm, based on a first order Taylor approximation of the equations. Since the method depends on such approximation, a confidence parameter τ is defined to measure the reliability of the estimated error. Field evaluations showed that by applying this information to an improved weighted-based auto-localization algorithm (WLAL), the standard deviation of the inter-beacon distances can be improved by more than 30% on average with respect to the original LAL method.

1. Introduction

There are many location-aware applications that require the estimation of the position of persons or objects in indoor environments. Since the Global Positioning System (GPS) is not available inside buildings, several localization systems have been designed to work in indoor environments, which are commonly known as Local Positioning Systems (LPSs) [1]. These systems require the installation of several nodes at fixed positions (called beacons) in the indoor environment, whose positions must be known in advance. The location of a mobile node (attached to the person or object to locate) can be obtained by the trilateration method using the measured distances between the beacons and the mobile node.

The manual determination of the beacons’ positions (e.g., using measuring tapes or ultrasonic/laser rangers) is a cumbersome error-prone method. Therefore various techniques have been proposed to address the problem of obtaining the position of the beacons [211], sometimes called the auto-calibration or auto-localization problem.

Typical auto-localization solutions are based on measuring distances from the beacons to a group of localized nodes, or a mobile node at several known positions, and performing an inverse trilateration (locating the beacons using the mobile node) [2]. In [4] four different mobile node positions are used to locate the beacons of a 3D LPS. In [3] three nodes with known positions are required plus a group of nodes with unknown positions. In [5] only the relative distances between four nodes mounted on a mobile robot are required. These methods, however, require an external localization system to obtain the position of the mobile node at each static location, which is not always available in indoor environments.

Assuming that no information on the position of any node is known, the only available data are the distance measurements between the beacons and the mobile node. Using this information Duff and Muller [6] proposed a nonlinear least-square optimization algorithm, where the objective functions were the distance equations and the variables the coordinates of all nodes. The initial conditions for the algorithm (a first position estimation of all nodes) were obtained by a trial and error method, by randomly generating those conditions and choosing the best solution. In [7] and [8] an Extended Kalman Filter and an H-infinite filter are used, respectively, to obtain the position of the beacons and the mobile node. In both cases an initial position estimation for the beacons is obtained using dead-reckoning information of the mobile node. In [9], a distance matrix is formed with the range measurements between beacons and the mobile node at different locations. With that information a rough approximation of the inter-beacon ranges is obtained using an interpolation scheme.

In [10] a solution of the auto-localization problem, based on the linearization of the trilateration equations, was presented. The method, known as the Linearized Auto-Localization (LAL) algorithm, neither requires a trial and error approximation (i.e., randomly generated positions) nor any external positioning information (such as dead-reckoning data). The LAL algorithm is used when all beacons are located in a plane (e.g., when all the beacons are located on the ceiling) parallel to the plane containing the mobile node trajectory. By using these conditions, all nonlinear terms of the trilateration equations can be grouped in additional variables in order to obtain linear auto-localization equations.

In the present paper we propose a method to estimate the solution’s standard deviation of the LAL algorithm based on a first order Taylor approximation. Since the method depends on the assumption that the solution has a Gaussian distribution, a confidence parameter τ is defined to measure the reliability of the estimated standard deviation. The information obtained from the solution’s precision is later used to further improve the LAL algorithm.

This paper is organized as follows. Section 2 describes the LAL algorithm [10]. In Section 3 the equations to estimate the solution’s standard deviation and their reliability are proposed. These are evaluated by simulation in Section 4 and also experimentally on an ultrasonic 3D LPS system in Section 6. Section 7 presents our conclusions.

2. The LAL Algorithm

In [10] a new solution for the auto-localization problem of 3D LPSs was presented based on the distance measurements between beacons and the mobile node. While other techniques require an external positioning system (e.g., a second LPS or an odometric sensor) to estimate the positions of the mobile node, the proposed method uses only the measurements available to the LPS being used for mobile node location. The method is based on the linearization of the trilateration equations by grouping nonlinear terms in additional variables. The LAL algorithm is defined for the case where all the beacons are located in a plane parallel to the plane containing the mobile node trajectory. For this particular configuration, a solvable initial subset of three beacon nodes is obtained based on [11]. Figure 1 shows the auto-localization configuration for n = 3 beacon nodes Ni, i = {1, 2, 3} and m measurement points Nj, j = {4, . . ., m + 3} from the mobile node path. For n = 3 beacons nodes a minimum of m = 6 measurements points are required [11], though more measurement points can be added to improve the estimated solution. From now on, the measurement points of the mobile node will be referred to as virtual nodes.

For a configuration of n = 3 beacons and m = 6 virtual nodes, the trilateration equations can be written as a function of two groups of distances: the unknown inter-beacon distances DU = {d12, d13, d23} (the variables to be estimated) and the measured distances DM = {d14, d24, d34, . . ., d39} between beacon nodes and virtual nodes (the available data). The final equation is expressed in the linear form [10]:

AX = B
where:
A = [ d 34 2 d 39 2 d 24 2 d 29 2 d 14 2 d 19 2 E 324 E 214 E 329 E 219 E 314 E 214 + E 319 E 219 d 35 2 d 39 2 d 25 2 d 29 2 d 15 2 d 19 2 E 325 E 215 E 329 E 219 E 315 E 215 + E 319 E 219 d 38 2 d 39 2 d 28 2 d 29 2 d 18 2 d 19 2 E 328 E 218 E 329 E 219 E 318 E 218 + E 319 E 219 ]
E ikj = d ij 2 d kj 2 ,     i , k = { 1 , 2 , 3 }     j = { 4 , 5 , , 9 }
X = 1 d 12 2   [ d 12 2 ( d 13 2 + d 23 2 d 12 2 ) d 13 2 ( d 12 2 + d 23 2 d 13 2 ) d 23 2 ( d 12 2 + d 13 2 d 23 2 ) d 13 2 d 23 2 ]
B = [ E 314   E 324 E 319   E 329 E 315   E 325 E 319   E 329 E 318   E 328 E 319   E 329 ]

Matrix A ∈ ℝ(m−1)×5 and vector B ∈ ℝ(m−1)×1 are composed of the distance measurements from the virtual nodes DM. The vector X ∈ ℝ5×1 includes the unknown inter-beacon distances. Once X is obtained the inter-beacon distances solution DU can be calculated with:

d 12 = 1 2   X 2 X 4 + X 3 X 5 d 13 = 1 2   X 1 + X 3 X 5 d 23 = 1 2   X 1 + X 2 X 4

For more than three beacon nodes an incremental procedure is used. The mobile node is moved on a plane trying to obtain at least six measurements shared by subsets of three beacon nodes. The goal is to obtain all possible inter-beacon solutions DU applying the linearized Equation (1) on such subsets. In large areas, or if the range of the nodes is limited, multiple paths can be required. Once enough inter-beacon distances have been obtained, a 2D localization algorithm is used to locate every beacon position. In the LAL method an algorithm based on the Multidimensional Scaling (MDS) is used to obtain all the beacons’ positions [12]. The MDS algorithm requires the distance matrix which represents the pairwise distances between all the beacons. Since for large localization areas it is difficult to move the mobile node in a path where all beacons are in range, some inter-beacon distances must be estimated by other means like the nearest path between those beacons. However, this estimation is a rough approximation of the real inter-beacons distances and it has a negative impact in the precision of the MDS method. To avoid this problem, the LAL algorithm uses a modification of the MDS proposed in [13] known as LaMSM. This method locates only subsets of fully connected beacons forming local maps that are later merged in a global map containing all beacons.

3. Precision of the Auto-Localization Method

Consider that the distance measurements DM, between beacons and virtual nodes, are given by:

d ij = d ij ¯ + ɛ ij
where d ij ¯ is the true distance between beacon i = {1, 2, 3} and virtual node j = {4, . . ., m + 3}, and ɛij is a zero mean Gaussian distribution with variance σ ij 2. An estimation of the precision of the inter-beacon distances DU obtained by the LAL algorithm presented in Section 2 would be very useful in the posterior determination of the beacons’ location. However, since it is clear that the resultant error distribution of DU is not Gaussian, a direct evaluation of such error can be a complex task. A good estimation can be obtained by using a first order approximation of the nonlinear equations based on a technique called differential sensitivity analysis [14].

3.1. Differential Sensitivity Analysis

Differential Sensitivity Analysis (DSA) is a technique used to evaluate the error of a given function originated by the uncertainty present on its variables. This technique uses a first order Taylor expansion to obtain an approximation of the function’s variance.

Consider l variables yu = gu(Z), u = {1, .., l}, each one function of p uncorrelated variables Z = {z1, . . ., zp} normally distributed with mean k and variance σ z k 2. The variance V (yu) of each variable yu can be expressed as:

V ( y u ) = k = 1 p ( g u z k ) 2   σ z k 2
and the covariance C(yu, yv) between variables yu and yv:
C ( y u ,   y v ) = k = 1 p ( g u z k )     ( g v z k )     σ z k 2

Equations (8) and (9) can be written in matrix form to obtain the covariance matrix C(Y) of vector Y = [y1, . . ., yl]T [15]:

C ( Y ) = 𝒢 C ( Z ) 𝒢 T
where C(Z) is the covariance matrix of Z and 𝒢 is the Jacobian matrix of G = {g1, . . ., gl} with respect to vector Z.

The approximation obtained with the DSA is only valid if the real distribution of Y is close to a Gaussian distribution, that is, if the higher terms of the Taylor series can be neglected. It is important to verify the range of application of the method, otherwise the estimated parameters could be significantly different from the true ones. The range of validity of this approximation will be further discussed in the present section.

In the next subsection we will apply the DSA to obtain the variance of the least squares solution of Equation (1).

3.2. Error Perturbation in Least Squares Solutions

For the estimation of the error in the least squares solution, we will use an approximation proposed by Stewart [16]. Similarly to the differential technique previously used, it is based on the assumption that the error can be modelled as a Gaussian distribution. Although in [16] this approximation was obtained for the case when the components of matrix A and vector B are row-wise independent, we extend it to a more general case where all elements of A and B can be dependent.

Let us define the observed matrices A′ = A+ à and B′ = B + , where à and represent the error perturbation on the true matrices A and B, originated by the noisy distance measurements. Matrices A′ and B′ are the ones used to obtain the least squares solution, since we do not know the actual values A and B of such matrices. Equation (1) can be rewritten as:

( A A ˜ ) X = B B ˜
left multiplying by the pseudo-inverse A+ = (A′T A′)−1 AT and rearranging the equation we obtain:
X = X + A + ( B ˜ A ˜ X )
where X′ = A+ B′ is the calculated solution we want to evaluate. The pseudo-inverse A+ can be expressed as a series function of matrix A [16]:
X = X + A + ( B ˜ A ˜ X ) + F ( B ˜ A ˜ X )
where F represents the higher order terms of A+.

Finally we can simplify Equation (13) to obtain the approximation:

X = X + A + ( B ˜ A ˜ X )
where the step from Equation (13) to Equation (14) is based on the hypothesis that A+F. Notice that Equation (14) is not a first order approximation of Equation (13), but rather a more restricted one.

Using Equation (14) we can obtain the desired covariance matrix C(X′) of the least squares solution X′, which is given by:

C ( X ) = C ( A + ( B ˜ A ˜ X ) ) = A + C ( B ˜ A ˜ X ) A + T = A + C ( D M ) T A + T
where the is the Jacobian matrix of R = BAX with respect to DM.

3.3. Covariance Matrix of the Inter-Beacon Distances

The covariance matrix C(DU) of the inter-beacon distances DU is obtained by applying Equation (10) to the distance Equation (6):

C ( D U ) = C ( X ) T
where represents the Jacobian matrix of Equation (6) with respect to X.

Replacing Equation (15) in Equation (16) we obtain the final variance error of the LAL method:

C ( D U ) = A + C ( D M ) T A + T T

Considering the simplest case, where all the measured distances DM have the same variance σ ij 2 = σ M 2, i = {1, 2, 3} j = {4, . . ., m + 3} we have that Equation (17) can be simplified as:

C ( D U ) = σ M 2 A + T A + T T = σ M 2   M DOP
where the effect of matrix MDOP ∈ ℝ3×3 is equivalent to the effect of the geometric dilution of precision (GDOP) in GPS [17]. For our case, we define the Distance Dilution of Precision (DDOP) as the amplification of the standard deviation σM of the measurement errorsDM onto the inter-beacon distances DU. This amplification only depends on the relative position of the nodes (beacons and virtuals).

For each inter-beacon distance dik i, k = {1, 2, 3}, between beacons i and k, a respective distance dilution of precision ddopik can be obtained by:

[ ddop 12 ddop 13 ddop 23 ] = 1 σ M   [ σ 12 σ 13 σ 23 ] = [ m DOP 11 m DOP 22 m DOP 33 ]
where σik is the standard deviation of the distance dik and mDOPuv denotes the (uv) element of matrix MDOP .

3.4. Reliability of the Variance Estimation

As discussed in Section 3.1, the variance estimation of the errors of DU is only valid if its distribution can be approximated to a zero mean Gaussian distribution. Since such condition depends on several variables, a parameter to measure the validity of the estimation is necessary. For example, a second order Taylor series approximation could be used to evaluate higher order moments of the error distribution, such as the skewness and kurtosis [18]. In practice, however, we found that evaluating only the error perturbation on the least squares solution (Section 3.2) is enough to verify the reliability of the estimated variance. This is because a more restricted approximation is used for the pseudo-inverse linearization than for the other non-linear equations.

The first order Taylor series of the pseudo-inverse of matrix A′ = A + Ã can be expressed as [19]:

A + A + A +   A ˜ A + + ( A T   A ) 1 A ˜ T   P A + + F
where P = IAA+ is the complementary orthogonal projection onto the orthogonal space of A. Here F represents the effect of noise over matrix A′. Since the approximation used in Equation (14) is valid when A+F, a comparison between these terms is proposed in [16] by using a parameter τ that meets the condition:
τ F S A + F
where ||·||F represents the Frobenius norm and ||·||S the Stochastic norm as defined in [19]. The parameter τ represents the ratio between the error perturbation F, associated by the noisy measurements, and the real value of matrix A+. As the τ value increases, the perturbation became increasingly important and can not be disregarded.

In Appendix 7 it is shown that, assuming the general case where all elements of A and B can be dependent, τ can be expressed as:

τ = max { m 1 , 2 × 5 }   || Q 1 2   A + || 2
where the elements quv of matrix Q ∈ ℝ5×5 are obtained from:
q uv = k = 1 m 1 C ( a ˜ ku ,   a ˜ kv )
being C (ãku, ãkv) the covariance between elements u and v of the kth row of matrix Ã. The covariance of these elements can be calculated by applying Equation (9) to the elements of matrix A defined in Equation (2).

4. Evaluation of the DSA Method

In order to evaluate the performance of the DSA method to predict the precision of the LAL inter-beacon solution, a LPS was simulated based on the node configuration shown in Figure 2. The measurement points are distributed along a circular path under the beacons, plus one measurement point at its center. In the following simulations, unless otherwise stated, the parameters are: radius r = 3 m, number of virtual nodes m = 12 and height h = 2m. The ranging data was generated with an additive white Gaussian noise with zero mean and a standard deviation of 0.01 m. All simulations were performed 1,000 times.

To avoid any confusion, we here define some terms used to refer to the different types of distance errors present in this section. We will refer as input errors those associated with the distance measurements DM. These errors have a known distribution (zero mean Gaussian noise) and a known standard deviation. The output errors are associated with the unknown inter-beacon distances DU obtained using the LAL algorithm. We call simulated values those obtained by simulating the auto-localization problem 1,000 times in order to calculate the output error statistics. In contrast the estimated values are those obtained using the methods and equations defined in this paper (Equations (18) and (22)). Finally offline estimates are those solutions obtained with the actual values of the distance measurements (i.e., when we know the exact position of all nodes). The online estimates are solutions obtained with noisy distance measurements and the estimated position of the nodes obtained by the auto-localization method.

The offline estimation allows us to evaluate the limits of the DSA method under ideal conditions (i.e., using the noise-free values of matrices A and X and the vector B). The estimation will be reliable as long as the output distribution of distances DU remains approximately Gaussian. As proposed in Section 3.4, the τ value should be sufficient to verify this reliability. The online method allows us to estimate “on-the-fly” the standard deviation of distances DU and use this information to improve the beacons’ localization. This method will present the same limits than the offline estimation, plus the effects of using the noisy matrices A′ and X′ and vector B′. The effect of the noisy matrices is evaluated by performing 1,000 times the online method on each test. The τ value will also be sufficient to verify the reliability of this method.

4.1. Standard Deviation Estimation Analysis

To evaluate the DSA method we calculated the quadratic mean output standard deviation σ mean = 1 3 σ 12 2 + σ 13 2 + σ 23 2, of the inter-beacon distances DU, obtained when shifting the center of the virtual nodes path of the LPS configuration (Figure 2). The resultant standard deviation maps obtained by simulation (Figure 3(a)) and by offline estimation (Figure 3(b)) match almost exactly. The advantage of the offline estimation over the simulation is that the calculation processes of the maps are computationally more efficient. The evaluation of the standard deviation maps can be very useful to analyse the ideal route path of the virtual nodes. For example, in Figure 3(b) it is shown that the ideal position of the center of the path is near the central point of the beacons (1.5 m,1.3 m). If the center of the path is far from the central point, σmean can increase up to 8 cm, that is, eight times the input standard deviation.

Figure 3(c) shows the estimated standard deviation obtained by the offline and online methods compared to the one obtained by simulation. As expected the online estimates’ errors are higher than the ones obtained by the offline method. Also, since the calculation of the vector X is directly related with the output standard deviation (Equation (16)), the online estimates worsens as the output standard deviation increases.

4.2. Reliability of the Variance Estimation

In Section 3.4 we proposed that the parameter τ can be used to weight the validity of the estimated output noise. To evaluate this premise we perform a simulation changing the two parameters that affect the output variance: the DDOP and the input variance. To change the DDOP we use the same LPS configuration as before but change the radius to the values r = {0.8 m, 1.2 m, 2 m} which produce a ddop12 = {7.79, 3.57, 1.52}. For the input variance a standard deviation ranging from 0 to 0.21 m is used.

In Figure 4(a) the simulated and estimated output standard deviation of distance d12 is shown with respect to the input standard deviation. For the online estimates the mean and the 95% and 5% percentile are shown. The graph corresponds to a ddop12 of 1.52. For an input noise less than 0.1 m, the output standard deviation increments almost linearly with the input standard deviation, as predicted by Equation (18). For higher values the simulation shows that the simulated standard deviation is always higher than the obtained from Equation (18). This happens because for high input noises the assumption that the error distribution of the output distances is close to a zero mean Gaussian distribution is no longer valid. This can be verified in Figure 4(b), where it is shown how the skewness of the distance d12 increases with the input noise. Since the skewness calculated for d12 is positive, the estimated standard deviation will be always lower than the one obtained in the simulation. This shows that the estimated offline standard deviation should not be used for a high input noise. For example, a 0.1 m input noise limit can be chose for this particular LPS configuration. Figure 4(a) also shows, as expected, that the estimated online standard deviation presents a higher error with the increment of the input noise.

Figure 5 shows the offline and 95% percentile online standard deviation error in percentage (compared to the one obtained by simulation) as a function of the τ value. These results are obtained for a ddop12 = {7.79, 3.57, 1.52} and a standard deviation input error ranging from 0 to 0.21 m. The graphics are limited to values of 0 ≤ τ ≤ 1.5 to focus within a reasonable range, since for values of τ > 1.5 the errors are much higher than the ones shown in the graphics.

Figure 5(a) shows the standard deviation error in percentage obtained with the offline method. This error is always under 10% for τ values inferior to 0.7, and within this range it seems almost constant regardless of the value of τ. For τ > 0.7 the error increases with τ . It is clear that, at this point, the effects of the skewness of the output distances begin to be appreciable. Notice that the value of τ where the standard deviation error begins to increase seems to be independent of the ddop12 value. This shows that the τ value is sufficient to verify the reliability of Equation (18) for the offline method, regardless of the input errors or the geometry of the LPS system.

A similar analysis used for the offline estimation can be used for the online case. Figure 5(b) shows the 95% percentile of the standard deviation error obtained with the online method in percentage. Since for every of the 1,000 iterations a slightly different value of τ is obtained (due to noisy matrices) we grouped τ in fixed intervals of 0.05 values (e.g., any value of τ between 1 and 1.05 is represented in the graph as a value of τ = 1). It can be seen that for τ < 0.7 the standard deviation error is always under 40%, though within this range the error is not constant and always increases with the value of τ . This effect is generated by the noise perturbation over the matrices A′ and B′. Still, we can use τ to limit the maximum error obtained when estimating the output standard deviation. For example, to obtain errors not higher than 20% a value of τ ≤ 0.3 can be chosen.

In this case, the value τ where the standard deviation error is lower than a given limit is not totally independent of the dop12 value, so it is possible that for higher values of DDOP or input noise the error could be higher. However, in practical applications, the values of the DDOP, the input noise, and τ will mostly be lower than the ones here simulated. For example, the value of τ obtained in the tests performed in Section 4.1 was always of τ ≤ 0.13

5. Weighted Linearized Auto-Localization Algorithm (WLAL)

In this section a new algorithm of auto-localization is presented based on the LAL algorithm and the error estimation techniques exposed in Section 3. As presented in [10], and reviewed in Section 2, the LAL algorithm calculates the distances between the unlocated beacons that are later used to obtain their relative position. The algorithm calculates the inter-beacon distances of every subset of three beacons, therefore more than one distance estimation is usually obtained between any pair of beacons. The final distance between two beacons is obtained by calculating the mean of all the available distances. We propose to rather use a weighted mean based on the online estimated standard deviations of the distances using the method presented in Section 3. If for each inter-beacon distance dij we obtain k distances estimations, the respective weighted mean distance ij is obtained by:

d ¯ ij = k = 1 n ( d ij ) k / ( σ ij 2 ) k k = 1 n 1 / ( σ ij 2 ) k
where (dij)k is the k-th distance dij estimation and (σij)k the respective standard deviation that can be calculated from Equation (19).

By using the estimated standard deviation for each possible subset, we can give more weight to solutions with the best subset of beacons. With the WLAL algorithm we expect an improvement of the resultant standard deviation of the inter-beacon distances. The expected variance σ̄ij for each interbeacon distance ij obtained by applying the weighted mean is:

σ ¯ ij 2 = 1 k = 1 n 1 / ( σ ij 2 ) k

6. Evaluation on an Ultrasonic LPS

We apply the DSA method on a real ultrasonic LPS, where the assumptions of co-planarity between beacons and the presence of only Gaussian noise in the measurements are not completely true. We want to evaluate the performance of the method under non-ideal conditions. The test is performed on the 3DLocus system [4], shown in Figure 6(a), which is an acoustic LPS composed by n = 7 beacons deployed on a cell of 2.8 m × 2.8 m × 2.8 m. The calculated standard deviation of the distance measurements obtained with the 3DLocus is 0.23 mm, which is much more accurate than the one used in the simulations in Section 4. The height between the beacons and the virtual nodes is 1.4 m.

The linearized auto-localization method is compared with the inverse positioning method presented in [2]. Since the inverse positioning method requires the exact position of the mobile node, a Staübli Unimation industrial robotic arm with a 50 μm accuracy was used for positioning such node. An approximately circular path with a radius of 0.8 m of m = 9 virtual nodes were used in this test, and a total of 100 measurements were made on each point. From these measurements we obtain 100 individual estimations of the positions of the beacons. The nodes’ configuration used in this test is shown in Figure 6(b).

6.1. Inter-Beacon Distances Estimation

Figure 7(a) shows the estimated online and measured distance dilution of precision ddop12 obtained with different subsets of beacons. The measured ddop12 was calculated by dividing the standard deviation of the calculated distance d12 by the standard deviation of the distance measurements DM. In most cases the estimated value is slightly lower than the actual ddop12 obtained in the 3DLocus, though the difference is always below 30%. One reason for this difference is that we are using a mean standard deviation of 0.23 mm for all the 3DLocus measurements while the actual value of the measurements’ errors depends on various factors such as the distance and angle between nodes. In Figure 7(b) the histogram of the calculated distance d12 obtained using the beacons’ subset {1, 2, 3} is shown. As can be observed, the error distribution of the calculated distance resembles a Gaussian distribution. In order to verify the assumption of a Gaussian error distribution on all the calculated inter-beacon distances, we run a Lilliefors test for normality [20]. The test established with a p-value of 0.05 that the error distribution on the calculated inter-beacon distances approximates a normal distribution (i.e., there is a 5% probability that the normality is a false positive).

In Table 1 a comparison of the output standard deviation σl obtained with the LAL algorithm and the standard deviation σw obtained with the WLAL algorithm are shown. In every case σw was lower than σl. For example, for the distance d12 the standard deviation obtained by the WLAL algorithm is 65% less than the one obtained with the LAL algorithm. Since we know that the ideal position of the virtual nodes path is under the center of the triangle formed by the subset of three beacons, it is clear that the solution obtained with the subset {1, 2, 3} will have a higher error than the one obtained with the subset {1, 2, 5}, when estimating the d12 distance. On average, an improvement of 32.8% was obtained using the WLAL algorithm.

6.2. Beacons’ Position Estimation

In Figure 8 the beacons’ position estimation obtained over one hundred trials with the LAL method, WLAL algorithm and the inverse positioning method is shown on the X-Y plane. A coordinate system was defined using the beacon 1 as origin and the beacon 2 as the X axis. As expected the standard deviation of the beacons’ position is improved by using the weighted mean instead of a simple mean of the estimated inter-beacon distances. On average, an improvement of 22% was obtained on the standard deviation on axis X and Y using the WLAL algorithm.

Figure 8 also shows bias errors between the auto-localization methods. Using the inverse positioning solution as the ground truth location of the beacons, a mean RMS error of the beacons’ positions of 0.76 cm and 0.92 cm is obtained with the LAL and WLAL algorithms respectively. The greatest bias observed is for the fifth beacon (1.3 cm). The presence of a bias between the different methods could be originated by several causes. First, although the cell structure of the 3DLocus is positioned parallel to the floor, it seems that there is a small inclination that causes a height variation on the beacons. The inverse positioning method shows a mean height variation of up to 8 mm between the beacons. Since the LAL and WLAL methods assume that the beacons are on a plane, the height difference will generate a discrepancy between the inverse positioning and the linearized methods. A second source of the observed bias could be the noisy distances used in the solution of the inter-beacon distances. In this paper we are only taking into account the influence of the noise using a first order approximation, however, it is known that the higher terms of the noise can also generate a bias [21].

To evaluate these hypotheses, we simulated a LPS using the positions obtained with the inverse positioning as the real location of the beacons, but positioning them in the same plane (by setting the Z coordinate to be zero on all beacons) and using simulated measurements with the same standard deviation. The simulation showed no bias between the methods in any beacon, as can be seen for the fifth beacon in Figure 9, showing that the effect of the higher terms of the noise is negligible. When the beacons were not placed in the same plane, the simulation showed a bias between all the auto-localization methods. Using the beacon’s altitude obtained with the inverse positioning, the simulation showed a bias of 5 mm for the fifth beacon. Finally, systematic range errors originated from the actual LPS can also add to the bias observed [22] (e.g., the range measurement is affected by the orientation between the beacon and the mobile node). Any error originated by this bias and also the small height variation between beacons can be compensated using an optimization algorithm such as the ones used in [68]. The WLAL solution could be used as a first estimation of the beacons position, since with a 1.5 cm offset any optimization algorithm will easily converge to a more accurate solution.

7. Conclusions

In this paper a method to estimate the error of the solutions obtained with the Linearized Auto-Localization (LAL) method was presented. The method is based on the differential sensitivity analysis that uses a first order Taylor approximation to obtain the function’s error variance. Since the method depends on such approximation, a confidence parameter τ was defined to measure the reliability of the estimated error distribution. The differential sensitivity analysis showed that the standard deviation of the solution, obtained by the LAL method, is proportional to the standard deviation of the distance measurements and a matrixMDOP dependent of the geometry of the localization system, which is similar to the geometric dilution of precision on localization problems.

Two versions of the error estimation were evaluated: (a) the offline method that requires the location of all nodes and (b) the online method where all the nodes’ positions are unknown. The first method is useful to evaluate possible path routes and strategies that maximize the accuracy of the auto-localization algorithm. The advantage of this method is that the calculation processes are computationally more efficient than an evaluation based on a simulation process. Simulated test showed that τ can be used to limit the maximum error obtained by the offline method regardless of the measured distances noise and the DDOP. The online method can be used to estimate “on-the-fly” the standard deviation of the inter-beacon distances obtained by the auto-localization algorithm, and use this information to improve its accuracy. It was also shown by simulation that the parameter τ can be used to limit the maximum error obtained with the online method.

Finally, a modification of the LAL method was used to evaluate the online method on an acoustic LPS. For the calculation of the inter-beacons distances, a weighted mean rather than a simple mean was proposed by using the output error statistics obtained by the online method. By using this Weighted Linearized Auto-Localization (WLAL) algorithm, the standard deviation of the calculated inter-beacon distances showed an improvement of 32.8% on average. The beacons’ position estimation also presented a improvement of 22% of the standard deviation on axes X and Y.

Acknowledgments

The authors would like to thank the financial support provided by projects LEMUR (TIN2009-14114-C04-03) and LOCA (CSIC-PIE Ref.200450 E430).

References

  1. Hightower, J.; Borriello, G. Location systems for ubiquitous computing. Computer 2001, 34, 57–66. [Google Scholar]
  2. Mahajan, A.; Figueroa, F. An automatic self-installation and calibration method for a 3D position sensing system using ultrasonics. Robot. Auton. Syst 1999, 28, 281–294. [Google Scholar]
  3. Ruiz, F.D.; Ureña, J.; Villadangos, J.M.; Gude, I.; García, J.J.; Hernández, Á.; Jiménez, A. Optimal Test-Point Positions for Calibrating an Ultrasonic LPS System. Proceedings of the IEEE International Conference on Emerging Technologies and Factory Automation (ETFA ’08), Hamburg, Germany, 15–18 September 2008; pp. 338–344.
  4. Prieto, J.C.; Jiménez, A.R.; Guevara, J.; Ealo, J.; Seco, F.; Roa, J.; Ramos, F. Performance evaluation of 3D-LOCUS advanced acoustic LPS. IEEE Trans. Instrum. Meas 2009, 58, 2385–2395. [Google Scholar]
  5. Nishitani, A.; Nishida, Y.; Hori, T.; Mizoguchi, H. Portable Ultrasonic 3D Tag System Based on a Quick Calibration Method. Proceedings of the 2004 IEEE International Conference on Systems, Man and Cybernetics, The Hague, The Netherlands, 10–13 October 2004; pp. 1561–1568.
  6. Duff, P.; Muller, H. Autocalibration Algorithm for Ultrasonic Location Systems. Proceedings of the 7th IEEE International Symposium on Wearable Computers (ISWC ’03), White Plains, NY, USA, 21–23 October 2003; pp. 62–68.
  7. Olson, E.; Leonard, J.J.; Teller, S. Robust range-only beacon localization. IEEE J. Ocean. Eng 2006, 31, 949–958. [Google Scholar]
  8. Ruiz, D.; Ureña, J.; García, J.C.; Hernández, Á.; García, E.; Aparicio, J. Simultaneous Mobile Robot Positioning and LPS Self-Calibration in A Smart Space. Proceedings of the 2010 IEEE International Symposium on Industrial Electronics (ISIE ’10), Bari, Italy, 4–7 July 2010; pp. 2865–2870.
  9. Mautz, R.; Ochieng, W. A robust indoor positioning and auto-localisation algorithm. J. Glob. Position. Syst 2007, 6, 38–46. [Google Scholar]
  10. Guevara, J.; Jiménez, A.R.; Prieto, J.C.; Seco, F. Auto-localization algorithm for local positioning systems. Ad Hoc Netw 2012. in press.. [Google Scholar]
  11. Guevara, J.; Jiménez, A.R.; Morse, A.S.; Fang, J.; Prieto, J.C.; Seco, F. Auto-Localization in Local Positioning Systems: A Closed-Form Range-Only Solution. Proceedings of the 2010 IEEE International Symposium on Industrial Electronics (ISIE ’10), Bari, Italy, 4–7 July 2010; pp. 2834–2840.
  12. Ji, X.; Zha, H. Sensor Positioning in Wireless Ad-Hoc Sensor Networks Using Multidimensional Scaling. Proceedings of the 23rd Annual Joint Conference of the IEEE Computer and Communications Societies (INFOCOM ’04), Hong Kong, China, 7–11 March 2004; 4, pp. 2652–2661.
  13. Kim, E.; Woo, S.; Kim, C.; Kim, K. LaMSM: Localization Algorithm with Merging Segmented Maps for Underwater Sensor Networks. Proceedings of the 2007 Conference on Emerging Direction in Embedded and Ubiquitous Computing (EUC ’07), Taipei, Taiwan, 17–20 December 2007; pp. 445–454.
  14. Hamby, D.M. A review of techniques for parameter sensitivity analysis of environmental models. Environ. Monit. Assess 1994, 32, 135–154. [Google Scholar]
  15. Haensch, T.; Leschiutta, S.; Wallard, A.J. Metrology and Fundamental Constants; IOS Press: Amsterdam, The Netherlands, 2007. [Google Scholar]
  16. Stewart, G.W.; Perturbation, Theory. Least Squares with Errors in the Variables. In Contemporary Mathematics 112: Statistical Analysis of Measurement Error Models and Applications; American Mathematical Society: Washington, DC, USA, 1990; pp. 171–181. [Google Scholar]
  17. Yarlagadda, R.; Ali, I.; Al-Dhahir, N.; Hershey, J. GPS GDOP metric. IEE Proc. Radar Sonar Navig 2000, 147, 259–264. [Google Scholar]
  18. Mekid, S.; Vaja, D. Propagation of uncertainty: Expressions of second and third order uncertainty with third and fourth moments. Measurement 2008, 41, 600–609. [Google Scholar]
  19. Stewart, G.W. Stochastic perturbation theory. SIAM Rev 1990, 32, 579–610. [Google Scholar]
  20. Lilliefors, H.W. On the kolmogorov-smirnov test for normality with mean and variance unknown. J. Am. Stat. Assoc 1967, 62, 399–402. [Google Scholar]
  21. Ji, Y.; Yu, C.; Anderson, B.D.O. Bias-Correction in Localization Algorithms. Proceedings of the the 2009 Global Telecommunications Conference (GLOBECOM ’09), Honolulu, HI, USA, 30 November–4 December 2009; pp. 1–7.
  22. Prieto, J.C.; Jiménez, A.R.; Guevara, J.I.; Ealo, J.L.; Seco, F.A.; Roa, J.O.; Ramos, F.X. Subcentimeter-Accuracy Localization Through Broadband Acoustic Transducers. Proceedings of the IEEE International Symposium on Intelligent Signal Processing (WISP ’07), 1–6 October 2007; pp. 929–934.

Appendix—Calculation of the τ Value

To obtain the value of τ we begin with the first order expansion of F derived from Equation (20):

F = A + A ˜ A + + ( A T   A ) 1 A ˜ T   P

Let ||F||S be the stochastic norm defined as:

F S = 𝔼 ( || F || F 2 )
where 𝔼 represents the expected value. The correspondent squared stochastic norm of F using Equation (26) is:
|| F || S 2 = || A + A ˜ A + + ( A T   A ) 1 A ˜ T   P || S 2
which based on [19] (Section 3.1.3) can be written as:
|| F || S 2 = || A + A ˜ A + || S 2 + || ( A T   A ) 1 A ˜ T   P || S 2

Beginning with the first term of Equation (29) we have:

|| A + A ˜   A + || S 2 = 𝔼   ( || A + A ˜   A + || F 2 ) 𝔼   ( || A + || F 2   || A ˜   A + || F 2 ) || A + || F 2   𝔼   ( || A ˜   A + || F 2 )

The Frobenius norm of the second term in Equation (30) can be simplified as:

𝔼   ( || A ˜   A + || F 2 ) = 𝔼   ( trace   ( ( A ˜   A + ) T   A ˜   A + ) ) = 𝔼   ( trace ( A + T   A ˜ T   A ˜   A + ) ) = trace   ( A + T   𝔼 ( A ˜ T   A ˜ )   A + ) = trace   ( A + T   Q A + )
where the elements quv of matrix Q ∈ ℝ5×5 can be calculated by:
q uv = 𝔼   ( row   ( A ˜ T ) u × column   ( A ˜ ) v ) = 𝔼   ( column   ( A ˜ ) u T × column   ( A ˜ ) v ) = 𝔼   ( a ˜ 1 u × a ˜ 1 v + a ˜ 2 u × a ˜ 2 v + + a ˜ ( m 1 ) u × a ˜ ( m 1 ) v )

Since we assume that all elements of the random matrix à have mean zero, we have that:

𝔼   ( a ˜ ku × a ˜ kv ) = C ( a ˜ ku ,   a ˜ kv )
therefore Equation (32) can be simplified as:
q uv = k = 1 ( m 1 ) C ( a ˜ ku ,   a ˜ kv )

Since matrix Q is symmetric and also positive semidefinite (since it is the product of (ÃTÃ)), we can define Q1/2 as the symmetric square root of Q where:

( Q 1 / 2 ) T   Q 1 / 2 = Q

Replacing Equation (35) in Equation (31) we obtain:

𝔼 ( || A ˜   A + || F 2 ) = trace   ( A + T   ( Q 1 / 2 ) T   Q 1 / 2   A + ) = || Q 1 / 2   A + || F 2

So the first term of the squared stochastic norm of F Equation (29) can be written as:

|| A +   A ˜   A + || S 2 || A + || F 2   || Q 1 / 2   A + || F 2

A similar procedure can be used to simplify the second term of Equation (29), obtaining:

|| ( A T   A ) 1   A ˜ T   P || S 2 || Q 1 / 2 ( A T   A ) 1 || F 2   || P || F 2

Hence:

|| F || S 2 || A + || F 2   || Q 1 / 2   A + || F 2 + || Q 1 / 2   ( A T   A ) 1 || F 2   || P || F 2

Based on [16] (Section 7), from Equation (39) we obtain the expression of τ as:

|| F || S || A + || F max { m 1 , 2 × 5 }   || Q 1 2   A + || 2 τ
since A ∈ ℝ(m−1)×5. From which we obtain Equation (22).

Figure 1. Solvable node subset composed by three beacon nodes and six virtual nodes on a plane.
Figure 1. Solvable node subset composed by three beacon nodes and six virtual nodes on a plane.
Sensors 12 02561f1 1024
Figure 2. LPS configuration used for simulation composed by 3 beacons and a circular path with radius r, height h and m virtual nodes.
Figure 2. LPS configuration used for simulation composed by 3 beacons and a circular path with radius r, height h and m virtual nodes.
Sensors 12 02561f2 1024
Figure 3. Output quadratic mean standard deviation σmean obtained when changing the center of the virtual nodes path: (a) obtained by simulation (1,000 times); (b) estimated offline. In (c) estimated σmean obtained with the offline and online methods compared to simulation.
Figure 3. Output quadratic mean standard deviation σmean obtained when changing the center of the virtual nodes path: (a) obtained by simulation (1,000 times); (b) estimated offline. In (c) estimated σmean obtained with the offline and online methods compared to simulation.
Sensors 12 02561f3 1024
Figure 4. Evaluation of d12 when increasing the input standard deviation from 0 to 0.21 m. The used LPS configuration has a ddop12 of 1.52. In (a) the simulated, estimated offline and online output standard deviation. For the online estimates the mean and the 95% and 5% percentile are shown; In (b) the skewness of the distance d12 is shown.
Figure 4. Evaluation of d12 when increasing the input standard deviation from 0 to 0.21 m. The used LPS configuration has a ddop12 of 1.52. In (a) the simulated, estimated offline and online output standard deviation. For the online estimates the mean and the 95% and 5% percentile are shown; In (b) the skewness of the distance d12 is shown.
Sensors 12 02561f4 1024
Figure 5. Calculated standard deviation errors of distance d12 using a LPS configuration with ddop12 = {7.79, 3.57, 1.52} and a standard deviation input error ranging from 0 to 0.21 m. The correlation between the value τ and the estimated standard deviation error is shown for offline estimation in (a) and for online estimates (the 95% percentile) in (b).
Figure 5. Calculated standard deviation errors of distance d12 using a LPS configuration with ddop12 = {7.79, 3.57, 1.52} and a standard deviation input error ranging from 0 to 0.21 m. The correlation between the value τ and the estimated standard deviation error is shown for offline estimation in (a) and for online estimates (the 95% percentile) in (b).
Sensors 12 02561f5 1024
Figure 6. 3DLocus beacons and virtual nodes configuration used in the test. (a) 3DLocus acoustic localization system; (b) Nodes configuration used in the test.
Figure 6. 3DLocus beacons and virtual nodes configuration used in the test. (a) 3DLocus acoustic localization system; (b) Nodes configuration used in the test.
Sensors 12 02561f6 1024
Figure 7. Experimental evaluation of the estimated distance d12 obtained with different beacons’ subsets. In (a) the estimated online and measured ddop12; In (b) the histogram of the calculated distance d12 obtained using the beacons’ subset {1, 2, 3}.
Figure 7. Experimental evaluation of the estimated distance d12 obtained with different beacons’ subsets. In (a) the estimated online and measured ddop12; In (b) the histogram of the calculated distance d12 obtained using the beacons’ subset {1, 2, 3}.
Sensors 12 02561f7 1024
Figure 8. X-Y coordinates of 3DLocus beacons’ positions obtained with the LAL algorithm, the WLAL algorithm and the inverse positioning method. The solid lines represent the 90% confidence ellipses of the estimated positions. (a) Beacon 2; (b) Beacon 3; (c) Beacon 4; (d) Beacon 5; (e) Beacon 6; (f) Beacon 7.
Figure 8. X-Y coordinates of 3DLocus beacons’ positions obtained with the LAL algorithm, the WLAL algorithm and the inverse positioning method. The solid lines represent the 90% confidence ellipses of the estimated positions. (a) Beacon 2; (b) Beacon 3; (c) Beacon 4; (d) Beacon 5; (e) Beacon 6; (f) Beacon 7.
Sensors 12 02561f8 1024
Figure 9. X-Y coordinates of beacon 5 obtained by simulation. The simulation does not presents the bias observed in the experimental test with the 3Dlocus LPS.
Figure 9. X-Y coordinates of beacon 5 obtained by simulation. The simulation does not presents the bias observed in the experimental test with the 3Dlocus LPS.
Sensors 12 02561f9 1024
Table 1. Standard deviation of the inter-beacon distances dij obtained using the mean σl and the weighted mean σw.
Table 1. Standard deviation of the inter-beacon distances dij obtained using the mean σl and the weighted mean σw.
dijσl[cm]σw[cm]dijσl[cm]σw[cm]dijσl[cm]σw[cm]
d120.1180.041d240.0850.058d370.0600.047
d130.0910.064d250.0940.046d450.0540.041
d140.1150.103d260.1040.102d460.0740.030
d150.0900.036d270.0580.041d470.0480.029
d160.0710.049d340.1160.053d560.0600.047
d170.0530.043d350.1650.144d570.0440.033
d230.0760.055d360.0990.036d670.0460.040

Share and Cite

MDPI and ACS Style

Guevara, J.; Jiménez, A.R.; Prieto, J.C.; Seco, F. Error Estimation for the Linearized Auto-Localization Algorithm. Sensors 2012, 12, 2561-2581. https://doi.org/10.3390/s120302561

AMA Style

Guevara J, Jiménez AR, Prieto JC, Seco F. Error Estimation for the Linearized Auto-Localization Algorithm. Sensors. 2012; 12(3):2561-2581. https://doi.org/10.3390/s120302561

Chicago/Turabian Style

Guevara, Jorge, Antonio R. Jiménez, Jose Carlos Prieto, and Fernando Seco. 2012. "Error Estimation for the Linearized Auto-Localization Algorithm" Sensors 12, no. 3: 2561-2581. https://doi.org/10.3390/s120302561

APA Style

Guevara, J., Jiménez, A. R., Prieto, J. C., & Seco, F. (2012). Error Estimation for the Linearized Auto-Localization Algorithm. Sensors, 12(3), 2561-2581. https://doi.org/10.3390/s120302561

Article Metrics

Back to TopTop