Journal of Sensors

Volume 2013 (2013), Article ID 197090, 10 pages

http://dx.doi.org/10.1155/2013/197090

## Low-Cost MEMS-Based Pedestrian Navigation Technique for GPS-Denied Areas

Department of Geomatics Engineering, The University of Calgary, 2500 University Drive N.W., Calgary, AB, Canada T2N 1N4

Received 10 April 2013; Accepted 22 July 2013

Academic Editor: Kai-Wei Chiang

Copyright © 2013 Abdelrahman Ali and Naser El-Sheimy. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The progress in the micro electro mechanical system (MEMS) sensors technology in size, cost, weight, and power consumption allows for new research opportunities in the navigation field. Today, most of smartphones, tablets, and other handheld devices are fully packed with the required sensors for any navigation system such as GPS, gyroscope, accelerometer, magnetometer, and pressure sensors. For seamless navigation, the sensors’ signal quality and the sensors availability are major challenges. Heading estimation is a fundamental challenge in the GPS-denied environments; therefore, targeting accurate attitude estimation is considered significant contribution to the overall navigation error. For that end, this research targets an improved pedestrian navigation by developing sensors fusion technique to exploit the gyroscope, magnetometer, and accelerometer data for device attitude estimation in the different environments based on quaternion mechanization. Results indicate that the improvement in the traveled distance and the heading estimations is capable of reducing the overall position error to be less than 15 m in the harsh environments.

#### 1. Introduction

Personal navigation requires technologies that are immune to signal obstructions and fading. One of the major challenges is obtaining a good heading solution in different environments and for different user positions without external absolute reference signals. Part of this challenge arises from the complexity and freedom of movement of a typical handheld user where the heading observability considerably degrades in low-speed walking, making this problem even more challenging. However, for short periods, the relative attitude and heading information is quite reliable. Self-contained systems requiring minimal infrastructure, for example, inertial measurement units (IMUs), stand as a viable option, since pedestrian navigation is not only focused on outdoor navigation but also on indoor navigation.

Nowadays, most of the smartphones are programmable and equipped with self-contained, low cost, small size, and power-efficient sensors, such as magnetometers, gyroscopes, and accelerometers. Hence, integrating IMUs navigation solution with a magnetometer-based heading can play an important role in pedestrian navigation in all environments. In the current state of the art in MEMS technology, the accuracy of gyroscopes is not good enough for deriving an absolute heading or relative heading over longer durations of time. However, for short periods, the relative attitude information is quite reliable. Magnetometers, on the other hand, provide absolute heading information once calibrated. However, they can easily be disturbed by ferrous objects nearby, making them unreliable for brief intervals. This calls for the investigation of possible sources of heading error in complementary sensors such as a gyroscope and a magnetometer and improving the accuracy of the result based on an improved Kalman filter design.

Much research towards the heading estimation for personal positioning applications has been conducted in the recent years. Some approaches use magnetometers exclusively for heading estimation [1] while others integrate it tightly with an IMU [2, 3]. One commercially available personal locator system based on this principle is the Dead Reckoning Module DRM-4000 made by Honeywell [4]. A quaternion-based method to integrate IMU with magnetometer is presented by [5]. Three body angular rates and four quaternion elements were used to express attitude and were selected as the states of the Kalman filter. The method needs to model the angular motion of the body. In [6], a linear system error model based on the Euler angles errors expressing the local frame errors is developed, and the corresponding system observation model is derived. The proposed method does not need to model the system angular motion and also avoids the nonlinear problem which is inherent in the customarily used methods. A similar technique is proposed by [7] where the angular rates were modeled to be a constant. A nonlinear derivative equation for the Euler angle integration kinematics is investigated in [8]. Work in [9, 10] presented an Euler angle error based method to integrate IMU with magnetometer data where three Euler angle errors and three gyroscope biases were used as states for the Kalman filter. The estimated states were used to correct the Euler angles and to compensate gyroscope drifts, respectively. The work at [11] presented a mathematical model for compass deviation by creating an a priori look-up table for heading corrections. A Kalman filtering approach was investigated by [12] to estimate the angular rotation from the input of a magnetometer compass and three gyroscopes. References [13, 14] presented a least squares technique with improvement which is used for the estimation of the compass deviation model. In addition, much research has been conducted to use the 3D magnetometer-based heading for personal navigation applications in the recent years [15].

