CN104050675A - 基于三角形描述的特征点匹配方法 - Google Patents

基于三角形描述的特征点匹配方法 Download PDF

Info

Publication number
CN104050675A
CN104050675A CN201410303630.6A CN201410303630A CN104050675A CN 104050675 A CN104050675 A CN 104050675A CN 201410303630 A CN201410303630 A CN 201410303630A CN 104050675 A CN104050675 A CN 104050675A
Authority
CN
China
Prior art keywords
triangle
stack features
features point
supporting zone
remaining
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
CN201410303630.6A
Other languages
English (en)
Other versions
CN104050675B (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.)
Henan University of Technology
Original Assignee
Henan University of Technology
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 Henan University of Technology filed Critical Henan University of Technology
Priority to CN201410303630.6A priority Critical patent/CN104050675B/zh
Publication of CN104050675A publication Critical patent/CN104050675A/zh
Application granted granted Critical
Publication of CN104050675B publication Critical patent/CN104050675B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)

Abstract

本发明涉及一种基于三角形描述的特征点匹配方法，包括步骤：采集图像并输入计算机，利用已有特征点检测算子提取图像特征点，将图像中的特征点组合成三角形，确定三角形的支撑区域，对支撑区域进行子区域划分并计算描述子，计算三角形之间的相似性距离，建立三角形匹配，由三角形匹配确立特征点匹配。相比于已有基于单个特征点构造描述子的方法，本发明提供的方法构造的匹配描述子能够克服尺度误差及方向误差且在纹理简单或重复时稳定性更好，因而比已有方法具有更优的匹配性能。

Description

