Hindawi Publishing Corporation International Journal of Navigation and Observation Volume 2012, Article ID 603254, 15 pages doi:10.1155/2012/603254

Research Article

Pedestrian Tracking Solution Combining

an Impulse Radio Handset Transmitter with an Ankle-Mounted Inertial Measurement Unit

Joe Youssef, Benoit Denis, Christelle Godin, and Suzanne Lesecq

CEA-Leti Minatec Campus, 17 rue des Martyrs, 38054 Grenoble Cedex 09, France Correspondence should be addressed to Christelle Godin, christeHe.godin@cea.fr Received 16 December 2011; Accepted 17 May 2012 Academic Editor: Amadou Idrissa Bokoye

Copyright © 2012 Joe Youssef et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We address the indoor tracking problem by combining an Impulse Radio-Ultra-Wideband handset with an ankle-mounted Inertial Measurement Unit embedding an accelerometer and a gyroscope. The latter unit makes possible the detection of the stance phases to overcome velocity drifts. Regarding radiolocation, a time-of-arrival estimator adapted to energy-based receivers is applied to mitigate the effects of dense multipath profiles. A novel quality factor associated with this estimator is also provided as a function of the received signal-to-noise ratio, enabling us to identify outliers corresponding to obstructed radio links and to scale the covariance matrix of radiolocation measurements. Finally, both radio and inertial subsystems are loosely-coupled into one single navigation solution relying on a specific extended Kalman filter. In the proposed fusion strategy, processed inertial data control the filter state prediction whereas Combined Time Differences Of Arrival are formed as input observations. These combinations offer low computational complexity as well as a unique filter structure over time, even after removing outliers. Experimental results obtained in a representatively harsh indoor environment emphasize the complementarity of the two technologies and the relevance of the chosen fusion method while operating with low-cost, noncollocated, asynchronous, and heterogeneous sensors.

1. Introduction

For the last past years, new location and tracking (LT) needs have been gradually introduced into a wide variety of applications, such as security, health care, rescue, logistics; or house automation. A growing interest has been more particularly expressed in location-dependent indoor services, which require seamless pedestrian navigation capabilities in harsh environments where satellite-based solutions cannot operate. In this context, alternative technologies are currently under investigation, based on for example, location-enabled wireless networks [1,2].

On their own, most of modern wireless networks can indeed retrieve the positions of mobile radio devices relative to the known position of reference anchors or base stations (BS). The radiolocation functionality simply relies on the measurement of radio metrics, which depend on the distance traveled in the air by transmitted signals. For instance, when a mobile radio device is synchronized with a BS (e.g., using

n-way ranging protocol transactions [3]), range information can be derived from the time of arrival (TOA) of the received signal. Several TOA-based range measurements collected (with respect to fixed BSs with known locations) can hence feed positioning or tracking algorithms to solve out a circular (resp., spherical) location estimation problem in 2D (resp., 3D).

Alternatively, if the surrounding BSs are strictly synchronized (i.e., independently of the clock of a mobile transmitter), the time difference of arrival (TDOA) can be considered, leading to a hyperbolic problem formulation.

Benefiting from unprecedented resolution and synchronization capabilities for the acquisition of such temporal radiolocation measurements, the Impulse Radio-Ultra-Wide band (IR-UWB) technology is today viewed as a credible solution to address LT applications through short-range and low-data-rate (LDR) links [4-9], a fortiori in the context of ad hoc wireless sensor networks (WSN). However, despite the claimed fine properties of IR-UWB signals, obstructed

radio links notoriously introduce additional biases and high dispersion onto measurements, for example, when direct paths are blocked, or when transmitted paths are severely attenuated and shifted in time [3]. Overall, IR-UWB systems claim to provide very fine location and tracking performances in line of sight (LOS) (i.e., typically within submetric precision), but they can hardly guarantee such a precision level under generalized non-line-of-sight (NLOS) conditions (e.g., affecting several links with respect to distinct BSs simultaneously), what is a rather common situation in indoor environments. Hence, one first challenge is to derive robust TOA estimators with practical receiver architectures. One more point would be to evaluate the instantaneous quality of TOA estimates, or even to efficiently remove measurement outliers due to NLOS to assist and enhance the fed LT algorithms.

But it is also well known that radiolocation solutions operating in harsh radio environments could benefit from external means, like assisting inertial navigation systems (INSs) based on inertial measurement units (IMUs). In a navigation scenario, one INS can deliver relevant information out of raw inertial measurements for example, the pedestrian displacement amplitude, velocity, or heading [10-12]. Such information can then be used in addition to radiolocation measurements, as proposed by recent works in the field of hybrid data fusion for example, with GNSS pseudoranges in [13], with IR-UWB TOA in [14, 15], with IR-UWB TDOA and AOA in [16, 17], or even with WiFi received signal strength indicators (RSSI) fingerprints in [18].

Besides strict performance considerations, fusion is all the more relevant, since mobile devices (e.g., personal terminals) are expected in the near future to physically integrate multiple standards and sensors and/or to cooperate with other systems or networks, making heterogeneous modalities naturally available on the user side. Accordingly, several system architectures and configurations involving either collocated or noncollocated radio devices and IMUs can be considered for navigation purposes. As an illustrating example, smart clothes comprising distributed IMUs might form a body area network (BAN) using one first radio access technology and interacting with a portable handset displaying or relaying inertial information to an external access point, to the infrastructure equipping the building, or to another wearable network in the vicinity (potentially through another radio technology). Finally, coupling low-power and low-cost technologies particularly makes sense for the perennial and massive deployment of such systems. This implies the use of adequate radio technologies with energy-efficient transceiver design, reasonably simple inertial units (i.e., with a limited number of embedded sensors), and low-computational complexity for further postprocessing (including data fusion tasks).

In this context, we address herein the pedestrian navigation problem in indoor environments that present dense multipath profiles and magnetic perturbations, by loosely coupling an IR-UWB handset transmitter with a shoe-mounted IMU endowed with a 3-axis accelerometer and a 3-axis gyroscope. A specific extended Kalman filter (EKF) is defined to optimally hybridize the radio and inertial data

and to cope with the specific system architecture constrains (i.e., operating with noncollocated sensors). Adapting and extending previous results from [19], this formulation also combines UWB measurements (obtained through the method proposed in [20]) into new observations defined as combined TDOA (CTDOA). While reducing filter complexity (as a function of the number of available measurements), the proposed solution allows us to remove outlier measurements, without reconfiguring the whole filter structure and without omitting relevant measurement information. The selection of nonoutlier measurements is performed by monitoring the instantaneous quality of TOA estimates, based on a new practical SNR-dependent indicator. As for inertial data, we take benefit from the IMU place to detect stance phases and to reset the foot velocity, hence mitigating in turn the drift that usually affects the estimated INS velocity. The obtained pedestrian heading and average body velocity subsequently feed the fusion filter to control the state prediction phase. Indoor experiments were carried out to validate the proposed fusion approach, as well as to draw intermediary statistical UWB channel parameters useful to TOA estimation.

The paper is structured as follows. In Section 2, we present a selection of techniques and concerns from the recent state of the art in the fields of IR-UWB TOA-based ranging, inertial navigation systems, and the fusion of both modalities. We also try to position the main contributions of this paper in comparison with existing solutions. Then in Section 3, we briefly discuss the specificities of the statistical UWB channel models, and we estimate channel parameters from real indoor measurements. Then a robust TOA estimator adapted to energy detection receivers is recalled. We also show how to practically evaluate the quality of this estimator and to remove outlier measurements, from a filter-oriented perspective. In Section 4, we present our INS, which basically consists of an ankle/foot-mounted magnetometer-free IMU. In particular, we show how to use detected stance phases to remove the drift on INS velocity and infer the average body velocity. In Section 5 we detail and justify further our loosely-coupled fusion strategy, along with the corresponding filter structure. Then we account in Section 6 for experimental results, which were obtained in a typical indoor environment offering mixed operating conditions in terms of signal-to-noise ratio (SNR), radio links obstructions, geometric dilution of precision (GDOP), and magnetic disturbances. Finally, Section 7 concludes the paper.

2. Related State-of-the-Art and Paper Contributions

2.1. Impulse Radio-Ultra-Wideband Time-of-Arrival Estimation. Many solutions have been described in the literature to cope with IR-UWB TOA estimation, including sophisticated algorithms inspired by former high-resolution channel estimation solutions [6] requiring high sampling rates. More recently, various other techniques adapted to the low-complexity LDR context have been proposed and compared (e.g., [7]). A specific focus is usually made on noncoherent

receivers like energy detectors (ED), for which one simple approach consists in comparing the energy collected in consecutive time bins with an appropriate detection threshold. The index of the first bin exceeding this threshold is then associated with the TOA estimate [7, 21]. Unfortunately, within realistic indoor channels, threshold-based ED suffers from overlapping multipath components (MPC) and poor signal-to-noise ratio (SNR) conditions, introducing significant estimation errors and biases. One weakness of these methods is that they do not take benefit from the whole MPC profile (though conveying constructive information), but they only depend on marginal and independent energy terms. Therefore, new TOA estimators have been proposed very recently in [20], assuming realistic path amplitude statistics in compliance with IEEE 802.15.4a recommendations [22], and considering the whole observed energy profile before making a decision on the estimated TOA. These estimators were shown to exhibit low-estimation dispersion around the actual channel leading edge, over a wide range of practical SNRs and channels.

But one more challenge is to allow the real-time prediction of TOA estimation uncertainty or dispersion. As regards to this measurement quality assessment, for most of the TOA estimators proposed in the literature, one could estimate offline the uncertainty based on the received signal power and conditioned on the environment category (e.g., indoor industrial, indoor residential, outdoor, etc.), on the channel configuration (e.g., LOS, NLOS, severe NLOS, etc.), and/or on the actual distance. A relationship could then be established a priori between the error affecting the measured TOA and the SNR through simulations, experimental campaigns, or even theoretical analysis (e.g., [23]). However, this kind of method can hardly benefit from the specificities of the received signal at each instant (e.g., of the current multipath energy profile) under mobility. As an example, under given SNR conditions, the received signal could have either sparse or dense multipath profile, which directly impacts the reliability of the estimated TOA. Moreover, offline characterization is usually mostly intended for ranging performance assessment, but still remains unex-ploited for online tracking purposes.

Hence a new practical method is still required to associate the instantaneous TOA estimation quality with practical estimators (e.g., [20]).

2.2. Inertial Navigation Systems. INSs based on integrated IMUs are more and more used for navigation purposes due to their low cost, low weight, and low consumption. Moreover, as stand-alone autonomous solutions, they can be used indifferently in indoor or outdoor situations, in the lack of surrounding means or infrastructure. They can also be considered for dead reckoning navigation (DRN) (i.e., estimation of the current position from a known starting point) when any absolute positioning system or GNSS is not available, ensuring the navigation service continuity.

A typical IMU consists of a combination of low-cost microelectromechanical sensors (MEMS) for example, acceler-ometers, magnetometers, or gyroscopes. Such units, available as commercial devices now [24], are intended for numerous

applications such as motion capture, unmanned vehicular control, antenna stabilization, video gaming, and pedestrian navigation. Those IMUs provide metrics related to acceleration, orientation, magnetic field, and angular rate. The raw measurements must be processed and analyzed further in order to get relevant navigation information, such as the direction, displacement, speed, and so forth.

In the specific pedestrian navigation context, DRN systems rely on step detection and, for each detected step, on length and heading estimation [25]. Examples of step detection methods are given in [26] or [27]. For step length estimation, two kinds of approaches can be considered. One approach is based on an empirical model linking the step length with some parameters extracted from the sensors measurement during the step [28, 29]. The advantage of this approach is that the models can be adapted to any placement of the IMU. The disadvantage is that a calibration is needed for each pedestrian. A second approach uses the integration of the gyroscope and the accelerometer raw data. Nevertheless, even small errors on the measurements would lead, after integration, to a large error cumulated over time. The drift of the estimated velocity can be reduced with a foot-mounted IMU using zero velocity update (ZUPT) [10], which enables to reset the velocity to zero when the foot is on the floor. Other methods enable to limit the drift of foot-mounted INSs [30]. As no calibration is needed and drifts are limited, this kind ofmethod leads to better performances than empirical methods. However, the IMU has to be placed on the foot in order to benefit from zero velocity. For heading estimation, magnetometers can be used. But they can be subject to local magnetic disturbances induced by pieces of furniture, buried or on-body metallic materials. Those perturbation are particularly present near walls and floors, making hazardous the use of foot-mounted magnetometers in typical indoor environments.

2.3. Hybrid Data Fusion Strategies. Three different fusion strategies can be applied to merge heterogeneous data in a tracking context.

The first one is the so-called decoupled strategy, for example [13]: the mobile node position is estimated separately and independently by each subsystem and then the set of estimated positions (delivered by independent systems) is fused into one final solution. This strategy is preferably applied when both location systems are not subject to drift, or when raw data are not available.

The second strategy is tightly coupled, for example [15]: some raw data are available at each subsystem, which are processed at once to track the location. Hence, it necessitates a global model that relates the heterogeneous measurements of both systems to a correction model accounting for defected measurements (i.e., orientation drift, outlier ranging measurements, magnetic disturbance). The tightly coupled fusion seems to be the best method at first sight since it enables to jointly optimize the output estimate given all the available data. However, these solutions require that the involved subsystems are physically collocated or rigidly connected, what is neither necessarily relevant nor practical for inertial-based pedestrian navigation. Finally, the

synchronization issue between subsystems, which is inherent to any fusion strategy, is all the more critical within these tightly coupled solutions.

The last fusion strategy is loosely coupled, for example [16, 17]: the inputs to the location estimator can be raw and/or preprocessed data issued from sub-systems. This strategy enables to mitigate the previously mentioned drawbacks of the decoupled and tightly coupled fusion strategies. It tolerates that the radio part includes a controlling module to remove outlier measurements and/or to adjust the assumed quality of radio measurements. Finally, this loosely-coupled strategy enables different placements of the involved subsystems.

Besides the coupling strategy itself, the fusion of heterogeneous data also requires the use of specific estimation tools and advanced filtering techniques, such as a particle filter in [14], a combination of a particle filter and an extended Kalman filter (EKF) in [15, 18], or a backward/forward Kalman filter with a recording/smoothing unit in [31]. But another feature of navigation fusion filters concerns the way the IMU data are exploited. In a first approach, the inertial data are used in the prediction phase of the filter [17]. In the second one, they are integrated as observations in the correction phase [16].

In the following, we will provide detailed justifications for retaining a loosely-coupled scheme based on EKF and using inertial data in the prediction phase of the filter, given our system architecture constraints.

2.4. Summary of the Main Paper Contributions. Overall, the main paper contributions are as follows:

(i) application of an energy-based Bayesian TOA estimator under realistic channel parameters;

(ii) proposal of a new quality indicator for TOA estimates, feeding the tracking filter;

(iii) proposal of a specific fusion-oriented filter admitting:

(a) IMU-based pedestrian heading and velocity measurements as control inputs into the filter state prediction;

(b) combined temporal radiolocation measurements as filter observations, reducing computational complexity and mitigating local harmful effects due to EKF linearization;

(iv) evaluation of the previous items within one single unified tracking scheme through representative experiments in a typical indoor environment.

3. TOA Estimation with UWB Energy Detection Receivers

3.1. Received Signal Statistics. The received IR-UWB mul-tipath signal r(t) can be modeled as a sum of weighted

BPF (■)2

Figure 1: Typical block diagram of an energy detection (ED) receiver.

replicas of an energy-normalized pulse p(t) of duration Tp, as follows:

r(t) = ^Ep Y^o-jp^ - Tj - Ttoa) + n(t), j=i

