CN109238084B - 一种微型圆孔测量的自动导引方法 - Google Patents

一种微型圆孔测量的自动导引方法 Download PDF

Info

Publication number
CN109238084B
CN109238084B CN201810988563.4A CN201810988563A CN109238084B CN 109238084 B CN109238084 B CN 109238084B CN 201810988563 A CN201810988563 A CN 201810988563A CN 109238084 B CN109238084 B CN 109238084B
Authority
CN
China
Prior art keywords
image
ellipse
round hole
probe
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.)
Active
Application number
CN201810988563.4A
Other languages
English (en)
Other versions
CN109238084A (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.)
Hefei Luyang Technology Innovation Group Co.,Ltd.
Original Assignee
Hefei 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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201810988563.4A priority Critical patent/CN109238084B/zh
Publication of CN109238084A publication Critical patent/CN109238084A/zh
Application granted granted Critical
Publication of CN109238084B publication Critical patent/CN109238084B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B5/00Measuring arrangements characterised by the use of mechanical techniques
    • G01B5/004Measuring arrangements characterised by the use of mechanical techniques for measuring coordinates of points
    • G01B5/008Measuring arrangements characterised by the use of mechanical techniques for measuring coordinates of points using coordinate measuring machines
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/14Measuring arrangements characterised by the use of optical techniques for measuring distance or clearance between spaced objects or spaced apertures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
    • G01B21/04Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness by measuring coordinates of points
    • G01B21/047Accessories, e.g. for positioning, for tool-setting, for measuring probes
    • 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/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
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image

Abstract

本发明公开了一种微型圆孔测量的自动导引方法,利用左相机和右相机针对探针和待测圆孔进行拍摄获取左源图和右源图,通过边缘提取生成图像A和图像B;分别对图像A及图像B进行椭圆检测,获取每张图像中探针及圆孔的椭圆轮廓。依据四个椭圆获取四个包络矩形,再依据四个包络矩形截取左源图及右源图中探针和圆孔的影像放大生成四张图像进行边缘提取,保留最外层边缘点得到四个椭圆心坐标,依此获得探针及待测圆孔中心在左源图和右源图中的坐标,最后利用标定的相机参数,获取探针与待测微型圆孔中心的三维距离,驱动三坐标测量机完成导引。本发明实现了探针向微型圆孔中心的自动导引,能极大地提高微型圆孔测量的测量效率。

Description

