CN109993800B - 一种工件尺寸的检测方法、装置及存储介质 - Google Patents

一种工件尺寸的检测方法、装置及存储介质 Download PDF

Info

Publication number
CN109993800B
CN109993800B CN201910201974.9A CN201910201974A CN109993800B CN 109993800 B CN109993800 B CN 109993800B CN 201910201974 A CN201910201974 A CN 201910201974A CN 109993800 B CN109993800 B CN 109993800B
Authority
CN
China
Prior art keywords
image
standard
pixel
sub
matched
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.)
Active
Application number
CN201910201974.9A
Other languages
English (en)
Other versions
CN109993800A (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.)
Wuyi University
Original Assignee
Wuyi University
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 Wuyi University filed Critical Wuyi University
Priority to CN201910201974.9A priority Critical patent/CN109993800B/zh
Publication of CN109993800A publication Critical patent/CN109993800A/zh
Application granted granted Critical
Publication of CN109993800B publication Critical patent/CN109993800B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/02Affine transformations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration

Landscapes

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

Abstract

本发明公开了一种工件尺寸的检测方法、装置及存储介质,包括如下步骤:获取工件的候选图像信息进行相机的标定;创建工件的标准图像和图像金字塔,将所述标定坐标输入到图像金字塔中进行筛选;利用Canny算法和双线性插值法分别计算所述标准图像和匹配图像的亚像素边缘;获取标准亚像素边缘和匹配亚像素边缘之间最小的Hausdorff距离的匹配图像;计算所述标准图像和拟合图像之间的平移矩阵和旋转矩阵,对拟合图像进行仿射变换;将所述校正图像进行边缘分切与拟合;比较所述标准图像和标准拟合图像之间的尺寸差异,在允许误差的范围内选取合格样品;该检测方法及装置能够实时检测工件的质量,检测的速度快、精确度高,提高了产品的合格率。

Description

