Next Article in Journal
Low-Temperature Storage Improves the Over-Time Stability of Implantable Glucose and Lactate Biosensors
Next Article in Special Issue
UWB Localization with Battery-Powered Wireless Backbone for Drone-Based Inventory Management
Previous Article in Journal
Ripeness Prediction of Postharvest Kiwifruit Using a MOS E-Nose Combined with Chemometrics
Previous Article in Special Issue
Indoor 3-D Localization Based on Received Signal Strength Difference and Factor Graph for Unknown Radio Transmitter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Indoor Positioning System Based on Chest-Mounted IMU

1
Graduate School of Information Science and Electrical Engineering, Kyushu University, Fukuoka 819-0395, Japan
2
Library, Kyushu University, Fukuoka 819-0395, Japan
3
Faculty of Information Science and Electrical Engineering, Kyushu University, Fukuoka 819-0395, Japan
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(2), 420; https://doi.org/10.3390/s19020420
Submission received: 4 January 2019 / Revised: 15 January 2019 / Accepted: 16 January 2019 / Published: 21 January 2019
(This article belongs to the Special Issue Sensor Fusion and Novel Technologies in Positioning and Navigation)

Abstract

:
Demand for indoor navigation systems has been rapidly increasing with regard to location-based services. As a cost-effective choice, inertial measurement unit (IMU)-based pedestrian dead reckoning (PDR) systems have been developed for years because they do not require external devices to be installed in the environment. In this paper, we propose a PDR system based on a chest-mounted IMU as a novel installation position for body-suit-type systems. Since the IMU is mounted on a part of the upper body, the framework of the zero-velocity update cannot be applied because there are no periodical moments of zero velocity. Therefore, we propose a novel regression model for estimating step lengths only with accelerations to correctly compute step displacement by using the IMU data acquired at the chest. In addition, we integrated the idea of an efficient map-matching algorithm based on particle filtering into our system to improve positioning and heading accuracy. Since our system was designed for 3D navigation, which can estimate position in a multifloor building, we used a barometer to update pedestrian altitude, and the components of our map are designed to explicitly represent building-floor information. With our complete PDR system, we were awarded second place in 10 teams for the IPIN 2018 Competition Track 2, achieving a mean error of 5.2 m after the 800 m walking event.

1. Introduction

