三维重建(7)--运动恢复结构SfM系统解析

2024-01-26 06:04

本文主要是介绍三维重建(7)--运动恢复结构SfM系统解析,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

目录

一、SfM系统(两视图)

1、特征提取

2、特征匹配

3、RANSAC求解基础矩阵F 

4、完整的欧式结构恢复算法流程

二、基于增量法的SfM系统(以OpenMVG为例)

1、预处理 

2、图像特征点提取与匹配

3、两视图重构点云

4、增加新视图,多视图重构


一、SfM系统(两视图)

        对于欧式结构恢复的两视图问题,需要获得三维场景的m张图像的像坐标作为已知条件,求解三维场景结构(即三维点坐标),m个摄像机的外参数R和T。所以现在的问题在于如何标注m张图像的像坐标的对应关系,将摄像机拍摄的照片转换成已知条件。

        图像的像坐标对应关系可以转换为特征提取特征匹配两个部分。

1、特征提取

         特征提取:通过输入m张不同摄像机拍摄同一场景的图片,输出具有尺度不变性的特征点,即尺度不变特征转换(SIFT转换),

        SIFT算法特点:具有尺度不变性、能够适应旋转图像,改变图像亮度,移动拍摄位置的变化,能在一定程度上不受视角变化、仿射变换、噪声的干扰。

        

        SIFT算法实现特征提取的流程:提取尺度不变的区域(高斯微分函数),进行尺寸归一化和旋转归一化,对特征点区域计算特征方向,进行两两比较找出相互匹配的若干对特征点,建立两个像平面之间特征点的对应关系。

2、特征匹配

        对于二视图的特征匹配而言,首先使用特征提取器对两幅图像进行特征提取,提取出的特征点可能不同,对于下图,假设我们找右图特征点在左图的对应点。

        首先,选择一个右图中的特征点i,并计算左图所有特征点的一个相似度度量(比如余弦相似度度量),并选择这其中top2距离小的两个左图特征点j_1,j_2,并记录j_1,j_2与特征点i之间的距离d_1,d_2

        计算距离比d_1/d_2(指定小的数除以大的数),如果其小于某个给定阈值(如0.6),则认为右图特征点i与左图特征点j_1是一对对应点,这步算法的设计,如果说d_1/d_2较大,趋于1,则说明左图两个特征点之间较近,越说明左图两个特征点与右图特征点i建立的两组对应关系很难辨认,极为接近,所以我们不作为匹配对象,这更容易引入噪声。

        我们应该注意的是,特征匹配是为了后续求解基础矩阵F进行准备的。

3、RANSAC求解基础矩阵F 

         RANSAC求解基础矩阵的方法是对归一化八点法估计基础矩阵F的扩充,为了进一步去减少离群点和噪声,我们在求解基础矩阵时,可以使用RANSAC方法来充分排除离群点对求解基础矩F的干扰,提高估计的准确性和鲁棒性。

4、完整的欧式结构恢复算法流程

        对于欧式结构恢复问题,一般假设摄像机1与世界坐标系不存在旋转平移关系,即外参数为[I \hspace{0.2cm} 0],摄像机2与世界坐标系存在R,T关系,即外参数为[R \hspace{0.2cm}T],换言之摄像机2与摄像机1之间存在[R \hspace{0.2cm}T]关系。

        算法流程: 

(1)对应点计算(特征提取+特征匹配)

(2)求解基础矩阵F(RANSAC+归一化八点法)

(3)求解本质矩阵E

(4)分解本质矩阵E\rightarrow R,T\rightarrow M_2

(5)三角化

二、基于增量法的SfM系统(以OpenMVG为例)

        基于增量法的SfM系统流程:图像特征点提取与匹配->两视图重构初始点云->增加视图到系统实现多视图重构

        算法流程如下: 

1、预处理 

        预处理工作:图像特征点提取与近邻匹配,基于RANSAC的基础矩阵或单应矩阵的估计。