一种工件尺寸的检测方法、装置及存储介质
技术领域
本发明涉及工件检测领域,具体涉及一种工件尺寸的检测方法、装置及存储介质。
背景技术
随着智能制造技术的发展,工件的制造工艺逐步走向智能化,从原材料到产品出库,一条无人化的生产线能满足所有的工艺需求。在生产过程中,对质量的监控尤为重要,目前针对工件的尺寸检测,主要通过使用米尺进行手工测量,这样的测量方法误差大、精度低、效率下,因此现代工业生产迫切需要新的技术出现,来代替人工的尺寸测量。
发明内容
为解决上述问题,本发明的目的在于提供一种工件尺寸的检测方法及装置,能够实时检测工件的质量,检测的速度快、精确度高,提高了产品的合格率。
本发明解决其问题所采用的技术方案是:第一方面,本发明实施例提出了一种工件尺寸的检测方法,包括如下步骤:
获取工件的候选图像信息,对候选图像信息中的坐标进行相机的标定,得到候选图像的标定坐标;
创建工件的标准图像和用于与标准图像匹配的图像金字塔,将所述候选图像的标定坐标输入到图像金字塔中进行筛选,得到匹配图像;
利用Canny算法和双线性插值法分别计算所述标准图像和匹配图像的亚像素边缘,得到标准图像的标准亚像素边缘,和匹配图像的匹配亚像素边缘;
计算所述标准亚像素边缘和匹配亚像素边缘之间的Hausdorff距离,获取Hausdorff距离最小的匹配图像,得到拟合图像;
计算所述标准图像和拟合图像之间的平移矩阵和旋转矩阵,对拟合图像进行仿射变换,得到校正图像;
将所述校正图像进行边缘分切与拟合,得到标准拟合图像;
比较所述标准图像和标准拟合图像之间的尺寸差异,在允许误差的范围内选取合格样品。
进一步,获取工件的候选图像信息,对候选图像信息中的坐标进行相机的标定,得到候选图像的标定坐标,包括如下步骤:
采集候选图像,获取候选图像中各点的世界坐标,将世界坐标转换为像素坐标;
对所述像素坐标进行畸变矫正,得到标定坐标。
进一步,采集候选图像包括:采用红光隧道型光源采集候选图像。
进一步,创建标准图像和用于与标准图像匹配的图像金字塔,将所述标定坐标输入到图像金字塔中进行筛选,得到匹配图像,包括如下步骤:
采集标准模板并计算标准模板的标定坐标,将标准模板的标定坐标输入到扰动算法,生成标准图像;
提取所述标准图像和候选图像的中心坐标和水平夹角;
将所述标准图像的中心坐标和水平夹角作为图像金字塔的匹配条件,将候选图像的中心坐标和水平夹角输入到图像金字塔进行筛选,得到匹配图像。
进一步,利用Canny算法和双线性插值法分别计算所述标准图像和匹配图像的亚像素边缘,得到标准图像的标准亚像素边缘,和匹配图像的匹配亚像素边缘,包括如下步骤:
将所述标准图像和匹配图像输入到Canny算法,分别计算标准图像和匹配图像的像素级边缘,得到标准图像的像素级边缘的标准像素点,和匹配图像的像素级边缘的匹配像素点;
将所述标准像素点和匹配像素点各自组合为四个互相关联的关联像素点,得到标准像素点的标准关联像素点,和匹配像素点的匹配关联像素点;
将标准关联像素点和匹配关联像素点输入到双线性插值算法中,得到标准图像的标准亚像素边缘和匹配图像的匹配亚像素边缘。
进一步,计算标准图像和拟合图像之间的平移矩阵和旋转矩阵,对拟合图像进行仿射变换,包括如下步骤:
计算拟合图像的中心坐标和水平夹角的大小,得到拟合图像相对于标准图像的平移矩阵和旋转矩阵;
将所述的平移矩阵和旋转矩阵作用于拟合图像上进行仿射变换,得到校正图像。
进一步,计算标准图像和拟合图像之间的平移矩阵和旋转矩阵,对拟合图像进行仿射变换,还包括如下步骤:
将标准图像和校正图像内相邻但未闭合的区域用线段首尾连接,分别计算对比标准图像和校正图像的闭合区域内的亚像素边缘的亚像素单元数量,在亚像素单元数量允许的误差范围内,筛选出校正图像的闭合区域,得到保留真实边缘的校正图像。
进一步,将所述校正图像进行边缘分切与拟合,得到标准拟合图像,包括如下步骤:
采用区域分裂算法将校正图像分割为若干轮廓曲线;
将所述轮廓曲线各自首尾相连,得到连接线段;
设定距离阈值dset,比较轮廓曲线到连接线段的最大距离dmax与阈值dset的值,将最大距离dmax大于阈值dset的轮廓曲线在所在点处分割;
重复上述所述的将轮廓曲线首尾相连,以及分割轮廓曲线的步骤,直到所有的轮廓曲线均满足最大距离dmax小于阈值dset为止;
计算分割完成的轮廓曲线的亚像素边缘的全局特征值,根据全局特征值将轮廓线段拟合成为几何图像;
对所述的几何图像进行最小二乘法的拟合计算,得到标准拟合图像。
第二方面,本发明实施例还提出了一种工件尺寸的检测装置,包括至少一个控制处理器和用于与所述至少一个控制处理器通信连接的存储器;所述存储器存储有可被所述至少一个控制处理器执行的指令,所述指令被所述至少一个控制处理器执行,以使所述至少一个控制处理器能够执行如以上任一项所述的一种工件尺寸的检测方法。
第三方面,本发明实施例还提出了一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机可执行指令,所述计算机可执行指令用于使计算机执行如以上任一项所述的一种工件尺寸的检测方法。
本发明实施例中提供的技术方案,至少具有如下有益效果:通过运用图像金字塔、Canny算法、双线性插值法和Hausdorff距离的结合,实现对候选图像的识别,提高了匹配的精度和速度;对拟合图像进行仿射变换与边缘的分切、拟合,得到的标准拟合图像的拟合程度高,误差率底,有效的提高了产品的合格率。
附图说明
下面结合附图和实例对本发明作进一步说明:
图1是本发明的工件尺寸的检测方法的一个实施例的整体流程图;
图2是本发明的工件尺寸的检测方法的关联像素点一个实施例的示意图;
图3是本发明的工件尺寸的检测方法的校正图像进行边缘分切的一个实施例的示意图。
具体实施方式
在生产过程中,对质量的监控尤为重要,目前针对工件的尺寸检测,主要通过使用米尺进行手工测量,这样的测量方法误差大、精度低、效率下,因此现代工业生产迫切需要新的技术出现,来代替人工的尺寸测量。
基于此,本发明提供了一种工件尺寸的检测方法及装置,能够实时检测工件的质量,检测的速度快、精确度高,提高了产品的合格率。
下面结合附图,对本发明实施例作进一步阐述。
参照图1,本发明的一个实施例提供了一种工件尺寸的检测方法,包括如下步骤:
步骤S100:获取工件的候选图像信息,对候选图像信息中的坐标进行相机的标定,得到候选图像的标定坐标;
步骤S200:创建工件的标准图像和用于与标准图像匹配的图像金字塔,将所述候选图像的标定坐标输入到图像金字塔中进行筛选,得到匹配图像;
步骤S300:利用Canny算法和双线性插值法分别计算所述标准图像和匹配图像的亚像素边缘,得到标准图像的标准亚像素边缘,和匹配图像的匹配亚像素边缘;
步骤S400:计算所述标准亚像素边缘和匹配亚像素边缘之间的Hausdorff距离,获取Hausdorff距离最小的匹配图像,得到拟合图像;
步骤S500:计算所述标准图像和拟合图像之间的平移矩阵和旋转矩阵,对拟合图像进行仿射变换,得到校正图像;
步骤S600:将所述校正图像进行边缘分切与拟合,得到标准拟合图像;
步骤S700:比较所述标准图像和标准拟合图像之间的尺寸差异,在允许误差的范围内选取合格样品。
在本实施例中,步骤S100对候选图像信息进行标定,使候选图像信息更加准确;步骤S200中,图像金字塔的匹配条件是与标准图像的相似程度,例如是标准图像的标定坐标、面积等,通过图像金字塔对候选图像的标定坐标的筛选,可以加快匹配的速度、提高搜索效率;步骤S300利用Canny算法和双线性插值法的结合,计算标准图像的标准亚像素边缘,和匹配图像的匹配亚像素边缘,能够有效的提高执行速度和插值精度。
在步骤S400中,Hausdorff距离是标准亚像素边缘与匹配亚像素边缘间的最大距离h(A,B),以及匹配亚像素边缘与标准亚像素边缘的最大距离h(B,A)的最大值H(A,B):
Figure BDA0001997751930000051
其中,Hausdorff距离值越小,标准亚像素边缘与匹配亚像素边缘就越接近,那么标准图像和匹配图像就越相似,因此拟合图像选取Hausdorff距离差异较小的匹配图像,即H(A,B)的值差异较小,并且Hausdorff距离的计算过程简单,操作方便,拟合程度高。
在步骤S500中,对拟合图像进行仿射变换,使拟合图像通过旋转和平移变换到标准图像的位置,实现了对拟合图像的精确变换与定位,提高了拟合的程度;步骤S600和步骤S700,对校正图像进行边缘分切与拟合,拟合程度高,误差率底,并且在允许误差的范围内选取合格样品,能够有效的提高产品的合格率,优选地,标准拟合图像的允许误差范围为标准图像规格的±2.5%。
进一步地,本发明的另一个实施例还提供了一种工件尺寸的检测方法,其中,获取采集到的候选图像信息,将候选图像信息中的坐标进行相机的标定,包括如下步骤:
步骤S110:采集候选图像,获取候选图像中各点的世界坐标,将世界坐标转换为像素坐标;
步骤S120:对所述像素坐标进行畸变矫正,得到标定坐标。
在本实施例中,步骤S110通过摄像机获取候选图像中各点的世界坐标,将候选图像的世界坐标(XW,YW,ZW)转换为像素坐标(u,v):
Figure BDA0001997751930000061
其中,
Figure BDA0001997751930000062
为摄像机内参矩阵,
Figure BDA0001997751930000063
为摄像机外参矩阵,外参矩阵主要包括一个3*3的旋转矩阵R和3*1的平移矩阵T。在步骤S120中,由于摄像机镜头畸变的存在,会使像素坐标发生变化,例如是径向畸变和切向畸变,因此需要对像素坐标进行畸变矫正:
Figure BDA0001997751930000064
其中,径向半径
Figure BDA0001997751930000065
(k1,k2,k3)是径向畸变参数,(p1,p2)是切向畸变参数,(xcorrected,ycorrected)为畸变校正后的标定坐标。
步骤S110和步骤S120将候选图像的坐标转换像素坐标,并对像素进行了畸变校正,使候选图像的坐标信息更加准确、对候选图像进行比较筛选的过程更加方便,简化了操作流程。
进一步地,本发明的另一个实施例还提供了一种工件尺寸的检测方法,其中,获取候选图像包括:采用红光隧道型光源采集候选图像。
在本实施例中,摄像机在采集候选图像时,使用红光隧道型光源进行采集,因为红光的波长最高,并且接近传感器敏感度的最高值,能够有效避免工件表面成像的反高光问题,使候选图像的清晰度更高。
进一步地,本发明的另一个实施例还提供了一种工件尺寸的检测方法,其中,创建标准图像和用于与标准图像匹配的图像金字塔,将所述标定坐标输入到图像金字塔中进行筛选,得到匹配图像,包括如下步骤:
步骤S210:获取标准模板的标定坐标,将标定坐标输入到扰动算法,生成标准图像;
步骤S220:提取所述标准图像和候选图像的中心坐标和水平夹角;
步骤S230:将所述标准图像的中心坐标和水平夹角作为图像金字塔的匹配条件,将候选图像的中心坐标和水平夹角输入到图像金字塔进行筛选,得到匹配图像。
在本实施例中,步骤S210的标准模板选取形状固定,规格大小在所需范围,且灰度差异较大的图像,即轮廓比较明显的图像作为标准模板,并计算标准模板的标定坐标;扰动算法是线性方法,在标定坐标附近进行迭代,直到迭代结果满足误差允许范围,例如是百万分之一,因此将标定坐标输入到扰动算法中进行计算,生成的标准图像能够减少外界因素对标准模板数据的影响,提高了准确率。在步骤S220中,根据标准图像的标定坐标和候选图像的标定坐标,分别计算标准图像的中心坐标和水平夹角,以及候选图像的中心坐标和水平夹角;步骤S230将标准图像的中心坐标和水平夹角作为图像金字塔的匹配条件,每层的图像金字塔设置的筛选条件为候选图像的中心坐标和水平夹角与标准图像的中心坐标和水平夹角的差异程度,将候选图像的中心坐标和水平夹角输入到图像金字塔中进行筛选,最后经过筛选得到差异范围在较低程度的匹配图像的集合。图像金字塔的应用加快了对候选图像的匹配速度、提高搜索效率。
进一步地,本发明的另一个实施例还提供了一种工件尺寸的检测方法,其中,利用Canny算法和双线性插值法分别计算所述标准图像和匹配图像的亚像素边缘,得到标准图像的标准亚像素边缘,和匹配图像的匹配亚像素边缘,包括如下步骤:
步骤S310:将所述标准图像和匹配图像输入到Canny算法,分别计算标准图像和匹配图像的像素级边缘,得到标准图像的像素级边缘的标准像素点,和匹配图像的像素级边缘的匹配像素点;
步骤S320:将所述标准像素点和匹配像素点各自组合为四个互相关联的关联像素点,得到标准像素点的标准关联像素点,和匹配像素点的匹配关联像素点;
步骤S330:将标准关联像素点和匹配关联像素点输入到双线性插值算法中,得到标准图像的标准亚像素边缘和匹配图像的匹配亚像素边缘。
在本实施例中,步骤S310的Canny算法具有良好的信噪比,对单一边缘具有唯一的响应,高低阈值的选取能优化边缘提取的精度。将标准图像和匹配图像输入到Canny算法,首先利用高斯滤波器对图像进行平滑消噪,二维高斯滤波函数为:
Figure BDA0001997751930000081
其对应的梯度矢量为:
Figure BDA0001997751930000082
其中δ为高斯滤波器的尺度参数,控制着图像的平滑程度,δ越小,边缘定位精度越高,但图像的信噪比降低;δ越大,噪声抑制效果越好,但边缘的定位精度降低。
然后,利用一阶偏导的有限差分来计算Canny算法的梯度幅值和方向:
Figure BDA0001997751930000083
其中
Figure BDA0001997751930000084
梯度幅值A(i,j)反映了图像上(i,j)的边缘强度;
Figure BDA0001997751930000085
为图像点(i,j)的法向矢量,且垂直于边缘的方向。通过判断每一点的梯度方向上A(i,j)是否为其领域的最大值来确定该点是否为边缘点,实现抑制梯度方向上的所有非屋脊峰值的梯度幅值,有效地消除了虚假的边缘点,得到细化的边缘。用高阈值TH和低阈值TL,分别作用于非极大值抑制后的图像,得到两个边缘图像TH(i,j)和TL(i,j),TH(i,j)是用高阈值得到的边缘阵列,含有较少的假边缘,TL(i,j)是用低阈值得到的边缘阵列,保留了较多的细节信息;根据边缘的连通性,以边缘图像阵列TH(i,j)为基础,在边缘图像阵列TL(i,j)中搜索边缘点进行连接,得到既消除了较多假边缘又保留了比较重要边缘信息的像素级边缘,即可得到标准图像的像素级边缘的标准像素点,和匹配图像的像素级边缘的匹配像素点。
参照图2,在步骤S320中,获取标准图像的i行j列的标准像素点(i1,j1),则对于标准像素点(i1,j1)的四个互相关联的标准关联像素点为(i1,j1),(i1+1,j1),(i1,j1+1),(i1+1,j1+1);步骤S330通过双线性插值算法计算标准关联像素点的标准亚像素点(i.j):
f(i,j)=(i1,j1)(j1+1-j)(i1+1-i)+(i1+1,j1)(j1+1-j)(i-i1)+(i1,j1+1)(j-j1)(i1+1-i)+(i1+1,j1+1)(j-j1)(i-i1)。
对匹配图像的匹配像素点同样运行步骤S320和步骤S330,即可得到匹配亚像素点。双线性插值算法的插值精度较高,因此,标准关联像素点和匹配关联像素点经过双线性插值算法运算后,各自联合所有得到的标准亚像素点和匹配亚像素点,即可得到高精度的标准亚像素边缘和匹配亚像素边缘。
进一步地,本发明的另一个实施例还提供了一种工件尺寸的检测方法,其中,计算标准图像和拟合图像之间的平移矩阵和旋转矩阵,对拟合图像进行仿射变换,包括如下步骤:
步骤S510:计算拟合图像的中心坐标和水平夹角的大小,得到拟合图像相对于标准图像的平移矩阵和旋转矩阵;
步骤S520:将所述的平移矩阵和旋转矩阵作用于拟合图像上进行仿射变换,得到校正图像。
在本实施例中,步骤S510通过计算拟合图像的中心坐标(Row,Column)和水平夹角phi的大小,结合标准图像的中心坐标(RefRow,RefColumn)和水平夹角Refphi,得到拟合图像的平移矩阵在X轴方向上的平移效果(RefRow-Row)、Y轴方向上的平移效果(RefColumn-Column),以及旋转矩阵的旋转效果(Refphi-phi),则平移矩阵和旋转矩阵分别为:
Figure BDA0001997751930000091
步骤S520将平移矩阵和旋转矩阵作用于拟合图像上进行仿射变换,使拟合图像通过旋转和平移变换到标准图像的位置,实现了对拟合图像变换的精确定位,提高了拟合的精确度。
进一步地,本发明的另一个实施例还提供了一种工件尺寸的检测方法,其中,计算标准图像和拟合图像之间的平移矩阵和旋转矩阵,对拟合图像进行仿射变换,还包括如下步骤:
步骤S530:将标准图像和校正图像内相邻但未闭合的区域用线段首尾连接,分别计算对比标准图像和校正图像的闭合区域内的亚像素边缘的亚像素单元数量,在亚像素单元数量允许的误差范围内,筛选出校正图像的闭合区域,得到保留真实边缘的校正图像。
在本实施例中,步骤S530通过计算标准图像和校正图像的闭合区域内的亚像素边缘的亚像素单元数量,在允许的误差范围内,筛选出在相同位置的亚像素单元数量近似的闭合区域,剔除数量差异较大的闭合区域,得到保留真实边缘的校正图像,保证了校正图像的准确性。优选地,允许的误差范围为校正图像的亚像素单元数量与标准图像的亚像素单元数量的差异在±1%内。
进一步地,本发明的另一个实施例还提供了一种工件尺寸的检测方法,其中,将所述校正图像进行边缘分切与拟合,得到标准拟合图像,包括如下步骤:
步骤S610:采用区域分裂算法将校正图像分割为若干轮廓曲线;
步骤S620:将所述轮廓曲线各自首尾相连,得到连接线段;
步骤S630:设定距离阈值dset,比较轮廓曲线到连接线段的最大距离dmax与阈值dset的值,将最大距离dmax大于阈值dset的轮廓曲线在所在点处分割;
步骤S640:重复上述所述的将轮廓曲线首尾相连,以及分割轮廓曲线的步骤,直到所有的轮廓曲线均满足最大距离dmax小于阈值dset为止;
步骤S650:计算分割完成的轮廓曲线的亚像素边缘的全局特征值,根据不同的全局特征值将轮廓线段拟合成不同的几何图像;
步骤S660:对所述的几何图像进行最小二乘法的拟合计算,得到标准拟合图像。
在本实施例中,步骤S610利用区域分裂算法对图像中区域特征不一致的区域进行分裂,形成四个相等的子区域,既得到区域特征一致的若干轮廓曲线。
图3为本实施例的边缘分切的示意图,参照图3(a)和图3(b),步骤S620和步骤S630将轮廓曲线到连接线段的最大距离dmax与阈值dset进行比较,将最大距离dmax大于阈值dset的轮廓曲线在所在点处分割,形成两段轮廓曲线,其中,所在点为分割点;参照图3(c)和图3(d),步骤S640为重复步骤S620和步骤S630的操作,即重复将分割的轮廓曲线首尾连接,以及重新分割轮廓曲线,直到所有的轮廓曲线的最大距离dmax均小于阈值dset为止。步骤S650对分割完成的轮廓曲线计算亚像素边缘的全局特征值cont_approx,将cont_approx=1的轮廓曲线用直线拟合,cont_approx=0的轮廓曲线用椭圆形弧拟合,将cont_approx=-1的轮廓曲线用圆拟合,得到初步完成拟合的轮廓曲线,且轮廓曲线拟合程度较高。
步骤S660对不同的曲线轮廓分别进行最小二乘法的拟合,首先根据曲线轮廓假设拟合的曲线轮廓的函数关系,例如是直线、圆形或椭圆的函数关系,然后对函数关系进行最小二乘法处理,在计算最小二乘法的参数时,需要满足拟合值的残差最小,使得实际点与最小二乘法的拟合点之间的距离的平方和达到最小,即可得到关于曲线轮廓的标准轮廓段方程,进而得到标准拟合图像。最小二乘法的计算过程简单,拟合值的精确度高,得到的标准拟合图像的拟合程度也较高,很好地还原了图像的真实度。
另外,本发明的另一个实施例还提供了一种工件尺寸的检测方法,该识别方法包括如下步骤:
步骤S110:采用红光隧道型光源采集候选图像,获取候选图像中各点的世界坐标,以及摄像机的内参矩阵和外参矩阵,将世界坐标转换为像素坐标;
步骤S120:对所述像素坐标进行畸变矫正,得到标定坐标;
步骤S210:获取标准模板的标定坐标,将标定坐标输入到扰动算法,生成标准图像;
步骤S220:提取所述标准图像和候选图像的中心坐标和水平夹角;
步骤S230:将所述标准图像的中心坐标和水平夹角作为图像金字塔的匹配条件,将候选图像的中心坐标和水平夹角输入到图像金字塔进行筛选,得到匹配图像;
步骤S310:将所述标准图像和匹配图像输入到Canny算法,分别计算标准图像和匹配图像的像素级边缘,得到标准图像的像素级边缘的标准像素点,和匹配图像的像素级边缘的匹配像素点;
步骤S320:将所述标准像素点和匹配像素点各自组合为四个互相关联的关联像素点,得到标准像素点的标准关联像素点,和匹配像素点的匹配关联像素点;
步骤S330:将标准关联像素点和匹配关联像素点输入到双线性插值算法中,得到标准图像的标准亚像素边缘和匹配图像的匹配亚像素边缘;
步骤S400:计算所述标准亚像素边缘和匹配亚像素边缘之间的Hausdorff距离,获取Hausdorff距离最小的匹配图像,得到拟合图像;
步骤S510:计算拟合图像的中心坐标和水平夹角的大小,得到拟合图像相对于标准图像的平移矩阵和旋转矩阵;
步骤S520:将所述的平移矩阵和旋转矩阵作用于拟合图像上进行仿射变换,得到校正图像;
步骤S530:将标准图像和校正图像内相邻但未闭合的区域用线段首尾连接,分别计算对比标准图像和校正图像的闭合区域内的亚像素边缘的亚像素单元数量,在亚像素单元数量允许的误差范围内,筛选出校正图像的闭合区域,得到保留真实边缘的校正图像;
步骤S610:采用区域分裂算法将校正图像分割为若干轮廓曲线;
步骤S620:将所述轮廓曲线各自首尾相连,得到连接线段;
步骤S630:设定距离阈值dset,比较轮廓曲线到连接线段的最大距离dmax与阈值dset的值,将最大距离dmax大于阈值dset的轮廓曲线在所在点处分割;
步骤S640:重复上述所述的将轮廓曲线首尾相连,以及分割轮廓曲线的步骤,直到所有的轮廓曲线均满足最大距离dmax小于阈值dset为止;
步骤S650:计算分割完成的轮廓曲线的亚像素边缘的全局特征值,根据不同的全局特征值将轮廓线段拟合成不同的几何图像;
步骤S660:对所述的几何图像进行最小二乘法的拟合计算,得到标准拟合图像;
步骤S700:比较所述标准图像和标准拟合图像之间的尺寸差异,在允许误差的范围内选取合格样品。
在本实施例中,步骤S110通过红光隧道型光源的摄像机获取候选图像中各点的世界坐标,红光隧道型光源的灵敏度高,能够有效避免工件表面成像的反高光问题,使候选图像的清晰度更高。将候选图像的世界坐标(XW,YW,ZW)转换为像素坐标(u,v)的公式为:
Figure BDA0001997751930000131
其中,
Figure BDA0001997751930000132
为摄像机内参矩阵,
Figure BDA0001997751930000133
为摄像机外参矩阵,外参矩阵主要包括一个3*3的旋转矩阵R和3*1的平移矩阵T。在步骤S120中,由于摄像机镜头畸变的存在,会使像素坐标发生变化,例如是径向畸变和切向畸变,因此需要对像素坐标进行畸变矫正:
Figure BDA0001997751930000134
其中,径向半径
Figure BDA0001997751930000135
(k1,k2,k3)是径向畸变参数,(p1,p2)是切向畸变参数,(xcorrected,ycorrected)为畸变校正后的标定坐标。
步骤S110和步骤S120将候选图像的坐标转换像素坐标,并对像素进行了畸变校正,使候选图像的坐标信息更加准确、对候选图像进行比较筛选的过程更加方便,简化了操作流程。
步骤S210的标准模板选取形状固定,规格大小在所需范围,且灰度差异较大的图像,即轮廓比较明显的图像作为标准模板,并计算标准模板的标定坐标;扰动算法是线性方法,在标定坐标附近进行迭代,直到迭代结果满足误差允许范围,例如是百万分之一,因此将标定坐标输入到扰动算法中进行计算,生成的标准图像能够减少外界因素对标准模板数据的影响,提高了准确率。在步骤S220中,根据标准图像的标定坐标和候选图像的标定坐标,分别计算标准图像的中心坐标和水平夹角,以及候选图像的中心坐标和水平夹角;步骤S230将标准图像的中心坐标和水平夹角作为图像金字塔的匹配条件,每层的图像金字塔设置的筛选条件为候选图像的中心坐标和水平夹角与标准图像的中心坐标和水平夹角的差异程度,将候选图像的中心坐标和水平夹角输入到图像金字塔中进行筛选,最后经过筛选得到差异范围在较低程度的匹配图像的集合。图像金字塔的应用加快了对候选图像的匹配速度、提高搜索效率。
步骤S310的Canny算法具有良好的信噪比,对单一边缘具有唯一的响应,高低阈值的选取能优化边缘提取的精度。将标准图像和匹配图像输入到Canny算法,首先利用高斯滤波器对图像进行平滑消噪,二维高斯滤波函数为:
Figure BDA0001997751930000141
其对应的梯度矢量为:
Figure BDA0001997751930000142
其中δ为高斯滤波器的尺度参数,控制着图像的平滑程度,δ越小,边缘定位精度越高,但图像的信噪比降低;δ越大,噪声抑制效果越好,但边缘的定位精度降低。
然后,利用一阶偏导的有限差分来计算Canny算法的梯度幅值和方向:
Figure BDA0001997751930000143
其中
Figure BDA0001997751930000144
梯度幅值A(i,j)反映了图像上(i,j)的边缘强度;
Figure BDA0001997751930000145
为图像点(i,j)的法向矢量,且垂直于边缘的方向。通过判断每一点的梯度方向上A(i,j)是否为其领域的最大值来确定该点是否为边缘点,实现抑制梯度方向上的所有非屋脊峰值的梯度幅值,有效地消除了虚假的边缘点,得到细化的边缘。用高阈值TH和低阈值TL,分别作用于非极大值抑制后的图像,得到两个边缘图像TH(i,j)和TL(i,j),TH(i,j)是用高阈值得到的边缘阵列,含有较少的假边缘,TL(i,j)是用低阈值得到的边缘阵列,保留了较多的细节信息;根据边缘的连通性,以边缘图像阵列TH(i,j)为基础,在边缘图像阵列TL(i,j)中搜索边缘点进行连接,得到既消除了较多假边缘又保留了比较重要边缘信息的像素级边缘,即可得到标准图像的像素级边缘的标准像素点,和匹配图像的像素级边缘的匹配像素点。
参照图2,在步骤S320中,获取标准图像的i行j列的标准像素点(i1,j1),则对于标准像素点(i1,j1)的四个互相关联的标准关联像素点为(i1,j1),(i1+1,j1),(i1,j1+1),(i1+1,j1+1);步骤S330通过双线性插值算法计算标准关联像素点的标准亚像素点(i.j):
f(i,j)=(i1,j1)(j1+1-j)(i1+1-i)+(i1+1,j1)(j1+1-j)(i-i1)+(i1,j1+1)(j-j1)(i1+1-i)+(i1+1,j1+1)(j-j1)(i-i1)。
对匹配图像的匹配像素点同样运行步骤S320和步骤S330,即可得到匹配亚像素点。双线性插值算法的插值精度较高,因此,标准关联像素点和匹配关联像素点经过双线性插值算法运算后,各自联合所有得到的标准亚像素点和匹配亚像素点,即可得到高精度的标准亚像素边缘和匹配亚像素边缘。
在步骤S400中,计算标准亚像素边缘与匹配亚像素边缘间的Hausdorff距离,Hausdorff距离是标准亚像素边缘与匹配亚像素边缘间的最大距离h(A,B),以及匹配亚像素边缘与标准亚像素边缘的最大距离h(B,A)的最大值H(A,B):
Figure BDA0001997751930000151
其中,Hausdorff距离值越小,标准亚像素边缘与匹配亚像素边缘就越接近,那么标准图像和匹配图像就越相似,因此拟合图像选取Hausdorff距离差异较小的匹配图像,即H(A,B)的值差异较小,并且Hausdorff距离的计算过程简单,操作方便,拟合程度高。
步骤S510通过计算拟合图像的中心坐标(Row,Column)和水平夹角phi的大小,结合标准图像的中心坐标(RefRow,RefColumn)和水平夹角Refphi,得到拟合图像的平移矩阵在X轴方向上的平移效果(RefRow-Row)、Y轴方向上的平移效果(RefColumn-Column),以及旋转矩阵的旋转效果(Refphi-phi),则平移矩阵和旋转矩阵分别为:
Figure BDA0001997751930000161
步骤S520将平移矩阵和旋转矩阵作用于拟合图像上进行仿射变换,使拟合图像通过旋转和平移变换到标准图像的位置,实现了对拟合图像变换的精确定位,提高了拟合的精确度。
步骤S530通过计算标准图像和校正图像的闭合区域内的亚像素边缘的亚像素单元数量,在允许的误差范围内,筛选出在相同位置的亚像素单元数量近似的闭合区域,剔除数量差异较大的闭合区域,得到保留真实边缘的校正图像,保证了校正图像的准确性。优选地,允许的误差范围为校正图像的亚像素单元数量与标准图像的亚像素单元数量的差异在±1%内。
步骤S610利用区域分裂算法对图像中区域特征不一致的区域进行分裂,形成四个相等的子区域,既得到区域特征一致的若干轮廓曲线。
图3为本实施例的边缘分切的示意图,参照图3(a)和图3(b),步骤S620和步骤S630将轮廓曲线到连接线段的最大距离dmax与阈值dset进行比较,将最大距离dmax大于阈值dset的轮廓曲线在所在点处分割,形成两段轮廓曲线,其中,所在点为分割点;参照图3(c)和图3(d),步骤S640为重复步骤S620和步骤S630的操作,即重复将分割的轮廓曲线首尾连接,以及重新分割轮廓曲线,直到所有的轮廓曲线的最大距离dmax均小于阈值dset为止。步骤S650对分割完成的轮廓曲线计算亚像素边缘的全局特征值cont_approx,将cont_approx=1的轮廓曲线用直线拟合,cont_approx=0的轮廓曲线用椭圆形弧拟合,将cont_approx=-1的轮廓曲线用圆拟合,得到初步完成拟合的轮廓曲线,且轮廓曲线拟合程度较高。
步骤S660对不同的曲线轮廓分别进行最小二乘法的拟合,首先根据曲线轮廓假设拟合的曲线轮廓的函数关系,例如是直线、圆形或椭圆的函数关系,然后对函数关系进行最小二乘法处理,在计算最小二乘法的参数时,需要满足拟合值的残差最小,使得实际点与最小二乘法的拟合点之间的距离的平方和达到最小,即可得到关于曲线轮廓的标准轮廓段方程,进而得到标准拟合图像。最小二乘法的计算过程简单,拟合值的精确度高,得到的标准拟合图像的拟合程度也较高,很好地还原了图像的真实度。
步骤S700在允许误差的范围内选取合格样品,能够有效的提高产品的合格率,优选地,标准拟合图像的允许误差范围为标准图像规格的±2.5%。
根据上述步骤,特别地选取下水槽来进行实验,已知下水槽的外边缘尺寸为500mm*440mm,内槽边缘尺寸为400mm*400mm,针对下水槽的已知尺寸,对下水槽进行人工多次测量和基于工件的尺寸检测方法进行测量,测量结果如下表1和表2所示:
表1:人工测量下水槽尺寸结果
Figure BDA0001997751930000171
表2:基于本发明的下水槽尺寸测量结果
Figure BDA0001997751930000172
Figure BDA0001997751930000181
根据表1和表2的数据可知,人工测量的精度较低,特别是在对角线的测量时,由于存在悬空测量的过程,极易产生超出误差允许范围的情况,导致产品的误检。相比较于本文的基于工件的尺寸检测方法,该方法只在成像平面进行测量,对对角线和对边缘的测量精度是一致的,所以平均误差值较低,可靠性强,精度较高。并且,人工测量下水槽的一条边缘的时间大概是5s,测量一条对角线大概是9s,那么对一个下水槽的完整的测量时间至少需要80s,而工件的尺寸检测方法中,对下水槽的所有边缘和对角线的测量几乎同时进行,因此对一个下水槽的完整的测量时间只需要1.2s。综上所述,工件的尺寸检测方法能够有效地提高检测的精度和速率,很好地提高了产品的合格率。
此外,本发明的另一个实施例还提供了一种工件尺寸的检测装置,包括至少一个控制处理器和用于与所述至少一个控制处理器通信连接的存储器;所述存储器存储有可被所述至少一个控制处理器执行的指令,所述指令被所述至少一个控制处理器执行,以使所述至少一个控制处理器能够执行如上的任一项所述的一种工件尺寸的检测方法。
在本实施例中,检测装置包括:一个或多个控制处理器和存储器,控制处理器和存储器可以通过总线或者其他方式连接。
存储器作为一种非暂态计算机可读存储介质,可用于存储非暂态软件程序、非暂态性计算机可执行程序以及模块,如本发明实施例中的检测方法对应的程序指令/模块。控制处理器通过运行存储在存储器中的非暂态软件程序、指令以及模块,从而执行检测装置的各种功能应用以及数据处理,即实现上述方法实施例的检测方法。
存储器可以包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需要的应用程序;存储数据区可存储根据检测装置的使用所创建的数据等。此外,存储器可以包括高速随机存取存储器,还可以包括非暂态存储器,例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施方式中,存储器可选包括相对于控制处理器远程设置的存储器,这些远程存储器可以通过网络连接至该检测装置。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
所述一个或者多个模块存储在所述存储器中,当被所述一个或者多个控制处理器执行时,执行上述方法实施例中的检测方法,例如,执行以上描述检测方法步骤S100至S700、S110至S120、S210至S230、S310至S320、S510至S530,以及S610至S660的功能。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机可执行指令,该计算机可执行指令被一个或多个控制处理器执行,例如,一个控制处理器执行,可使得上述一个或多个控制处理器执行上述方法实施例中的检测方法,例如,执行以上描述的图1中的方法步骤S100至S700、S110至S120、S210至S230、S310至S320、S510至S530,以及S610至S660的功能。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
通过以上的实施方式的描述,本领域技术人员可以清楚地了解到各实施方式可借助软件加通用硬件平台的方式来实现。本领域技术人员可以理解实现上述实施例方法中的全部或部分流程是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(ReadOnly Memory,ROM)或随机存储记忆体(Random Access Memory,RAM)等。
以上是对本发明的较佳实施进行了具体说明,但本发明并不局限于上述实施方式,熟悉本领域的技术人员在不违背本发明精神的前提下还可作出种种的等同变形或替换,这些等同的变形或替换均包含在本申请权利要求所限定的范围内。

