Moving target localization in indoor wireless sensor networks mixed with LOS/NLOS situations

Indoor target localization is an essential and fundamental issue for wireless sensor networks (WSN). However, it is rather difficult for WSN to maintain the localization accuracy in line-of-sight (LOS) and non-line-of-sight (NLOS) mixed environment. NLOS propagation always leads to larger ranging error than LOS does. When the target moves in the rooms and corridors, the signal transmission state will switch frequently between LOS and NLOS. It is a challenging task to deal with this situation because the ranging error characteristics under LOS and NLOS conditions are quite different. In this paper, we propose an interacting multiple model-extended Kalman filter (IMM-EKF) algorithm to improve the localization accuracy for moving target in indoor environment. In the IMM structure, two Kalman filters (KF) are adopted in parallel to accurately smoothen the distance measurement. The proposed algorithm can adapt to the dynamically changing condition between LOS and NLOS due to the two KFs' interaction so that large NLOS ranging errors are further reduced. Once the estimated ranges are obtained, the EKF is employed to estimate the target's location. Empirical measurement results are obtained from typical office environment to verify the effectiveness of the proposed algorithm. Experimental results illustrate that the IMM smoother can efficiently mitigate the NLOS effects on ranging errors and achieve high localization accuracy.


Introduction
In recent years, indoor target localization and tracking technologies are undoubtedly booming in various fields such as robotics and emergency systems. GPS localization is not suitable for indoor environment due to the building's occlusion [1]. Since it owns many advantages such as reliable data and high positioning accuracy, wireless sensor networks (WSN) is widely used for indoor target localization.
According to the relative position between anchor node (AN) and target node (TN), the propagation condition of wireless signal can be divided into such two situations as line-of-sight (LOS) and non-line-of-sight (NLOS). As shown in Figure 1, LOS link represents the direct signal propagation between AN and TN without obstruction. However, in WSN, especially in indoor environment, due to the obstruction such as building structure, furniture, and human activity [2][3][4], the signal propagation will fall into reflection, diffraction, scattering, etc. This indirect signal propagation is defined as NLOS. Whether measured by time-of-arrival (TOA) or received signal strength (RSS), the distances obtained under NLOS link always contain obvious positive bias because the signal travels longer path or time than direct link. Furthermore, the measurement error under NLOS is often larger than that under LOS. With chirp spread spectrum (CSS) nodes [5,6], the different measurement noises under LOS and NLOS condition are shown in Figure 2. Generally, a mixture of LOS and biased NLOS measurements will lead to large error, which is a significant challenge in practice. This paper focuses on how to maintain the accuracy of the mobile localization in indoor environment mixed with LOS/NLOS conditions. Several techniques [7][8][9][10][11][12] have been proposed to mitigate the bad effect of NLOS measurement error. The basic approach is to identify the NLOS condition and then remove associated measurements [7,8]. In this way, the adverse effect of NLOS propagation is eliminated and traditional localization methods can be used to obtain accurate positions. Constrained optimization techniques have also been applied to reduce NLOS errors [9,10]. Furthermore, some researchers even try to explore new methods without identification of LOS/NLOS, but these methods always need either abundant prior knowledge [11] or complex training process [12].
Indoor environment is very complex and should be considered carefully due to the presence of obstructions such as walls, doors, and corridors. When the moving target (MN) moves through rooms and corridors, the propagation conditions change between LOS and NLOS alternatively over a short distance. The range estimation conducted with the NA5TR1 CSS nodes shows that the ranging error characteristics are quite different between LOS and NLOS measurements. This frequent LOS/ NLOS transition will cause a serious measurement error. Therefore, it is significant to cope with the dynamically changing propagation for mobile localization in indoor environment mixed with LOS and NLOS situations. For this case, the transmission channel between the MN and ANs can be treated as a switching mode. We need to determine the LOS/NLOS probabilities of the propagation path between the MN and each AN, respectively. Some proper strategies, such as first-order Markov [13] and IMM [14,15], were proposed to deal with this switching mode system. IMM estimator has been demonstrated as one of the most effective methods to estimate the state of a hybrid dynamic system with several switching modes under uncertain conditions.
In this paper, we focus on how to achieve high accuracy under the dynamic changing LOS/NLOS conditions in indoor environment. We use a Kalman filtering-based interacting multiple model (KF-IMM) smoother to filter the distance measurement in indoor LOS/NLOS situations. Through the LOS and NLOS error analysis on off-line measurements, we get the error parameters in two propagation conditions. After that, we respectively set up two channel error parameters of the Kalman filter model for IMM structure. These two KF filters smooth the distance measurements at the same time, and the model probabilities are calculated adaptively based on the error parameters. Then, the measurement error can be mitigated at the output combination stage. As the filtered distances are obtained, the EKF is employed to achieve the target position. To validate the performance of the proposed method, we conducted the experiment on the fifth floor of the Information Building at Northeastern University.
The main contributions of this paper are listed as follows: (1)We propose a new idea to deal with the frequent transition between NLOS and LOS conditions. Based on the ranging information, the proposed IMM filtering algorithm can adaptively calculate the probabilities of LOS model and NLOS model. Then, the two models interact through the Markov chain automatically. Therefore, the range smoother can response in time when the LOS/NLOS conditions changes. This is significant for range filtering in LOS/NLOS mixed environment. (2)By our proposed algorithm, the filtered results of LOS and NLOS models can be combined automatically by IMM structure. According to the different error characteristics of LOS and NLOS, we adopt two parallel KFs. The range filtering results of the two KFs are combined based on the model probabilities. In this way, it is not only feasible for severe NLOS condition, but also viable to various degrees of NLOS. algorithm. Experimental results show that positive NLOS error is significantly mitigated and moving target localization can be accurately achieved.
The remainder of this paper is organized as follows: Section 2 gives a brief overview of the localization techniques in NLOS environments. The system model is described in Section 3. In Section 4, the IMM smoother and EKF-based localization algorithm are proposed. The performance of our algorithm in practical indoor environment is presented in Section 5. Finally, the conclusion is drawn in Section 6.

