ORB-SLAM-Mono论文要点总结(六)——— 附录-非线性优化 Nonlinear Optimizations

非线性优化

  1. BA优化[1]
    地图点的3-D位置坐标为 X w , j ∈ R 3 X_{w,j}\in R^3 Xw,jR3,关键帧位姿变换 T i w ∈ S E ( 3 ) T_{iw}\in SE(3) TiwSE(3)(变换空间)。 w w w表示世界参考系,优化使得地图点和与之匹配的特征点 x i , j ∈ R 2 x_{i,j}\in R^2 xi,jR2间的重投影误差最小。在关键帧 i i i中一个地图点 j j j的观测误差为:
    e i , j = x i , j − π i ( T i w , X w , j ) e_{i,j} = x_{i,j} - \pi_i({T_{iw},X_{w,j}}) eij=xi,jπi(TiwXw,j)
    这里 π \pi π为投影函数:
    ORB-SLAM-Mono论文要点总结(六)——— 附录-非线性优化 Nonlinear Optimizations_第1张图片
    图中世界系到 i i i关键帧的旋转矩阵 R i w ∈ S O ( 3 ) R_{iw}\in SO(3) RiwSO(3)旋转矩阵李群空间,位移变换 t i w ∈ R 3 t_{iw}\in R^3 tiwR3,分为为变换矩阵 T i w T_{iw} Tiw的一部分。 ( f i , u , f i , v ) (f_{i,u},f_{i,v}) fi,ufi,v ( c i , u , c i , v ) (c_{i,u},c_{i,v}) (ci,uci,v)分别表示与相机 i i i对应的焦距长度和主点坐标。最小化代价函数为:
    C = ∑ i , j ρ h ( e i T Ω i , j − 1 e i , j ) C = \sum_{i,j}\rho_h(e_i^T\Omega_{i,j}^{-1}e_{i,j}) C=i,jρh(eiTΩi,j1ei,j)上式中 ρ h \rho_h ρh是huber鲁棒损失函数, Ω i , j = σ i , j 2 I 2 ∗ 2 \Omega_{i,j} = \sigma_{i,j}^2I_{2*2} Ωi,j=σi,j2I22为检测到的特征点与尺度相关的协方差矩阵。
    下面是基于以上重投影误差原理和建立的代价函数在不同工作模块工作时,对应的优化方式,固定不变的状态,以及每种优化方式待优化的状态:
工作模块 优化方式 固定项 待优化项
地图初始化 full BA 第一个关键帧 所有地图点和关键帧
局部建图 local BA 关键帧 局部地图里所有的地图点
跟踪 位姿优化/motion-only BA 所有地图点 相机位姿
  1. 基于相似变换空间Sim3约束的位姿图优化[6]
    基于两个边的位姿图(见VII-D), 可以定义一边的误差为:
    e i , j = l o g S i m ( 3 ) ( S i j S j w S i w − 1 ) e_{i,j} = log_{Sim(3)}(S_{ij}S_{jw}S_{iw}^{-1}) ei,j=logSim(3)(SijSjwSiw1)上式中 S i j S_{ij} Sij是图优化之前基于SE(3)位姿计算的两个关键帧的相对显示变换,令尺度系数为1?。在回环边,基于Horn[42]方法计算相对变换。 l o g S i m ( 3 ) log_{Sim(3)} logSim(3)[48]转换到切空间为了将误差转换为7维向量?。优化相似变换空间的的关键帧位姿时最小化的代价函数为:
    C = ∑ i , j ( e i , j T Λ i , j e i , j ) C = \sum_{i,j}(e_{i,j}^T\Lambda_{i,j}e_{i,j}) C=i,j(ei,jTΛi,jei,j)上式中 Λ i , j \Lambda_{i,j} Λi,j如[48]中所述的信息矩阵,作者将其设为单位阵。为了固定规范7自由度,固定回环关键帧。尽管此方法是一个大致的全局BA(full BA),在实验VIII-E中证明比BA能够更快更好的收敛。
  2. 关键帧间的相似变换优化
    基于关键帧1和关键帧2间给定的一组n个匹配 i ⇒ j i\Rightarrow j ij(特征点和它们对应的3D地图点),为了优化相对位姿 S 12 S_{12} S12(见VII-B),可以最小化两帧图像间的重投影误差为:
    ORB-SLAM-Mono论文要点总结(六)——— 附录-非线性优化 Nonlinear Optimizations_第2张图片
    代价函数为:
    在这里插入图片描述
    上式中 Ω 1 , i \Omega_{1,i} Ω1,i Ω 2 , j \Omega_{2,j} Ω2,j为与图像1和图像2中特征点尺度相关的协方差度量。在此优化过程中,地图点时固定的。

你可能感兴趣的:(SLAM)