Claims (7)

1.一种工件尺寸的检测方法,其特征在于:包括如下步骤:
获取工件的候选图像信息,对候选图像信息中的坐标进行相机的标定,得到候选图像的标定坐标;
采集标准模板并计算标准模板的标定坐标,将标准模板的标定坐标输入到扰动算法,生成标准图像;
提取所述标准图像和候选图像的中心坐标和水平夹角;
创建图像金字塔,将所述标准图像的中心坐标和水平夹角作为图像金字塔的匹配条件,将候选图像的中心坐标和水平夹角输入到图像金字塔进行筛选,得到匹配图像;
利用Canny算法和双线性插值法分别计算所述标准图像和匹配图像的亚像素边缘,得到标准图像的标准亚像素边缘,和匹配图像的匹配亚像素边缘;
计算所述标准亚像素边缘和匹配亚像素边缘之间的Hausdorff距离,获取Hausdorff距离最小的匹配图像,得到拟合图像;
计算拟合图像的中心坐标和水平夹角的大小,得到拟合图像相对于标准图像的平移矩阵和旋转矩阵;
将所述的平移矩阵和旋转矩阵作用于拟合图像上进行仿射变换,得到校正图像;
采用区域分裂算法将校正图像分割为若干轮廓曲线;
首尾相连线段:将所述轮廓曲线各自首尾相连,得到连接线段;
分割轮廓曲线:设定距离阈值dset,比较轮廓曲线到连接线段的最大距离dmax与阈值dset的值,将最大距离dmax大于阈值dset的轮廓曲线在所在点处分割;
重复上述首尾相连线段和分割轮廓曲线的步骤,直到所有的轮廓曲线均满足最大距离dmax小于阈值dset为止;
计算分割完成的轮廓曲线的亚像素边缘的全局特征值,根据全局特征值将轮廓线段拟合成为几何图像;
对所述的几何图像进行最小二乘法的拟合计算,得到标准拟合图像;
比较所述标准图像和标准拟合图像之间的尺寸差异,在允许误差的范围内选取合格样品。
2.根据权利要求1所述的一种工件尺寸的检测方法,其特征在于:获取工件的候选图像信息,对候选图像信息中的坐标进行相机的标定,得到候选图像的标定坐标,包括如下步骤:
采集候选图像,获取候选图像中各点的世界坐标,将世界坐标转换为像素坐标;
对所述像素坐标进行畸变矫正,得到标定坐标。
3.根据权利要求2所述的一种工件尺寸的检测方法,其特征在于:采集候选图像包括:采用红光隧道型光源采集候选图像。
4.根据权利要求1所述的一种工件尺寸的检测方法,其特征在于:利用Canny算法和双线性插值法分别计算所述标准图像和匹配图像的亚像素边缘,得到标准图像的标准亚像素边缘,和匹配图像的匹配亚像素边缘,包括如下步骤:
将所述标准图像和匹配图像输入到Canny算法,分别计算标准图像和匹配图像的像素级边缘,得到标准图像的像素级边缘的标准像素点,和匹配图像的像素级边缘的匹配像素点;
将所述标准像素点和匹配像素点各自组合为四个互相关联的关联像素点,得到标准像素点的标准关联像素点,和匹配像素点的匹配关联像素点;
将标准关联像素点和匹配关联像素点输入到双线性插值算法中,得到标准图像的标准亚像素边缘和匹配图像的匹配亚像素边缘。
5.根据权利要求1所述的一种工件尺寸的检测方法,其特征在于:计算标准图像和拟合图像之间的平移矩阵和旋转矩阵,对拟合图像进行仿射变换,还包括如下步骤:
将标准图像和校正图像内相邻但未闭合的区域用线段首尾连接,分别计算对比标准图像和校正图像的闭合区域内的亚像素边缘的亚像素单元数量,在亚像素单元数量允许的误差范围内,筛选出校正图像的闭合区域,得到保留真实边缘的校正图像。
6.一种工件尺寸的检测装置,其特征在于,包括至少一个控制处理器和用于与所述至少一个控制处理器通信连接的存储器;所述存储器存储有可被所述至少一个控制处理器执行的指令,所述指令被所述至少一个控制处理器执行,以使所述至少一个控制处理器能够执行如权利要求1-5任一项所述的一种工件尺寸的检测方法。
7.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机可执行指令,所述计算机可执行指令用于使计算机执行如权利要求1-5任一项所述的一种工件尺寸的检测方法。
CN201910201974.9A 2019-03-18 2019-03-18 一种工件尺寸的检测方法、装置及存储介质 Active CN109993800B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910201974.9A CN109993800B (zh) 2019-03-18 2019-03-18 一种工件尺寸的检测方法、装置及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910201974.9A CN109993800B (zh) 2019-03-18 2019-03-18 一种工件尺寸的检测方法、装置及存储介质

