SLAM++ 增量式BA优化库

incSLAM++是3DV2017的最佳论文,亮点在于其增量式Schur补更新加速BA求解,针对稀疏变化有效缩短计算时间。此外,文章介绍了一种高效协方差恢复技术,降低边缘协方差计算成本。文中通过与g2oSLAM++和Ceres比较,突出了稀疏性和外点抑制的优势。

incSLAM++是3DV 2017最佳论文,源文件网址 https://sourceforge.net/p/slam-plus-plus/wiki/Home/

incSLAM++有两个创新点:
一是增量式Schur补更新。
如果delta更新是稀疏的,则采用增量式Schur补来求解非线性最小二乘问题,即增量式BA问题,会缩短计算时间。

求解不同关键帧之间的位姿和每帧上的点即求解增量式BA(解不是全局的),就是求解如下非线性最小二乘问题
在这里插入图片描述

在每次迭代时,在初始点附近利用泰勒展开进行线性化,得到如下最小二乘问题

在这里插入图片描述

在使用最小二乘求解过程中将相机位姿和路标点分离,则LS求解的系统矩阵和右边为
在这里插入图片描述

因为在线性化点附近小的改变可以忽略 所以系统矩阵用增量表达更新

在这里插入图片描述

增量式Schur补变为

在这里插入图片描述


在这里插入图片描述

如果所有的变量都更新的话,P和deltaP-1就有同样的秩,这时问题就退化到批处理了,但是如果没有更新所有变量,当前的方法就有加速作用。

二是一种高效的协方差恢复技术。
可以使边缘协方差(marginal covariance)计算和NLS的计算时间在一个数量级上。

系统矩阵的逆就是协方差。这里采用Woodbury formula 求得了协方差,将Schur补的求逆改为Schur补的迭代。

文章强调了所提方法的优越性在于更新的稀疏性,当变量更新小于50%时候,计算代价更小。文章还采用了损失函数来削弱外点的影响。在与g2o slam++ ceres作对比的时候,g2o与slam++均采用LM solver,incSlam++与ceres均采用dogleg算法。效果对比
在这里插入图片描述

<think>嗯,用户问的是ORB-SLAM中的BA优化实现和改进方法。首先,我需要回忆一下ORB-SLAM的基本结构。ORB-SLAM主要分为跟踪、局部建图和闭环检测三个线程,而BA优化主要在局部建图和闭环检测中使用。 BA(Bundle Adjustment)在ORB-SLAM中应该是在局部地图优化和全局优化时进行的。记得ORB-SLAM使用的是g2o作为优化后端,所以需要确认这一点。用户可能想了解具体的实现细节,比如如何构建图优化模型,使用的节点和边是什么类型的。 在实现部分,应该分步骤说明:数据结构准备、图模型构建、优化策略、参数设置。比如,关键帧和地图点作为顶点,观测关系作为边。需要提到使用的鲁棒核函数,比如Huber,以及优化后的外点处理机制。 改进方法方面,ORB-SLAM3可能对BA有改进,比如加入IMU或视觉惯性BA。此外,可能提到使用其他优化如Ceres,或者增量式BA、滑动窗口法来提升效率。还有边缘化策略和自适应参数调整,比如根据运动状态调整迭代次数。 用户提供的引用中,引用[1]提到了两阶段的全局-局部图优化,这可能和ORB-SLAM的全局BA与局部BA有关,可以引用这一点。引用[3]提到图形学中的算法,可能和BA中的优化算法或数值方法有关联,但不太确定是否需要引用。需要检查引用内容是否相关,再决定是否加入。 还要注意生成相关问题,可能涉及BA的作用、与其他SLAM的区别、实时性挑战、鲁棒性处理等。确保问题覆盖用户可能关心的其他方面。 需要验证提到的优化策略是否准确,比如ORB-SLAM是否真的使用Huber核函数,以及边缘化的具体应用。可能还需要确认ORB-SLAM3中的改进点,比如视觉惯性BA的具体实现方式。 最后,确保回答结构清晰,分点说明,并且引用正确。可能需要调整语言更口语化,但按照用户要求,使用中文回答,并符合格式要求。</think>### ORB-SLAM中的BA优化实现与改进方法 #### 一、ORB-SLAMBA优化实现 1. **BA的作用与层级** - **局部BA**:优化当前关键帧及其共视关键帧,以及这些关键帧观测到的地图点,减少局部累积误差[^1]。 - **全局BA**:在闭环检测后优化所有关键帧和地图点,确保全局一致性。 - 数学形式可表示为最小化重投影误差: $$ \min_{\mathbf{T}_i, \mathbf{P}_j} \sum_{i,j} \rho \left( \|\pi(\mathbf{T}_i \mathbf{P}_j) - \mathbf{u}_{ij} \|^2 \right) $$ 其中$\mathbf{T}_i$为位姿,$\mathbf{P}_j$为3D点,$\rho$为Huber鲁棒核函数。 2. **实现步骤** - **数据结构准备**:提取当前激活的关键帧集合及其观测的地图点。 - **图模型构建**: - 关键帧位姿作为SE(3)节点 - 地图点坐标作为3D点节点 - 观测关系构建一元边(单目)或二元边(双目) - **优化策略**: - 使用g2o的LM(Levenberg-Marquardt)优化器 - 固定非活跃关键帧的位姿 - 采用两级优化:先优化位姿再联合优化位姿与地图点 3. **代码级实现(简化逻辑)** ```cpp void Optimizer::LocalBundleAdjustment(KeyFrame* pKF, bool* pbStopFlag) { // 1. 获取局部关键帧集合 vector<KeyFrame*> vLocalKeyFrames = pKF->GetVectorCovisibleKeyFrames(); // 2. 构建优化器 g2o::SparseOptimizer optimizer; g2o::BlockSolver_6_3::LinearSolverType* linearSolver = new g2o::LinearSolverEigen<g2o::BlockSolver_6_3::PoseMatrixType>(); g2o::OptimizationAlgorithmLevenberg* solver = new g2o::OptimizationAlgorithmLevenberg( new g2o::BlockSolver_6_3(linearSolver)); // 3. 添加顶点与边 for(KeyFrame* pKFi : vLocalKeyFrames) { g2o::VertexSE3Expmap* vSE3 = new g2o::VertexSE3Expmap(); vSE3->setEstimate(Converter::toSE3Quat(pKFi->GetPose())); optimizer.addVertex(vSE3); // 添加地图点观测边 for(MapPoint* pMP : pKFi->GetMapPoints()) { g2o::EdgeSE3ProjectXYZ* e = new g2o::EdgeSE3ProjectXYZ(); e->setVertex(0, dynamic_cast<g2o::OptimizableGraph::Vertex*>(vSE3)); // ... 设置其他参数 optimizer.addEdge(e); } } // 4. 执行优化 optimizer.initializeOptimization(); optimizer.optimize(5); // 迭代5次 } ``` #### 二、改进方法 1. **算法层面的改进** - **视觉惯性BA**(ORB-SLAM3新增): $$ E = \sum \rho(\| \mathbf{r}_{IMU} \|_{\Sigma_{IMU}}^2 ) + \sum \rho(\| \mathbf{r}_{proj} \|_{\Sigma_{img}}^2 ) $$ 联合优化相机位姿、速度、IMU偏置和地图点[^1]。 - **滑动窗口BA**:限制优化关键帧数量,保持实时性。 - **边缘化策略**:使用Schur补消除旧关键帧的状态,保留其约束信息。 2. **工程优化** - **增量式BA**:仅优化受新观测影响的部分变量 - **多线程加速**:将雅可比矩阵计算与线性求解步骤并行化 - **自适应参数调整**:根据运动剧烈程度动态调整迭代次数 3. **鲁棒性增强** - **动态权重分配**:基于特征匹配质量调整信息矩阵权重 - **多假设BA**:维护多个优化假设,选择最一致的解 #### 三、典型改进案例 DiSCo-SLAM提出的两阶段优化框架[^1]: 1. **第一阶段**:分布式局部BA优化 - 各机器人独立优化局部地图 - 采用scan context进行快速数据关联 2. **第二阶段**:全局一致性优化 - 整合各局部优化结果 - 构建全局位姿图进行轻量级优化
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

大江东去浪淘尽千古风流人物

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值