where Ep is the overall channel energy received per transmitted pulse, Tj and aj = fije'°j are, respectively, the arrival time and the normalized complex gain of the jth path (with Xj E[fij ] = 1 without loss of generality), if fij and Oj denote the path amplitude and phase, respectively, and TTOA is the TOA of the first path. The received signal is filtered with an ideal Bandpass filter (BPF) of bandwidth W. Accordingly, a zero-mean additive white Gaussian noise (AWGN) process with double-sided power spectral density N0/2 is also filtered in the band of the transmitted signal into n(t), with the resulting variance a2 = N0W. Finally, we define SNR = Ep/No.

Considering a typical ED receiver [20, 21, 32] (see Figure 1), the overall observation time T is divided into K time slots of length A « Tp. As we assume nonoverlapping replicas ofthe transmitted pulse hereafter, the noise samples taken in different slots are considered as statistically independent. The slots are numbered starting from slot 1 (i.e., for t e [0, A]) up to slot K (i.e., for t e [T - A, T]). The output sample associated with the kth time slot of the ED can be written as

Pk = I I r( t) |2 dt.

J(k-1)A

Let kTOA be the slot index to be estimated, which is associated with the bin that contains TTOA. Then the slots with indexes {kTOA,..., K} correspond to the multipath region. According to [33], under the AWGN assumption the normalized energy samples = ^k/N0 at the integrator output of the ED follow a Chi-square distribution. For k = {1,...,K}, is either a central Chi-square distributed (C^2) random variable (r.v.) with V = 2WA degrees of freedom or a noncentral Chi-square distributed r.v. (NC^2), still with V degrees of freedom, but with a noncentrality parameter Ek = SNRfi that accounts for the energy of the useful signal in slot k. fik is usually assumed as a Nakagami-m distributed r.v. with parameters mk and E[fik] = Ak. Finally, in [34] the average power delay profile (APDP) is modeled by a single exponential decay, with a decay constant e.In [22], the APDP also follows an exponential decay, but with further multipath clustering effects.

To verify these model assumptions and characterize key channel parameters, a channel measurement campaign was carried out in the indoor environment described in

(k - £TOa)A (ns)

Ak-kTOA

---Pp (k - kTOA)

Figure 2: Fitted average power decay and probability of path presence out of measured channel profiles, as two functions of the excess delay.

Section 6.1 for our tracking experiments. Averaging over all the measured channel profiles, an exponential model was fitted to the empirical average power decay, and an empirical function was drawn for the path presence probability, depending on the channel excess delay (or time index, equivalently). To get a better fit of the power decay to our measurements, we preliminarily conditioned received energy upon the presence of a path in each slot so that the slots with no detected path energy were not taken into account in the average. In order to decide whether & contains only noise energy, we simply set a threshold vn equal to 10 dB above the noise power. Finally, the probability of path presence Pp was computed over all the observed channel profiles, as follows:

Pp(k) = Pr(£k > vn). (3)

According to [22], each received profile realization is affected by a distinct shadowing effect. Thus, the measured channel profiles were normalized in power before averaging. Figure 2 represents the exponential power decay model Ak fitted to the measured channel profiles, along with the empirical path presence probability Pp(k), both of which decrease as a function of the excess slot index (i.e., the time index after channel leading edge). For the exponential power decay, the best fit was obtained for e = 40 ns, which is rather close to the mean value e = 39.8 ns found in [34].

Since fa is Nakagami-m distributed with a mean value Ak, then Xk = fik is Gamma distributed with a shape parameter equal to mk and a scale parameter equal to Ak/mk. For the sake of simplicity, we assume that all the received paths have the same parameter mk = m (for all k corresponding to the part of the received signal comprising the useful multipath energy), like in [21]. As shown in Figure 3, for m = 2 the theoretical cumulative density function (CDF) of the normalized path energy satisfactorily fits to the empirical one computed out of our experiments. In the following, the exponential power decay and Nakagami-m parameters, as well as the probability of path presence as a function of the excess time index, are used as prior information to feed our TOA estimator. Similarly, in [21, 33],

- Experimental CDF ---Pr(m, 1/m)> m = 2

----Pr(m, 1/m), m = 1 ...... Pr(m, 1/m), m = 3

Figure 3: Empirical and fitted CDFs of the normalized path energy, for different Nakagami-m path parameter assumptions.

such channel characteristics are assumed available a priori to set optimal detection thresholds.

3.2. TOA Estimation. We consider the minimum mean square error (MMSE) estimator proposed in [20]. Accordingly, the sample index associated with the estimated TOA is defined as

£toa = X kTOA p | c krO\ ), (4)

kTOA=1 V^,..., W

where p( •) is now a conditional probability density function (pdf). Complete analytical developments with their intermediary results can be found in [20]. The TOA estimate is finally obtained as TOA = A fcTOA, where A is the ED sampling period.

3.3. TOA Covariance Estimation. We now intend to provide an indicator reflecting the quality of energy-based TOA estimation. One step beyond, the underlying idea is to feed a tracking filter with further information to dynamically adjust the covariance matrices used in the filter correction step, and/or even to help the detection of outlier TOA-based observations, for instance based on a filter innovation test [35].

For this purpose, we rely here on the instantaneous SNR of the received signal. First, the noise power spectral density N0/2 can be rather straightforwardly estimated in the absence of transmitted signal. Thus we assume that it is available on the receiver side. One ideal solution would consist in jointly estimating the SNR and the TOA. However, for the sake of practicability and simplicity, it is chosen to estimate the SNR first out of all the collected normalized energy samples i^k}k=i,...,K. Unfortunately, there is no close form for the maximum likehood estimator (MLE) of the corresponding

(NCx2) noncentrality parameters {Ek}k=i,...,K. However, it is at least shown in [36] that the following estimator:

Ek = max{& - V,0}

has a lower mean square error (MSE) than the MLE for a single observation and for V > 0.5, which is practically verified in our case. Now, since the TOA is a priori unknown before applying (4), then the mean energy value Ak = E[Ak ] of the multipath components within each time slot k is also unknown. Hence, we use Er = Xk=1 Ek to determine the SNR. Here, two possibilities are available, as follows:

.Erl = X Ek, k=1

Er2 = max|x& - KV,oL

where XK=1 is an NC^2 distributed r.v. with KV degrees of freedom and a noncentrality parameter equal to Er. Both estimators in (6) provide the same estimation result if > V, for all k e {1 : K}, which is more probable at high SNR values. Otherwise, through simulations, we found out that Er2 gives lower MSE than Er1. Moreover, for & > V, for all k e {1 : K}, Er2 is unbiased, and the standard deviation of the relative error (Er2 - Er )/Er is equal to "J4Er + 2V/Er, which decreases when Er increases (i.e., at high SNR).

Let Er = SNR Ar, where Ar is the sum of K Gamma distributed r.v. {Ak}k=1,...,K, each with a probability Pp(k), then Ar is also a r.v. whose characteristic function (CF) (t) is given by

fa (t) = n(pp(k)(1 - J A^t) m + 1 - Pp(k)), (7)

where ±j is the square root of -1, and still assuming mk = m (for all k corresponding to the part of the received signal comprising the useful multipath energy). Suppose now that Er is given, then a nonbiased estimator of the SNR can be built as

SNR = Er

E[Ar ]'

where E[Ar ] is the expected value of Ar, which can be

practically computed as (d^A(t)/dt)| 0 = Sit=1 AkPp(k),

that is to say, as the derivative of (7) with respect to t evaluated at t = 0.

Using the experimental Ak and Pp(k) found in Section 3.1, it turns out that ag^R ~ 0.1SNR practically, which is fairly acceptable.

Finally, SNR is estimated at the real receiver using the approximation Er from (6) instead of Er in (8), as follows:

E[Ar ]'

2 4 6 8 10 12 14 16 18 20 SNR (dB)

- RMSE of the estimated TOA

Figure 4: RMSE of the proposed ED-based TOA estimator, as a function of the SNR value approximated at the receiver side SNR.

