CN109993787A - 一种图像法螺纹中径求解方法 - Google Patents

一种图像法螺纹中径求解方法 Download PDF

Info

Publication number
CN109993787A
CN109993787A CN201910126315.3A CN201910126315A CN109993787A CN 109993787 A CN109993787 A CN 109993787A CN 201910126315 A CN201910126315 A CN 201910126315A CN 109993787 A CN109993787 A CN 109993787A
Authority
CN
China
Prior art keywords
point
row
pixel
image
thread
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
CN201910126315.3A
Other languages
English (en)
Other versions
CN109993787B (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.)
Shaanxi University of Technology
Original Assignee
Shaanxi 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 Shaanxi University of Technology filed Critical Shaanxi University of Technology
Priority to CN201910126315.3A priority Critical patent/CN109993787B/zh
Publication of CN109993787A publication Critical patent/CN109993787A/zh
Application granted granted Critical
Publication of CN109993787B publication Critical patent/CN109993787B/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
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30164Workpiece; Machine component

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种图像法螺纹中径求解方法,具体包括如下步骤:步骤1,获取螺纹图像,并建立图像坐标系xoy;步骤2,逐行遍历式扫描螺纹图像,寻找螺纹各牙侧边中径点的位置;步骤3,根据步骤2找到的螺纹中径点计算螺纹图像选择区域总像素A;步骤4,根据步骤3所得的螺纹图像选择区域总像素A求螺纹中径d2。本发明解决了现有螺纹测量方法中难以确定螺纹中径的问题。

Description