Publications (2)

Publication Number Publication Date
CN109993800A CN109993800A (zh) 2019-07-09
CN109993800B true CN109993800B (zh) 2023-01-03

Family

ID=67129536

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910201974.9A Active CN109993800B (zh) 2019-03-18 2019-03-18 一种工件尺寸的检测方法、装置及存储介质

Country Status (1)

Country Link
CN (1) CN109993800B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111445431B (zh) * 2018-12-28 2023-10-20 Tcl科技集团股份有限公司 一种图像分割方法、终端设备及计算机可读存储介质
CN110544276A (zh) * 2019-08-19 2019-12-06 西安交通大学 最小二乘法椭圆拟合活塞裙部最大点处尺寸测量方法
CN110675393B (zh) * 2019-09-29 2024-03-26 上海呈彧智能科技有限公司 一种基于机器视觉的坯料规格检测方法
CN111524099A (zh) * 2020-04-09 2020-08-11 武汉钢铁有限公司 评价试样截面几何形状参数的方法
CN111860501B (zh) * 2020-07-14 2021-02-05 哈尔滨市科佳通用机电股份有限公司 基于形状匹配的高铁高度调整杆脱出故障图像识别方法
CN113114830B (zh) * 2021-06-11 2021-08-24 成都新西旺自动化科技有限公司 对位位置的校正量获取方法、调节方法、贴合方法及介质
CN113514007B (zh) * 2021-06-30 2022-04-19 微见智能封装技术(深圳)有限公司 一种工件尺寸的测量方法及装置
CN113554616A (zh) * 2021-07-21 2021-10-26 深圳市蓝蓝科技有限公司 基于数控机床的在线测量引导方法及系统
CN113674214A (zh) * 2021-07-26 2021-11-19 常州市新创智能科技有限公司 一种基于双目相机拼接的直线度检测方法和装置
CN114331951B (zh) * 2021-09-30 2024-09-24 腾讯科技(深圳)有限公司 图像检测方法、装置、计算机、可读存储介质及程序产品
CN114022569B (zh) * 2021-11-18 2024-06-07 湖北中烟工业有限责任公司 一种基于视觉测量箱体方正度的方法及装置
CN115018835B (zh) * 2022-08-08 2022-11-04 聊城市宁泰电机有限公司 一种汽车起动机齿轮检测方法
CN116519688B (zh) * 2023-04-26 2024-05-14 中国科学院植物研究所 一种浆果表型特征高通量获取及自动分析方法和系统
CN116580022B (zh) * 2023-07-07 2023-09-29 杭州鄂达精密机电科技有限公司 工件尺寸检测方法、装置、计算机设备及存储介质
CN118470021B (zh) * 2024-07-12 2024-09-27 青岛珞宾通信有限公司 一种基于人工智能的施工质量评价方法

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010097436A (ja) * 2008-10-16 2010-04-30 Keyence Corp 画像処理におけるパターンモデルの位置決め方法、画像処理装置、画像処理プログラム及びコンピュータで読み取り可能な記録媒体
CN102107374A (zh) * 2011-02-01 2011-06-29 吉林大学 轴盘类零件直径尺寸在线检测方法
CN104346788A (zh) * 2013-07-29 2015-02-11 展讯通信(上海)有限公司 一种拼接图像的方法和装置
CN104359403A (zh) * 2014-11-21 2015-02-18 天津工业大学 基于亚像素边缘算法的平面零件尺寸测量方法
CN105740899A (zh) * 2016-01-29 2016-07-06 长安大学 一种机器视觉图像特征点检测与匹配复合的优化方法
CN105930858A (zh) * 2016-04-06 2016-09-07 吴晓军 一种带旋转、缩放的快速高精度几何模板匹配方法
DE102016114337A1 (de) * 2015-08-06 2017-02-09 Cognex Corporation System und verfahren zum verknüpfen von koordinatenräumen maschinellen sehens in einer umgebung angeleiteten zusammenbaus
CN106824806A (zh) * 2016-12-27 2017-06-13 聊城大学 基于机器视觉的小模数塑料齿轮的检测系统及检测方法
CN106969706A (zh) * 2017-04-02 2017-07-21 聊城大学 基于双目立体视觉的工件检测与三维测量系统及检测方法
CN107301661A (zh) * 2017-07-10 2017-10-27 中国科学院遥感与数字地球研究所 基于边缘点特征的高分辨率遥感图像配准方法
CN107388991A (zh) * 2017-07-03 2017-11-24 中国计量大学 一种端面多圆角轴类零件圆角半径测量方法
CN107588721A (zh) * 2017-08-28 2018-01-16 武汉科技大学 一种基于双目视觉的零件多尺寸的测量方法及系统
CN108648205A (zh) * 2018-05-07 2018-10-12 广州大学 一种亚像素边缘检测方法
CN108986126A (zh) * 2018-06-15 2018-12-11 哈尔滨工业大学 基于高斯拟合亚像素边缘检测及改进ransac算法的圆心检测法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7349922B2 (en) * 2001-11-14 2008-03-25 Yeda Research And Development Co. Ltd. Method and apparatus for data clustering including segmentation and boundary detection
US8175412B2 (en) * 2004-02-17 2012-05-08 Yeda Research & Development Co. Ltd. Method and apparatus for matching portions of input images
JP5301239B2 (ja) * 2008-08-09 2013-09-25 株式会社キーエンス 画像処理におけるパターンモデルの位置決め方法、画像処理装置、画像処理プログラム及びコンピュータで読み取り可能な記録媒体
JP5685665B1 (ja) * 2014-04-28 2015-03-18 ヴィスコ・テクノロジーズ株式会社 画像認識装置
JP6570370B2 (ja) * 2015-08-21 2019-09-04 キヤノン株式会社 画像処理方法、画像処理装置、プログラム及び記録媒体
JP2017162391A (ja) * 2016-03-11 2017-09-14 東芝メモリ株式会社 画像処理方法及び画像処理プログラム
CN107194959A (zh) * 2017-04-25 2017-09-22 北京海致网聚信息技术有限公司 基于切片进行图像配准的方法和装置
EP3404583A1 (en) * 2017-05-19 2018-11-21 MVTec Software GmbH System and method for model adaptation
CN108537831B (zh) * 2018-03-09 2021-06-15 中北大学 对增材制造工件进行ct成像的方法及装置

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010097436A (ja) * 2008-10-16 2010-04-30 Keyence Corp 画像処理におけるパターンモデルの位置決め方法、画像処理装置、画像処理プログラム及びコンピュータで読み取り可能な記録媒体
CN102107374A (zh) * 2011-02-01 2011-06-29 吉林大学 轴盘类零件直径尺寸在线检测方法
CN104346788A (zh) * 2013-07-29 2015-02-11 展讯通信(上海)有限公司 一种拼接图像的方法和装置
CN104359403A (zh) * 2014-11-21 2015-02-18 天津工业大学 基于亚像素边缘算法的平面零件尺寸测量方法
DE102016114337A1 (de) * 2015-08-06 2017-02-09 Cognex Corporation System und verfahren zum verknüpfen von koordinatenräumen maschinellen sehens in einer umgebung angeleiteten zusammenbaus
CN105740899A (zh) * 2016-01-29 2016-07-06 长安大学 一种机器视觉图像特征点检测与匹配复合的优化方法
CN105930858A (zh) * 2016-04-06 2016-09-07 吴晓军 一种带旋转、缩放的快速高精度几何模板匹配方法
CN106824806A (zh) * 2016-12-27 2017-06-13 聊城大学 基于机器视觉的小模数塑料齿轮的检测系统及检测方法
CN106969706A (zh) * 2017-04-02 2017-07-21 聊城大学 基于双目立体视觉的工件检测与三维测量系统及检测方法
CN107388991A (zh) * 2017-07-03 2017-11-24 中国计量大学 一种端面多圆角轴类零件圆角半径测量方法
CN107301661A (zh) * 2017-07-10 2017-10-27 中国科学院遥感与数字地球研究所 基于边缘点特征的高分辨率遥感图像配准方法
CN107588721A (zh) * 2017-08-28 2018-01-16 武汉科技大学 一种基于双目视觉的零件多尺寸的测量方法及系统
CN108648205A (zh) * 2018-05-07 2018-10-12 广州大学 一种亚像素边缘检测方法
CN108986126A (zh) * 2018-06-15 2018-12-11 哈尔滨工业大学 基于高斯拟合亚像素边缘检测及改进ransac算法的圆心检测法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于HALCON的双目立体视觉工件尺寸测量;张如如等;《计算机测量与控制》;20180125(第01期);第59-63页 *
微装配中显微视觉关键技术研究;黄海亭;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20190215(第02期);第B022-687页 *
铝合金零件尺寸误差自动检测方法研究;冯淼等;《机电工程》;20180420(第04期);第380-384页 *