At this point, deriving the analytical expression for the root mean square error (RMSE) of the estimated TOA as a function of SNR still remains challenging. Consequently, we proceed by simulations to draw the required relationship. In these simulations, the true TOA is a uniformly distributed r.v. in the time interval [0, T] (i.e., in the entire observation window), and the received signal is simulated according to the statistical parameters given in Section 3.1. Figure 4 shows the RMSE of the estimated TOA as a function of SNR, the latter being computed on the receiver side based on (9). This monotonic evolution can be tabulated as a function of SNR. In an online tracking context, the instantaneous standard deviation aTOAl. of the current TOA estimate with respect to the z'th BS can then be approximated by the previous tabulated RMSE function, assuming unbiased TOA estimate (i.e., <rTOA(- = RMSE ra aTOA.). In the following, irToA- is used in Sections 5 and 6 to determine the instantaneous measurement covariance matrix at each update of a tracking filter (out of the collected energy samples).

3.4. Measurement Outliers Detection and Discarding. In the previous section, we have proposed a dispersion estimator ffTOA, assuming unbiased estimated TOAs. Unfortunately, most of NLOS situations result in biased estimates.

From a tracking filter perspective, a first approach to mitigate harmful NLOS effects consists in increasing the diagonal elements of the measurement covariance matrix. Different solutions can also be used to detect NLOS outliers, for instance based on innovation tests [15, 16]. In this case, the innovation is related to the difference between the current observation and the predicted one, computed under an unbiased hypothesis. However, the complete measurement set would be discarded, even if only one single measurement was biased. Moreover, in some scenarios, severely obstructed NLOS situations might also coincide with situations where the location information is adversely conditioned, for instance due to a poor geometrical dilution of precision (GDOP). In this case, the values in the predicted state covariance matrix are so large (and so are the values in the covariance matrix of predicted measurements) that outliers can hardly be detected through simple innovation

tests. Another method is to preliminarily identify the channel state LOS/NLOS based for example, on collected channel energy [37, 38]. Then, identified NLOS channel conditions are associated with TOA measurement outliers. Whatever the detection method, outliers can be discarded before the location estimation step [15] independently of the latest state estimate. In our context, as ED receivers naturally provide access to the channel energy profile, the last approach is considered.

One indicator of the channel status (and hence indirectly of potential TOA estimation biases) is conditioned on the received signal energy as follows:

Ic = ^ + A2 i (fcToA - M2 p(^V)

kT OA = 1

Practically, this indicator reflects the variance of the estimated TOA conditioned on the energy samples [£k }k=i,...,K of the received signal, averaging over all the possible true kTOA, supposed to be uniformly distributed in the observation window. Since the multipath energy spread of NLOS channels is notoriously larger [22], the hypothesis kTOA = kTOA conditioned on the observed energy samples [£k}k=1,...,K in (10) is likely equiprobable for all the slot indexes. Therefore, the indicator will exhibit larger values for NLOS channels. Note that computing Ia is rather convenient, since all the required quantities are already available for the calculation of fcTOA (see Section 3.2).

Practically, to decide if one TOA measurement is biased (i.e., if it is an outlier) or not, sJI^ is compared with a threshold, which is set at three times above the LOS ffTOA shown in Figure 4. This threshold ensures correctly identifying 97% of the LOS channels under the default centered Gaussian assumption.

4. Inertial Navigation System

As already pointed out in Section 2.2, it is better to place the INS on the foot in order to benefit from ZUPT resetting methods [10] or other update methods [30]. Moreover, we do not consider the use of a magnetometer as magnetic disturbances are too important near the floor. Then we assume a foot/ankle-mounted IMU with a triaxis gyrometer and a triaxis accelerometer only.

First, two main frames are defined:

(i) the body frame (BF) is the frame in which raw data are measured. BF data are referred to their definition frame with (•

(ii) the navigation frame (NF) is the reference frame, which is linked to the earth frame. The NF data are referred to this frame with (• • • )n. Since no magnetometer is used, it is not possible to estimate the orientation of the sensor with respect to the North. Hence the navigation frame includes a rotation around the vertical axis with respect to the North frame defined by the East, North, and Up orientations. This rotation is fixed at the beginning

of the walk, using the projection of the sensor axis on the horizontal plane instead of East and North axis. Vertical axis (orthogonal to the horizontal plane) is given during the stance phase by the measurement of the accelerometer.

Let Rbn be the rotation matrix related to the unitary quaternion qbn, denoting the orientation of the NF in the BF using the relation given in the appendix. Initial orientation with respect to the vertical axis is given, at the beginning of the walk and during the stance phase, by the accelerometer measurement. Considering our definition of the NF, with respect to horizontal axis, initial orientation is set to zero. Then, as the gyroscope measures the angular rate, qbn and then Rbn can be estimated by integrating its raw data. This rotation matrix also transforms the coordinates of a vector in the NF into its coordinates in the BF. Then the raw data issued at the accelerometer can be written as

ab(t) = Rbn(t)an(t) = Rbn(t)(anp(t) - gn), (11)

where gn and anp denote, respectively, the gravity field and the proper acceleration in the NF. The proper acceleration anp corresponds to the derivative of the velocity and to the second derivative of the position. The velocity can hence be estimated by integrating this proper acceleration. In order to limit the drift due to noise or bias integration, it is necessary to detect correctly the stance phase and suppose the velocity is null during this phase, as proposed within the ZUPT method [10].

4.1. Stance Phase Detection. Many step detection methods are proposed in the literature. For instance, in [39] the step detection method is based on the Fourier transform through counting zero-crossing points over a threshold of the accelerometer output. In [26], the pedestrian step pattern is first detected. Then, the beginning and the end of the step can be defined according to this detected pattern. In [27] several other methods are also compared.

In the scenario considered here, the pedestrian is assumed to be walking in a building, moving from one room to another one with some stops and a direction that may change after a few steps. The pedestrian cannot walk very quickly, thus the stance phase of his foot is large enough to be easily detectable. A simple method used to detect the stance phase consists in comparing the acceleration variance to a threshold. Note that even if the ankle trajectory looks reproducible over distinct steps, the pattern ofthe IMU measurements will depend on the sensor orientation, see (11). To overcome this problem, as ||ab(t)|| = ||an(t)|| does not depend on the instantaneous rotation of the IMU, one can simply use the amplitude ||ab (t)|| of the 3D accelerometer measurement [26]. During the stance phase, this amplitude is mostly constant and equal to the amplitude of the gravity vector. We then compute the variance of ||ab (t)|| in a sliding window of temporal length equal to 0.25 s. When the variance value falls under a fixed threshold, a stance phase is detected. During our experiments presented in Section 6, this very simple method has efficiently detected all the steps.

4.2. IMU Proper Acceleration. Let [tf, tf] be the swing phase time interval related to the Kth detected step, and let l be the index of inertial data within this time interval (i.e., {l1 : l2}). Then the acceleration has to be integrated only during this time interval, whereas during the remaining time, the velocity is set to zero. Thus, the velocity of the pedestrian ankle is computed step by step, whereas the orientation of the NF in the BF at time index l is continuously determined with the quaternion qb„ • • • qbP results from the integration of the 3D angular rate measurements tb(l) = [wX® wby!l> wZ®]r, as follows:

(l+1) _ b(l) (l) qbn = q-t 0 qbn,

where 0 denotes the quaternion product and

q-J = [cos(\\wb('-> Ts\\/2) - sin(\\wb('-> Ts\\/2)q^]T,

where Ts = 1/fs is the time sampling interval of the IMU and qt® = tb(l)/||tb® ||. Once the body orientation is determined, the accelerometer measurement terms are simply rotated into the NF as follows:

n(i) _ d№ab(i)

= Rnba ,

where rtnb = Rbn

Then we have

J"* an(t)dt = £ [a"p(t) - gn) dt = £ (anp(t) - gn) dt

= (vn(if) - v?(ti)) - (t$ - ti)gn,

where v; (t) denotes the IMU velocity. Since at t = tf and t = tK the foot is on the floor, then v^ti) = v^tf) = 0, and consequently

an(t)dt.

(tf - t1) Jti

Finally, the proper acceleration is given by anp (t) = (an(t) + gn),

where gn in (17) compensates at the same time the gravity and the velocity bias.

4.3. Inertial Support to Pedestrian Navigation. Our aim is to continuously compute the ankle velocity of the pedestrian together with his heading. Note that the heading is defined as the angle given by the direction of the walk in the horizontal plane of the Navigation Frame. Since during the stance phase the velocity is set to zero, then the velocity is computed separately for each step with

v"(t) = aP (t)dt.

Unlike in [10], the velocity here is continuous between different gait phases after centering the acceleration in (17) during the swing phase of each step. Thus, the computed IMU velocity in (18) starts at zero and ends up at zero.

We consider that the pedestrian walks on a flat floor. The pieces of horizontal information of the ankle velocity and heading are then given as follows:

f? = atanfv| I

\ Vx,l

During a walk with a constant speed, the ankles alternate between stance and swing phases. Hence, the pedestrian waist experiences almost a constant velocity, whereas the ankles experience high velocity variations alternating between maximum and null speed values. As the horizontal waist velocity of the pedestrian does not vary as much as the ankle velocity, then we compute the waist velocity vW using a smoothed version of . In the following, the horizontal velocity vW of the pedestrian and the heading fb are the two processed inertial data simultaneously incorporated into the tracking filter, contrarily to the loosely-coupled fusion strategies in [16,17], where the IMU is used only for heading and step detection. Furthermore, in our proposal, there is no need to have the same pedestrian heading and displacement orientation. For instance, pedestrian side walk and back walk are freely enabled since the real displacement of the foot is estimated. Moreover, there is no need to estimate the step length as in [16, 17] or to calibrate the leg length as in [40].

5. Tracking Problem Statement and Fusion Filter Formulation

In the considered scenario, we remind that inertial sensors are attached on the ankle, whereas the pedestrian holds in his hands an IR-UWB transmitter, and NA known reference receivers are disseminated in the environment. Hence, out of {r;};=1,...,NA, {0;};=1,...,NA, Vw, and f the problem here is to track at least the unknown pedestrian position p = [x y]T, where x and y refer to the Cartesian coordinates. For this purpose, we consider using an Extended Kalman Filter (EKF), which is widely used in mobile tracking applications. After recalling the general EKF formulation in Section 5.1,, we will then justify our overall fusion strategy in Section 5.2. Finally, the state and observation models will be detailed in Sections 5.3 and 5.4.

5.1. Generic Extended Kalman Filter Formulation. We start formulating the tracking problem with the following generic state equation

s(l) = f(s(l-1), u(l-1)

where l is the time index, s(l), u(l) and n(' denote respectively the state vector that contains all the parameters to be

Estimated TOA and related quality

Combined TDOA and scaled observation covariance

Figure 5: Global inputs from building subsystems and fusion filter structure.

estimated, the input control to the dynamic system, and an

1.1.d. process noise sequence. f is the function that relates two consecutive states. As for the observation function h, it binds the state vector to an observation vector r composed of measurements, as follows:

r® = h(s®) + n®, (21)

where n® is an i.i.d. process noise sequence.

The EKF is a popular linearized version of the Kalman filter that can handle nonlinear functions for both estimated states and observations [35]. It enables to estimate the state vector s® and the related state covariance matrix P(l) at time increment l from the previous estimates s(l-1) and P(l-1), the control input u(l-1) obtained at l - 1, and the current observations r(l).

The estimation of s(l) with the EKF is typically split in two steps, namely the state prediction and the state update, which rely, respectively, on the dynamic state equation in (20) and the observation equation in (21).

(a) Prediction:

s<l/l-1) = f (s^-1), u(l-1))

p(l/l-1) = pp(l-1)pr + Q.

(b) Update:

K(l) = P(l/l-1)Hr (HP(l/l-1)Hr + R)-1

¿(0 = ¿(M-1) + K(l) (r(l) - hfs^-1')) (23)

P® = (i - K(°HP(l/l-1),

where F is the Jacobian of f with respect to s computed in ¿(l-1), and H is the Jacobian of h computed in ¿(l/l-1). Q® and R®, which are assumed to be known, are the covariance matrices of n® and n®.

5.2. Overall Fusion Strategy and Filter Structure. In Section 2.3, we presented "decoupled," "tightly coupled," and "loosely coupled" options as three different fusion strategies to merge radio and inertial subsystems into a single navigation solution.

We choose a "loosely coupled" strategy here because it authorizes different placements of the subsystems. For instance, the pedestrian can hold the UWB transmitter in his hand, avoiding many near-floor obstacles, and increasing visibility with respect to anchor nodes, whereas a foot-mounted IMU can take advantage of ZUPT to overcome the INS velocity drift. Moreover, with the "loosely coupled" strategy, each drift-free observation issued from the UWB subsystem can be used to limit the INS drifts.

In our application, we choose to integrate the IMU data in the prediction phase ofthe filter for the following reasons:

(i) it is hard to quantify the errors affecting the inertial data model due to for example, step missing and pedestrian walk behavior (variations of walk speed or unusual pedestrian movements);

(ii) in an alternative approach, the fusion filter update would have to run at the rate of the inertial data instead of the UWB rate (e.g., fs = 200 Hz instead of3-4 Hz in our experiments), while the update part of the filter includes matrix inversions/manipulations and a few matrix multiplications, which is computationally demanding for real-time applications;

(iii) the first proposed approach naturally supports different sampling rates for the prediction and the correction phases. For instance, in our application, the UWB rate is applied to the correction phase, whereas the inertial rate is applied to the prediction phase.

Figure 5 shows the global structure of the proposed EKF fusion filter. The raw inertial data are preprocessed as described in Section 4 to get pedestrian velocity and heading as inputs to the filter prediction phase. In parallel, IR-UWB TOA measurements are delivered at each receiver with their estimated variance, as computed in Section 3.3. In the update part of the filter, specific combinations of such measurements are performed, as introduced in [19], enabling us to remove outliers using the indicator presented in Section 3.4 without changing the filter structure (e.g., the size of involved matrices).

5.3. IMU-Based State Prediction. Generally speaking, at time index l, in the absence of state noise, s® is expressed as a function of the previous state s(l-1), the pedestrian velocity

v(l 1), and the nonbiased or biased compensated heading, p(l-1) as follows:

s(0 = f^-1), v^1), p^1)).

For instance if s = p, then f is fully accounted by

¿0 = x(l-1) +

v(l-1)cos( p(l-1)

yd) = y(l-1) + .

■(^sini rn(l"1)

V(l ^sml p

with fs the inertial refreshment rate.

As already pointed out, the INS on its own suffers from heading drift because ofthe integration ofthe biased angular rate measurement from the gyroscope. Hence, the heading angle fb determined in Section 4.3 is biased by bp. Thus, as one further variable to estimate, we incorporate the latter bias in the modeled EKF state s = [pT bp]T, whereas the control input is composed of the data issued from the INS, namely, u = [vW fb]T (i.e., assuming vW is sufficiently representative for v). Therefore, in our filter the state prediction at time step l is given by

,,(l-i)

£(W-1) = £(l-1) + .

J(l/l-1) = y>(l-1) +

bPl/l-1) = bPl-1),

cos(yjl-1) - S<l-1))

with fs the inertial refreshment rate and hence the Jacobian F of f computed at the predicted state s(l/l-1) is as follows:

vrw phl-1) - bpl/l-1)

,,(l-1),

Jl-1) ftl/l-1)

cos ^ ph - hp

Note that the state covariance prediction also runs at the inertial rate fs. The state covariance matrix Q = Q(l) is supposed to be time invariant. The update part of the EKF is achieved only when new radiolocation measurements are available. Thus of' changes only at the update step, and holds the same value during the prediction step.

5.4. Filter Update with Combined T(D)OA Measurements. At the z'th base station (BS) (i = 1,...,NA), TOA-based measurements r. = c TOA. are performed through ED (see Section 3.2), along with an approximation a. = cOToAi for the corresponding standard deviation a. = c aTOAi. (see Section 3.3). It is also assumed that the NA receivers are

synchronized, but independently of the mobile transmitter, so that all the pseudo-TOA measurements are biased by a common unknown delay.