In outdoor positioning systems, GPS and magnetometers are assembled to estimate the position and orientation of a system with respect to Earth’s frame. With a standard GPS on mobile, positioning error can be around 10 m in well-conditioned outdoor environments, such as an open-space area. However, accuracy largely becomes less reliable in indoor environments because not only do the walls of the building block the GPS signals, but also additional magnetic fields from electrical devices produce noises. Therefore, localization using such devices generally works only in outdoor situations.
Nowadays, demand for indoor navigation systems has been rapidly increasing for location-based services, such as navigation to shops, and advertisements. From a technical point of view, indoor navigation systems can be separated into two categories: infrastructure-based systems, and infrastructure-free ones. The first category requires to install devices in the environment, and the position of the pedestrian is computed by using those devices. Radio wave-based systems, such as radio-frequency identification (RFID) [1,2], bluetooth low energy (BLE) [3], and WiFi [4,5], belong in this category. These normally need many resources to calibrate the devices so that all device positions are known. The second category, which does not require any additional hardware to be installed, can be more cost-effective. Camera-based simultaneous localization and mapping (SLAM) and inertial measurement unit (IMU)-based pedestrian dead reckoning (PDR) [6] systems belong in this class, and can provide the relative pose from a starting location. Compared with the first category, they need the initial pose be known by other means, whereas both the sampling rate and the accuracy of the relative pose estimation are quite high. With the recent advances of cameras and SLAM algorithms, SLAM-based navigation systems can work in real time. However, a drawback of camera-based systems is requiring both bright illumination and a static environment, which sometimes cannot be fulfilled in a real environment, such as a shopping mall. Although accuracy is not as high as the SLAM systems, IMU-based PDR systems have no impact on the environment, and have been researched for years.
In the literature of the IMU-based PDR, the IMU can be mounted on different body parts to adapt to various conditions. As the most accurate way, the IMU was mounted on a foot to utilize the framework of the zero-velocity update (ZUPT) [7,8]. The velocity computed by sequentially integrating acceleration normally drifts away from actual value over time due to the accumulation of sensor noise. The ZUPT was proposed to reduce this error as follows. While walking, the velocity of a foot can be considered as zero when the foot is attached on the ground. By detecting these periods and resetting the velocity to zero, the error on the velocity can be suppressed. However, the ZUPT cannot be used for other body parts because there is normally no periodical zero-velocity moment.
A specific motion model was normally proposed according to the body installation position of the IMU. For example, handheld IMU-based systems can be implemented in smartphones [9]. Head-mounted systems can be implemented in smartglasses [10]. Some systems can also recognize where the IMU is mounted and apply different processes according to positions [11,12]. However, it is not easy to accurately compute walking displacement from the IMU attached onto such body parts because there are many unexpected movements in addition to pure walking movement, such as head rotation and hand shaking. Unexpected movements can obviously introduce errors and decrease positioning accuracy.
In this paper, we propose an indoor positioning system that uses a chest-mounted IMU with an efficient map-matching algorithm for multifloor navigation tasks. We propose to select the chest as a novel installation position because the chest has the fewest unexpected movements in the upper body while walking. In other words, the chest has lower degrees of freedom for movement compared with other parts, such that it normally faces the moving direction. Therefore, our motivation is to create a body-suit-type PDR system to achieve navigation in a multifloor building. The input data of our system are the temporal sequence of acceleration, angular velocity, and barometric pressure from the IMU, and the floor map for map matching. As a device, we used NGIMU (http://x-io.co.uk/ngimu/), of which the specifications are listed in Table 1, to develop a prototype system, and attached it at the chest with a stretching band, as illustrated in Figure 1. To adapt the PDR system to the chest-mounted IMU, we propose a novel regression model based on step-length estimation. Different from methods for other upper-body-based PDR, our approach can accurately compute not only step length but also step direction. To reduce error accumulation on the position and heading of the pedestrian in the building, we integrated the idea of an efficient map-matching algorithm based on particle filtering into our system. The components of our map were carefully designed, and can explicitly represent building-floor information. As onsite performance evaluation, we were awarded second place in 10 teams for the IPIN 2018 Competition (http://ipin2018.ifsttar.fr/competition/about/) Track 2, achieving a mean error of 5.2 m after walking 800 m.
In summary, the contribution of our paper is threefold as follows.
  • A novel chest-mounted IMU-based PDR for body-suit-type systems is proposed.
  • Step-length estimation for a chest-mounted IMU is proposed.
  • A complete system for multifloor navigation tasks was implemented and the code is open (https://github.com/rairyuu/PDR-with-Map-Matching).

2. Related Work

The framework of IMU-based PDR systems generally consists of four parts: orientation update, step detection, step-length estimation, and map matching. Orientation update and map matching are common techniques in navigation systems [13], while step detection and step-length estimation are dedicated to human motion because they use the constraints of a pedestrian step. In this section, we review the existing methods on these three parts.
As the cornerstone of navigation systems, orientation update has been researched for decades. Orientation accuracy heavily affects positioning accuracy. Generally, IMU orientation is computed from gyroscope data. Therefore, its accuracy depends on the accuracy of the gyroscope. In practice, the error on gyroscope data consists of two parts, the high-frequency noise and the slow-changing temperature drift [14]. Many studies have explored methods to reduce noise [15,16]. To suppress the temperature drift, Shen et al. proposed a model based on a genetic algorithm and the Elman neural network [17]. However, the computation cost of removing sensor noise and the temperature drift is relatively high. To achieve a fast orientation update, Madgwick et al. proposed an efficient orientation filter that is robust against the error on sensor data [18]. By utilizing the accelerometer data, the Madgwick filter achieves a high-accuracy orientation update. Moreover, its computation cost is low. Therefore, in our system, we used the Madgwick filter to update the orientation.
The input data for standard IMU-based PDR systems are the temporal sequence of acceleration and angular velocity expressed in the device frame. By computing the gravity direction and using the double integration of the acceleration with the initial velocity after removing gravity acceleration, it is possible to compute the displacement on the 2D map. However, the computed displacement is not accurate over time owing to error accumulation caused by sensor noise. To avoid error accumulation, the process of step detection was proposed. For foot-mounted IMU-based systems, Ojeda et al. proposed two empirical rules to separate steps [7,8]. Madgwick et al. proposed a step-detection algorithm based on filtering and thresholding, which only required the norm of the acceleration [18,19]. For handheld IMU-based systems, step detection is more difficult because hand motion is not always correlated with step motion. To deal with this problem, Susi et al. proposed an adaptive algorithm that first recognized the motion type, then applied a proper method according to the motions [20]. For head-mounted IMU-based systems, a filtering and thresholding-based algorithm is still a good choice, as implemented by Zhang et al. [10].
Step-length estimation is the process of computing the displacement of each step. Displacement consists of two parts: length and direction. This means that this process can be classified into two categories: computing the length and direction, separately or not. The use of a fixed step length computed from, for instance, body height, may work for computing the displacement of normal walking. However, it is generally not accurate for other motions such as running or slow steps. To improve the accuracy, Weinberg et al. proposed an empirical equation to compute the step length with a hip (or upper-body)-mounted IMU [21]. Step direction can be computed from the IMU poses. For foot-mounted IMU-based systems, the double integration of the acceleration can produce both length and direction together. However, this straightforward computation cannot fit all motions. To settle this problem, Shin et al. proposed an adaptive step-length estimation algorithm [22]. For handheld IMU systems, Renaudin et al. proposed an adaptive algorithm based on motion mode classification [23]. Yan et al. proposed a regression model to compute step displacement [11].
Utilizing both step detection and step-length estimation can avoid errors in each step. However, position and heading errors in the global or world frame still accumulate over time. To solve this problem, a map-matching algorithm was proposed by making the constraints of the positions in the global map frame [24,25,26]. The particle filter-based approach is particularly good for dealing with nonlinear and non-Gaussian estimation problems [27]. By implementing the particle filter for a map-matching algorithm, positioning accuracy can be significantly improved. Davidson et al. showed how the particle filter improved the performance of indoor navigation systems [28]. By combining the backtracking particle filter with different levels of a 2D building map, Widyawan et al. achieved highly accurate indoor navigation [29]. Bojja et al. proposed a 3D map-matching algorithm in order to navigate the 3D space [30].

3. Proposed System

As illustrated in Figure 2, our system mainly consists of four parts: orientation update, step detection, step-length estimation, and map matching. This framework is basically the same as the standard used in IMU-based PDR systems. Particularly, we propose novel step-length estimation with a chest-mounted IMU, and we implemented a complete system for multifloor navigation tasks.
In the orientation update, IMU orientation with respect to the world frame, which is the initial IMU sensor frame, was computed. Step detection is the process of separating sequential data into steps. Once one step is detected, the displacement of the step is computed in step-length estimation. Since only the PDR with the IMU contains the positioning and heading errors, consistency between the estimated trajectory and the map is optimized in map matching. In the remainder of this section, we introduce each process in detail.

3.1. Orientation Update

The data acquired from the IMU are expressed in the IMU sensor frame. To compute the current position in the world frame, we first transform the 3D acceleration data from the sensor frame to the world frame. Normally, we defined the initial IMU sensor frame as the world frame, as illustrated in Figure 3. For instance, the z axis can be parallel to the gravity direction when the IMU is placed on the ground at the beginning of the process. In our system, the world frame is defined as follows: the x axis is the initial moving direction, the y axis is perpendicular to the x axis, and the z axis is parallel to the gravity direction. The transformation matrix is simply computed with the algorithm by Madgwick et al. [18], which uses Quaternion representation for three degrees-of-freedom (DoF) poses. After this process, we can have acceleration in the world frame to remove gravitational acceleration on the z axis, and finally have the acceleration of the device itself.
It should be noted that the orientation computed with the IMU gradually contains errors on the azimuth angle in the world frame, whereas the elevation angle can be compensated with gravity acceleration. For instance, the x axis, which is initially defined as the heading, does not match with the actual heading. To compensate this, we use a map-matching algorithm to correct the heading, as explained in Section 3.4. Since we additionally parameterize the correction value against the orientation, we can use any orientation-update method as an external orientation estimator. This means that our system is based on a loosely coupled approach between the PDR and map matching.

3.2. Step Detection

To reduce the error accumulation caused by sensor noise, we apply the step-detection method to separate sequential data into steps so that we can independently process the data of each step. As illustrated in Figure 4, we propose to use the filtering and thresholding approach for the acceleration norm to detect steps, similarly to References [18,19] because the peaks in the acceleration-data sequence captured at the chest has similar tendencies as other body parts.
First, we compute the 3D acceleration norm in the sequence, as illustrated in Figure 4a. To remove noise, we apply a low-pass filter to the norm sequence. In Figure 4b, the peaks of the result correspond to the moments when the foot hits the ground, where acceleration is maximum. As illustrated in Figure 5, we define one step as the interval from one peak to the next one. An example result of the step detection is illustrated in Figure 4c.

3.3. Step-Length Estimation

In the literature of the PDR with an IMU attached on an upper-body part, step-length estimation was achieved by designing an empirical motion model derived from IMU data [21]. Therefore, we followed the same way to derive the motion model for a chest-mounted IMU.
By analyzing the relationship between displacement and the IMU data, we empirically formulated that the displacement of one step can be computed by:
D = K · a ( t ) d t d t
where a ( t ) is the 3D acceleration at time t, the range of the integral is the duration of one step, and K is a scalar parameter that needs to be calibrated for each user. To avoid drift errors, velocity is not sequentially computed, and Equation (1) takes the initial zero velocity for distance computation. Since the double integration of accelerations is used in our approach, both the length and direction of one step can be computed at the same time.
The motivation of formulating the Equation (1) is as follows. The 2D acceleration data of one step in the x y co-ordinates of the world frame is illustrated in Figure 6a. In this case, the moving-forward direction is along the x axis. From the figure, the decelerating phase is obviously more observed than the accelerating phase. We found that this can be a chest-acceleration characteristic. Therefore, the displacement computed by the double integration of the sequential 2D accelerations is opposite the forward direction, as illustrated in Figure 6b. To transform it to correct the direction and length, a user parameter K is required.
As a calibration method, K is computed for each pedestrian by:
K = D g D d
where D g is a given distance in which the pedestrian is asked to walk straightly, and D d is the total length computed from the double integration of accelerations in one step for all of the steps. To decrease the error on the measurement, we tested five times, and used the mean as the final K. For calibration, the longer D g is, the better for accuracy it is. In our experiment, we empirically set D g to 20 m. For a pedestrian whose height is about 170 cm, K was normally in the range of [ 5.0 , 3.0 ] . After calibration, the error on step length was generally around 10% in our experiments. This error range is used for map matching, as described in Section 3.4.
It should be noted that displacement on the y axis was not zero, as illustrated in Figure 6b. Normally, steps are alternately made by the left foot and right foot. Therefore, the upper body tends to alternately lean to the left side and right side. This is the reason why there is side displacement in one step. Since the step made by the left or right foot is usually followed by a step made by the other foot, the overall movement is forward.
For PDR systems based on an upper-body-mounted IMU, it is usually difficult to compute the step direction. Most of the systems require that the step direction is equivalent to the IMU heading. This means that it is assumed that the IMU heading is equivalent to the moving direction. This constraint makes the systems inconvenient to use, and it is not valid when unexpected movements occur. In contrast, our approach does not depend on the actual IMU heading because step displacement is computed from accelerations in the world frame. Therefore, our system has no restriction on step direction with respect to the IMU heading, which is more user-friendly and accurate under several motions.

3.4. Map Matching

By applying step detection to the IMU data, the error accumulation caused by the sensor noise can be suppressed compared with the sequential double integration of the accelerations. However, the positioning and heading errors in the world frame still grow over time. To improve accuracy, we integrated the idea of an efficient map-matching algorithm based on particle filtering [24] into our system.
As illustrated in Figure 7, we propose to define six components in our map for a multifloor building. l i n e corresponds to a wall that cannot be crossed. In z o n e _ n o _ p a r t i c l e where pedestrians cannot come into, the existing particles are deleted, and new particles cannot be generated. Only in z o n e _ u p d a t e _ a l t i t u d e , where altitude can be changed, such as stairs, escalators, and elevators, will the altitude be updated according to the movement. a r r o w is used for explicitly describing the direction of the possible movement at such locations. a r r o w _ u p connects to the upper floor, a r r o w _ d o w n connects to the lower floor, and a r r o w _ b o t h connects to both floors.
In our map-matching algorithm, each particle has three parameters: position P i , scale s i , and heading correction h i for the i-th particle. Scale s i is a scalar value to have error tolerance of the step length, which is randomized from the range [ 1.0 r , 1.0 + r ] , where r controls the range. r can be determined with the possible error range from preliminary calibration experiments. For instance, r = 0.1 when error is up to 10%, as described in Section 3.3. Since the heading in the world frame may drift away from its actual value, as discussed in Section 3.1, we propose to parameterize heading correction h i as a rotation angle to rectify the current heading computed from the IMU orientation into the actual heading in the x y co-ordinate system, as illustrated in Figure 8.
For each particle, position P i is updated by:
P i t + 1 = P i t + R ( D t · s i , h i )
where D t is the 2D displacement of one current step computed from the IMU, R ( X , y ) is a function that rotates vector X with angle y on the x y co-ordinates. h i and s i are determined and fixed once a new particle is generated. Final position P and heading correction h of the pedestrian are finally computed by:
P = 1 n · i = 1 n P i
h = 1 n · i = 1 n h i
where n is the number of particles.
The procedure of generating new particles is as follows. After updating P i with Equation (3) at every step, our system applies collision detection with the map to the particles. For instance, if the line from P i t to P i t + 1 intersects with the l i n e , the ith particle is deleted, and our system generates new particles for the next step, as described in Algorithm 1 and Table 2. The position of N e w P a r t i c l e is randomized from the circle with center P a r t i c l e and radius R a d i u s . The heading correction of N e w P a r t i c l e is randomized around final heading correction h. The scale of N e w P a r t i c l e is randomized from the range of [ 1.0 r , 1.0 + r ] . Once N e w P a r t i c l e is generated, its heading correction and scale are fixed.
Algorithm 1: Generating new particles.
P a r t i c l e N u m b e r M a x = 256
T r y T i m e M a x = 8
R a d i u s = 3.0
P a r t i c l e N u m b e r = n
while P a r t i c l e N u m b e r P a r t i c l e N u m b e r M a x :
      P a r t i c l e = S e l e c t R a n d o m P a r t i c l e ( )
      T r y T i m e = 1
     while T r y T i m e T r y T i m e M a x :
           N e w P a r t i c l e = P r o p o s e N e w P a r t i c l e ( P a r t i c l e , R a d i u s )
          if B a c k t r a c k i n g T e s t ( N e w P a r t i c l e ) = = P A S S :
                A p p e n d P a r t i c l e ( N e w P a r t i c l e )
                P a r t i c l e N u m b e r = P a r t i c l e N u m b e r + 1
               break
           T r y T i m e = T r y T i m e + 1
When generating new particles, the particles may be generated in inappropriate positions. This problem decreases positioning accuracy, and sometimes even crashes the map-matching module. Therefore, inspired by the backtracking particle filter proposed by Widyawan et al. [31], we applied a backtracking test to new particles to select only appropriate particles. In Figure 9, blue arrows correspond to recent steps. First, new particles are randomly proposed around a randomly selected particle. Then, new particles go back recent steps and apply collision detection to the particles, as illustrated in Figure 9c. If a particle hits the wall, it is deleted. As illustrated in Figure 9d, only particles that pass this test are left. By using the backtracking test, most new particles are produced around the correct position. This algorithm largely allows our system to achieve higher accuracy. Moreover, this algorithm could be used when map matching failed to track the user. By applying the backtracking test to randomly selected points around the last estimated position, it is possible to find the correct current position and restore the system. In our experiment, considering real-time processing and accuracy, we used the 32 most recent steps for the backtracking test.
For multifloor-building tasks, we use the barometer to update the altitude of the pedestrian. Delta altitude Δ a can be computed by:
Δ a = Δ p · 0.09 m / Pa
where Δ p is the delta pressure. As illustrated in Figure 10, pressure decreases when going upstairs. However, barometric pressure varies over time even in indoor situations. To correctly update the altitude, we propose to use z o n e _ u p d a t e _ a l t i t u d e as a component of the map. Only in this zone can the altitude be updated. All a r r o w s are defined in this zone. Since the pedestrian does not stay in this zone for a long time, this approach usually works accurately enough. When delta altitude Δ a is closed to the floor height, our system matches the current position to the export of the nearest a r r o w , and produces new particles. Every time the floor changes, the positioning error of our system can be reset.

3.5. Map Editor

To easily generate the map, we developed a map editor so that the map input to our system can be manually edited. The input to the editor is the map images, and the output is a set of files that contain the multifloor-building information. As illustrated in Figure 11, the basic component of our map is a point. By linking the points, we can create any component. Our map editor is available online (https://github.com/rairyuu/PDR-with-Map-Matching).

4. Evaluation

In this section, we present the quantitative-experiment results of our system.

4.1. Evaluation Criterion

We followed the evaluation procedure used in IPIN 2018 Track 2 as follows. In each experiment, the pedestrian was asked to walk along a given route. Several keypoints, where the positions were measured in order to be known, were set in the route before the experiment. The pedestrian was asked to slightly stop at each keypoint so that the estimated position output by our system was recorded.
Error e i of the i-th keypoint is defined by the following equations:
e i = D x y ( P g i , P e i ) + D z ( P g i , P e i )
D x y ( P g i , P e i ) = ( P g i · x P e i · x ) 2 + ( P g i · y P e i · y ) 2
D z ( P g i , P e i ) = | P g i · z P e i · z | · p z
where P g i is the ground-truth position, and P e i is the estimated position. P i has three elements, x, y, and z. x and y are expressed in meters, whereas element z is a scalar that represents the current floor. p z is the penalty on floor error, which is set to 15.

4.2. System Initialization

In each experiment, the initial position and heading were given to set up the system as follows.
  • Open the IMU, launch our system software, and connect the IMU to the computer.
  • Calibrate the IMU pose with the known initial heading in the world frame of the map, as illustrated in Figure 12.
  • Attach the IMU on the pedestrian chest, and run the system.

4.3. Lab Experiment

In the first experiment, we evaluated our system in the campus building with different configurations. As illustrated in Figure 13a, the pedestrian was asked to start at the pink point, go along the orange arrows, walk around the green route for three times, and finally walk back to the starting point along the blue arrows. There was no open-space area in the route, and the width of corridors was about 3 m. In this experiment, we used 256 particles. The travelled distance was 432.22 m. The red points represent the keypoints for the evaluation. Keypoint positions were given by the map of the building.
The estimated routes of our system are illustrated as Figure 13c,d. At the beginning of the route, the system without map matching was accurate enough. However, the drift on heading was growing over time, as can be seen in the figure. Therefore, the result obviously became worse. In our system with map matching, since the estimated position and heading were calibrated dynamically, the result was much more accurate and stable. The detailed results are illustrated in Figure 14 and Table 3.

4.4. IPIN 2018 Competition Experiment

In the second experiment, we evaluated our system by participating in the IPIN 2018 Competition Track 2. The competition was held in a huge shopping mall that contained several open-space areas. The evaluation route crossed three floors, and consisted of indoor and outdoor spaces. In the indoor-space areas, the width range of corridors was from 10 to 20 m. In this experiment, we used 512 particles that were optimized by preliminary experiments. During the competition, the shopping mall was crowded.
The competition route, our estimated result, and the errors are illustrated in Figure 15. The error at each keypoint is illustrated in Figure 16. In outdoor-space areas, since the map matching cannot work, our result was not so accurate. When we went into an indoor area, the result improved, where error was generally less than 5 m. As described in Table 4, the mean error of our system was 5.2 m. The IPIN competition evaluated the 75th percent error. Therefore, our final score was 5.7 m, which was 0.2 m worse than the champion. In the 10 teams who participated in this track, we were awarded second place. The detail of the competition can be found on the IPIN 2018 homepage (http://ipin2018.ifsttar.fr/competition/about/). We also uploaded a video (https://youtu.be/Ti7w-r3dB2E) that provides a more detailed result and shows how our map matching worked in the competition.

5. Conclusions

In this paper, we proposed a novel indoor positioning system with map matching. In our system, the IMU is mounted on the chest where movement is simplest and most stable. To reduce error accumulation on position and heading, we designed an efficient map-matching algorithm based on particle filtering. Different from most existing algorithms, our map-matching algorithm was designed for 3D navigation, which can be used in multifloor buildings. With carefully designed components (wall, stair, escalator, elevator), our map can clearly represent building-floor information. With the proposed algorithm implemented, our navigation system was awarded second place in the IPIN 2018 Competition Track 2, achieving a mean error of 5.2 m after 800 m of walking. In our future work, we are improving the accuracy of the step-length estimation. As mentioned in our motivation, we are also designing a body suit to support the indoor work.

Author Contributions

All authors conceived and designed the study. C.L. designed and implemented the system, developed the tool programs, and drafted most of the manuscript under the supervision of and with suggestions from H.U., D.T., A.S., and R.-i.T. C.L. and H.U. participated in the IPIN 2018 Competition Track 2. All authors contributed to the result analysis and discussions. All authors approved the submitted manuscript.

Funding

Part of this work was supported by JSPS KAKENHI, Grant Number JP18H04125.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; and in the decision to publish the results.

References

  1. Finkenzeller, K. RFID Handbook: Fundamentals and Applications in Contactless Smart Cards, Radio Frequency Identification and Near-Field Communication; John Wiley & Sons: Hoboken, NJ, USA, 2010. [Google Scholar]
  2. Ni, L.M.; Liu, Y.; Lau, Y.C.; Patil, A.P. LANDMARC: indoor location sensing using active RFID. In Proceedings of the First IEEE International Conference on Pervasive Computing and Communications, Fort Worth, TX, USA, 26–26 March 2003; pp. 407–415. [Google Scholar]
  3. Newman, N. Apple iBeacon technology briefing. J. Direct Data Digit. Mark. Pract. 2014, 15, 222–225. [Google Scholar] [CrossRef] [Green Version]
  4. Yang, C.; Shao, H.R. WiFi-based indoor positioning. IEEE Commun. Mag. 2015, 53, 150–157. [Google Scholar] [CrossRef]
  5. Woo, S.; Jeong, S.; Mok, E.; Xia, L.; Choi, C.; Pyeon, M.; Heo, J. Application of WiFi-based indoor positioning system for labor tracking at construction sites: A case study in Guangzhou MTR. Autom. Constr. 2011, 20, 3–13. [Google Scholar] [CrossRef]
  6. Jimenez, A.R.; Seco, F.; Prieto, C.; Guevara, J. A comparison of pedestrian dead-reckoning algorithms using a low-cost MEMS IMU. In Proceedings of the 2009 IEEE International Symposium on Intelligent Signal Processing, Budapest, Hungary, 26–28 August 2009; pp. 37–42. [Google Scholar]
  7. Ojeda, L.; Borenstein, J. Non-GPS navigation for security personnel and first responders. J. Navig. 2007, 60, 391–407. [Google Scholar] [CrossRef]
  8. Ojeda, L.; Borenstein, J. Personal dead-reckoning system for GPS-denied environments. In Proceedings of the 2007 IEEE International Workshop on Safety, Security and Rescue Robotics, Rome, Italy, 27–29 September 2007; pp. 1–6. [Google Scholar]
  9. Antigny, N.; Servières, M.; Renaudin, V. Pedestrian Track Estimation with Handheld Monocular Camera and Inertial-Magnetic Sensor for Urban Augmented Reality. In Proceedings of the 2017 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Sapporo, Japan, 18–21 September 2017. [Google Scholar]
  10. Zhang, Y.; Hu, W.; Xu, W.; Wen, H.; Chou, C.T. NaviGlass: Indoor Localisation Using Smart Glasses. In Proceedings of the 2016 International Conference on Embedded Wireless Systems and Networks, EWSN, Graz, Austria, 15–17 February 2016; pp. 205–216. [Google Scholar]
  11. Yan, H.; Shan, Q.; Furukawa, Y. RIDI: Robust IMU Double Integration. arXiv, 2017; arXiv:1712.09004. [Google Scholar]
  12. Elhoushi, M.; Georgy, J.; Korenberg, M.; Noureldin, A. Robust motion mode recognition for portable navigation independent on device usage. In Proceedings of the Position, Location and Navigation Symposium-PLANS 2014, Monterey, CA, USA, 5–8 May 2014; pp. 158–163. [Google Scholar]
  13. Quddus, M.A.; Ochieng, W.Y.; Noland, R.B. Current map-matching algorithms for transport applications: State-of-the art and future research directions. Transp. Res. Part C Emerg. Technol. 2007, 15, 312–328. [Google Scholar] [CrossRef] [Green Version]
  14. Shen, C.; Yang, J.; Tang, J.; Liu, J.; Cao, H. Note: Parallel processing algorithm of temperature and noise error for micro-electro-mechanical system gyroscope based on variational mode decomposition and augmented nonlinear differentiator. Rev. Sci. Instrum. 2018, 89, 076107. [Google Scholar] [CrossRef] [PubMed]
  15. Algrain, M.C.; Ehlers, D.E. Novel Kalman filtering method for the suppression of gyroscope noise effects in pointing and tracking systems. Opt. Eng. 1995, 34, 3016–3031. [Google Scholar] [CrossRef]
  16. Lo, P.H.; Halatyn, P.; Geoca, E.; Archibald, J.B. Gyroscope Noise Reduction and Drift Compensation. U.S. Patent 5,795,988, 18 August 1998. [Google Scholar]
  17. Chong, S.; Rui, S.; Jie, L.; Xiaoming, Z.; Jun, T.; Yunbo, S.; Jun, L.; Huiliang, C. Temperature drift modeling of MEMS gyroscope based on genetic-Elman neural network. Mech. Syst. Signal Process. 2016, 72, 897–905. [Google Scholar] [CrossRef]
  18. Madgwick, S.O.; Harrison, A.J.; Vaidyanathan, R. Estimation of IMU and MARG orientation using a gradient descent algorithm. In Proceedings of the 2011 IEEE International Conference on Rehabilitation Robotics, Zurich, Switzerland, 29 June–1 July 2011; pp. 1–7. [Google Scholar]
  19. Carberry, J.; Hinchly, G.; Buckerfield, J.; Tayler, E.; Burton, T.; Madgwick, S.; Vaidyanathan, R. Parametric design of an active ankle foot orthosis with passive compliance. In Proceedings of the 2011 24th International Symposium on Computer-Based Medical Systems (CBMS), Bristol, UK, 27–30 June 2011; pp. 1–6. [Google Scholar]
  20. Susi, M.; Renaudin, V.; Lachapelle, G. Motion mode recognition and step detection algorithms for mobile phone users. Sensors 2013, 13, 1539–1562. [Google Scholar] [CrossRef] [PubMed]
  21. Weinberg, H. Using the ADXL202 in pedometer and personal navigation applications. Analog Devices AN-602 Appl. Note 2002, 2, 1–6. [Google Scholar]
  22. Shin, S.; Park, C.; Kim, J.; Hong, H.; Lee, J. Adaptive step length estimation algorithm using low-cost MEMS inertial sensors. In Proceedings of the Sensors Applications Symposium, San Diego, CA, USA, 6–8 February 2007; pp. 1–5. [Google Scholar]
  23. Renaudin, V.; Susi, M.; Lachapelle, G. Step length estimation using handheld inertial sensors. Sensors 2012, 12, 8507–8525. [Google Scholar] [CrossRef] [PubMed]
  24. Ristic, B.; Arulampalam, S.; Gordon, N. Beyond the Kalman Filter: Particle Filters for Tracking Applications; Artech House: Norwood, MA, USA, 2003. [Google Scholar]
  25. Xiao, Z.; Wen, H.; Markham, A.; Trigoni, N. Lightweight map matching for indoor localisation using conditional random fields. In Proceedings of the 13th International Symposium on Information Processing in Sensor Networks, Berlin, Germany, 15–17 April 2014; pp. 131–142. [Google Scholar]
  26. Shin, S.H.; Park, C.G.; Choi, S. New Map-Matching Algorithm Using Virtual Track for Pedestrian Dead Reckoning. ETRI J. 2010, 32, 891–900. [Google Scholar] [CrossRef]
  27. Arulampalam, M.S.; Maskell, S.; Gordon, N.; Clapp, T. A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking. IEEE Trans. Signal Process. 2002, 50, 174–188. [Google Scholar] [CrossRef] [Green Version]
  28. Davidson, P.; Collin, J.; Takala, J. Application of particle filters for indoor positioning using floor plans. In Proceedings of the Ubiquitous Positioning Indoor Navigation and Location Based Service (UPINLBS), Kirkkonummi, Finland, 14–15 October 2010; pp. 1–4. [Google Scholar]
  29. Klepal, M.; Beauregard, S.; Widyawan. A backtracking particle filter for fusing building plans with PDR displacement estimates. In Proceedings of the 2008 5th Workshop on Positioning, Navigation and Communication, Hannover, Germany, 27 March 2008; pp. 207–212. [Google Scholar]
  30. Bojja, J.; Kirkko-Jaakkola, M.; Collin, J.; Takala, J. Indoor localization methods using dead reckoning and 3D map matching. J. Signal Process. Syst. 2014, 76, 301–312. [Google Scholar] [CrossRef]
  31. Klepal, M.; Beauregard, S.; Widyawan. A novel backtracking particle filter for pattern matching indoor localization. In Proceedings of the First ACM International Workshop on Mobile Entity Localization and Tracking in GPS-Less Environments, Orlando, FL, USA, 30 September 2008; pp. 79–84. [Google Scholar]
Figure 1. A prototype system of the chest-mounted intertial measurement unit (IMU)-based pedestrian dead reckoning (PDR).
Figure 1. A prototype system of the chest-mounted intertial measurement unit (IMU)-based pedestrian dead reckoning (PDR).
Sensors 19 00420 g001
Figure 2. System overview.
Figure 2. System overview.
Sensors 19 00420 g002
Figure 3. Sensor and world frames.
Figure 3. Sensor and world frames.
Sensors 19 00420 g003
Figure 4. Step-detection process.
Figure 4. Step-detection process.
Sensors 19 00420 g004
Figure 5. Definition of one step.
Figure 5. Definition of one step.
Sensors 19 00420 g005
Figure 6. Data of one step.
Figure 6. Data of one step.
Sensors 19 00420 g006
Figure 7. Map components.
Figure 7. Map components.
Sensors 19 00420 g007
Figure 8. Definition of heading correction h i .
Figure 8. Definition of heading correction h i .
Sensors 19 00420 g008
Figure 9. Process of backtracking test. Blue arrows are recent steps.
Figure 9. Process of backtracking test. Blue arrows are recent steps.
Sensors 19 00420 g009
Figure 10. Pressure variations when going upstairs.
Figure 10. Pressure variations when going upstairs.
Sensors 19 00420 g010
Figure 11. Map-editing process.
Figure 11. Map-editing process.
Sensors 19 00420 g011
Figure 12. Calibrate IMU pose with initial heading.
Figure 12. Calibrate IMU pose with initial heading.
Sensors 19 00420 g012
Figure 13. Experiment 1: system evaluation with different configurations.
Figure 13. Experiment 1: system evaluation with different configurations.
Sensors 19 00420 g013
Figure 14. Experiment 1: errors at each keypoint.
Figure 14. Experiment 1: errors at each keypoint.
Sensors 19 00420 g014
Figure 15. Experiment 2: evaluation in IPIN 2018 Competition Track 2.
Figure 15. Experiment 2: evaluation in IPIN 2018 Competition Track 2.
Sensors 19 00420 g015
Figure 16. Experiment 2: errors at each keypoint.
Figure 16. Experiment 2: errors at each keypoint.
Sensors 19 00420 g016
Table 1. NGIMU specifications.
Table 1. NGIMU specifications.
RangeResolutionSampling Rate
Accelerometer ± 16 g490 μ g400 Hz
Gyroscope ± 2000 /s 0.06 /s400 Hz
Barometer 30 110 kPa 0.18 Pa25 Hz
Size 56 × 39 × 18 mm
Weight 46 g
Table 2. Definition of variables and functions in Algorithm 1.
Table 2. Definition of variables and functions in Algorithm 1.
P a r t i c l e N u m b e r M a x Maximal number of particles in our system.
P a r t i c l e N u m b e r Number of current existing particles.
nNumber of existing particles before generating new particles.
T r y T i m e M a x Maximal trying time on proposing a new particle.
T r y T i m e Counter of trying time on proposing a new particle.
P a r t i c l e Randomly selected particle.
R a d i u s Second input parameter of function P r o p o s e N e w P a r t i c l e .
S e l e c t R a n d o m P a r t i c l e ( ) Randomly select one particle from existing particles and return it.
P r o p o s e N e w P a r t i c l e ( _ p a r t i c l e , _ r a d i u s ) Propose a new particle around _ p a r t i c l e and return it.
B a c k t r a c k i n g T e s t ( _ p a r t i c l e ) Apply backtracking test to _ p a r t i c l e . Return P A S S if _ p a r t i c l e passed the test.
A p p e n d P a r t i c l e ( _ p a r t i c l e ) Append _ p a r t i c l e to existing particles.
Table 3. Experiment 1: error distribution.
Table 3. Experiment 1: error distribution.
Travelled distance: 432.22 m
Total number of keypoints: 25
Our system
Mean0.78 m
Median0.51 m
75th percent0.76 m
Standard deviation0.92 m
Table 4. Experiment 2: error distribution.
Table 4. Experiment 2: error distribution.
Travelled distance: 792.49 m
Keypoint number: 70
Our system
Mean5.2 m
Median3.6 m
75th percent5.7 m
Standard deviation5.0 m

Share and Cite

MDPI and ACS Style

Lu, C.; Uchiyama, H.; Thomas, D.; Shimada, A.; Taniguchi, R.-i. Indoor Positioning System Based on Chest-Mounted IMU. Sensors 2019, 19, 420. https://doi.org/10.3390/s19020420

AMA Style

Lu C, Uchiyama H, Thomas D, Shimada A, Taniguchi R-i. Indoor Positioning System Based on Chest-Mounted IMU. Sensors. 2019; 19(2):420. https://doi.org/10.3390/s19020420

Chicago/Turabian Style

Lu, Chuanhua, Hideaki Uchiyama, Diego Thomas, Atsushi Shimada, and Rin-ichiro Taniguchi. 2019. "Indoor Positioning System Based on Chest-Mounted IMU" Sensors 19, no. 2: 420. https://doi.org/10.3390/s19020420

APA Style

Lu, C., Uchiyama, H., Thomas, D., Shimada, A., & Taniguchi, R. -i. (2019). Indoor Positioning System Based on Chest-Mounted IMU. Sensors, 19(2), 420. https://doi.org/10.3390/s19020420

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