CN108830868A - 一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法 - Google Patents
一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法 Download PDFInfo
- Publication number
- CN108830868A CN108830868A CN201810354070.5A CN201810354070A CN108830868A CN 108830868 A CN108830868 A CN 108830868A CN 201810354070 A CN201810354070 A CN 201810354070A CN 108830868 A CN108830868 A CN 108830868A
- Authority
- CN
- China
- Prior art keywords
- edge
- point
- formula
- iteration
- snake model
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 241000270295 Serpentes Species 0.000 title claims abstract description 28
- 238000005457 optimization Methods 0.000 claims abstract description 8
- 230000008569 process Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000005452 bending Methods 0.000 claims description 4
- 230000000295 complement effect Effects 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 10
- 238000001514 detection method Methods 0.000 abstract description 7
- 238000007781 pre-processing Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 238000003708 edge detection Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000003455 independent Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/181—Segmentation; Edge detection involving edge growing; involving edge linking
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30164—Workpiece; Machine component
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Quality & Reliability (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法,用于定位轴类工件图像中的R角中心和求解R角半径,属于机器视觉精密测量技术研究领域。对框选的多边形图像感兴趣区域进行图像预处理,消除原始图像的孤立噪声点,提高图像的对比度;通过Canny算子进行边缘提取,实现边缘的初定位;利用Snake模型滤除边缘上较大的毛刺产生的伪边缘点,并在稀疏区域插值边缘点,即对原始边缘进行Snake优化;采用Zernike正交矩对边缘进行精定位,获取精确的亚像素级边缘;通过迭代极性转换回归方法对亚像素边缘点进行圆拟合,获得精确的R角的圆心位置和半径。该检测方法解决了由于附着物干扰降低轴类工件R角检测的精确度问题,以及短圆弧拟合的可靠性问题。
Description
技术领域
本发明涉及一种机器视觉的测量技术及轴类工件测量系统,具体涉及一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法,属于机器视觉精密测量技术研究领域。
背景技术
轴类工件是工业中极为常见也是非常重要的旋转类工件,其综合形状误差精度不仅直接影响机械的运动性能、使用寿命,而且对减少能源消耗和环境污染等有重要影响。R角是轴类工件误差检测项中的重点和难点,原因在于轴类工件的R角是一个曲面,即使在远心镜头和平行背光的配合下,也无法使曲面上的每一个点完全平行地投射到工业相机的CCD感光芯片上,再加上工件本身的附着物干扰,使R角定位的精确度和准确度受到极大的影响。
R角是指通过机器视觉获取的二维图像上的一段圆弧,工业中通常用接触式的测量方式检测其参数,如R规测量法,其原理是通过光隙法测量圆弧半径,测量时必须使R规的测量面与工件圆弧完全的紧密接触,当测量面与工件的圆弧中间没有间隙时,工件的圆弧尺寸则为此时R规上对应的数字。由于该方法是目测,故准确度不高。基于图像的非接触式无损测量方式,测量结果更为精确,也更易实现工业自动化。R角是圆的一部分,所以R角定位与圆定位的算法是一致,而圆定位的精确与否,主要由边缘提取的准确性和圆拟合的精确性两方面决定。边缘检测的算法主要分为差分边缘检测算法和基于微分算子的边缘检测算法,其中,基于微分算子的Canny边缘检测算子在数字图像处理中得到广泛应用。传统的圆拟合算法主要有3类:基于Hough变换的圆拟合算法、基于不变矩的方法和基于最小二乘法。其中,当随机误差为正态分布时,最小二乘法是由最大似然法推出的一个最优估计技术,可使测量误差的平方和最小,因此也被认为是从一组测量值中获得未知量的理想方法之一。
发明内容
本发明的目的是针对轮廓附着物对视觉测量中R角点检测的干扰以及短圆弧拟合不稳定的问题,设计一种Snake模型和迭代极性转换回归结合的圆弧拟合方法,具体是指一种通过工业线阵相机和平行背光同步运动采集轴类工件图像,经图像处理技术精确定位R角中心与求解R角半径的方法,能够可靠且精确提取感兴趣区域内的R角中心的方法,以此测量轴类工件的部分尺寸参数。
本发明的技术方案如下:
一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法,其特征在于,包含以下步骤:
步骤1、对图像进行预处理;
步骤2、利用Canny算子对该预处理区域进行边缘提取;
步骤3、利用Snake模型对所述区域的边缘优化,去除伪边缘点,并插补邻近边缘点;
步骤4、利用Zernike正交矩获取亚像素级边缘;
步骤5、通过迭代极性转换回归实现R角中心的精确定位。
其进一步的技术方案为,所述步骤3包括:
步骤3.1、通过对能量函数迭代,得到最优轮廓;具体的有:
利用Snake模型滤除边缘上毛刺产生的伪边缘点;参数化为v(s)=(x(s),y(s))的轮廓的Snake模型表示为下式:
在参数v(s)=(x(s),y(s))中,x(s)和y(s)分别表示每个控制点在图像中的坐标位置,s是以傅立叶变换形式描述边界的自变量;式(1)中,v′(s)为弹性能;v″(s)为弯曲能;Eext为外部能;外部能的表达式为:
Eext(v(s))=-|▽[Gσ(x,y)*I(x,y)]|2 (2)
式(2)中,I(x,y)为二维灰度图像,*表示卷积,Gσ(x,y)为标准差为σ的二维高斯函数,▽为梯度算子;
选择参数α、β,使式(1)中能量E趋于最小,即E→min,则使得轮廓全局最优地逼近真实边缘点;式(1)的连续曲线能量函数写为下面的离散化表示形式:
迭代计算边缘点邻域点的能量值,并将选择能量值最小的邻域点作为真实点,由此获得Snake模型下的最优轮廓;
步骤3.2、以最优轮廓为基础,消除原始边缘上的伪边缘点;分别计算原始边缘点到所有轮廓点的距离平方n表示边缘点个数,i表示轮廓点个数;求解距离平方的最小值设定阈值T,判断是否有Dmin>T;如果是,则删除该边缘点,如果否,则保留该边缘点;
步骤3.3、对步骤3.2获得的边缘的稀疏区域进行最邻域插值;当相邻边缘点间距大于边缘点平均距离的2倍,在两边缘点间利用最邻近插值法进行补点。
其进一步的技术方案为,所述步骤5包括:
步骤5.1、首先将步骤4得到的边缘点转换至坐标系的第一象限:xi=xi-min(xi),yi=yi-min(yi);然后通过下式计算圆心坐标a,b和两个候选半径R和Rpol:
式(12)(13)中有,
式(12)~式(20)中,公式中的x和y表示边缘点的坐标,n表示边缘点的个数;
然后通过式(12)~(18)计算圆心坐标a,b和两个候选半径R和Rpol,并通过式(19)计算两个候选半径下圆拟合的均方根误差RMSER和
步骤5.2、保留RMSE值小的那组参数(a*,b*,R);然后利用迭代的思想降低圆弧拟合的均方差,提高定位精度;具体的:判断是否有RMSEm-1<RMSEm,如果是,则m=m+1;m表示迭代次数,其初始值为2;每次迭代通过改变ρ*的值;其中是经过验证的黄金数,取从而通过式(20)改变tanθ的值,进而通过式(12)~(15)再次计算圆心和两个候选半径,如此循环,直至RMSEm-1≥RMSEm,停止迭代;
步骤5.3、翻转变换得到圆心:a=a*+min(xi),b=b*+min(yi),并输出(a,b,R)。
本发明的有益技术效果是:
本发明针对轮廓附着物对视觉测量中圆弧拟合的干扰以及噪声干扰下短圆弧段拟合结果不稳定的问题,提出一种基于Snake模型和迭代极性转换回归的圆弧拟合方法,首先通过Snake模型对边缘优化,去除伪边缘点,并插补邻近边缘点,然后通过迭代极性转换回归精确定位圆弧的中心,并求解圆的半径。该方法不仅解决了附着物干扰问题,而且提高了短圆弧定位的可靠性和精确度。
具体的,本发明提供了一种精确、可靠的定位受附着物干扰的轴类工件检测图像中的R角中心的方法,并将其应用在轴类工件的尺寸测量上。对框选的多边形图像感兴趣区域进行图像预处理,消除原始图像的孤立噪声点,提高图像的对比度,便于之后的边缘提取;通过Canny算子进行边缘提取,实现边缘的初定位;利用Snake模型滤除边缘毛刺产生的伪边缘点,并在稀疏区域插值边缘点,即对原始边缘进行优化;采用Zernike正交矩对边缘进行精定位,获取精确的亚像素级边缘;最后通过迭代极性转换回归方法对亚像素边缘点进行圆拟合,获得精确的R角的圆心位置和半径。该检测方法解决了由于附着物干扰降低轴类工件R角检测的精确度问题,以及短圆弧拟合的可靠性问题。
附图说明
图1是本发明整体流程图。
图2是Snake模型优化边缘的流程图。
图3是迭代极性转换圆拟合算法流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合具体实施例,并参照附图,对本发明作进一步详细说明。
本发明的目的是针对轮廓附着物对视觉测量中R角检测干扰的问题以及短圆弧拟合不稳定问题,设计一种能够可靠且精确定位感兴趣区域内R角中心的方法,流程分为五个部分:图像预处理、边缘粗定位、边缘优化、亚像素定位,圆弧中心定位,具体流程如图1所示。
具体步骤为:
步骤1、对感兴趣区域进行预处理。预处理的步骤具体为:
步骤1.1、用双边滤波器有效地保边去噪;
步骤1.2、然后用高斯滤波器,对感兴趣区域的边缘进行平滑,改善图像质量;
步骤1.3、最后利用灰度变换提高图像的对比度。
步骤2、边缘粗定位。
利用Ganny算子进行边缘检测,获得一个二值边缘图。
步骤3、边缘优化。边缘优化的具体流程如图2所示。首先手动获取Snake模型的控制点(xi,yi);
步骤3.1、通过对能量函数迭代求解方法得到最佳轮廓;具体的有:
利用Snake模型滤除边缘上毛刺产生的伪边缘点。Snake模型是能够使能量最小化的样条曲线,通过内外力的共同作用将轮廓点收缩到图像,提取图像轮廓。一条参数化为v(s)=(x(s),y(s))的轮廓的Snake模型表示为下式:
在参数v(s)=(x(s),y(s))中,x(s)和y(s)分别表示每个控制点在图像中的坐标位置,s是以傅立叶变换形式描述边界的自变量。v′(s)为弹性能,能够抑制轮廓伸长,v″(s)为弯曲能,能够抑制轮廓弯曲,两者统称为内部能,两处内能的综合作用是使轮廓平滑地收缩。Eext为外部能,外部能的作用是使轮廓收缩到一定的图像特征上,针对不同的特征,具有不同的表达形式,为了使轮廓能够收缩到图像边缘,外部能的表达式可写为:
Eext(v(s))=-|▽[Gσ(x,y)*I(x,y)]|2 (2)
式(2)中,I(x,y)为二维灰度图像,*表示卷积,Gσ(x,y)为标准差为σ的二维高斯函数,▽为梯度算子。
选择合适的参数α、β,可使得E→min,就可以使轮廓全局最优地逼近真实边缘点。式(1)的连续曲线能量函数重写为下面的离散化表示形式:
迭代计算边缘点邻域点的能量值,并将选择能量值最小的邻域点作为真实点,由此获得Snake模型下的最优轮廓。
步骤3.2、以最优轮廓为基础,消除原始边缘上的伪边缘点。分别计算原始边缘点到所有轮廓点的距离平方n表示边缘点个数,i表示轮廓点个数,求解距离平方的最小值设定阈值T为0.5,判断是否有如果是,则删除该边缘点,如果否,则保留该边缘点
步骤3.3、对步骤3.2获得的边缘的稀疏区域进行最邻域插值。当相邻边缘点间距大于边缘点平均距离的2倍,那么在两边缘点间利用最邻近插值法进行补点。
步骤4、亚像素定位。
首先图像中某一点的m阶n次的Zernike矩的定义为:
式(4)中,f(x,y)为离散的二维数字图像,x、y为图像中点的坐标位置,n为正整数或0,m为正负整数,且满足n-|m|为偶数和|m|≤n的条件限制,ρ表示原点到象素(x,y)的向量的距离;θ表示按照逆时针方向向量ρ跟x轴之间的夹角,x2+y2≤1,*表示共轭,复数多项式Vnm(ρ,θ)表示为:
当图像旋转角度后,Zernike矩和旋转前关系:
利用此方法进行边缘定位时需要用到Z00,Z11,Z20三个Zernike矩,根据旋转前后的关系可以得出以下三个关系:
Z′00=Z00 (7)
Z′20=Z20 (9)
当图像旋转角度后,图像关于x轴上下对称,此时的Z11虚部为零,即
采用Zernike矩对边缘点进行细分。算法的具体流程如下:
根据模板大小为N*N,计算每一点的Zernike矩Z11和Z20,Im[Z11]为Z11的虚部,Re[Z11]为Z11的实部。
计算旋转角度
计算
计算l=Z20/Z′11;
根据亚像素坐标(x′,y′)与像素坐标(x,y)关系:
计算亚像素坐标。
步骤5、圆弧中心定位。输入n个坐标点(xi,yi);
步骤5.1、利用极性转换回归定位圆弧中心。首先将步骤4得到的边缘点转换至坐标系的第一象限,xi=xi-min(xi),yi=yi-min(yi);
式(18)中有,
式(12)~式(20)中,公式中的x和y表示边缘点的坐标,n表示边缘点的个数。
然后通过式(12)~(18)计算圆心坐标a,b和两个候选半径R和Rpol,并通过式(19)计算两个候选半径下圆拟合的均方根误差RMSER和
步骤5.2、保留RMSE值小的那组参数(a*,b*,R);然后利用迭代的思想降低圆弧拟合的均方差,提高定位精度。主要思想是将数据点转换为到原点的距离,然后搜索一个RMSE减小量一直迭代,直至当前RMSE值小于前一次RMSE值,停止迭代。
具体的,判断是否有RMSEm-1<RMSEm,如果是,则m=m+1,m表示迭代次数,随着迭代次数的增加,会在每次迭代中通过改变ρ*的值,从而每次都刷新圆心坐标,因此每次迭代会改变RMSE值。在中参数是经过验证的黄金数,取从而通过式(20)改变tanθ的值。进一步的,通过式(12)~(15)再次计算圆心和两个候选半径,如此循环,直至RMSEm-1≥RMSEm,停止迭代。
步骤5.3、翻转变换得到圆心:a=a*+min(xi),b=b*+min(yi),并输出(a,b,R)。
以上所述的仅是本发明的优选实施方式,本发明不限于以上实施例。可以理解,本领域技术人员在不脱离本发明的精神和构思的前提下直接导出或联想到的其他改进和变化,均应认为包含在本发明的保护范围之内。
Claims (3)
1.一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法,其特征在于,包含以下步骤:
步骤1、对图像进行预处理;
步骤2、利用Canny算子对该预处理区域进行边缘提取;
步骤3、利用Snake模型对所述区域的边缘优化,去除伪边缘点,并插补邻近边缘点;
步骤4、利用Zernike正交矩获取亚像素级边缘;
步骤5、通过迭代极性转换回归实现R角中心的精确定位。
2.如权利要求1所述的一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法,其特征在于,所述步骤3包括:
步骤3.1、通过对能量函数迭代,得到最优轮廓;具体的有:
利用Snake模型滤除边缘上毛刺产生的伪边缘点;参数化为v(s)=(x(s),y(s))的轮廓的Snake模型表示为下式:
在参数v(s)=(x(s),y(s))中,x(s)和y(s)分别表示每个控制点在图像中的坐标位置,s是以傅立叶变换形式描述边界的自变量;式(1)中,v′(s)为弹性能;v″(s)为弯曲能;Eext为外部能;外部能的表达式为:
式(2)中,I(x,y)为二维灰度图像,*表示卷积,Gσ(x,y)为标准差为σ的二维高斯函数,为梯度算子;
选择参数α、β,使式(1)中能量E趋于最小,即E→min,则使得轮廓全局最优地逼近真实边缘点;式(1)的连续曲线能量函数写为下面的离散化表示形式:
迭代计算边缘点邻域点的能量值,并将选择能量值最小的邻域点作为真实点,由此获得Snake模型下的最优轮廓;
步骤3.2、以最优轮廓为基础,消除原始边缘上的伪边缘点;分别计算原始边缘点到所有轮廓点的距离平方n表示边缘点个数,i表示轮廓点个数;求解距离平方的最小值设定阈值T,判断是否有Dmin>T;如果是,则删除该边缘点,如果否,则保留该边缘点;
步骤3.3、对步骤3.2获得的边缘的稀疏区域进行最邻域插值;当相邻边缘点间距大于边缘点平均距离的2倍,在两边缘点间利用最邻近插值法进行补点。
3.如权利要求1所述的一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法,其特征在于,所述步骤5包括:
步骤5.1、首先将步骤4得到的边缘点转换至坐标系的第一象限:xi=xi-min(xi),yi=yi-min(yi);然后通过下式计算圆心坐标a,b和两个候选半径R和Rpol:
式(12)(13)中有,
式(12)~式(20)中,公式中的x和y表示边缘点的坐标,n表示边缘点的个数;
然后通过式(12)~(18)计算圆心坐标a,b和两个候选半径R和Rpol,并通过式(19)计算两个候选半径下圆拟合的均方根误差RMSER和
步骤5.2、保留RMSE值小的那组参数(a*,b*,R);然后利用迭代的思想降低圆弧拟合的均方差,提高定位精度;具体的:判断是否有RMSEm-1<RMSEm,如果是,则m=m+1;m表示迭代次数,其初始值为2;每次迭代通过改变ρ*的值;其中是经过验证的黄金数,取从而通过式(20)改变tanθ的值,进而通过式(12)~(15)再次计算圆心和两个候选半径,如此循环,直至RMSEm-1≥RMSEm,停止迭代;
步骤5.3、翻转变换得到圆心:a=a*+min(xi),b=b*+min(yi),并输出(a,b,R)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810354070.5A CN108830868B (zh) | 2018-04-19 | 2018-04-19 | 一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810354070.5A CN108830868B (zh) | 2018-04-19 | 2018-04-19 | 一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108830868A true CN108830868A (zh) | 2018-11-16 |
CN108830868B CN108830868B (zh) | 2021-07-16 |
Family
ID=64154848
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810354070.5A Active CN108830868B (zh) | 2018-04-19 | 2018-04-19 | 一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108830868B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111724366A (zh) * | 2020-06-16 | 2020-09-29 | 武汉工程大学 | 一种激光腔识别方法及装置 |
CN112381850A (zh) * | 2020-12-04 | 2021-02-19 | 亿嘉和科技股份有限公司 | 柜体表面圆形目标自动分割方法、系统、装置及存储介质 |
CN112990207A (zh) * | 2019-12-02 | 2021-06-18 | 中移智行网络科技有限公司 | 车牌图像提取方法、装置、系统、存储介质和计算机设备 |
CN113436212A (zh) * | 2021-06-22 | 2021-09-24 | 广西电网有限责任公司南宁供电局 | 一种断路器静触头齿合状态图像检测内轮廓的提取方法 |
CN114004981A (zh) * | 2021-11-10 | 2022-02-01 | 柳州沪信汽车科技有限公司 | 一种非完整点云条件下的车身r角视觉检测方法及系统 |
CN115063431A (zh) * | 2022-08-19 | 2022-09-16 | 山东远盾网络技术股份有限公司 | 一种基于图像处理的汽车零部件质量溯源方法 |
CN116823925A (zh) * | 2023-08-30 | 2023-09-29 | 苏州聚视兴华智能装备有限公司 | 高精度o型橡胶圈内径外径测量方法、装置及电子设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040186369A1 (en) * | 2003-03-21 | 2004-09-23 | Scimed Life Systems, Inc. | Scan conversion of medical imaging data |
CN104070523A (zh) * | 2013-03-27 | 2014-10-01 | 深圳市生命之泉科技发展有限公司 | 基于空间坐标转换的工业机器人实时圆弧插补实现方法 |
CN104570928A (zh) * | 2013-10-29 | 2015-04-29 | 中国科学院沈阳自动化研究所 | 基于共形参数化的网格曲面上数控加工轨迹规划方法 |
CN107341802A (zh) * | 2017-07-19 | 2017-11-10 | 无锡信捷电气股份有限公司 | 一种基于曲率与灰度复合的角点亚像素定位方法 |
CN107830792A (zh) * | 2016-09-16 | 2018-03-23 | Nm数字模拟有限公司 | 用于确定位置测量系统的位置指示器的位置的方法 |
-
2018
- 2018-04-19 CN CN201810354070.5A patent/CN108830868B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040186369A1 (en) * | 2003-03-21 | 2004-09-23 | Scimed Life Systems, Inc. | Scan conversion of medical imaging data |
CN104070523A (zh) * | 2013-03-27 | 2014-10-01 | 深圳市生命之泉科技发展有限公司 | 基于空间坐标转换的工业机器人实时圆弧插补实现方法 |
CN104570928A (zh) * | 2013-10-29 | 2015-04-29 | 中国科学院沈阳自动化研究所 | 基于共形参数化的网格曲面上数控加工轨迹规划方法 |
CN107830792A (zh) * | 2016-09-16 | 2018-03-23 | Nm数字模拟有限公司 | 用于确定位置测量系统的位置指示器的位置的方法 |
CN107341802A (zh) * | 2017-07-19 | 2017-11-10 | 无锡信捷电气股份有限公司 | 一种基于曲率与灰度复合的角点亚像素定位方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112990207A (zh) * | 2019-12-02 | 2021-06-18 | 中移智行网络科技有限公司 | 车牌图像提取方法、装置、系统、存储介质和计算机设备 |
CN111724366A (zh) * | 2020-06-16 | 2020-09-29 | 武汉工程大学 | 一种激光腔识别方法及装置 |
CN111724366B (zh) * | 2020-06-16 | 2023-08-25 | 武汉工程大学 | 一种激光腔识别方法及装置 |
CN112381850A (zh) * | 2020-12-04 | 2021-02-19 | 亿嘉和科技股份有限公司 | 柜体表面圆形目标自动分割方法、系统、装置及存储介质 |
CN113436212A (zh) * | 2021-06-22 | 2021-09-24 | 广西电网有限责任公司南宁供电局 | 一种断路器静触头齿合状态图像检测内轮廓的提取方法 |
CN114004981A (zh) * | 2021-11-10 | 2022-02-01 | 柳州沪信汽车科技有限公司 | 一种非完整点云条件下的车身r角视觉检测方法及系统 |
CN115063431A (zh) * | 2022-08-19 | 2022-09-16 | 山东远盾网络技术股份有限公司 | 一种基于图像处理的汽车零部件质量溯源方法 |
CN115063431B (zh) * | 2022-08-19 | 2022-11-11 | 山东远盾网络技术股份有限公司 | 一种基于图像处理的汽车零部件质量溯源方法 |
CN116823925A (zh) * | 2023-08-30 | 2023-09-29 | 苏州聚视兴华智能装备有限公司 | 高精度o型橡胶圈内径外径测量方法、装置及电子设备 |
CN116823925B (zh) * | 2023-08-30 | 2023-11-17 | 苏州聚视兴华智能装备有限公司 | 高精度o型橡胶圈内径外径测量方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN108830868B (zh) | 2021-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108830868A (zh) | 一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法 | |
CN107341802B (zh) | 一种基于曲率与灰度复合的角点亚像素定位方法 | |
CN106204610A (zh) | 基于图像处理的轴孔装配自动对中系统 | |
CN108986126A (zh) | 基于高斯拟合亚像素边缘检测及改进ransac算法的圆心检测法 | |
CN103425988A (zh) | 一种具有圆弧几何基元的实时定位与匹配方法 | |
CN111047588A (zh) | 一种轴型小零件尺寸的图像化测量方法 | |
CN108416787A (zh) | 应用于机器视觉检测的工件直线边缘定位方法 | |
Yang et al. | Polygon-invariant generalized Hough transform for high-speed vision-based positioning | |
CN111161260A (zh) | 一种基于深度学习的热轧带钢表面缺陷检测方法及装置 | |
CN113077392B (zh) | 一种针对指针表模糊照片的高准确度自动读数方法 | |
CN111311618A (zh) | 一种基于高精度几何基元提取的圆弧工件匹配与定位方法 | |
CN104715491A (zh) | 一种基于一维灰度矩的亚像素边缘检测方法 | |
Li et al. | Adaptive bilateral smoothing for a point-sampled blade surface | |
CN111027530A (zh) | 基于轮胎压印字符识别的预处理方法 | |
CN101840513A (zh) | 一种图像形状特征提取方法 | |
Jiang et al. | A novel centerline extraction algorithm for a laser stripe applied for turbine blade inspection | |
CN103839274A (zh) | 一种基于几何比例关系的扩展目标跟踪方法 | |
Teutsch et al. | Real-time detection of elliptic shapes for automated object recognition and object tracking | |
Liu et al. | A comparative study of different corner detection methods | |
CN112529869B (zh) | 一种阀套节流方孔检测方法 | |
CN106780562A (zh) | 一种基于像差图图像处理的主轴回转误差新算法 | |
Rao et al. | A real-time auto-recognition method for pointer-meter under uneven illumination | |
Zhang et al. | Non-contact high-precision aperture detection based on distortion compensation and adaptive mean blur algorithm | |
Zheng et al. | Pointer instrument recognition algorithm based on Haar-like feature and polar expansion | |
Zou et al. | Contour Primitive Recognition of Planar Thin Sheet Metal Based on Feature Point |
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 |