The magnetometer cannot be used as standalone source for heading information in the harsh environments, especially indoor [16]. In addition, it is required to have knowledge about the preexisted magnetic anomalies resulted from some of the man-made infrastructure [17]. Using magnetic field measurements in heading estimation for indoor navigation also has some limitations as the magnetic field signal needs to be strong enough. Also, the mobile navigation device should be away from any source of disturbances to avoid any perturbation effect [18]. Besides that, the magnetic field during the indoor environment is not completely constant due to the presence of the electronic and electrical devices everywhere. To avoid the problem of magnetometer anomaly, arising out of ferrous materials in the vicinity of the magnetometers, a perturbation detection technique is required. In such scenario, the filter works only in the propagation mode without any update for the attitude. Also the gyroscope bias drifts with time and temperature can be compensated by magnetometers. In this paper, a method is presented to obtain seamless attitude information by integrating the heading outputs based on magnetometer, accelerometer, and gyroscope measurements using the Kalman filter (KF).

#### 2. Pedestrian Dead Reckoning (PDR)

Pedestrian dead reckoning is a relative means of positioning where the initial position and heading of the user are supposed to be known. The basic concept and components of the proposed PDR algorithm are shown in Figure 1. Generally, steps of the user are detected based on the accelerometer data. To get the travelled distance, the total number of steps is multiplied by the step length. With known heading and reference point, the user can be tracked by successive steps count.

Step detection is a basic step for any PDR algorithm. Usually, the accelerometer signal is used to detect the steps of the person. Once the step is detected, the total number of steps for a pedestrian can be counted. As a result, the total travelled distance can be estimated by multiplying the step length with the total number of steps. Using travelled distance and heading, user can be located during a typical trip. Step detection algorithm can be performed based on different kinds of sensors, that is, not only accelerometers but even gyroscopes and magnetometers. However, our step event detection scheme is based on using the accelerometer signal. The norm of the three accelerometers is used as in (1), where it is possible to clearly identify the steps by observing, for example, the signal over time: Steps are detected as peaks in the resulting norm, where the step is the highest local maximum in the norm acceleration between the current peak and the previous step peak.

#### 3. Sensors Performance

The used device in the test, Samsung Galaxy Nexus smartphone, is shown in Figure 2 with axes definition.

Besides other sensors, the device is equipped with triad magnetometer (M), triad gyroscope (G), and triad accelerometer (A). The manufactures and the ranges of the main used sensors are listed in Table 1.

##### 3.1. Gyroscope Drift

Gyroscopes are mainly used to determine device attitude. The output of such sensor is rotational rate, and performing a single integration on the gyroscopes outputs is necessary to obtain a relative change in angle. Due to the integration process which is very sensible to the systematic errors of the gyroscopes, the bias introduces a quadratic error in the velocity and a cubic error in the position [19]. Gyroscopes measurements can generally be described using where is the measured angular rate, is the true angular rate, is the gyroscope bias, is the linear scale factor matrix, is the nonorthogonality matrix, and is the sensor noise. With integration, the gyroscope bias will introduce an angle error in pitch or roll proportional to time; that is, ; this small angle will cause misalignment of the IMU. Therefore, when projecting the acceleration from for example the gravity vector , from the body frame to the local-level frame, the acceleration vector will be incorrectly projected due to this misalignment error. This will introduce an error in one of the horizontal acceleration; that is, . Consequently, this leads to an error in velocity and in position . To overcome the problem of error drift, a bias compensation of gyroscopes is required. When the device is in stationary mode, the deterministic bias of the gyroscope can be estimated by calculating the average gyroscope output during that time.

##### 3.2. Magnetometer Perturbation

During operation, especially inside buildings, a magnetometer is subject to many external disturbances such as large metal objects [20]. Other objects like steel structures and electromagnetic power lines can affect the solution of the magnetometer. These kinds of disturbance lead to unpredictable performance of the magnetometer which is a major drawback of using magnetic sensors. In the meanwhile, the magnetic field parameters such as strength, horizontal and vertical magnetic field, and change in the inclination angle can be checked to detect the perturbation in magnetometer measurement by comparing to the reference values which can be found at [21]. Figure 3 shows an example for harsh environment as the magnetometer is totally disturbed due to the steel constructions. The figure shows the total magnetic field in a perturbed area during a walking test compared to the reference value which is 570 mGauss for Calgary [21].

#### 4. Multisensors Heading Fusion Filter