Related works
Traditional localization algorithm based on WSN mainly includes trilateration localization, triangulation localization, and least squares. These algorithms directly measure the ranges or angles, which are used for the localization. Under the LOS environment, these algorithms can achieve high accuracy. But for NLOS environment, they cannot obtain desired accuracy because of the significant NLOS measurement error [16,17].
Nowadays, researches toward target localization in NLOS environments mainly focus on LOS/NLOS identification and NLOS error mitigation. The basic idea of NLOS identification algorithm can be divided into two types, including measurement-based identification [18,19] and the channel statistics-based identification [20][21][22]. An angle of arrival (AOA)-based NLOS identification for multiple-antenna system is proposed in [18]. On the basis of the phase difference across two antenna elements, the hypothesis test on K-factor is adopted to identify the NLOS state. However, the limitation lies in that AOA-based measurement needs additional hardware support. Reference [19] presents a TOA-based NLOS identification method with prior knowledge such as the standard deviation of measurement noise. RSS-based identification algorithms are also widely used, but the ideal accuracy can only be achieved in ultra wideband condition [20]. In [21], the statistics of radio propagation channel metrics along with binary hypothesis testing are applied to identify the NLOS states. In [22], the RSS and time characteristic are utilized to estimate the NLOS situations, but the analysis process is relatively complicated.
Most of the NLOS mitigation technologies are dependent upon the NLOS identification. Residual test algorithm, presented in [23], can identify LOS ANs to avoid NLOS impact by discarding the NLOS measurement. However, it might be unfeasible when the number of available range measurements is limited. In [11] and [24], a non-parametric approach is proposed for NLOS identification based on support vector machines, which can attain good accuracy at the cost of high computation. In [25], a linear programming (LP) approach is adopted for NLOS mitigation. Some researchers propose soft-decision algorithms without NLOS identification and discarding [26][27][28]. In [26,27], Chen uses a residual weighting (RWgh)-based algorithm to alleviate the NLOS errors. On the basis of Chen's research, Hammes and Zoubir [28] suggest a data fusion algorithm to achieve higher localization accuracy. However, RWgh algorithms cannot perform well when the number of LOS nodes is deficient. But, the computational complexity will grow rapidly when the number of ANs increases. In [29], vehicle velocity and heading direction measurements are exploited to make constrained optimization for NLOS error mitigation. Other researchers propose some model matching and database linked methods, such as fingerprinting algorithm, which are always time-consuming [12]. Some researchers begin to explore new methods without NLOS identification. Typical algorithms include least square multi-lateration [30], optimized residual weighting [31], support vector machine classifier [11], and machine learning [32]. These methods often require ideal synchronization or massive experimental data.
As for frequent LOS/NLOS transition, researchers have made good progress with interacting multiple mode. Yang et al. [33] propose a location estimation scheme using fuzzy-based IMM smoother for mobile location. IMM is employed as a switch between the LOS and NLOS states, which are considered to be a Markov process with two interactive modes. Associated with data fusion, this scheme can efficiently mitigate the NLOS effects on the measurement. In [15], interacting multiple model is proposed for moving target location estimation in cellular network. In the rough wireless environment such as urban area, the proposed method can efficiently mitigate the NLOS effects on the ranging error. However, the urban structure is quite different from indoor environment. The NLOS situation in urban area is caused by obstruction due to large buildings, while the indoor NLOS always considers the signal attenuated by walls or doors. Taking this into consideration, the measurement model should be different between these two environments.

System model
The WSN-based indoor moving target localization system contains several ANs and a moving node (MN). MN communicates with ANs to acquire distance measurements, with which the position of MN at a particular time step can be estimated.
At the time step t, the distances obtained through the communication between MN and M ANs by CSS system are used to build the distance vector: where d mes m t ð Þ is the range measurement between the mth AN and the MN. Our experiment is operated in the office building as shown in Figure 3. The channel conditions between MN and ANs change dynamically when the target moves. The statistics of channel conditions between the moving MN and ANs are shown in Figure 4. The distances between AN and MN are directly acquired by NA5TR1 CSS node ranging, which is defined in the standard IEEE 802.15.4a [5,6]. CSS ranging is based on the time of flight (TOF) of radio frequency signals. In the CSS system, the TOF is calculated by symmetric double-sided two-way ranging (SDS-TWR) [34]. In this paper, we assume that the error parameters, which are calculated with offline measurements, of LOS and NLOS follow the normal distribution.
At time step t, the range state vector between the mth AN and the MN is defined as where d m (t) represents the distances between the mth AN and the MN, _ d m t ð Þ is the velocity of the MN according to mth AN, and M is the number of the ANs. The state-space model can be expressed as linear dynamic equations Ti represents the sampling time interval, ω d (t − 1) denotes the process noise, which is assumed as independently and identically distributed (i.i.d.) Gaussian noise, with zero mean and covariance matrix Q d .
The MN's state vector at each time step t is defined as a two-dimensional position and velocity vector of the MN.
where (x(t), y(t)) denotes the coordinate of the MN at time step t, while _ x t ð Þ; _ y t ð Þ ð Þ represents the velocity of the MN in x and y directions. The MN's state updates over time according to the random force model. The range measurement matrix has been defined in Equation (1). Then, let where h(X(t)) = [h 1 (X(t)), h 2 (X(t)), ⋯, h M (X(t))] T describes the Euclidean distance between MN and the ANs; the subelements are defined as  error μ d and variance σ 2 d . Here, μ d and σ 2 d are assumed to be i.i.d. Gaussian with different parameters according to the LOS and NLOS states

Target localization based on IMM algorithm
Since the indoor construction is complicated and the communication signals are obstructed frequently, the NLOS propagation condition will be caused between MN and ANs. However, the measurement errors under LOS and NLOS conditions are quite different. Since the uniform measurement filtering method cannot fit both LOS and NLOS conditions, it is necessary to introduce a mixture of LOS and biased NLOS range filter. Therefore, we adopt the IMM algorithm to filter the measurements under the LOS/NLOS conditions. Then, EKF is employed to estimate the MN's location.

