CN105678763A - 一种稳定的直线/圆特征检测方法 - Google Patents

一种稳定的直线/圆特征检测方法 Download PDF

Info

Publication number
CN105678763A
CN105678763A CN201610006030.2A CN201610006030A CN105678763A CN 105678763 A CN105678763 A CN 105678763A CN 201610006030 A CN201610006030 A CN 201610006030A CN 105678763 A CN105678763 A CN 105678763A
Authority
CN
China
Prior art keywords
detection
point
straight line
line
circle
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
CN201610006030.2A
Other languages
English (en)
Other versions
CN105678763B (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.)
Bozhong Suzhou Precision Industry Technology Co Ltd
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201610006030.2A priority Critical patent/CN105678763B/zh
Publication of CN105678763A publication Critical patent/CN105678763A/zh
Application granted granted Critical
Publication of CN105678763B publication Critical patent/CN105678763B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20112Image segmentation details
    • G06T2207/20164Salient point detection; Corner detection

Landscapes

  • Image Analysis (AREA)

Abstract

本发明提供了一种直线/圆特征检测方法,首先设定基准点及初始检测区域、通过模板匹配对检测区域进行重新定位、再利用位置补正调整检测区域,通过提取测量线方向数据、投影线范围内测量线局部噪声尺度计算和测量线上局部尺度的滤波操作计算得到噪声的方差估计值,然后计算出一维数据的边缘点,再通过拟合内点及噪声点分割与筛选剔除检测到点中的噪声点,保留直线/圆上的内点,最后利用最小二乘方法对检测到内点拟合出直线/圆。本发明的方法能够在强噪声图像中稳定地提取边缘特征点并拟合为直线/圆,由于通过模板匹配进行检测区域及ROI区域设定,直线/圆只在ROI区域进行,具有非常好的计算实时性,在工业机器视觉的测量技术中将有重要的应用。

Description

