CN110111391B - 基于改进的下逼近非负矩阵对极几何估计方法 - Google Patents
基于改进的下逼近非负矩阵对极几何估计方法 Download PDFInfo
- Publication number
- CN110111391B CN110111391B CN201910418346.6A CN201910418346A CN110111391B CN 110111391 B CN110111391 B CN 110111391B CN 201910418346 A CN201910418346 A CN 201910418346A CN 110111391 B CN110111391 B CN 110111391B
- Authority
- CN
- China
- Prior art keywords
- model
- matrix
- negative matrix
- point
- hypothesis
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/2433—Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
- G06T7/74—Determining position or orientation of objects or cameras using feature-based methods involving reference images or patches
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Image Analysis (AREA)
Abstract
基于改进的下逼近非负矩阵对极几何估计方法,涉及计算机视觉技术。提供伴随着离群点和模型假设修剪技术的基于改进的下逼近非负矩阵对极几何估计方法。首先,使用误匹配修剪技术分析匹配对之间的关系剔除离群点(误匹配点)的影响。接着,使用模型假设修剪技术来选择有意义的模型假设。然后,引入空间约束项(空间上相邻的数据点更可能属于相同的模型假设)和稀疏约束项(稀疏非负元素更能体现数据点对模型的一致特征)到下逼近非负矩阵。最后,使用交替迭代法求解下逼近非负矩阵的u和v,从u中提取多结构模型。
Description
技术领域
本发明涉及计算机视觉技术,尤其是涉及一种基于改进的下逼近非负矩阵对极几何估计方法。
背景技术
随着科技的进步,计算机视觉已经融入到人类生活的方方面面。计算机视觉的研究任务是通过模拟人类的视觉功能,把自然场景信息高效地转换为数字信息。而对极几何估计是通过分析两幅图像之间的约束关系,减缓噪声和离群点对真实数据的污染,极大地提高了模型拟合的鲁棒性。模型拟合是计算机视觉的一项基础性研究任务,其在图像拼接、三维重建、运动分割、平面检测和消失点检测等方面都具有广泛的应用。
针对多结构模型拟合问题,一些基于分解的方法(NMF)(比如[1]Magri,Luca,andAndrea Fusiello,"Robust Multiple Model Fitting with Preference Analysis andLow-rank Approximation,"British Machine Vision Conference,Vol.20.No.1-20,2015;[2]Cai,Deng and He,Xiaofei and Wang,Xuanhui and Bao,Hujun and Han,Jiawei,"Locality preserving nonnegative matrix factorization,"InternationalJoint Conference on Artificial Intelligence,2009)通过执行矩阵分解技术,使用类似M估计样本与一致策略来处理多模型拟合问题,使用鲁棒主成分分析和对称非负矩阵分解实现多模型的提取。然而,NMF难以在合理的时间内找到全局最优值,它通常应用奇异值分解(SVD)来求解最优秩一的NMF。这种情况下,输入矩阵与其秩一近似之间的差异不一定是非负的,这可能导致SVD因素分解出现负值。一些基于双聚类的因素分析方法(比如[3]Denitto,Matteo and Magri,Luca and Farinelli,Alessandro and Fusiello,Andreaand Bicego,Manuele,"Multiple structure recovery via probabilisticbiclustering,"Joint IAPR International Workshops on Statistical Techniques inPattern Recognition and Structural and Syntactic Pattern Recognition,2016;[4]Tepper,Mariano,and Guillermo Sapiro,"Fast l1-nmf for multiple parametricmodel estimation,"Tehnical Report,2016)常用于分离出偏好矩阵中具有一致“行为”的数据。这些方法的基本思想是将数据按层分解,每个层次对应一个不同的双集类。FABIA方法(比如[4])则是在偏好矩阵上使用l1进行非负分解,然后采用偏好分析或一致性分析推导出特征向量,最后采用变分期望最大化算法来估计模型参数。虽然这种方法在性能上有较大的提升,但它们需要增加后处理步骤以获得该信息的双聚类因素。
下逼近非负矩阵(NMU)(比如[6]Gillis,Nicolas,andGlineur,"Usingunderapproximations for sparse nonnegative matrix factorization,"PatternRecognition,Vol.43.No.1676-1687,2010)提供了一种高效的算法来处理参数估计问题。为了在合理的时间内找到最优的非负秩一近似解,它采用了拉格朗日松弛法,并允许用户通过停止递归过程以找到NMU的近似解。最近,下逼近非负矩阵分解方法(比如[5]Tepper,Mariano,and Guillermo Sapiro,"Nonnegative matrix underapproximation forrobust multiple model fitting,"Computer Vision and Pattern Recognition,2017)被提出用于偏好矩阵分析以解决多结构模型拟合问题。这种方法基于一种交替方向乘法器(ADMM)来求解下逼近非负矩阵,再使用鲁棒统计方法提取模型实例。然而它所使用的多高斯采样方法非常耗时,尤其是在数据包含有大量离群点的情况下。
发明内容
本发明的目的在于提供伴随着离群点和模型假设修剪技术的基于改进的下逼近非负矩阵对极几何估计方法。
本发明包括以下步骤:
2)使用误匹配修剪技术移除离群点;
在步骤2)中,所述误匹配修剪技术的具体步骤可为:
2.2把一个假定为正确匹配的点[α,β,1]T加入到Ax和Ay(其中,(α,β,1)表示齐次坐标),得到:
2.4.使用最大重投影误差分离内点和离群点,如下所示:
其中,γt为一个阈值,θ为内点集,||·||F为F-范数操作。
3)使用模型假设修剪技术选择有意义的模型假设;
在步骤3)中,所述模型假设修剪技术的具体步骤可为:
那么,ξi的先验概率ρi可由下式计算:
然后,先验概率的熵E计算如下:
其中,log(·)表示对数函数。
其中,log(·)表示对数函数。
其中,exp(·)为指数函数,η通常设置为2.5;
5)使用改进的下逼近非负矩阵进行模型拟合;
在步骤5)中,所述改进的下逼近非负矩阵进行模型拟合方法的具体步骤可为:
5.1首先,下逼近非负矩阵(NMU)可以被定义为:
其中,Ψ为一个满足下逼近约束拉格朗日乘子法的矩阵;
5.2将空间约束项φ||Ωu||1和稀疏约束项ψ||u||1融入到下逼近非负矩阵中:
其中,uT(M-Ψ)v为最小二乘残差,φ和ψ分别为空间约束项和稀疏约束项的惩罚系数;Ω为数据点的邻居矩阵,这种融合方法能提取到更多局部信息和一致特征,提升模型拟合的鲁棒性。
6)应用交替迭代法求解u和v。u的每一行和列分别表示数据点和模型假设;
7)从u中提取模型实例,输出多结构模型。
本发明提出一种基于改进的下逼近非负矩阵对极几何估计方法。首先,该方法使用误匹配修剪技术分析匹配对之间的关系剔除离群点(误匹配点)的影响。接着,使用模型假设修剪技术来选择有意义的模型假设。然后,引入空间约束项(空间上相邻的数据点更可能属于相同的模型假设)和稀疏约束项(稀疏非负元素更能体现数据点对模型的一致特征)到下逼近非负矩阵。最后,使用交替迭代法求解下逼近非负矩阵的u和v,从u中提取多结构模型。
附图说明
图1为本发明实施例的流程示意图。
图2为基于改进后的下逼近非负矩阵对二进制数值图像进行重构的示意图。
图3为本发明在合成数据上的直线、圆和平面的拟合结果。在图3中,(a)为合成数据上的直线,(b)为合成数据上的圆,(c)为合成数据上的平面。
图4为本发明在AdelaideRMF数据集上的单应性估计结果。在图4中,(a)为Hartley,(b)为Sene,(c)Elderhalla。
图5为本发明在AdelaideRMF数据集上的基础矩阵分割结果。在图5中,(a)为Biscuitbook,(b)为Breadcube,(c)为Carchipscube。
图6为本发明在York Urban数据集上的消失点估计结果。在图6中,(a)为P1040822,(b)为P1020861,(c)为P1040812。
具体实施方式
下面结合附图和实施例对本发明的方法作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了实施方式和具体操作过程,但本发明的保护范围不限于下述的实施例。
参见图1,本发明实施例包括以下步骤:
S2.使用误匹配修剪技术移除离群点具体包括:
S2-2.把一个假定为正确匹配的点[α,β,1]T加入到Ax和Ay(其中,(α,β,1)表示齐次坐标),可得到:
S2-4.使用最大重投影误差分离内点和离群点,如下所示:
其中,γt为一个阈值,θ为内点集,||·||F为F-范数操作。
S3.使用模型假设修剪技术选择有意义的模型假设具体包括:
那么,ξi的先验概率ρi可由下式计算:
然后,先验概率的熵E计算如下:
其中,log(·)表示对数函数。
其中,log(·)表示对数函数。
其中,exp(·)为指数函数,η通常设置为2.5。
S5.使用改进的下逼近非负矩阵进行模型拟合方法具体包括:
S5-1.首先,下逼近非负矩阵(NMU)可以被定义为:
其中,Ψ为一个满足下逼近约束拉格朗日乘子法的矩阵。
S5-2.将空间约束项φ||Ωu||1和稀疏约束项ψ||u||1融入到下逼近非负矩阵中:
其中,uT(M-Ψ)v为最小二乘残差,φ和ψ分别为空间约束项和稀疏约束项的惩罚系数,Ω为数据点的邻居矩阵。这种融合方法能提取到更多局部信息和一致特征,提升模型拟合的鲁棒性。
基于改进后的下逼近非负矩阵对二进制数值图像进行重构的示意图参见图2。
S6.应用交替迭代法求解u和v。u的每一行和列分别表示数据点和模型假设。
S7.从u中提取模型实例,输出多结构模型。
本发明在合成数据集、AdelaideRMF数据集(包含19个单应性结构和19个基础矩阵结构)和York Urban数据集(包含102个室内和室外场景)进行多结构数据估计的部分结果如图3~6所示。图3为直线、圆形和平面拟合结果,图4为单应性估计结果,图5为基础矩阵分割结果,图6为消失点估计结果。
由实验结果所示,本发明所提出的对极几何估计方法获得了较准确的拟合结果。
Claims (3)
1.基于改进的下逼近非负矩阵对极几何估计方法,其特征在于包括以下步骤:
2)使用误匹配修剪技术移除离群点;
3)使用模型假设修剪技术选择有意义的模型假设;
其中,exp(·)为指数函数,η通常设置为2.5;
5)使用改进的下逼近非负矩阵进行模型拟合,具体步骤为:
5.1首先,下逼近非负矩阵被定义为:
满足u,v≥0,M-uvT≥0
满足Θ=uT(M-Ψ)v,
||u||2≤1,||v||2≤1,
其中,Ψ为一个满足下逼近约束拉格朗日乘子法的矩阵;
5.2将空间约束项φ||Ωu||1和稀疏约束项ψ||u||1融入到下逼近非负矩阵中:
其中,uT(M-Ψ)v为最小二乘残差,φ和ψ分别为空间约束项和稀疏约束项的惩罚系数;Ω为数据点的邻居矩阵,这种融合方法能提取到更多局部信息和一致特征,提升模型拟合的鲁棒性;
6)应用交替迭代法求解u和v,u的每一行和列分别表示数据点和模型假设;
7)从u中提取模型实例,输出多结构模型。
2.如权利要求1所述基于改进的下逼近非负矩阵对极几何估计方法,其特征在于在步骤2)中,所述误匹配修剪技术的具体步骤为:
2.2把一个假定为正确匹配的点[α,β,1]T加入到Ax和Ay,得到:
2.4.使用最大重投影误差分离内点和离群点,如下所示:
其中,γt为一个阈值,θ为内点集,||·||F为F-范数操作。
3.如权利要求1所述基于改进的下逼近非负矩阵对极几何估计方法,其特征在于在步骤3)中,所述模型假设修剪技术的具体步骤为:
那么,ξi的先验概率ρi由下式计算:
然后,先验概率的熵E计算如下:
其中,log(·)表示对数函数;
其中,log(·)表示对数函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910418346.6A CN110111391B (zh) | 2019-05-20 | 2019-05-20 | 基于改进的下逼近非负矩阵对极几何估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910418346.6A CN110111391B (zh) | 2019-05-20 | 2019-05-20 | 基于改进的下逼近非负矩阵对极几何估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110111391A CN110111391A (zh) | 2019-08-09 |
CN110111391B true CN110111391B (zh) | 2021-01-12 |
Family
ID=67491109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910418346.6A Active CN110111391B (zh) | 2019-05-20 | 2019-05-20 | 基于改进的下逼近非负矩阵对极几何估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110111391B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111767941B (zh) * | 2020-05-15 | 2022-11-18 | 上海大学 | 一种基于对称非负矩阵分解的改进谱聚类及并行化方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104050653A (zh) * | 2014-07-07 | 2014-09-17 | 西安电子科技大学 | 基于非负结构稀疏的高光谱图像超分辨率算法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150269314A1 (en) * | 2014-03-20 | 2015-09-24 | Rudjer Boskovic Institute | Method and apparatus for unsupervised segmentation of microscopic color image of unstained specimen and digital staining of segmented histological structures |
CN104392207B (zh) * | 2014-11-10 | 2017-07-14 | 天津大学 | 一种用于数字图像内容识别的特征编码方法 |
CN104408723B (zh) * | 2014-11-26 | 2017-04-12 | 浙江大学 | 一种基于非负矩阵临近的拉曼光谱图像解混方法 |
US9946767B2 (en) * | 2016-01-19 | 2018-04-17 | Conduent Business Services, Llc | Smoothed dynamic modeling of user traveling preferences in a public transportation system |
CN108875686A (zh) * | 2018-07-02 | 2018-11-23 | 扬州大学 | 一种基于非负稀疏的鉴别性低秩表示分类方法 |
-
2019
- 2019-05-20 CN CN201910418346.6A patent/CN110111391B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104050653A (zh) * | 2014-07-07 | 2014-09-17 | 西安电子科技大学 | 基于非负结构稀疏的高光谱图像超分辨率算法 |
Also Published As
Publication number | Publication date |
---|---|
CN110111391A (zh) | 2019-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | Text-guided graph neural networks for referring 3d instance segmentation | |
Chen et al. | Denoising hyperspectral image with non-iid noise structure | |
Xiong et al. | Fast algorithms to implement N-FINDR for hyperspectral endmember extraction | |
CN111738143B (zh) | 一种基于期望最大化的行人重识别方法 | |
CN111325165B (zh) | 考虑空间关系信息的城市遥感影像场景分类方法 | |
CN107977661B (zh) | 基于fcn与低秩稀疏分解的感兴趣区域检测方法 | |
CN108510013B (zh) | 基于低秩核心矩阵的改进稳健张量主成分分析的背景建模方法 | |
Galvanin et al. | Extraction of building roof contours from LiDAR data using a Markov-random-field-based approach | |
CN110838095B (zh) | 一种基于循环密集神经网络的单幅图像去雨方法和系统 | |
CN114187520B (zh) | 一种建筑物提取模型的构建及应用方法 | |
CN111046868A (zh) | 基于矩阵低秩稀疏分解的目标显著性检测方法 | |
Rasouli et al. | Visual saliency improves autonomous visual search | |
CN110111391B (zh) | 基于改进的下逼近非负矩阵对极几何估计方法 | |
Huynh et al. | Outlier detection in video sequences under affine projection | |
Ling et al. | Human object inpainting using manifold learning-based posture sequence estimation | |
Mittal et al. | Generalized projection based m-estimator: Theory and applications | |
CN112967210A (zh) | 一种基于全卷积孪生网络的无人机图像去噪方法 | |
CN114998647B (zh) | 基于注意力多实例学习的乳腺癌全尺寸病理图像分类方法 | |
CN111126123A (zh) | 一种基于压缩的增量核零空间变换行人再识别方法 | |
CN113763474B (zh) | 一种基于场景几何约束的室内单目深度估算方法 | |
CN115082337A (zh) | 一种基于双重全变分的高光谱混合噪声去除方法 | |
CN115063672A (zh) | 基于改进yolox模型的小目标检测方法 | |
Cai et al. | Building facade completion using semantic-synchronized GAN | |
Xu et al. | 6D-Diff: A Keypoint Diffusion Framework for 6D Object Pose Estimation | |
CN114495163A (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 |