Measurement filtering based on IMM algorithm
Single model is insufficient to capture the measurements under LOS/NLOS mixed environment, because the measurement errors are quite different between LOS and NLOS conditions. The IMM structure employs two parallel Kalman filters to smoothen the measurement errors in LOS and NLOS separately. Then, the state estimates of both filters are combined automatically according to the model probabilities. As the MN moves, due to the construction and furniture, the communication conditions between MN and ANs frequently alternate between LOS and NLOS. The LOS/NLOS transition channel between AN and MN is considered as a switching-mode system, which can be seen as a firstorder time-homogeneous Markov chain [15] depicted in Figure 5. Markov chain consists of two states. We let i and j (i, j = 1, 2) denote the model variable, where i, j = 1 represents the LOS model, and i, j = 2 represents the NLOS model. The transition probabilities p ij denote the conditional probability of states transition from status i at time step t − 1 to status j at time step t. Let D i (t) denote the input of the state model i at the time step t, then where represent the normalized transition probabilities, and p 11 + p 12 = 1.
The framework of the IMM algorithm is illustrated in Figure 6. The IMM algorithm consists of four steps: input interaction, model filtering, model probability update, and output combination.

Input interaction
Since the IMM algorithm is recursive, the inputs of both KFs in the IMM structure at time step t come from the two KFs' filtering results at time step t − 1. The mixing probability μ m,i|j (t|t − 1) representing the mixing probability of the ith KF filtering result inputs to jth KF filter can be calculated as where μ m,i (t − 1) represents the ith model probability of the mth AN at time step t − 1, p ij denotes the transition probability defined in Markov chain, and c m;j is a normalized variable, which is defined as With the input mixing probabilities, the state estimation for mode inputD m;0 j (t|t − 1) and a priori covariance matrix for the jth mode-matched filter of mth AN can be computed aŝ where P m,i (t − 1|t − 1) denotes ith model's posterior probability at time step t − 1.

Model filtering
As mentioned above, we adopt two parallel KFs in this step. Based on the a priori knowledge that the measurement noises are quite different between LOS and NLOS conditions, one KF is designed for LOS measurement filtering and another is for NLOS. The KF can be described in two steps as Prediction: Update:

Model probability update
When the estimated states are obtained, the model likelihoods and probabilities are needed to be calculated to prepare for the output combination.
where Λ m,j (t) denotes the Gaussian density function of residual error V m, j (t), with zero mean and covariance S m,j (t).

Output combination
According to the model probability, the estimated distance can be updated as The value of μ d,i refers to Equation (9), i = 1 means LOS condition, i = 2 means NLOS condition.
The final distance estimation can be calculated aŝ d m t ð Þ ¼ GD m t tÞ: j ð ð 28Þ

Mobile localization
We use EKF to deal with non-linear dynamic problems of the MN localization. The basic idea of EKF is to convert the non-linear system according to the one-order Taylor function expansion into the linear system. The process of EKF-based MN localization is described as Prediction: Update: whereD Arr t ð Þ ¼d 1 t ð Þ;d 2 t ð Þ; ⋯;d M t ð Þ h i T denotes the estimated measurement vector obtained from the IMM filtering. Then, the Jacobian matrix H can be defined as The measured distances from LOS or NLOS have been mitigated by IMM. Furthermore, the EKF has smoothed the estimated error at the localization stage. Combining the IMM and EKF, we can get a more accurate estimation.