Filter complexity may be critical with respect to both hardware and software capabilities, especially when many available measurements must be processed in realtime. It is hence worth noting that the most significant part of EKF complexity results from filter gain computations, which directly depend on the size of the observation vector. Hence we consider using the combination-based observations proposed in [19] to reduce the number of observation functions from the number NA of anchors (or measurements) down to the dimension K of the location problem with no precision degradation.

For instance, let W(p) = [w1 w2 ] be an NA X K combination matrix generating K = 2 observation functions used in filter update (out of NA = 4 measurements), as follows:

MP = [hi (p)

Mp)f = W(p)T d(p),

where d(p) = [d1(p),..., d4(p)] is a vector composed of the pseudodistances between the mobile node and the 4 BSs, and it is reminded p = [x y]T. Thus, the same position accuracy as that of the standard TDOA formulation can be preserved using CTDOA, as follows:

vW-1)sin(phl-1) - bPl-1)) (26) where

W(p) = R-1/2[A'^(p) (p)],

AXMp) = AX(p) - (l'LiAX(p)) 1Na,I,

a/ (p) = A;(p) - 1'Na,iA;(p) 1Na,i

AX(p) = R-12 ^di1

ay (p) = R-1/2 ^

with 1NA,1 = r-1/21na,1/Vtr(R-1), if 1na,1 is an Na X 1 vector of ones and Rd the diagonal covariance matrix of TOA-based measurements, whose diagonal is practically composed of of (i = 1,...,Na). Note that similar expressions are available in [19] for the simpler CTOA formulation.

Using such combinations, it is still possible to remove outliers with the indicator Ia., while maintaining a constant filter architecture as a function of time. Within a more classical filter formulation, this would lead to change the size of the filter gain matrix, whereas this size is constant in our proposal and only dependent on K. Practically, so as to discard the ith measurement, the only operation consists in replacing in R-1/2 the ith diagonal element by 0 when computing W(p). Finally, the problem of choosing the best referential BS among the BSs providing nonoutlier measurements does not exist anymore here, unlike in the classical TDOA formulation.

R-UWB Tx handset

mounted 3A3G IMU

Departure Arrival

8 10 m

12 14 16 18

Figure 6: Experimental scenario: pedestrian equipped with an IR-UWB Tx handset and an ankle-mounted 3A3G IMU (a), back-forth real trajectory and layout of the indoor scene, including 4 isochronous IR-UWB Rx (BSs) in room A, with walls separating the 2 Rooms and the corridor.

6. Indoor Experiments

6.1. Experimental Setup. Real-life tracking experiments were carried out in a typical indoor environment at CEA-Leti Minatec premises, as shown on Figure 6. A pedestrian followed a 100m-long round-trip path in two rooms and a corridor, referring to visual markers on the floor. The pedestrian was holding an IR-UWB transmitter in his hand and an IMU was attached to one of his ankles, in compliance with the fusion scenario and system architecture considered so far. We used inertial data from a 3A3G IMU at the sampling rate fs = 200 Hz to estimate the heading and the waist average velocity. The latter was obtained by smoothing the online ankle velocity with a low-pass filter whose impulse response is 1.5s-width rectangular. Without loss of generality with respect to TOA resolution or fusion concepts, the considered IR-UWB transmitter emitted in the band [0.5,1.1] GHz for implementation convenience at the pulse repetition period (PRP) of 1 ^s. UWB antennas were connected to the four synchronous channels of a 6 GHz-bandwidth digital storage oscilloscope (DSO) using cables with N/SMA connectors, serving as four surrounding BSs in Room A. The DSO, enabling signal acquisition at 10 Gsps, integrated a PC for postprocessing. Relying on this setup, after averaging the incoming signal over sequences of 40 successive pulses to increase SNR before energy integration, the received signal was stored for consecutive time intervals of T = 250 ns (i.e., for the duration of the ED observation window). Then the ED-based TOA estimation method described in Section 3.2 was applied to the acquired signals. At the emulated receivers, the observation window was divided into K = 125 resolvable time bins of length A = 2 ns (i.e., approximately the pulse width). The four DSO Rx channels were triggered synchronously but affected by the same unknown delay, which could vary from one acquisition to the next, hence requiring CTDOA as observations in the tracking filter, as

previously mentioned. Finally, the UWB acquisition rate was not constant, and reached 3 to 4 acquisitions per second, which is much lower than the inertial sampling rate.

6.2. Evaluation Procedure and Algorithms Benchmark. Five different estimator settings, depicted as so-called "scenarios" in the following, were tested and compared:

(i) scenario 1: EKF tracking with IR-UWB only, using in the observation vector 2 CTDOA combinations of the 4 available pseudo-TOA measurements, with a timeinvariant covariance matrix R^1) = R^ for all l (i.e., l' being the time step whenever an UWB observation is available, with l a multiple of l');

(ii) scenario 2: idem as scenario 1, except the covariance matrix R^1 ) that is adjusted depending on

