On-Manifold Preintegration for Real-Time Visual--Inertial Odometry

Christian ForsterLuca CarloneFrank DellaertDavide Scaramuzza

article2015IEEE Transactions on robotics1,355 citations2017 King-Sun Fu Transactions on Robotics Best Paper Award

Develops an inertial measurement preintegration formulation on the rotation group with closed-form bias updates, allowing real-time visual-inertial odometry to achieve high accuracy within factor graph optimization frameworks.

Listen

Autonomous navigation and motion tracking in environments without GPS rely heavily on combining monocular cameras with inertial measurement units (IMUs), a technique known as visual-inertial odometry (VIO). However, existing systems face an ongoing trade-off between estimation accuracy and real-time computation: traditional filtering methods are fast but suffer from accumulated linearization errors, while full optimization methods are highly accurate but computationally overwhelmed by the high rate of IMU data and the growing trajectory history.

The main objective of the article is to demonstrate an efficient and accurate VIO framework that achieves real-time, maximum a posteriori motion estimation. It evaluates a mathematical formulation for preintegrating high-rate inertial measurements directly on the rotation manifold and combining them with landmark-free visual factors in an incremental optimization pipeline.

The researchers developed a theoretical model that combines multiple IMU measurements between keyframes into single relative motion constraints on the 3D rotation manifold, avoiding repeated integrations when sensor biases are updated. They integrated this preintegrated IMU model into a factor graph using a structureless vision approach—which analytically eliminates 3D landmark points—and solved the resulting optimization incrementally using the iSAM2 algorithm. The system was validated using Monte Carlo simulations over 50 runs alongside real-world evaluations on a 430-meter indoor trajectory with motion-capture ground truth and multi-floor outdoor datasets, benchmarking against established methods such as OKVIS, MSCKF, and Google Tango.

The findings confirm substantial performance and accuracy gains over existing approaches. First, in real-world indoor testing, the proposed method accumulated only 0.3 meters of average drift over 360 meters of traveled distance, outperforming both OKVIS and MSCKF, which each accumulated roughly 0.7 meters (a drift reduction of more than 50%). Second, the framework achieved steady real-time execution on standard hardware, taking roughly 10 milliseconds for backend state updates and 3 milliseconds for visual feature tracking. Third, outdoor trials showed superior end-to-end accuracy, recording 1.0 meter and 0.5 meter drift on respective paths compared to 2.2 meters and 1.4 meters for Google Tango. Finally, statistical consistency evaluations demonstrated that the estimator properly tracks sensor biases and avoids overconfidence along unobservable degrees of freedom, completely removing the mathematical singularities and coordinate-dependent errors inherent to older Euler-angle preintegration methods.

These results show that engineering teams do not need to sacrifice trajectory accuracy to achieve real-time tracking on computationally constrained mobile platforms, such as drones or augmented reality devices. By avoiding redundant recalculations and landmark state explosion, the approach lowers computing overhead and minimizes risk of catastrophic tracking drift in GPS-denied environments.

Decision-makers and engineering leads should consider adopting on-manifold IMU preintegration and structureless vision factors for future autonomous navigation stacks, leveraging the authors' open-source implementation within the GTSAM toolbox. Future deployment planning should evaluate integrating loop-closure detection to eliminate long-term residual drift and explore higher-order integration schemes if deploying with lower-frequency IMUs.

Confidence in the system's performance is high based on the consistency of simulated Monte Carlo tests and physical validation against ground truth. Key operational considerations include sensitivity to initial sensor calibration (such as camera-IMU extrinsic alignment and time synchronization) and the assumption of sufficient visual feature tracking from the front-end pipeline, which may degrade in textureless or highly dynamic environments.

  • Paper: ORB-SLAM: A Versatile and Accurate Monocular SLAM System, Raul Mur-Artal et al. (2015). Introduces the keyframe-based tracking and local optimization architecture that underpins modern real-time visual estimation frameworks.
  • Paper: MonoSLAM: Real-Time Single Camera SLAM, Andrew J. Davison et al. (2007). Establishes real-time probabilistic visual SLAM using recursive filtering, providing fundamental context for subsequent graph-based and optimization-driven VIO methods.
  • Paper: An Introduction to the Kalman Filter, Greg Welch et al. (1995). Presents foundational state estimation and covariance propagation principles essential for formulating inertial sensor noise models and kinematic state propagation.
  • Paper: Estimating uncertain spatial relationships in robotics, Randall Smith et al. (1986). Develops the foundational probabilistic representation of uncertain spatial coordinate transformations used across robotics state estimation.
Cover for On-Manifold Preintegration for Real-Time Visual--Inertial Odometry

Abstract

Current approaches for visual-inertial odometry (VIO) are able to attain highly accurate state estimation via nonlinear optimization. However, real-time optimization quickly becomes infeasible as the trajectory grows over time, this problem is further emphasized by the fact that inertial measurements come at high rate, hence leading to fast growth of the number of variables in the optimization. In this paper, we address this issue by preintegrating inertial measurements between selected keyframes into single relative motion constraints. Our first contribution is a \emph{preintegration theory} that properly addresses the manifold structure of the rotation group. We formally discuss the generative measurement model as well as the nature of the rotation noise and derive the expression for the \emph{maximum a posteriori} state estimator. Our theoretical development enables the computation of all necessary Jacobians for the optimization and a-posteriori bias correction in analytic form. The second contribution is to show that the preintegrated IMU model can be seamlessly integrated into a visual-inertial pipeline under the unifying framework of factor graphs. This enables the application of incremental-smoothing algorithms and the use of a \emph{structureless} model for visual measurements, which avoids optimizing over the 3D points, further accelerating the computation. We perform an extensive evaluation of our monocular \VIO pipeline on real and simulated datasets. The results confirm that our modelling effort leads to accurate state estimation in real-time, outperforming state-of-the-art approaches.

Table of Contents

  • I Introduction
  • II Related Work
  • II-A Filtering
  • II-B Fixed-lag Smoothing
  • II-C Full Smoothing
  • III Preliminaries
  • III-A Notions of Riemannian geometry
  • III-B Uncertainty Description in SO⁡(3)\mathrm{SO}(3)
  • III-C Gauss-Newton Method on Manifold
  • IV Maximum a Posteriori Visual-Inertial State Estimation
  • IV-A The State
  • IV-B The Measurements
  • IV-C Factor Graphs and MAP Estimation
  • V IMU Model and Motion Integration
  • VI IMU Preintegration on Manifold
  • VI-A Preintegrated IMU Measurements
  • VI-B Noise Propagation
  • VI-C Incorporating Bias Updates
  • VI-D Preintegrated IMU Factors
  • VI-E Bias Model
  • VII Structureless Vision Factors
  • VIII Experimental Analysis
  • VIII-A Simulation Experiments
  • VIII-A1 Pose Estimation Accuracy and Timing
  • VIII-A2 Consistency
  • VIII-A3 Bias Estimation Accuracy
  • VIII-A4 First-Order Bias Correction
  • VIII-A5 Advantages over the Euler-angle-based formulation
  • VIII-B Real Experiments
  • VIII-B1 Implementation
  • VIII-B2 Indoor Experiments
  • VIII-B3 Outdoor Experiments
  • IX Conclusion
  • IX-A Iterative Noise Propagation
  • IX-B Bias Correction via First-Order Updates
  • IX-C Jacobians of Residual Errors
  • IX-C1 Jacobians of 𝐫Δ​𝐩i​j\mathbf{r}_{\Delta\mathbf{p}_{ij}}
  • IX-C2 Jacobians of 𝐫Δ​𝐯i​j\mathbf{r}_{\Delta\mathbf{v}_{ij}}
  • IX-C3 Jacobians of 𝐫Δ​𝚁i​j\mathbf{r}_{\Delta\mathtt{R}_{ij}}
  • IX-D Structureless Vision Factors: Null Space Projection
  • IX-E Rotation Rate Integration Using Euler Angles
  • References

