CN109215118B - 一种基于图像序列的增量式运动结构恢复优化方法 - Google Patents
一种基于图像序列的增量式运动结构恢复优化方法 Download PDFInfo
- Publication number
- CN109215118B CN109215118B CN201811086146.7A CN201811086146A CN109215118B CN 109215118 B CN109215118 B CN 109215118B CN 201811086146 A CN201811086146 A CN 201811086146A CN 109215118 B CN109215118 B CN 109215118B
- Authority
- CN
- China
- Prior art keywords
- image
- images
- point
- points
- point cloud
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种基于图像序列的增量式运动结构恢复优化方法,涉及三维重建领域。包括:获取所需重建物体各个角度的图像;获取两两图像之间特征点的关系及描述符;过滤特征点的异常值;排除异解;形成初始点云模型;丰富点云模型;得到最终的点云模型;对点云模型进行优化。本方法在初始图像对选择过程中使用自适应阈值估计法,使其能够自动的调整过滤异常特征点的值,并引入约束进一步过滤,获得更精确的初始图像对,确保能够进行鲁棒的初始重建,提高整体三维重建质量;在增量添加迭代过程中使用EPNP法求解,提高迭代的解算效率及精度,从而提高本方法的性能;通过BA法,进一步提高了模型优化的精确度,解决了增量迭代过程中误差累计导致的问题。
Description
技术领域
本发明涉及三维重建技术领域,具体涉及一种基于图像序列的增量式运动结构恢复优化方法。
背景技术
随着计算机相关技术的不断进步,虚拟现实技术得到了迅猛的发展。在虚拟场景的构建过程中,大量的三维模型将被使用。而传统工作流程中,三维模型的建立大多依靠繁琐的手工建模软件,大大延长了开发周期,同时显著提高了开发成本。为此,我们寻求数据源更易获取且更为快速的三维建模方法。
图像,是三维空间的二维映射,其中包含着三维空间中相机位置,焦距等信息,利用已知的图像中的这些信息恢复为三维立体结构,可以看做是相机成像的逆过程。在一组对同一物体不同角度拍摄的图像序列中,图像与图像之间的一些特征点存在着对应关系,利用现有的方法
运动结构恢复算法,是实现三维重建中稀疏重建步骤的关键算法。其中增量式运动结构算法是目前国内外研究的主流算法。它通过几何、摄影测量、摄影成像等知识,以增量迭代为核心思路,从图像序列中由初始重建,不断添加,将二维图像上的像素点通过图像之间的联系还原为三维点云中的结构。传统增量式运动结构恢复算法存在着效率和重建质量等多方面的问题,对其的优化是一个重要的研究方向。
发明内容
针对现有技术存在的问题,本发明提供一种基于图像序列的增量式运动结构恢复优化方法,可以加快增量式运动结构恢复算法的计算效率,并提高生成点云模型的质量。
为了实现上述目的,一种基于图像序列的增量式运动结构恢复优化方法,包括以下步骤:
步骤1:获取所需重建物体各个角度的图像;
步骤2:利用SFIT特征点提取和匹配方法对图像进行处理,得到两两图像之间特征点的关系及描述符;
步骤3:使用反向理论研究模型自适应估计阈值对图像之间特征点的异常值即外点进行过滤,并通过约束对图像之间特征点的异常值进行进一步过滤;
步骤4:采用改进的Nister的五点法,利用点到极线距离最小原理排除异解;
步骤5:选择符合约束且经过异常值过滤后异常值占比最低的两幅图像作为初始图像对,并对初始图像中匹配的点进行初始重建,形成初始点云模型,具体步骤如下:
步骤5.1:选择符合约束且经过异常值过滤后异常值占比最低的两幅图像作为初始图像对;
步骤5.2:计算初始化图像对的本质矩阵E;
步骤5.3:使用奇异值分解法求解旋转矩阵R和平移向量T;
步骤5.4:根据本质矩阵E、旋转矩阵R以及平移向量T确定初始图像对的相机位置;
步骤5.5:根据相机位置,利用双视图重建原理进行初始重建;
步骤6:加入新的图像,根据新增图像与点云模型中已有的图像的轨迹关系,计算新增图像的旋转矩阵和平移向量,继续对平面图像中匹配的二维点进行三角化从而计算出三维空间中的位置点,并添加到点云模型中;
步骤7:重复步骤6,直到所有图像中的所有已匹配特征点均已被还原为三维空间中的位置点,得到最终的点云模型;
步骤8:基于捆绑调整算法对点云模型进行优化,具体步骤如下:
步骤8.1:使用全局捆绑调整策略,引入Huber函数作为损失函数,改良捆绑调整模型,并利用改良后的捆绑调整模型对计算出的所有三维空间中的位置点进行修正;
步骤8.2:进行图形完整性检查,在所有的观察射线对中引入最小三角测量角并对轨迹进行重三角测量;
步骤8.3:重复步骤7.1至步骤7.2两次,得到最终重建的的点云模型。
进一步地,所述步骤3中反向理论研究模型自适应估计阈值的计算公式如下:
式中,NFA(M,k)为反向理论研究模型自适应估计阈值,M为输入的图像数据,k为假定的内点联系的数量,Nout为可从RANSAC对应样本中采样估计的模型数目,n为总的关系数量,Nsample为RANSAC样本的主体,ek(M)为输入的图像数据M在n个对应关系中的第k个最小误差,d为误差维度,α0为具有1像素误差的随机对应关系的概率。
进一步地,所述步骤3中通过约束对异常值进行进一步过滤,包括如下约束:
约束1:使用五点法核验内点,将通过核验的图像根据内点数量进行排序,排在前n%的图像优先考虑选择为可以用来作为初始图像对的候选图像;
约束2:选择宽基线作为约束条件,弥补短基线匹配技术在深度信息估计方面的不足对重建效果的影响;
所述宽基线约束具体是:两幅图像经特征提取和匹配之后所形成的描述相机位置关系的一条直线,它可以弥补短基线在深度信息还原方面不足所带来对重建效果的影响。
进一步地,所述步骤4中利用点到极线距离最小原理排除异解的具体原理如下:
正确的相机姿态和位置元素对应的本质矩阵经过SVD奇异值分解计算,将第一幅图像上的点投影到第二幅图像上时将该点变为一条极线,第二幅图像上的对应点到该极线的距离最小,第二幅图像到第一幅图像同理,如果是误解,这两个距离的和则不是最小的;
所述相机姿态使用3x3的旋转矩阵R与1x3的平移向量T共同描述,由计算出的本质矩阵E分解得到。
本发明的有益效果:
本发明提出一种基于图像序列的增量式运动结构恢复优化方法,在初始图像对选择过程中使用自适应阈值估计方法,使算法能够自动的根据图像数据的状态调整过滤异常特征点(外点)的值,并引入约束进一步过滤,从而可以获得更精确的初始图像对,以便于确保能够进行鲁棒的初始重建,能够为最终提高整体三维重建质量打下基础;在增量添加迭代过程中使用高效n点透视法(EPNP)方法求解,能有效的提高增量迭代过程的解算效率,提高解算的精度,在一定程度上提高了整个方法的性能;通过优化捆绑调整(BA)方法,进一步提高了模型优化的精确度,解决了增量迭代过程中误差累计导致的问题。
附图说明
图1为本发明实施例中基于图像序列的增量式运动结构恢复优化方法的流程图。
具体实施方式
为了使本发明的目的、技术方案及优势更加清晰,下面结合附图和具体实施例对本发明做进一步详细说明。此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
一种基于图像序列的增量式运动结构恢复优化方法,流程如图1所示,包括以下步骤:
步骤1:获取所需重建物体各个角度的图像;
本实施例中,通过数码相机、手机或无人机等其他设备,对所需重建的物体进行拍摄,获取物体各个角度的图像。
步骤2:利用SFIT特征点提取和匹配方法对图像进行处理,得到两两图像之间特征点的关系及描述符;
步骤3:使用反向理论研究模型自适应估计阈值对图像之间特征点的异常值即外点进行过滤,并通过约束对图像之间特征点的异常值进行进一步过滤;
所述反向理论研究模型自适应估计阈值的计算公式如下:
式中,NFA(M,k)为反向理论研究模型自适应估计阈值,M为输入的图像数据,k为假定的内点联系的数量,Nout为可从RANSAC对应样本中采样估计的模型数目,n为总的关系数量,Nsample为RANSAC样本的主体,ek(M)为输入的图像数据M在n个对应关系中的第k个最小误差,d为误差维度,α0为具有1像素误差的随机对应关系的概率。
本实施例中,为了求解上式,还需要引入关于模型M的约束条件,公式如下:
式中,ε通常设为1,结合自适应估计阈值计算公式、约束条件和已知参数,即可计算出自适应估计阈值。
所述通过约束对异常值进行进一步过滤,包括如下约束:
约束1:使用五点法核验内点,将通过核验的图像根据内点数量进行排序,排在前n%的图像优先考虑选择为可以用来作为初始图像对的候选图像;
约束2:选择宽基线作为约束条件,弥补短基线匹配技术在深度信息估计方面的不足对重建效果的影响;
所述宽基线约束具体是:两幅图像经特征提取和匹配之后所形成的描述相机位置关系的一条直线,它可以弥补短基线在深度信息还原方面不足所带来对重建效果的影响。
步骤4:采用改进的Nister的五点法,利用点到极线距离最小原理排除异解;
所述利用点到极线距离最小原理排除异解的具体原理如下:
正确的相机姿态和位置元素对应的本质矩阵经过SVD奇异值分解计算,将第一幅图像上的点投影到第二幅图像上时将该点变为一条极线,第二幅图像上的对应点到该极线的距离最小,第二幅图像到第一幅图像同理,如果是误解,这两个距离的和则不是最小的;
所述相机姿态使用3x3的旋转矩阵R与1x3的平移向量T共同描述,由计算出的本质矩阵E分解得到。
本实施例中,所述利用点到极线距离最小原理排除异解的公式如下:
式中,xi为第一幅图像中第i个匹配点,xi’为第二幅图像中与第一幅图像中第i个匹配点对应的匹配点,N为图像中匹配点的总个数,E为本质矩阵,T为转置符号,d(*,*)2表示两点之间的距离。
步骤5:选择符合约束且经过异常值过滤后异常值占比最低的两幅图像作为初始图像对,并对初始图像中匹配的点进行初始重建,形成初始点云模型,具体步骤如下:
步骤5.1:选择符合约束且经过异常值过滤后异常值占比最低的两幅图像作为初始图像对:
步骤5.2:计算初始化图像对的本质矩阵E;
步骤5.3:使用奇异值分解法求解旋转矩阵R和平移向量T;
步骤5.4:根据本质矩阵E、旋转矩阵R以及平移向量T确定初始图像对的相机位置;
步骤5.5:根据相机位置,利用双视图重建原理进行初始重建;
步骤6:加入新的图像,根据新增图像与点云模型中已有的图像的轨迹关系,计算新增图像的旋转矩阵和平移向量,继续对平面图像中匹配的二维点进行三角化从而计算出三维空间中的位置点,并添加到点云模型中;
本实施例中,所述三角化的公式如下:
式中,表示未知的点坐标,{(wi)}i=1n表示未知点的第i个投影参数,表示已知的第i个二维点在图像平面上的坐标,其中在第三行添加1是为了便于进行矩阵计算,fu和fv均表示相机焦距参数,uc和vc表示相机主点相对二维图像平面的位置,αij表示灰度质心齐次坐标,由公式的第三行可知,根据未知和已知参数,该等式可分解为两个线性系统并进行求解,公式如下:
根据如上公式结合点云模型中已有的三维空间位置点,即可求解新添加入点云模型中二维平面图像上未被还原为三维空间中位置的点。
步骤7:重复步骤6,直到所有图像中的所有已匹配特征点均已被还原为三维空间中的位置点,得到最终的点云模型;
步骤8:基于捆绑调整算法对点云模型进行优化,具体步骤如下:
步骤8.1:使用全局捆绑调整策略,引入Huber函数作为损失函数,改良捆绑调整模型,并利用改良后的捆绑调整模型对计算出的所有三维空间中的位置点进行修正;
本实施例中,所述对计算出的所有三维空间中的位置点进行修正的公式如下:
式中,δij表示相机i是否观察到场景点j,若δij=1表示相机i观察到场景点j,若δij=0表示相机i没有观察到场景点j,Ki表示相机i的内参矩阵,Ri表示相机i的旋转矩阵,Ci表示相机i的中心,γ(Ki,Ri,Ci,Xj)表示投影函数,xij表示确定的二维图像点位置。
步骤8.2:进行图形完整性检查,在所有的观察射线对中引入最小三角测量角并对轨迹进行重三角测量;
本实施例中,所述三角测量角公式如下:
式中,Xab是已三角化的点。
经过本步骤可以过滤具有大的重投影误差的轨迹,更好的提高计算重建点的深度的准确性,降低误差累计并增强重建的完整性,使计算值更贴近于真实值。
步骤8.3:重复步骤7.1至步骤7.2两次,得到最终重建的的点云模型。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解;其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;因而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (6)
1.一种基于图像序列的增量式运动结构恢复优化方法,其特征在于,包括以下步骤:
步骤1:获取所需重建物体各个角度的图像;
步骤2:利用SFIT特征点提取和匹配方法对图像进行处理,得到两两图像之间特征点的关系及描述符;
步骤3:使用反向理论研究模型自适应估计阈值对图像之间特征点的异常值即外点进行过滤,并通过约束对图像之间特征点的异常值进行进一步过滤;
步骤4:采用改进的Nister的五点法,利用点到极线距离最小原理排除异解;
步骤5:选择符合约束且经过异常值过滤后异常值占比最低的两幅图像作为初始图像对,并对初始图像中匹配的点进行初始重建,形成初始点云模型;
步骤6:加入新的图像,根据新增图像与点云模型中已有的图像的轨迹关系,计算新增图像的旋转矩阵和平移向量,继续对平面图像中匹配的二维点进行三角化从而计算出三维空间中的位置点,并添加到点云模型中;
步骤7:重复步骤6,直到所有图像中的所有已匹配特征点均已被还原为三维空间中的位置点,得到最终的点云模型;
步骤8:基于捆绑调整算法对点云模型进行优化。
3.根据权利要求1所述的基于图像序列的增量式运动结构恢复优化方法,其特征在于,所述步骤3中通过约束对异常值进行进一步过滤,包括如下约束:
约束1:使用五点法核验内点,将通过核验的图像根据内点数量进行排序,排在前n%的图像优先考虑选择为可以用来作为初始图像对的候选图像;
约束2:选择宽基线作为约束条件,弥补短基线匹配技术在深度信息估计方面的不足对重建效果的影响;
所述宽基线约束具体是:两幅图像经特征提取和匹配之后所形成的描述相机位置关系的一条直线,它可以弥补短基线在深度信息还原方面不足所带来对重建效果的影响。
4.根据权利要求1所述的基于图像序列的增量式运动结构恢复优化方法,其特征在于,所述步骤4中利用点到极线距离最小原理排除异解的具体原理如下:
正确的相机姿态和位置元素对应的本质矩阵经过SVD奇异值分解计算,将第一幅图像上的点投影到第二幅图像上时将该点变为一条极线,第二幅图像上的对应点到该极线的距离最小,第二幅图像到第一幅图像同理,如果是误解,这两个距离的和则不是最小的;
所述相机姿态使用3x3的旋转矩阵R与1x3的平移向量T共同描述,由计算出的本质矩阵E分解得到。
5.根据权利要求1所述的基于图像序列的增量式运动结构恢复优化方法,其特征在于,所述步骤5包括如下步骤:
步骤5.1:选择符合约束且经过异常值过滤后异常值占比最低的两幅图像作为初始图像对;
步骤5.2:计算初始化图像对的本质矩阵E;
步骤5.3:使用奇异值分解法求解旋转矩阵R和平移向量T;
步骤5.4:根据本质矩阵E、旋转矩阵R以及平移向量T确定初始图像对的相机位置;
步骤5.5:根据相机位置,利用双视图重建原理进行初始重建。
6.根据权利要求1所述的基于图像序列的增量式运动结构恢复优化方法,其特征在于,所述步骤8包括如下步骤:
步骤8.1:使用全局捆绑调整策略,引入Huber函数作为损失函数,改良捆绑调整模型,并利用改良后的捆绑调整模型对计算出的所有三维空间中的位置点进行修正;
步骤8.2:进行图形完整性检查,在所有的观察射线对中引入最小三角测量角并对轨迹进行重三角测量;
步骤8.3:重复步骤7.1至步骤7.2两次,得到最终重建的的点云模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811086146.7A CN109215118B (zh) | 2018-09-18 | 2018-09-18 | 一种基于图像序列的增量式运动结构恢复优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811086146.7A CN109215118B (zh) | 2018-09-18 | 2018-09-18 | 一种基于图像序列的增量式运动结构恢复优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109215118A CN109215118A (zh) | 2019-01-15 |
CN109215118B true CN109215118B (zh) | 2022-11-29 |
Family
ID=64984562
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811086146.7A Active CN109215118B (zh) | 2018-09-18 | 2018-09-18 | 一种基于图像序列的增量式运动结构恢复优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109215118B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11166022B2 (en) * | 2019-06-04 | 2021-11-02 | Google Llc | Quantization constrained neural image coding |
CN111144483B (zh) * | 2019-12-26 | 2023-10-17 | 歌尔股份有限公司 | 一种图像特征点过滤方法以及终端 |
CN111815765B (zh) * | 2020-07-21 | 2022-07-05 | 西北工业大学 | 一种基于异构数据融合的图像三维重建方法 |
CN114003190B (zh) * | 2021-12-30 | 2022-04-01 | 江苏移动信息系统集成有限公司 | 一种适应多场景和多设备的增强现实方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012044217A1 (en) * | 2010-10-01 | 2012-04-05 | Saab Ab | Method and apparatus for optimization and incremental improvement of a fundamental matrix |
CN106600686A (zh) * | 2016-12-06 | 2017-04-26 | 西安电子科技大学 | 一种基于多幅未标定图像的三维点云重建方法 |
WO2017219391A1 (zh) * | 2016-06-24 | 2017-12-28 | 深圳市唯特视科技有限公司 | 一种基于三维数据的人脸识别系统 |
-
2018
- 2018-09-18 CN CN201811086146.7A patent/CN109215118B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012044217A1 (en) * | 2010-10-01 | 2012-04-05 | Saab Ab | Method and apparatus for optimization and incremental improvement of a fundamental matrix |
WO2017219391A1 (zh) * | 2016-06-24 | 2017-12-28 | 深圳市唯特视科技有限公司 | 一种基于三维数据的人脸识别系统 |
CN106600686A (zh) * | 2016-12-06 | 2017-04-26 | 西安电子科技大学 | 一种基于多幅未标定图像的三维点云重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109215118A (zh) | 2019-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107945220B (zh) | 一种基于双目视觉的重建方法 | |
CN110363858B (zh) | 一种三维人脸重建方法及系统 | |
CN109215118B (zh) | 一种基于图像序列的增量式运动结构恢复优化方法 | |
CN111462329B (zh) | 一种基于深度学习的无人机航拍影像的三维重建方法 | |
CN112102458A (zh) | 基于激光雷达点云数据辅助的单镜头三维图像重构方法 | |
CN114399554B (zh) | 一种多相机系统的标定方法及系统 | |
CN107507277B (zh) | 三维点云重构方法和装置、服务器及可读存储介质 | |
CN112381847B (zh) | 管路端头空间位姿测量方法及系统 | |
CN111160298A (zh) | 一种机器人及其位姿估计方法和装置 | |
CN112630469B (zh) | 一种基于结构光和多光场相机的三维检测方法 | |
CN113838191A (zh) | 一种基于注意力机制和单目多视角的三维重建方法 | |
CN113160335A (zh) | 一种基于双目视觉的模型点云及三维表面重建方法 | |
CN113313740B (zh) | 一种基于平面连续性的视差图和表面法向量联合学习方法 | |
CN114119987A (zh) | 基于卷积神经网络的特征提取和描述子生成方法及系统 | |
JP2009186287A (ja) | 平面パラメータ推定装置、平面パラメータ推定方法及び平面パラメータ推定プログラム | |
KR100965838B1 (ko) | 항공용 라이다로부터 추출된 빌딩폴리곤의 내재적기하정규화 방법 | |
Guan et al. | Efficient recovery of multi-camera motion from two affine correspondences | |
KR101166719B1 (ko) | 평행사변형 기반의 무한 호모그래피 산출 방법과 이를 이용한 건물 구조 복원 방법 | |
CN110570473A (zh) | 一种基于点线融合的权重自适应位姿估计方法 | |
CN114399547B (zh) | 一种基于多帧的单目slam鲁棒初始化方法 | |
CN113486928B (zh) | 一种基于有理多项式模型可微分张量表达的多视影像对齐方法 | |
CN114998532A (zh) | 一种基于数字图像重建的三维影像视觉传达优化方法 | |
JP2019032660A (ja) | 撮像システム、撮像方法 | |
Yang et al. | 3D reconstruction through measure based image selection | |
CN114004848A (zh) | 一种多视角遥感地形影像点云重建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |