Next Article in Journal
A Review of Infrared Thermography for Delamination Detection on Infrastructures and Buildings
Previous Article in Journal
Non-Destructive Detection of Damaged Strawberries after Impact Based on Analyzing Volatile Organic Compounds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Constrained Kalman Filter for Wi-Fi-Based Indoor Localization with Flexible Space Organization

by
Vincent Sircoulomb
* and
Houcine Chafouk
IRSEEM, ESIGELEC, UNIROUEN, Normandie University, 76000 Rouen, France
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(2), 428; https://doi.org/10.3390/s22020428
Submission received: 22 November 2021 / Revised: 17 December 2021 / Accepted: 29 December 2021 / Published: 6 January 2022
(This article belongs to the Section Navigation and Positioning)

Abstract

:
This paper presents a constrained Kalman filter for Wi-Fi-based indoor localization. The contribution of this work is to introduce constraints on the object speed and to provide a numerically optimized form for fast computation. The proposed approach is suitable to flexible space organization, as in warehouses, and when objects can be spun around, for example barcode readers in a hand. We experimented with the proposed technique using a robot and three devices, on five different journeys, in a 6000 m2 warehouse equipped with six Wi-Fi access points. The results highlight that the proposed approach provides a 19% improvement in localization accuracy.

1. Introduction

Nowadays, a wide-scale proliferation of wireless devices can be observed. This phenomenon is due to multiple factors, such as more and more efficient, small and low energy consumption electronic components and widespread wireless communication technologies. Among the services allowed by this technological progress, LBS (Localization Based Service) is facing a continuously increasing interest, notably due to dedicated smartphone applications. According to the statistics, 80% of LBS are for indoor purposes [1], for example for guidance of persons with mobility problems, emergency evacuation or goods localization in warehouses.
In the past decade, plenty of research has been conducted about LBS [2,3], in particular using wireless techniques [4,5,6] such as ultra-wideband, Bluetooth, or Wi-Fi. Wi-Fi has the advantage to be cost-advantageous. Nowadays, it is also very widespread, which avoids the installation of new antennas or additional transponders as for RFID (Radio Frequency IDentification) for example. Position information can be obtained through the wireless signal RSSI (Received Signal Strength Indication) [7] or the waves’ AoA (Angle of Arrival) [8], knowing the position of the emitters. However, such information is generally insufficient to provide accurate localization, especially using Wi-Fi. Thus, wireless localization is traditionally coupled to additional techniques such as fingerprinting [9]; the use of space grids [10,11]; or of additional sensors like smartphone ones, namely accelerometers, gyroscopes, or magnetometers [12]. This last technique is usually named inertial navigation and includes a variant called PDR (Pedestrian Dead Reckoning).
To fuse information coming from Wi-Fi and from one or several additional technique(s), Kalman filtering is the most popular approach [13]. Many improvements to the baseline Kalman Filter (KF) have been used, such as the extended (EKF) or unscented KF (UKF) to deal with the orientation nonlinearity, multi-stage, and adaptive or robust versions. In [14], landmarks were used as inputs of a limited complexity KF to improve the accuracy. In [15], two EKF were used to fuse the measurements and to update a fingerprinting database. In [16], an adaptive UKF was used to improve the PDR with a map-matching technique. In [17], particle swarm optimization was used to improve the correction step accuracy. In [18], several movement scenarios were considered with an adaptive and robust filter. In [19], a two-stage version was introduced to improve the robustness by separately dealing with attitude determination and tracking.
Another improvement in KF is state constraints. State constraints are all additional equations that do not fit into the structure of a KF [20]. For example, when estimating a quaternion, it is known that its Euclidian norm is equal to one. Taking into account such additional knowledge in the estimation process contributes to improving the localization and tracking accuracy.

1.1. Related Work

Extensive research has been conducted to incorporate state constraint equations into the structure of a state estimator—through model reparameterization or model reduction, pseudo-observation methods, gain modification techniques, projection approaches, etc. If both the system and constraints are linear, all these approaches result in the same state estimate [21]. For LBS purposes, constraints about the heading were taken into account in [1,22], allowing for a reduction of the state space-tracking model to a four-order linear one. In [23], the heading was constrained according to the type of environment: narrow passages like corridors versus wide scenes with an azimuth wheel. In [24], position constraints were considered in an UKF when at most four anchors contributed to the localization. In [25], nonlinear inequality constraints were considered for the measurements, in particular RSSI and AoA. However, it should be noted that no work has been done to consider speed constraints in KF-based LBS.