The attitude of the device is commonly estimated using the inertial sensor. There are three main approaches for attitude representation: DCM, Euler angle, and quaternion. Among the three techniques quaternion algebra is the preferred. However, the estimated attitude from the gyroscope is very noisy leading to unbounded growth in the heading errors. An integration scheme for the gyroscope, accelerometer, and magnetometer data is proposed to estimate the device attitude and the gyroscope bias. The proposed scheme is a quaternion-based KF as shown in Figure 4.

In order to use the KF-based estimator for quaternion parameters and gyroscope biases estimation for a device which is carried by a pedestrian, the required model for the states and measurements and their respective system and measurement error models are presented in this section.

##### 4.1. Quaternion Mechanization

The relationship between the direct cosine matrix (DCM) and Euler angles is given in (3) with the sequence of azimuth, pitch, and roll () or () [22]: A quaternion is a four-dimensional vector which is defined based on a vector and a rotation angle. The vector is given as The DCM matrix in terms of quaternion vector components can be obtained from using where represents the DCM matrix in terms of the quaternion vector.

The matrix transforms from the body frame to the local level (navigation) frame. The roll, pitch, and azimuth values can be obtained by using the math function on the values of the propagated inside the sensors navigation equations

##### 4.2. Filter States

Basically, the target of the proposed filter is to estimate the device attitude based on the quaternion technique. Consequently, any improvement in the quaternion estimate leads to improving the estimated attitude values. The implementation of the KF is optimal for linear systems driven by additive white Gaussian noise (AWGN). The state model can be written in the following form: where is the state vector, is the state transition matrix, and represents the covariance matrix of the applied state model. The measurement system can be represented by a linear equation of the following form: where is the vector of measurement updates, is the design (observation) matrix that relates the measurements to the state vector, and is the measurement noise.

The nonlinear form of the system model in the absence of the known input can be written as where is now a nonlinear function describing the time evolution of the states. Consider a nominal trajectory, , related to the actual trajectory, , as where is a perturbation from nominal trajectory. performing a Taylor series expansion equation (10) about the nominal trajectory yields where is now the dynamic matrix for a system with state vector which consists of the perturbed states, . Thus, the main states to be estimated are the errors in the quaternion parameters given by: where is the error in the th quaternion parameter.

The quaternion parameters are primarily determined using the angular rates obtained from gyroscopes’ measurements. The deterministic errors associated with the gyroscope can be compensated using data from static interval at the beginning of the test while the stochastic errors in biases are given by where , , and are the gyroscope biases.

The complete state vector is defined as a 7-dimensional vector with the first four components being errors in the elements of the quaternion and the last three being the elements of the gyroscope biases.

##### 4.3. The State Transition Model

The angular rate is linked to the quaternion parameters as in the following: where is the attitude quaternion , , and represent angular rate measurements in the sensor frame obtained using the rate gyroscopes. Quaternion is used to represent attitude in the filter design because it does not have the singularity problem associated with Euler angles. The previous equation can be rewritten as The Taylor series expansion to first order is shown in the following: where is the state vector is the error in the state. Thus, the quaternion parameters error can be obtained as A general equation for the 1st order Gauss-Markov model is given as where is the reciprocal of the correlation time and is the variance of the gyroscope signal. The different parameters of the Gauss-Markov model can be determined as shown in Figure 5.

According to (19), the gyroscopes bias can be modeled as The complete state model can be written as the following: where Using the dynamic matrix in (21), the state transition matrix can be defined as

##### 4.4. Measurement Models

The magnetometer measurements along with the accelerometer data are used as main source of update. The measured magnetic field is tested for perturbation. Once the magnetic field is free of disturbances, the geomagnetic heading is estimated from the calibrated data. Also, the accelerometer measurements are experienced to be noisy as low cost MEMS sensors are used and usually measured at higher rate. Therefore, the average of the measured data over the step time is used to estimate the roll and pitch values. The roll, pitch, and heading estimates are used to calculate the quaternion parameters.

The only source of update is the quaternion parameters. Thus, the design matrix can be set as

##### 4.5. Modeling of Process and Measurement Noises

In order to complete the design of the KF, it is necessary to define the noise covariance matrices, the process noise covariance matrix and the measurement noise covariance matrix . These matrices reflect the confidence in the system model and the measurements, respectively. The covariance of is often called the process noise matrix, , and can be computed as The matrix is a 7-dimension square matrix which can be computed using as follows [23, 24]:

The measurement noise covariance matrix is also known as the covariance matrix for . The matrix represents the level of confidence placed in the accuracy of the measurements and is given by The matrix is a 4-dimension diagonal square matrix. The diagonal elements are the variances of the individual measurements, which can be determined experimentally using measurement data from the used sensors.

##### 4.6. Filter State Initialization

The state vector should be initialized at the beginning of the process. For the gyroscope bias states, all biases are initialized as zeros. The quaternion states can be initialized from the DCM matrix using the Euler angles. The mean of the accelerometer calibrated data during a stationary period can be used to estimate the initial roll and pitch using the following relationships [25]: where is the mean of the accelerometer data. During the same interval, the roll and pitch estimates are used for leveling the magnetometer data to be in the navigation frame. The calibrated magnetometer data is used to estimate the initial azimuth as The DCM is calculated using the initial Euler angles values as in (3). Then, the following relation between the quaternion and the DCM [22] is used to calculate the initial quaternion vector: The initial quaternion vector, , is calculated using the initial Euler angles values (, , ). Equation (30) is also used to get the updated quaternion parameters.

#### 5. Results and Discussion

In this section, the performance of the proposed attitude algorithm is assessed. The device is held in the texting mode. All tests start with around 30 seconds of stationary period to calibrate for the gyroscope while the magnetometer is calibrated by moving the device in the 3D space. The first scenario is conducted to assess the performance of the proposed technique in the indoor environments while the second scenario is performed in the downtown area.

##### 5.1. Environment Change Test

In this scenario, the test is started outside the building close to the Olympic Oval at the University of Calgary. The PDR solution is initialized using initial position from GPS and initial heading from the magnetometer. The device is held in stationary at the beginning of the test for about 40 s to calibrate for the gyroscope deterministic bias and to estimate the initial orientation, roll and pitch, of the device. The device is held in the compass, texting mode, and the user keeps going through the first floor with a long corridor inside the building. The trajectory is ended outside the McEwan Student Center in front of Taylor Family Digital Library (TDFL). The time for this trajectory was about 5.5 minutes taken in 594 steps with a total travelled distance of 461 m. Various reasons made this place a candidate for the test. (i)It is a popular and attractive place for students’ activities.(ii)It has a long corridor which is a challenging area for magnetometer.(iii)At the starting point, in front of the building, there is a huge metal structure which can affect the magnetometer performance.(iv)This place is full of students and simulates the normal walking scenario of a smartphone user. The derived heading form the magnetometer is totally affected at the beginning due to the presence of a huge metal object in the surrounding area as shown in Figure 6. The total magnetic field is distorted yielding incorrect heading estimation from magnetometer during the first 40 seconds. In addition, the calibration process is performed based on data taken outside the building, and so, once the user is moved to be inside the building, the distribution of the magnetic field is different from that outside as shown in Figure 6. Thus, magnetometer-based heading is still not accurate. However, after around 2 minutes the heading from the magnetometer starts to be the main source of update for the attitude filter. As a result, during the perturbation interval the attitude filter does not perform the update stage and keeps propagating the attitude of the device based on the gyroscopes’ measurements.

The magnified part in Figure 7 shows the magnetometer heading during a perturbation area. The figure shows that the heading is diverted and scattered when the magnetic field is distorted in contrast to the attitude filter heading during the same interval.

To evaluate the overall accuracy of the PDR algorithm, a reference trajectory plotted on the map of the test site is used. A reference trajectory for the actual direction of the test is plotted on the map of the first floor of the University of Calgary. As shown in Figure 8, the maximum error is less than 5 m during the test. Also, the trajectory is finished at the correct point with position drift of less than 4 m. The error in the distance is about 1.1% of the total traveled distance.

##### 5.2. Downtown Test

Downtown is an attractive place for tourists in addition to most people in the city. It has the major attractive sites such as shopping centers, administration offices, museums, restaurants, cafes, and theaters. So, it is important to have a good navigation system which is able to help the pedestrians to their destinations. Some conventional techniques such as GPS suffer from the multipath and signal destruction due to the high buildings. The test is conducted in downtown Calgary to assess the performance of the proposed algorithm in the harsh environments. The performance of the magnetometer is totally affected due to the distortion of the magnetic field in the downtown area. The high buildings, cars, and traffic signals add more complications for the heading estimation based on the magnetometer. The test is started at the intersection of the 7th street SW and the 8th avenue SW downtown Calgary as shown in Figure 9. The selected trajectory is a square starting in the east direction followed by north, west, and south directions. The length of the trajectory is 490 m taking about 6 minutes of walking.

Although the magnetometer gets good maneuvering at the start of the test the presence of a strong perturbation source affected the heading estimate. As shown in Figure 10, the estimated heading from magnetometer is totally unused for major part of the test due to the distortion in the magnetic field. It shows that the heading update is available only for two minutes during the interval of 140 s to 260 s.

Figure 11 shows the PDR solution using the attitude filter compared to the same solution based on the heading from the magnetometer stand-alone, gyroscope heading stand-alone, and GPS solution. As shown in the figure, there is no accurate stand-alone solution for the test trajectory. (i)The gyroscope stand-alone solution gives the accurate direction with position drift due to the uncompensated bias. (ii)The magnetometer stand-alone solution is diverted at many parts of the test; however, in certain parts it performs well and provides the correct heading. (iii)The GPS is considered the poorest solution as it does not provide any correct information at all during the test due to the satellites unavailability. However, the PDR algorithm based on the attitude filter, in corporation with the magnetometer anomaly detection technique, provides an acceptable solution.

The overall accuracy of the PDR algorithm is evaluated by comparing the PDR solution with a reference trajectory on Google map. As shown in Figure 11, the maximum error in position happens at the north direction side of the trajectory where no update is available yet, and it was around 10 m. Also, the trajectory is finished at the correct place with position drift of less than 4 m. The overall error in distance is about 1% of the total traveled distance.

#### 6. Conclusion

An enhanced attitude estimation technique based on two different attitude sensors is presented. The algorithm is based on the quaternion mechanization to estimate the device attitude. The filter is working as a complementary technique; the Earth’s magnetic field and the angle rate are integrated to estimate the device heading. The filter is propagated in the prediction mode using the gyroscope measurements while the magnetometer and accelerometer measurements are used for the update stage. The improvement in the attitude estimation leads to an improved PDR result. The PDR-based estimated trajectory results are compared to reference trajectories to show the accuracy of the proposed technique. The results show that the presented algorithm is able to provide the necessary navigation information accurately even in the harsh environments.

#### Acknowledgments

This work was supported in part by research funds from TECTERRA Commercialization and Research Centre, the Canada Research Chairs Program, and the Natural Science and Engineering Research Council of Canada (NSERC).

#### References