Knowls

  1. Knowl 1 — On-Manifold IMU Preintegration Model

    model/method

    In visual-inertial state estimation, an Inertial Measurement Unit (IMU) outputs high-rate angular velocity ω~(t)=ω(t)+bg(t)+ηg(t)\tilde{\omega}(t) = \omega(t) + b^g(t) + \eta^g(t) and linear acceleration a~(t)=R(t)T(a(t)−g)+ba(t)+ηa(t)\tilde{a}(t) = R(t)^T(a(t) - g) + b^a(t) + \eta^a(t), where ω(t)∈R3\omega(t) \in \mathbb{R}^3 is instantaneous angular velocity, a(t)∈R3a(t) \in \mathbb{R}^3 is sensor acceleration, g∈R3g \in \mathbb{R}^3 is the gravity vector in world coordinates, R(t)∈SO(3)R(t) \in SO(3) is the orientation of the sensor frame relative to the world frame, bg(t),ba(t)∈R3b^g(t), b^a(t) \in \mathbb{R}^3 are sensor biases, and ηg(t),ηa(t)\eta^g(t), \eta^a(t) are additive white Gaussian noise processes.

    To avoid reintegrating measurements when the state estimate at keyframe ii changes during nonlinear optimization, inertial measurements between consecutive keyframes ii and jj at discrete times k∈{i,…,j−1}k \in \{i, \dots, j-1\} with time step Δt\Delta t are integrated into relative motion increments that are independent of the initial pose (Ri,pi)∈SE(3)(R_i, p_i) \in SE(3) and velocity vi∈R3v_i \in \mathbb{R}^3:

    ΔRij≐RiTRj=∏k=ij−1Exp((ω~k−big−ηkgd)Δt)\Delta R_{ij} \doteq R_i^T R_j = \prod_{k=i}^{j-1} \text{Exp}\left( (\tilde{\omega}_k - b_i^g - \eta_k^{gd}) \Delta t \right)

    Δvij≐RiT(vj−vi−gΔtij)=∑k=ij−1ΔRik(a~k−bia−ηkad)Δt\Delta v_{ij} \doteq R_i^T (v_j - v_i - g \Delta t_{ij}) = \sum_{k=i}^{j-1} \Delta R_{ik} (\tilde{a}_k - b_i^a - \eta_k^{ad}) \Delta t

    Δpij≐RiT(pj−pi−viΔtij−12gΔtij2)=∑k=ij−1(ΔvikΔt+12ΔRik(a~k−bia−ηkad)Δt2)\Delta p_{ij} \doteq R_i^T \left( p_j - p_i - v_i \Delta t_{ij} - \frac{1}{2} g \Delta t_{ij}^2 \right) = \sum_{k=i}^{j-1} \left( \Delta v_{ik} \Delta t + \frac{1}{2} \Delta R_{ik} (\tilde{a}_k - b_i^a - \eta_k^{ad}) \Delta t^2 \right)

    where Δtij≐∑k=ij−1Δt\Delta t_{ij} \doteq \sum_{k=i}^{j-1} \Delta t, ΔRik≐RiTRk\Delta R_{ik} \doteq R_i^T R_k, Δvik≐RiT(vk−vi−gΔtik)\Delta v_{ik} \doteq R_i^T (v_k - v_i - g \Delta t_{ik}), and ηkgd,ηkad\eta_k^{gd}, \eta_k^{ad} are discrete-time noise terms. The operator Exp(ϕ)≐exp⁡(ϕ∧)\text{Exp}(\phi) \doteq \exp(\phi^\wedge) denotes the exponential map from R3\mathbb{R}^3 to SO(3)SO(3) using the skew-symmetric matrix hat operator (⋅)∧(\cdot)^\wedge.

  2. Knowl 2 — Preintegrated IMU Measurement Model and Residual Errors

    equation

    Given keyframe states xi=[Ri,pi,vi,bi]x_i = [R_i, p_i, v_i, b_i] and xj=[Rj,pj,vj,bj]x_j = [R_j, p_j, v_j, b_j], where (R,p)∈SE(3)(R, p) \in SE(3), v∈R3v \in \mathbb{R}^3, and b=[bg,ba]∈R6b = [b^g, b^a] \in \mathbb{R}^6, the preintegrated IMU measurements computed at a linearization bias estimate bˉi=[bˉig,bˉia]\bar{b}_i = [\bar{b}_i^g, \bar{b}_i^a] are defined as:

    ΔR~ij≐∏k=ij−1Exp((ω~k−bˉig)Δt)\Delta \tilde{R}_{ij} \doteq \prod_{k=i}^{j-1} \text{Exp}\left( (\tilde{\omega}_k - \bar{b}_i^g) \Delta t \right)

    Δv~ij≐∑k=ij−1ΔR~ik(a~k−bˉia)Δt\Delta \tilde{v}_{ij} \doteq \sum_{k=i}^{j-1} \Delta \tilde{R}_{ik} (\tilde{a}_k - \bar{b}_i^a) \Delta t

    Δp~ij≐∑k=ij−1(Δv~ikΔt+12ΔR~ik(a~k−bˉia)Δt2)\Delta \tilde{p}_{ij} \doteq \sum_{k=i}^{j-1} \left( \Delta \tilde{v}_{ik} \Delta t + \frac{1}{2} \Delta \tilde{R}_{ik} (\tilde{a}_k - \bar{b}_i^a) \Delta t^2 \right)

    The generative measurement model relating the preintegrated measurements to the states and random noise vector ηijΔ=[δϕijT,δvijT,δpijT]T∼N(09×1,Σij)\eta_{ij}^\Delta = [\delta \phi_{ij}^T, \delta v_{ij}^T, \delta p_{ij}^T]^T \sim \mathcal{N}(0_{9\times 1}, \Sigma_{ij}) is:

    ΔR~ij=RiTRjExp(δϕij)\Delta \tilde{R}_{ij} = R_i^T R_j \text{Exp}(\delta \phi_{ij})

    Δv~ij=RiT(vj−vi−gΔtij)+δvij\Delta \tilde{v}_{ij} = R_i^T (v_j - v_i - g \Delta t_{ij}) + \delta v_{ij}

    Δp~ij=RiT(pj−pi−viΔtij−12gΔtij2)+δpij\Delta \tilde{p}_{ij} = R_i^T \left( p_j - p_i - v_i \Delta t_{ij} - \frac{1}{2} g \Delta t_{ij}^2 \right) + \delta p_{ij}

    Incorporating first-order bias updates δbi=[δbig,δbia]\delta b_i = [\delta b_i^g, \delta b_i^a] with respect to bˉi\bar{b}_i, the preintegrated IMU factor residual error rIij=[rΔRijT,rΔvijT,rΔpijT]T∈R9r_{\mathcal{I}_{ij}} = [r_{\Delta R_{ij}}^T, r_{\Delta v_{ij}}^T, r_{\Delta p_{ij}}^T]^T \in \mathbb{R}^9 is given by:

    rΔRij=Log((ΔR~ij(bˉig)Exp(∂ΔRˉij∂bgδbig))TRiTRj)r_{\Delta R_{ij}} = \text{Log}\left( \left( \Delta \tilde{R}_{ij}(\bar{b}_i^g) \text{Exp}\left( \frac{\partial \Delta \bar{R}_{ij}}{\partial b^g} \delta b_i^g \right) \right)^T R_i^T R_j \right)

    rΔvij=RiT(vj−vi−gΔtij)−[Δv~ij(bˉi)+∂Δvˉij∂bgδbig+∂Δvˉij∂baδbia]r_{\Delta v_{ij}} = R_i^T (v_j - v_i - g \Delta t_{ij}) - \left[ \Delta \tilde{v}_{ij}(\bar{b}_i) + \frac{\partial \Delta \bar{v}_{ij}}{\partial b^g} \delta b_i^g + \frac{\partial \Delta \bar{v}_{ij}}{\partial b^a} \delta b_i^a \right]

    rΔpij=RiT(pj−pi−viΔtij−12gΔtij2)−[Δp~ij(bˉi)+∂Δpˉij∂bgδbig+∂Δpˉij∂baδbia]r_{\Delta p_{ij}} = R_i^T \left( p_j - p_i - v_i \Delta t_{ij} - \frac{1}{2} g \Delta t_{ij}^2 \right) - \left[ \Delta \tilde{p}_{ij}(\bar{b}_i) + \frac{\partial \Delta \bar{p}_{ij}}{\partial b^g} \delta b_i^g + \frac{\partial \Delta \bar{p}_{ij}}{\partial b^a} \delta b_i^a \right]

    where Log(⋅)≐(log⁡(⋅))∨\text{Log}(\cdot) \doteq (\log(\cdot))^\vee maps an element in SO(3)SO(3) to its tangent vector in R3\mathbb{R}^3.

  3. Knowl 3 — Iterative Covariance Propagation for Preintegrated IMU Measurements

    equation

    The preintegrated measurement noise vector ηijΔ=[δϕijT,δvijT,δpijT]T∈R9\eta_{ij}^\Delta = [\delta \phi_{ij}^T, \delta v_{ij}^T, \delta p_{ij}^T]^T \in \mathbb{R}^9 is propagated in discrete time from the raw discrete IMU noise vector ηkd=[(ηkgd)T,(ηkad)T]T∼N(06×1,Ση)\eta_k^d = [(\eta_k^{gd})^T, (\eta_k^{ad})^T]^T \sim \mathcal{N}(0_{6\times 1}, \Sigma_\eta) with covariance Ση=1Δtdiag(Cov(ηg),Cov(ηa))\Sigma_\eta = \frac{1}{\Delta t} \text{diag}(\text{Cov}(\eta^g), \text{Cov}(\eta^a)).

    The linear recurrence relation for the noise vector is:

    ηijΔ=Aj−1ηij−1Δ+Bj−1ηj−1d\eta_{ij}^\Delta = A_{j-1} \eta_{ij-1}^\Delta + B_{j-1} \eta_{j-1}^d

    with transition matrices:

    Aj−1=[ΔR~j−1jT03×303×3−ΔR~ij−1(a~j−1−bˉia)∧ΔtI3×303×3−12ΔR~ij−1(a~j−1−bˉia)∧Δt2ΔtI3×3I3×3]A_{j-1} = \begin{bmatrix} \Delta \tilde{R}_{j-1 j}^T & 0_{3\times 3} & 0_{3\times 3} \\ -\Delta \tilde{R}_{ij-1} (\tilde{a}_{j-1} - \bar{b}_i^a)^\wedge \Delta t & I_{3\times 3} & 0_{3\times 3} \\ -\frac{1}{2} \Delta \tilde{R}_{ij-1} (\tilde{a}_{j-1} - \bar{b}_i^a)^\wedge \Delta t^2 & \Delta t I_{3\times 3} & I_{3\times 3} \end{bmatrix}

    Bj−1=[Jrj−1Δt03×303×3ΔR~ij−1Δt03×312ΔR~ij−1Δt2]B_{j-1} = \begin{bmatrix} J_r^{j-1} \Delta t & 0_{3\times 3} \\ 0_{3\times 3} & \Delta \tilde{R}_{ij-1} \Delta t \\ 0_{3\times 3} & \frac{1}{2} \Delta \tilde{R}_{ij-1} \Delta t^2 \end{bmatrix}

    where Jrj−1≐Jr((ω~j−1−bˉig)Δt)J_r^{j-1} \doteq J_r((\tilde{\omega}_{j-1} - \bar{b}_i^g)\Delta t) is the right Jacobian of SO(3)SO(3) defined by:

    Jr(ϕ)=I−1−cos⁡∥ϕ∥∥ϕ∥2ϕ∧+∥ϕ∥−sin⁡∥ϕ∥∥ϕ∥3(ϕ∧)2J_r(\phi) = I - \frac{1 - \cos \|\phi\|}{\|\phi\|^2} \phi^\wedge + \frac{\|\phi\| - \sin \|\phi\|}{\|\phi\|^3} (\phi^\wedge)^2

    The preintegrated measurement covariance Σij∈R9×9\Sigma_{ij} \in \mathbb{R}^{9\times 9} is computed recursively starting from Σii=09×9\Sigma_{ii} = 0_{9\times 9} via:

    Σij=Aj−1Σij−1Aj−1T+Bj−1ΣηBj−1T\Sigma_{ij} = A_{j-1} \Sigma_{ij-1} A_{j-1}^T + B_{j-1} \Sigma_\eta B_{j-1}^T

  4. Knowl 4 — First-Order A-Posteriori Bias Correction for Preintegrated Measurements

    model/method

    When optimization updates the bias estimate from bˉi=[bˉig,bˉia]\bar{b}_i = [\bar{b}_i^g, \bar{b}_i^a] to bˉi+δbi\bar{b}_i + \delta b_i, repeating numerical integration across all intermediate IMU measurements is avoided by applying a first-order Taylor expansion on SO(3)SO(3) and R3\mathbb{R}^3:

    ΔR~ij(big)≈ΔR~ij(bˉig)Exp(∂ΔRˉij∂bgδbig)\Delta \tilde{R}_{ij}(b_i^g) \approx \Delta \tilde{R}_{ij}(\bar{b}_i^g) \text{Exp}\left( \frac{\partial \Delta \bar{R}_{ij}}{\partial b^g} \delta b_i^g \right)

    Δv~ij(big,bia)≈Δv~ij(bˉig,bˉia)+∂Δvˉij∂bgδbig+∂Δvˉij∂baδbia\Delta \tilde{v}_{ij}(b_i^g, b_i^a) \approx \Delta \tilde{v}_{ij}(\bar{b}_i^g, \bar{b}_i^a) + \frac{\partial \Delta \bar{v}_{ij}}{\partial b^g} \delta b_i^g + \frac{\partial \Delta \bar{v}_{ij}}{\partial b^a} \delta b_i^a

    Δp~ij(big,bia)≈Δp~ij(bˉig,bˉia)+∂Δpˉij∂bgδbig+∂Δpˉij∂baδbia\Delta \tilde{p}_{ij}(b_i^g, b_i^a) \approx \Delta \tilde{p}_{ij}(\bar{b}_i^g, \bar{b}_i^a) + \frac{\partial \Delta \bar{p}_{ij}}{\partial b^g} \delta b_i^g + \frac{\partial \Delta \bar{p}_{ij}}{\partial b^a} \delta b_i^a

    The sensitivity Jacobians are computed incrementally during the preintegration step:

    ∂ΔRˉij∂bg=−∑k=ij−1ΔR~k+1j(bˉi)TJrkΔt\frac{\partial \Delta \bar{R}_{ij}}{\partial b^g} = - \sum_{k=i}^{j-1} \Delta \tilde{R}_{k+1 j}(\bar{b}_i)^T J_r^k \Delta t

    ∂Δvˉij∂ba=−∑k=ij−1ΔRˉikΔt\frac{\partial \Delta \bar{v}_{ij}}{\partial b^a} = - \sum_{k=i}^{j-1} \Delta \bar{R}_{ik} \Delta t

    ∂Δvˉij∂bg=−∑k=ij−1ΔRˉik(a~k−bˉia)∧∂ΔRˉik∂bgΔt\frac{\partial \Delta \bar{v}_{ij}}{\partial b^g} = - \sum_{k=i}^{j-1} \Delta \bar{R}_{ik} (\tilde{a}_k - \bar{b}_i^a)^\wedge \frac{\partial \Delta \bar{R}_{ik}}{\partial b^g} \Delta t

    ∂Δpˉij∂ba=∑k=ij−1(∂Δvˉik∂baΔt−12ΔRˉikΔt2)\frac{\partial \Delta \bar{p}_{ij}}{\partial b^a} = \sum_{k=i}^{j-1} \left( \frac{\partial \Delta \bar{v}_{ik}}{\partial b^a} \Delta t - \frac{1}{2} \Delta \bar{R}_{ik} \Delta t^2 \right)

    ∂Δpˉij∂bg=∑k=ij−1(∂Δvˉik∂bgΔt−12ΔRˉik(a~k−bˉia)∧∂ΔRˉik∂bgΔt2)\frac{\partial \Delta \bar{p}_{ij}}{\partial b^g} = \sum_{k=i}^{j-1} \left( \frac{\partial \Delta \bar{v}_{ik}}{\partial b^g} \Delta t - \frac{1}{2} \Delta \bar{R}_{ik} (\tilde{a}_k - \bar{b}_i^a)^\wedge \frac{\partial \Delta \bar{R}_{ik}}{\partial b^g} \Delta t^2 \right)

    where Jrk=Jr((ω~k−bˉig)Δt)J_r^k = J_r((\tilde{\omega}_k - \bar{b}_i^g)\Delta t).

  5. Knowl 5 — Analytic Jacobians of Preintegrated IMU Residual Errors

    theoretical result

    Gauss-Newton optimization on manifolds lifts the residual errors rIij=[rΔRijT,rΔvijT,rΔpijT]Tr_{\mathcal{I}_{ij}} = [r_{\Delta R_{ij}}^T, r_{\Delta v_{ij}}^T, r_{\Delta p_{ij}}^T]^T using local retractions: Ri←RiExp(δϕi),pi←pi+Riδpi,vi←vi+δvi,δbi←δbi+δb~iR_i \leftarrow R_i \text{Exp}(\delta \phi_i), \quad p_i \leftarrow p_i + R_i \delta p_i, \quad v_i \leftarrow v_i + \delta v_i, \quad \delta b_i \leftarrow \delta b_i + \tilde{\delta b}_i Rj←RjExp(δϕj),pj←pj+Rjδpj,vj←vj+δvjR_j \leftarrow R_j \text{Exp}(\delta \phi_j), \quad p_j \leftarrow p_j + R_j \delta p_j, \quad v_j \leftarrow v_j + \delta v_j

    The analytic Jacobians of the preintegrated residuals with respect to these local perturbations are:

    For rotation residual rΔRijr_{\Delta R_{ij}}: ∂rΔRij∂δϕi=−Jr−1(rΔR(Ri))RjTRi,∂rΔRij∂δϕj=Jr−1(rΔR(Rj))\frac{\partial r_{\Delta R_{ij}}}{\partial \delta \phi_i} = -J_r^{-1}(r_{\Delta R}(R_i)) R_j^T R_i, \quad \frac{\partial r_{\Delta R_{ij}}}{\partial \delta \phi_j} = J_r^{-1}(r_{\Delta R}(R_j)) ∂rΔRij∂δb~ig=−Jr−1(rΔRij(δbig))Exp(rΔRij(δbig))TJr(∂ΔRˉij∂bgδbig)∂ΔRˉij∂bg\frac{\partial r_{\Delta R_{ij}}}{\partial \tilde{\delta b}_i^g} = -J_r^{-1}(r_{\Delta R_{ij}}(\delta b_i^g)) \text{Exp}(r_{\Delta R_{ij}}(\delta b_i^g))^T J_r\left( \frac{\partial \Delta \bar{R}_{ij}}{\partial b^g} \delta b_i^g \right) \frac{\partial \Delta \bar{R}_{ij}}{\partial b^g} ∂rΔRij∂δpi=∂rΔRij∂δvi=∂rΔRij∂δpj=∂rΔRij∂δvj=∂rΔRij∂δb~ia=03×3\frac{\partial r_{\Delta R_{ij}}}{\partial \delta p_i} = \frac{\partial r_{\Delta R_{ij}}}{\partial \delta v_i} = \frac{\partial r_{\Delta R_{ij}}}{\partial \delta p_j} = \frac{\partial r_{\Delta R_{ij}}}{\partial \delta v_j} = \frac{\partial r_{\Delta R_{ij}}}{\partial \tilde{\delta b}_i^a} = 0_{3\times 3}

    For velocity residual rΔvijr_{\Delta v_{ij}}: ∂rΔvij∂δϕi=(RiT(vj−vi−gΔtij))∧,∂rΔvij∂δvi=−RiT,∂rΔvij∂δvj=RiT\frac{\partial r_{\Delta v_{ij}}}{\partial \delta \phi_i} = \left( R_i^T (v_j - v_i - g \Delta t_{ij}) \right)^\wedge, \quad \frac{\partial r_{\Delta v_{ij}}}{\partial \delta v_i} = -R_i^T, \quad \frac{\partial r_{\Delta v_{ij}}}{\partial \delta v_j} = R_i^T ∂rΔvij∂δb~ig=−∂Δvˉij∂bg,∂rΔvij∂δb~ia=−∂Δvˉij∂ba\frac{\partial r_{\Delta v_{ij}}}{\partial \tilde{\delta b}_i^g} = -\frac{\partial \Delta \bar{v}_{ij}}{\partial b^g}, \quad \frac{\partial r_{\Delta v_{ij}}}{\partial \tilde{\delta b}_i^a} = -\frac{\partial \Delta \bar{v}_{ij}}{\partial b^a} ∂rΔvij∂δϕj=∂rΔvij∂δpi=∂rΔvij∂δpj=03×3\frac{\partial r_{\Delta v_{ij}}}{\partial \delta \phi_j} = \frac{\partial r_{\Delta v_{ij}}}{\partial \delta p_i} = \frac{\partial r_{\Delta v_{ij}}}{\partial \delta p_j} = 0_{3\times 3}

    For position residual rΔpijr_{\Delta p_{ij}}: ∂rΔpij∂δϕi=(RiT(pj−pi−viΔtij−12gΔtij2))∧,∂rΔpij∂δpi=−I3×3,∂rΔpij∂δpj=RiTRj\frac{\partial r_{\Delta p_{ij}}}{\partial \delta \phi_i} = \left( R_i^T \left( p_j - p_i - v_i \Delta t_{ij} - \frac{1}{2} g \Delta t_{ij}^2 \right) \right)^\wedge, \quad \frac{\partial r_{\Delta p_{ij}}}{\partial \delta p_i} = -I_{3\times 3}, \quad \frac{\partial r_{\Delta p_{ij}}}{\partial \delta p_j} = R_i^T R_j ∂rΔpij∂δvi=−RiTΔtij,∂rΔpij∂δb~ig=−∂Δpˉij∂bg,∂rΔpij∂δb~ia=−∂Δpˉij∂ba\frac{\partial r_{\Delta p_{ij}}}{\partial \delta v_i} = -R_i^T \Delta t_{ij}, \quad \frac{\partial r_{\Delta p_{ij}}}{\partial \tilde{\delta b}_i^g} = -\frac{\partial \Delta \bar{p}_{ij}}{\partial b^g}, \quad \frac{\partial r_{\Delta p_{ij}}}{\partial \tilde{\delta b}_i^a} = -\frac{\partial \Delta \bar{p}_{ij}}{\partial b^a} ∂rΔpij∂δϕj=∂rΔpij∂δvj=03×3\frac{\partial r_{\Delta p_{ij}}}{\partial \delta \phi_j} = \frac{\partial r_{\Delta p_{ij}}}{\partial \delta v_j} = 0_{3\times 3}

    where the inverse right Jacobian is: Jr−1(ϕ)=I+12ϕ∧+(1∥ϕ∥2−1+cos⁡∥ϕ∥2∥ϕ∥sin⁡∥ϕ∥)(ϕ∧)2J_r^{-1}(\phi) = I + \frac{1}{2} \phi^\wedge + \left( \frac{1}{\|\phi\|^2} - \frac{1 + \cos \|\phi\|}{2 \|\phi\| \sin \|\phi\|} \right) (\phi^\wedge)^2

  6. Knowl 6 — Structureless Vision Factor via Null Space Projection

    model/method

    To eliminate the need to include 3D landmark positions ρl∈R3\rho_l \in \mathbb{R}^3 (l=1,…,Ll = 1, \dots, L) in the state vector, visual reprojection errors are marginalized out analytically at each Gauss-Newton iteration.

    For a landmark ll observed in a set of nln_l keyframes X(l)\mathcal{X}(l) with image measurements zil∈R2z_{il} \in \mathbb{R}^2, the lifted and linearized reprojection error is written in stacked matrix form as:

    ∥FlδTX(l)+Elδρl−bl∥2\| F_l \delta T_{\mathcal{X}(l)} + E_l \delta \rho_l - b_l \|^2

    where δTX(l)\delta T_{\mathcal{X}(l)} stacks the pose perturbations [δϕiT,δpiT]T∈R6[\delta \phi_i^T, \delta p_i^T]^T \in \mathbb{R}^6 for i∈X(l)i \in \mathcal{X}(l), El∈R2nl×3E_l \in \mathbb{R}^{2n_l \times 3} is the stacked landmark Jacobian, Fl∈R2nl×6nlF_l \in \mathbb{R}^{2n_l \times 6n_l} is the stacked pose Jacobian, and bl∈R2nlb_l \in \mathbb{R}^{2n_l} is the residual error vector normalized by the measurement noise standard deviation.

    Minimizing this quadratic term with respect to the landmark position correction yields:

    δρl∗=−(ElTEl)−1ElT(FlδTX(l)−bl)\delta \rho_l^* = -(E_l^T E_l)^{-1} E_l^T (F_l \delta T_{\mathcal{X}(l)} - b_l)

    Substituting δρl∗\delta \rho_l^* back eliminates the landmark variable, producing the structureless vision factor:

    ∥(El⊥)T(FlδTX(l)−bl)∥2\| (E_l^\perp)^T (F_l \delta T_{\mathcal{X}(l)} - b_l) \|^2

    where El⊥∈R2nl×(2nl−3)E_l^\perp \in \mathbb{R}^{2n_l \times (2n_l - 3)} is an orthonormal basis for the null space of ElE_l computed via Singular Value Decomposition (SVD), satisfying (El⊥)TEl⊥=I(E_l^\perp)^T E_l^\perp = I and El⊥(El⊥)T=I−El(ElTEl)−1ElTE_l^\perp (E_l^\perp)^T = I - E_l(E_l^T E_l)^{-1} E_l^T. This formulation avoids matrix inversions and reduces the optimization problem to only camera poses.

  7. Knowl 7 — Maximum A Posteriori Visual-Inertial Odometry on Factor Graphs

    model/method

    The visual-inertial odometry problem is formulated as Maximum A Posteriori (MAP) estimation over the trajectory states Xk={xi}i∈Kk\mathcal{X}_k = \{x_i\}_{i \in \mathcal{K}_k} at keyframe times Kk\mathcal{K}_k, where each state is xi=[Ri,pi,vi,bi]x_i = [R_i, p_i, v_i, b_i].

    The full negative log-posterior cost function optimized by the back-end is:

    Xk∗=arg⁡min⁡Xk(∥r0∥Σ02+∑(i,j)∈Kk∥rIij∥Σij2+∑(i,j)∈Kk∥rbij∥Σbij2+∑l=1L∥(El⊥)T(FlδTX(l)−bl)∥2)\mathcal{X}_k^* = \arg\min_{\mathcal{X}_k} \left( \| r_0 \|_{\Sigma_0}^2 + \sum_{(i,j) \in \mathcal{K}_k} \| r_{\mathcal{I}_{ij}} \|_{\Sigma_{ij}}^2 + \sum_{(i,j) \in \mathcal{K}_k} \| r_{b_{ij}} \|_{\Sigma_{b_{ij}}}^2 + \sum_{l=1}^L \| (E_l^\perp)^T (F_l \delta T_{\mathcal{X}(l)} - b_l) \|^2 \right)

    where:

    1. ∥r0∥Σ02\| r_0 \|_{\Sigma_0}^2 is the prior factor on the initial state x0x_0.
    2. ∥rIij∥Σij2\| r_{\mathcal{I}_{ij}} \|_{\Sigma_{ij}}^2 is the preintegrated IMU factor between keyframes ii and jj weighted by covariance Σij\Sigma_{ij}.
    3. ∥rbij∥Σbij2=∥bjg−big∥Σbgd2+∥bja−bia∥Σbad2\| r_{b_{ij}} \|_{\Sigma_{b_{ij}}}^2 = \| b_j^g - b_i^g \|_{\Sigma_{bgd}}^2 + \| b_j^a - b_i^a \|_{\Sigma_{bad}}^2 represents the random-walk IMU bias drift factors between keyframes, with covariances Σbgd=ΔtijCov(ηbg)\Sigma_{bgd} = \Delta t_{ij} \text{Cov}(\eta^{bg}) and Σbad=ΔtijCov(ηba)\Sigma_{bad} = \Delta t_{ij} \text{Cov}(\eta^{ba}) derived from continuous white noise densities ηbg,ηba\eta^{bg}, \eta^{ba}.
    4. The final term represents structureless vision factors constraining keyframes observing landmark ll.

    Inference is executed incrementally using the iSAM2 algorithm over the Bayes tree, which updates only the subset of variables affected by new measurements.

  8. Knowl 8 — SO(3) Manifold Preintegration Invariance and Singularity Avoidance over Euler Angles

    theoretical result

    Preintegration directly on the SO(3)SO(3) manifold overcomes three structural flaws inherent in Euler-angle-based IMU preintegration:

    1. Integration Exactness: Integrating angular rates using the matrix exponential map on SO(3)SO(3) is exact for constant angular velocity across the interval Δt\Delta t. In contrast, Euler angle integration θk+1=θk+[E′(θk)]−1(ω~k−ηkg)Δt\theta_{k+1} = \theta_k + [E'(\theta_k)]^{-1}(\tilde{\omega}_k - \eta_k^g)\Delta t (where E′(θk)E'(\theta_k) is the conjugate Euler angle rate matrix) is only a first-order approximation whose error accumulates rapidly with large sampling periods Δt\Delta t or high angular rates ω~\tilde{\omega}.

    2. Statistical Fairness and Coordinate Invariance: The negative log-likelihood on SO(3)SO(3), L(R)=12∥Log(R~−1R)∥Σ2\mathcal{L}(R) = \frac{1}{2}\| \text{Log}(\tilde{R}^{-1} R) \|_\Sigma^2, is left-invariant under rigid body coordinate transformations (fair metric). The Euler angle likelihood L(θ)=12∥θ~−θ∥Σ2\mathcal{L}(\theta) = \frac{1}{2}\| \tilde{\theta} - \theta \|_\Sigma^2 is not invariant, causing state estimation results to depend on the choice of the reference world frame.

    3. Singularity Avoidance: Euler angle parametrizations (such as ZYXZYX) suffer from gimbal lock singularities at pitch angles θ=±π/2+nπ\theta = \pm \pi/2 + n\pi (n∈Zn \in \mathbb{Z}). In the vicinity of these points, Euler angle noise propagation severely degrades, producing diverging Kullback-Leibler (KL) divergence against the ground-truth covariance. SO(3)SO(3) preintegration avoids representation singularities entirely.

  9. Knowl 9 — Real-Time Benchmark Performance and Odometry Accuracy on Indoor Benchmark

    empirical result

    The proposed monocular VIO pipeline (combining the SVO front-end, on-manifold preintegrated IMU factors, structureless vision factors, and iSAM2 back-end) was evaluated on a 430 m indoor trajectory recorded with a forward-looking VI-Sensor (20 Hz camera, 800 Hz ADIS16448 IMU) and ground truth from a Vicon motion-capture system.

    1. Accuracy vs. State-of-the-Art: Over trajectory segments of lengths {10,40,90,160,250,360}\{10, 40, 90, 160, 250, 360\} m, the proposed system achieved an average translation drift of 0.3 m over 360 m traveled distance, compared to 0.7 m average drift for both OKVIS (keyframe-based sliding-window optimization) and MSCKF (filtering). The approach exhibited significantly lower yaw drift while pitch and roll errors remained low and constant across all methods due to gravity observability.

    2. Outdoor Trajectory Accuracy: On a 300 m outdoor loop around a building, the proposed system accumulated 1.0 m end-to-end drift, compared to 2.2 m drift for Google Tango Peanut (mapper v3.15). On a 160 m multi-floor trajectory across three building floors, the proposed system accumulated 0.5 m drift compared to 1.4 m for Google Tango.

    3. Execution Runtime: The SVO front-end required ~3 ms per frame for feature tracking on an Intel i7 2.4 GHz laptop. The back-end iSAM2 incremental update optimized keyframes with 10 iterations in ~10 ms per keyframe on average.

  10. Knowl 10 — Estimation Consistency and Observability-Preserving Covariance Bounds

    empirical result

    The statistical consistency of the on-manifold VIO estimator was verified using 50 Monte Carlo simulation runs over a 120 m 3D sinusoidal circular trajectory (corrupted with continuous gyro noise σg=0.0007 rad/(sHz)\sigma^g = 0.0007\text{ rad}/(\text{s}\sqrt{\text{Hz}}), accel noise σa=0.019 m/(s2Hz)\sigma^a = 0.019\text{ m}/(\text{s}^2\sqrt{\text{Hz}}), gyro bias noise σbg=0.0004 rad/(s2Hz)\sigma^{bg} = 0.0004\text{ rad}/(\text{s}^2\sqrt{\text{Hz}}), accel bias noise σba=0.012 m/(s3Hz)\sigma^{ba} = 0.012\text{ m}/(\text{s}^3\sqrt{\text{Hz}}), and σpx=1 pixel\sigma_{px} = 1\text{ pixel} image noise).

    1. NEES Hypothesis Test: The average Normalized Estimation Error Squared (NEES) across N=50N = 50 runs is ηˉk=1N∑i=1Nϵk(i)TΣ^k−1ϵk(i)\bar{\eta}_k = \frac{1}{N} \sum_{i=1}^N \epsilon_k^{(i)T} \hat{\Sigma}_k^{-1} \epsilon_k^{(i)}, where ϵk=[Log(R^kTRkgt)T,(R^kT(p^k−pkgt))T]T∈R6\epsilon_k = [\text{Log}(\hat{R}_k^T R_k^{gt})^T, (\hat{R}_k^T(\hat{p}_k - p_k^{gt}))^T]^T \in \mathbb{R}^6. For significance level α=2.5%\alpha = 2.5\% and n=6×50=300n = 6 \times 50 = 300 degrees of freedom, the two-sided χ3002\chi_{300}^2 acceptance region is ηˉk∈[5.0,7.0]\bar{\eta}_k \in [5.0, 7.0]. The proposed estimator's average pose NEES remained strictly within this interval and below 7.0 throughout the entire trajectory, demonstrating that the estimator does not become overconfident.

    2. Observability Properties: The estimation uncertainty and errors for roll, pitch, and velocity remained strictly bounded due to gravity observability, whereas yaw orientation and 3D position errors grew slowly over time, matching the theoretical 4 unobservable degrees of freedom in monocular visual-inertial odometry without spurious information gain.