一种图像法螺纹中径求解方法
技术领域
本发明属于螺纹测量技术领域,涉及一种图像法螺纹中径求解方法。
背景技术
机器视觉等图像法测量螺纹具有效率高的特点,其可通过一次性图像采集获得多项螺纹参数。随着图像传感器分辨率的不断提高,和曝光时间的不断缩小,图像法螺纹测量开始越来越多地应用于生产实践。然而,由于受螺纹制造过程中牙形误差,和图像法螺纹测量时图像坐标系与待测螺纹轴线很难完全重合这一客观条件的影响,使得利用图像法求解螺纹中径变得十分困难。对图像法螺纹测量而言,找到一种实效的螺纹中径求解方法就变得非常迫切。
发明内容
本发明的目的是提供一种图像法螺纹中径求解方法,解决了现有螺纹测量方法中难以确定螺纹中径的问题。
本发明所采用的技术方案是,一种图像法螺纹中径求解方法,具体包括如下步骤:
步骤1,在保证螺纹轴线oo′与螺纹图像坐标轴x最小夹角θ不超过条件下,获取螺纹图像,对螺纹图像进行二值化后,保证黑色像素总数超过图像总像素80%为条件裁剪螺纹图像,并建立图像坐标系xoy;
步骤2,逐行遍历式扫描螺纹图像,寻找螺纹各牙侧边中径点的位置;
步骤3,根据步骤2找到的螺纹中径点计算螺纹图像选择区域总像素A;
步骤4,根据步骤3所得的螺纹图像选择区域总像素A求螺纹中径d2
本发明的特点还在于,
步骤2的具体过程如下:
按行从第0行至M/2行开始扫描图像,获取像素灰度变化最大的行的行号im和灰度变化数m;继续按行扫描图像,获取另外一个灰度变化数为m的最大行号,记为jm;在经历扫描行数大于图像总行数一半M/2后,仍然按行扫描图像,获取像素灰度变化最大的行的行号i′m和灰度变化数m′,再继续按行扫描图像,获取另外一个灰度变化数为m′的行号最大行号j′m;然后取
并当m≠m′时,取二者较小者为灰度变化数m″,即
m″=min(m,m′) (26);
当m″为偶数时,取螺纹牙数
N0=m″/2 (27);
否则取
以n=2N0为上限,从第0列至第n-1列对第i行检索出左侧像素灰度为255、右侧为0或左侧像素灰度为0、右侧为255的点A1、A2、A3、……、An,n=1,2,3,……,2N0,以列、行号为坐标值存入二维数组A[xk,yk],对第j行检索出左侧像素灰度为255、右侧为0或左侧像素灰度为0、右侧为255的点E1、E2、E3、……、En以列、行号为坐标值存入二维数组E[xk,yk];对第i′,j′检索出左侧像素灰度为255、右侧为0或左侧像素灰度为0、右侧为255的点A′1、A′2、A′3、……、A′n和E′1、E′2、E′3、……、E′n,依次将这些点的列、行号存入二维数组A′[xk,yk]和E′[xk,yk]中(其中xk为第k个对应点的列号,yk为第k个对应点的行号,k=1,2,3,……,n),用A[xk]表示第i行上点Ak列坐标,A′[xk]表示第i′行上点A′k列坐标,E[xk]表示第j行上点Ek列坐标,E′[xk]表示第j′行上点E′k列坐标;
取p=(1,2,3,……,n-1),令ap为第i行Ap+1点与Ap点列像素之差,即
ap=A[xp+1]-A[xp] (29);
bp为第j行Ep+1点与Ep点列像素之差,即
bp=E[xp+1]-E[xp] (30);
分别形成n-1个数构成的数列a1、a2、a3、…、ap、…、an-1和b1、b2、b3、…、bp、…、bn-1
R=N0-1 (31);
也即
并令q=1,2,3,……,R,对于第i,j上的点A2q-1、A2q、A2q+1、E2q-1、E2q、E2q+1以q为下标序列,代入求得的ap、bp数列的p下标,用式
求得y1、y2、y3、……、yR,以y1、y2、y3、……、yR为行号按其左侧像素为黑色,右侧像素为白色,检索其对应列坐标,并依次得到列坐标值xCq,q=1,2,3,……、R,令yCq=yq,进而形成中径上点C1、C2、C3、……、CR的列、行坐标值(xCq,yCq);
对于第i,j上的点A2q、A2q+1、A2q+2、E2q、E2q+1、E2q+2以q为下标序列,代入求得的ap、bp数列的p下标,用式
求得y′1、y′2、y′3、……、y′R,以y′1、y′2、y′3、……、y′R为行号按其右侧像素为黑色,左侧像素为白色,检索其对应列坐标值 进而形成中径上点C’1、C’2、C’3、……、C’R的列、行坐标值
对于第i′,j′上的点A′1、A′2、A′3、……、A′n和E′1、E′2、E′3、……、E′n,及形成的二维数组A′[xk,yk]和E′[xk,yk],取p=1,2,3,……,n-1,令a′p为第i′行A′p+1点与A′p点列像素之差,即
a′p=A′[xp+1]-A′[xp] (35);
b′p为第j′行E′p+1点与E′p点列像素之差,即
b′p=E′[xp+1]-E′[xp] (36);
分别形成n-1个数构成的数列a′1、a′2、a′3、…、a′p、…、a′n-1和b′1、b′2、b′3、…、b′p、…、b′n-1
R=N0-1 (37);
并令q=1,2,3,……,R,对于第i′,j′上的点A′2q-1、A′2q、A′2q+1、E′2q-1、E′2q、E′2q+1以q为下标序列,代入求得的a′p、b′p数列的p下标,用式
以这些行号按其左侧像素为黑色,而右侧像素为白色检索其对应列坐标进而形成下半部分螺纹中径上点F1、F2、F3、……、FR的列、行坐标值
对于第i′,j′上的点A′2q、A′2q+1、A′2q+2、E′2q、E′2q+1、E′2q+2以q为下标代替a′p、b′p数列的p下标,用式
求得y″′1、y″′2、y″′3、……、y″′R,以这些行号按其左侧像素为白色,右侧像素为黑色,检索其对应列坐标进而,形成下半部分螺纹中径上点F’1、F’2、F3、……、F’R的列、行坐标值
步骤3的具体过程如下:
对步骤2中形成的上半部分中径点C1、C2、C3、……、CR及C’1、C’2、C’3、……、C’R,分别以(xCq,yCq)和(xC′q,yC′q)表示其坐标值;对步骤3中形成的下半部分中径点F1、F2、F3、……、FR和F’1、F’2、F’3、……、F’R,则以(xFq,yFq)和(xF′q,yF′q)表示其坐标值,下标q=1,2,3……,R;对用上、下相邻两点Cq、C’q、Fq、F’q点构成的四边形,用列、行坐标按高斯四边形面积公式
和下列公式(40):
求得S1,S2,S3,……,SR之和SS1
对用上、下相邻两点C’q、Cq+1、F’q、Fq+1点构成的四边形,用列、行坐标按高斯四边形面积公式:
和下列公式(41):
求得S’1,S’2,S’3,……,S’R之和SS2
A=SS1+SS2=S1+S2=S1+S2+……+SR+S′1+S′2+……+S′R (20)。
步骤4的具体过程如下:
以C1和C’R的坐标求得由点C1至C’R距离L0,则可求得螺纹中径d2
式中,d2为螺纹中径,单位:mm;
A为螺纹图像选择区域总像素,单位:平方像素;
L0为螺纹图像选择区域单边像素数,单位:像素;
Ky为y坐标方向上的像素标定系数,单位:mm/像素。
本发明的有益效果如下:
(1)这种方法符合精密测量技术中的平均读数原理,可以自然消除毛刺、微小油珠或水珠以及个别牙廓边缘缺陷的影响;
(2)直接采用图像像素做为计算依据,减少了如图像旋转等中间环节,计算过程简单;
(3)采用牙形中部作为计算依据,规避了牙顶毛刺和牙根油污滞留对中径计算的影响;
(4)降低了图像法螺纹测量时被测螺纹工件安装要求,简化了图像处理过程,提高了图像法螺纹测量的适应性。
附图说明
图1(a)、(b)是螺纹轴线与图像坐标系之间的关系示意图;
图2(a)~(e)是本发明一种图像法螺纹中径求解方法中螺纹中径点分析过程示意图;
图3是本发明一种图像法螺纹中径求解方法中螺纹中径求解过程示意图;
图4是本发明一种图像法螺纹中径求解方法中建立的螺纹图像坐标示意图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
对图像方法进行螺纹测量的采集系统而言,按采集的螺纹图像与图像坐标系关系划分,可以分为两种情况。第一种情况见图1(a)所示,采集的图像中螺纹轴线O’O’与图像坐标系某坐标轴平行,其经二值化后,螺纹图像可以视为由M行N列的M×N的黑、白像素阵列。按螺纹中径定义:螺纹中径是指母线通过螺纹牙型上沟槽和凸起宽度相等的假想圆柱的直径。在这种情况下,螺纹中径求解可按相邻牙侧边线上点像素数目相等为条件,通过计算机程序遍历图像中所有0行至M-1行,寻找相邻边线黑像素与白像素相等的行。比如图1(a)中,若在图像上部分存在第i行,相邻牙侧黑色像素数等于紧挨着的相邻白色像素数,即P1=P2,则可以确定第i行是一条中径边,同样,在图像下半部分也可以同样条件遍历得到第j行是另一条中径边,此时可由下列关系式得到中径d2计算式:
d2=(j-i)Ky (1);
式中,d2为螺纹中径;
i,j为中径边所在行的行号;
Ky为y坐标方向上的像素标定系数;
第二中情况,当图像中的螺纹轴线与图像坐标系坐标轴不平行时,如图1(b)所示,就不能使用上述求解方法了。螺纹工作面为螺纹牙两个侧面,因而螺纹在实际制造和使用过程中,对牙顶和牙根没有提出较高要求,毛刺等缺陷较多,加上油污的影响,牙顶和牙根就不能作为其自身定位基准,而牙侧边则由于螺纹制造过程中产生的牙型角误差,使得中径的求解无法得到相对统一的结论。
本发明一种图像法螺纹中径求解方法,该方法可就第一种和第二种情况,给出统一解决方案。计算原理如下:
在螺纹轴线与其图像坐标轴最小夹角不超过条件下(一般定位装置均能达到),对螺纹图像按照图2(a)所示自上而下逐行扫描,得到螺纹上部轮廓扫描行i、j和下部轮廓扫描行i’、j’,i、j分别为图像中上部螺纹靠近牙顶和牙根部分的行号,i’、j’为分别为下部螺纹靠近牙根和牙顶部分的行号;当i、j’足够靠近牙顶,而j、i’足够靠近牙根时,展开几何分析如图2(b)所示;i,j是与图像坐标轴平行的靠近牙顶和牙根部任意两行行号,第i行和第j行还分别在牙凸形成总像素数为S的图2(b)所示阴影面积,在牙槽部形成总像素数为S’的图2(b)所示阴影面积;若i行上被牙凸截得像素数为a1,即图2(b)中A1A2,被其相邻牙槽截得像素数为a2,即图2(b)中A2A3;而靠近牙根的j行上被牙凸截得像素数为b1,即图2(b)中E1E2,被其相邻牙槽截得像素数为b2,即图2(b)中E2E3;若C1点为螺纹中径与螺纹牙侧边A2E2交点,TT为过C1点的中径线,过C1点做平行于牙凸侧边A2E2平行线BD,由图中几何关系可以得到
和四边形A1A2E1E2面积与四边形A2A3E2E3面积之差等于三角形ΔA2BC1面积与三角形ΔC1DE2面积之差,即
进而得到:
y1为图像坐标系中螺纹中径线TT与牙侧边A2E2交点C1的行号。
如图2(c)所示,对相邻牙侧A3E3中径点C’1的行号y’1,也可由四边形A2A3E2E3面积(即图中S’)与四边形A3A4E3E4面积(即图中S″)等于三角形ΔA3B’C’1面积与三角形ΔC’1D’E3面积之差,即
进而得到:
采用上述相同的方法得到,即
同理,对图2(a)中螺纹图像的下部轮廓进行同样的牙廓中径点行号计算,如图2(d)、(e)所示;在图2(d)中,若靠近牙根的扫描行i’与牙侧边交点分别为A’1、A’2,靠近牙顶的扫描行j’与牙侧边交点分别为E’1、E’2,若设i’行上被牙槽截得像素数为a’1,即图2中A’1A’2,被其相邻牙凸截得像素数为a’2,而靠近牙顶的j’行上被牙凸截得像素数为b’1,即图2(d)中E’1E’2,被其相邻牙槽截得像素数为b’2,对牙侧边A’2E’2中径点F1的行号y”1,可由四边形A’1A’2E’1E’2面积(即图中S0)与四边形A’2A’3E’2E’3面积(即图中S’0)等于三角形ΔB”A’2F1面积与三角形ΔF1D”E’2面积之差,即
进而得到:
同样由四边形A’2A’3E’2E’3面积(即图中S’0)与四边形A’3A’4E’3E’4面积(即图中S”0)等于三角形ΔB”’A’3F’1面积与三角形ΔF’1D”’E’3面积之差,
进而得到:
如图3所示,依照这种方法,用牙凸1与其右侧牙槽求得牙凸1右侧边中径点C1,用牙凸2与其左侧牙槽求得牙凸2左侧边中径点C’1,用牙凸2与其右侧牙槽求得牙凸2右侧边中径点C2,用牙凸3与其左侧牙槽求得牙凸3左侧边中径点C’2,用牙凸3与其右侧牙槽求得牙凸3右侧边中径点C3,等等,至牙凸R,依次求得中径点C1、C2、C3、……、CR和C’1、C’2、C’3、……、C’R的行、列坐标;以同样方法求得螺纹另一投影牙侧中径点F1、F2、F3、……、FR和F’1、F’2、F’3、……、F’R的列、行坐标;当用以(xCq,yCq)、(xC′q,yC′q)、(xFq,yFq)和(xF′q,yF′q)分别表示Cq、C’q、Fq、F’q各点列、行坐标,按照高斯多边形面积公式
在q=1,2,3,……,R时,依次求得S1,S2,S3,……,SR
SS1=S1+S2+S3,+……+SR (17);
同理,对C’q、Cq+1、F’q、Fq+1点构成的四边形,用列、行坐标按高斯四边形面积公式
在q=1,2,3,……,R时,求得S’1,S’2,S’3,……,S’R
SS2=S’1+S’2+S’3+……+S’R; (19);
A=SS1+SS2=S1+S2=S1+S2+……+SR+S′1+S′2+……+S′R (20);
以C1和C’R的坐标求得由点C1至C’R距离L0(也即C1C’R连线段的总像素数),则可求得中径
式中,Ky为y坐标方向上的像素标定系数。
本发明一种图像法螺纹中径求解方法,具体包括如下步骤:
步骤1,根据安装条件,保证螺纹轴线oo′与其图像坐标轴x最小夹角θ不超过条件下获取螺纹图像,对螺纹图像进行二值化后,以保证黑色像素总数超过图像总像素80%为条件裁剪螺纹图像,并建立图像坐标系xoy,如图4所示;
步骤2,按行从第0行至M/2行开始扫描图像,获取像素灰度变化最大的行的行号im和灰度变化数m;继续按行扫描图像,获取另外一个灰度变化数为m的最大行号,记为jm;在经历扫描行数大于图像总行数一半M/2后,仍然按行扫描图像,获取像素灰度变化最大的行的行号i′m和灰度变化数m′,再继续按行扫描图像,获取另外一个灰度变化数为m′的行号最大行号j′m;然后取
并当m≠m′时,取二者较小者为灰度变化数m″,即
m″=min(m,m′)(26);
当m″为偶数时,取螺纹牙数
N0=m″/2(27);
否则取
以n=2N0为上限,从第0列至第n-1列对第i行检索出左侧像素灰度为255、右侧为0或左侧像素灰度为0、右侧为255的点A1、A2、A3、……、An(n=1,2,3,……,2N0),以列、行号为坐标值存入二维数组A[xk,yk],对第j行检索出左侧像素灰度为255、右侧为0或左侧像素灰度为0、右侧为255的点E1、E2、E3、……、En以列、行号为坐标值存入二维数组E[xk,yk];同样方法,对第i′,j′检索出左侧像素灰度为255、右侧为0或左侧像素灰度为0、右侧为255的点A′1、A′2、A′3、……、A′n和E′1、E′2、E′3、……、E′n,依次将这些点的列、行号存入二维数组A′[xk,yk]和E′[xk,yk]中(其中xk为第k个对应点的列号,yk为第k个对应点的行号,k=1,2,3,……,n),用A[xk]表示第i行上点Ak列坐标,A′[xk]表示第i′行上点A′k列坐标,E[xk]表示第j行上点Ek列坐标,E′[xk]表示第j′行上点E′k列坐标。
取p=(1,2,3,……,n-1),令ap为第i行Ap+1点与Ap点列像素之差,即
ap=A[xp+1]-A[xp] (29);
bp为第j行Ep+1点与Ep点列像素之差,即
bp=E[xp+1]-E[xp] (30);
分别形成n-1个数构成的数列a1、a2、a3、…、ap、…、an-1和b1、b2、b3、…、bp、…、bn-1
R=N0-1 (31);
也即
并令q=1,2,3,……,R,对于第i,j上的点A2q-1、A2q、A2q+1、E2q-1、E2q、E2q+1以q为下标序列,代入求得的ap、bp数列的p下标,用式
求得y1、y2、y3、……、yR,以y1、y2、y3、……、yR为行号按其左侧像素为黑色(比如灰度值为0)右侧像素为白色(比如灰度值为255)检索其对应列坐标xCq(下标q=1,2,3,……、R),为统一起见,令yCq=yq,进而形成中径上点C1、C2、C3、……、CR的列、行坐标值(xCq,yCq);
对于第i,j上的点A2q、A2q+1、A2q+2、E2q、E2q+1、E2q+2以q为下标序列,代入求得的ap、bp数列的p下标,用式
求得y′1、y′2、y′3、……、y′R,以y′1、y′2、y′3、……、y′R为行号按其左侧像素为黑色(比如灰度值为0)右侧像素为白色(比如灰度值为255)检索其对应列坐标值(下标q=1,2,3,……、R),为统一起见,令进而形成中径上点C’1、C’2、C’3、……、C’R的列、行坐标值
同样方法,对于第i′,j′上的点A′1、A′2、A′3、……、A′n和E′1、E′2、E′3、……、E′n,及形成的二维数组A′[xk,yk]和E′[xk,yk],取p=1,2,3,……,n-1,令a′p为第i′行A′p+1点与A′p点列像素之差,即
a′p=A′[xp+1]-A′[xp] (35);
b′p为第j′行E′p+1点与E′p点列像素之差,即
b′p=E′[xp+1]-E′[xp] (36);
分别形成n-1个数构成的数列a′1、a′2、a′3、…、a′p、…、a′n-1和b′1、b′2、b′3、…、b′p、…、b′n-1
R=N0-1 (37);
并令q=1,2,3,……,R,对于第i′,j′上的点A′2q-1、A′2q、A′2q+1、E′2q-1、E′2q、E′2q+1以q为下标序列,代入求得的a′p、b′p数列的p下标,用式
以这些行号按其左侧像素为黑色,而右侧像素为白色检索其对应列坐标进而形成下半部分螺纹中径上点F1、F2、F3、……、FR的列、行坐标值
对于第i′,j′上的点A′2q、A′2q+1、A′2q+2、E′2q、E′2q+1、E′2q+2以q为下标代替a′p、b′p数列的p下标,用式
求得y″′1、y″′2、y″′3、……、y″′R,以这些行号按其左侧像素为白色,右侧像素为黑色,检索其对应列坐标进而,形成下半部分螺纹中径上点F’1、F’2、F3、……、F’R的列、行坐标值
步骤3:如图3所示,对步骤3中形成的上半部分中径点C1、C2、C3、……、CR及C’1、C’2、C’3、……、C’R,分别以(xCq,yCq)和(xC′q,yC′q)表示其坐标值;对步骤3中形成的下半部分中径点F1、F2、F3、……、FR和F’1、F’2、F’3、……、F’R,则以(xFq,yFq)和(xF′q,yF′q)表示其坐标值,下标q=1,2,3……,R;对用上、下相邻两点Cq、C’q、Fq、F’q点构成的四边形,用列、行坐标按高斯四边形面积公式
和下列公式(31):
求得S1,S2,S3,……,SR之和SS1
对用上、下相邻两点C’q、Cq+1、F’q、Fq+1点构成的四边形,用列、行坐标按高斯四边形面积公式
和下列公式(41):
求得S’1,S’2,S’3,……,S’R之和SS2
A=SS1+SS2=S1+S2=S1+S2+……+SR+S′1+S′2+……+S′R (20);
步骤4,
以C1和C’R的坐标求得由点C1至C’R距离L0(也即C1C’R连线段的总像素数),则可求得螺纹中径d2
式中,d2为螺纹中径,单位:mm;
A为螺纹图像选择区域总像素,单位:平方像素;
L0为螺纹图像选择区域单边像素数,单位:像素;
Ky为y坐标方向上的像素标定系数,单位:mm/像素。

Claims (4)

1.一种图像法螺纹中径求解方法,其特征在于:具体包括如下步骤:
步骤1,在保证螺纹轴线oo′与螺纹图像坐标轴x最小夹角θ不超过条件下,获取螺纹图像,对螺纹图像进行二值化后,保证黑色像素总数超过图像总像素80%为条件裁剪螺纹图像,并建立图像坐标系xoy;
步骤2,逐行遍历式扫描螺纹图像,寻找螺纹各牙侧边中径点的位置;
步骤3,根据步骤2找到的螺纹中径点计算螺纹图像选择区域总像素A;
步骤4,根据步骤3所得的螺纹图像选择区域总像素A求螺纹中径d2
2.根据权利要求1所述的一种图像法螺纹中径求解方法,其特征在于:所述步骤2的具体过程如下:
按行从第0行至M/2行开始扫描图像,获取像素灰度变化最大的行的行号im和灰度变化数m;继续按行扫描图像,获取另外一个灰度变化数为m的最大行号,记为jm;在经历扫描行数大于图像总行数一半M/2后,仍然按行扫描图像,获取像素灰度变化最大的行的行号i′m和灰度变化数m′,再继续按行扫描图像,获取另外一个灰度变化数为m′的行号最大行号j′m;然后取
并当m≠m′时,取二者较小者为灰度变化数m″,即
m″=min(m,m′) (26);
当m″为偶数时,取螺纹牙数
N0=m″/2 (27);
否则取
以n=2N0为上限,从第0列至第n-1列对第i行检索出左侧像素灰度为255、右侧为0或左侧像素灰度为0、右侧为255的点A1、A2、A3、……、An,n=1,2,3,……,2N0,以列、行号为坐标值存入二维数组A[xk,yk],对第j行检索出左侧像素灰度为255、右侧为0或左侧像素灰度为0、右侧为255的点E1、E2、E3、……、En以列、行号为坐标值存入二维数组E[xk,yk];对第i′,j′检索出左侧像素灰度为255、右侧为0或左侧像素灰度为0、右侧为255的点A′1、A′2、A′3、……、A′n和E′1、E′2、E′3、……、E′n,依次将这些点的列、行号存入二维数组A′[xk,yk]和E′[xk,yk]中(其中xk为第k个对应点的列号,yk为第k个对应点的行号,k=1,2,3,……,n),用A[xk]表示第i行上点Ak列坐标,A′[xk]表示第i′行上点A′k列坐标,E[xk]表示第j行上点Ek列坐标,E′[xk]表示第j′行上点E′k列坐标;
取p=(1,2,3,……,n-1),令ap为第i行Ap+1点与Ap点列像素之差,即
ap=A[xp+1]-A[xp] (29);
bp为第j行Ep+1点与Ep点列像素之差,即
bp=E[xp+1]-E[xp] (30);
分别形成n-1个数构成的数列a1、a2、a3、…、ap、…、an-1和b1、b2、b3、…、bp、…、bn-1
R=N0-1 (31);
也即
并令q=1,2,3,……,R,对于第i,j上的点A2q-1、A2q、A2q+1、E2q-1、E2q、E2q+1以q为下标序列,代入求得的ap、bp数列的p下标,用式
求得y1、y2、y3、……、yR,以y1、y2、y3、……、yR为行号按其左侧像素为黑色,右侧像素为白色,检索其对应列坐标,并依次得到列坐标值xCq,q=1,2,3,……、R,令yCq=yq,进而形成中径上点C1、C2、C3、……、CR的列、行坐标值(xCq,yCq);
对于第i,j上的点A2q、A2q+1、A2q+2、E2q、E2q+1、E2q+2以q为下标序列,代入求得的ap、bp数列的p下标,用式
求得y′1、y′2、y′3、……、y′R,以y′1、y′2、y′3、……、y′R为行号按其右侧像素为黑色,左侧像素为白色,
检索其对应列坐标值xC′q,q=1,2,3,……、R,令yC′q=y′q,进而形成中径上点C’1、C’2、C’3、……、C’R的列、行坐标值(xC′q,yC′q);
对于第i′,j′上的点A′1、A′2、A′3、……、A′n和E′1、E′2、E′3、……、E′n,及形成的二维数组A′[xk,yk]和E′[xk,yk],取p=1,2,3,……,n-1,令a′p为第i′行A′p+1点与A′p点列像素之差,即
a′p=A′[xp+1]-A′[xp] (35);
b′p为第j′行E′p+1点与E′p点列像素之差,即
b′p=E′[xp+1]-E′[xp] (36);
分别形成n-1个数构成的数列a′1、a′2、a′3、…、a′p、…、a′n-1和b′1、b′2、b′3、…、b′p、…、b′n-1
R=N0-1 (37);
并令q=1,2,3,……,R,对于第i′,j′上的点A′2q-1、A′2q、A′2q+1、E′2q-1、E′2q、E′2q+1以q为下标序列,代入求得的a′p、b′p数列的p下标,用式
以这些行号按其左侧像素为黑色,而右侧像素为白色检索其对应列坐标xFq,q=1,2,3,……、R,令yFq=y″q,进而形成下半部分螺纹中径上点F1、F2、F3、……、FR的列、行坐标值(xFq,yFq);
对于第i′,j′上的点A′2q、A′2q+1、A′2q+2、E′2q、E′2q+1、E′2q+2以q为下标代替a′p、b′p数列的p下标,用式
求得y″′1、y″′2、y″′3、……、y″′R,以这些行号按其左侧像素为白色,右侧像素为黑色,检索其对应列坐标xF′q,q=1,2,3,……、R,令yF′q=y″′q,进而,形成下半部分螺纹中径上点F’1、F’2、F3、……、F’R的列、行坐标值(xF′q,yF′q)。
3.根据权利要求2所述的一种图像法螺纹中径求解方法,其特征在于:所述步骤3的具体过程如下:
对步骤2中形成的上半部分中径点C1、C2、C3、……、CR及C’1、C’2、C’3、……、C’R,分别以(xCq,yCq)和(xC′q,yC′q)表示其坐标值;对步骤3中形成的下半部分中径点F1、F2、F3、……、FR和F’1、F’2、F’3、……、F’R,则以(xFq,yFq)和(xF′q,yF′q)表示其坐标值,下标q=1,2,3……,R;对用上、下相邻两点Cq、C’q、Fq、F’q点构成的四边形,用列、行坐标按高斯四边形面积公式
和下列公式(40):
求得S1,S2,S3,……,SR之和SS1
对用上、下相邻两点C’q、Cq+1、F’q、Fq+1点构成的四边形,用列、行坐标按高斯四边形面积公式:
和下列公式(41):
求得S’1,S’2,S’3,……,S’R之和SS2
A=SS1+SS2=S1+S2=S1+S2+……+SR+S′1+S′2+……+S′R (20)。
4.根据权利要求3所述的一种图像法螺纹中径求解方法,其特征在于:所述步骤4的具体过程如下:
以C1和C’R的列坐标求得由点C1至C’R距离L0,则可求得螺纹中径d2
式中,d2为螺纹中径,单位:mm;
A为螺纹图像选择区域总像素,单位:平方像素;
L0为螺纹图像选择区域单边像素数,单位:像素;
Ky为y坐标方向上的像素标定系数,单位:mm/像素。
CN201910126315.3A 2019-02-20 2019-02-20 一种图像法螺纹中径求解方法 Active CN109993787B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910126315.3A CN109993787B (zh) 2019-02-20 2019-02-20 一种图像法螺纹中径求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910126315.3A CN109993787B (zh) 2019-02-20 2019-02-20 一种图像法螺纹中径求解方法

Publications (2)

Publication Number Publication Date
CN109993787A true CN109993787A (zh) 2019-07-09
CN109993787B CN109993787B (zh) 2022-10-04

Family

ID=67130218

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910126315.3A Active CN109993787B (zh) 2019-02-20 2019-02-20 一种图像法螺纹中径求解方法

Country Status (1)

Country Link
CN (1) CN109993787B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110260825A (zh) * 2019-07-10 2019-09-20 陕西理工大学 一种图像法螺纹牙型角求解方法
CN110849279A (zh) * 2019-11-27 2020-02-28 陕西理工大学 一种机器视觉测量螺纹中径的补偿方法
CN111462066A (zh) * 2020-03-30 2020-07-28 华南理工大学 一种基于机器视觉的螺纹参数检测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003087524A (ja) * 2001-09-07 2003-03-20 Seiko Epson Corp 画像読み取り装置
CN109003258A (zh) * 2018-06-15 2018-12-14 广东工业大学 一种高精度亚像素圆形零件测量方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003087524A (ja) * 2001-09-07 2003-03-20 Seiko Epson Corp 画像読み取り装置
CN109003258A (zh) * 2018-06-15 2018-12-14 广东工业大学 一种高精度亚像素圆形零件测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈淑娴等: "基于改进SUSAN算法的卡钳排气螺钉参数辨识", 《中国测试》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110260825A (zh) * 2019-07-10 2019-09-20 陕西理工大学 一种图像法螺纹牙型角求解方法
CN110260825B (zh) * 2019-07-10 2020-12-25 陕西理工大学 一种图像法螺纹牙型角求解方法
CN110849279A (zh) * 2019-11-27 2020-02-28 陕西理工大学 一种机器视觉测量螺纹中径的补偿方法
CN110849279B (zh) * 2019-11-27 2021-05-25 陕西理工大学 一种机器视觉测量螺纹中径的补偿方法
CN111462066A (zh) * 2020-03-30 2020-07-28 华南理工大学 一种基于机器视觉的螺纹参数检测方法
CN111462066B (zh) * 2020-03-30 2023-04-21 华南理工大学 一种基于机器视觉的螺纹参数检测方法

Also Published As

Publication number Publication date
CN109993787B (zh) 2022-10-04

Similar Documents

Publication Publication Date Title
CN109993787A (zh) 一种图像法螺纹中径求解方法
EP3278302B1 (de) Bewegungsmesssystem einer maschine und verfahren zum betreiben des bewegungsmesssystems
CN111750789B (zh) 小模数齿轮视觉测量中齿距偏差和齿廓偏差评定方法
CN109409290B (zh) 一种温度表检定读数自动识别系统及方法
CN101067548A (zh) 基于直线相位的光栅自校正方法
DE102016205382A1 (de) Randerkennungsabweichungskorrekturwertberechnungsverfahren, Randerkennungsabweichungskorrekturverfahren und Randerkennungsabweichungskorrekturprogramm
DE112014002505T5 (de) Verfahren zum Selbstkalibrieren eines Drehgebers
EP2388556B1 (de) Optischer Positionsgeber
CN110672037A (zh) 基于相移法的线性光源光栅投影三维测量系统及方法
CN106408616A (zh) 一种ct成像中投影图背景不一致的校正方法
WO2016199439A1 (ja) 三次元計測装置
DE102008055158B4 (de) Verfahren zur 3D-Vermessung der Oberfläche eines Objekts, insbesondere für zahnmedizinische Zwecke
CN113689451A (zh) 一种碳板边界提取方法、装置、存储介质和电子设备
DE112015004528T5 (de) Formmessvorrichtung und Formmessverfahren
CN110926340A (zh) 一种新型光栅尺位移测量装置
CN112634375B (zh) 一种ai智能检测中的平面标定和三维重建方法
CN116664554B (zh) 一种基于图像处理的螺栓螺纹缺陷检测方法
CN115289997B (zh) 一种双目相机三维轮廓扫描仪及其使用方法
CN112146593A (zh) 一种基于机器视觉的外螺纹检测方法及实时检测系统
DE102004058655A1 (de) Verfahren und Anordnung zum Messen von Geometrien eines Objektes mittels eines Koordinatenmessgerätes
CN110260825B (zh) 一种图像法螺纹牙型角求解方法
CN1648600A (zh) 利用单台非量测数码相机进行钣金件三维量测的方法
CN108205811A (zh) 一种植物叶片最小外接矩形计算方法
CN112001379A (zh) 基于机器视觉的汽车仪表固定视点读取仪表的修正算法
CN117788569B (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