DTAM: Dense Tracking and Mapping in Real-Time, Robot vision Group
14 likes5,872 views
The document details DTAM, a real-time camera tracking and dense mapping system that utilizes dense, pixel-based methods instead of traditional feature extraction. It outlines system components, dense mapping and tracking methodologies, and presents evaluations demonstrating its performance against other methods. Future work is suggested to improve robustness to lighting variations and depth assumptions.
![Key observation
● Data term can be globally optimized by
exhaustive search(point-wise optimization)
● Convex regularization term can be solved
efficiently by convex optimization algorithms
● And we can approximate the energy functional
by decoupling data and regularity term following
the approach described in [1][2]
[1]F.Steinbrucker et.al: Large displacement optical flow computation without warping
[2]A.Chambolle et.al: An Algorithm for Total Variation Minimization and Applications](https://image.slidesharecdn.com/dtam-gm-lihang-robotvison-140528043940-phpapp01/85/DTAM-Dense-Tracking-and-Mapping-in-Real-Time-Robot-vision-Group-24-320.jpg){kind=tracking}
![Primal-Dual Approach
● General class of energy minimization problems:
● Can obtain dual form by replacing F(Kx) by its convex conjugate
F*(y)
● Use duality principles to arrive at the primal-dual form of
following [1][2][3]
* Usually regularization term
* Often a norm:||Kx|| * Data term
g(u)‖∇ξ(u)‖ϵ+Q(u)
[1] J.-F. Aujol. Some first-order algorithms for total variation based image restoration
[2] A. Chambolle and T. Pock. A first-order primal-dual algorithm for convex problems
with applications to imaging
[3] M.Zhu. Fast numerical algorithms for total variation based image restoration](https://image.slidesharecdn.com/dtam-gm-lihang-robotvison-140528043940-phpapp01/85/DTAM-Dense-Tracking-and-Mapping-in-Real-Time-Robot-vision-Group-29-320.jpg){kind=tracking}
![Algorithm
● Compute partial derivatives
-
-
● For fixed a, gradient ascent w.r.t q and gradient descent
w.r.t d is performed
● For fixed d, exhaustive search w.r.t a is performed
● is decremented
● Until >
θ
θn
θend
[1] A. Chambolle and T. Pock. A first-order primal-dual algorithm for convex problems
with applications to imaging](https://image.slidesharecdn.com/dtam-gm-lihang-robotvison-140528043940-phpapp01/85/DTAM-Dense-Tracking-and-Mapping-in-Real-Time-Robot-vision-Group-35-320.jpg){kind=tracking}
![SSD optimization
● Problem:
- Align template image T(x) with input image I(x)
● Formulation:
- Find the transform that best maps the pixels of the templates into
the ones of the current image minimizing:
- are the displacement parameters to be optimized
● Hypothesis:
- Known a coarse approximation of the template position
Σ
x
[I(W(x;p))−T(x)]2
W(x; p)
p=( p1 ,..., pn)T
(p0)](https://image.slidesharecdn.com/dtam-gm-lihang-robotvison-140528043940-phpapp01/85/DTAM-Dense-Tracking-and-Mapping-in-Real-Time-Robot-vision-Group-40-320.jpg){kind=tracking}
![SSD optimization
● Problem:
- minimize
- The current estimation of p is iteratively updated to reach the minimum of the function.
● Formulations:
- Direct additional
- Direct compositional
-Inverse
Σ
x
[I(W(x;p))−T(x)]2
Σ
x
[I(W(x;p+Δ p))−T(x)]2
Σ
x
[I(W(W(x;Δ p);p))−T(x)]
2
Σ
x
[I(W(x;Δ p))−I(W(x;p))]
2](https://image.slidesharecdn.com/dtam-gm-lihang-robotvison-140528043940-phpapp01/85/DTAM-Dense-Tracking-and-Mapping-in-Real-Time-Robot-vision-Group-41-320.jpg){kind=tracking}
![SSD optimization
●
Example: Direct additive method
- Minimize:
- First order Taylor expansion:
- Solution:
- with:
Σ
x
[I(W(x;p+Δ p))−T(x)]2
Σ
x
[I(W(x;p))+∇ I
∂W
∂p
Δ p−T(x)]
2
Δ p=Σ
x
H
−1
[∇ I
∂W
∂p
]
T
[T(x)−I(W(x;p))]
H=Σ
x
[∇ I
∂W
∂p
]
T
[∇ I
∂W
∂p
]](https://image.slidesharecdn.com/dtam-gm-lihang-robotvison-140528043940-phpapp01/85/DTAM-Dense-Tracking-and-Mapping-in-Real-Time-Robot-vision-Group-42-320.jpg){kind=tracking}
![SSD robustified
● Formulation:
● Problem: In case of occlusion, the occluded pixels
cause the optimum of the function to be changed.
The occluded pixels have to be ignored from the optimization
● Method
- Only the pixels with a difference
lower than a threshold are selected
- Threshold is iteratively updated to get more
selective as the optimization reaches the optimum
Δ p=Σ
x
H
−1
[∇ I
∂W
∂p
]
T
[T(x)−I(W(x;p))]
[T (x)−I (W (x; p))]](https://image.slidesharecdn.com/dtam-gm-lihang-robotvison-140528043940-phpapp01/85/DTAM-Dense-Tracking-and-Mapping-in-Real-Time-Robot-vision-Group-43-320.jpg){kind=tracking}
DTAM: Dense Tracking and Mapping in Real-Time, Robot vision Group
- 1. DTAM:Dense Tracking and Mapping in Real-Time NewCombe, Lovegrove & Davision ICCV11
- 2. Outline ➢ Introduction ➢ Related Work ➢ System Overview ➢ Dense Mapping ➢ Dense Tracking ➢ Evaluation and Results ➢ Conclusions and Future Work
- 3. Outline ➢ Introduction ➢ Related Work ➢ System Overview ➢ Dense Mapping ➢ Dense Tracking ➢ Evaluation and Results ➢ Conclusions and Future Work
- 4. Introduction ● Dense Tracking and Mapping in Real-Time ● DTAM is a system for real-time camera tracking and reconstruction which relies not on feature extraction but dense, every pixel methods. ● Simultaneous frame-rate Tracking and Dense Mapping
- 5. Outline ➢ Introduction ➢ Related Work ➢ System Overview ➢ Dense Mapping ➢ Dense Tracking ➢ Evaluation and Results ➢ Conclusions and Future Work
- 6. Related Work ● Real-time SFM(Structure from Motion) ● PTAM(G. Klein and D. W. Murray. ISMAR 2007) ● Improving the agility of keyframe-based SLAM(G. Klein and D. W. Murray. ECCV 2008) ● Live dense reconstruction with a single moving camera(R. A. Newcombe and A. J. Davison. CVPR 2010) ● Real-time dense geometry from a handheld camera(J. Stuehmer et.al. 2010)
- 7. Outline ➢ Introduction ➢ Related Work ➢ System Overview ➢ Dense Mapping ➢ Dense Tracking ➢ Evaluation and Results ➢ Conclusions and Future Work
- 8. System Overview ● Input -Single hand held RGB Camera ● ● Objective: -Dense Mapping -Dense Tracking Input Imgage 3D Dense Map
- 10. Outline ➢ Introduction ➢ Related Work ➢ System Overview ➢ Dense Mapping ➢ Dense Tracking ➢ Evaluation and Results ➢ Conclusions and Future Work
- 11. Dense Mapping ● Estimate inverse depth map from bundles of images
- 12. Photometric error ● Total cost ● Photometric error ● Where: - K ...intrinsic matrix - Tmr ...transformation from m to r - -
- 13. Depth map estimation ● Principle: - S depth hypothesis are considered for each pixel of the reference image Ir - Each corresponding 3D point is projected onto a bundle of images Im - Keep the depth hypothesis that best respects the color consistency from the reference to the bundle of images ● Formulation: - :pixel position and depth hypothesis - :number of valid reprojection of the pixel in the bundle - :photometric error between reference and current image
- 14. Depth map estimation Example reference image pixel Reprojection of depth Hypotheses on one image of bundle Rerojectionin ImagebundlePhotomtricerror Depth Hypotheses
- 15. Inverse Depth Map Computation ● Inverse depth map can be computed by minimizing the photometric error(exhaustive search ove the volume): ● But featureless regions are prone to false minima
- 16. Inverse Depth Map Computation
- 17. Depth map filtering approach ● Problem: - Uniform regions in reference image do not give discriminative enough photometric error ● Idea: - Assume the depth is smooth on uniform regions - Use total variational approach where depth map is the functional to optimize: *photometric error defines the data term *the smoothness constraint defines the regularization
- 18. Inverse Depth Map Computation ● Featureless regions are prone to false minima ● Solution:Regularization term - We want to penalize deviation from spatially smooth solution - But preserve edges and discontinuities
- 19. Depth map filtering approach ● Formulation of the variational approach - First term: regularization constraint, g is defined as 0 for image gradients and 1 for uniform regions. So that gradient on depth map is penalized for uniform regions - Second term: data term defined by the photometric error - Huber norm: differentiable replacement to L1 norm that better preserve discontinuities compared to L2
- 20. Energy Functional ● Regularised cost Huber norm Weight Regularization term Photometric cost term
- 21. Total Variation(TV) Regularization ● L1 penalization of gradient magnitudes - Favors sparse, piecewise-constant solutions - Allows sharp discontinuities in the solution ● Problem - Staircasing - Can be reduced by using quadratic penalization for small gradient magnitudes
- 22. Energy Functional Analysis Convex Not convex
- 23. Energy Minimization ● Composition of both terms is non-convex fuction ● Possible solution - Linearize the cost volume to get a convex approximation of the data term - Solve approximation iteratively within coarse-to-fine warping scheme * Can lead in loss of the reconstruction details ● Better solution?
- 24. Key observation ● Data term can be globally optimized by exhaustive search(point-wise optimization) ● Convex regularization term can be solved efficiently by convex optimization algorithms ● And we can approximate the energy functional by decoupling data and regularity term following the approach described in [1][2] [1]F.Steinbrucker et.al: Large displacement optical flow computation without warping [2]A.Chambolle et.al: An Algorithm for Total Variation Minimization and Applications
- 25. Alternating two Global Optimizations *Drives original and aux. Variables together *Minimizing functional above equivalent to minimizing original formulation as θ -> 0 *Data and regularity terms are decoupled via aux. Variable α *Optimization process is split into two sub-prolems α=Ω→ℝ
- 26. Alternating two Global Optimizations ● Energy functional can be globally minimized w.r.t ξ * Since it is convex in ξ * E.g. gradient descent ● Energy functional can be globally minimized w.r.t α * Not convex w.r.t α, but trivially point-wise optimizable * Exhaustive search
- 27. Algorithm ● Initialization - Compute = = - θ = large_value ● Until > - Compute * Minimize with fixed α * Use convex optimization tools, e.g. gradient descent - Compute * Minize with fixed ξ * Exhaustive search – Decrement θ θn θend αu 0 ξu 0 mind C(u,d) ξu n αu n Eξ, α Eξ, α
- 28. Even better ● Problem - optimization badly conditioned as (uniform regions) - expensive when doing exhaustive search - accuracy is not good enough ● Solution - Primal-Dual approach for convex optimization step - Acceleration of non-convex search - Sub-pixel accuracy ∇u→0
- 29. Primal-Dual Approach ● General class of energy minimization problems: ● Can obtain dual form by replacing F(Kx) by its convex conjugate F*(y) ● Use duality principles to arrive at the primal-dual form of following [1][2][3] * Usually regularization term * Often a norm:||Kx|| * Data term g(u)‖∇ξ(u)‖ϵ+Q(u) [1] J.-F. Aujol. Some first-order algorithms for total variation based image restoration [2] A. Chambolle and T. Pock. A first-order primal-dual algorithm for convex problems with applications to imaging [3] M.Zhu. Fast numerical algorithms for total variation based image restoration
- 30. Primal-Dual Approach ● General problem formulation: ● By definition(Legendre-Fenchel transform): ● Dual Form(Saddle-point problem):
- 31. Primal-Dual Approach ● Conjugate of Huber norm(obtained via Legendre-Fenchel transform)
- 32. Minimization ● Solving a saddle point problem now! ● Condition of optimality met when ● Compute partial derivatives - - ● Perform gradient descent - Ascent on y(maximization) - Descent on x(minimization)
- 33. Discretisation ● First some notation: - Cost volume is discretized in M X N X S array * M X N ... reference image resolution * S ... number of points linearly sampling the inverse depth range - Use MN X 1 stacked rasterised column vector * d ... vector version of ξ * a ... vector version of α * g ... MN X 1 vector with per-pixel weights * G=diag(g) ... element-wise weighting matrix - Ad computes 2MN X 1 gradient vector
- 34. Implementation ● Replace the weighted Huber regularizer by its conjugate ● Saddle-point problem - Primal variable d and dual variable q - Coupled with data term * Sum of convex and non-convex functions F(AGd) F*(q) F*(q)
- 35. Algorithm ● Compute partial derivatives - - ● For fixed a, gradient ascent w.r.t q and gradient descent w.r.t d is performed ● For fixed d, exhaustive search w.r.t a is performed ● is decremented ● Until > θ θn θend [1] A. Chambolle and T. Pock. A first-order primal-dual algorithm for convex problems with applications to imaging
- 36. Outline ➢ Introduction ➢ Related Work ➢ System Overview ➢ Dense Mapping ➢ Dense Tracking ➢ Evaluation and Results ➢ Conclusions and Future Work
- 37. Dense Tracking ● Inputs: - 3D texture model of the scene - Pose at previous frame ● Tracking as a registration problem - First inter-frame rotation estimation: the previous image is aligned on the current image to estimate a coarse inter-frame rotation - Estimated pose is used to project the 3D model into 2.5D image - The 2.5D image is registered with the current frame to find the current pose Template matching problem
- 38. Tracking Strategy and Algorithm ● Based on image alignment against dense model ● Coarse-to-fine strategy - Pyramid hierarchy of images ● Lucas-Kanade algorithm - Estimate “warp” between images - Iterative minimization of a cost function - Parameters of warp correspond to dimensionality of search space
- 39. Tracking in Two Stages ● Two stages - Constrained rotation estimation * Use coarser scales * Rough estimate of pose - Accurate 6-DOF pose refinement * Set virtual camera at location Project dense model to the virtual camera Image ,inverse depth image * Align live image and to estimate * Final pose estimation T ^ wl ν Tw ν=T ^ wl Iν ξν Il Iν Tlv Twl=Tw νTlν
- 40. SSD optimization ● Problem: - Align template image T(x) with input image I(x) ● Formulation: - Find the transform that best maps the pixels of the templates into the ones of the current image minimizing: - are the displacement parameters to be optimized ● Hypothesis: - Known a coarse approximation of the template position Σ x [I(W(x;p))−T(x)]2 W(x; p) p=( p1 ,..., pn)T (p0)
- 41. SSD optimization ● Problem: - minimize - The current estimation of p is iteratively updated to reach the minimum of the function. ● Formulations: - Direct additional - Direct compositional -Inverse Σ x [I(W(x;p))−T(x)]2 Σ x [I(W(x;p+Δ p))−T(x)]2 Σ x [I(W(W(x;Δ p);p))−T(x)] 2 Σ x [I(W(x;Δ p))−I(W(x;p))] 2
- 42. SSD optimization ● Example: Direct additive method - Minimize: - First order Taylor expansion: - Solution: - with: Σ x [I(W(x;p+Δ p))−T(x)]2 Σ x [I(W(x;p))+∇ I ∂W ∂p Δ p−T(x)] 2 Δ p=Σ x H −1 [∇ I ∂W ∂p ] T [T(x)−I(W(x;p))] H=Σ x [∇ I ∂W ∂p ] T [∇ I ∂W ∂p ]
- 43. SSD robustified ● Formulation: ● Problem: In case of occlusion, the occluded pixels cause the optimum of the function to be changed. The occluded pixels have to be ignored from the optimization ● Method - Only the pixels with a difference lower than a threshold are selected - Threshold is iteratively updated to get more selective as the optimization reaches the optimum Δ p=Σ x H −1 [∇ I ∂W ∂p ] T [T(x)−I(W(x;p))] [T (x)−I (W (x; p))]
- 44. Template matching in DTAM ● Inter-frame rotation estimation - the template is the previous image that is matched with current image. Warp is defined on . The initial estimate of p is identity. ● Full pose estimation - template is 2.5D, warp is defined by full 3D motion estimation, which is on - The initial pose is given by the pose estimated at the previous frame and the inter-frame rotation estimation SO(3) SE(3)
- 45. 6 DOF Image Alignment ● Gauss-Newton gradient descent non-linear optimization ● Non-linear expression linearized by first-order Taylor expansion F(ψ)= 1 2 Σ u∈Ω (fu(ψ)) 2 fu(ψ)=Il(π(KTlν(ψ)π−1 (u,ξυ(u))))−Iυ(u) Tlv (ψ)=exp( Σ i=1 6 ψi gen SE(3) i) ψ∈R6 BelongstoLieAlgebraςϱ3
- 46. Outline ➢ Introduction ➢ Related Work ➢ System Overview ➢ Dense Mapping ➢ Dense Tracking ➢ Evaluation and Results ➢ Conclusions and Future Work
- 47. Evaluation and Results ● Runs in real-time - NVIDIA GTX 480 GPU - i7 quad-core GPU - Grey Flea2 camera * Resolution 640*480 * 30Hz ● Comparison with PTAM - a challenging high acceleration back-and-forth trajectory close to a cup - with DTAM's relocaliser disabled
- 48. Evaluation and Results ● Unmodelled objects ● Camera defocus
- 49. Outline ➢ Introduction ➢ Related Work ➢ System Overview ➢ Dense Mapping ➢ Dense Tracking ➢ Evaluation and Results ➢ Conclusions and Future Work
- 50. Conclusions ● First live full dense reconstruction system ● Significant advance in real-time geometrical vision ● Robust - rapid motion - cemera defocus ● Dense modelling and dense tracking make the system beat any point-based method with modelling and tracking performance
- 51. Future Work ● Short comings - Brightness constancy assumption * often violated in real-world * not robust to global illumination changes - Smoothness assumption on depth ● Possible solutions - integrate a normalized cross correlation measure into the objective function for more robustness to local and global lighting changes - joint modelling of the dense lighting and reflectance properties of the scene to enable moe accurate photometric cost functions(the authors are more interested in this approach)
- 52. Thank You!