1.2. Contribution

In this paper, the contribution is to develop a Kalman filter for Wi-Fi-based indoor LBS with constraints on the moving object speed. This contribution arises from the kind of indoor environment under consideration in this work: a warehouse with flexible topology and space organization, according to the products stored. The nature of the products can badly affect Wi-Fi coverage, impacting localization. Moreover, small IMU chips with a limited cost usually have larger drifts and errors, which do not contribute to improving the positioning. State constraints are thus considered to compensate for this accuracy decrease. The flexible topology implies that the paths can evolve with time, making position constraints unsuitable. Moreover, the objects to localize can be spun around anywhere, preventing the consideration of heading constraints. We thus introduce speed constraints according to the limited knowledge we have about the flexible organization. A corresponding Kalman LBS algorithm is developed in this paper and the constraint stage is detailed to provide fast computation with a minimal amount of scalar additions and multiplications.

1.3. Outline

This paper is organized as follows. The problem to solve is first stated in Section 2. Then, the proposed approach to deal with the problem under consideration is depicted in Section 3. The results obtained on several practical case studies in a warehouse are given in terms of plots and figures in Section 4. These results are then analyzed and discussed in Section 5. Finally, concluding remarks and the outlook of this work are given in Section 6.

2. Problem Statement

Let us consider a device moving into an indoor environment. The movement is supposed to be horizontal. The device is equipped with an IMU (Inertial Measurement Unit). The objective is to track the device, that is to say to estimate its position in real time. The device movement can be modeled by Equation (1).
r x , k r y , k V x , k V y , k θ k = 1 0 T s 0 0 0 1 0 T s 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 r x , k 1 r y , k 1 V x , k 1 V y , k 1 θ k 1 + 0 0 0 0 0 0 cos θ k 1 sin θ k 1 0 sin θ k 1 cos θ k 1 0 0 0 1 Δ V x , k Δ V y , k Δ θ k ,
where k * is the discrete time; T s is the sampling time; r x , k ; r y , k and V x , k ; V y , k denote the 2D position and velocity, respectively, at time k ; and where θ k is the device heading (see Figure 1). Δ V x , k , Δ V y , k , and Δ θ k are the IMU measurements, that is to say the velocity increments between discrete times k 1 and k along the G , x m and G , y m axes and the angle increment around the G ,   z axis. Writing down on the first hand:
x = r x r y V x V y θ ,     u = Δ V x Δ V y Δ θ ,     F = 1 0 T s 0 0 0 1 0 T s 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 ,     G x = 0 0 0 0 0 0 cos θ sin θ 0 sin θ cos θ 0 0 0 1 ,
and considering process noises in the second-hand, Equation (1) can be rewritten as follows:
x k = F x k 1 + G x k 1 u k + s k + w k .
where s k is the additive noise affecting the IMU measurements and w k is introduced to take into account the model approximations (discretization, neglecting the earth rotation effect, etc.).
To compensate for the mid- and long-term effects of the stochastic processes s k and w k , the device position is also measured according to Equation (4):
RSSI i = 10 n log 10 d i + A i ,
where RSSI i is the Received Signal Strength Indication of the Wi-Fi signals exchanged with the access point i   1 i N , n is a factor depending on the environment ( n = 2 for outdoor and 4 for indoor situations), d i is the distance between the device and the access point, and A i is a calibration coefficient. Introducing the following:
y k = RSSI 1 , k RSSI N , k ,     γ k = γ 1 , k 0 0 0 0 0 0 γ N , k     with :   γ i , k = 1   if   RSSI i   is   detected   at   time   k 0   otherwise
Equation (4) can be rewritten in the following compact form:
y k = γ k h x k + v k ,
where v k is the measurement noise. The problem to solve is to estimate as accurately as possible x k for any time k using the set of measurements y 1 , , y k and their existence γ 1 , , γ k .
Using Equations (3) and (6), the estimation can recursively be done using an EKF, because function h is easy to linearize and does not contain discontinuities. The EKF Equations (7)–(8) for the prediction step and Equations (9)–(11) for the Joseph’s form correction step are given below:
x ^ k = F x ^ k 1 + G x ^ k 1 u k ,
Σ k = F Σ k 1 F T + G x ^ k 1 P G T x ^ k 1 + Q ,
K k = γ k Σ k H k T H k Σ k H k T + R 1 ,
x ^ k + = x ^ k + K k y k h x ^ k
Σ k + = I K k H k Σ k I K k H k T + K k R K k T
where I is the identity matrix of appropriate dimension; the superscript T denotes the transposition operation; and P , Q ,   R , Σ k , and Σ k + are the following covariances:
P = E s k s k T ,           Q = E w k w k T ,           R = E v k v k T , Σ k = E x k x ^ k x k x ^ k T ,       Σ k + = E x k x ^ k + x k x ^ k + T .
The filter is initialized with x ^ 0 and Σ 0 . It should be noted that at a given time k , if no RSSI is detected, then γ k is null and so x ^ k + = x ^ k and Σ k + = Σ k . As output of the EKF, the estimated state and error covariance are as follows: x ^ k x ^ k + and Σ k Σ k + . The matrix H k is obtained by linearizing h , as follows:
H k = h x x = x ^ k = 10 n log 2 10 x ^ k 1 X b 1 d 1 x ^ k 2 Y b 1 d 1 0 0 0 x ^ k 1 X b N d N x ^ k 2 Y b N d N 0 0 0 ,
where X b i ; Y b i are the coordinates of the access point i   and where d i = x ^ k 1 X b i 2 + x ^ k 2 Y b i 2 is the distance between the object predicted position x ^ k 1 ; x ^ k 2 and the access point i .
Localizations using Wi-Fi RSSI, IMU, and EKF usually give pretty good results, with the accuracy depending on the number of access points. However, in warehouses, the Wi-Fi coverage is often depreciated because of metallic shelving and harsh stored products such as pallets of water bottles, conserves, or chemical products. Moreover, the topology of warehouses can frequently change because of permanently changing stocks. Thus, additional state constraints are introduced in the baseline EKF to improve the localization accuracy.

3. Proposed Approach

Let us consider an indoor environment such as a warehouse. It includes several rectangular storage areas containing shelves, as illustrated in Figure 2. Each storage area can have its own organization changed, i.e., the shelves’ position and quantity can be modified according to the products stored. Consequently, the corridors between the shelves are not fixed. Outside the storage areas, the movements are free. For circulation convenience, all shelves are parallel to the O , y axis.
Because the shelves’ localization into an area are subject to possible modifications, position constraints are not suitable. Moreover, the objects to localize can be spun around, for example a barcode reader in a hand when scanning some products put on the shelves. Thus, no additional constraint about the heading could be taken into account. However, as the displacements between the shelves are along the O , y axis due to the area arrangement, the speed along the O , x axis is close to zero. This constraint is taken into account according to a smooth additional measurement, as in [26]:
D x k + ζ k = 0     with :   D = 0   0   1   0   0 .
Constraint (14) is equivalent to: V x , k = ζ k , where ζ k is a zero-mean white random process with a standard deviation σ . It imposes the speed along the O , x axis to be close to zero, but not strictly, because the object can nevertheless laterally navigate between the left and right shelves. Equation (14) must be considered only if the object is inside a storage area. Writing down X i min , X i max , Y i min , Y i max , the limits of area i (see Figure 2), and considering the uncertainty of the position estimate x ^ k + 1 ; x ^ k + 2 , the condition is as follows:
X i min + 3 Σ k + 1 , 1 x ^ k + 1 X i max 3 Σ k + 1 , 1 Y i min + 3 Σ k + 2 , 2 x ^ k + 2 Y i max 3 Σ k + 2 , 2 ,
where Σ k + i , j is the coefficient at the i th row and j th column of Σ k + . Considering the noise distribution is Gaussian, the probability of wrong detection of a storage area is less than 0.1%. When condition (15) is satisfied for one of the storage areas, taking into account constraint (14) leads to the following constraint update, similar to the correction step of Equations (9)–(11):
L k = Σ k + D T D Σ k + D T + σ 2 1 ,
x ^ k = x ^ k + + L k 0 D x ^ k + = I L k D x ^ k + ,
Σ k = I L k D Σ k + ,
where x ^ k and Σ k are the constrained estimated state and covariance, respectively. Considering that D = 0   0   1   0   0 , Equation (16) becomes the following:
L k = 1 Σ k + 3 , 3 + σ 2 Σ k + 1 , 3 Σ k + 5 , 3 .
This result allows for writing down the following:
M k = I L k D = 1 0 α 1 k 0 0 0 1 α 2 k 0 0 0 0 β k 0 0 0 0 α 4 k 1 0 0 0 α 5 k 0 1     with : T k = 1 Σ k + 3 , 3 + σ 2 α i k = T k Σ k + i , 3     i = 1 , 2 , 4 , 5 β k = T k σ 2
Finally, the constraints in scalar form for fast computation become the following:
x ^ k i = x ^ k + i α i k x ^ k + 3     i = 1 , 2 , 4 , 5 x ^ k 3 = β k x ^ k + 3 Σ k i , j = Σ k + i , j α i k Σ k + 3 , j     i = 1 , 2 , 4 , 5 ,   j = 1 , , 5 Σ k 3 , j = β k Σ k + 3 , j       j = 1 , , 5    
Equations (20) and (21) can be optimally computed with 25 additions (1 for Equation (20), 4 for x ^ k , and 4 × 5 for Σ k ) and 36 multiplications (6 for Equation (20), 5 for x ^ k , and 4 × 5 + 5 for Σ k ).
The flowchart of the constrained LBS is summarized in Figure 3.

4. Results

Now, let us apply the proposed approach to a 6000 m 2 warehouse. This warehouse contains m = 2 storage areas and is equipped with N = 6 Wi-Fi access points (reference: Cambium CnPilot E510 with firmware 3.11.4.1-r3). A 425 m journey is done with a robot (see Figure 4), according to the route drawn in Figure 5. The robot is a WiFibot Lab V4 equipped with an Intel Celeron quad core SBC running Linux Ubuntu 18.04 LTS. It embeds for the wireless communication an Atheros AR9280 wireless-N dual band half mini-card. The robot is also equipped with a YEI 3-space IMU, whose noise characteristics are 99 μ g / Hz for the accelerometers and 9 mdeg / s for the gyroscopes, leading to the following P with sample time T s = 0.1   s :
P = 10 8 9.424 0 0 0 9.424 0 0 0 0.247 .    
To set the parameters A i , n   and the variance R , the robot was placed at different known fixed positions all over the warehouse. About 2000 RSSI values coming from the six access points were acquired. Remarking that Equation (4) is linear in parameters A i and n , the least squares method can be applied to compute their value. The variances of the computed A i that can further be obtained are by definition equal to R . The results are n = 3.6885 , R = 1   dB 2 , and:
Acc .  pt   i 1 2 3 4 5 6 A i   ( dB ) 11.66 14.04 17.16 14.52 12.53 16.43
Finally, the additional process noise and velocity constraint covariances were set to:
Q = 10 6 0.01 0 0 0 0 0 0.01 0 0 0 0 0 25 0 0 0 0 0 25 0 0 0 0 0 0.04 ,       σ 2 = 1 m 2 s 2 ,    
and the EKF was initialized with Equation (25).
x ^ 0 = 30 106 0 0 π / 3 x 0     Σ 0 = 3 2 0 0 0 0 0 3 2 0 0 0 0 0 0.1 2 0 0 0 0 0 0.1 2 0 0 0 0 0 5 2 .    
The localization results are given in Figure 6 for the position with respect to the time and in Figure 7 for the 2D trajectory.
The same journey was also done five times at different speeds, with three different devices: two barcode readers and one tablet. During these experiments, the RSSI and IMU measurements were sent to a log server through a 4G connection and a 4G Teltonika RUT950 router, to create a database. The constrained and unconstrained LBS algorithms were then applied using this database. The results in terms of RMS (root mean square) are summarized in Table 1. The length of each journey and the total number of RSSI detected during the route are also given for information.

5. Discussion

Let us now discuss about the results presented in the previous section. It can clearly be seen in Figure 7 that taking into account the speed constraint leads to a better tracking of the trajectory. It is true in the storage areas, but it also has a positive impact outside, because the heading estimate is more accurate. This can especially be observed in the last part of the route, when the robot is leaving the area 2. It can also be noticed that the constraint allows for a faster realignment of the estimated trajectory to the real one after each abrupt 90° turn. Such situations are usually tricky because of the difficulty to accurately estimate the heading.
In Figure 6, the impact of the constraint can be seen in the confidence interval ± 3 × STD . An increase of this interval stands for a localization uncertainty increase, due to the low rate of RSSI detected. It is equivalent to an observability deterioration. Such a phenomenon is from far less important in the constrained case (Figure 6b), more especially along the X axis. It highlights the positive effect of the constraint.
Looking to the figures in Table 1, the first thing to point out is that the worst results were obtained in the shortest journey (the second) and the best ones were in the longest (fourth and fifth). This makes sense because for a given distance, the slower the journey, the longer it is, which involves more RSSI detections. The second thing to notice is that the global RMS over the five journeys and three devices are 1.48 m versus 1.20 m without/with the constraint. The improvement is thus 19%. The best obtained result (device 2, journey 5) is a mean error smaller than 70 cm, which is pretty good considering the warehouse harshness.
To have a point of comparison, our approach was also experimented with the dataset [27] presented in [28]. This dataset included several journeys done by an object (a Nexus 4 running Android 4.4). The data included the real position and the RSSI of the different access points. The route is drawn in Figure 8, extracted from [28].
Unlike in [28], we considered only one MAC address for each access point, i.e., that six RSSI were available at each scan. We also did not consider markers. According to our approach, four areas can be considered, as detailed in Table 2.
The constraint equations suitable with these areas are:
D x k + ζ k = 0     with :   D = 0   0   0   1   0   for   areas   1   and   3 0   0   1   0   0   for   areas   2   and   4 .
The results in Table 3 highlight that even if our approach is slightly less efficient than the Soft Range Limited K-Nearest Neighbors (SRL-KNN) algorithm introduced in [27], it has smaller average errors than all of the other methods. The estimated 2D trajectory is plot and compared with the real one in Figure 9.