- S. Y. Cho, K. W. Lee, C. G. Park, and J. G. Lee, “A personal navigation system using low-cost MEMS/GPS/Fluxgate,” in
*Proceedings of the 59th Annual Meeting of the Institute of Navigation and CIGTF 22nd Guidance Test Symposium*, Albuquerque, NM, USA, 2003. - C. Aparicio,
*Implementation of a Quaternion-Based Kalman Filter for Human Body Motion Tracking Using MARG Sensors*, California Naval Postgraduate School, Monterey, Calif, USA, 2004. - X. Yun and E. R. Bachmann, “Design, implementation, and experimental results of a quaternion-based Kalman filter for human body motion tracking,”
*IEEE Transactions on Robotics*, vol. 22, no. 6, pp. 1216–1227, 2006. View at Publisher · View at Google Scholar · View at Scopus - Honeywell, “DRM 4000 Dead Reckoning Module,” 2009, http://www51.honeywell.com/aero/common/documents/myaerospacecatalog-documents/Missiles-Munitions/DRM4000.pdf.
- J. L. Marins, X. Yun, E. R. Bachmann, R. B. McGhee, and M. J. Zyda, “An extended Kalman filter for quaternion-based orientation estimation using MARG sensors,” in
*Proceedings of IEEE/RSJ International Conference on Intelligent Robots and Systems*, vol. 4, pp. 2003–2011, November 2001. View at Scopus - S. Han and J. Wang, “A novel method to integrate IMU and magnetometers in attitude and heading reference systems,”
*Journal of Navigation*, vol. 64, no. 4, pp. 727–738, 2011. View at Publisher · View at Google Scholar · View at Scopus - S. Emura and S. Tachi, “Compensation of time lag between actual and virtual spaces by multi-sensor integration,” in
*Proceedings of IEEE International Conference on Multisensor Fusion and Integration for Intelligent Systems (MFI '94)*, pp. 463–469, October 1994. View at Scopus - J. M. Cooke, M. J. Zyda, D. R. Pratt, and R. B. McGhee, “NPSNET: flight simulation dynamic modeling using quaternions,”
*Presence*, vol. 1, pp. 404–420, 1992. View at Google Scholar - E. Foxlin, “Inertial head-tracker sensor fusion by a complementary separate-bias Kalman filter,” in
*Proceedings of the IEEE Virtual Reality Annual International Symposium*, pp. 185–194, April 1996. View at Scopus - P. Setoodeh, A. Khayatian, and E. Farjah, “Attitude estimation by separate-bias Kalman filter-based data fusion,”
*Journal of Navigation*, vol. 57, no. 2, pp. 261–273, 2004. View at Publisher · View at Google Scholar · View at Scopus - S.-W. Liu, Z.-N. Zhang, and J. C. Hung, “High accuracy magnetic heading system composed of fluxgate magnetometers and a microcomputer,” in
*Proceedings of the IEEE National Aerospace and Electronics Conference (NAECON '89)*, pp. 148–152, May 1989. View at Scopus - B. Hoff and R. Azuma, “Autocalibration of an electronic compass in an outdoor augmented reality system,” in
*Proceedings of IEEE and ACM International Symposium on Augmented Reality*, pp. 159–164, Munich, Germany, October 2000. - D. Gebre-Egziabher, G. H. Elkaim, J. D. Powell, and B. W. Parkinson, “A non-linear, two step estimation algorithm for calibrating solid-state Strapdown magnetometers,” in
*Proceedings of the 8th International St. Petersburg Conference on Navigation Systems*, IEEE/AIAA, St. Petersburg, Russia, 2001. - C. L. Tsai,
*Motion-based differential carrier phase positioning: heading determination and its applications [Ph.D. dissertation]*, National Taiwan University of Science and Technology, 2006. - S. Kwanmuang, L. Ojeda, and J. Borenstein, “Magnetometer-enhanced personal locator for tunnels and GPS-denied outdoor environments,” in
*Sensors, and Command, Control, Communications, and Intelligence (C3I) Technologies for Homeland Security and Homeland Defense X*, Proceedings of SPIE, Orlando, Fla, USA, April 2011. View at Publisher · View at Google Scholar - L. Xue, W. Yuan, H. Chang, and C. Jiang, “MEMS-based multi-sensor integrated attitude estimation technology for MAV applications,” in
*Proceedings of the 4th IEEE International Conference on Nano/Micro Engineered and Molecular Systems (NEMS '09)*, pp. 1031–1035, January 2009. View at Publisher · View at Google Scholar · View at Scopus - W. F. Storms and J. F. Raquet, “Magnetic field aided vehicle tracking,” in
*Proceedings of the 22nd International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS '09)*, pp. 1767–1774, September 2009. View at Scopus - E. R. Bachmann, X. Yun, and C. W. Peterson, “An investigation of the effects of magnetic variations on inertial/magnetic orientation sensors,” in
*Proceedings of IEEE International Conference on Robotics and Automation (ICRA '04)*, pp. 1115–1122, May 2004. View at Scopus - N. El-Sheimy,
*Inertial Techniques and INS/DGPS Integration*, Lecture Notes ENGO 623, Department of Geomatics Engineering, the University of Calgary, Calgary, Canada, 2012. - A. S. Ali, S. Siddharth, Z. Syed, C. L. Goodall, and N. El-Sheimy, “An efficient and robust maneuvering mode to calibrate low cost magnetometer for improved heading estimation for pedestrian navigation,”
*Journal of Applied Geodesy*, vol. 7, no. 1, pp. 65–73, 2013. View at Publisher · View at Google Scholar - C. C. Finlay, S. Maus, C. D. Beggan et al., “International geomagnetic reference field: the eleventh generation,”
*Geophysical Journal International*, vol. 183, no. 3, pp. 1216–1230, 2010. View at Publisher · View at Google Scholar · View at Scopus - J. B. Kuipers,
*Quaternions and Rotation Sequences*, Princeton University Press, Princeton, NJ, USA, 1999. - R. G. Brown and P. Y. Hwang,
*Introduction to Random Signals and Applied Kalman Filtering*, Wiley, New York, NY, USA, 1992. - A. Gelb,
*Applied Optimal Estimation*, MIT Press, 1974. - H. J. Luinge and P. H. Veltink, “Measuring orientation of human body segments using miniature gyroscopes and accelerometers,”
*Medical and Biological Engineering and Computing*, vol. 43, no. 2, pp. 273–282, 2005. View at Publisher · View at Google Scholar · View at Scopus