CN102592292A - 基于内积运算的对称图形中心定位方法 - Google Patents

基于内积运算的对称图形中心定位方法 Download PDF

Info

Publication number
CN102592292A
CN102592292A CN2011104318950A CN201110431895A CN102592292A CN 102592292 A CN102592292 A CN 102592292A CN 2011104318950 A CN2011104318950 A CN 2011104318950A CN 201110431895 A CN201110431895 A CN 201110431895A CN 102592292 A CN102592292 A CN 102592292A
Authority
CN
China
Prior art keywords
point
image
local maximum
symmetric
symmetry
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
CN2011104318950A
Other languages
English (en)
Other versions
CN102592292B (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 CN201110431895.0A priority Critical patent/CN102592292B/zh
Publication of CN102592292A publication Critical patent/CN102592292A/zh
Application granted granted Critical
Publication of CN102592292B publication Critical patent/CN102592292B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

本发明涉及数字图像中一种基于内积运算的对称图形中心定位方法,包括:采集图像并输入计算机、计算各像素点的梯度、计算图像的边缘图、基于内积运算计算各像素点的对称性能量获得图像的对称性能量分布图、在阈值约束下在对称性能量分布图上检测局部极大值点,最后,剔除非中心点的局部极大值点,输出对称图形的中心点。本发明提供的方法能够简单、高效的对图像中对称图形的中心进行准确定位,较已有方法简单,具有广泛的适用性。

Description

基于内积运算的对称图形中心定位方法
技术领域
本发明涉及计算机视觉中的图像特征自动检测领域,特别是数字图像中对称图形中心定位方法。
背景技术
图形中心定位是图像处理领域的一个重要问题,被广泛应用于车轮、X射线衍射图、同心圆、星云等图像中心的检测,为驾驶员、天文学家等提供重要的结构信息。目前,用于图形中心定位的方法有多种。Dammer等[1]利用Hough变换实现了X射线衍射图像中心点的确定。Wong等[2]提出一种基于运动信息的中心定位算法,从视频序列中获取出圆形和螺旋形物体的中心。而关于人眼中心的定位,多种方法被提出[3-4]。上述方法主要面向特定的研究对象,Cauchie等[5]利用自适应Hough变换,采用逐渐精确的方式提出了通用的静态图片中心定位算法;Chung[6]等在该算法的基础上,利用随机化方法,首先利用任意三个不同梯度方向的边缘点组合确定可能的中心点,然后再剔除错误的中心,实现图像中心点的准确测定。该算法的检测准确性较好,但由于采用边缘点的随机组合方法,运算较为复杂,计算量大。针对现有检测方法中存在的问题,本发明提出一种基于内积运算的对称图形中心定位方法,该方法引入数学中的内积运算,构造对称性描述子,获取图像的对称性能量分布图,进而实现图形中心的准确定位。该方法不需要已知图像的相关信息,运算简单,计算量小,易于实现。
参考文献:
[1]C.Dammer,P.Leleux,D.Villers,M.Dosiere,Use of the Hough transform to determine the center of digitizedX-ray diffraction patterns,Nuclear Instruments and Methods in Physics Research Section B:BeamInteractions with Materials and Atoms 132(1)(1997)214-220.
[2]K.Y.Wong,C.L.Yip,Identifying centers of circulating and spiraling vector field patterns and its applications,Pattern Recognition 42(7)(2009)1371-1387.
[3]N.T.H.B.N.Huan,H.Kim,Location of iris based on circular and linear filters,in:International Conference onControl,Automation,Robotics and Vision,Hanoi,2008,pp.412-416.
[4]T.R.B.S.E.-R.V.Kith,M.EI-Sharkawy,S.Elnoubi,A feature and appearance based method for eye detectionon gray intensity face images,in:International Conference on Computer Engineering and Systems,Hahaii,2008,pp.41-47.
[5]J.Cauchie,V.Fioletand,D.Villers,Optimization of an Hough transform algorithm for the search of a center,Pattern Recognition 41(2)(2008)567-574.
[6]Kuo-Liang Chung,Yong-HuaiHuang,Jyun-PinWang,Ting-Chin Chang,Hong-Yuan Mark Liao,Fastrandomized algorithm for center-detection.Pattern Recognition 43(2010)2659-2665.
发明内容
本发明针对数字图像中的中心定位问题,提供一种基于内积运算的对称图形中心定位方法。为了实现本目的,本发明基于内积运算的对称图形中心定位方法,包括以下步骤:
步骤S1:采集图像并输入计算机;
步骤S2:利用高斯梯度模板计算各像素点的梯度;
步骤S3:利用Canny边缘检测算子计算图像的边缘图;
步骤S4:基于内积运算计算图像中各像素点的对称性能量,获得图像的对称性能量分布图;
步骤S5:在阈值约束下,在对称性能量分布图上检测局部极大值点;
步骤S6:验证步骤S5获得的局部极大值点,剔除不合理的局部极大值点,保留的局部极大值点即为对称图形的中心点。
本发明提供的基于内积运算的对称图形中心定位方法,主要利用图形轮廓点关于中心点对称的特性(即图形轮廓上某一点的梯度与该点关于中心点的对称点的梯度反向平行),引入数学中的内积运算,构造对称性描述子计算各像素点的对称性能量,从而获得图像的对称性能量分布图;然后在阈值约束下,在对称性能量分布图上检测局部极大值点,并利用约束条件剔除非中心点的局部极大值点,获得图形中心的准确定位。本发明提供的方法运算简单、易于实现,且计算效率较高。
附图说明
图1为本发明基于内积运算的对称图形中心定位方法流程图。
图2a为实施例中使用的原始图像。图2b为利用Canny算子在图2a上获得的边缘图;图2c为本发明提供的方法在图2a上获得的对称性能量分布图;图2d为在图2c上进行局部极大值检测获取的局部极大值点;图2e为利用本方法获得的对称图形中心定位结果。
具体实施方式
如图1所示为本发明基于内积运算的对称图形中心定位方法流程图,包括:采集图像并输入计算机、计算各像素点的梯度、计算图像的边缘图、计算各像素点的对称性能量获得图像的对称性能量分布图、在阈值约束下在对称性能量分布图上检测局部极大值点,最后,剔除非中心点的局部极大值点,输出对称图形的中心点。各步骤的具体实施细节如下:
步骤S1:采集图像并输入计算机。
步骤S2:利用高斯梯度模板计算各像素点的梯度。
步骤S3:利用Canny边缘检测算子计算图像的边缘图。
步骤S4:基于内积运算计算图像中各像素点的对称性能量,获得图像的对称性能量分布图。给定一个检测半径r,对于图像中任意一点X,利用公式 E ( X ) = &Sigma; P i , P i &prime; &Element; G ( X , r ) e - d i &CenterDot; abs ( < P i , P i &prime; > mag ( P i ) &CenterDot; mag ( P i &prime; ) ) 计算点X处的对称性能量,从而获得图像的对称性能量分布图;其中G(X,r)表示边缘图中以点X为圆心,r为半径的一个圆形区域;Pi和P′i为G(X,r)内的一个对称点对,记点Pi关于点X的对称点为点Pis,则点P′i为以点Pis为中心的5×5邻域内距离点Pis最近的边缘点;i=1,2,...,NoX,其中NoX为G(X,r)内关于点X对称的边缘点对的个数;<Pi,P′i>表示点Pi与P′i的内积运算;为权重系数,di为点P′i和点Pis之间的距离;mag(·)表示梯度幅值;abs(·)表示取绝对值。
步骤S5:在阈值约束下,在对称性能量分布图上检测局部极大值点。计算阈值T=k·Mean(E),其中Mean(E)表示步骤S4获得的对称性能量分布图E的均值,比例系数k的取值范围为2~5,在对称性能量分布图上检测对称性能量大于T且在5×5邻域内为局部极大值的点。
步骤S6:验证步骤S5获得的局部极大值点,剔除不合理的局部极大值点,保留的局部极大值点即为对称图形的中心点。对于步骤S5获得的任一局部极大值点P,记圆形区域G(P,r)内关于点P对称的对称点对的个数为NoP,将不满足条件NoP>Ts·2πr的局部极大值点剔除,其中Ts为比例系数,一般取0.4~0.8,则被保留的局部极大值点即为图像中对称图形的中心点。
实施例
如图2a所示一幅原始灰度图像,图像大小为61×120,图像中包含2个大小相同的圆形硬币。下面介绍使用本发明提供的方法定位对称图形中心的具体实施步骤:
步骤S1:采集图像并输入计算机。
步骤S2:利用高斯梯度模板计算各像素点的梯度。
步骤S3:利用Canny边缘检测算子计算图像的边缘图。利用Canny算子进行边缘检测时的高斯尺度设置为1.0,进行连接时的高低阈值参数分别设置为0.2与0.3,如图2b所示为使用Canny算子对图2a进行边缘检测后获得的边缘图。
步骤S4:基于内积运算计算图像中各像素点的对称性能量,获得图像的对称性能量分布图。给定一个检测半径r=35,对于图像中任意一点X,利用公式 E ( X ) = &Sigma; P i , P i &prime; &Element; G ( X , r ) e - d i &CenterDot; abs ( < P i , P i &prime; > mag ( P i ) &CenterDot; mag ( P i &prime; ) ) 计算点X处的对称性能量,从而获得图像的对称性能量分布图;其中G(X,r)表示边缘图中以点X为圆心,r为半径的一个圆形区域;Pi和P′i为G(X,r)内的一个对称点对,记点Pi关于点X的对称点为点Pis,则点P′i为以点Pis为中心的5×5邻域内距离点Pis最近的边缘点;i=1,2,...,NoX,其中NoX为G(X,r)内关于点X对称的边缘点对的个数;<Pi,P′i>表示点Pi与P′i的内积运算;
Figure BSA00000641871900052
为权重系数,di为点P′i和点Pis之间的距离;mag(·)表示梯度幅值;abs(·)表示取绝对值。如图2c所示为图2a的对称性能量分布图。
步骤S5:在阈值约束下,在对称性能量分布图上检测局部极大值点。设置k=2,Mean(E)=2.4086,计算阈值T=4.8173,在对称性能量分布图上检测对称性能量大于T且在5×5邻域内为局部极大值的点,共计获得6个局部极大值点,如图2d为在图2c上进行局部极大值检测获取的局部极大值点。
步骤S6:验证步骤S5获得的局部极大值点,剔除不合理的局部极大值点,保留的局部极大值点即为对称图形的中心点。设置比例系数Ts=0.4,对于获得的6个局部极大值点Pi(i=1,2,...,6),将不满足条件
Figure BSA00000641871900053
的局部极大值点剔除,则被保留的3个局部极大值点即为图像中对称图形的中心点,如图2e所示。
可以看出,本发明提供的方法利用内积这一数学运算,结合对称图形的特征,构造对称性能量描述子,获取图像的对称性能量分布图,实现图像对称性结构的表示,其中,对称性能量较大的点对应于图形的中心点。
该方法不需要已知图像的相关信息,构造描述子后,不需要其他复杂的计算,运算简单且易于实现。在计算时间方面,均采用C++语言进行编程实现,利用本发明提供的方法在实施例图像上进行检测耗时0.5s,而利用随机化方法对图2a进行检测分别耗时0.8s,本发明提供的方法计算效率更高。

Claims (1)

1.一种基于内积运算的对称图形中心定位方法,其特征在于,包括步骤:
步骤S1:采集图像并输入计算机;
步骤S2:利用高斯梯度模板计算各像素点的梯度;
步骤S3:利用Canny边缘检测算子计算图像的边缘图;
步骤S4:基于内积运算计算图像中各像素点的对称性能量,获得图像的对称性能量分布图;具体方法为:给定一个检测半径r,对于图像中任意一点X,利用公式 E ( X ) = &Sigma; P i , P i &prime; &Element; G ( X , r ) e - d i &CenterDot; abs ( < P i , P i &prime; > mag ( P i ) &CenterDot; mag ( P i &prime; ) ) 计算点X处的对称性能量,从而获得图像的对称性能量分布图;其中G(X,r)表示边缘图中以点X为圆心,r为半径的一个圆形区域;Pi和P′i为G(X,r)内的一个对称点对,记点Pi关于点X的对称点为点Pis,则点P′i为以点Pis为中心的5×5邻域内距离点Pis最近的边缘点;i=1,2,...,NoX,其中NoX为G(X,r)内关于点X对称的边缘点对的个数;<Pi,P′i>表示点Pi与P′i的内积运算;
Figure FSA00000641871800012
为权重系数,di为点P′i和点Pis之间的距离;mag(·)表示梯度幅值;abs(·)表示取绝对值;
步骤S5:在阈值约束下,在对称性能量分布图上检测局部极大值点;具体方法为:计算阈值T=k·Mean(E),其中Mean(E)表示步骤S4获得的对称性能量分布图E的均值,比例系数k的取值范围为2~5,在对称性能量分布图上检测对称性能量大于T且在5×5邻域内为局部极大值的点;
步骤S6:验证步骤S5获得的局部极大值点,剔除不合理的局部极大值点,保留的局部极大值点即为对称图形的中心点;具体方法为:对于步骤S5获得的任一局部极大值点P,记圆形区域G(P,r)内关于点P对称的对称点对的个数为NoP,将不满足条件NoP>Ts·2πr的局部极大值点剔除,其中Ts为比例系数,一般取0.4~0.8,则被保留的局部极大值点即为图像中对称图形的中心点。
CN201110431895.0A 2011-12-12 2011-12-12 基于内积运算的对称图形中心定位方法 Expired - Fee Related CN102592292B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110431895.0A CN102592292B (zh) 2011-12-12 2011-12-12 基于内积运算的对称图形中心定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110431895.0A CN102592292B (zh) 2011-12-12 2011-12-12 基于内积运算的对称图形中心定位方法

Publications (2)

Publication Number Publication Date
CN102592292A true CN102592292A (zh) 2012-07-18
CN102592292B CN102592292B (zh) 2014-07-16

Family

ID=46480882

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110431895.0A Expired - Fee Related CN102592292B (zh) 2011-12-12 2011-12-12 基于内积运算的对称图形中心定位方法

Country Status (1)

Country Link
CN (1) CN102592292B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102831610A (zh) * 2012-08-13 2012-12-19 中国科学院自动化研究所 基于内积与距离分布的椭圆快速检测方法
CN103034998A (zh) * 2012-12-04 2013-04-10 中国科学院自动化研究所 旋转对称图形的中心与旋转角检测方法及装置
CN108876845A (zh) * 2018-06-14 2018-11-23 歌尔股份有限公司 菲涅尔纹中心的确定方法和装置
CN110031014A (zh) * 2019-03-27 2019-07-19 浙江亚特电器有限公司 基于图案识别的视觉定位方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6019151A (en) * 1997-01-07 2000-02-01 Eastman Kodak Company Printing onto discs such as compact discs and the like
US6859555B1 (en) * 2000-09-19 2005-02-22 Siemens Corporate Research, Inc. Fast dominant circle detection through horizontal and vertical scanning
CN102096820A (zh) * 2011-01-12 2011-06-15 河南理工大学 数字图像中基于距离分布的正方形检测方法
CN102201118A (zh) * 2010-12-30 2011-09-28 天津大学 应用于轮胎x光图像胎冠部位带束层边界定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6019151A (en) * 1997-01-07 2000-02-01 Eastman Kodak Company Printing onto discs such as compact discs and the like
US6859555B1 (en) * 2000-09-19 2005-02-22 Siemens Corporate Research, Inc. Fast dominant circle detection through horizontal and vertical scanning
CN102201118A (zh) * 2010-12-30 2011-09-28 天津大学 应用于轮胎x光图像胎冠部位带束层边界定位方法
CN102096820A (zh) * 2011-01-12 2011-06-15 河南理工大学 数字图像中基于距离分布的正方形检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LI-QIN JIA 等: "A Fast Randomized Circle Detection Algorithm", 《2011 4TH INTERNATIONAL CONGRESS ON IMAGE AND SIGNAL PROCESSING (CISP)》 *
刘红敏 等: "基于基元表示的多边形检测方法", 《自动化学报》 *
王志衡 等: "内积能量与边缘检测", 《计算机学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102831610A (zh) * 2012-08-13 2012-12-19 中国科学院自动化研究所 基于内积与距离分布的椭圆快速检测方法
CN103034998A (zh) * 2012-12-04 2013-04-10 中国科学院自动化研究所 旋转对称图形的中心与旋转角检测方法及装置
CN108876845A (zh) * 2018-06-14 2018-11-23 歌尔股份有限公司 菲涅尔纹中心的确定方法和装置
CN108876845B (zh) * 2018-06-14 2022-02-25 歌尔股份有限公司 菲涅尔纹中心的确定方法和装置
CN110031014A (zh) * 2019-03-27 2019-07-19 浙江亚特电器有限公司 基于图案识别的视觉定位方法
CN110031014B (zh) * 2019-03-27 2024-01-26 浙江亚特电器股份有限公司 基于图案识别的视觉定位方法

Also Published As

Publication number Publication date
CN102592292B (zh) 2014-07-16

Similar Documents

Publication Publication Date Title
CN104063702B (zh) 一种基于遮挡修复和局部相似性匹配的三维步态识别方法
CN106251355B (zh) 一种融合可见光图像与对应夜视红外图像的检测方法
CN102592292B (zh) 基于内积运算的对称图形中心定位方法
CN110443879B (zh) 一种基于神经网络的透视误差补偿方法
CN103903246A (zh) 物体检测方法和装置
CN107742289A (zh) 一种基于机器视觉回转体工件检测方法
CN107292923B (zh) 基于深度图挖掘的后向传播图像视觉显著性检测方法
CN104050666A (zh) 基于分割的脑部mr图像配准方法
CN105096300A (zh) 对象检测方法和设备
CN113093216A (zh) 一种激光雷达和相机融合的不规则物体测量方法
CN106446833B (zh) 一种用于复杂场景图像识别的多通道仿生视觉方法
CN110070503A (zh) 基于卷积神经网络的尺度校准方法、系统及介质
CN104899852A (zh) 一种基于遗传算法及支持向量机的微小颗粒检测系统标定方法
CN103514599B (zh) 一种基于邻域总变分的影像最佳分割尺度选择方法
CN102074011B (zh) 数字图像中任意三角形的检测方法
CN113902779B (zh) 一种基于张量投票方法的点云配准方法
CN116883945B (zh) 一种融合目标边缘检测和尺度不变特征变换的人员识别定位方法
CN101393607B (zh) 利用快速边界跟踪实现脑血管图像识别的方法
Wijenayake et al. Automatic detection and decoding of photogrammetric coded targets
CN102096820B (zh) 数字图像中基于距离分布的正方形检测方法
CN102072707B (zh) 数字图像中圆的中心与半径快速检测方法
Cai et al. Deep representation and stereo vision based vehicle detection
CN115935767B (zh) 基于统计力学的三维船舶交通复杂度识别感知方法及装置
CN103955929B (zh) 图像局部边缘模式与非边缘模式判断方法及判断装置
CN105740869A (zh) 一种基于多尺度多分辨率的方形算子边缘提取方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Liu Hongmin

Inventor after: Wang Zhiheng

Inventor after: Jiang Guoquan

Inventor after: Xu Huilin

Inventor after: Wu Jun

Inventor before: Liu Hongmin

Inventor before: Wang Zhiheng

Inventor before: Jiang Guoquan

Inventor before: Li Jun

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: LIU HONGMIN WANG ZHIHENG JIANG GUOQUAN LI JUN TO: LIU HONGMIN WANG ZHIHENG JIANG GUOQUAN XU HUILIN WU JUN

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140716

Termination date: 20161212