Coverage note — None was omitted; all primary theoretical foundations (on-manifold preintegration, iterative noise propagation, a-posteriori bias corrections, residual Jacobians, structureless factors), graphical model formulations, and empirical evaluations are represented.

References

  1. 1.A. Martinelli. Vision and IMU data fusion: Closed-form solutions for attitude, speed, absolute scale, and bias determination. IEEE Trans. Robotics, 28(1):44–60, 2012.
  2. 2.T. Lupton and S. Sukkarieh. Visual-inertial-aided navigation for high-dynamic motion in built environments without initial conditions. IEEE Trans. Robotics, 28(1):61–76, Feb 2012.
  3. 3.M. Kaess, H. Johannsson, R. Roberts, V. Ila, J. Leonard, and F. Dellaert. iSAM2: Incremental smoothing and mapping using the Bayes tree. Intl. J. of Robotics Research, 31:217–236, Feb 2012.
  4. 4.L. Carlone, Z. Kira, C. Beall, V. Indelman, and F. Dellaert. Eliminating conditionally independent sets in factor graphs: A unifying perspective based on smart factors. In IEEE Intl. Conf. on Robotics and Automation (ICRA), 2014.
  5. 5.A.I. Mourikis and S.I. Roumeliotis. A multi-state constraint Kalman filter for vision-aided inertial navigation. In IEEE Intl. Conf. on Robotics and Automation (ICRA), pages 3565–3572, April 2007.
  6. 6.C. Forster, L. Carlone, F. Dellaert, and D. Scaramuzza. IMU preintegration on manifold for efficient visual-inertial maximum-a-posteriori estimation. In Robotics: Science and Systems (RSS), 2015.
  7. 7.Frank Dellaert. Factor graphs and GTSAM: A hands-on introduction. Technical Report GT-RIM-CP&R-2012-002, Georgia Institute of Technology, September 2012.
  8. 8.K. J. Wu, A. M. Ahmed, G. A. Georgiou, and S. I. Roumeliotis. A square root inverse filter for efficient vision-aided inertial navigation on mobile devices. In Robotics: Science and Systems (RSS), 2015.
  9. 9.Bradley M. Bell and Frederick W. Cathey. The iterated kalman filter update as a gauss-newton method. IEEE Transactions on Automatic Control, 38(2):294–297, 1993.
  10. 10.A.J. Davison, I. Reid, N. Molton, and O. Stasse. MonoSLAM: Real-time single camera SLAM. IEEE Trans. Pattern Anal. Machine Intell., 29(6):1052–1067, Jun 2007.
  11. 11.M. Bloesch, S. Omari, M. Hutter, and R. Siegwart. Robust visual inertial odometry using a direct EKF-based approach. In IEEE/RSJ Intl. Conf. on Intelligent Robots and Systems (IROS). IEEE, 2015.
  12. 12.E.S. Jones and S. Soatto. Visual-inertial navigation, mapping and localization: A scalable real-time causal approach. Intl. J. of Robotics Research, 30(4), Apr 2011.
  13. 13.S. I. Roumeliotis and J. W. Burdick. Stochastic cloning: A generalized framework for processing relative state measurements. In IEEE Intl. Conf. on Robotics and Automation (ICRA). IEEE, 2002.
  14. 14.K. Tsotsos, A. Chiuso, and S. Soatto. Robust inference for visual-inertial sensor fusion. In IEEE Intl. Conf. on Robotics and Automation (ICRA), 2015.
  15. 15.A. Martinelli. Observability properties and deterministic algorithms in visual-inertial structure from motion. Foundations and Trends in Robotics, pages 1–75, 2013.
  16. 16.Dimitrios G. Kottas, Joel A. Hesch, Sean L. Bowman, and Stergios I. Roumeliotis. On the consistency of vision-aided inertial navigation. In Intl. Sym. on Experimental Robotics (ISER), 2012.
  17. 17.G.P. Huang, A.I. Mourikis, and S.I. Roumeliotis. A first-estimates jacobian EKF for improving SLAM consistency. In Intl. Sym. on Experimental Robotics (ISER), 2008.
  18. 18.J.A. Hesch, D.G. Kottas, S.L. Bowman, and S.I. Roumeliotis. Camera-imu-based localization: Observability analysis and consistency improvement. Intl. J. of Robotics Research, 33(1):182–201, 2014.
  19. 19.J. Hernandez, K. Tsotsos, and S. Soatto. Observability, identifiability and sensitivity of vision-aided inertial navigation. In IEEE Intl. Conf. on Robotics and Automation (ICRA), 2015.
  20. 20.A.I. Mourikis and S.I. Roumeliotis. A dual-layer estimator architecture for long-term localization. In Proc. of the Workshop on Visual Localization for Mobile Platforms at CVPR, Anchorage, Alaska, June 2008.
  21. 21.G. Sibley, L. Matthies, and G. Sukhatme. Sliding window filter with application to planetary landing. J. of Field Robotics, 27(5):587–608, 2010.
  22. 22.T-C. Dong-Si and A.I. Mourikis. Motion tracking with fixed-lag smoothing: Algorithm consistency and analysis. In IEEE Intl. Conf. on Robotics and Automation (ICRA), 2011.
  23. 23.S. Leutenegger, P. Furgale, V. Rabaud, M. Chli, K. Konolige, and R. Siegwart. Keyframe-based visual-inertial slam using nonlinear optimization. In Robotics: Science and Systems (RSS), 2013.
  24. 24.S. Leutenegger, S. Lynen, M. Bosse, R. Siegwart, and P. Furgale. Keyframe-based visual-inertial slam using nonlinear optimization. Intl. J. of Robotics Research, 2015.
  25. 25.P. Maybeck. Stochastic Models, Estimation and Control, volume 1. Academic Press, New York, 1979.
  26. 26.G.P. Huang, A.I. Mourikis, and S.I. Roumeliotis. An observability-constrained sliding window filter for SLAM. In IEEE/RSJ Intl. Conf. on Intelligent Robots and Systems (IROS), pages 65–72, 2011.
  27. 27.S-H. Jung and C.J. Taylor. Camera trajectory estimation using inertial sensor measurements and structure fom motion results. In IEEE Conf. on Computer Vision and Pattern Recognition (CVPR), 2001.
  28. 28.D. Sterlow and S. Singh. Motion estimation from image and inertial measurements. 2004.
  29. 29.M. Bryson, M. Johnson-Roberson, and S. Sukkarieh. Airborne smoothing and mapping using vision and inertial sensors. In IEEE Intl. Conf. on Robotics and Automation (ICRA), pages 3143–3148, 2009.
  30. 30.V. Indelman, S. Wiliams, M. Kaess, and F. Dellaert. Information fusion in navigation systems via factor graph based incremental smoothing. Robotics and Autonomous Systems, 61(8):721–738, August 2013.
  31. 31.A. Patron-Perez, S. Lovegrove, and G. Sibley. A spline-based trajectory representation for sensor fusion and rolling shutter cameras. Intl. J. of Computer Vision, February 2015.
  32. 32.H. Strasdat, J.M.M. Montiel, and A.J. Davison. Real-time monocular SLAM: Why filter? In IEEE Intl. Conf. on Robotics and Automation (ICRA), pages 2657–2664, 2010.
  33. 33.G. Klein and D. Murray. Parallel tracking and mapping on a camera phone. In IEEE and ACM Intl. Sym. on Mixed and Augmented Reality (ISMAR), 2009.
  34. 34.E.D. Nerurkar, K.J. Wu, and S.I. Roumeliotis. C-KLAM: Constrained keyframe-based localization and mapping. In IEEE Intl. Conf. on Robotics and Automation (ICRA), 2014.
  35. 35.G. Klein and D. Murray. Parallel tracking and mapping for small AR workspaces. In IEEE and ACM Intl. Sym. on Mixed and Augmented Reality (ISMAR), pages 225–234, Nara, Japan, Nov 2007.
  36. 36.M. Kaess, A. Ranganathan, and F. Dellaert. iSAM: Incremental smoothing and mapping. IEEE Trans. Robotics, 24(6):1365–1378, Dec 2008.
  37. 37.V. Indelman, S. Wiliams, M. Kaess, and F. Dellaert. Factor graph based incremental smoothing in inertial navigation systems. In Intl. Conf. on Information Fusion, FUSION, 2012.
  38. 38.V. Indelman, A. Melim, and F. Dellaert. Incremental light bundle adjustment for robotics navigation. In IEEE/RSJ Intl. Conf. on Intelligent Robots and Systems (IROS), November 2013.
  39. 39.S. Shen. Autonomous Navigation in Complex Indoor and Outdoor Environments with Micro Aerial Vehicles. PhD Thesis, University of Pennsylvania, 2014.
  40. 40.N. Keivan, A. Patron-Perez, and G. Sibley. Asynchronous adaptive conditioning for visual-inertial SLAM. In Intl. Sym. on Experimental Robotics (ISER), 2014.
  41. 41.J. Hornegger and C. Tomasi. Representation issues in the ML estimation of camera motion. In Intl. Conf. on Computer Vision (ICCV), pages 640–647, Kerkyra, Greece, September 1999.
  42. 42.M. Moakher. Means and averaging in the group of rotations. SIAM Journal on Matrix Analysis and Applications, 24(1):1–16, 2002.
  43. 43.G. S. Chirikjian. Stochastic Models, Information Theory, and Lie Groups, Volume 2: Analytic Methods and Modern Applications (Applied and Numerical Harmonic Analysis). Birkhauser, 2012.
  44. 44.Y. Wang and G.S. Chirikjian. Nonparametric second-order theory of error propagation on motion groups. Intl. J. of Robotics Research, 27(11–12):1258–1273, 2008.
  45. 45.P.A. Absil, C.G. Baker, and K.A. Gallivan. Trust-region methods on Riemannian manifolds. Foundations of Computational Mathematics, 7(3):303–330, 2007.
  46. 46.T. D. Barfoot and P. T. Furgale. Associating uncertainty with three-dimensional poses for use in estimation problems. IEEE Trans. Robotics, 30(3):679–693, 2014.
  47. 47.Y. Wang and G.S. Chirikjian. Error propagation on the euclidean group with applications to manipulator kinematics. IEEE Trans. Robotics, 22(4):591–602, 2006.
  48. 48.S. T. Smith. Optimization techniques on Riemannian manifolds. Hamiltonian and Gradient Flows, Algorithms and Control, Fields Inst. Commun., Amer. Math. Soc., 3: 113–136, 1994.
  49. 49.J.A. Farrell. Aided Navigation: GPS with High Rate Sensors. McGraw-Hill, 2008.
  50. 50.J. H. Manton. Optimization algorithms exploiting unitary constraints. IEEE Trans. Signal Processing, 50:63–650, 2002.
  51. 51.F.R. Kschischang, B.J. Frey, and H-A. Loeliger. Factor graphs and the sum-product algorithm. IEEE Trans. Inform. Theory, 47(2), February 2001.
  52. 52.F. Dellaert. Square Root SAM: Simultaneous localization and mapping via square root information smoothing. Technical Report GIT-GVU-05-11, Georgia Institute of Technology, 2005.
  53. 53.R.M. Murray, Z. Li, and S. Sastry. A Mathematical Introduction to Robotic Manipulation. CRC Press, 1994.
  54. 54.P. E. Crouch and R. Grossman. Numerical integration of ordinary differential equations on manifolds. J. Nonlinear Sci., 3:1–22, 1993.
  55. 55.H. Munthe-Kaas. Higher order runge-kutta methods on manifolds. Appl. Numer. Math., 29(1):115–127, 1999.
  56. 56.J. Park and W.-K. Chung. Geometric integration on euclidean group with application to articulated multibody systems. IEEE Trans. Robotics, 2005.
  57. 57.M. S. Andrle and J. L. Crassidis. Geometric integration of quaternions. Journal of Guidance, Control, and Dynamics, 36(6):1762–1757, 2013.
  58. 58.J.L. Crassidis. Sigma-point Kalman filtering for integrated GPS and inertial navigation. IEEE Trans. Aerosp. Electron. Syst., 42(2):750–756, 2006.
  59. 59.P. Furgale, J. Rehder, and R. Siegwart. Unified temporal and spatial calibration for multi-sensor systems. In IEEE/RSJ Intl. Conf. on Intelligent Robots and Systems (IROS), 2013.
  60. 60.M. Li and A.I. Mourikis. Online temporal calibration for camera-imu systems: Theory and algorithms. Intl. J. of Robotics Research, 33(6), 2014.
  61. 61.R. I. Hartley and A. Zisserman. Multiple View Geometry in Computer Vision. Cambridge University Press, second edition, 2004.
  62. 62.M. Kaess, V. Ila, R. Roberts, and F. Dellaert. The Bayes tree: Enabling incremental reordering and fluid relinearization for online mapping. Technical Report MIT-CSAIL-TR-2010-021, Computer Science and Artificial Intelligence Laboratory, MIT, Jan 2010.
  63. 63.Y. Bar-Shalom, X. R. Li, and T. Kirubarajan. Estimation with Applications To Tracking and Navigation. John Wiley and Sons, 2001.
  64. 64.J. Hornegger. Statistical modeling of relations for 3-D object recognition. In Intl. Conf. Acoust., Speech, and Signal Proc. (ICASSP), volume 4, pages 3173–3176, Munich, April 1997.
  65. 65.C. Forster, M. Pizzoli, and D. Scaramuzza. SVO: Fast Semi-Direct Monocular Visual Odometry. In IEEE Intl. Conf. on Robotics and Automation (ICRA), 2014. doi: 10.1109/ICRA.2014.6906584.
  66. 66.J. Nikolic, J., Rehderand M. Burri, P. Gohl, S. Leutenegger, P. Furgale, and R. Siegwart. A Synchronized Visual-Inertial Sensor System with FPGA Pre-Processing for Accurate Real-Time SLAM. In IEEE Intl. Conf. on Robotics and Automation (ICRA), 2014.
  67. 67.F.C. Park and B.J. Martin. Robot sensor calibration: Solving AX=XB on the euclidean group. 10(5), 1994.
  68. 68.A. Geiger, P. Lenz, and R. Urtasun. Are we ready for autonomous driving? the KITTI vision benchmark suite. In IEEE Conf. on Computer Vision and Pattern Recognition (CVPR), pages 3354–3361, Providence, USA, June 2012.
  69. 69.Tango IMU specifications. URL http://ae-bst.resource.bosch.com/media/products/dokumente/bmx055/BST-BMX055-FL000-00 2013-05-07 onl.pdf.
  70. 70.Adis IMU specifications. URL http://www.analog.com/media/en/technical-documentation/data-sheets/ADIS16448.pdf.
  71. 71.C.D. Meyer. Matrix Analysis and Applied Linear Algebra. SIAM, 2000.
  72. 72.J. Diebel. Representing attitude: Euler angles, unit quaternions, and rotation vectors. Technical report, Stanford University, 2006.