为什么要估计基础矩阵的同时估计单应矩阵?

        由于我们不能确定所取的特征点在世界坐标系下有没有很多是同一平面的点,如果同一平面的点较多,那么使用单应矩阵的估计后去计算本质矩阵E的误差要小很多,受到噪声的影响也很小。所以利用RANSAC来估计基础矩阵和单应矩阵,如果単应矩阵求解更容易拟合则使用単应矩阵。

2、图像特征点提取与匹配

(1)计算对应点的轨迹t

        对应点的轨迹:多个像平面下对应特征点的连线,在OpenMVG中只保留关联至少三张图片的特征点的连线。

 (2)计算连通图G

        连通图:这里的连通图指的是以图片作为结点,两个图片之间匹配特征数量达到某一阈值则作为一条边,未达到阈值则不构成边。

3、两视图重构点云

(1)在G中选取一条边e

        在连通图G中的多条边里,选取一条特殊的边,这条边的两个节点对应的两个图像满足:这两个图像的对应点三角化时的射线夹角中位数不大于60度且不小于3度,此时能保证这两个摄像机即不会离的太近容易同时丢失太多信息,也不会离的过远而导致对应点对偏少。

(2)后续过程

        接下来的三步:鲁棒估计本质矩阵E,分解本质矩阵E,三角化,都是两视图重构的流程,在这里不再重复介绍。       

        删除G中的边e:这一步也很好理解,由于我们使用了这个边进行两视图重构,那么后续将不再用这条边。

4、增加新视图,多视图重构

        当G中存在其他边时,进行多视图重构。

        选择一条G中的边e,满足该边的track(e)∩ {已重建3D点}最大化,即这条边所关联的两个图像的对应点对,为所有剩余边中最大的,优先选择该边。

        用PnP方法,对边e所对应两个图像中,未进行重建摄像机位姿的图像一方进行重建摄像机位姿。

        然后,三角化新的tracks,删除G中的边e,并执行捆绑调整操作。

        每一次引入新的视图,将执行一次捆绑调整操作,可以对世界坐标的三维点进行不断拟合。

Reference:https://www.bilibili.com/video/BV1DQ4y1e7x6/?spm_id_from=333.999.0.0&vd_source=2152dec20715e478285cc87cc31201ed

这篇关于三维重建(7)--运动恢复结构SfM系统解析的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/645853

相关文章

网页解析 lxml 库--实战

lxml库使用流程 lxml 是 Python 的第三方解析库,完全使用 Python 语言编写,它对 XPath表达式提供了良好的支 持,因此能够了高效地解析 HTML/XML 文档。本节讲解如何通过 lxml 库解析 HTML 文档。 pip install lxml lxm| 库提供了一个 etree 模块,该模块专门用来解析 HTML/XML 文档,下面来介绍一下 lxml 库

不懂推荐算法也能设计推荐系统

本文以商业化应用推荐为例,告诉我们不懂推荐算法的产品,也能从产品侧出发, 设计出一款不错的推荐系统。 相信很多新手产品,看到算法二字,多是懵圈的。 什么排序算法、最短路径等都是相对传统的算法(注:传统是指科班出身的产品都会接触过)。但对于推荐算法,多数产品对着网上搜到的资源,都会无从下手。特别当某些推荐算法 和 “AI”扯上关系后,更是加大了理解的难度。 但,不了解推荐算法,就无法做推荐系

基于人工智能的图像分类系统

目录 引言项目背景环境准备 硬件要求软件安装与配置系统设计 系统架构关键技术代码示例 数据预处理模型训练模型预测应用场景结论 1. 引言 图像分类是计算机视觉中的一个重要任务,目标是自动识别图像中的对象类别。通过卷积神经网络(CNN)等深度学习技术,我们可以构建高效的图像分类系统,广泛应用于自动驾驶、医疗影像诊断、监控分析等领域。本文将介绍如何构建一个基于人工智能的图像分类系统,包括环境

水位雨量在线监测系统概述及应用介绍

在当今社会,随着科技的飞速发展,各种智能监测系统已成为保障公共安全、促进资源管理和环境保护的重要工具。其中,水位雨量在线监测系统作为自然灾害预警、水资源管理及水利工程运行的关键技术,其重要性不言而喻。 一、水位雨量在线监测系统的基本原理 水位雨量在线监测系统主要由数据采集单元、数据传输网络、数据处理中心及用户终端四大部分构成,形成了一个完整的闭环系统。 数据采集单元:这是系统的“眼睛”,

使用SecondaryNameNode恢复NameNode的数据

1)需求: NameNode进程挂了并且存储的数据也丢失了,如何恢复NameNode 此种方式恢复的数据可能存在小部分数据的丢失。 2)故障模拟 (1)kill -9 NameNode进程 [lytfly@hadoop102 current]$ kill -9 19886 (2)删除NameNode存储的数据(/opt/module/hadoop-3.1.4/data/tmp/dfs/na

电脑桌面文件删除了怎么找回来?别急,快速恢复攻略在此

在日常使用电脑的过程中,我们经常会遇到这样的情况:一不小心,桌面上的某个重要文件被删除了。这时,大多数人可能会感到惊慌失措,不知所措。 其实,不必过于担心,因为有很多方法可以帮助我们找回被删除的桌面文件。下面,就让我们一起来了解一下这些恢复桌面文件的方法吧。 一、使用撤销操作 如果我们刚刚删除了桌面上的文件,并且还没有进行其他操作,那么可以尝试使用撤销操作来恢复文件。在键盘上同时按下“C

嵌入式QT开发:构建高效智能的嵌入式系统

摘要: 本文深入探讨了嵌入式 QT 相关的各个方面。从 QT 框架的基础架构和核心概念出发,详细阐述了其在嵌入式环境中的优势与特点。文中分析了嵌入式 QT 的开发环境搭建过程,包括交叉编译工具链的配置等关键步骤。进一步探讨了嵌入式 QT 的界面设计与开发,涵盖了从基本控件的使用到复杂界面布局的构建。同时也深入研究了信号与槽机制在嵌入式系统中的应用,以及嵌入式 QT 与硬件设备的交互,包括输入输出设

JAVA智听未来一站式有声阅读平台听书系统小程序源码

智听未来,一站式有声阅读平台听书系统 🌟 开篇:遇见未来,从“智听”开始 在这个快节奏的时代,你是否渴望在忙碌的间隙,找到一片属于自己的宁静角落?是否梦想着能随时随地,沉浸在知识的海洋,或是故事的奇幻世界里?今天,就让我带你一起探索“智听未来”——这一站式有声阅读平台听书系统,它正悄悄改变着我们的阅读方式,让未来触手可及! 📚 第一站:海量资源,应有尽有 走进“智听

【C++】_list常用方法解析及模拟实现

相信自己的力量,只要对自己始终保持信心,尽自己最大努力去完成任何事,就算事情最终结果是失败了,努力了也不留遗憾。💓💓💓 目录   ✨说在前面 🍋知识点一:什么是list? •🌰1.list的定义 •🌰2.list的基本特性 •🌰3.常用接口介绍 🍋知识点二:list常用接口 •🌰1.默认成员函数 🔥构造函数(⭐) 🔥析构函数 •🌰2.list对象

usaco 1.3 Mixing Milk (结构体排序 qsort) and hdu 2020(sort)

到了这题学会了结构体排序 于是回去修改了 1.2 milking cows 的算法~ 结构体排序核心: 1.结构体定义 struct Milk{int price;int milks;}milk[5000]; 2.自定义的比较函数,若返回值为正,qsort 函数判定a>b ;为负,a<b;为0,a==b; int milkcmp(const void *va,c