基于三角形描述的特征点匹配方法
技术领域
本发明涉及计算机视觉中的图像特征自动匹配领域，特别是数字图像中特征点的自动匹配方法。
背景技术
特征匹配技术在图像检索、物体识别、视频跟踪以及增强现实等领域有重要应用，目前已有图像特征点匹配方法主要有以下两类：
第一类是基于局部纹理信息构造匹配描述子的方法，主要代表有SIFT、SURF、GLOH等基于梯度特征匹配描述子（详见文献K. Mikolajczyk and C. Schmid, A performance evaluation of local descriptors. IEEE Transaction on Pattern Analysis and Machine Intelligence, 2005, 27(10): 1615–1630）与基于局部灰度序特征的描述方法（Fan B, Wu F and Hu Z. Rotationally invariant descriptors using intensity order pooling. IEEE Transaction on Pattern Analysis and Machine Intelligence, 2012, 34(10): 2031-2045），该类方法主要利用特征局部纹理信息进行特征匹配，首先为特征确定支撑区域，然后将支撑区域划分为一系列子区域，最后将每个子区域内包含的图像纹理信息表示为向量并组成匹配描述子。在图像形变存在时，该类方法构造的匹配描述子分辨力下降较大，获得的匹配数量较少。
第二类是基于几何一致性优化的匹配方法（如文献J.H. Lee and C.H. Won. Topology Preserving Relaxation Labeling for Non-rigid Point Matching. IEEE Transactions on Pattern Analysis and Machine Intelligence. 2011, 33(2): 427–432.），这类方法利用特征点之间的几何约束整体构造匹配目标优化函数，从而将匹配问题转化为一个目标函数优化问题。该类方法利用的特征之间的几何约束信息比第一类方法在图像形变下更为稳定。但该类方法需要解决的几何优化问题运算一般非常复杂，且运算量随着特征点数量增加呈几何级数增长，导致该类方法运算效率很低。
第一类方法由于简单且易于实现，在实际应用中使用较多。但该类法存在如下误差问题：（1）为实现尺度不变性而进行尺度选择导致存在尺度误差；（2）为实现旋转不变性而基于局部梯度直方图进行主方向估计产生方向误差；尺度误差与方向误差的存在影响了描述子的分辨能力；另外，在图像纹理比较简单或者重复纹理较多时，基于单个特征点描述的方法不太稳定而效果较差。
发明内容
本发明针对第一类特征点匹配方法存在的尺度误差与方向误差问题，提供一种不需要进行尺度选择与方向估计的具有更优性能的特征点匹配方法。为了实现本目的，本发明提供的基于三角组合描述的特征点匹配方法包括以下步骤：
步骤S1：从不同角度拍摄同一场景两幅不同图像并输入计算机；
步骤S2：利用已有特征点检测算子提取图像中特征点，如使用Harris算子；
步骤S3：利用Delaunay三角化技术将图像中的特征点组合成三角形；
步骤S4：通过计算外接圆确定每个三角形的支撑区域；
步骤S5：分别以圆心与三个顶点确定的矢量为主方向对支撑区域进行子区域划分并计算描述子；
步骤S51: 对支撑区域进行子区域划分；
步骤S52: 计算旋转不变特征向量；
步骤S53: 构造三角形描述子；
步骤S6：通过计算描述子之间的欧式距离计算三角形之间的相似性距离；
步骤S7：基于双向匹配准则进行三角形匹配；
步骤S8：由三角形匹配确立特征点匹配。
本发明提出的方法将单个特征点组合为三角形代替单个特征点作为基本匹配单位，相比于传统基于单个特征点构造描述子的方法，本发明提供的方法具有如下优势：（1）直接根据三角形的外接圆确定支撑区域可避免传统尺度选择环节引入的尺度误差；（2）直接利用圆心与顶点确定主方向可避免传统进行主方向估计引入的方向误差；（3）图像纹理比较简单或者重复较多时，基于单个特征点描述的方法不太稳定，本发明将单个特征点组合成三角形后更加稳定。
附图说明
图1所示为本发明基于三角形描述的特征点匹配方法流程图。
图2所示为确定三角形支撑区域示意图。
图3所示为计算三角形的3个描述子示意图。 
具体实施方式
如图1所示为本发明基于三角形描述的特征点匹配方法流程图，包括：采集图像并输入计算机，利用已有特征点检测算子提取图像特征点，将图像中的特征点组合成三角形，确定三角形的支撑区域，对支撑区域进行子区域划分并计算三角形的描述子，计算三角形之间的相似性距离，建立三角形匹配，由三角形匹配确立特征点匹配。各步骤的具体实施细节如下：
步骤S1：从不同角度拍摄同一场景两幅不同图像并输入计算机；
步骤S2：利用已有技术进行特征点检测，如使用Harris算子；
步骤S3：利用Delaunay三角化技术将图像中的特征点组合成三角形，具体方式为，记步骤S2获得的第1幅图像中特征点集合为                                                ，第2幅图像中的特征点集合为，其中为图像中特征点个数，利用Delaunay三角化技术将进行三角化获得个三角形，将进行三角化获得个三角形；
步骤S4：通过计算外接圆确定每个三角形的支撑区域，具体方式为，如图2所示，对于步骤S3获得的三角形，其中，利用如下公式计算其外接圆圆心和半径：
，
，
,,，
，
定义所述外接圆覆盖的图像区域为三角形的支撑区域，记为；
步骤S5：分别以圆心与三个顶点确定的矢量为主方向对支撑区域进行子区域划分并计算描述子，具体方式为，
步骤S51：对支撑区域进行子区域划分,具体方式为,如图3所示，对于步骤S4确定三角形的支撑区域，以圆心与顶点的连线向量为主方向将支撑区域在角度方向上等分为8个方向，然后再沿着径向方向将支撑区域等分为2部分，共计获得16个子区域；
步骤S52：计算旋转不变特征向量,具体方式为,记向量及其顺时针旋转90度方向对应的单位向量分别为与，记支撑区域内任一点的高斯梯度，将分别投影到、方向上获得点的旋转不变特征向量，其中、 ；
步骤S53：构造三角形描述子,具体方式为,对于三角形任一子区域，记内有个像素点及它们对应特征向量，通过计算特征向量的均值与标准差获得子区域的2维均值描述向量与2维标准差描述向量；将16个子区域对应的均值描述向量与标准差描述向量分别组成一个向量，可获得支撑区域的32维均值描述向量与32维标准差描述向量，分别进行归一化后获得区域的64维均值标准差描述向量，最后再进行归一化处理获得支撑区域的64维描述子；
记以向量代替重复执行步骤S51、S52、S53获得的描述子为；记以向量代替重复执行步骤S51、S52、S53获得的描述子为；最终每个三角形获得三个描述子；
步骤S6：通过计算描述子之间的欧式距离确定三角形之间的相似性距离，具体方式为，记第1幅图像中三角形经过步骤S5获得的3个描述子分别为,第2幅图像中三角形经过步骤S5获得的3个描述子分别为，则三角形与的相似性距离定义为；
步骤S7：基于双向匹配准则进行三角形匹配,具体方式为，对于第1幅图像中待匹配三角形和第2幅图像中的待匹配三角形，记中与相似性距离最小的三角形为,如果同时也是中与相似性距离最小的三角形，则 为一组三角形匹配；
步骤S8：由三角形匹配确立特征点匹配，具体方式为，对于步骤S7获得的描述子分别为与匹配三角形与，假设两个三角形的描述子之间的欧式距离最小,则可首先确立第1组特征点匹配，然后保持次序方式依次确立另两组特征点匹配，具体方式为，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,。
本发明提出的方法将单个特征点组合为三角形代替单个特征点作为基本匹配单位，相比于传统基于单个特征点构造描述子的方法，本发明提供的方法具有如下优势：（1）直接根据三角形的外接圆确定支撑区域可避免传统尺度选择环节引入的尺度误差；（2）直接利用圆心与顶点确定主方向可避免传统进行主方向估计引入的方向误差；（3）图像纹理比较简单或者重复较多时，基于单个特征点描述的方法不太稳定，本发明将单个特征点组合成三角形后更加稳定。

Claims (1)

1.一种数字图像中基于三角形描述的特征点匹配方法，其特征在于，包括步骤：
步骤S1：从不同角度拍摄同一场景两幅不同图像并输入计算机；
步骤S2：利用已有技术进行特征点检测，如使用Harris算子；
步骤S3：利用Delaunay三角化技术将图像中的特征点组合成三角形，具体方式为，记步骤S2获得的第1幅图像中特征点集合为                                                ，第2幅图像中的特征点集合为，其中为图像中特征点个数，利用Delaunay三角化技术将进行三角化获得个三角形，将进行三角化获得个三角形；
步骤S4：通过计算外接圆确定每个三角形的支撑区域，具体方式为，对于步骤S3获得的三角形，其中，利用如下公式计算其外接圆圆心和半径：
，
  ，
,,，
，
定义所述外接圆覆盖的图像区域为三角形的支撑区域，记为；
步骤S5：分别以圆心与三个顶点确定的矢量为主方向对支撑区域进行子区域划分并计算描述子，具体方式为，
步骤S51：对支撑区域进行子区域划分,具体方式为，对于步骤S4确定三角形的支撑区域，以圆心与顶点的连线向量为主方向将支撑区域在角度方向上等分为8个方向，然后再沿着径向方向将支撑区域等分为2部分，共计获得16个子区域；
步骤S52：计算旋转不变特征向量,具体方式为,记向量及其顺时针旋转90度方向对应的单位向量分别为与，记支撑区域内任一点的高斯梯度，将分别投影到、方向上获得点的旋转不变特征向量，其中、 ；
步骤S53：构造三角形描述子,具体方式为,对于三角形任一子区域，记内有个像素点及它们对应特征向量，通过计算特征向量的均值与标准差获得子区域的2维均值描述向量与2维标准差描述向量；将16个子区域对应的均值描述向量与标准差描述向量分别组成一个向量，可获得支撑区域的32维均值描述向量与32维标准差描述向量，分别进行归一化后获得区域的64维均值标准差描述向量，最后再进行归一化处理获得支撑区域的64维描述子；
记以向量代替重复执行步骤S51、S52、S53获得的描述子为；记以向量代替重复执行步骤S51、S52、S53获得的描述子为；最终每个三角形获得三个描述子；
步骤S6：通过计算描述子之间的欧式距离确定三角形之间的相似性距离，具体方式为，记第1幅图像中三角形经过步骤S5获得的3个描述子分别为,第2幅图像中三角形经过步骤S5获得的3个描述子分别为，则三角形与的相似性距离定义为；
步骤S7：基于双向匹配准则进行三角形匹配,具体方式为，对于第1幅图像中待匹配三角形和第2幅图像中的待匹配三角形，记中与相似性距离最小的三角形为,如果同时也是中与相似性距离最小的三角形，则 为一组三角形匹配；
步骤S8：由三角形匹配确立特征点匹配，具体方式为，对于步骤S7获得的描述子分别为与的匹配三角形与，假设两个三角形的描述子之间的欧式距离最小,则可首先确立第1组特征点匹配，然后保持次序方式依次确立另两组特征点匹配，具体方式为，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,，
如第1组特征点匹配为则剩下两组特征点匹配分别为,。
CN201410303630.6A 2014-06-30 2014-06-30 基于三角形描述的特征点匹配方法 Active CN104050675B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410303630.6A CN104050675B (zh) 2014-06-30 2014-06-30 基于三角形描述的特征点匹配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410303630.6A CN104050675B (zh) 2014-06-30 2014-06-30 基于三角形描述的特征点匹配方法

Publications (2)

Publication Number Publication Date
CN104050675A true CN104050675A (zh) 2014-09-17
CN104050675B CN104050675B (zh) 2017-02-08

Family

ID=51503458

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410303630.6A Active CN104050675B (zh) 2014-06-30 2014-06-30 基于三角形描述的特征点匹配方法

Country Status (1)

