CN107564006A - 一种利用Hough变换的圆形目标检测方法 - Google Patents
一种利用Hough变换的圆形目标检测方法 Download PDFInfo
- Publication number
- CN107564006A CN107564006A CN201711000106.1A CN201711000106A CN107564006A CN 107564006 A CN107564006 A CN 107564006A CN 201711000106 A CN201711000106 A CN 201711000106A CN 107564006 A CN107564006 A CN 107564006A
- Authority
- CN
- China
- Prior art keywords
- ballot
- circle
- circular target
- image
- hough transform
- 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
Links
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种利用Hough变换的圆形目标检测方法,包括如下步骤:步骤一,灰度转换与增强;步骤二,图像二值化;步骤三,Hough投票;步骤四,搜索峰值对应的投票距离ρm;步骤五,计算圆心坐标;步骤六,计算圆半径。本发明的圆形目标检测方法,通过在图像空间中进行粗采样,通过选用1D的累加数组,进一步降低了Hough变换的内存需求和时间消耗;通过求解方程组的解,得到圆形目标的圆心,提高圆心检测的精度。同时,直接在二值图像上进行圆形目标检测,不需要对图像进行边缘提取。
Description
技术领域
本发明涉及数字图像处理和模式识别领域,特别是在数字图像处理中利用Hough变换快速检测圆形目标的方法。
背景技术
圆形目标检测在计算机视觉和模式识别领域有着广泛的应用。例如,在工业生产线上有圆形的工件,道路上有圆形的交通信号,生物信息领域中的圆形虹膜等。在这些应用场景中,往往需要对圆形目标进行快速、准确地识别和定位,以便进一步分析处理。
常用的圆形目标检测的方法有:基于最大估计的拟合方法、基于Hough变换的投票方法。拟合方法受噪声影响大;Hough变换方法对图像噪声、目标不完整性的鲁棒性好、抗干扰能力强。
传统的利用Hough变换的圆检测算法将图像空间中一个边缘点变换到三维参数空间中的一个圆锥,所有边缘点在参数空间对应的圆锥簇相交于一峰值点。进而将圆检测问题转化为参数空间的峰值搜索问题。由于圆的投票方程为:(x-x0)2+(y-y0)2=r2,标准Hough变换定义了一个三维的参数空间,需要更多的时间消耗和大量的内存开销。
为此,出现了各种改进的Hough变换方法用于圆检测。
第一类是运用随机采样方法【Chung K L,Huang Y H,Shen S M,et al.Efficientsampling strategy and refinement strategy for randomized circle detection[J].Pattern Recognition,2012,45(1):252-263.】。随机Hough变换从图像空间中随机选择三个点,变换到参数空间中的一个点。虽然随机Hough变换降低了算法的时间开销,但是仍然采用的是3D的参数空间。当处理复杂图像时,由于无效采样而引入无效累积,导致算法性能大大降低。
第二类方法是运用低维参数空间的两步策略。首先利用圆边缘点的梯度方向信息【Jain N,Jain N.Coin Recognition Using Circular Hough Transform[J].International Journal of Electronics Communication&Computer Technology,2012,2(3):101-104】、或者圆弦的中垂线【Kim,HeungSoo,Kim,et al.A two-step circledetection algorithm from the intersecting chords[J].Pattern RecognitionLetters,2001,22(6):787-798.】,向一个2D累加器投票求得圆心坐标。然后运用1D的直方图投票估计圆的半径。这类算法采用2D+1D的参数空间,提高了Hough变换的执行效率,但是与传统Hough变换相比受噪声影响大,当出现同心圆时效果并不好。
相关的专利文献有:
【黄海清.利用改进霍夫变换进行类圆形检测的图像处理方法:CN,CN102568007A[P].2012.】、【CN201410594409.0一种基于Hough变换的圆检测方法】利用边缘点梯度方向信息降低了对内存的需求,提高了算法效率。而边缘梯度的提取对图像噪声非常敏感。
【张翼,孙永荣,刘晓俊,等.一种基于图像的圆形小目标精确检测方法:CN,CN103150730 A[P].2013.】、【CN201410598274.5一种圆形图形的检测方法】通过Hough变换和最小二乘拟合进行二次检测。而拟合算法对图像噪声敏感,并且受Hough检测结果的影响。
【Poyil B,Sripadarao A M.Hough transform for circles:,EP2738711[P].2014.】、【B·波伊勒,A·M·斯里帕达拉奥.Hough transform for circles:,CN103854279A[P].2014.】运用软、硬件结合的方案,提高Hough圆检测的效率。但是需要提供圆半径的范围。
【齐敏,辛红娟,吴志超,等.面向圆形标志符的两级优化霍夫变换快速检测定位方法:,CN104463920A[P].2015.】在平移Hough变换的基础上,通过两级Hough变换计算圆形标志符的圆心坐标。节约了投票计算的时间开销,但并没有降低空间需求。
【胡德志,孙碧亮,袁超飞.一种基于Hough变化的圆形识别方法及系统:,CN106780599 A[P].2017.】直接调用OpenCV函数库中的cvHoughCircles函数进行圆形检测。而该函数在实际应用中受限于参数选择等一系列困难。
另外,这些圆检测方法是基于边缘图像的,需要事先运用Sobel或Canny边缘检测技术计算机边缘点的梯度。如果梯度幅度大于设定的阈值,就认为是边缘点,同时记录该点的梯度方向。
发明内容
为了进一步降低Hough变换用于圆检测时的时间消耗和内存需求、提高算法效率,降低Hough变换圆检测的计算和空间复杂度,同时保持检测的精度,本发明提供一种利用Hough变换的圆形目标检测方法,忽略边缘检测,直接在二值图像上识别和定位圆形目标。
本发明的技术方案如下:
一种利用Hough变换的圆形目标检测方法,包括如下步骤:
步骤一,灰度转换与增强;
步骤二,图像二值化;
步骤三,Hough投票;
步骤四,搜索峰值对应的投票距离ρm;
步骤五,计算圆心坐标;
步骤六,计算圆半径。
进一步,所述步骤一,对输入的彩色图像进行变换,转换成灰度图像;并做图像平滑处理,降低图像噪声。
进一步,所述步骤二,运用otsu算法计算阈值,对经过步骤一处理后的灰度图像进行二值化处理,得到二值图像;如果像素点的灰度值大于该阈值,赋予该点灰度值255,否则,赋予该点灰度值0。
进一步,所述步骤三,从图像空间向参数空间进行投票时,在图像空间进行粗采样Δx=Δy>1,选用直线的投票方程:ρ=x·cosθ+y·sinθ;随机选择两个投票角度θ1和θ2;对图像中的每个黑像素点,运用直线投票方程,分别计算对应的投票距离ρ1和ρ2;分别向两个一维的累加数组A1(ρ)和A2(ρ)投票。
进一步,所述步骤四,分别在两个1D累加数组中,搜索最大投票值所对应的投票距离ρ;对于每个投票角θ,所有圆形目标都向累加数组投票后,当投票距离ρ经过圆心坐标时,投票值最大;设此时投票距离为ρm,并且最大投票值等于圆形目标的直径。
进一步,所述步骤五,圆心坐标(x0,y0)对应的投票距离为ρm;基于ρ1 m和ρ2 m构成了两个方程:ρ1 m=x0·cosθ1+y0·sinθ1和ρ2 m=x0·cosθ2+y0·sinθ2,解方程组得到圆心坐标(x0,y0)。
进一步,所述步骤六,在每个累加数组中,最大投票值为圆形目标的直径;当投票距离为ρm时,投票值A1(ρ1 m)和A2(ρ2 m)最大;根据在ρm处的最大投票值A(ρm)的一半,计算圆半径。
本发明的有益效果如下:
本发明提出了一种利用Hough变换的圆形目标检测方法,通过在图像空间中进行粗采样,通过选用1D的累加数组,进一步降低了Hough变换的内存需求和时间消耗;通过求解方程组的解,得到圆形目标的圆心,提高圆心检测的精度。同时,直接在二值图像上进行圆形目标检测,不需要对图像进行边缘提取。
附图说明
图1是本发明的检测圆形目标的流程图。
图2是从图像空间到参数空间的投票分析图。
图3是本发明实例虹膜图像中圆形目标的检测结果。
图4是本发明实例仿真图像中圆形目标的检测结果。
具体实施方式
下面结合附图对本发明的实施例作详细说明。
本发明的目的在于提供一种利用Hough变换的圆形目标检测方法,进一步降低Hough变换圆检测的计算和空间复杂度。其特征在于,按如下步骤实现:
如图1所示为本发明的流程图,算法以一幅数字图像为输入,输出圆形目标的圆心和半径。实施例包括如下步骤:
步骤一,图像的灰度转换与增强。首先对输入的原始彩色图像进行灰度转换,变换成灰度图像。然后,对灰度图像进图像平滑处理,降低图像噪声。
步骤二,图像二值化。运用otsu算法计算阈值,对经过步骤一处理后的灰度图像进行处理,得到二值图像。如果像素点的灰度值大于该阈值,就赋予该点灰度值为255,否则,赋予该点灰度值为0。此处,也可以根据实际应用,人工选取阈值。
步骤三,Hough投票。当从图像空间向参数空间进行投票时,(1)在图像空间进行粗采样Δx=Δy=2。此处,也可以选择更大的采样间隔Δx和Δy,这样图像中只有一小部分像素向参数空间进行投票。(2)选用直线的投票方程:ρ=x·cosθ+y·sinθ。选取大的采样间隔,这样图像中只有一小部分像素向参数空间进行投票。随机选择两个投票角度θ1和θ2。对图像中的每个黑像素点,运用直线投票方程,分别计算对应的投票距离ρ1和ρ2。定义两个1D的累加数组A1(ρ)和A2(ρ)。并分别向两个一维的累加数组A1(ρ)和A2(ρ)投票。
从图像空间向参数空间投票过程如图2所示。对于任一投票角θ,投票距离在一定范围内,并且投票值A(ρ)从0到最大,再减到0。
步骤四,搜索峰值对应的投票距离ρm。当图像中所有圆形目标像素点都向两个累加数组A1(ρ)和A2(ρ)投票后,分别在两个1D累加数组中,搜索最大投票值所对应的投票距离。对于每个投票角θ,所有圆形目标都向累加数组投票后,当投票距离ρ经过圆心坐标时,投票值最大。并且最大投票值等于圆形目标的直径。设最大投票值对应的投票距离为:ρ1 m和ρ2 m。也就是说,当ρ=ρ1 m时,A1(ρ)最大;当ρ=ρ2 m时,A2(ρ)最大。
步骤五,计算圆心坐标。因为当投票距离经过圆心(x0,y0)时,投票值最大。也就是圆心坐标(x0,y0)对应的投票距离为ρm。这样,基于第四步得到的ρ1 m和ρ2 m,构成了两个方程:
ρ1 m=x0·cosθ1+y0·sinθ1
ρ2 m=x0·cosθ2+y0·sinθ2 (1)
解上述方程组,就得到圆心坐标。
x0=(ρ2 m·sinθ1-ρ1 m·sinθ2)/sin(θ1-θ2)
y0=(ρ1 m·cosθ2-ρ2 m·cosθ1)/sin(θ1-θ2) (2)
步骤六,计算圆半径。因为在每个累加数组中,最大投票值为圆形目标的直径。并且当投票距离为ρ1 m和ρ2 m时,投票值A1(ρ)和A2(ρ)最大。所以,圆半径为:
r=A1(ρ1 m)/2或r=A2(ρ2 m)/2 (3)
检测的圆心坐标和半径如图3所示。
以下实例采用仿真实例进行检验。
构造一幅仿真图像,包含一个圆心在(184,97),半径为59的一个圆形对象,并额外添加10%的图像椒盐噪声。
运用本发明利用Hough变换检测圆形目标方法对仿真图像进行处理,得到圆形目标的圆心坐标(x0,y0)为(183.65,98.68),半径为r=60。如图4所示。检测的圆形目标圆心和半径与真值非常接近。圆心误差1.7像素,半径误差1.7%。
本实例在进行圆形目标检测时,一方面降低了算法对内存的需求,提高了算法效率。另一方面,直接在二值图像上检测圆形目标,不需要进行图像边缘检测。
以上是本发明的实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (7)
1.一种利用Hough变换的圆形目标检测方法,包括如下步骤:
步骤一,灰度转换与增强;
步骤二,图像二值化;
步骤三,Hough投票;
步骤四,搜索峰值对应的投票距离ρm;
步骤五,计算圆心坐标;
步骤六,计算圆半径。
2.根据权利要求1所述的一种利用Hough变换的圆形目标检测方法,其特征在于:所述步骤一,对输入的彩色图像进行变换,转换成灰度图像;并做图像平滑处理,降低图像噪声。
3.根据权利要求1所述的一种利用Hough变换的圆形目标检测方法,其特征在于:所述步骤二,运用otsu算法计算阈值,对经过步骤一处理后的灰度图像进行二值化处理,得到二值图像;如果像素点的灰度值大于该阈值,赋予该点灰度值255,否则,赋予该点灰度值0。
4.根据权利要求1所述的一种利用Hough变换的圆形目标检测方法,其特征在于:所述步骤三,从图像空间向参数空间进行投票时,在图像空间进行粗采样Δx=Δy>1,选用直线的投票方程:ρ=x·cosθ+y·sinθ;随机选择两个投票角度θ1和θ2;对图像中的每个黑像素点,运用直线投票方程,分别计算对应的投票距离ρ1和ρ2;分别向两个一维的累加数组A1(ρ)和A2(ρ)投票。
5.根据权利要求1所述的一种利用Hough变换的圆形目标检测方法,其特征在于:所述步骤四,分别在两个1D累加数组中,搜索最大投票值所对应的投票距离ρ;对于每个投票角θ,所有圆形目标都向累加数组投票后,当投票距离ρ经过圆心坐标时,投票值最大;设此时投票距离为ρm,并且最大投票值等于圆形目标的直径。
6.根据权利要求1所述的一种利用Hough变换的圆形目标检测方法,其特征在于:所述步骤五,圆心坐标(x0,y0)对应的投票距离为ρm;基于ρ1 m和ρ2 m构成了两个方程:ρ1 m=x0·cosθ1+y0·sinθ1和ρ2 m=x0·cosθ2+y0·sinθ2,解方程组得到圆心坐标(x0,y0)。
7.根据权利要求1所述的一种利用Hough变换的圆形目标检测方法,其特征在于:所述步骤六,在每个累加数组中,最大投票值为圆形目标的直径;当投票距离为ρm时,投票值A1(ρ1 m)和A2(ρ2 m)最大;根据在ρm处的最大投票值A(ρm)的一半,计算圆半径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711000106.1A CN107564006B (zh) | 2017-10-24 | 2017-10-24 | 一种利用Hough变换的圆形目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711000106.1A CN107564006B (zh) | 2017-10-24 | 2017-10-24 | 一种利用Hough变换的圆形目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107564006A true CN107564006A (zh) | 2018-01-09 |
CN107564006B CN107564006B (zh) | 2020-07-31 |
Family
ID=60986022
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711000106.1A Active CN107564006B (zh) | 2017-10-24 | 2017-10-24 | 一种利用Hough变换的圆形目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107564006B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109365877A (zh) * | 2018-12-20 | 2019-02-22 | 长安大学 | 一种数控机床零件加工自动定位方法 |
CN109658428A (zh) * | 2018-12-14 | 2019-04-19 | 南京工程学院 | 基于Hough变换实现圆的快速检测方法 |
CN110060259A (zh) * | 2019-04-04 | 2019-07-26 | 东华大学 | 一种基于霍夫变换的鱼眼镜头有效区域提取方法 |
CN110097596A (zh) * | 2019-04-30 | 2019-08-06 | 湖北大学 | 一种基于opencv的目标检测系统 |
CN110197488A (zh) * | 2019-05-31 | 2019-09-03 | 中国人民解放军96962部队 | 一种基于ct图像的定位分割系统和方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331885A (zh) * | 2014-10-29 | 2015-02-04 | 电子科技大学 | 一种基于投票线聚类的圆形目标检测方法 |
CN105631846A (zh) * | 2014-10-31 | 2016-06-01 | 陕西高华知本化工科技有限公司 | 一种圆形图形的检测方法 |
-
2017
- 2017-10-24 CN CN201711000106.1A patent/CN107564006B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331885A (zh) * | 2014-10-29 | 2015-02-04 | 电子科技大学 | 一种基于投票线聚类的圆形目标检测方法 |
CN105631846A (zh) * | 2014-10-31 | 2016-06-01 | 陕西高华知本化工科技有限公司 | 一种圆形图形的检测方法 |
Non-Patent Citations (1)
Title |
---|
杨全银: "基于Hough变换的图像形状特征检测", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109658428A (zh) * | 2018-12-14 | 2019-04-19 | 南京工程学院 | 基于Hough变换实现圆的快速检测方法 |
CN109365877A (zh) * | 2018-12-20 | 2019-02-22 | 长安大学 | 一种数控机床零件加工自动定位方法 |
CN110060259A (zh) * | 2019-04-04 | 2019-07-26 | 东华大学 | 一种基于霍夫变换的鱼眼镜头有效区域提取方法 |
CN110097596A (zh) * | 2019-04-30 | 2019-08-06 | 湖北大学 | 一种基于opencv的目标检测系统 |
CN110097596B (zh) * | 2019-04-30 | 2023-06-09 | 湖北大学 | 一种基于opencv的目标检测系统 |
CN110197488A (zh) * | 2019-05-31 | 2019-09-03 | 中国人民解放军96962部队 | 一种基于ct图像的定位分割系统和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107564006B (zh) | 2020-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107564006A (zh) | 一种利用Hough变换的圆形目标检测方法 | |
CN105046252B (zh) | 一种人民币冠字码识别方法 | |
CN107392963B (zh) | 一种用于软式自主空中加油的仿鹰眼运动目标定位方法 | |
CN104167003A (zh) | 一种遥感影像的快速配准方法 | |
CN106909941A (zh) | 基于机器视觉的多表字符识别系统及方法 | |
CN104568986A (zh) | 基于surf算法的遥控器面板印刷缺陷自动检测方法 | |
CN106709500B (zh) | 一种图像特征匹配的方法 | |
CN106339707B (zh) | 一种基于对称性特征的仪表指针图像识别方法 | |
CN106022315A (zh) | 一种用于虹膜识别的瞳孔中心定位方法 | |
CN106156684B (zh) | 一种二维码识别方法及装置 | |
CN103426186A (zh) | 一种改进的surf快速匹配方法 | |
CN106407924A (zh) | 基于路面特征的双目道路识别检测方法 | |
CN103400384A (zh) | 结合区域匹配和点匹配的大视角图像匹配方法 | |
CN104951940A (zh) | 一种基于掌纹识别的移动支付验证方法 | |
CN101303768A (zh) | 圆形标志点在摄像机透视投影变换时圆心偏差的修正方法 | |
CN108256467B (zh) | 一种基于视觉注意机制和几何特征的交通标志检测方法 | |
CN107506688B (zh) | Harris角点检测影像金字塔掌纹ROI识别方法 | |
CN107092876A (zh) | 基于Retinex与S‑SIFT特征结合的低照度车型识别方法 | |
CN102938147A (zh) | 一种基于快速鲁棒特征的低空无人机视觉定位方法 | |
CN106327534A (zh) | 一种基于定位块的轮胎内壁纹理识别方法 | |
CN106503748A (zh) | 一种基于s‑sift特征和svm训练器的车型识别算法 | |
CN110245626B (zh) | 一种在虹膜图像中精确检测眼睫毛图像的方法 | |
CN106127258A (zh) | 一种目标匹配方法 | |
CN110443295A (zh) | 改进的图像匹配与误匹配剔除算法 | |
CN110909772B (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 |