Citation

MLA
Forster, C., et al. “On-Manifold Preintegration for Real-Time Visual--Inertial Odometry”. IEEE Transactions on Robotics, vol. 33, no. 1, 2017, pp. 1–1, https://doi.org/10.1109/TRO.2016.2597321.
APA
Forster, C., Carlone, L., Dellaert, F., & Scaramuzza, D. (2017). On-Manifold Preintegration for Real-Time Visual--Inertial Odometry. IEEE Transactions on Robotics, 33(1), 1–21. https://doi.org/10.1109/TRO.2016.2597321
Chicago
Forster, C., L. Carlone, F. Dellaert, and D. Scaramuzza. 2017. “On-Manifold Preintegration for Real-Time Visual--Inertial Odometry”. IEEE Transactions on Robotics 33 (1): 1–21. https://doi.org/10.1109/TRO.2016.2597321.
Harvard
Forster, C. et al. (2017) “On-Manifold Preintegration for Real-Time Visual--Inertial Odometry”, IEEE Transactions on Robotics, 33(1), pp. 1–21. Available at: https://doi.org/10.1109/TRO.2016.2597321.
Vancouver
1. Forster C, Carlone L, Dellaert F, Scaramuzza D (2017) On-Manifold Preintegration for Real-Time Visual--Inertial Odometry. IEEE Transactions on Robotics 33:1–21

BibTeX

@article{Forster_2017, title={On-Manifold Preintegration for Real-Time Visual--Inertial Odometry}, volume={33}, ISSN={1941-0468}, url={http://dx.doi.org/10.1109/TRO.2016.2597321}, DOI={10.1109/tro.2016.2597321}, number={1}, journal={IEEE Transactions on Robotics}, publisher={Institute of Electrical and Electronics Engineers (IEEE)}, author={Forster, Christian and Carlone, Luca and Dellaert, Frank and Scaramuzza, Davide}, year={2017}, month=Feb, pages={1–21} }
Metadata:Crossref

Source Code

This paper has an official code repository available. Click below to access the source code.

View Repository

Access the Paper

This paper is available from its original source. Click below to access the PDF.

Open PDF