Country Link
CN (1) CN104050675B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105184786A (zh) * 2015-08-28 2015-12-23 大连理工大学 一种浮点型三角形特征描述方法
CN105678814A (zh) * 2016-01-05 2016-06-15 武汉大学 一种结合相位相关分析的建筑物立面影像重复纹理检测方法
CN106355577A (zh) * 2016-09-08 2017-01-25 武汉科技大学 基于特征状态与全局一致性的快速图像匹配方法及系统
CN106485264A (zh) * 2016-09-20 2017-03-08 河南理工大学 基于梯度序划分及映射策略的曲线描述及匹配方法
CN108537287A (zh) * 2018-04-18 2018-09-14 北京航空航天大学 基于图模型的图像闭环检测方法及装置
CN109263811A (zh) * 2018-09-29 2019-01-25 刘愉强 具有无磁低功耗的方向识别与定位功能的浮标及方法
CN109636790A (zh) * 2018-12-13 2019-04-16 北京理工大学 一种管路结构的识别方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103208003B (zh) * 2013-04-12 2015-03-11 大连理工大学 一种基于几何图形特征点形状描述子的方法
CN103473766B (zh) * 2013-09-03 2016-08-17 浙江工商大学 一种基于视觉感知的二维形状特征点匹配方法

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105184786A (zh) * 2015-08-28 2015-12-23 大连理工大学 一种浮点型三角形特征描述方法
CN105184786B (zh) * 2015-08-28 2017-10-17 大连理工大学 一种浮点型三角形特征描述方法
CN105678814A (zh) * 2016-01-05 2016-06-15 武汉大学 一种结合相位相关分析的建筑物立面影像重复纹理检测方法
CN105678814B (zh) * 2016-01-05 2018-03-13 武汉大学 一种结合相位相关分析的建筑物立面影像重复纹理检测方法
CN106355577A (zh) * 2016-09-08 2017-01-25 武汉科技大学 基于特征状态与全局一致性的快速图像匹配方法及系统
CN106355577B (zh) * 2016-09-08 2019-02-12 武汉科技大学 基于特征状态与全局一致性的快速图像匹配方法及系统
CN106485264A (zh) * 2016-09-20 2017-03-08 河南理工大学 基于梯度序划分及映射策略的曲线描述及匹配方法
CN108537287A (zh) * 2018-04-18 2018-09-14 北京航空航天大学 基于图模型的图像闭环检测方法及装置
CN109263811A (zh) * 2018-09-29 2019-01-25 刘愉强 具有无磁低功耗的方向识别与定位功能的浮标及方法
CN109263811B (zh) * 2018-09-29 2024-03-08 国家海洋局南海调查技术中心 具有无磁低功耗的方向识别与定位功能的浮标及方法
CN109636790A (zh) * 2018-12-13 2019-04-16 北京理工大学 一种管路结构的识别方法及装置
CN109636790B (zh) * 2018-12-13 2021-07-27 北京理工大学 一种管路结构的识别方法及装置

Also Published As

Publication number Publication date
CN104050675B (zh) 2017-02-08

Similar Documents

Publication Publication Date Title
CN104050675A (zh) 基于三角形描述的特征点匹配方法
CN111028277B (zh) 基于伪孪生卷积神经网络的sar和光学遥感图像配准方法
Li et al. Robust feature matching via support-line voting and affine-invariant ratios
Zhao et al. SPHORB: A fast and robust binary feature on the sphere
Buch et al. Pose estimation using local structure-specific shape and appearance context
CN102968777B (zh) 一种基于重叠区域sift特征点的图像拼接方法
Fan et al. Line matching leveraged by point correspondences
CN101833765B (zh) 基于双向匹配和三角形约束的特征匹配方法
CN102592281B (zh) 图像匹配方法
CN102663733B (zh) 基于特征组对的特征点匹配方法
CN103679702A (zh) 一种基于图像边缘矢量的匹配方法
Zhou et al. Progressive large scale-invariant image matching in scale space
Chen et al. Robust affine-invariant line matching for high resolution remote sensing images
CN104616297A (zh) 一种用于图像篡改取证的改进型sift算法
Li et al. 4FP-structure: A robust local region feature descriptor
Cordes et al. High-resolution feature evaluation benchmark
CN102122359A (zh) 一种图像配准方法及装置
Varytimidis et al. W α SH: weighted α-shapes for local feature detection
Cao et al. Fast and robust local feature extraction for 3D reconstruction
Zhang et al. Incorporating depth information into few-shot semantic segmentation
CN105631860A (zh) 基于局部排序方向直方图描述子的图像同名点提取方法
CN103336964A (zh) 一种基于模值差镜像不变性的sift图像匹配方法
CN105184327B (zh) 一种基于顶点三分策略的遥感图像特征点匹配方法
CN107229935B (zh) 一种三角形特征的二进制描述方法
Noury et al. How to overcome perceptual aliasing in ASIFT?

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