一种微型圆孔测量的自动导引方法
技术领域
本发明涉及图像处理领域,具体涉及一种适用于用三坐标测量机测量的情况下,利用接触式探头测量微型圆孔场合的自动导引方法。
背景技术
随着微纳米技术的发展,科学仪器朝向精密化发展,零部件向着微型化发展,其中圆孔在零部件中的应用十分广泛。但是微型圆孔的加工,尤其是孔径1mm以下的圆孔加工是十分困难的,例如油泵,油嘴及磨具等。因此,如何准确地检测微型圆孔的加工质量及相关参数是重要的技术问题。
近年来,科研工作者利用结构光、激光干涉及激光三角法等方式测量微型圆孔的相关参数,但对于高深宽比的圆孔的加工质量及参数,这些测量方法难以测量。所以,其中有一些测量方式会剖开圆孔,再对圆孔进行相关测量,但这些测量方式会对待测工件产生破坏性影响。
利用接触式测量探头对微型圆孔进行测量,不仅测量精度满足要求且可以测量不同深度的圆孔加工质量及相关参数。但是,利用接触式探头进行测量时,需要将探头顶部的探针置入圆孔之内,由于圆孔的孔径小,其孔径仅仅比探针直径稍大,肉眼难以完成导引操作。虽然工业上有利用相机对孔位及探针进行拍摄,人们通过拍摄的影像进行判断,导引探针进入圆孔的实例,但是此方法的测量效率低下、自动化程度不高。
发明内容
本发明是为避免上述现有技术的不足,提供一种适用于接触式探头测量微型圆孔的自动导引方法,以期在利用接触式探头测量高深宽比的圆孔时,能够快速、准确地到达待测圆孔的中心以开始自动测量。
本发明为解决技术问题采用如下技术方案:
本发明微型圆孔测量的自动导引方法,所述微型圆孔测量是指采用三坐标测量机,并利用接触式测量探头前端的探针探入待测圆孔中进行微孔测量;其特点是:在所述三坐标测量机的一侧基座上固定设置型号相同、且光轴平行的左相机和右相机,所述探针和待测圆孔均处在左相机和右相机的视场中,其中相机所拍摄图像均以像素为单位,以图像左上角顶点为原点建立坐标系,像素的横坐标与纵坐标分别是在其图像数组中所在的列数与所在行数;所述微型圆孔测量的自动导引方法是按如下步骤进行:
步骤1、利用左相机拍摄获取左源图,所述左源图中包含有探针和待测圆孔的完整影像;利用右相机拍摄获取右源图,所述右源图中包含有探针和待测圆孔的完整影像;
步骤2、根据设定的边缘阈值对所述左源图进行边缘提取并生成图像A,对所述右源图进行边缘提取并生成图像B;
步骤3、对所述图像A和图像B分别进行自动阈值的霍夫椭圆检测,针对所述图像A分别获得由所述探针的边缘轮廓生成的椭圆C1,以及由待测圆孔的边缘轮廓生成的椭圆C3;针对所述图像B分别获得由所述探针的边缘轮廓生成的椭圆C2,以及由待测圆孔的边缘轮廓生成的椭圆C4;且椭圆C1的离心率小于椭圆C3的离心率,椭圆C2的离心率小于椭圆C4的离心率;针对自动阈值的霍夫椭圆检测,若不满足条件一,则调整霍夫椭圆检测阈值直到从所述图像A和图像B检测出的椭圆满足条件一时继续执行步骤4;若是调整霍夫椭圆检测阈值达到最大值时仍不满足条件一,则退出导引;
条件一:椭圆C1和椭圆C3之间的相对位置与椭圆C2和椭圆C4之间的相对位置为相同;
步骤4、获取探针和待测圆孔的影像生成图像:
获取椭圆C1在所述图像A中的包络矩形R1,根据所述包络矩形R1各顶点在所述图像A中的坐标,在所述左源图中截取图像并按倍数N1进行图像放大,获得探针的影像生成图像P1;获取椭圆C2在所述图像B中的包络矩形R2,根据所述包络矩形R2各顶点在所述图像B中的坐标,在所述右源图中截取图像并按倍数N2进行图像放大,获得探针的影像生成图像P2;获取椭圆C3在所述图像A中的包络矩形R3,根据所述包络矩形R3各顶点在所述图像A中的坐标,在所述左源图中截取图像并按倍数N3进行图像放大,获得待测圆孔的影像生成图像P3;获取椭圆C4在所述图像B中的包络矩形R4,根据所述包络矩形R4各顶点在所述图像B中的坐标,在所述右源图中截取图像并按倍数N4进行图像放大,获得待测圆孔的影像生成图像P4;所述图像放大是以不产生影像形变的方式进行图像放大;
步骤5、针对所述图像P1、图像P2、图像P3和图像P4分别进行边缘提取,去除内部边缘点,只保留最外层的边缘点并拟合椭圆;针对图像P1、图像P2、图像P3和图像P4预设距离阈值,若图像Pi中存在一边缘点M到所述拟合椭圆上与所述边缘点M对应的切点的距离dij大于拟合椭圆相应的距离阈值Di,则剔除过于离散的边缘点M并重新拟合椭圆,直至所有边缘点均满足要求;由此对应获得图像Pi中拟合椭圆的椭圆心Ei的坐标,i=1,2,3,4,j=1,2,3……n,n表示所述图像Pi中的边缘点个数;
步骤6、按如下方式分别获得各源图中探针和待测圆孔的中心点坐标:
根据所述图像P1中椭圆心E1、包络矩形R1、图像A及倍数N1,获得左源图中探针的中心点坐标;根据所述图像P2中椭圆心E2、包络矩形R2、图像B及倍数N2,获得右源图中探针的中心点坐标;根据所述图像P3中椭圆心E3、包络矩形R3、图像A及倍数N3,获得左源图中待测圆孔的中心点坐标;根据所述图像P4中椭圆心E4、包络矩形R4、图像B及倍数N4,获得右源图中待测圆孔的中心点坐标;
利用通过标定得到的左相机和右相机的内部参数、外部参数和畸变参数,经过矩阵运算获得所述探针的中心点与待测圆孔的中心点的三维距离,依据所述三维距离驱动三坐标测量机完成导引。
本发明微型圆孔测量的自动导引方法的特点也在于,在所述步骤5中针对图像P1、图像P2、图像P3和图像P4分别进行边缘提取,按如下方式获得最外层的边缘点:针对边缘提取后的图像P1、图像P2、图像P3和图像P4,分别以图像上边界向图像下边界、图像下边界向图像上边界、图像左边界到图像右边界以及图像右边界到图像左边界四个方向分别遍历图像边缘点,且仅保留每个方向第一个遍历到的边缘点作为所述最外层的边缘点。
本发明微型圆孔测量的自动导引方法的特点也在于,所述步骤3中自动阈值的霍夫椭圆检测是按如下方法进行:首先预设最小圆心距离以防止在图像A与图像B中针对探针的轮廓及待测圆孔的轮廓检测出过多的椭圆;再以低霍夫椭圆检测阈值检测图像A,然后加大霍夫椭圆检测阈值直到图像A中只存在两个椭圆;以低霍夫椭圆检测阈值检测图像B,然后加大霍夫椭圆检测阈值直到图像B中只存在两个椭圆。
本发明微型圆孔测量的自动导引方法的特点也在于,在所述步骤3中,针对所述条件一是按如下方式进行判断:
若是:在图像A中,椭圆C1的圆心位于椭圆C3的圆心的下方;且在图像B中,椭圆C2的圆心位于椭圆C4的圆心的下方;或是:在图像A中,椭圆C1的圆心位于椭圆C3的圆心的上方,且在图像B中,椭圆C2的圆心位于椭圆C4的圆心的上方,均为满足条件一。
本发明微型圆孔测量的自动导引方法的特点也在于,在所述步骤5中,所述距离dij由式(1)计算得出:
Figure BDA0001780256040000031
其中,xij表示图像Pi中边缘点M在图像Pi中的横坐标,yij表示图像Pi中边缘点M在图像Pi中的纵坐标,fij表示边缘点M在图像Pi的拟合椭圆上的相应切点在图像Pi中的横坐标,gij表示边缘点M在图像Pi的拟合椭圆上的相应切点在图像Pi中的纵坐标。
本发明微型圆孔测量的自动导引方法的特点也在于,所述步骤4中所述图像放大倍数Ni由式(2)计算得出:
Figure BDA0001780256040000041
其中,Ni为图像Pi的放大倍数,wi为包络矩形Ri的宽度,hi为包络矩形Ri的高度,W为对应源图的宽度,H为对应源图的高度。
本发明微型圆孔测量的自动导引方法的特点也在于,按式(3)分别获得探针和待测圆孔在各源图中的中心点坐标:
Figure BDA0001780256040000042
Ui表示探针的中心点或待测圆孔的中心点在对应源图中的横坐标;
Vi表示探针的中心点或待测圆孔的中心点在对应源图中的纵坐标;
X1i表示包络矩形Ri的左上角顶点在所述对应源图中的横坐标,Y1i表示包络矩形Ri的左上角顶点在对应源图中的纵坐标,X2i表示椭圆心Ei在所述图像Pi中的横坐标,Y2i表示椭圆心Ei在所述图像Pi中的纵坐标。
本发明微型圆孔测量的自动导引方法的特点也在于,所述包络矩形的方向与所述左源图或右源图的方向保持一致,即所述包络矩形与对应的左源图或右源图的水平轴的夹角为0°。
与已有技术相比,本发明有益效果体现在:
1、本发明针对三坐标测量机、利用接触式探头测量圆孔的应用场合,自动导引探头进入待测圆孔从而实现自动化测量,将定位及导引整合为一体,极大地提高了圆孔测量实用性和测量效率;
2、本发明方法针对接触式探头的球形探针及待测圆孔区域分别进行提取,利用提取结果的相对位置关系进行提取准确与否的判断,能有效减少误导引;在物像不发生形变的情况下进行放大,能减小定位误差,有效提高了自动导引的准确性;
3、本发明方法在拟合接触式探头的球形探针及待测圆孔边界的椭圆心时,预先根据边缘点的位置对外部边缘点进行保留,对内部边缘点进行去除,再根据距离阈值剔除过于离散的坏点,提高了拟合椭圆心的准确性。
附图说明
图1为本发明方法流程图;
图2为应用在本发明方法中的测量系统示意图;
具体实施方式
参见图2,本实施例中微型圆孔测量是指采用三坐标测量机,并利用接触式测量探头前端的探针2探入待测圆孔3中进行微孔测量;三坐标测量机由基座、X轴6、Y轴7及Z轴8组成,在三坐标测量机的一侧基座上固定设置型号相同、且光轴平行的左相机4和右相机5,接触式测量探头1悬挂于三坐标测量机的Z轴8上,在接触式测量探头1的顶端安装有球形探针2,探针2和待测圆孔3均处在左相机4和右相机5的视场中,其中相机所拍摄图像均以像素为单位,以图像左上角顶点为原点建立坐标系,像素的横坐标与纵坐标分别是在其图像数组中所在的列数与所在行数。
参见图1,本实施例中微型圆孔测量的自动导引方法是按如下步骤进行:
步骤1、利用左相机4拍摄获取左源图,左源图中包含有探针2和待测圆孔3的完整影像;利用右相机5拍摄获取右源图,右源图中包含有探针2和待测圆孔3的完整影像。
步骤2、根据设定的边缘阈值对左源图进行边缘提取并生成图像A,对右源图进行边缘提取并生成图像B;边缘阈值满足以该阈值提取后的边缘不损失探针及待测微型圆孔的边缘信息,同时没有太多的噪声干扰。
步骤3、对图像A和图像B分别进行自动阈值的霍夫椭圆检测,针对图像A分别获得由探针2的边缘轮廓生成的椭圆C1,以及由待测圆孔3的边缘轮廓生成的椭圆C3;针对图像B分别获得由探针2的边缘轮廓生成的椭圆C2,以及由待测圆孔3的边缘轮廓生成的椭圆C4;且椭圆C1的离心率小于椭圆C3的离心率,椭圆C2的离心率小于椭圆C4的离心率。由于相机与待测圆孔3所在平面成一定角度,所以待测圆孔3在图像中呈现出椭圆形轮廓。而由于光照及图像插值的影响,所以探针2呈现出椭圆形轮廓而非圆形。但是探针2本身是球形,所以探针2的椭圆形轮廓的离心率比待测圆孔3的椭圆形轮廓的离心率小。
为了降低椭圆误检测的几率,针对自动阈值的霍夫椭圆检测,若不满足条件一,则调整霍夫椭圆检测阈值直到从图像A和图像B检测出的椭圆满足条件一时继续执行步骤4;若是调整霍夫椭圆检测阈值达到最大值时仍不满足条件一,则退出导引。
条件一:椭圆C1和椭圆C3之间的相对位置与椭圆C2和椭圆C4之间的相对位置为相同。
步骤4、获取探针2和待测圆孔3的影像生成图像:
获取椭圆C1在图像A中的包络矩形R1,根据包络矩形R1各顶点在图像A中的坐标,在左源图中截取图像并按倍数N1进行图像放大,获得探针2的影像生成图像P1;获取椭圆C2在图像B中的包络矩形R2,根据包络矩形R2各顶点在图像B中的坐标,在右源图中截取图像并按倍数N2进行图像放大,获得探针2的影像生成图像P2;获取椭圆C3在图像A中的包络矩形R3,根据包络矩形R3各顶点在图像A中的坐标,在左源图中截取图像并按倍数N3进行图像放大,获得待测圆孔3的影像生成图像P3;获取椭圆C4在图像B中的包络矩形R4,根据包络矩形R4各顶点在图像B中的坐标,在右源图中截取图像并按倍数N4进行图像放大,获得待测圆孔3的影像生成图像P4;包络矩形在包含探针2或待测圆孔3完整影像的基础上,其宽度与高度应尽可能的小;图像放大是以不产生影像形变的方式进行图像放大;包络矩形的方向与左源图或右源图的方向保持一致,即包络矩形与对应的左源图或右源图的水平轴的夹角为0°。
步骤5、针对图像P1、图像P2、图像P3和图像P4分别进行边缘提取,去除内部边缘点,只保留最外层的边缘点并拟合椭圆;针对图像P1、图像P2、图像P3和图像P4预设距离阈值,若图像Pi中存在一边缘点M到拟合椭圆上与边缘点M对应的切点的距离dij大于拟合椭圆相应的距离阈值Di,则剔除过于离散的边缘点M并重新拟合椭圆,直至所有边缘点均满足要求;由此对应获得图像Pi中拟合椭圆的椭圆心Ei的坐标,i=1,2,3,4,j=1,2,3……n,n表示图像Pi中的边缘点个数。由于三坐标测量机的检测台、待测工件除待测圆孔之外的表面边缘等影响,对图像的边缘检测会产生过于离散的边缘点,这些边缘点会降低拟合椭圆的准确度,剔除后有助于提升椭圆心坐标的精度。
步骤6、按如下方式分别获得各源图中探针2和待测圆孔3的中心点坐标:
根据图像P1中椭圆心E1、包络矩形R1、图像A及倍数N1,获得左源图中探针2的中心点坐标;根据图像P2中椭圆心E2、包络矩形R2、图像B及倍数N2,获得右源图中探针2的中心点坐标;根据图像P3中椭圆心E3、包络矩形R3、图像A及倍数N3,获得左源图中待测圆孔3的中心点坐标;根据图像P4中椭圆心E4、包络矩形R4、图像B及倍数N4,获得右源图中待测圆孔3的中心点坐标;
利用通过标定得到的左相机4和右相机5的内部参数、外部参数和畸变参数,经过矩阵运算获得探针2的中心点与待测圆孔3的中心点的三维距离,依据三维距离驱动三坐标测量机完成导引。
具体实施中,相应措施也包括:
在步骤5中针对图像P1、图像P2、图像P3和图像P4分别进行边缘提取,按如下方式获得最外层的边缘点:针对边缘提取后的图像P1、图像P2、图像P3和图像P4,分别以图像上边界向图像下边界、图像下边界向图像上边界、图像左边界到图像右边界以及图像右边界到图像左边界四个方向分别遍历图像边缘点,且仅保留每个方向第一个遍历到的边缘点作为最外层的边缘点。
在步骤3中自动阈值的霍夫椭圆检测是按如下方法进行:首先预设最小圆心距离以防止在图像A与图像B中针对探针2的轮廓及待测圆孔3的轮廓检测出过多的椭圆;再以低霍夫椭圆检测阈值检测图像A,然后加大霍夫椭圆检测阈值直到图像A中只存在两个椭圆;以低霍夫椭圆检测阈值检测图像B,然后加大霍夫椭圆检测阈值直到图像B中只存在两个椭圆。
在步骤3中,针对条件一是按如下方式进行判断:
若是:在图像A中,椭圆C1的圆心位于椭圆C3的圆心的下方;且在图像B中,椭圆C2的圆心位于椭圆C4的圆心的下方;或是:在图像A中,椭圆C1的圆心位于椭圆C3的圆心的上方,且在图像B中,椭圆C2的圆心位于椭圆C4的圆心的上方,均为满足条件一。
在步骤5中,距离dij由式(1)计算得出:
Figure BDA0001780256040000071
其中,xij表示图像Pi中边缘点M在图像Pi中的横坐标,yij表示图像Pi中边缘点M在图像Pi中的纵坐标,fij表示边缘点M在图像Pi的拟合椭圆上的相应切点在图像Pi中的横坐标,gij表示边缘点M在图像Pi的拟合椭圆上的相应切点在图像Pi中的纵坐标。
在步骤4中图像放大倍数Ni由式(2)计算得出:
Figure BDA0001780256040000072
其中,Ni为图像Pi的放大倍数,wi为包络矩形Ri的宽度,hi为包络矩形Ri的高度;W为对应源图的宽度,H为对应源图的高度。此方式能够以不产生物像形变且最大程度地放大探针及待测圆孔的物像,提高获取的探针中心坐标及待测圆孔中心坐标的准确度。
按式(3)分别获得探针2和待测圆孔3在各源图中的中心点坐标:
Figure BDA0001780256040000073
Ui表示探针2的中心点或待测圆孔3的中心点在对应源图中的横坐标;
Vi表示探针2的中心点或待测圆孔3的中心点在对应源图中的纵坐标,
X1i表示包络矩形Ri的左上角顶点在对应源图中的横坐标,Y1i表示包络矩形Ri的左上角顶点在对应源图中的纵坐标,X2i表示椭圆心Ei在图像Pi中的横坐标,Y2i表示椭圆心Ei在图像Pi中的纵坐标。
实验验证:
S1、利用两台型号相同且分辨率为2048×1536的工业相机同时拍摄,分别获取左源图和右源图;
S2、对左源图进行边缘提取生成图像A,对右源图进行边缘提取生成图像B;
S3、对图像A、图像B进行自动阈值的霍夫椭圆检测,一共获取四个椭圆,对应椭圆的相对位置满足条件一,继续执行;
S4、根据四个椭圆圆心的图像坐标,获取对应包络矩形,依据包络矩形坐标截取并放大左源图和右源图中的探针及圆孔影像,图像A中探针的包络矩形的左上角顶点坐标为(1456,816),包络矩形宽度52像素,高度为50像素,则图像放大倍数为30;图像A中待测圆孔的包络矩形的左上角顶点坐标为(1050,788),包络矩形宽度60像素,高度为54像素,则放大倍数为28;图像B中探针的包络矩形的左上角顶点坐标为(892,822),包络矩形宽度为40像素,高度为38像素,则放大倍数为40;图像B中待测圆孔的包络矩形的左上角顶点坐标为(536,790),包络矩形宽度为48像素,高度为32像素,则放大倍数为42;
S5、对图像P1、图像P2、图像P3、图像P4进行边缘提取,之后保留外部边缘点,去除内部边缘点,对探针图像P1、图像P2剔除坏点后拟合椭圆,对待测圆孔图像P3、图像P4剔除坏点后拟合椭圆。图像P1中椭圆心坐标为(614,750),则探针中心在左源图的坐标为(1476.5,841);图像P2中椭圆心坐标为(556,610),则探针中心在右源图的坐标为(1069.9,809.8);图像P3中椭圆心坐标为(742,718),则待测圆孔中心在左源图的坐标为(910.6,840.0);图像P4中椭圆心坐标为(906,806),则待测圆孔中心在右源图的坐标为(557.6,809.2);
S6、根据预先标定得到的两台工业相机内部参数、外部参数及畸变参数,探针中心的世界坐标为(24.0458,3.97679,192.132),待测微型圆孔中心的世界坐标为(2.12379,2.44826,209.639),其中探针直径为1.00mm,圆孔直径为1.30mm,根据相机坐标系及三坐标测量机运动轴之间的平移矩阵及旋转矩阵,X轴行进距离为13.5172mm,Y轴行进距离为-21.7002mm,Z轴行进距离为-11.6518mm,将接触式探头的探针导引进入圆孔从而完成自动测量。

