iSAM: Incremental Smoothing and Mapping

iSAM: Incremental
Smoothing and Mapping
Michael Kaess (Student Member IEEE)
Ananth Ranganathan(Student Member IEEE)
Frank Dellaert(Member IEEE)
Created By: Akanksha, October 2015
Key Idea

iSAM performs fast incremental updates of the information matrix thus
avoiding unnecessary calculations.

Periodic variable reordering prevents unnecessary fill-in in the square root
factor.

Online data association, hence relevant estimation uncertainties are
retrieved from incrementally updated square root factor exploiting on
sparsity of the full covariance matrix.
Created By: Akanksha, October 2015
Related Work

Square root SAM: Simultaneous localization and mapping via square root
information smoothing by F. Dellaert and M. Kaess (2006)

Square Root SAM: Simultaneous location and mapping via square root
information smoothing by F. Dellaert (2005)

Stochastic Models, Estimation and Control by P. Maybeck (1079)

Factorization method for discrete sequential estimation by G. Bierman (1977)
Created By: Akanksha, October 2015
SAM Exposition
𝑀
𝑃 𝑋, 𝐿, 𝑈, 𝑍 ∝ 𝑃(𝑥0 )
𝑃(𝑥𝑖 |𝑥𝑖−1 , 𝑢𝑖 )
𝑖=1
Created By: Akanksha, October 2015
𝐾
𝑃(𝑧𝑘 |𝑥𝑖𝑘 , 𝑙𝑗𝑘 )
𝑘=1
Assumption – Gaussian Measurement Models

𝑥𝑖 = 𝑓𝑖 𝑥𝑖−1 , 𝑢𝑖 + 𝑤𝑖
where 𝑤𝑖 is process noise as Ν(0, Λ𝑖 )

𝑧𝑘 = ℎ𝑘 𝑥𝑖𝑘 , 𝑙𝑗𝑘 + 𝜈𝑖
where 𝑣𝑖 is process noise as Ν(0, Γ𝑖 )
Created By: Akanksha, October 2015
SLAM as Least Square Problem
maximum a posteriori (MAP) estimate
𝑋 ∗ , 𝐿∗ = argmax 𝑃(𝑋, 𝐿, 𝑈, 𝑍) = argmin − log 𝑃(𝑋, 𝐿, 𝑈, 𝑍)
𝑋,𝐿
𝑋,𝐿
In terms of process and measurement models
𝑀
∥ 𝑓𝑖 𝑥𝑖−1 , 𝑢𝑖 − 𝑥𝑖 ∥2Λ𝑖 +
𝑋 ∗ 𝐿∗ = 𝑎𝑟𝑔min{
𝑋,𝐿
𝐾
𝑖=1
∥ ℎ𝑘 𝑥𝑖𝑘 , 𝑙𝑗𝑘 − 𝑧𝑘 ∥2Γ𝑖 }
𝑘=1
Where ∥ 𝑒 ∥Σ = 𝑒 𝑇 Σ −1 𝑒 for the squared Mahalanobis distance
Linearization of the system : 𝜃 ∗ = argmin ∥ 𝐴𝜃 − 𝑏 ∥2
𝜃
Created By: Akanksha, October 2015
QR Factorization
Substitute 𝐴 = 𝑄
𝑅
in linearized least square problem:
0
𝐴𝜃 − 𝑏 2 = 𝑅𝜃 − 𝑑 2 + 𝑒 2
Where 𝑑, 𝑒 𝑇 ≔ 𝑄𝑇 𝑏 as derived in the paper.
The above is minimum if and only if 𝑅𝜃 = 𝑑, thus it can be implied that
𝑅𝜃 ∗ = 𝑑
Created By: Akanksha, October 2015
iSAM: Incremental Smoothing and Mapping

Matrix Factorization by Givens Rotations
Created By: Akanksha, October 2015
iSAM: Incremental Smoothing and Mapping

Incremental Updating
Created By: Akanksha, October 2015
iSAM: Incremental Smoothing and Mapping

Incremental SAM
Maximum number of Givens Rotations needed for adding a new measurement
row is 𝑛. Due to sparsity of R and the new measurement row, constant Givens
Rotations are needed in this method.
Created By: Akanksha, October 2015
Loops and Periodic Variable Reordering
Created By: Akanksha, October 2015
Data Association

Maximum likelihood Data Association (ML)

Marginal Covariances
Created By: Akanksha, October 2015
Experimental Results and Discussion

Landmark based iSAM
Created By: Akanksha, October 2015
Experimental Results and Discussion

Pose Constraint-based iSAM
Created By: Akanksha, October 2015
Experimental Results and Discussion

Sparsity of square root factor
Created By: Akanksha, October 2015
Questions?
Created By: Akanksha, October 2015
Thank you!
Created By: Akanksha, October 2015