摘要
大家好，今天为大家带来的文章是 OpenREALM: Real-time Mapping for Unmanned Aerial Vehicles
本文提出了一种基于无人机(UAV)的实时建图框架OpenREALM。利用无人机的机载计算机上的照机来获取目标感兴趣区域的高分辨率图像。不同的操作模式允许OpenREALM在近似平面地面下执行简单的拼接，或者完全恢复复杂的三维表面信息，以提取高程图和几何校正的正射影像。此外，还利用无人机的全局位置对数据进行地理推理。在所有模式下，地面操作员都可以现场查看结果地图的建图进度。获得的最新的表面信息将推动各种无人机应用。
Github：https://github.com/laxnpander/OpenREALM
文章：https://arxiv.org/abs/2009.10492v1
主要工作与贡献
目前大多数框架只在无人机飞行过后通过捕获的所有图像，进行建图和测绘，这是显然不是最优的。在任务期间根据所获得的数据进行干预是不可能的，尽管理论上所有的数据都是可用的。从任务开始到最终确定全地图的总时间可能需要长达几个小时，在搜索和救援场景的情况下，这是不可接受的。
因此，我们提出了一个框架，它能够读取无人机稳定摄像机获得的连续图像流，并且它提供了三种不同的操作模式：
1)基于GNSS定位和实时图像拼接生成2D地图；
2)基于GNSS定位和实时图像生成2D地图和视觉姿态估计；
3)基于GNSS定位和GPU加速立体重建的视觉姿态估计生成2.5D数字高程模型。
模式（1)被设计用于现成的硬件，如大疆产品。计算出的地图仅依赖于GNSS和标题信息来将图像对齐成一个全局地图。然而，它也适合于创建一个无特征或高度动态的区域的全局地图。模式（2)和（3）需要一个具有高帧率的相机，适合于视觉姿态估计。模式（3)中的输出包括摄像机运动的估计，每个视图的深度图，高程图（正射影像）以及观测场景的密集三维点云图。
(a)结果是一个实时更新的经过几何校正的正射影像图
(b)地面表示为由三维点云数据生成的高程图。利用概率方法识别具有高噪声的数据区域并删除。
算法流程
1.相关工作
算法ROS框架
图2.显示了三个不同层的总体设计。
第一层是通过机器人操作系统(ROS)实现的“传输层”。在航空映射场景中，至少有两种不同的处理设备，无人机的相机和地面站组成，所捕获的数据由带GNSS信息的图像组成。
第二层是“阶段层”，其中包含过程控制算法。每个阶段都被封装起来，交换只通过其相应的ROS节点在一个方向上产生。这样，一个阶段的输出就会成为以下阶段的输入。
姿态估计
姿态估计阶段的工作流程如图3所示。开始时，它为支持的框架创建一个可视化SLAM接口。这目前可以是ORBSLAM2、DSO或SVO，这取决于传递的初始参数，输入帧由接口到实际的SLAM实现从而估计摄像机姿态矩阵M。M被定义为从照相机到世界帧的转换。如果跟踪成功，该模块试图利用视觉和GNSS位置来识别从局部、视觉到全局、地理坐标系统的转换。由于一组初始测量是必要的，传入帧可能会等待，直到估计误差低于一定的阈值。在求解任意尺度并对齐视觉和GNSS轨迹后，所有的帧都被发布。
如果视觉SLAM框架由于水或平面场而无法跟踪当前帧，状态切换为“丢失”，不设置视觉姿势。为了避免在此场景中映射过程完全失败，可以使用计算备用解决方案：
φ为无人机的自带的姿态,t为GNSS提供的无人机坐标,
只有当相机面对地面时，才能假设这种替代姿势是稳定的。
稠密建图
在前一阶段，当前输入帧的摄像机姿态是在一个地理坐标系中计算的。这个姿态可以直观地估计，也只能基于GNSS和视觉信息。前者具有较高的精度，但在无特征区域内缺乏鲁棒性。另一方面，后者通常总是可计算的，但不确定性很。在致密化阶段，根据图3中的工作流程，只使用具有视觉估计姿态的帧来重建观测表面密集的三维点云。
首先，检查输入帧的适用性。如果将姿态识别为视觉估计，则初始化深度地图创建。根据所选择的SLAM端口，将一组帧传递给重建接口。这反过来又提供了，类似于在姿态估计阶段的可视化SLAM界面，集成最先进的重建框架的可能性。目前只合并了PSL。在稠密重建后，将深度图投影到三维点云中，并覆盖任何先前存在的稀疏点。
纹理生成
在前面的阶段中，对所有处理过的帧估计了一个地理参考的姿态。这可以是基于使用视觉SLAM的图像数据，或基于先验的GNSS位置和姿态。然后，只对那些具有视觉估计的姿态深度图的帧生成并投影到稠密的点云中。因此，后续阶段必须处理可能发生的三种不同类型的帧：
1）GNSS位置帧、固定姿态、无点云（如果视觉姿态估计失败）帧；
2）视觉估计帧、精确姿态和稀疏点云（仅稠密失败）；
3）视觉估计帧、精确姿态和密集点云帧（如果之前所有阶段都成功的话）；
曲面生成阶段对当前输入帧的所有数据进行了评估，并提出了一种数字曲面模型。并用一个简单的平面或一个2.5D高程图来描述观察到的场景。但人们关注的不仅是海拔高度，表面法线和观测角度通常也可以在同一过程中重建。为此目的，一个有效的结构是有益的，因为它将一个特定的地理位置与收集到的所有信息紧密结合起来。欣兹曼等人提出了由范克豪瑟[18]开发的开源库“网格地图”。它是由一个感兴趣的区域和一个地面采样距离(GSD)来定义的。可以堆叠多层数据，使网格中的每个单元都由一个多维的信息向量组成。对于OpenREALM，我们采用了这个想法，但重新实现了一些模块。
图2 纹理生成算法伪代码
正射校正
正射校正阶段的目标是利用先前估计的表面模型和相机姿态来修正由视角和表面结构引起的图像的视觉失真。在最好的情况下，由此产生的正射影是高分辨率的，所以感兴趣的点（例如人类，汽车，…）可以很容易地检测到。Hinzmann等人提出了两种不同的方法来实现这种校正：
1)基于点云的正投影（正投影）
2）基于网格的正投影（向后投影）
方法（1)的计算时间最低，对重建的稠密点云有很强的依赖性。特别是由于点云上的孔洞而造成的区域覆盖范围较小。方法（2)通过将观测场景的高程保存为2.5D网格图，可以将其视为具有浮点数据而不是强度值的常规单通道图像。因此，它也可以有效地调整大小到任何必要的分辨率，就像一个常规的图像一样。总之，空间分辨率和纹理分辨率可以作为两个独立的参数来处理。空间分辨率主要取决于致密化和表面生成阶段，纹理是独立的，可以设置为一个仅受原始图像分辨率限制的选择值。如果映射的输入图像明显大于多视图重建所能处理的图像，这就特别有用。因此，为了实现，选择了“基于网格的正射校正”技术。
图4.在“正射校正”期间，初始网格图被调整到最大期望的GSD。
图像与点云拼接
通过从正位校正阶段接收第一帧，全局映射被初始化。之后，新的帧被称为“地图更新”，可以被划分为没有先验信息的区域或重叠的区域。前者被直接写入到全局地图中。然而，后者的提取为两个子地图，都只描述各自数据的重叠（全局地图重叠区域或地图更新重叠区域）。在下一步中，将为子映射的每个网格元素计算一个最能描述每个层表面的混合值。这个混合区域最终被写成了全拼接。
对于网格单元的混合，可以应用各种策略。本文选择了一种概率方法。潜在的问题陈述可以总结为：“如果一个网格单元的高度存在两个不同的假设，那么选择哪一个？”。为了决定这个基本问题，在全局地图上又添加了三层，即“海拔方差”、“海拔假说”和“观测数”。一旦地图更新到达，对于重叠区域的每个单元，将计算新高程的临时浮动平均值。
图5.所获取的自定义数据集的概述。在A、B和C中显示了无人机视角图像，由于缺乏或模糊，对视觉算法具有挑战性。
实验结果
1.视觉SLAM测姿的准确度
（a）绘制了ORBSLAM2（蓝色）和
Agisoft Metashape（虚线）的轨迹。整体对齐效果很好，特别是在x轴和y轴上。差异最大的是za轴，它表示重建的深度。
（b）在整个数据集上，ORBSLAM2方法的相对和绝对姿态误差(RPE/APE)。
首先，应该对姿态估计进行评估。高姿态精度是很重要的，因为所有连续的处理阶段都将继承不确定性。作为地面事实，我们使用经典的离线摄影测量软件Agisoft Metashape计算了所有3276张图像的轨迹。作为可视化的SLAM框架，我们将分析限制在ORBSLAM2上。结果如图所示。在a)中，显示了每个轴的平移误差。特别是在xy对齐时，几乎没有可见的位移。另一方面，Z轴等于估计的场景深度，它是使用多个视图几何图形重建的。因此，这个轴应该是最不确定的轴，而且确实可以测量到对地面真相的轻微偏差。在60-100秒之间，无人机的第一个转折点发生，这也是图7中的视觉障碍B。该区域的特征数量较少，明显降低了整体精度。图8b)中的相对姿态误差(RPE)支持这一观点，因为它在转弯前有最高的峰值，并且在转弯之后迅速下降。绝对位姿误差(APE)为0.53m，整体性能良好。
2.纹理表面重构质量
图6.通过OpenREALM重建的网格显示在左侧。为了评价地面真相，用Agisoft计算，并对得到的表面偏差进行颜色编码。边界区域中的红色区域是数据缺失的结果。
我们分析了三维表面重建问题，它被用来生成地面纹理。为此目的，数据集的前端重叠被减少到80%，因此计算时间保持合理(1：45h)。选择了各个阶段的设置，这样处理仍然在硬件上运行，但空间分辨率尽可能高。最终地图中产生的密集云被导出到开源软件“CloudCompare”中。在那里，使用迭代最近点算法(ICP)将其与地面真值对齐，以减少地理参考的影响。然后，将实时创建的密集云的每个点投影到参考网格上。生成的距离由颜色值编码，并显示在其在网格上的各自位置。如上图所示。
使用PSL对提出的实现进行稠密处理。获得的GSD为0.15m/cell。图5c所示的数据集中的混叠效应似乎对重建过程有严重影响。这所房子在3D图像中几乎无法被识别出来。这并不奇怪，因为PSL试图计算每个像素的深度值，尽管三角测量非常不确定。跑道与地面情况很吻合。在地图左上角的平面场，就在转折点之前，有最高的偏差(-1.5到0.5m)。这也是姿态估计显示出显著偏移的区域。因此，第一阶段的初始误差很可能已经传播到交叉处。总而言之，结果是可以改进的，但重要的是要记住，OpenREALM的处理时间是实时的，相比之下，就算快的软件也需要花几个小时。
图7.OpenREALM的不同操作模式与Agisoftmetashape提供的地面真值（左）相比。第一个左图使用ORBSLAM2进行对比视觉姿态估计。然而，没有进行三维重建。右图显示了本文方法，它提供了一个2.5D高程图以及真正的矫形影像。。最后一个图显示了纯粹基于无人机的GNSS位置和航向对齐获得的图像所生成的地图。
3.总结
我们提出了一种无人机的实时测绘框架。不同的操作模式使用户能够执行基于GNSS或视觉SLAM的图像拼接，或者完全重建三维表面并提取几何校正的正射影像。今后工程的表面质量有待进一步提高。目前，SLAM和三维重建的实现并没有为航空测绘场景而设计和优化。因此，在这一领域的进一步研究应大大提高研究结果。此外，为基准测试建图框架提供一个公共数据集可能会推动未来的进展。
转自：测绘学术资讯
文章来源于泡泡机器人SLAM
如有侵权，请联系本站删除！
