In recent years, microelectromechanical system (MEMS) inertial sensors (3D accelerometers and 3D gyroscopes) have become widely available due to their small size and low cost. Inertial sensor measurements are obtained at high sampling rates and can be integrated to obtain position and orientation information. These estimates are accurate on a short time scale, but suffer from integration drift over longer time scales. To overcome this issue, inertial sensors are typically combined with additional sensors and models. In this tutorial we focus on the signal processing aspects of position and orientation estimation using inertial sensors.We discuss different modeling choices and a selected number of important algorithms. The algorithms include optimizationbased smoothing and filtering as well as computationally cheaper extended Kalman filter and complementary filter implementations. The quality of their estimates is illustrated using both experimental and simulated data.

Anomalies in the ambient magnetic field can be used as features in indoor positioning and navigation. By using Maxwell’s equations, we derive and present a Bayesian nonparametric probabilistic modeling approach for interpolation and extrapolation of the magnetic field. We model the magnetic field components jointly by imposing a Gaussian process (GP) prior to the latent scalar potential of the magnetic field. By rewriting the GP model in terms of a Hilbert space representation, we circumvent the computational pitfalls associated with GP modeling and provide a computationally efficient and physically justified modeling tool for the ambient magnetic field. The model allows for sequential updating of the estimate and time-dependent changes in the magnetic field. The model is shown to work well in practice in different applications. We demonstrate mapping of the magnetic field both with an inexpensive Raspberry Pi powered robot and on foot using a standard smartphone.

The use of inertial measurement units (IMUs) has gained popularity for the estimation of lower limb kinematics. However, implementations in clinical practice are still lacking. The aim of this review is twofold—to evaluate the methodological requirements for IMU-based joint kinematic estimation to be applicable in a clinical setting, and to suggest future research directions. Studies within the PubMed, Web Of Science and EMBASE databases were screened for eligibility, based on the following inclusion criteria: (1) studies must include a methodological description of how kinematic variables were obtained for the lower limb, (2) kinematic data must have been acquired by means of IMUs, (3) studies must have validated the implemented method against a golden standard reference system. Information on study characteristics, signal processing characteristics and study results was assessed and discussed. This review shows that methods for lower limb joint kinematics are inherently application dependent. Sensor restrictions are generally compensated with biomechanically inspired assumptions and prior information. Awareness of the possible adaptations in the IMU-based kinematic estimates by incorporating such prior information and assumptions is necessary, before drawing clinical decisions. Future research should focus on alternative validation methods, subject-specific IMU-based biomechanical joint models and disturbed movement patterns in real-world settings.

We present a computationally efficient algorithm for using variations in the ambient magnetic field to compensate for position drift in integrated odometry measurements (dead-reckoning estimates) through simultaneous localization and mapping (SLAM). When the magnetic field map is represented with a reduced-rank Gaussian process (GP) using Laplace basis functions defined in a cubical domain, analytic expressions of the gradient of the learned magnetic field become available. An existing approach for magnetic field SLAM with reduced-rank GP regression uses a Rao-Blackwellized particle filter (RBPF). For each incoming measurement, training of the magnetic field map using an RBPF has a computational complexity per time step of $O\left(N_p{N}_m^2\right)$, where $N_p$ is the number of particles, and $N_m$ is the number of basis functions used to approximate the Gaussian process. Contrary to the existing particle filter-based approach, we propose applying an extended Kalman filter based on the gradients of our learned magnetic field map for simultaneous localization and mapping. Our proposed algorithm only requires training a single map. It, therefore, has a computational complexity at each time step of $O\left(N_m^2\right)$. We demonstrate the workings of the extended Kalman filter for magnetic field SLAM on an open-source data set from a foot-mounted sensor and magnetic field measurements collected onboard a model ship in an indoor pool. We observe that the drift compensating abilities of our algorithm are comparable to what has previously been demonstrated for magnetic field SLAM with an RBPF.

The ability to capture joint kinematics in outside-laboratory environments is clinically relevant. In order to estimate kinematics, inertial measurement units can be attached to body segments and their absolute orientations can be estimated. However, the heading part of such orientation estimates is known to drift over time, resulting in drifting joint kinematics. This study proposes a novel joint kinematic estimation method that tightly incorporates the connection between adjacent segments within a sensor fusion algorithm, to obtain drift-free joint kinematics. Drift in the joint kinematics is eliminated solely by utilizing common information in the accelerometer and gyroscope measurements of sensors placed on connecting segments. Both an optimization-based smoothing and a filtering approach were implemented. Validity was assessed on a robotic manipulator under varying measurement durations and movement excitations. Standard deviations of the estimated relative sensor orientations were below 0.89° in an optimization-based smoothing implementation for all robot trials. The filtering implementation yielded similar results after convergence. The method is proven to be applicable in biomechanics, with a prolonged gait trial of 7 minutes on 11 healthy subjects. Three-dimensional knee joint angles were estimated, with mean RMS errors of 2.14°, 1.85°, 3.66° in an optimization-based smoothing implementation and mean RMS errors of 3.08°, 2.42°, 4.47° in a filtering implementation, with respect to a golden standard optical motion capture reference system.

Hyperbolic localization measures the time difference of arrivals (TDOAs) of signals to determine the location of a wireless source or receiver. Traditional methods depend on precise clock synchronization between nodes so that time measurements from independent devices can be meaningfully compared. Imperfect synchronization is often the dominant source of error. We propose two new message-based TDOA equations for hyperbolic localization which require no synchronization and meet or exceed state-of-the-art accuracy. Our approaches leverage anchor nodes that observe each other’s packet arrival times and a novel reformulation of the TDOA equation to reduce the effect of clock drift error. Closed-form equations are derived for computing TDOA in both self-localization and source-localization modes of operation along with bounds on maximum clock drift error. Three experiments are performed, including a clock drift simulation, a nonline-of-sight (NLOS) simulation, and an indoor validation experiment on custom ultra wideband (UWB) hardware all of which involved eight anchor nodes and one localizing node in a 128- $\text{m}^{3}$ capture volume. Our source-localization approach achieved unprecedented accuracy with lower cost equipment and trivial setup. Our self-localization matched state-of-the-art accuracy but with infinite scalability and high privacy. These results could enable economical and infinite density indoor navigation and dramatically reduce the economic cost and increase the accuracy of implementing industrial and commercial tracking applications.