{SNRj ure 4);

}/=1,...,4 (using the tabulated function in Fig-

(iii) scenario 3: idem as scenario 2, with additional outliers detection with {I^ }i=1,...,4;

(iv) scenario 4: tracking with a stand-alone INS only;

(v) scenario 5: EKF tracking fusing IR-UWB and INS data as in Section 5, that is, idem as scenario 3, with INS outputs controlling the filter state prediction (with predictions of the IMU frequency and updates whenever UWB measurements are available).

For scenarios 1, 2, and 3 (i.e., UWB only) the state vector of the corresponding EKF just contains the 2D cartesian

.Room A Room B Room A-circle Room B Room A Corridor

Table 1: Average estimated traveled distance per actually traveled meter, for different tracking scenarios.

100 150

Time (s)

■- vw

Figure 7: Pedestrian velocity vW obtained by smoothing the horizontal velocity of the pedestrian ankle vI, along the traveled trajectory.

coordinates and velocities s = [xyvxvy ] , with a classical linear transformation matrix

1 0 dt(l'-1) 0

0 1 0 dt(l'-1) 0 0 1 0 0 0 0 1

where df17'-1' is the time elapsed between samples I' -1 and I' at the UWB subsystem. Scenario 4 being without UWB (i.e., stand-alone INS), only the state prediction step is applied, and the heading bias is systematically set to zero.

6.3. Performance Indicators. To assess and benchmark the performances ofthe proposed tracking schemes, two kinds of metrics are used. First, we consider an error distance between each estimated position p® issued at the filter output at time step I and its nearest orthogonal projection p® onto the actual trajectory

= llp(l) - P(l)l

This instantaneous distance error is finally averaged into Id, over specific portions of the trajectory or geographic areas (i.e., over selected but contiguous sequences of I).

As for the second performance indicator, we consider the difference between the real and estimated traveled distance (over the entire back and forth trajectory), normalized by the real traveled distance. This second relative performance indicator gives an idea about the uncertainty on the overall trajectory length, what could be interesting in several applications besides navigation (e.g., sports analysis, activity monitoring in physical rehab or as dietetics support, etc.).

6.4. Results and Discussion. Relying uniquely on the heading angle fb and on the pedestrian velocity vW (i.e., after smoothing the ankle velocity vi) shown on Figure 7 as a function of time, the tracking performance of a standalone dead-reckoning INS (i.e., scenario 4) is illustrated in Figure 8(a). The drifts affecting both the estimated heading and position clearly justify the use of a side IR-UWB subsystem here.

Figure 8(b) shows the estimated trajectory obtained within scenario 5, where the location drift is now significantly reduced. The dashed part of the trajectory refers to

Scenarios 1 2 3 4 5

D (/1m) 1.8911 2.1088 1.5712 0.8245 0.9466

the first part of the walk (i.e., the one-way portion of the trajectory forth). Noting that BS1 and BS4 are in NLOS configurations in this first portion, thus the mobile could not properly correct the position drift until it gets sufficiently good pseudorange estimates. Consequently, the mobile gets closer to the real trajectory in the middle of the scene, and it even sticks to the real trajectory for the remaining part of the walk (see e.g., the trip back in straight lines).

For each scenario, Figure 9 shows the average location error Id over the entire trajectory and in each room separately. The error is particularly large in Room B for a single IR-UWB radiolocation system (i.e., scenarios 1 to 3). This is due to the combined harmful effects of generalized NLOS links (with respect to the four BSs simultaneously), to body shadowing, and to poor GDOP conditions. Comparing Id for the three first scenarios in Room B specifically, one can notice that using an adaptive observation covariance matrix R(1) (scenario 2) and removing further outlier measurements (scenario 3) clearly help to reduce the error. In more favorable areas, for example in Room A (with at least one LOS links systematically) or in the corridor with light NLOS conditions (i.e., through plasterboard walls), even in the first scenario, the IR-UWB system alone would slightly outperform the INS system alone. As expected, fusing the two subsystems (scenario 5) reduces systematically the overall error, even if the enhancement is far more spectacular in Room B in comparison with both scenarios 3 and 4. These results open the floor to parsimonious fusion schemes, where one could switch from a stand-alone subsystem into the complete fusion-oriented system on demand, depending on the operating conditions, hence saving energy and complexity at the price of slight performance degradations.

Finally, Table 1 shows the average estimated distance D per actually traveled meter. This takes into account the distance between consecutively estimated locations. Due to NLOS situations, a tracking system only based on IR-UWB would tend to overestimate the traveled distance mostly because of occasional but strongly biased TOA-based measurements, which lead to nonstraight and more erratic estimated paths. Omitting measurement outliers, the estimated traveled distance is significantly reduced, but still rather large. The use of INS then enables us to reduce significantly this error.

The previous results illustrate the complementarity of the two subsystems and the potential of the proposed fusion scheme, under the architectural constraints of noncollo-cated and asynchronous sensors. The IR-UWB part of the system tends to correct the heading drift and resolves the growing error of the INS in time, whereas the INS part helps IR-UWB radiolocation in generalized NLOS situations and/or in penalizing mobile locations that would experience bad geometrical configurations with respect to the BSs

8 6 4 2 0

Room B

❖ Corridor

10 12 14 16

Figure 8: Indoor layout of tracking experiments, with 4 BSs in room A, the 2 rooms and the corridor being separated by walls. Actual (red) and estimated (blue) trajectories with an ankle-mounted 3A3G IMU only (a), or a 3A3G IMU loosely coupled with IR-UWB (b).

_ 5 J= 4 IS? 3 2 1 0

All scene

■ Scenario 1

■ Scenario 2 □ Scenario 3

Room A

Room B

■ Scenario 4

■ Scenario 5

Corridor

Figure 9: Average location errors for the whole trajectory and for each room in different tracking scenarios.

(e.g., in room B). Finally coupling both systems enables reliable and robust pedestrian tracking with an almost uniform quality all over the scene.

7. Conclusion

In this paper, we have addressed the problem of pedestrian tracking by coupling an IR-UWB transmitter handset with an ankle-mounted INS device. One motivation to couple these subsystems was to overcome their respective limitations in harmful operating indoor environments, while benefiting from their complementary capabilities.

The raw measurements of each subsystem have been carefully studied, and new preprocessing techniques have been proposed before applying hybrid data fusion techniques. Regarding IR-UWB first, TOAs are estimated at low complexity energy-based receiver following a Bayesian approach. A new practical criterion predicting the standard deviation of estimated TOAs and a methodology to identify outlier measurements have been proposed. Real channel measurements have been carried out and exploited to validate a few statistical multipath parameters. As for INS, in order to limit the drift due to noise integration and

avoid magnetic disturbances, we have considered one ankle-mounted IMU with a 3-axis gyroscope and a 3-axis accelerometer, whose measurements are processed to determine the pedestrian average velocity and biased heading. This ankle-mounted INS (when considered with a UWB Tx handset) obviously imposes further challenging constraints in terms of system architecture, while operating with non-collocated and asynchronous heterogeneous devices, hence impacting in turn the fusion strategy.

To merge the heterogeneous data from both subsystems, we have proposed a specific EKF filter formulation. Combined time difference of arrival [19] of nonoutlier TOA measurements are used as observation inputs to the filter, the observation covariance matrix being dynamically scaled thanks to the new practical TOA estimation quality indicator. The use of CTDOA enables to remove outliers without changing the low complexity structure of the tracking filter. The velocity and heading estimates issued at the IMU are taken into account into the filter during the prediction phase. This option enables us to take benefit from zero velocity information during the stance phase, as well as to operate under different sampling/refreshment rates for the IMU and IR-UWB subsystems.

Experimental results illustrate the complementarity of the two subsystems and the efficiency of the proposed fusion scheme. In particular, the IR-UWB part corrects the heading drift and resolves the growing error of the INS as a function of time, whereas the INS part advantageously assists IR-UWB radiolocation in generalized NLOS situations and/or in penalizing mobile locations that would experience bad geometrical configurations with respect to the BSs. Overall, coupling both systems enables reliable and robust tracking with uniform quality of service over the scene.

Finally, the possibility to apply parsimonious fusion schemes, hence saving energy and complexity at the price of slight performance degradations, has been pointed out and briefly discussed. Accordingly, one could switch from a stand-alone subsystem into the complete fusion-oriented system on demand, depending on the operating conditions (i.e., while experiencing generalized NLOS or generalized LOS, poor or favorable GDOP, etc.).

Appendix

Let qxy = [qxy,o qxy,1 qxy,2 qxyr3] be the quaternion that denotes the orientation of the "X" frame relative to the " Y"

Acknowledgments

This work has been performed in the framework of the WHERE (ICT-217033) and WHERE2 (ICT-248894) projects, which are partly funded by the European Union. The authors would like to thank also the people from Movea (http://www.movea.com/) for their precious technological support as regards to inertial systems. The authors also thank the reviewers for their valuable comments and suggestions.

References

[1] A. H. Sayed, A. Tarighat, and N. Khajehnouri, "Network-based wireless location: challenges faced in developing techniques for accurate wireless location information," IEEE Signal Processing Magazine, vol. 22, no. 4, pp. 24-40, 2005.

[2] K. Yu, I. Sharp, and Y. J. Guo, Ground-Based Wireless Positioning, IEEE Series on Digital & Mobile Communication, John Wiley & Sons, 2009, 2012 DRAFT.

[3] B. Denis, J. B. Pierrot, and C. Abou-Rjeily, "Joint distributed synchronization and positioning in UWB Ad Hoc networks using TOA," IEEE Transactions on Microwave Theory and Techniques, vol. 54, no. 4, part 2, pp. 1896-1910, 2006.

[4] "PART 15.4: Wireless Medium Access Control (MAC) and Physical Layer (PHY) Specifications for Low-Rate Wireless Personal Area Networks (LR-WPANs)," IEEE Std 802.15.4a-2007, Amendment 1, Add Alternate PHYs, 2007.

[5] D. Lachartre, D. Morche, L. Ouvry et al., "A 1.1nJ/bit 802.15.4a-compliant fully integrated UWB transceiver in 0.13,wm CMOS," in Proceedings of the IEEE International SolidState Circuits Conference (ISSCC '09), San Francisco, Calif, USA, February 2009.

[6] J.-Y. Lee and R. A. Scholtz, "Ranging in a dense multipath environment using an UWB radio link," IEEE Journal on Selected Areas in Communications, vol. 20, no. 9, pp. 16771683, 2002.

[7] I. Guvenc, S. Gezici, and Z. Sahinoglu, "Ultra-wideband range estimation: theoretical limits and practical algorithms," in Proceedings of the IEEE International Conference on UltraWideband (ICUWB '08), pp. 93-96, September 2008.

[8] S. Gezici, Z. Tian, G. B. Giannakis et al., "Localization via ultra-wideband radios: a look at positioning aspects of future sensor networks," IEEE Signal Processing Magazine, vol. 22, no. 4, pp. 70-84, 2005.

[9] I. Guvenc, S. Gezici, and Z. Sahinoglu, "Ultra-wideband range estimation: theoretical limits and practical algorithms," in Proceedings of the IEEE International Conference on UltraWideband (ICUWB '08), pp. 93-96, September 2008.

frame, then the associated rotation matrix Rxy is given in the following equation:

[10] E. Foxlin, "Pedestrian tracking with shoe-mounted inertial sensors," IEEE Computer Graphics and Applications, vol. 25, no. 6, pp. 38-46, 2005.

[11] C. W. Tan and S. Park, "Design of accelerometer-based inertial navigation systems," IEEE Transactions on Instrumentation and Measurement, vol. 54, no. 6, pp. 2520-2530, 2005.

[12] L. Fang, P. J. Antsaklis, L. A. Montestruque et al., "Design of a wireless assisted pedestrian dead reckoning system—the NavMote experience," IEEE Transactions on Instrumentation and Measurement, vol. 54, no. 6, pp. 2342-2358, 2005.

[13] F. Caron, E. Duflos, D. Pomorski, and P. Vanheeghe, "GPS/ IMU data fusion using multisensor Kalman filtering: introduction of contextual aspects," Information Fusion, vol. 7, no. 2, pp. 221-230, 2006.

[14] D. B. Jourdan, J. J. Deyst, M. Z. Win, and N. Roy, "Monte Carlo localization in dense multipath environments using UWB ranging," in Proceedings of the IEEE International Conference on Ultra-Wideband (ICU'05), pp. 314-319, September 2005.

[15] J. D. Hol, F. Dijkstra, H. Luinge, and T. B. Schony, "Tightly coupled UWB/IMU pose estimation," in Proceedings of the IEEE International Conference on Ultra-Wideband (ICUWB '09), pp. 688-692, September 2009.

[16] S. Pittet, V. Renaudin, B. Merminod, and M. Kasser, "UWB and MEMS based indoor navigation," Journal of Navigation, vol. 61, no. 3, pp. 369-384, 2008.

[17] V. Renaudin, B. Merminod, and M. Kasser, "Optimal data fusion for pedestrian navigation based on UWB and MEMS," in Proceedings of the IEEE/ION Position, Location and Navigation Symposium (PLANS '08), pp. 341-349, May 2008.

[18] F. Evennou and F. Marx, "Advanced integration of WiFi and inertial navigation systems for indoor mobile positioning," Eurasip Journal on Applied Signal Processing, vol. 2006, Article ID 86706, 2006.

[19] J. Youssef, B. Denis, C. Godin, and S. Lesecq, "Reducing the complexity order of position estimators with combined radiolocation measurements," in Proceedings of the 6th Workshop on Positioning, Navigation and Communication (WPNC '09), pp. 217-222, March 2009.

[20] J. Youssef, B. Denis, C. Godin, and S. Lesecq, "New TOA estimators within energy-based receivers under realistic UWB channel statistics," in Proceedings of the IEEE 71st Vehicular Technology Conference (VTC '10), pp. 1-5, May 2010.

[21] D. Dardari, C. C. Chong, and M. Z. Win, "Threshold-based time-of-arrival estimators in UWB dense multipath channels," IEEE Transactions on Communications, vol. 56, no. 8, pp. 13661378,2008.

[22] A. F. Molisch, D. Cassioli, C. Chong et al., "A comprehensive standardized model for ultrawideband propagation channels," IEEE Transactions on Antennas and Propagation, vol. 54, no. 11, part 1,pp. 3151-3166, 2006.

2 qxy,0 — 1 + 2 qxy,1 2 qxy,1 qxy,2 — 2 qxy,0 qxy,3 2 qxy,1 qxy,3 + 2 qxy,0 qxy,2

2 qxy,1 qxy,2 + 2 qxy,0 qxy,3 2 qxy,0 — 1 + 2 qxy,2 2 qxy,2 qxy,3 — 2 qxy,0 qxy,1 2 qxy,1 qxy,3 — 2 qxy,0 qxy,2 2 qxy ,2 qxy,3 + 2 qxy,0 qxy,1 2 qxy,0 — 1 + 2 qxy,3

[23] D. Dardari, C.-C. Chong, and M. Z. Win, "Improved lower bounds on time-of-arrival estimation error in realistic UWB channels," in Proceedings of the IEEE International Conference on Ultra-Wideband (ICUWB '06), pp. 531-538, September 2006.

[24] http://www.movea.com/.

[25] A. Croci, M. De Agostino, and A. M. Manzino, "A GNSS/INS-based architecture for rescue team monitoring," in Proceedings of the International Conference on Indoor Positioning and Indoor Navigation (IPIN '10), pp. 827-841, Zurich, Switzerland, September 2010.

[26] H.-J. Jang, J. W. Kim, and D.-H. Hwang, "Robust step detection method for pedestrian navigation systems," Electronics Letters, vol. 43, no. 14, pp. 749-751, 2007.

[27] I. Skog, J.-O. Nilsson, and P. Handel, "Evaluation of zero-velocity detectors for foot-mounted inertial navigation systems," in Proceedings of the International Conference on Indoor Positioning and Indoor Navigation (IPIN '10), pp. 959-969, Zurich, Switzerland, September 2010.

[28] S. H. Shin, C. G. Park, J. W. Kim, H. S. Hong, and J. M. Lee, "Adaptive step length estimation algorithm using low-cost MEMS inertial sensors," in Proceedings of the IEEE Sensors Applications Symposium (SAS '07), San Diego, Calif, USA, February 2007.

[29] S. Miyazaki, "Long-term unrestrained measurement of stride length and walking velocity utilizing a piezoelectric gyroscope," IEEE Transactions on Biomedical Engineering, vol. 44, no. 8, pp. 753-759, 1997.

[30] N. Castaneda and S. Lamy-Perbal, "An improved shoe-mounted inertial navigation system," in Proceedings of the International Conference on Indoor Positioning and Indoor Navigation(IPIN '10), pp. 934-944, Zurich, Switzerland, September 2010.

[31] P. Coronel, S. Furrer, W. Schott, and B. Weiss, "Indoorlocation tracking using inertial navigation sensors and radio beacons," Lecture Notes in Computer Science, vol. 4952, pp. 325-340, 2008.

[32] I. Guvenc, Z. Sahinoglu, and P. V. Orlik, "TOA estimation for IR-UWB systems with different transceiver types," IEEE Transactions on Microwave Theory and Techniques, vol. 54, no. 4, pp. 1876-1886, 2006.

[33] I. Guvenc and Z. Sahinoglu, "Threshold-based TOA estimation for impulse radio UWB systems," in Proceedings of the IEEE International Conference on Ultra-Wideband (ICU '05), pp. 420-425, Zurich, Switzerland, September 2005.

[34] D. Cassioli, M. Z. Win, and A. F. Molisch, "The ultra-wide bandwidth indoor channel: from statistical model to simulations," IEEE Journal on Selected Areas in Communications, vol. 20, no. 6, pp. 1247-1257, 2002.

[35] Y. Bar-Shalom and T. E. Fortmann, Tracking and Data Association, Academic Press Professional, 1987.

[36] K. M.-L. Saxena and K. Alam, "Estimation of the non-centrality parameter of a chi squared distribution," The Annals of Statistics, vol. 10, no. 3, pp. 1012-1016, 1982.

[37] S. Gezici, H. Kobayashi, and H. V. Poor, "Non-parametric non-line-of-sight identification," in Proceedings of the IEEE 58th Vehicular Technology Conference (VTC-Fall '03), pp. 2544-2548, Orlando, Fla, USA, October 2003.

[38] A. Rabbachin, I. Oppermann, and B. Denis, "GML ToA estimation based on low complexity UWB energy detection," in Proceedings of the IEEE 17th International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC 06), Helsinki, Finland, September 2006.

[39] Q. Ladetto and B. Merminod, "In step with INS: navigation for the blind, tracking emergency crews," GPS World, vol. 13, no. 10, pp. 30-38, 2002.

[40] V. Renaudin, O. Yalak, P. Tom, and B. Merminod, "Indoor navigation of emergency agents," European Journal of Navigation, vol. 5, no. 3, pp. 36-45, 2007.

Copyright of International Journal of Navigation & Observation is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.