6. Conclusions

In this paper, a constrained Kalman filter for Wi-Fi-based indoor localization was developed. This approach is suitable when the space organization is flexible, as in warehouses. The constraint is about the object speed and only applies in some areas of the warehouse. The object to track can be spun around, for example a barcode reader in a hand. The proposed approach was experimented with a robot and three devices on five different journeys. The results highlight that the proposed approach improves the localization accuracy by 19%. As an outlook of this work, the proposed approach can be extended to the 3D case. The use of an estimator such as the UKF and nonlinear smoothing will also be considered. Finally, the regional classifications analysis is another prospect of this work.

Author Contributions

Conceptualization, V.S. and H.C.; methodology, V.S.; software, V.S.; validation, V.S.; formal analysis, V.S.; investigation, V.S.; resources, H.C.; data curation, V.S.; writing—original draft preparation, V.S.; writing—review and editing, V.S.; visualization, V.S. and H.C..; supervision, H.C..; funding acquisition, H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Normandy Region and the European Regional Development Fund within the “GéoWiLog” project.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This work utilized part of the open data [27], which can be found at: https://ieee-dataport.org/open-access/wifi-rssi-indoor-localization (accessed on 21 November 2021).

Acknowledgments

This work is financed within the “GéoWiLog” project. The program is cofounded the Normandy Region and the European Union (Europe is involved in Normandy through the European Funds for Regional Development-ERDF). The authors also gratefully acknowledge the association Twelve Solution and the companies Circoé and Bolloré Logistics for their valuable help in the experimental part of this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, Z.; Zhao, L.; Qin, C.; Wang, Y. Wi-Fi/PDR integrated navigation with robustly constrained Kalman filter. Meas. Sci. Technol. 2020, 31, 84002. [Google Scholar] [CrossRef]
  2. Mautz, R. Indoor Positioning Technologies; Habilitation Thesis; ETH Zurich: Zürich, Switzerland, 2012. [Google Scholar]
  3. Syberfeldt, A.; Ayani, M.; Holm, M.; Wang, L.; Lindgren-Brewster, R. Localizing operators in the smart factory: A review of existing techniques and systems. In Proceedings of the International Symposium on Flexible Automation (ISFA), Cleveland, OH, USA, 1–3 August 2016; pp. 179–185. [Google Scholar]
  4. Zafari, F.; Gkelias, A.; Leung, K.K. A survey of indoor localization systems and technologies. IEEE Commun. Surv. Tutor. 2019, 21, 2568–2599. [Google Scholar] [CrossRef] [Green Version]
  5. Farid, Z.; Nordin, R.; Ismail, M. Recent advances in wireless indoor localization techniques and system. J. Comput. Netw. Commun. 2013, 2013, 185138. [Google Scholar] [CrossRef]
  6. Basri, C.; El Khadimi, A. Survey on indoor localization system and recent advances of Wi-Fi fingerprinting technique. In Proceedings of the 5th International Conference on Multimedia Computing and Systems, Marrakech, Morocco, 29 September–1 October 2016; pp. 253–259. [Google Scholar]
  7. Aly, H.; Youssef, M. New insights into WiFi-based device-free localization. In Proceedings of the 2013 ACM Conference on Pervasive and Ubiquitous Computing Adjunct Publication, Zurich, Switzerland, 8–12 September 2013; pp. 541–548. [Google Scholar]
  8. Tian, Z.; Wang, Z.; Li, Z.; Zhou, M. RTIL: A Real-Time Indoor Localization System by Using Angle of Arrival of Commodity WiFi Signal. In Proceedings of the 11th International Conference on Wireless Communications and Signal Processing (WCSP), Xi’an, China, 23–25 October 2019; pp. 1–6. [Google Scholar]
  9. Shu, Y.; Huang, Y.; Zhang, J.; Coué, P.; Cheng, P.; Chen, J.; Shin, K.G. Gradient-based fingerprinting for indoor localization and tracking. IEEE Trans. Ind. Electron. 2015, 63, 2424–2433. [Google Scholar] [CrossRef]
  10. Xiong, J.; Jamieson, K. Arraytrack: A fine-grained indoor location system. In Proceedings of the 10th {USENIX} Symposium on Networked Systems Design and Implementation ({NSDI}, Lombard, IL, USA, 2–5 April 2013; pp. 71–84. [Google Scholar]
  11. Xu, W.; Liu, L.; Zlatanova, S.; Penard, W.; Xiong, Q. A pedestrian tracking algorithm using grid-based indoor model. Autom. Constr. 2018, 92, 173–187. [Google Scholar] [CrossRef]
  12. Yu, J.; Na, Z.; Liu, X.; Deng, Z. Wi-Fi/PDR-integrated indoor localization using unconstrained smartphones. EURASIP J. Wirel. Commun. Netw. 2019, 41. [Google Scholar] [CrossRef]
  13. Zhang, H.; Xia, Y.; Liu, K.; Jin, F.; Chen, C.; Liao, Y. A Kalman Filter Based Indoor Tracking System via Joint Wi-Fi/PDR Localization. In Proceedings of the IEEE SmartWorld, Ubiquitous Intelligence & Computing, Advanced & Trusted Computing, Scalable Computing & Communications, Cloud & Big Data Computing, Internet of People and Smart City Innovation (SmartWorld/SCALCOM/UIC/ATC/CBDCom/IOP/SCI), Guangzhou, China, 8–12 October 2018; pp. 1444–1449. [Google Scholar]
  14. Chen, Z.; Zou, H.; Jiang, H.; Zhu, Q.; Soh, Y.C.; Xie, L. Fusion of Wi-Fi, smartphone sensors and landmarks using the Kalman filter for indoor localization. Sensors 2015, 15, 715–732. [Google Scholar] [CrossRef]
  15. Zhan, M.; Xi, Z.H. Indoor Location Method of Wi-Fi/PDR Fusion Based on Extended Kalman Filter Fusion. J. Phys. Conf. Ser. 2020, 1601, 042004. [Google Scholar]
  16. Wang, J.; Hu, A.; Li, X.; Wang, Y. An improved PDR/magnetometer/floor map integration algorithm for ubiquitous positioning using the adaptive unscented Kalman filter. ISPRS Int. J. Geo-Inf. 2015, 4, 2638–2659. [Google Scholar] [CrossRef] [Green Version]
  17. Xujian, H.; Hao, W. Wi-Fi indoor positioning algorithm based on improved Kalman filtering. In Proceedings of the International Conference on Intelligent Transportation, Big Data & Smart City (ICITBS), Changsha, China, 17–18 December 2016; pp. 349–352. [Google Scholar]
  18. Li, Z.; Liu, C.; Gao, J.; Li, X. An improved Wi-Fi/PDR integrated system using an adaptive and robust filter for indoor localization. ISPRS Int. J. Geo-Inf. 2016, 5, 224. [Google Scholar] [CrossRef] [Green Version]
  19. Cui, Y.; Zhang, Y.; Huang, Y.; Wang, Z.; Fu, H. Novel Wi-Fi/MEMS integrated indoor navigation system based on two-stage EKF. Micromachines 2019, 10, 198. [Google Scholar] [CrossRef] [Green Version]
  20. Zhang, Y.M.; Sircoulomb, V.; Langlois, N. Robust State Estimation in Networked Control Systems under Linear Hard Equality Constraints. Int. J. Model. Identif. Control 2015, 3, 167–175. [Google Scholar] [CrossRef]
  21. Simon, D. Kalman Filtering with State Constraints: A Survey of Linear and Nonlinear Algorithms. IET Control Theory Appl. 2010, 4, 1303–1318. [Google Scholar] [CrossRef] [Green Version]
  22. Zhou, M.; Dolgov, M.; Liu, Y.; Wang, Y. Wi-Fi/PDR integrated system for 3D indoor localization. In Proceedings of the International Conference on Machine Learning and Intelligent Communications, Hangzhou, China, 6–8 July 2018; pp. 451–459. [Google Scholar]
  23. Yuan, D.; Zhang, J.; Wang, J.; Cui, X.; Liu, F.; Zhang, Y. Robustly Adaptive EKF PDR/UWB Integrated Navigation Based on Additional Heading Constraint. Sensors 2021, 21, 4390. [Google Scholar] [CrossRef]
  24. Zhao, Y.; Li, X.; Wang, Y.; Xu, C.Z. Biased constrained hybrid Kalman filter for range-based indoor localization. IEEE Sens. J. 2017, 18, 1647–1655. [Google Scholar] [CrossRef]
  25. Gao, Z.; Guo, H.; Xie, Y.; Lu, H.; Zhang, J.; Diao, W.; Xu, R. An improved localization method in cyber-social environments with obstacles. Comput. Electr. Eng. 2020, 86, 106694. [Google Scholar] [CrossRef]
  26. Porrill, J. Optimal combination and constraints for geometrical sensor data. Int. J. Robot. Res. 1988, 7, 66–77. [Google Scholar] [CrossRef]
  27. Hoang, M.T.; Dong, X.; Lu, T.; Yuen, B.; Westendorp, R. WiFi RSSI Indoor Localization. IEEE Dataport, 30 November 2019. [Google Scholar] [CrossRef]
  28. Hoang, M.T.; Zhu, Y.; Yuen, B.; Rees, T.; Dong, X.; Lu, T.; Westendorp, R.; Xie, M. A soft range limited K-nearest neighbours algorithm for indoor localization enhancement. IEEE Sens. 2018, 18, 10208–10216. [Google Scholar] [CrossRef] [Green Version]
  29. Zou, H.; Jin, M.; Jiang, H.; Xie, L.; Spanos, C.J. WinIPS: WiFibased non-intrusive indoor positioning system with online radio map construction and adaptation. IEEE Trans. Wirel. Commun. 2017, 16, 8118–8130. [Google Scholar] [CrossRef]
  30. Xie, Y.; Wang, Y.; Nallanathan, A.; Wang, L. An improved K-Nearest-Neighbor indoor localization method based on Spearman distance. IEEE Signal Process. Lett. 2016, 23, 351–355. [Google Scholar] [CrossRef] [Green Version]
  31. Kushki, A.; Plataniotis, K.N.; Venetsanopoulos, A.N. Kernel based positioning in wireless local area networks. IEEE Trans. Mob. Comput. 2007, 6, 689–705. [Google Scholar] [CrossRef] [Green Version]
  32. Au, A.W.S.; Feng, C.; Valaee, S.; Reyes, S.; Sorour, S.; Markowitz, S.N.; Gold, D.; Gordon, K.; Eizenman, M. Indoor tracking and navigation using received signal strength and compressive sensing on a mobile device. IEEE Trans. Mob. Comput. 2013, 12, 2050–2062. [Google Scholar] [CrossRef]
Figure 1. Moving device to track into an indoor environment.
Figure 1. Moving device to track into an indoor environment.
Sensors 22 00428 g001
Figure 2. Varying organization of the warehouse, with m storage areas. An example of shelves’ organization is drawn for area m .
Figure 2. Varying organization of the warehouse, with m storage areas. An example of shelves’ organization is drawn for area m .
Sensors 22 00428 g002
Figure 3. Flowchart of the proposed constrained LBS approach.
Figure 3. Flowchart of the proposed constrained LBS approach.
Sensors 22 00428 g003
Figure 4. WiFibot moving into the warehouse.
Figure 4. WiFibot moving into the warehouse.
Sensors 22 00428 g004
Figure 5. Warehouse under consideration with the 425 m journey done by the robot and the localization of the different Wi-Fi access points and storage areas.
Figure 5. Warehouse under consideration with the 425 m journey done by the robot and the localization of the different Wi-Fi access points and storage areas.
Sensors 22 00428 g005
Figure 6. Position with respect to the time: (a) without application of constraints; (b) with application of the constraints. The right part of plots (a,b) are the error (difference between real and estimated position) and the confidence interval ± 3 × STD , where STD is the standard deviation: STD = Σ k 1 , 1 for the X-axis and STD = Σ k 2 , 2 for the Y-axis.
Figure 6. Position with respect to the time: (a) without application of constraints; (b) with application of the constraints. The right part of plots (a,b) are the error (difference between real and estimated position) and the confidence interval ± 3 × STD , where STD is the standard deviation: STD = Σ k 1 , 1 for the X-axis and STD = Σ k 2 , 2 for the Y-axis.
Sensors 22 00428 g006
Figure 7. Comparison of the 2D positions with and without constraint applications.
Figure 7. Comparison of the 2D positions with and without constraint applications.
Sensors 22 00428 g007
Figure 8. Experiment in the engineering office wing at University of Victoria, considering open datasets [27]. The six access points were located at the AP circles.
Figure 8. Experiment in the engineering office wing at University of Victoria, considering open datasets [27]. The six access points were located at the AP circles.
Sensors 22 00428 g008
Figure 9. Real (ground truth) trajectory and the estimated one with our approach.
Figure 9. Real (ground truth) trajectory and the estimated one with our approach.
Sensors 22 00428 g009
Table 1. RMS localization errors (in m) with and without constraint applications.
Table 1. RMS localization errors (in m) with and without constraint applications.
Device 1Device 2Device 3
Journey 1
Length: 8′09″
Nbr of RSSI detect.: 105
RMS unconstr.: 1.86
RMS constrained: 1.05
Nbr of RSSI detect.: 155
RMS unconstr.: 1.31
RMS constrained: 1.28
Nbr of RSSI detect.: 107
RMS unconstr.: 1.53
RMS constrained: 1.36
Journey 2
Length: 6′34″
Nbr of RSSI detect.: 127
RMS unconstr.: 1.87
RMS constrained: 1.43
Nbr of RSSI detect.: 94
RMS unconstr.: 2.26
RMS constrained: 1.68
Nbr of RSSI detect.: 114
RMS unconstr.: 1.59
RMS constrained: 1.58
Journey 3
Length: 8′25″
Nbr of RSSI detect.: 115
RMS unconstr.: 1.76
RMS constrained: 1.04
Nbr of RSSI detect.: 151
RMS unconstr.: 0.93
RMS constrained: 0.89
Nbr of RSSI detect.: 244
RMS unconstr.: 1.01
RMS constrained: 0.94
Journey 4
Length: 10′38″
Nbr of RSSI detect.: 90
RMS unconstr.: 1.43
RMS constrained: 1.01
Nbr of RSSI detect.: 105
RMS unconstr.: 1.23
RMS constrained: 0.96
Nbr of RSSI detect.: 142
RMS unconstr.: 1.56
RMS constrained: 1.38
Journey 5
Length: 10′34″
Nbr of RSSI detect.: 150
RMS unconstr.: 1.31
RMS constrained: 1.12
Nbr of RSSI detect.: 123
RMS unconstr.: 0.77
RMS constrained: 0.66
Nbr of RSSI detect.: 119
RMS unconstr.: 0.89
RMS constrained: 0.75
Table 2. Constraint areas in Figure 8.
Table 2. Constraint areas in Figure 8.
Area 1
(Including AP3)
Area 2
(Including AP4)
Area 3
(Including AP5)
Area 4
(Including AP6)
X min 6 m16.5 m6 m0 m
X max 16.5 m21 m16.5 m6 m
Y min 2 m−11 m−14 m−11 m
Y max 4 m2 m−11 m2 m
Table 3. Comparison of the average errors.
Table 3. Comparison of the average errors.
Our approach: 0.89 mSRL-KNN [28]: 0.66 m to 1.2 mSTI-WKNN [29]: 1.09 m
Spearman rank [30]: 1.45 mKernel method [31]: 1.07 mKalman [32]: 0.96 m
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sircoulomb, V.; Chafouk, H. A Constrained Kalman Filter for Wi-Fi-Based Indoor Localization with Flexible Space Organization. Sensors 2022, 22, 428. https://doi.org/10.3390/s22020428

AMA Style

Sircoulomb V, Chafouk H. A Constrained Kalman Filter for Wi-Fi-Based Indoor Localization with Flexible Space Organization. Sensors. 2022; 22(2):428. https://doi.org/10.3390/s22020428

Chicago/Turabian Style

Sircoulomb, Vincent, and Houcine Chafouk. 2022. "A Constrained Kalman Filter for Wi-Fi-Based Indoor Localization with Flexible Space Organization" Sensors 22, no. 2: 428. https://doi.org/10.3390/s22020428

APA Style

Sircoulomb, V., & Chafouk, H. (2022). A Constrained Kalman Filter for Wi-Fi-Based Indoor Localization with Flexible Space Organization. Sensors, 22(2), 428. https://doi.org/10.3390/s22020428

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop