CN102385750A - 基于几何关系的直线匹配方法及*** - Google Patents

基于几何关系的直线匹配方法及*** Download PDF

Info

Publication number
CN102385750A
CN102385750A CN2011101694533A CN201110169453A CN102385750A CN 102385750 A CN102385750 A CN 102385750A CN 2011101694533 A CN2011101694533 A CN 2011101694533A CN 201110169453 A CN201110169453 A CN 201110169453A CN 102385750 A CN102385750 A CN 102385750A
Authority
CN
China
Prior art keywords
straight line
image
coordinate
point
same image
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.)
Granted
Application number
CN2011101694533A
Other languages
English (en)
Other versions
CN102385750B (zh
Inventor
陈茂银
周东华
侯强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN 201110169453 priority Critical patent/CN102385750B/zh
Publication of CN102385750A publication Critical patent/CN102385750A/zh
Application granted granted Critical
Publication of CN102385750B publication Critical patent/CN102385750B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

本发明公开了基于几何关系的直线匹配方法及***。该直线匹配方法包括:步骤1,对若干个图像进行直线检测;步骤2,对从同一个图像检测到的直线进行投影变换;步骤3,求取同一个图像检测到的直线在投影面上形成的圆弧的交点;步骤4,生成同一个图像检测到的直线的匹配特征;步骤5,利用从各图像检测到的直线对应的匹配特征进行直线匹配。本发明无需预先求取基本矩阵,避免了点匹配精度对直线匹配的精度的影响,提高了匹配算法的效率。

Description

基于几何关系的直线匹配方法及***
技术领域
本发明涉及图像处理领域,尤其涉及基于几何关系的直线匹配方法及***。
背景技术
直线是三维世界中构成物体中重要几何元素。从三维场景中获得的二维图像中物体的直线特征仍然保留。不同视角下对同一场景拍摄的图像中一部分直线保持不变或直线间具有一定的相关性,找到不同图像中的对应的直线,利用直线之间的对应关系确定两幅图像间的匹配及变换关系有重要的意义。三维场景的恢复主要利用点及线将三维的物体目标等再现出来,匹配不同图像中的直线是三维重建的重要基础。
直线匹配的方法主要采用极几何约束进行初匹配,再利用相关法进行精匹配。极几何约束方法是利用两视角下获得的图像的几何关系,建立了图像点间的对应几何变换关系,在一幅图像中的点在另一视角的图像中的几何约束可通过基本矩阵表示,获得图像间的基本矩阵利用直线的两个端点可得到另一视角图像中满足几何约束的区域,属于几何约束的区域直线作为候选直线,完成粗匹配。精匹配采用相关的方法,选取直线邻域内像素的灰度变化作为相似度实现匹配。
采用极几何约束的方法需要预先求取基本矩阵,基本矩阵的求取是采用匹配点对求线性方程的方法实现,在直线匹配前需要进行点匹配,直线匹配的精度受点匹配精度影响,匹配算法效率不高。
发明内容
针对现有技术中存在的上述问题,本发明提供了基于几何关系的直线匹配方法及***。
本发明提供了基于几何关系的直线匹配方法,包括:
步骤1,对若干个图像进行直线检测;
步骤2,对从同一个图像检测到的直线进行投影变换;
步骤3,求取同一个图像检测到的直线在投影面上形成的圆弧的交点;
步骤4,生成同一个图像检测到的直线的匹配特征;
步骤5,利用从各图像检测到的直线对应的匹配特征进行直线匹配。
在一个示例中,步骤1中,对图像的边缘进行检测从而检测出图像中的直线。
在一个示例中,步骤2中,以图像
Figure BDA0000070251860000021
的几何中心作为坐标变换的中心点,投影中心为单位球面的中心点,其中L1,L2,L3,L L,Ln为同一个图像中检测到的直线。
在一个示例中,步骤2中,设坐标变换中心点在图像中的坐标为(qx0 qy0),同一图像中检测到的直线上任一点的坐标为(qx qy),则投影变换关系为:
q x - q x 0 q y - q y 0 - 1 | | q x - q x 0 q y - q y 0 - 1 | | .
在一个示例中,步骤3中,设同一图像中检测到的两条直线对应的第一圆弧上和第二圆弧上的两点的坐标向量分别为
Figure BDA0000070251860000023
Figure BDA0000070251860000024
Figure BDA0000070251860000025
则圆弧的交点为
Figure BDA0000070251860000026
在一个示例中,步骤4包括:
步骤41,选取弧线的交点所对应的在直线上的两点之间距离最远的坐标点pmax1,pmax2,并以
Figure BDA0000070251860000027
的中点作为新的坐标变换中心点,并将同一图像中检测到的直线进行新的投影变换;
Figure BDA0000070251860000028
是坐标点pmax1,pmax2在X轴的中点,
Figure BDA0000070251860000031
是坐标点pmax1,pmax2在Y轴的中点;
步骤42,求取坐标点pmax1,pmax2在新的投影下所在圆弧的单位法向量,并将该单位法向量设置为参考向量;
步骤43,计算同一图像中检测到的直线在新投影下的坐标向量与参考向量的夹角;
步骤44,将夹角进行排序后作为直线的匹配特征。
在一个示例中,步骤5中,利用距离最近准则对直线的匹配特征进行匹配。
本发明提供了一种实现基于几何关系的直线匹配方法的***。
本发明无需预先求取基本矩阵,避免了点匹配精度对直线匹配的精度的影响,提高了匹配算法的效率。
附图说明
下面结合附图来对本发明作进一步详细说明,其中:
图1是本发明的基于几何关系的直线匹配方法流程图;
图2是投影变换模型图;
图3是圆弧交点求取原理图。
具体实施方式
本发明提出了一种新的直线匹配方法,利用图像中直线间的几何关系进行匹配判断。图像中的直线与直线之间有平行和相交关系,不同视角下获得的图像直线间的相交或平行关系在一定条件下是不变的,从一条直线出发确定出直线和其他直线的几何关系,利用几何关系的一直性判别是否为同一直线。将直线投影到球面上后,每两条直线在球面上存在两个交点,该两交点为球面大圆的两个极点,选取极点作为几何关系的表示量,将每一条直线同其他直线的几何关系采用向量表示,利用向量相关性作为匹配准则进行匹配。
本发明提供的方法如图1所示,其使用了2个图像(不限于2个图像),即图像1和图像2,但对每个图像的处理流程相同。本发明提供的方法包括:
步骤101、对图像1和图像2进行直线检测。进行直线检测时,通过对图像边缘检测,进而从图像边缘中检测出直线。原两幅图像分别记为:f1(x y)和f2(x y),检测出图像f1(x y)中的直线n条分别记为
Figure BDA0000070251860000041
图像f2(x y)中的直线m条分别记为
Figure BDA0000070251860000042
属于直线上
Figure BDA0000070251860000043
上的第j个点表示为
Figure BDA0000070251860000044
仅有直线点的图像记为∏1,∏2,表示为
Figure BDA0000070251860000045
Figure BDA0000070251860000046
步骤102、投影变换。经过投影变换图像上同一直线上的点在投影面上处于同一大圆弧上。不同圆弧间必然存在两个交点,且两交点的连线通过球心。图像中直线的几何及关系不同对应的交点的位置就不同,例如图像中两直线平行则球面上的圆弧交点靠近投影模型的赤道,两直线在图像中相交则球面上的圆弧交点偏离赤道。球面圆弧间的交点位置信息反应了直线间的几何关系。投影变换模型如图2所示,以图像f1(x y)中的直线
Figure BDA0000070251860000047
几何关系特征向量求取为例说明投影变换的过程。
Figure BDA0000070251860000048
为图像∏1中的直线,选取∏1的几何中心
Figure BDA0000070251860000049
为坐标变换的中心点,投影中心C选为球面的中心,选取单位球面,连接直线上点和投影中心C得投影线,投影线和球面的交点
Figure BDA00000702518600000411
Figure BDA00000702518600000412
在球面上的投影点,将∏1上的点采用投影变换关系投影到球面上。记:图像∏1
Figure BDA00000702518600000413
取值为1点为,
Figure BDA00000702518600000414
坐标变换中心点在图像中的坐标为
Figure BDA00000702518600000415
则投影点cp1坐标为:
q x 1 - q x 0 1 q y 1 - q y 0 1 - 1 | | q x 1 - q x 0 1 q y 1 - q y 0 1 - 1 | |
步骤103、求取几何交点。同一大圆弧及投影中心可确定唯一平面,两两相交圆弧的交点的连线是过圆弧及投影中心的两平面相交的交线。通过两相交平面的法向量可求出交线的方向向量。如前所示交线过球心,交点位于球面上,交点的坐标可通过交线的方向向量求出。以直线
Figure BDA0000070251860000051
和直线
Figure BDA0000070251860000052
为例,图像∏1中直线上两点在球面上的投影点
Figure BDA0000070251860000054
Figure BDA0000070251860000055
的坐标向量为
Figure BDA0000070251860000056
直线
Figure BDA0000070251860000057
上两点在球面上的投影点
Figure BDA0000070251860000058
Figure BDA0000070251860000059
的坐标向量
Figure BDA00000702518600000510
Figure BDA00000702518600000511
投影圆弧与
Figure BDA00000702518600000512
投影圆弧的交点为:
两交点关于球心对称,取向量第一分量不小于零的值记为求出其他圆弧同
Figure BDA00000702518600000515
投影圆弧的交点,在∏1上的对应点
Figure BDA00000702518600000517
同理依次求取各圆弧间的交点,即图像∏1中直线
Figure BDA00000702518600000518
Figure BDA00000702518600000519
投影圆弧的交点为
Figure BDA00000702518600000520
在原图像中的对应点为
Figure BDA00000702518600000521
步骤104、生成匹配特征。该步骤中,利用投影变换的逆变换求出圆弧交点对应的图像点,这些点即为图像中的直线交点。选取直线交点中距离最远两点的中点作为新的变换中心,将所有交点以新的变换中心重新投影到球面上,距离最远两点所在的圆弧的单位法向量作为参考向量。在新投影下交点的坐标向量和参考向量的夹角(或内积绝对值)作为交点的表示量,该表示量作为直线特征的一个元素。将属于同一直线上交点的表示量按从小到大排列,生成特征向量。
作为示例,选出直线交点中距离最远两点,
Figure BDA00000702518600000522
Figure BDA00000702518600000523
两点中点:
p mid 1 : p mid x 1 q mid y 1 = p max 1 x 1 + p max 2 x 1 2 p max 1 y 1 + p ma x 2 y 1 2 .
Figure BDA00000702518600000525
为变换中心将交点
Figure BDA00000702518600000526
投影到球面上求新的投影坐标。
直线交点
Figure BDA00000702518600000527
在新投影下投影坐标向量为:
Figure BDA00000702518600000528
Figure BDA00000702518600000529
在新投下的投影点的坐标向量:
Figure BDA00000702518600000530
过该两点圆弧的单位法向量为: C ‾ V max 1 = C ‾ p max 1 1 × C ‾ P max 1 1 | | C ‾ p max 1 1 × C ‾ p max 1 1 | | , 为参考向量
交点投影坐标向量
Figure BDA0000070251860000063
和参考向量
Figure BDA0000070251860000064
的内积:
Figure BDA0000070251860000065
直线
Figure BDA0000070251860000066
特征向量生成:
Figure BDA0000070251860000067
上各交点投影坐标向量
各交点投影坐标向量同参考向量的夹角: ar cos ( C ‾ V max 1 · [ C ‾ P 12 1 ; C ‾ P 13 1 ; ΛΛ ; C ‾ P 1 n - 1 1 ] T ) .
ar cos ( C ‾ V max 1 · [ C ‾ P 12 1 ; C ‾ P 13 1 ; ΛΛ ; C ‾ P 1 n - 1 1 ] T ) 按大小排序结果
sort ( ar cos ( C ‾ V max 1 · [ C ‾ P 12 1 ; C ‾ P 13 1 ; ΛΛ ; C ‾ P 1 n - 1 1 ] T ) )
则匹配特征表示为:
MG 1 = 1 sort ( ar cos ( C ‾ V max 1 · [ C ‾ P 12 1 ; C ‾ P 13 1 ; ΛΛ ; C ‾ P 1 n - 1 1 ] T ) ) .
生成所有直线的几何关系匹配特征,
Figure BDA00000702518600000613
采用同样步骤生成图像∏2的特征
Figure BDA00000702518600000614
步骤105、直线匹配。利用距离最近准则匹配直线的特征向量。以上所述仅为本发明的优选实施方式,但本发明保护范围并不局限于此。任何本领域的技术人员在本发明公开的技术范围内,均可对其进行适当的改变或变化,而这种改变或变化都应涵盖在本发明的保护范围之内。

Claims (8)

1.基于几何关系的直线匹配方法,其特征在于,包括:
步骤1,对若干个图像进行直线检测;
步骤2,对从同一个图像检测到的直线进行投影变换;
步骤3,求取同一个图像检测到的直线在投影面上形成的圆弧的交点;
步骤4,生成同一个图像检测到的直线的匹配特征;
步骤5,利用从各图像检测到的直线对应的匹配特征进行直线匹配。
2.如权利要求1所述的基于几何关系的直线匹配方法,其特征在于,步骤1中,对图像的边缘进行检测从而检测出图像中的直线。
3.如权利要求1或2所述的基于几何关系的直线匹配方法,其特征在于,步骤2中,以图像
Figure FDA0000070251850000011
的几何中心作为坐标变换的中心点,投影中心为单位球面的中心点,其中L1,L2,L3,L L,Ln为同一个图像中检测到的直线。
4.如权利要求3所述的基于几何关系的直线匹配方法,其特征在于,步骤2中,设坐标变换中心点在图像中的坐标为(qx0 qy0),同一图像中检测到的直线上任一点的坐标为(qx qy),则投影变换关系为:
q x - q x 0 q y - q y 0 - 1 | | q x - q x 0 q y - q y 0 - 1 | | .
5.如权利要求4所述的基于几何关系的直线匹配方法,其特征在于,步骤3中,设同一图像中检测到的两条直线对应的第一圆弧上和第二圆弧上的两点的坐标向量分别为
Figure FDA0000070251850000013
Figure FDA0000070251850000014
Figure FDA0000070251850000015
Figure FDA0000070251850000016
则圆弧的交点为
Figure FDA0000070251850000017
6.如权利要求5所述的基于几何关系的直线匹配方法,其特征在于,步骤4包括:
步骤41,选取弧线的交点所对应的在直线上的两点之间距离最远的坐标点pmax1,pmax2,并以
Figure FDA0000070251850000021
的中点作为新的坐标变换中心点,并将同一图像中检测到的直线进行新的投影变换;
Figure FDA0000070251850000022
是坐标点pmax1,pmax2在X轴的中点,
Figure FDA0000070251850000023
是坐标点pmax1,pmax2在Y轴的中点;
步骤42,求取坐标点pmax1,pmax2在新的投影下所在圆弧的单位法向量,并将该单位法向量设置为参考向量;
步骤43,计算同一图像中检测到的直线在新投影下的坐标向量与参考向量的夹角;
步骤44,将夹角进行排序后作为直线的匹配特征。
7.如权利要求6所述的基于几何关系的直线匹配方法,其特征在于,步骤5中,利用距离最近准则对直线的匹配特征进行匹配。
8.一种实现如权利要求1-7任意一项所述的基于几何关系的直线匹配方法的***。
CN 201110169453 2011-06-22 2011-06-22 基于几何关系的直线匹配方法及*** Active CN102385750B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110169453 CN102385750B (zh) 2011-06-22 2011-06-22 基于几何关系的直线匹配方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110169453 CN102385750B (zh) 2011-06-22 2011-06-22 基于几何关系的直线匹配方法及***

Publications (2)

Publication Number Publication Date
CN102385750A true CN102385750A (zh) 2012-03-21
CN102385750B CN102385750B (zh) 2013-07-10

Family

ID=45825139

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110169453 Active CN102385750B (zh) 2011-06-22 2011-06-22 基于几何关系的直线匹配方法及***

Country Status (1)

Country Link
CN (1) CN102385750B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103345642A (zh) * 2013-06-28 2013-10-09 华中科技大学 一种基于点线对偶的图像匹配方法
CN105160311A (zh) * 2015-08-26 2015-12-16 清华大学 直线段匹配方法及装置
CN105719309A (zh) * 2016-01-27 2016-06-29 大连理工大学 一种基于射影不变量的直线匹配方法
CN106296645A (zh) * 2015-06-25 2017-01-04 株式会社理光 图像处理方法以及图像处理装置
CN106934788A (zh) * 2015-12-30 2017-07-07 中国科学院沈阳自动化研究所 一种直线倾斜角的快速提取方法
CN109064440A (zh) * 2018-06-19 2018-12-21 广东工业大学 一种基于机器视觉的喇叭音圈焊线识别方法
WO2019024723A1 (zh) * 2017-08-01 2019-02-07 歌尔股份有限公司 特征点匹配结果处理方法和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1277079B1 (de) * 2000-04-25 2006-03-29 Rodenstock GmbH Verfahren zum berechnen eines progressiven brillenglases sowie verahren zur herstellung eines derartigen brillenglases
CN1941850A (zh) * 2005-09-29 2007-04-04 中国科学院自动化研究所 多摄像机下基于主轴匹配的行人跟踪方法
CN101561269A (zh) * 2009-05-26 2009-10-21 张征宇 近景摄影测量的特征线自动匹配方法
CN101621711A (zh) * 2009-07-23 2010-01-06 东南大学 采用两个相同圆进行摄像机标定的方法
CN101862177A (zh) * 2010-04-20 2010-10-20 中山大学中山眼科中心 一种视网膜裂孔三维定位的方法及其装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1277079B1 (de) * 2000-04-25 2006-03-29 Rodenstock GmbH Verfahren zum berechnen eines progressiven brillenglases sowie verahren zur herstellung eines derartigen brillenglases
CN1941850A (zh) * 2005-09-29 2007-04-04 中国科学院自动化研究所 多摄像机下基于主轴匹配的行人跟踪方法
CN101561269A (zh) * 2009-05-26 2009-10-21 张征宇 近景摄影测量的特征线自动匹配方法
CN101621711A (zh) * 2009-07-23 2010-01-06 东南大学 采用两个相同圆进行摄像机标定的方法
CN101862177A (zh) * 2010-04-20 2010-10-20 中山大学中山眼科中心 一种视网膜裂孔三维定位的方法及其装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
文贡坚: "一种基于特征编组的直线立体匹配全局算法", 《软件学报》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103345642A (zh) * 2013-06-28 2013-10-09 华中科技大学 一种基于点线对偶的图像匹配方法
CN103345642B (zh) * 2013-06-28 2016-05-25 华中科技大学 一种基于点线对偶的图像匹配方法
CN106296645A (zh) * 2015-06-25 2017-01-04 株式会社理光 图像处理方法以及图像处理装置
CN105160311A (zh) * 2015-08-26 2015-12-16 清华大学 直线段匹配方法及装置
CN105160311B (zh) * 2015-08-26 2018-08-28 清华大学 直线段匹配方法及装置
CN106934788A (zh) * 2015-12-30 2017-07-07 中国科学院沈阳自动化研究所 一种直线倾斜角的快速提取方法
CN106934788B (zh) * 2015-12-30 2020-11-24 中国科学院沈阳自动化研究所 一种直线倾斜角的快速提取方法
CN105719309A (zh) * 2016-01-27 2016-06-29 大连理工大学 一种基于射影不变量的直线匹配方法
CN105719309B (zh) * 2016-01-27 2018-08-14 大连理工大学 一种基于射影不变量的直线匹配方法
WO2019024723A1 (zh) * 2017-08-01 2019-02-07 歌尔股份有限公司 特征点匹配结果处理方法和装置
US11100361B1 (en) 2017-08-01 2021-08-24 Goertek Inc. Method and apparatus for processing feature point matching result
CN109064440A (zh) * 2018-06-19 2018-12-21 广东工业大学 一种基于机器视觉的喇叭音圈焊线识别方法

Also Published As

Publication number Publication date
CN102385750B (zh) 2013-07-10

Similar Documents

Publication Publication Date Title
CN102385750B (zh) 基于几何关系的直线匹配方法及***
Tang et al. 3D mapping and 6D pose computation for real time augmented reality on cylindrical objects
US20090153577A1 (en) Method and system for texturing of 3d model in 2d environment
CN102043966B (zh) 基于部分主分量分析和姿态估计联合的人脸识别方法
CN109671110B (zh) 一种局部几何结构约束的城区宽基线影像特征点匹配方法
CN107886546A (zh) 利用球像及公共自极三角形标定抛物折反射摄像机的方法
CN101236653A (zh) 图像校正方法及***
CN105678687A (zh) 基于图像内容的立体图像拼接方法
CN101159017A (zh) 基于局部图像纹理特征匹配的焊缝自动识别视觉方法
CN102509342A (zh) 一种共享增强现实场景中协同式虚实遮挡处理方法
CN113362457A (zh) 一种基于散斑结构光的立体视觉测量方法及***
CN106203429A (zh) 基于双目立体视觉复杂背景下的遮挡目标检测方法
Zhang et al. RI-Fusion: 3D object detection using enhanced point features with range-image fusion for autonomous driving
CN105809615A (zh) 点云数据影像化方法
CN105243137A (zh) 一种基于草图的三维模型检索视点选择方法
CN101383051B (zh) 一种基于图像重投影的视图合成方法
CN106204717A (zh) 一种立体图像快速三维重建方法与装置
CN103945206B (zh) 一种基于相似帧比较的立体图像合成***
CN1917658B (zh) 一种从单目图像序列生成立体图像序列的方法
CN103824076A (zh) 一种图像尺度不变特征的检测提取方法和***
CN116091987A (zh) 一种面向工业场景的多策略图像异常样本生成方法
CN104156973A (zh) 基于立体匹配的实时三维视频监控方法
CN113313751A (zh) 图像数据目标检测倾斜框自动标注方法及***
Li et al. 3D reconstruction based on stereovision and texture mapping
JP2010237941A (ja) マスク画像生成装置、3次元物体モデル情報生成装置及びプログラム

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant