CN109612390A - 基于机器视觉的大尺寸工件自动测量系统 - Google Patents
基于机器视觉的大尺寸工件自动测量系统 Download PDFInfo
- Publication number
- CN109612390A CN109612390A CN201811541212.5A CN201811541212A CN109612390A CN 109612390 A CN109612390 A CN 109612390A CN 201811541212 A CN201811541212 A CN 201811541212A CN 109612390 A CN109612390 A CN 109612390A
- Authority
- CN
- China
- Prior art keywords
- image
- straight line
- size
- size workpiece
- edge
- 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
- 238000005259 measurement Methods 0.000 claims abstract description 68
- 238000001514 detection method Methods 0.000 claims abstract description 40
- 238000006243 chemical reaction Methods 0.000 claims abstract description 18
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 238000003708 edge detection Methods 0.000 claims description 30
- 238000000034 method Methods 0.000 claims description 25
- 230000006872 improvement Effects 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 6
- 238000005286 illumination Methods 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 3
- 230000033001 locomotion Effects 0.000 abstract description 9
- 238000004519 manufacturing process Methods 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 7
- 229910000831 Steel Inorganic materials 0.000 description 5
- 239000010959 steel Substances 0.000 description 5
- 238000000691 measurement method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000003705 background correction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000036651 mood Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- 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/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- 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/10—Image acquisition modality
- G06T2207/10024—Color image
-
- 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/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- 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/20—Special algorithmic details
- G06T2207/20092—Interactive image processing based on input by user
- G06T2207/20104—Interactive definition of region of interest [ROI]
-
- 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/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
-
- 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)
- Geometry (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于机器视觉的大尺寸工件自动测量系统,属于机器视觉测量领域。包括大尺寸工件自动测量硬件系统和大尺寸工件自动测量软件系统;所述大尺寸工件自动测量硬件系统,包括光源、相机和镜头以及高精度运动装置;所述大尺寸工件自动测量软件系统,包括图像预处理算法、图像数据转换建模、基于改进RANSAC的亚像素直线/圆弧检测算法和尺寸计算。本发明的系统能够针对大尺寸的工件进行测量;填补了机器视觉在大尺寸工件测量领域应用的空白。本发明中的软件系统能够对图像有效的去除噪声、无效点,并且还能对非均匀背景进行校正;提高了本发明的系统的生产效率、系统的稳定型、系统测量准确率。
Description
技术领域
本发明涉及一种基于机器视觉的大尺寸工件自动测量系统,属于视觉测量领域。
背景技术
近年来,随着航空航天、大型船舶制造、采矿及新能源工业的快速发展,大尺寸工件在市场上的需求日益增多,测量精度的要求也越来越高。目前,大尺寸工件的测量方式主要是人工测量,但是这种测量方式存在很多的问题:(1)人工测量精度低;(2)实时性差,效率低,产品的成本提高;(3)稳定性不高,人工测量精度往往会受到操作者本身情绪的影响。因此,基于机器视觉的测量方法得到了快速发展。
机器视觉是利用计算机对图像或视频进行处理,实现对二维或三维场景的识别、检测、感知和理解等功能。目前,机器视觉已广泛应用于工业测量领域,而且具有非接触、高效、可靠、通用性好、易于实现自动化等优点。因此,机器视觉在大批量、连续自动化生产流水线上具有很高的应用价值,能够对产品进行辨识和外形尺寸测量等。
现有的机器视觉测量工件的方法,大多只能满足小尺寸工件的测量,没有测量大尺寸工件的方法。并且现有的机器视觉测量方法对于图像的处理大部分都是先采用高斯滤波进行预处理,再采用通过特征点将图像拼接成一张大图的方法来进行图像拼接,最后对图像进行测量;这样的处理方法在预处理阶段只能消除高斯噪声,不能对非均匀背景起到校正,不能完全的去除无效点和干扰点对结果的影响;现有的图像拼接的方法效率低,测量误差大,尤其是在特征点不明显的情况下更加突出。
发明内容
本发明针对现有技术存在的上述不足,提出一种基于机器视觉的大尺寸工件自动测量系统。本发明的系统硬件稳定,算法鲁棒性好,能够解决目前测量方法效率低下、误差大的缺点,能够用于大尺寸工件的自动测量,具有较好的实用价值。
本发明是通过以下技术方案实现的:
一种基于机器视觉的大尺寸工件自动测量系统,包括大尺寸工件自动测量硬件系统和大尺寸工件自动测量软件系统。
所述大尺寸工件自动测量硬件系统,包括光源、相机和镜头以及高精度运动装置。
所述大尺寸工件自动测量软件系统,包括图像预处理算法、图像数据转换建模、基于改进RANSAC的亚像素直线/圆弧检测算法和尺寸计算。
所述的图像预处理算法包括:(1)灰度化操作,将彩色图像转化为灰度图;(2)对步骤(1)得到的灰度图进行高斯滤波,得到边缘保留完整且对比度高的图像;(3)通过Canny边缘检测得到待测区域的粗边缘信息。
所述的图像数据转换建模是通过坐标转换来获得所需尺寸的算法,包括:
(a)通过光电编码器精确定位相机位置,采集到多幅工件图像,以第一张图像为基准图像;
(b)利用式(2)和式(3)所示的坐标转换模型将其余图像的坐标转换为基准图像下的坐标;
i=i′+(in-i0)×xi (2)
j=j′+(jn-j0)×yj (3)
式中,每幅图像的位置索引为(i,j),基准图像的索引为(1,1),其中i、j是转换之后的坐标,i'、j'是待转换图像自身的像素坐标,in、jn是待转换图像的位置索引,i0、j0是基准图像的索引,xi、yj是图像在i和j方向上的像素个数。
所述的基于改进RANSAC的亚像素直线/圆弧检测算法,包括基于三次样条插值的亚像素边缘检测算法及基于改进RANSAC的直线/圆弧检测算法。
所述的基于三次样条插值的亚像素边缘检测算法包括:
①使用Canny边缘检测算子对所选的ROI区域进行边缘检测得到粗边缘;
②以粗边缘上某一点为中心沿垂直于边缘的左右方向或上下方向各取三个像素点并得到其灰度值;
③对步骤②中取得的六个点和中心点进行三次样条插值,得到边界点的光照强度函数,为六个三次多项式S(x);
④求解函数S(x)二阶导数的过零点,所述的过零点为亚像素级的边界点;然后对求得的亚像素边界点重新进行边缘拟合。
所述的基于改进RANSAC的直线/圆弧检测算法,包括:
Ⅰ.首先利用霍夫直线/圆弧检测对亚像素边缘进行检测得到多个直线/圆弧模型;
Ⅱ.选取尺寸最大的直线/圆弧模型作为最优模型;
Ⅲ.求取ROI区域中的点集Q(xi,yi)到所拟合出模型的距离,其中i=1,2…n;
Ⅳ.判断所得到的距离是否满足阈值要求,不满足要求的点直接舍弃;
Ⅴ.将步骤Ⅳ中满足阈值要求的点保存在内点集合P(xi,yi)中,其中i=1,2…n-k;k为所舍弃的点的个数;
Ⅵ.对内点集P进行最小二乘直线/圆弧拟合,得到精确的直线/圆弧方程。
所述的尺寸测量是对大尺寸工件进行二维尺寸测量,包括直线与直线之间的距离以及圆弧的半径。对于直线间距的测量,假定其中一条直线的一般式为Ax+By+C=0,另外一条直线上的点集P=(xi,yi),其中i=0,1,2……,n;则直线间距测量的数学模型可表示为:
其中,A、B、C是直线参数。
在一种实施方式中,所述的光源为环形光源,设置为前向光源方式。
在一种实施方式中,所述的镜头的焦距符合式(1)
f=wD/V (1)
其中,w为CMOS芯片的靶面尺寸;V为视场尺寸;D为工作距离。
在一种实施方式中,所述的镜头为Kowa品牌的12mm镜头,镜头焦距为12±3mm。
在一种实施方式中,所述的高精度运动装置包括高精度三维运动装置和辅助装置,能够完成镜头在整个工作范围的运动,实现大尺寸工件的完整图像采集。在一种实施方式中,所述的相机通过转接板安装在三维运动装置的单轴驱动器上,在伺服电机的驱动下随单轴驱动器一起运动。
本发明还提供一种基于机器视觉的大尺寸工件自动测量方法,所述的方法包括获取大尺寸工件图像、图像预处理、图像数据转换建模、基于改进RANSAC的亚像素直线/圆弧检测、尺寸计算。
一种基于机器视觉的大尺寸工件自动测量方法中的图像预处理包括:(1)灰度化操作,将彩色图像转化为灰度图;(2)对步骤(1)得到的灰度图进行高斯滤波,得到边缘保留完整且对比度高的图像;(3)通过Canny边缘检测得到待测区域的粗边缘信息。
一种基于机器视觉的大尺寸工件自动测量方法的图像数据转换建模是通过坐标转换来获得所需尺寸的算法,包括:(a)通过光电编码器精确定位相机位置,采集到多幅工件图像,以第一张图像为基准图像;(b)利用式(2)和式(3)所示的坐标转换模型将其余图像的坐标转换为基准图像下的坐标;
i=i′+(in-i0)×xi (2)
j=j′+(jn-j0)×yj (3)
式中,每幅图像的位置索引为(i,j),基准图像的索引为(1,1),其中i、j是转换之后的坐标,i'、j'是待转换图像自身的像素坐标,in、jn是待转换图像的位置索引,i0、j0是基准图像的索引,xi、yj是图像在i和j方向上的像素个数;
所述的基于改进RANSAC的亚像素直线/圆弧检测算法,包括基于三次样条插值的亚像素边缘检测算法及基于改进RANSAC的直线/圆弧检测算法;
所述的基于三次样条插值的亚像素边缘检测算法包括:
①使用Canny边缘检测算子对所选的ROI区域进行边缘检测得到粗边缘;
②以粗边缘上某一点为中心沿垂直于边缘的左右方向或上下方向各取三个像素点并得到其灰度值;
③对步骤②中取得的六个点和中心点进行三次样条插值,得到边界点的光照强度函数,为六个三次多项式S(x);
④求解函数S(x)二阶导数的过零点,所述的过零点为亚像素级的边界点;然后对求得的亚像素边界点重新进行边缘拟合;
所述的基于改进RANSAC的直线/圆弧检测算法,包括:
Ⅰ.首先利用霍夫直线/圆弧检测对亚像素边缘进行检测得到多个直线/圆弧模型;
Ⅱ.选取尺寸最大的直线/圆弧模型作为最优模型;
Ⅲ.求取ROI区域中的点集Q(xi,yi)到所拟合出模型的距离,其中i=1,2…n;
Ⅳ.判断所得到的距离是否满足阈值要求,不满足要求的点直接舍弃;
Ⅴ.将步骤Ⅳ中满足阈值要求的点保存在内点集合P(xi,yi)中,其中i=1,2…n-k;k为所舍弃的点的个数;
Ⅵ.对内点集P进行最小二乘直线/圆弧拟合,得到精确的直线/圆弧方程;
所述的尺寸测量是对大尺寸工件进行二维尺寸测量,包括直线与直线之间的距离以及圆弧的半径;对于直线间距的测量,假定其中一条直线的一般式为Ax+By+C=0,另外一条直线上的点集P=(xi,yi),其中i=0,1,2……,n;则直线间距测量的数学模型可表示为:
其中,A、B、C是直线参数。
本发明有益效果是:
(1)本发明提供的系统能够针对大尺寸的工件进行测量;填补了机器视觉在大尺寸工件测量领域应用的空白;
(2)本发明中使用了灰度转换、高斯滤波、Canny边缘检测来对图像做预处理,有效的去除了噪声、无效点,并且还能对非均匀背景进行校正;解决了图像预处理中干扰点对检测结果的影响;
(3)本发明中使用图像数据转换建模替代传统的图像拼接方法,显著提高了生产效率,提高了检测的稳定性,在特征点不明显的情况下也不会影响测量结果;进而提高了测量系统的生产效率、系统的稳定型、系统测量准确率;
(4)本发明使用了改进RANSAC检测算法和亚像素边缘检测算法,克服了传统的基于RANSAC思想的检测方法效率低下的问题,提高了检测效率,缩短了整个处理周期;引入了亚像素边缘检测算法克服了常规边缘检测算子无法对单个像素进行细分,检测出的边缘于实际边缘存在偏差的问题,进而提高了整个系统的生产效率和检测精度。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1系统机械结构;
图2坐标转换模型;
图3基于改进RANSAC的直线检测算法流程图;
图4基于改进RANSAC的直线检测结果;
图5直线间距测量流程图;
图6系统算法流程图;
图7钢板工件外形图。
具体实施方式
实施例一:
本发明的一种基于机器视觉的大尺寸工件自动测量系统,如图1所示,包括工作台,待测钢板工件,传动装置,相机镜头和光源。
所述的相机镜头设置在传动装置下方,所述的光源设置在相机镜头的下方,所述的待测钢板工件放置在工作台上;工作台放置在光源的下方。
光源一般可设置为背光光源和前向光源;背向光源通常能够获得较高的对比度,前向光源可以获得待测工件的表面信息;本实施例中的光源设置为前向光源的方式。而在图像采集过程中,如果没有合适的光源将会导致所拍摄的图像明暗不一,这会对后续的图像处理带来不必要的麻烦。针对待测工件的特点,本申请采用的是环形光源。
本实施例中采用成本低廉、原理简单的CMOS工业摄像机进行图像的采集。由于工件尺寸较大,本实施例选取了是大恒公司生产的MER-1520-13U3C-L型CMOS面阵摄像机。像素个数为4608*3288,CMOS传感器面积为1/2.3英寸,摄像机接口为Mini USB3.0。
本事实力中的镜头的基本光学性能有焦距、分辨率和景深等,其中,焦距是确定镜头最主要的参数;焦距f应符合公式(1):
f=wD/V (1)
式中,w为CMOS芯片的靶面尺寸;V为视场尺寸;D为工作距离。本实施例中,待测工件的尺寸为长度范围100mm×100mm-2000mm×2000mm,视场尺寸定为130mm×100mm。通过计算得到焦距为11.781mm,本实施例中选用了Kowa品牌的12mm镜头,该镜头焦距为12±3mm。
本实施例中的机械装置,具体是:由于待测工件尺寸较大,本申请采用相机运动而工件不动的方法实现图像的采集。系统的总体机械结构如图1所示,本申请主要涉及高精度三维运动装置和辅助装置,以完成传感器在整个工作范围的运动,从而实现大尺寸工件的完整图像采集。相机通过转接板安装在单轴驱动器上,在伺服电机的驱动下随单轴驱动器一起运动,其运动距离是已知的。通过x、y和z三个方向上的运动,可以按顺序采集到大尺寸工件的多张图像。
本实施例中的大尺寸工件自动测量软件系统,包括图像预处理算法、图像数据转换建模、基于改进RANSAC的亚像素直线/圆弧检测算法和尺寸计算。
所述的图像预处理算法,具体是:为了更直观地获得所要识别的特征,需要对采集到的图像进行预处理操作。本申请的图像预处理算法主要包括三个,首先进行灰度化操作,将彩色图像转化为灰度图;然后进行高斯滤波,得到边缘保留完整且对比度高的图像;最后通过Canny边缘检测得到待测区域的粗边缘信息。
所述的图像数据转换建模,具体是:一般来说,对于大尺寸工件的尺寸计算,常用的方法是通过特征点匹配将所有图像拼接成一张大图再对其进行测量。但是,这种方法处理起来速度比较慢,影响生产效率。因此,本申请设计了一种通过坐标转换来获得所需尺寸的算法。
首先,通过光电编码器精确定位相机位置,采集大尺寸工件的多张图像,以第一张图像为基准图像,利用发明中公式(2)和公式(3)所示的坐标转换模型将其余图像的坐标转换为基准图像下的坐标;然后,利用本申请设计的算法对其进行亚像素边缘检测以及后续的尺寸测量。本申请中坐标转换的基本思路如下:
通过光电编码器精确定位相机位置,采集到六幅工件图像。如图2所示,图中的六个矩形代表采集到的六幅图像,并给它们一个位置索引(i,j)。本申请将第一张图像(索引为(1,1))作为基准图像,其坐标原点作为基准原点,则剩余图像的坐标可以通过以下公式转换为基准图像坐标系下的坐标:
i=i′+(in-i0)×xi (2)
j=j′+(jn-j0)×yj (3)
式中,i、j是转换之后的坐标,i'、j'是待转换图像自身的像素坐标,in、jn是待转换图像的位置索引,i0、j0是基准图像的索引,xi、yj是图像在i和j方向上的像素个数。
所述的基于改进RANSAC的亚像素直线/圆弧检测算法,主要包括基于三次样条插值的亚像素边缘检测算法以及基于改进RANSAC的直线/圆弧检测算法。
所述的基于三次样条插值的亚像素边缘检测算法,具体是:由于光照等因素的影响,一般的边缘检测算子,例如Sobel算子、Canny算子和Roberts算子等,这些边缘检测算子只能判断边缘的整像素位置,无法对单个像素再进行细分,因此,所检测出的边缘不一定是图像实际的边缘,不能满足更高的测量精度要求。针对这一问题,本申请采用基于三次样条插值的亚像素边缘检测算法,实现图像边缘亚像素坐标点的提取。该算法的具体实现步骤如下:
(1)首先用Canny边缘检测算子对所选的ROI区域进行边缘检测得到粗边缘。
(2)以粗边缘上某一点为中心沿垂直于边缘的方向左右(上下)各取三个像素点并得到其灰度值。
(3)对这七个点进行三次样条插值,得到六个三次多项式S(x)(边界点的光照强度函数)。
(4)求解函数S(x)二阶导数的过零点,即为我们所需要的亚像素级的边界点,然后对这些亚像素边界点重新进行边缘拟合。
所述的基于改进RANSAC的直线/圆弧检测算法,其算法流程图如图3所示。具体是:对于通用的直线/圆弧检测算法,例如霍夫直线检测、随机抽样一致算法(比如RANSAC)、最小二乘法等,这些算法的测量精度低,计算复杂度高,不适用于工业自动化检测。因此,本申请设计了一种基于改进RANSAC的直线/圆弧检测算法。RANSAC算法是根据一组包含异常数据的样本数据集,计算出数据的数学模型参数,得到有效样本数据的算法。该算法的优点是能鲁棒地估计模型参数,缺点是计算量大,运算周期长。针对这一缺点,本申请所设计算法的基本思路如下:
(1)首先利用霍夫直线/圆弧检测对亚像素边缘进行检测得到多个直线/圆弧模型。
(2)选取尺寸最大的直线/圆弧模型作为最优模型。
(3)求取ROI区域中的点集Q(xi,yi)(i=1,2…n)到所拟合出模型的距离。
(4)判断所得到的距离是否满足阈值要求,不满足要求的点直接舍弃。
(5)将(4)中满足阈值要求的点保存在内点集合P(xi,yi)(i=1,2…n-k)(其中k为所舍弃的点的个数)中。
(6)对内点集P进行最小二乘直线/圆弧拟合,得到精确的直线/圆弧方程。
如图4所示,为使用本申请设计的基于改进RANSAC的直线检测算法对亚像素边缘进行检测得到的直线,黑色区域为检测出的直线。
所述的尺寸测量,即对大尺寸工件进行二维尺寸测量,包括直线与直线之间的距离以及圆弧的半径。对于直线间距的测量,假定其中一条直线的一般式为Ax+By+C=0,另外一条直线上的点集P=(xi,yi)(i=0,1,2……,n),则直线间距测量的数学模型可表示为:
式中,A、B、C是直线参数,直线间距测量的流程图如图5所示。
实施效果
如图6所示,为本申请的算法流程图。按照该流程,对图7所示的钢板工件进行检测,便可以得到钢板的二维尺寸。在本申请的实验中,利用所设计的大尺寸工件尺寸测量算法测得的结果平均误差为0.05mm,对单个尺寸测量的重复精度在0.01mm以内。而且,算法稳定性好,鲁棒性强,在工业测量领域具有较好的应用前景。
本实施例相比于人工检测,拥有更高的检测速度和精度,并且拥有良好的重复性,避免了人工测量带来的偶然误差;同时,该系统的软件接口能够很好地与计算机进行连接,发展和改良空间巨大,可以有效地代替人工测量。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于机器视觉的大尺寸工件自动测量系统,其特征在于,包括大尺寸工件自动测量硬件系统和大尺寸工件自动测量软件系统;
所述大尺寸工件自动测量硬件系统,包括光源、相机和镜头以及机械装置;
所述大尺寸工件自动测量软件系统,包括图像预处理算法、图像数据转换建模、基于改进RANSAC的亚像素直线/圆弧检测算法和尺寸计算。
2.根据权利要求1所述的基于机器视觉的大尺寸工件自动测量系统,其特征是,所述的图像预处理算法包括:(1)灰度化操作,将彩色图像转化为灰度图;(2)对步骤(1)得到的灰度图进行高斯滤波,得到边缘保留完整且对比度高的图像;(3)通过Canny边缘检测得到待测区域的粗边缘信息。
3.根据权利要求1所述的基于机器视觉的大尺寸工件自动测量系统,其特征在于,
所述的图像数据转换建模是通过坐标转换来获得所需尺寸的算法,包括:
(a)通过光电编码器精确定位相机位置,采集到多幅工件图像,以第一张图像为基准图像;
(b)利用式(2)和式(3)所示的坐标转换模型将其余图像的坐标转换为基准图像下的坐标;
i=i′+(in-i0)×xi (2)
j=j′+(jn-j0)×yj (3)
式中,每幅图像的位置索引为(i,j),基准图像的索引为(1,1),其中i、j是转换之后的坐标,i'、j'是待转换图像自身的像素坐标,in、jn是待转换图像的位置索引,i0、j0是基准图像的索引,xi、yj是图像在i和j方向上的像素个数。
4.根据权利要求1所述的基于机器视觉的大尺寸工件自动测量系统,其特征在于,所述的基于改进RANSAC的亚像素直线/圆弧检测算法,包括基于三次样条插值的亚像素边缘检测算法及基于改进RANSAC的直线/圆弧检测算法。
5.根据权利要求4所述的基于机器视觉的大尺寸工件自动测量系统,其特征在于,
所述的基于三次样条插值的亚像素边缘检测算法包括:
①使用Canny边缘检测算子对所选的ROI区域进行边缘检测得到粗边缘;
②以粗边缘上某一点为中心沿垂直于边缘的左右方向或上下方向各取三个像素点并得到其灰度值;
③对步骤②中取得的六个点和中心点进行三次样条插值,得到边界点的光照强度函数,为六个三次多项式S(x);
④求解函数S(x)二阶导数的过零点,所述的过零点为亚像素级的边界点;然后对求得的亚像素边界点重新进行边缘拟合。
6.根据权利要求4所述的基于机器视觉的大尺寸工件自动测量系统,其特征在于,
所述的基于改进RANSAC的直线/圆弧检测算法,包括:
Ⅰ.首先利用霍夫直线/圆弧检测对亚像素边缘进行检测得到多个直线/圆弧模型;
Ⅱ.选取尺寸最大的直线/圆弧模型作为最优模型;
Ⅲ.求取ROI区域中的点集Q(xi,yi)到所拟合出模型的距离,其中i=1,2…n;
Ⅳ.判断所得到的距离是否满足阈值要求,不满足要求的点直接舍弃;
Ⅴ.将步骤Ⅳ中满足阈值要求的点保存在内点集合P(xi,yi)中,其中i=1,2…n-k;k为所舍弃的点的个数;
Ⅵ.对内点集P进行最小二乘直线/圆弧拟合,得到精确的直线/圆弧方程。
7.根据权利要求1所述的基于机器视觉的大尺寸工件自动测量系统,其特征在于,
所述的尺寸测量是对大尺寸工件进行二维尺寸测量,包括直线与直线之间的距离以及圆弧的半径;对于直线间距的测量,假定其中一条直线的一般式为Ax+By+C=0,另外一条直线上的点集P=(xi,yi),其中i=0,1,2……,n;则直线间距测量的数学模型可表示为:
其中,A、B、C是直线参数。
8.一种基于机器视觉的大尺寸工件自动测量方法,其特征在于,所述的方法包括获取大尺寸工件图像、图像预处理、图像数据转换建模、基于改进RANSAC的亚像素直线/圆弧检测、尺寸计算。
9.根据权利要求8所述的基于机器视觉的大尺寸工件自动测量方法,其特征在于,所述的图像预处理包括:(1)灰度化操作,将彩色图像转化为灰度图;(2)对步骤(1)得到的灰度图进行高斯滤波,得到边缘保留完整且对比度高的图像;(3)通过Canny边缘检测得到待测区域的粗边缘信息。
10.根据权利要求9所述的基于机器视觉的大尺寸工件自动测量方法,其特征在于,所述的图像数据转换建模是通过坐标转换来获得所需尺寸的算法,包括:(a)通过光电编码器精确定位相机位置,采集到多幅工件图像,以第一张图像为基准图像;
(b)利用式(2)和式(3)所示的坐标转换模型将其余图像的坐标转换为基准图像下的坐标;
i=i′+(in-i0)×xi (2)
j=j′+(jn-j0)×yj (3)
式中,每幅图像的位置索引为(i,j),基准图像的索引为(1,1),其中i、j是转换之后的坐标,i'、j'是待转换图像自身的像素坐标,in、jn是待转换图像的位置索引,i0、j0是基准图像的索引,xi、yj是图像在i和j方向上的像素个数;
所述的基于改进RANSAC的亚像素直线/圆弧检测算法,包括基于三次样条插值的亚像素边缘检测算法及基于改进RANSAC的直线/圆弧检测算法;
所述的基于三次样条插值的亚像素边缘检测算法包括:
①使用Canny边缘检测算子对所选的ROI区域进行边缘检测得到粗边缘;
②以粗边缘上某一点为中心沿垂直于边缘的左右方向或上下方向各取三个像素点并得到其灰度值;
③对步骤②中取得的六个点和中心点进行三次样条插值,得到边界点的光照强度函数,为六个三次多项式S(x);
④求解函数S(x)二阶导数的过零点,所述的过零点为亚像素级的边界点;然后对求得的亚像素边界点重新进行边缘拟合;
所述的基于改进RANSAC的直线/圆弧检测算法,包括:
Ⅰ.首先利用霍夫直线/圆弧检测对亚像素边缘进行检测得到多个直线/圆弧模型;
Ⅱ.选取尺寸最大的直线/圆弧模型作为最优模型;
Ⅲ.求取ROI区域中的点集Q(xi,yi)到所拟合出模型的距离,其中i=1,2…n;
Ⅳ.判断所得到的距离是否满足阈值要求,不满足要求的点直接舍弃;
Ⅴ.将步骤Ⅳ中满足阈值要求的点保存在内点集合P(xi,yi)中,其中i=1,2…n-k;k为所舍弃的点的个数;
Ⅵ.对内点集P进行最小二乘直线/圆弧拟合,得到精确的直线/圆弧方程;
所述的尺寸测量是对大尺寸工件进行二维尺寸测量,包括直线与直线之间的距离以及圆弧的半径;对于直线间距的测量,假定其中一条直线的一般式为Ax+By+C=0,另外一条直线上的点集P=(xi,yi),其中i=0,1,2……,n;则直线间距测量的数学模型可表示为:
其中,A、B、C是直线参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811541212.5A CN109612390B (zh) | 2018-12-17 | 2018-12-17 | 基于机器视觉的大尺寸工件自动测量系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811541212.5A CN109612390B (zh) | 2018-12-17 | 2018-12-17 | 基于机器视觉的大尺寸工件自动测量系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109612390A true CN109612390A (zh) | 2019-04-12 |
CN109612390B CN109612390B (zh) | 2019-11-08 |
Family
ID=66009455
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811541212.5A Active CN109612390B (zh) | 2018-12-17 | 2018-12-17 | 基于机器视觉的大尺寸工件自动测量系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109612390B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109631787A (zh) * | 2018-12-27 | 2019-04-16 | 长安大学 | 透射式靶标图像的光斑中心检测方法及桥梁挠度图像式检测装置 |
CN110288651A (zh) * | 2019-06-21 | 2019-09-27 | 湖南大学 | 大尺寸工件的公差视觉检测方法、装置和计算设备 |
CN110517276A (zh) * | 2019-07-12 | 2019-11-29 | 郑州磨料磨具磨削研究所有限公司 | 一种基于机器视觉的砂轮修整方法及装置 |
CN110675393A (zh) * | 2019-09-29 | 2020-01-10 | 上海呈彧智能科技有限公司 | 一种基于机器视觉的坯料规格检测方法 |
CN110879080A (zh) * | 2019-11-15 | 2020-03-13 | 武汉华中天经通视科技有限公司 | 一种高温锻件高精度智能测量仪和测量方法 |
CN111504188A (zh) * | 2020-04-28 | 2020-08-07 | 湖北文理学院 | 基于机器视觉的圆弧零件测量方法及装置 |
CN111862050A (zh) * | 2020-07-22 | 2020-10-30 | 无锡先导智能装备股份有限公司 | 一种物料检测系统、方法及设备 |
CN112418123A (zh) * | 2020-11-30 | 2021-02-26 | 西南交通大学 | 一种基于Hough变换的工程图图线及线型识别的方法 |
CN113432543A (zh) * | 2021-06-18 | 2021-09-24 | 浙江大学台州研究院 | 一种边缘零件直径尺寸测量方法 |
CN113510536A (zh) * | 2021-04-29 | 2021-10-19 | 厦门大学 | 一种加工中心的在机检测装置和方法 |
CN113643225A (zh) * | 2020-04-26 | 2021-11-12 | 北京配天技术有限公司 | 一种圆弧检测方法以及圆弧检测装置 |
CN113828948A (zh) * | 2021-11-23 | 2021-12-24 | 济南邦德激光股份有限公司 | 一种激光切割机的板材寻边方法、标定系统及寻边系统 |
CN114264236A (zh) * | 2021-12-22 | 2022-04-01 | 含光微纳科技(深圳)有限公司 | 基于机器视觉的工件几何尺寸自动化检测方法及系统 |
CN114383505A (zh) * | 2022-01-06 | 2022-04-22 | 江苏大学 | 一种用于短轴类零件尺寸的自动检测装置 |
CN114608444A (zh) * | 2020-10-15 | 2022-06-10 | 浙江大学台州研究院 | 一种基于评分系统的尺寸测量方法 |
CN114663427A (zh) * | 2022-04-25 | 2022-06-24 | 北京与子成科技有限公司 | 基于图像处理的锅炉零件尺寸检测方法 |
CN115235335A (zh) * | 2022-05-24 | 2022-10-25 | 武汉工程大学 | 一种高铁动车组走行部尺寸智能检测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1098268A2 (de) * | 1999-11-03 | 2001-05-09 | Institut für Neurosimulation und Bildtechnologien GmbH | Verfahren zur dreidimensionalen optischen Vermessung von Objektoberflächen |
CN101182990A (zh) * | 2007-11-30 | 2008-05-21 | 华南理工大学 | 一种基于机器视觉的大型在制工件几何测量系统 |
CN101349542A (zh) * | 2008-06-27 | 2009-01-21 | 东南大学 | 大尺寸零件的视觉测量装置 |
CN102722705A (zh) * | 2012-06-12 | 2012-10-10 | 武汉大学 | 一种基于ransac算法的多车道线检测方法 |
CN104197860A (zh) * | 2014-07-04 | 2014-12-10 | 丽水学院 | 大尺寸工件的三维表面形貌测量方法 |
CN104359403A (zh) * | 2014-11-21 | 2015-02-18 | 天津工业大学 | 基于亚像素边缘算法的平面零件尺寸测量方法 |
CN105678763A (zh) * | 2016-01-05 | 2016-06-15 | 吴晓军 | 一种稳定的直线/圆特征检测方法 |
JP2017049050A (ja) * | 2015-08-31 | 2017-03-09 | 株式会社ミツトヨ | 画像測定装置、その制御プログラム及び測定装置 |
CN108489394A (zh) * | 2018-04-17 | 2018-09-04 | 沈阳建筑大学 | 一种大尺寸薄板金属工件几何质量自动检测装置及方法 |
CN108844462A (zh) * | 2018-06-05 | 2018-11-20 | 广州视源电子科技股份有限公司 | 一种尺寸检测方法、装置、设备和系统 |
-
2018
- 2018-12-17 CN CN201811541212.5A patent/CN109612390B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1098268A2 (de) * | 1999-11-03 | 2001-05-09 | Institut für Neurosimulation und Bildtechnologien GmbH | Verfahren zur dreidimensionalen optischen Vermessung von Objektoberflächen |
CN101182990A (zh) * | 2007-11-30 | 2008-05-21 | 华南理工大学 | 一种基于机器视觉的大型在制工件几何测量系统 |
CN101349542A (zh) * | 2008-06-27 | 2009-01-21 | 东南大学 | 大尺寸零件的视觉测量装置 |
CN102722705A (zh) * | 2012-06-12 | 2012-10-10 | 武汉大学 | 一种基于ransac算法的多车道线检测方法 |
CN104197860A (zh) * | 2014-07-04 | 2014-12-10 | 丽水学院 | 大尺寸工件的三维表面形貌测量方法 |
CN104359403A (zh) * | 2014-11-21 | 2015-02-18 | 天津工业大学 | 基于亚像素边缘算法的平面零件尺寸测量方法 |
JP2017049050A (ja) * | 2015-08-31 | 2017-03-09 | 株式会社ミツトヨ | 画像測定装置、その制御プログラム及び測定装置 |
CN105678763A (zh) * | 2016-01-05 | 2016-06-15 | 吴晓军 | 一种稳定的直线/圆特征检测方法 |
CN108489394A (zh) * | 2018-04-17 | 2018-09-04 | 沈阳建筑大学 | 一种大尺寸薄板金属工件几何质量自动检测装置及方法 |
CN108844462A (zh) * | 2018-06-05 | 2018-11-20 | 广州视源电子科技股份有限公司 | 一种尺寸检测方法、装置、设备和系统 |
Non-Patent Citations (7)
Title |
---|
YAO ZHIYING: "Dynamic Detection System of Workpiece Based on Machine Vision", 《2010 INTERNATIONAL CONFERENCE ON INTELLIGENT COMPUTATION TECHNOLOGY AND AUTOMATION》 * |
伍济钢等: "薄片零件机器视觉图像亚像素边缘检测", 《中国机械工程》 * |
冯西等: "基于机器视觉的大尺寸工件自动测量系统", 《传感器与微系统》 * |
刘凌云等: "基于坐标变换的快速高精度图像拼接算法", 《自动化与仪表》 * |
李雪娇等: "基于机器视觉的大型零件尺寸测量技术研究", 《光电技术应用》 * |
秦煜等: "基于RANSAC的激光网格标记图像特征提取", 《计算机工程与科学》 * |
马文娟: "基于机器视觉的零件平面尺寸自动测量", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109631787A (zh) * | 2018-12-27 | 2019-04-16 | 长安大学 | 透射式靶标图像的光斑中心检测方法及桥梁挠度图像式检测装置 |
CN110288651A (zh) * | 2019-06-21 | 2019-09-27 | 湖南大学 | 大尺寸工件的公差视觉检测方法、装置和计算设备 |
CN110517276A (zh) * | 2019-07-12 | 2019-11-29 | 郑州磨料磨具磨削研究所有限公司 | 一种基于机器视觉的砂轮修整方法及装置 |
CN110675393A (zh) * | 2019-09-29 | 2020-01-10 | 上海呈彧智能科技有限公司 | 一种基于机器视觉的坯料规格检测方法 |
CN110675393B (zh) * | 2019-09-29 | 2024-03-26 | 上海呈彧智能科技有限公司 | 一种基于机器视觉的坯料规格检测方法 |
CN110879080A (zh) * | 2019-11-15 | 2020-03-13 | 武汉华中天经通视科技有限公司 | 一种高温锻件高精度智能测量仪和测量方法 |
CN113643225A (zh) * | 2020-04-26 | 2021-11-12 | 北京配天技术有限公司 | 一种圆弧检测方法以及圆弧检测装置 |
CN111504188A (zh) * | 2020-04-28 | 2020-08-07 | 湖北文理学院 | 基于机器视觉的圆弧零件测量方法及装置 |
CN111504188B (zh) * | 2020-04-28 | 2021-11-16 | 湖北文理学院 | 基于机器视觉的圆弧零件测量方法及装置 |
CN111862050A (zh) * | 2020-07-22 | 2020-10-30 | 无锡先导智能装备股份有限公司 | 一种物料检测系统、方法及设备 |
CN114608444A (zh) * | 2020-10-15 | 2022-06-10 | 浙江大学台州研究院 | 一种基于评分系统的尺寸测量方法 |
CN114608444B (zh) * | 2020-10-15 | 2024-01-26 | 浙江大学台州研究院 | 一种基于评分系统的尺寸测量方法 |
CN114705120B (zh) * | 2020-10-15 | 2023-06-23 | 浙江大学台州研究院 | 一种基于掩膜信息的尺寸测量方法及评分系统 |
CN114705120A (zh) * | 2020-10-15 | 2022-07-05 | 浙江大学台州研究院 | 一种基于掩膜信息的尺寸测量方法及评分系统 |
CN112418123A (zh) * | 2020-11-30 | 2021-02-26 | 西南交通大学 | 一种基于Hough变换的工程图图线及线型识别的方法 |
CN112418123B (zh) * | 2020-11-30 | 2021-08-03 | 西南交通大学 | 一种基于Hough变换的工程图图线及线型识别的方法 |
CN113510536B (zh) * | 2021-04-29 | 2022-07-29 | 厦门大学 | 一种加工中心的在机检测装置和方法 |
CN113510536A (zh) * | 2021-04-29 | 2021-10-19 | 厦门大学 | 一种加工中心的在机检测装置和方法 |
CN113432543A (zh) * | 2021-06-18 | 2021-09-24 | 浙江大学台州研究院 | 一种边缘零件直径尺寸测量方法 |
CN113828948A (zh) * | 2021-11-23 | 2021-12-24 | 济南邦德激光股份有限公司 | 一种激光切割机的板材寻边方法、标定系统及寻边系统 |
CN114264236A (zh) * | 2021-12-22 | 2022-04-01 | 含光微纳科技(深圳)有限公司 | 基于机器视觉的工件几何尺寸自动化检测方法及系统 |
CN114383505A (zh) * | 2022-01-06 | 2022-04-22 | 江苏大学 | 一种用于短轴类零件尺寸的自动检测装置 |
CN114383505B (zh) * | 2022-01-06 | 2024-06-07 | 江苏大学 | 一种用于短轴类零件尺寸的自动检测装置 |
CN114663427A (zh) * | 2022-04-25 | 2022-06-24 | 北京与子成科技有限公司 | 基于图像处理的锅炉零件尺寸检测方法 |
CN115235335A (zh) * | 2022-05-24 | 2022-10-25 | 武汉工程大学 | 一种高铁动车组走行部尺寸智能检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109612390B (zh) | 2019-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109612390B (zh) | 基于机器视觉的大尺寸工件自动测量系统 | |
CN110370286B (zh) | 基于工业机器人和单目相机的定轴运动刚体空间位置识别方法 | |
CN110672040B (zh) | 一种基于视觉的高精度旋转角度测量方法 | |
CN104266608B (zh) | 视觉传感器现场标定装置和标定方法 | |
CN105783723B (zh) | 基于机器视觉的精密模具表面加工精度检测装置及方法 | |
CN112629441B (zh) | 3d曲面玻璃轮廓扫描检测方法及系统 | |
CN107121093A (zh) | 一种基于主动视觉的齿轮测量装置及测量方法 | |
CN112614098B (zh) | 一种基于增强现实的毛坯定位与加工余量分析方法 | |
CN107392849B (zh) | 基于图像细分的靶标识别与定位方法 | |
CN111604598A (zh) | 一种机械臂进给式激光刻蚀系统的对刀方法 | |
CN103499297A (zh) | 一种基于ccd的高精度测量方法 | |
CN102853786A (zh) | 平整度检测装置和方法 | |
CN109211110A (zh) | 一种基于机器视觉的平板类工件尺寸检测系统及方法 | |
CN107358628B (zh) | 基于靶标的线阵图像处理方法 | |
CN102519400A (zh) | 基于机器视觉的大长径比轴类零件直线度误差检测方法 | |
CN110763136A (zh) | 一种高精度三维工件尺寸检测系统及方法 | |
CN111833392B (zh) | 标记点多角度扫描方法、系统及装置 | |
CN110608685A (zh) | 一种基于树莓派的物体尺寸快速测量方法 | |
JP7353757B2 (ja) | アーチファクトを測定するための方法 | |
Liu et al. | Measuring method for micro-diameter based on structured-light vision technology | |
CN113932710B (zh) | 一种复合式视觉刀具几何参数测量系统及方法 | |
CN108180825B (zh) | 一种基于线结构光的长方体形物体三维识别与定位方法 | |
CN105823430B (zh) | 成像视野小于柔性环形零件尺寸时的图像采集与拼接方法 | |
CN117537706A (zh) | 基于三轴坐标转换与单目视觉融合的钣金孔洞定位尺寸检测方法 | |
CN112747671A (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220517 Address after: 214000 No. 6, Huayun Road, Huazhuang street, Wuxi Economic Development Zone, Jiangsu Province Patentee after: Jiangsu Jumu Technology Co.,Ltd. Address before: No. 1800 Lihu Avenue, Wuxi City, Jiangsu Province Patentee before: Jiangnan University |
|
TR01 | Transfer of patent right |