一种稳定的直线/圆特征检测方法
技术领域
本发明涉及图像识别技术领域,尤其涉及一种直线/圆特征检测方法。
背景技术
直线/圆是图像中的重要特征,直线/圆检测是计算机视觉和机器视觉识别中的最为重要的工具之一,在机器视觉中有着重要的应用,例如尺寸测量、参考坐标系建立、网格特征识别、单据表格识别、刻度线检测、Mark点识别半径、直径测量等直线/圆类特征的检测。
在图像处理中,通常使用的方法集中在Hough变换和Radon变换,各种不同的算法重点在于解决检测的精度和速度问题(参见非专利文献1:JiqiangSong,MichaelRLyu.AHoughTransformbasedLineRecognitionMethodUtilizingbothParametersSpaceandImageSpace,PatternRecognition,38(4):539-552,2005.王立等,Radon变换在低信噪比图像中的线段检测,红外与激光工程,2003,32(2):163-166)。非专利文献2(史册等,实时图像处理中一种快速的直线检测算法,浙江大学学报,1999,33(5):482-486)提出了一种基于斑点的快速直线检测算法。非专利文献3(孙涵等,一种提出了一种基于链码的直线检测算法快速实用的直线检测算法,计算机应用研究,2006,2:256-260)。专利文献1(中国专利公开号CN101645172A)公开了一种基于Hough变换的直线检测算法;专利文献2(中国专利公开号CN103198443A)公开了一种图像处理装置和直线检测方法;专利文献3(中国专利公开号CN103914830A)公开了一种基于梯度直方图特征集的直线方法,根据梯度的局部直线特性,根据各个所述像素点的梯度直方图特征集确定待检测直线的至少一个备选方向;根据所述备选方向确定所述待检测直线的精确方向和位置。专利文献4(中国专利公开号CN104331876A)及专利文献5(中国专利公开号CN104408456A)公开了一种基于Hough变换的直线检测、图像处理的方法及装置;专利文献6(中国专利公开号CN102482027A)公开了一种基于Hough变换的圆检测方法,事先估算出圆的半径,在半径小范围利用Hough变换检测圆,该算法虽能提高效率,但稳定性不高;专利文献7(中国专利公开号CN104036514A)公开了一种基于直方图峰值搜索的圆检测方法,专利文献8(中国专利公开号CN1039032824A)公开了一种基于最小二乘的多圆检测方法,该专利中如何获得圆边缘点并没有公开;专利文献9(中国专利公开号CN103295227A)公开了一种基于梯度方向分割的圆检测方法。
在工业应用中,由于产品的形态多样,工厂环境复杂多变,或者由于产品表面的变化及反光等,使得产品图像质量变差,图像中含有大量的噪声、阴影等不稳定因素,采用已有的基于Hough变换、Radon变换或链码的方法,都不能稳定地检测出直线/圆特征,同时这类算法耗时较长,不能在工业机器视觉中进行实时检测。在实际机器视觉系统中通过依靠图像中的显著特征进行定位,然后利用直线/圆检测区域和显著特征区域的几何位置关系对后续的图像进行直线/圆检测区域的定位,在直线/圆检测区域中利用ROI区域进行直线/圆特征的检测,得到亚像素的边缘点后再利用最小二乘等方法拟合出直线/圆,但在强噪声图像中,经常会误检测到错误的边缘点,进行特征拟合时会对结果产生影响,造成检测误差。有的算法采用枚举的方式排出异常点(不在直线/圆上的点),但这种算法效率不高。
发明内容
本发明的目的在于提供一种在强噪声图像中稳定的直线/圆检测方法,为工业恶例环境中机器视觉测量功能提供了重要的工具。
为达上述目的,本发明通过以下技术方案实现:
一种直线/圆特征检测方法,包括以下步骤:步骤①:基准点及初始检测区域设置;步骤②:通过模板匹配对检测区域基准点重新定位;步骤③:参照基准点利用位置补正调整初始设置检测区域;步骤④:区域范围内待拟合特征检测;其中,步骤④分为:
⑤一维数据采样:选择二维图像中的一个区域并提取区域范围内的数据,在每个位置根据数据来源的特性,决定各位置滤波与否或是滤波的程度,通过分析数据变化的快慢判断某位置是否为边缘点,具体算法为:(1)提取测量线方向数据;(2)投影线范围内测量线局部噪声尺度计算;(3)测量线上局部尺度的滤波操作;
⑥一维边界点检测:边界点的定义是灰度变化强烈的位置,边缘滤波算法采用一阶滤波或二阶滤波算法,其中,一阶滤波算法的判断标准是选择梯度幅值极大值的位置,二阶滤波算法是选择导数过零点位置;
⑦拟合数据点提取:通过拟合内点及噪声点分割与筛选剔除检测到点中的噪声点,保留直线/圆上的内点;
⑧直线/圆拟合:利用最小二乘算法解决拟合参数的问题。
作为本发明的进一步改进,所述检测区域提供特征检测时的部分参数,包括:测量线方向、起始位置、长度、检测方向、检测段高度;同时在该区域还为当前图像Cimg的检测区域提供了基准位置。
作为本发明的进一步改进,直线拟合采用最小二乘拟合算法,圆拟合采用非线性的最小二乘拟合算法或牛顿迭代算法。
附图说明
图1是本发明的方法流程图;
图2(a)是直线检测区域设置示意图;
图2(b)是圆检测区域设置示意图;
图3是图像中显著特征定位示意图;
图4是模板匹配定位及位置补正得到新的特征检测区域示意图;
图5是一维边缘模型示意图;
图6是直线特征检测时一维数据采样示意图;
图7是测量线投影数据曲线示意图;
图8是一维测量线的数据不同处理方式结果示意图;
图9是传统的直线检测方法与本发明的直线检测方法的检测结果对比示意图,其中,图9(a)和图9(c)是传统的直线检测方法的检测点和直线拟合的结果示意图,图9(b)和图9(d)是本发明的直线检测方法的检测点和直线拟合的结果示意图。
具体实施方案
下面通过具体实施方式结合附图对本发明作进一步详细说明。
如附图1所示,本发明的一种稳定的直线/圆特征检测方法,包括以下步骤:步骤①:基准点及初始检测区域设置;步骤②:检测区域基准点重新定位;步骤③:参照基准点调整初始设置检测区域;步骤④:区域范围内待拟合特征检测;其中,步骤④分为:⑤一维数据采样、⑥一维边界点检测、⑦拟合数据点提取、⑧直线/圆拟合。
在直线/圆检测流程中步骤①中,首先在初始图像(参考图像)中设置特征检测区域,如附图2(a)和附图2(b)所示,检测区域提供特征检测时的部分参数,包括:测量线方向、起始位置、长度、检测方向、检测段高度。同时在该区域还为当前图像Cimg的检测区域提供了基准位置,根据步骤②和步骤③中的变换关系算出基准位置在图像Cimg中的新位置。
在步骤②中利用模板匹配方法得到图像中显著特征的位置FP,检测区域的基准点P设置在参考图像上,Fp和P点可以建立变换关系(R,t),如附图3所示。当图像改变时(当前图像),基准点的位置会发生改变,产生新的基准点P′。P′的计算方法在步骤③中得到。
对于后续输入的图像(当前图像)利用步骤③中位置补正关系将检测区域重新设定。基准点位置的改变伴随着产生坐标变换矩阵的改变,通过参考图像和当前图像的矩阵变换调整检测区域,得到检测区域在当前图像的位置及方向。如附图4所示,图中检测区域的位置和姿态都发生了变化,首先利用模板匹配定位方法找到图像中的显著特征及位置Fp,再利用参考图像中的变化关系(R,t)计算出新的检测区域的基准位置P′。
一维测量数据采样及处理
假设数字图像的组成公式为I(x,y)=G(x,y)-N(x,y),图像范围内各点都是由原始模拟图像G(x,y)及噪声图像N(x,y)叠加得到。相机拍摄时光照变化、相机芯片感光因素、不同芯片之间能量分布等各种因素都会体现在噪声图像上,因此噪声图像的成因一般没有简单的模型可以得到。边缘提取过程当中,噪声因素会干扰最终的结果。多数图像处理的分析都将高斯噪声模型替代噪声图像,通过处理范围内的高斯信号抑制该区域内的噪声图像。
一维边缘检测的数据模型是选择二维图像中的一个区域并提取区域范围内的数据,通过分析数据变化的快慢判断某位置是否为边缘点。一维投影数据来源于二维图像,在得到有效数据的同时也得到了噪声信号。附图5表示标准的一维边界点模型,边界点在位置R处。该模型没有受噪声影响,可以不使用滤波操作。使用高斯滤波或其它滤波(如双边滤波)处理附图5所示的模型时,在图像灰度阶跃处(位置R)会弱化边界点的梯度幅值。
一维数据生成过程中,在得到有效数据(附图5中R处附近)的同时会带来噪声信息。传统的滤波操作在平滑噪声区域数据及特征区域数据的同时,降低特征区域梯度强度并改变特征区域的数据走向。本发明所采用的算法的特点是在每个位置根据数据来源的特性,决定各位置滤波与否或是滤波的程度。
附图6的表示图像在生成一维数据时的步骤,粗线框是设定的检测范围,虚线是测量线(测量数据采样点),大空心箭头表示投影方向垂直于测量线方向。边界点检测方向和测量线方向平行。附图6中圆点左侧是强噪声区域,右侧是特征数据伴有弱噪声。本发明所采用的算法在维持特征信号的前提下,减弱噪声对结果的干扰。测量线上每个点的数据由检测范围内投影方向上的内容组成,决定滤波强弱及是否滤波的是投影线所处的环境。附图6中圆点左侧垂直测量线的投影线为噪声区域,且噪声强度较大,需要加强滤波;在圆点右侧附近是特征数据区域,为保持边界点的准确性,测量线保持原始数据;右侧偏右的弱噪声区域弱化滤波处理。
一维投影数据来源于二维图像,如同图像模型表示为i(x)=f(x)+n(x,σx)。σx是测量线上各点的噪声标准方差,由测量线上投影线范围内计算得到。噪声方差的估计值按照统计方法直接得到。本发明所采用的算法流程如下:
(1)提取测量线方向数据;
(2)投影线范围内测量线局部噪声尺度计算;
(3)测量线上局部尺度的滤波操作。
附图7是测量线投影数据曲线示意图,附图8是一维测量线的数据不同处理方式结果示意图,其中附图8中的点线即为附图7的曲线,虚线是全局滤波方式,实线是局部尺度滤波方式,滤波方法是双边滤波。局部尺度滤波在噪声波动强烈范围内较大尺度的平滑处理原始数据,在噪声较小区域或是无噪声范围内能更大程度的保持处理结果相对原始数据不变。
尺度滤波处理后的数据后期利用一阶导数操作,选择|f′g|≥fT设定为边界点。数据经过滤波操作降低噪声区的一阶导数幅值,更便于剔除噪声信号。
一维数据的边缘点检测
边缘点的定义是灰度变化强烈的位置,如附图5的R坐标处。边缘滤波算法主要分为一阶滤波和二阶滤波算法。一阶滤波算法的判断标准是选择梯度幅值极大值的位置,如sobel算法ei=|f′*h|i=|f*h′|i;二阶滤波算法是选择导数过零点位置,如Laplace算法pi=f″*h=f*h″,其中h是滤波函数,如高斯函数。当边界点的模糊程度较大时,一阶导数方式或二阶导数方法很难得到准确的单一边缘点信息。局部尺度的一阶算法及二阶算法能够针对不同模糊程度的边缘位置使用合适的尺度数值,得到单一准确数值。
高斯信号的一阶导数为 g 1 ′ ( x , σ 1 ) = - x 2 πσ 1 4 e - x 2 / 2 σ 1 2 , | r x | = | g x ′ ( x , σ 1 ) * I ( x ) | = r θ M , 尺度阈值计算方法如下。设参数αP为逐点意义数据,αI为全局意义数据,两者关系为αP=1-(1-αI)1/n,设定αI=0.05,n是图像尺寸。此外需使用命题1:假设一个线性变化sn是定义域的方差,线性变化的方差为SL=||L||2sn。设概率分布函数pU,V=f(U),得到如下概率公式:
p V ( v ) = p U ( f - 1 ( v ) ) | d d v f - 1 ( v ) | - - - ( 1 )
图像噪声设定为高斯噪声,标准方差为sn
P U ( u ) = 2 2 πs 1 e - u 2 / 2 s 1 2 , u ∈ [ 0 , ∞ ) - - - ( 2 )
其中s1=||g1(x,σ1||2sn(上述命题1设定)。由公式(1)f(u)=u2,推导并得公式:
p V ( v ) = 1 2 πvs 1 e - v / 2 s 1 2 - - - ( 3 )
P V 1 + V 2 ( v ) = p [ ( r θ M ) 2 = v ] = ∫ 0 v p V ( v ′ ) p V ( v - v ′ ) dv ′ - - - ( 4 )
积分结果如下:
p V 1 + V 2 ( v ) = 1 2 s 1 2 e - v / 2 s 1 2 - - - ( 5 )
为保证αP的阈值,设定积分下限值c1,该值满足如下等式:
∫ c 1 2 ∞ p V 1 + V 2 ( v ) d v = α p , c 1 = s 1 - 2 ln ( α p ) - - - ( 6 )
其中s1=||g1(x,σ1)||2αP=2.0×10-7,得到如下等式:
c 1 ( σ 1 ) = 1.1 s n σ 1 2 - - - ( 7 )
式(7)的图像噪声信号方差sn由图像内容计算,模糊信号方差(σ1)为未知数。计算一阶导数时也需要使用模糊信号方差。在应用中,为计算满足下列条件的模糊信号方差,在系列参数中枚举匹配符合要求数据。
σ ^ 1 ( x ) = i n f { σ 1 : r 1 θ M ( x , σ 1 ) > c 1 ( σ 1 ) } , σ 1 ∈ { 0.5 , 1 , 2 , 4 , 8 , 16 } - - - ( 8 )
同理,二阶导数的局部尺度数值使用上述方式计算得到,其中
拟合内点及噪声点分割与筛选
集合P:={pi|pi:=(xi,yi)}是待拟合的点集,拟合模型(如直线ax-by+c=0,圆(x-a)2+(y-b)2=r2)参数为点集在函数变化后的误差为:
r i = F ( p i , θ ^ ) - - - ( 9 )
其中F(·)是拟合函数相关的用来计算误差的函数。后期需要对误差集合{|ri|}从小到大排序。
设定内点分割尺度(比例)参数s,所有符合下面表达式的都属于内点集合:
| r i s | < E - - - ( 10 )
阈值参数E设置为2.5能包含98%的符合高斯模型误差的内点。
参数是第J次拟合参数(θJ)的第i个绝对误差,nJ是内点个数,尺度参数可以用迭代K阶尺度估计(IterativeScaleKEstimation,ISKE)
s ^ K J : = | r ~ i J | &Theta; - 1 ( 1 2 ( 1 + k J ) ) , k J : = K / n J - - - ( 11 )
其中是高斯密度函数的正态分布累积函数;K是一个预设定的固定参数,数值和输入点集的噪声信号相关,在0到1之间,式(11)的关键是如何估计内点个数nJ。迭代K阶尺度估计的算法设定nJ的初始值等于输入点集的势,根据s的特性迭代计算尺度数值,直至得到收敛的尺度数值。集合表示第J次拟合参数的第t次迭代,该集合收敛且单调递减,迭代K阶尺度估计算法描述如下所示:
K值的选择与输入点集合的特性有关:噪声点的比率较高时,设置较小K值;反之,设置较大的数据。默认K值设置为0.1。当该值设置较小时,存在原本属于内点的数据被分割到噪声集合的可能性。在得到内点集合后,针对噪声集合中的数据进一步做筛选,将绝对误差小且不会影响内点整体波动特性的数据重新分配到内点集合。
在挑选最佳拟合结果过程中,会选择不同组合的拟合参数,根据不同的组合得到初始拟合数据,最终的结果要选择权重最大的一组。权重的数据来源由拟合模型、内点集合的误差。
设定集合数据是第j次拟合参数的所有内点的误差,设计估计权重函数的定义如下:
f ^ K N , &theta; ^ j ( r ) : = 1 n &Sigma; i = 1 n 1 h ( &theta; ^ j ) K N &lsqb; r - r i ( &theta; j ) h ( &theta; ^ j ) &rsqb; - - - ( 12 )
式(12)中的叶帕涅奇尼科夫核函数(theEpanechnikovkernelKNE(r))及带宽函数的估计值的表达式如下:
KN E ( r ) : = 3 4 ( 1 - | | r | | 2 ) , | | r | | &le; 1 0 , | | r | | > 1 - - - ( 13 )
h ^ ( &theta; ^ j ) = &lsqb; 243 &Integral; - 1 1 K N ( r ) 2 d r 35 n &Integral; - 1 1 r 2 K N ( r ) d r &rsqb; 1 5 s ^ K ( &theta; ^ j ) - - - ( 14 )
得到权重函数如下:
w ^ j &Proportional; 1 n &Sigma; i = 1 n K N ( r i ( &theta; ^ j ) / h ) - - - ( 15 )
拟合的权重越大,表示拟合度越好。以上的权重注重的是内点的误差特性,会使某些波动较大的内点数据产生较大权重,线性度较好的内点集合权重较小。设内点集合的误差集合数据的方差是设定相关拟合权重的另一种模式为:
ω=ωσj(17)
最后的输出结果选择为权重最大的一组内点。
特征拟合的输入集是上述方式挑选出来的所有内点,利用最小二乘算法解决拟合参数的问题。直线拟合的最小二乘是线性方程,圆拟合采用非线性的最小二乘拟合算法或牛顿迭代算法。
为验证本发明的有效性进行了大量的测试,附图9为测试的部分实例,是实际工业机器视觉应用采集的图像,从图中可以看出图像质量较差,直线边缘部分模糊,特征不明显。附图9(a)和附图9(c)为传统的直线检测方法的检测点和拟合直线,附图9(b)和附图9(d)为本发明的直线检测方法的检测点及拟合直线,从图中可以看出本发明在图像质量很差时仍能够稳定地检测到需要的直线。
本发明的方法,首先设定基准点及初始检测区域、通过模板匹配对检测区域进行重新定位、再利用位置补正调整检测区域,通过提取测量线方向数据、投影线范围内测量线局部噪声尺度计算和测量线上局部尺度的滤波操作计算得到噪声的方差估计值,然后计算出一维数据的边缘点,再通过拟合内点及噪声点分割与筛选剔除检测到点中的噪声点,保留直线/圆上的内点,最后利用最小二乘方法对检测到内点拟合出直线/圆。本发明提出的方法能够在强噪声图像中稳定地提取边缘特征点并拟合为直线/圆的方法,由于通过模板匹配进行检测区域及ROI区域设定,直线/圆只在ROI区域进行,具有非常好的计算实时性,在工业机器视觉的测量技术中将有重要的应用。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (6)

1.一种直线/圆特征检测方法,其特征在于:所述方法包括以下步骤:步骤①:基准点及初始检测区域设置;步骤②:通过模板匹配对检测区域基准点重新定位;步骤③:参照基准点利用位置补正调整初始设置检测区域;步骤④:区域范围内待拟合特征检测;其中,步骤④分为:
⑤一维数据采样:选择二维图像中的一个区域并提取区域范围内的数据,在每个位置根据数据来源的特性,决定各位置滤波与否或是滤波的程度,通过分析数据变化的快慢判断某位置是否为边缘点,具体算法为:(1)提取测量线方向数据;(2)投影线范围内测量线局部噪声尺度计算;(3)测量线上局部尺度的滤波操作;
⑥一维边界点检测:边界点的定义是灰度变化强烈的位置,边缘滤波算法采用自适应的局部尺度参数实现一阶滤波或二阶滤波算法,其中,一阶滤波算法的判断标准是选择梯度幅值极大值的位置,二阶滤波算法是选择导数过零点位置;
⑦拟合数据点提取:通过拟合内点及噪声点分割与筛选剔除检测到点中的噪声点并在剔除的噪声点中取回符合要求的边界点,两者合并保留为直线/圆上的内点;
⑧直线/圆拟合:利用最小二乘算法解决拟合参数的问题。
2.根据权利要求1所述的方法,其特征在于:所述检测区域提供特征检测时的部分参数,包括:测量线方向、起始位置、长度、检测方向、检测段高度;同时在该区域还为当前图像Cimg的检测区域提供了基准位置。
3.根据权利要求1所述的方法,其特征在于:使用内点集合作输入,直线拟合采用最小二乘拟合算法,圆拟合采用非线性的最小二乘拟合算法或牛顿迭代算法。
4.一种直线/圆特征检测装置,其特征在于:所述装置包括:
设置模块,用于基准点及初始检测区域设置;
重定位模块,用于通过模板匹配对检测区域基准点重新定位;
调整模块,用于参照基准点利用位置补正调整初始设置检测区域;
拟合模块,用于区域范围内待拟合特征检测;其中,所述拟合模块包括:
一维数据采样子模块,用于选择二维图像中的一个区域并提取区域范围内的数据,在每个位置根据数据来源的特性,决定各位置滤波与否或是滤波的程度,通过分析数据变化的快慢判断某位置是否为边缘点,具体为:(1)提取测量线方向数据;(2)投影线范围内测量线局部噪声尺度计算;(3)测量线上局部尺度的滤波操作;
一维边界点检测子模块,用于边界点的定义是灰度变化强烈的位置,边缘滤波算法采用自适应的局部尺度参数实现一阶滤波或二阶滤波算法,其中,一阶滤波算法的判断标准是选择梯度幅值极大值的位置,二阶滤波算法是选择导数过零点位置;
拟合数据点提取子模块,用于通过拟合内点及噪声点分割与筛选剔除检测到点中的噪声点并在剔除的噪声点中取回符合要求的边界点,两者合并保留为直线/圆上的内点;
直线/圆拟合子模块,用于利用最小二乘算法解决拟合参数的问题。
5.根据权利要求4所述的装置,其特征在于:所述检测区域提供特征检测时的部分参数,包括:测量线方向、起始位置、长度、检测方向、检测段高度;同时在该区域还为当前图像Cimg的检测区域提供了基准位置。
6.根据权利要求4所述的装置,其特征在于:使用内点集合作输入,直线拟合采用最小二乘拟合算法,圆拟合采用非线性的最小二乘拟合算法或牛顿迭代算法。
CN201610006030.2A 2016-01-05 2016-01-05 一种稳定的直线/圆特征检测方法及装置 Active CN105678763B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610006030.2A CN105678763B (zh) 2016-01-05 2016-01-05 一种稳定的直线/圆特征检测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610006030.2A CN105678763B (zh) 2016-01-05 2016-01-05 一种稳定的直线/圆特征检测方法及装置

Publications (2)

Publication Number Publication Date
CN105678763A true CN105678763A (zh) 2016-06-15
CN105678763B CN105678763B (zh) 2018-11-27

Family

ID=56298953

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610006030.2A Active CN105678763B (zh) 2016-01-05 2016-01-05 一种稳定的直线/圆特征检测方法及装置

Country Status (1)

Country Link
CN (1) CN105678763B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107230212A (zh) * 2017-05-08 2017-10-03 武汉科技大学 一种基于视觉的手机尺寸的测量方法及系统
CN107909563A (zh) * 2017-10-09 2018-04-13 常熟理工学院 一种基于模板的快速Hough变换直线检测方法
CN109612390A (zh) * 2018-12-17 2019-04-12 江南大学 基于机器视觉的大尺寸工件自动测量系统
WO2020133538A1 (zh) * 2018-12-29 2020-07-02 深圳配天智能技术研究院有限公司 一种圆形特征检测方法、处理系统及具有存储功能的装置
CN112085759A (zh) * 2020-09-07 2020-12-15 凌云光技术股份有限公司 一种基于大数据的直线拟合方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101685251A (zh) * 2008-09-23 2010-03-31 中芯国际集成电路制造(上海)有限公司 镜头校验方法及系统
US20100098314A1 (en) * 2008-10-22 2010-04-22 Carestream Health, Inc. Tube detection in diagnostic images
CN103091331A (zh) * 2013-01-11 2013-05-08 华中科技大学 一种rfid天线毛刺和污点缺陷的视觉检测系统及方法
CN103425988A (zh) * 2013-07-03 2013-12-04 江南大学 一种具有圆弧几何基元的实时定位与匹配方法
CN104464079A (zh) * 2014-12-29 2015-03-25 北京邮电大学 基于模板特征点及其拓扑结构的多币种面值识别方法
CN104568986A (zh) * 2015-01-26 2015-04-29 中国科学院半导体研究所 基于surf算法的遥控器面板印刷缺陷自动检测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101685251A (zh) * 2008-09-23 2010-03-31 中芯国际集成电路制造(上海)有限公司 镜头校验方法及系统
US20100098314A1 (en) * 2008-10-22 2010-04-22 Carestream Health, Inc. Tube detection in diagnostic images
CN103091331A (zh) * 2013-01-11 2013-05-08 华中科技大学 一种rfid天线毛刺和污点缺陷的视觉检测系统及方法
CN103425988A (zh) * 2013-07-03 2013-12-04 江南大学 一种具有圆弧几何基元的实时定位与匹配方法
CN104464079A (zh) * 2014-12-29 2015-03-25 北京邮电大学 基于模板特征点及其拓扑结构的多币种面值识别方法
CN104568986A (zh) * 2015-01-26 2015-04-29 中国科学院半导体研究所 基于surf算法的遥控器面板印刷缺陷自动检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
范壮: "机械零件尺寸视觉检测系统的研究", 《中国优秀硕士学位论文全文数据库•信息科技辑》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107230212A (zh) * 2017-05-08 2017-10-03 武汉科技大学 一种基于视觉的手机尺寸的测量方法及系统
CN107230212B (zh) * 2017-05-08 2020-04-17 武汉科技大学 一种基于视觉的手机尺寸的测量方法及系统
CN107909563A (zh) * 2017-10-09 2018-04-13 常熟理工学院 一种基于模板的快速Hough变换直线检测方法
CN107909563B (zh) * 2017-10-09 2020-03-31 常熟理工学院 一种基于模板的快速Hough变换直线检测方法
CN109612390A (zh) * 2018-12-17 2019-04-12 江南大学 基于机器视觉的大尺寸工件自动测量系统
WO2020133538A1 (zh) * 2018-12-29 2020-07-02 深圳配天智能技术研究院有限公司 一种圆形特征检测方法、处理系统及具有存储功能的装置
CN111801709A (zh) * 2018-12-29 2020-10-20 深圳配天智能技术研究院有限公司 一种圆形特征检测方法、处理系统及具有存储功能的装置
CN111801709B (zh) * 2018-12-29 2023-12-22 深圳配天机器人技术有限公司 一种圆形特征检测方法、处理系统及具有存储功能的装置
CN112085759A (zh) * 2020-09-07 2020-12-15 凌云光技术股份有限公司 一种基于大数据的直线拟合方法及装置
CN112085759B (zh) * 2020-09-07 2023-11-10 凌云光技术股份有限公司 一种基于大数据的直线拟合方法及装置

Also Published As

Publication number Publication date
CN105678763B (zh) 2018-11-27

Similar Documents

Publication Publication Date Title
Jiang et al. HDCB-Net: A neural network with the hybrid dilated convolution for pixel-level crack detection on concrete bridges
CN105678763A (zh) 一种稳定的直线/圆特征检测方法
CN105913415A (zh) 一种具有广泛适应性的图像亚像素边缘提取方法
CN103759676A (zh) 一种工件表面粗糙度非接触式的检测方法
CN105787486B (zh) 一种基于图像处理的钢梁裂纹检测方法
Jiale et al. Reading recognition method of analog measuring instruments based on improved hough transform
CN106530271B (zh) 一种红外图像显著性检测方法
CN106228161A (zh) 一种指针式表盘自动读数方法
CN110751604A (zh) 一种基于机器视觉的钢管焊缝缺陷在线检测方法
CN109993154A (zh) 变电站单指针式六氟化硫型仪表智能识别方法
CN107092863A (zh) 一种变电站巡检机器人的指针式仪表读数识别方法
CN111127417B (zh) 一种基于sift特征匹配和改进ssd算法的印刷缺陷检测方法
CN110210448A (zh) 一种智能人脸皮肤老化程度的识别与评估方法
CN104778684A (zh) 钢材表面异质型缺陷的自动测量、表征分类方法及其系统
CN111368825A (zh) 一种基于语义分割的指针定位方法
CN105930852A (zh) 一种气泡图像识别方法
CN107704840A (zh) 一种基于深度学习的遥感图像道路检测方法
CN111415349A (zh) 一种基于图像处理技术的涤纶长丝毛丝检测方法
CN106056161A (zh) 一种针对平面旋转目标的视觉检测方法
CN108645490A (zh) 基于图像处理技术的冷水水表视觉检测系统及检测方法
CN110363196A (zh) 一种倾斜文本的文字精准识别的方法
CN113947583A (zh) 基于深度学习二维时间序列图像的焊缝无损检测方法
CN114092728A (zh) 一种基于深度学习的指针式仪表智能识别方法及系统
CN106022354A (zh) 基于svm的图像mtf测量方法
Aferi et al. Cotton texture segmentation based on image texture analysis using gray level co-occurrence matrix (GLCM) and Euclidean distance

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20160819

Address after: The 215200 West Jiangsu city of Suzhou province Wujiang economic and Technological Development Zone No. 666

Applicant after: Suzhou Bozhong Precision Industry Technology Co., Ltd.

Address before: 518000 building C, Innovation Research Institute, Nanshan District hi tech Zone, Guangdong, Shenzhen 1-6

Applicant before: Wu Xiaojun

CB02 Change of applicant information
CB02 Change of applicant information

Address after: 215200 No. 666 Hu Xin Road, Wujiang economic and Technological Development Zone, Suzhou, Jiangsu

Applicant after: Bo Seiko Polytron Technologies Inc

Address before: 215200 No. 666 Hu Xin Road, Wujiang economic and Technological Development Zone, Suzhou, Jiangsu

Applicant before: Suzhou Bozhong Precision Industry Technology Co., Ltd.

GR01 Patent grant
GR01 Patent grant