Claims (6)

1.一种微型圆孔测量的自动导引方法,所述微型圆孔测量是指采用三坐标测量机,并利用接触式测量探头(1)前端的探针(2)探入待测圆孔(3)中进行测量;其特征是,在所述三坐标测量机的一侧基座上固定设置型号相同、且光轴平行的左相机(4)和右相机(5),所述探针(2)和待测圆孔(3)均处在左相机(4)和右相机(5)的视场中,其中相机所拍摄图像均以像素为单位,以图像左上角顶点为原点建立坐标系,像素的横坐标与纵坐标分别是在其图像数组中所在的列数与所在行数;所述微型圆孔测量的自动导引方法是按如下步骤进行:
步骤1、利用左相机(4)拍摄获取左源图,所述左源图中包含有探针(2)和待测圆孔(3)的完整影像;利用右相机(5)拍摄获取右源图,所述右源图中包含有探针(2)和待测圆孔(3)的完整影像;
步骤2、根据设定的边缘阈值对所述左源图进行边缘提取并生成图像A,对所述右源图进行边缘提取并生成图像B;
步骤3、对所述图像A和图像B分别进行自动阈值的霍夫椭圆检测,针对所述图像A分别获得由所述探针(2)的边缘轮廓生成的椭圆C1,以及由待测圆孔(3)的边缘轮廓生成的椭圆C3;针对所述图像B分别获得由所述探针(2)的边缘轮廓生成的椭圆C2,以及由待测圆孔(3)的边缘轮廓生成的椭圆C4;且椭圆C1的离心率小于椭圆C3的离心率,椭圆C2的离心率小于椭圆C4的离心率;
针对自动阈值的霍夫椭圆检测,若不满足条件一,则调整霍夫椭圆检测阈值直到从所述图像A和图像B检测出的椭圆满足条件一时继续执行步骤4;若是调整霍夫椭圆检测阈值达到最大值时仍不满足条件一,则退出导引;
条件一:椭圆C1和椭圆C3之间的相对位置与椭圆C2和椭圆C4之间的相对位置为相同;
步骤4、获取探针(2)和待测圆孔(3)的影像生成图像:
获取椭圆C1在所述图像A中的包络矩形R1,根据所述包络矩形R1各顶点在所述图像A中的坐标,在所述左源图中截取图像并按倍数N1进行图像放大,获得探针(2)的影像生成图像P1
获取椭圆C2在所述图像B中的包络矩形R2,根据所述包络矩形R2各顶点在所述图像B中的坐标,在所述右源图中截取图像并按倍数N2进行图像放大,获得探针(2)的影像生成图像P2
获取椭圆C3在所述图像A中的包络矩形R3,根据所述包络矩形R3各顶点在所述图像A中的坐标,在所述左源图中截取图像并按倍数N3进行图像放大,获得待测圆孔(3)的影像生成图像P3
获取椭圆C4在所述图像B中的包络矩形R4,根据所述包络矩形R4各顶点在所述图像B中的坐标,在所述右源图中截取图像并按倍数N4进行图像放大,获得待测圆孔(3)的影像生成图像P4
所述图像放大是以不产生影像形变的方式进行图像放大;
步骤5、针对所述图像P1、图像P2、图像P3和图像P4分别进行边缘提取,去除内部边缘点,只保留最外层的边缘点并拟合椭圆;针对图像P1、图像P2、图像P3和图像P4预设距离阈值,若图像Pi中存在一边缘点M到所述拟合椭圆上与所述边缘点M对应的切点的距离dij大于拟合椭圆相应的距离阈值Di,则剔除过于离散的边缘点M并重新拟合椭圆,直至所有边缘点均满足要求;由此对应获得图像Pi中拟合椭圆的椭圆心Ei的坐标,i=1,2,3,4,j=1,2,3……n,n表示所述图像Pi中的边缘点个数;
所述距离dij由式(1)计算得出:
Figure FDA0002278295610000021
其中,xij表示图像Pi中边缘点M在图像Pi中的横坐标,yij表示图像Pi中边缘点M在图像Pi中的纵坐标,fij表示边缘点M在图像Pi的拟合椭圆上的相应切点在图像Pi中的横坐标,gij表示边缘点M在图像Pi的拟合椭圆上的相应切点在图像Pi中的纵坐标;
步骤6、按如下方式分别获得各源图中探针(2)和待测圆孔(3)的中心点坐标:
根据所述图像P1中椭圆心E1、包络矩形R1、图像A及倍数N1,获得左源图中探针(2)的中心点坐标;根据所述图像P2中椭圆心E2、包络矩形R2、图像B及倍数N2,获得右源图中探针(2)的中心点坐标;根据所述图像P3中椭圆心E3、包络矩形R3、图像A及倍数N3,获得左源图中待测圆孔(3)的中心点坐标;根据所述图像P4中椭圆心E4、包络矩形R4、图像B及倍数N4,获得右源图中待测圆孔(3)的中心点坐标;
利用通过标定得到的左相机(4)和右相机(5)的内部参数、外部参数和畸变参数,经过矩阵运算获得所述探针(2)的中心点与待测圆孔(3)的中心点的三维距离,依据所述三维距离驱动三坐标测量机完成导引。
2.根据权利要求1所述的微型圆孔测量的自动导引方法,其特征是,在所述步骤5中针对图像P1、图像P2、图像P3和图像P4分别进行边缘提取,并按如下方式获得最外层的边缘点:针对边缘提取后的图像P1、图像P2、图像P3和图像P4,分别以图像上边界向图像下边界、图像下边界向图像上边界、图像左边界到图像右边界以及图像右边界到图像左边界四个方向分别遍历图像边缘点,且仅保留每个方向第一个遍历到的边缘点作为所述最外层的边缘点。
3.根据权利要求1所述的微型圆孔测量的自动导引方法,其特征是,所述步骤3中自动阈值的霍夫椭圆检测是按如下方法进行:首先预设最小圆心距离以防止在图像A与图像B中针对探针(2)的轮廓及待测圆孔(3)的轮廓检测出过多的椭圆;再以低霍夫椭圆检测阈值检测图像A,然后加大霍夫椭圆检测阈值直到图像A中只存在两个椭圆;以低霍夫椭圆检测阈值检测图像B,然后加大霍夫椭圆检测阈值直到图像B中只存在两个椭圆。
4.根据权利要求1所述的微型圆孔测量的自动导引方法,其特征是,在所述步骤3中,针对所述条件一是按如下方式进行判断:若是:在图像A中,椭圆C1的圆心位于椭圆C3的圆心的下方;且在图像B中,椭圆C2的圆心位于椭圆C4的圆心的下方;或是:在图像A中,椭圆C1的圆心位于椭圆C3的圆心的上方,且在图像B中,椭圆C2的圆心位于椭圆C4的圆心的上方,均为满足条件一。
5.根据权利要求1所述的微型圆孔测量的自动导引方法,其特征是,所述步骤4中所述图像放大倍数Ni由式(2)计算得出:
Figure FDA0002278295610000031
其中,Ni为图像Pi的放大倍数,wi为包络矩形Ri的宽度,hi为包络矩形Ri的高度,W为对应源图的宽度,H为对应源图的高度。
6.根据权利要求5所述的微型圆孔测量的自动导引方法,其特征是,按式(3)分别获得探针(2)和待测圆孔(3)在各源图中的中心点坐标:
Figure FDA0002278295610000032
Ui表示探针(2)的中心点或待测圆孔(3)的中心点在对应源图中的横坐标;
Vi表示探针(2)的中心点或待测圆孔(3)的中心点在对应源图中的纵坐标;
X1i表示包络矩形Ri的左上角顶点在所述对应源图中的横坐标,Y1i表示包络矩形Ri的左上角顶点在对应源图中的纵坐标,X2i表示椭圆心Ei在所述图像Pi中的横坐标,Y2i表示椭圆心Ei在所述图像Pi中的纵坐标。
CN201810988563.4A 2018-08-28 2018-08-28 一种微型圆孔测量的自动导引方法 Active CN109238084B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810988563.4A CN109238084B (zh) 2018-08-28 2018-08-28 一种微型圆孔测量的自动导引方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810988563.4A CN109238084B (zh) 2018-08-28 2018-08-28 一种微型圆孔测量的自动导引方法

Publications (2)

Publication Number Publication Date
CN109238084A CN109238084A (zh) 2019-01-18
CN109238084B true CN109238084B (zh) 2020-04-14

Family

ID=65068744

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810988563.4A Active CN109238084B (zh) 2018-08-28 2018-08-28 一种微型圆孔测量的自动导引方法

Country Status (1)

Country Link
CN (1) CN109238084B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110360973B (zh) * 2019-08-28 2021-02-05 合肥工业大学 一种面向微型工件测量的自动引导方法
CN111030750B (zh) * 2019-10-09 2021-05-07 长飞光纤光缆股份有限公司 一种多模光纤dmd测试设备的探针配准方法及系统
CN112330599B (zh) * 2020-10-15 2024-02-20 浙江大学台州研究院 一种尺寸测量评分装置及调节方法和评分方法
CN112529869B (zh) * 2020-12-11 2023-07-21 中国航空工业集团公司金城南京机电液压工程研究中心 一种阀套节流方孔检测方法
CN113237425B (zh) * 2021-05-12 2022-07-12 合肥工业大学 单摄像头辅助坐标测量机测量微小工件的自动引导方法
CN114858056A (zh) * 2022-05-19 2022-08-05 华中科技大学 一种基于深孔测量装置自动大量测量小深孔孔深的方法
CN117470079B (zh) * 2023-12-26 2024-03-19 四川万圣通实业有限公司 一种管材表面的孔参数检测装置及孔参数检测方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001330428A (ja) * 2000-05-23 2001-11-30 Natl Inst Of Advanced Industrial Science & Technology Meti 3次元測定機の測定誤差評価方法及び3次元測定機用ゲージ
GB0508217D0 (en) * 2005-04-25 2005-06-01 Renishaw Plc Method for scanning the surface of a workpiece
JP2009047554A (ja) * 2007-08-20 2009-03-05 Tokai Kiyouhan Kk 三次元座標測定機
EP2698596A1 (en) * 2012-08-16 2014-02-19 Hexagon Technology Center GmbH Method and system for determining spatial coordinates with a mobile coordinate measuring machine

Also Published As

Publication number Publication date
CN109238084A (zh) 2019-01-18

Similar Documents

Publication Publication Date Title
CN109238084B (zh) 一种微型圆孔测量的自动导引方法
CN110118528B (zh) 一种基于棋盘靶标的线结构光标定方法
CN109612390B (zh) 基于机器视觉的大尺寸工件自动测量系统
CN103714535B (zh) 双目视觉测量系统中摄像机参数在线调整方法
CN107121093A (zh) 一种基于主动视觉的齿轮测量装置及测量方法
CN105783711B (zh) 三维扫描仪校正系统及其校正方法
CN110779797B (zh) 一种测量金属拉伸试验过程中塑性应变比的方法
CN109931869A (zh) 基于激光扫描成像的物料体积高精度检测方法
CN109751964B (zh) 一种高精度非接触式管径测量方法及装置
CN108562233A (zh) 利用二次曲线不变量的轴类零件直径尺寸在线测量方法
CN111127562B (zh) 一种单目面阵相机的标定方法及自动标定系统
CN109506629B (zh) 一种水下核燃料组件检测装置旋转中心标定的方法
CN112902869B (zh) 轨廓测量系统激光平面调整方法及装置
Lin et al. Real-time image-based defect inspection system of internal thread for nut
CN110657750A (zh) 一种用于刀具刃口钝化的检测系统和方法
CN114963981A (zh) 一种基于单目视觉的筒状零件对接非接触式测量方法
CN105928484A (zh) 基于双目视觉的电梯导轨自动测量系统
Iovenitti et al. Three-dimensional measurement using a single image
Fan et al. Dimensional inspecting system of shaft parts based on machine vision
CN111595560B (zh) 线激光器平面度的视觉检测方法
CN116091488B (zh) 一种发动机摇摆试验的位移测试方法及位移测试系统
CN104251687B (zh) 一种基于镜面图像处理的零件表面平整度检测方法
CN116817799B (zh) 用于超高层建筑物的垂直度测量系统
CN113983951B (zh) 三维目标的测量方法、装置、影像仪及存储介质
CN117329977B (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: 20220120

Address after: 230001 floor 6, block B, blue diamond Shangjie, No. 335, Suixi Road, Bozhou road street, Luyang District, Hefei City, Anhui Province

Patentee after: Hefei Luyang Technology Innovation Group Co.,Ltd.

Address before: Tunxi road in Baohe District of Hefei city of Anhui Province, No. 193 230009

Patentee before: Hefei University of Technology

TR01 Transfer of patent right