Also Published As

Publication number Publication date
CN109993800A (zh) 2019-07-09

Similar Documents

Publication Publication Date Title
CN109993800B (zh) 一种工件尺寸的检测方法、装置及存储介质
CN107798326B (zh) 一种轮廓视觉检测方法
CN110866924B (zh) 一种线结构光中心线提取方法及存储介质
CN105021124B (zh) 一种基于深度图的平面零件三维位置和法向量计算方法
Yang et al. Three-dimensional measurement of precise shaft parts based on line structured light and deep learning
CN108830868B (zh) 一种基于Snake模型和迭代极性转换回归结合的圆弧拟合方法
CN110906875B (zh) 一种孔径测量的视觉处理方法
CN111123242B (zh) 一种基于激光雷达和相机的联合标定方法及计算机可读存储介质
CN103886597A (zh) 一种基于边缘检测与拟合曲线聚类的圆检测方法
CN109448059B (zh) 一种快速x角点亚像素检测方法
CN118229673B (zh) 一种基于图像识别的管道内壁缺陷检测方法及系统
CN111667470A (zh) 一种基于数字图像的工业管道探伤内壁检测方法
Xia et al. A table method for coded target decoding with application to 3-D reconstruction of soil specimens during triaxial testing
CN116579955B (zh) 一种新能源电芯焊缝反光点去噪和点云补全方法及系统
Guo et al. A V-shaped weld seam measuring system for large workpieces based on image recognition
CN116819561A (zh) 一种点云数据匹配方法、系统、电子设备及存储介质
Ziqiang et al. Research of the algorithm calculating the length of bridge crack based on stereo vision
EP4430837A1 (en) Systems and methods for draft calculation
CN105260989B (zh) 基于多图配准复原三维图像的方法
CN109035335B (zh) 一种基于单目视觉的海底隧道渗水水位识别方法
CN209279912U (zh) 一种物体三维信息采集装置
JP6639369B2 (ja) テンプレートマッチング装置、方法、及びプログラム
CN118429345B (zh) 基于图像分割的石油管道异常状态数据提取方法
CN116521927B (zh) 基于网络地图瓦片的遥感影像匹配方法与系统
CN109584287B (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