Experimental result and analysis
To evaluate the performance of our proposed algorithm, practical experiments are performed on the fifth floor of the Information Building of Northeastern University. We adopt the NA5TR1 CSS nodes in the field experiments.
The experimental environment and deployment of the ANs are described in   [15] and the single model EKF localization algorithm. The IMM measurement filtering results are shown in Figure 7.
As the results of all ANs are similar, we just pick up AN1 as an example. Figure 7a shows the model probabilities of AN1. Contrast to the transition conditions of AN1, the probabilities of most of the time steps are reasonable. When AN1 is under the NLOS conditions, the NLOS model probability is larger than that of the LOS model, and when the condition turns into LOS, the probabilities will switch back timely. The detailed distance estimation errors along the path are shown in Figure 7b. It can be seen that at the time step t, the ranging errors are large due to the NLOS condition. However, with the IMM filtering, the errors are obviously decreased. The error statistics of all the ANs are shown in Table 1. The result demonstrates that the proposed IMM algorithm can effectively adapt the changes of the propagation conditions and significantly decrease the NLOS ranging errors.
With the filtered distance from the IMM algorithm, the MN's locations are calculated by the EKF. The IMM-EKF tracking result is compared with those of the IMM-only algorithm and EKF-only algorithm. Most of the related parameters in IMM are identical to IMM-EKF. For comparison, the EKF algorithm adopts the measurement error distribution as Gaussian distribution N(1,4). All visualized tracking results are shown in Figures 8, 9, and 10. The positioning errors are indicated by the lines connecting the estimated location (*) and real location (dots). The cumulative distribution function (CDF) of positioning error is shown in Figure 11.
As a valid algorithm for switching mode, the IMM algorithm has only been used successfully for LOS and NLOS environments in urban area. However, for wireless signal transmission, there is a large difference between outdoor and indoor environments. To compare the performance and effects of the different methods, the algorithm proposed in [15] is modified at some extent to be used for comparison in indoor LOS/NLOS environment. From the simulation results shown in Figures 8, 9, 10, and 11, the proposed IMM-EKF smoother can more accurately track the true range distances in indoor LOS/NLOS mixed situations. The main reason lies in that the mean noise error is introduced into the weighted fusion estimation at the stage of ranging filter, as shown in Equation (26). At the same time, the EKF algorithm goes a step further to depress the localization error. From the CDF of the positioning error in Figure 11, we can intuitively see the error distribution of EKF, IMM, and IMM-EKF algorithms. When the EKF and IMM are employed in the LOS/NLOS mixed environment, the localization procedure can achieve the accuracy of about 2.3 m. However, the IMM-EKF can achieve higher accuracy with 100% location error below 1.7 m, which is the best. Above all, the single model cannot fit the LOS/NLOS mixed environment very well.
The detailed statistical characters about the maximum, average, standard deviation, and 90% of the positioning error are summarized in Table 2. With the IMM-EKF algorithm, we can expect a localization performance of 90% error below 1.29 m. For general applications such as personal tracking and robot navigation, it is an acceptable result.
Experiments with different Markov chain probabilities are carried out to investigate the relationship between the Markov chain probabilities and localization performance. We assume that p 11 + p 12 = 1 and p 12 = p 21 . Therefore, once one of the transition probabilities is determined, the others are also fixed. So, we select p 11 as the variable to conduct the evaluation. Figure 12 shows the variation of the mean positioning errors due to the alteration of Markov chain probability. When the p 11 = 0.95, the IMM-EKF achieves the best localization accuracy.

Conclusions
This paper proposes a real-time localization algorithm, based on IMM-EKF, for moving target localization under LOS/NLOS mixed indoor environment. When the system obtains the measured distance vector between the MN and the ANs, the IMM is employed to mitigate the ranging errors either for the LOS or NLOS conditions.
With the filtered ranging result, the EKF algorithm is adopted to estimate the position of MN. In the range filtering stage, the measurement errors, especially the NLOS errors, are obviously weakened. By prediction and update, the EKF-based localization stage makes further improvement on the accuracy. Experimental result shows that the algorithm proposed in this paper achieves high accuracy even in complex environment where LOS and NLOS channel conditions switch randomly and frequently. This research is accomplished on the basis of a priori knowledge of the ranging error distribution for LOS and NLOS conditions. For further works, we will focus on the research to update the error parameters during the localization process.