CN110440726B - 双轴平移台与结构光扫描仪的栓接结合面共面度测量方法 - Google Patents
双轴平移台与结构光扫描仪的栓接结合面共面度测量方法 Download PDFInfo
- Publication number
- CN110440726B CN110440726B CN201910819433.2A CN201910819433A CN110440726B CN 110440726 B CN110440726 B CN 110440726B CN 201910819433 A CN201910819433 A CN 201910819433A CN 110440726 B CN110440726 B CN 110440726B
- Authority
- CN
- China
- Prior art keywords
- translation table
- axis
- structured light
- scanner
- target
- 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
Links
- 238000013519 translation Methods 0.000 title claims abstract description 128
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000005259 measurement Methods 0.000 claims abstract description 41
- 239000000919 ceramic Substances 0.000 claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000012216 screening Methods 0.000 claims abstract description 3
- 230000033001 locomotion Effects 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 5
- 230000003287 optical effect Effects 0.000 claims description 5
- 239000004831 Hot glue Substances 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000011900 installation process Methods 0.000 claims description 2
- 238000012795 verification Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 6
- 230000007547 defect Effects 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 229920000049 Carbon (fiber) Polymers 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 239000004917 carbon fiber Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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
- G01B11/30—Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明双轴平移台和结构光扫描仪的栓接结合面共面度测量方法属于视觉测量领域,涉及一种基于双轴平移台与结构光扫描仪的螺栓连接结合面共面度测量方法。该方法先分别操作平移台的两个轴,每次使扫描仪沿平移台的单一轴平动。通过拍摄标准陶瓷球靶标,采集靶球球心坐标与扫描仪平移距离并与靶标检定数据进行比对,分别求解平移台两个轴的标定参数。然后换下靶标安装待测工件,操作平移台平动扫描仪以拍摄工件的螺栓连接结合面,通过图像处理算法提取特征区域形貌点云,根据平动距离进行点云拼接,筛去粗差点,最终计算共面度。该方法具有设备简单、易于操作、自动化程度高、精度较高的特点,有效扩展了结构光扫描仪的测量应用范围。
Description
技术领域
本发明属于视觉测量领域,涉及一种基于双轴平移台与结构光扫描仪的螺栓连接结合面共面度测量方法。
背景技术
螺栓连接是机械工程中一种非常普遍的连接方法,具有结构简单、性能可靠、方便拆卸、适应性强等特点。在机械设计时为了确保连接强度并获取较好的力学性能,通常两个零部件的结合面处设计有多个螺栓组成的螺栓组。而所有螺栓的结合面的共面度对连接处的刚度和模态特性有着极大的影响。而在航空航天等高端制造装配领域当中,考虑到设备的极端工况和特殊指标,对螺栓组的结合面共面度的测量就有极高的研究意义和现实需求。当前对于表面几何形貌的测量,主要有接触式测量与非接触式测量两种方法。接触式测量,最典型的应用就是三坐标测量机。该方法精度极高,但缺点是设备庞大昂贵,效率较慢,采集点云数据较少,不适合于大型复杂曲面。这一方法通常不适合于加工装配过程中的在线实时测量。而作为非接触式测量领域中的代表,结构光测量方法就能较好地解决上述问题,并广泛地运用于工业的实时在线测量当中。然而结构光测量也存在单次拍摄视野过小的问题,因此也需要多位姿拍摄数据拼接的方法来测量大尺寸零部件的表面几何特征。而采用高运动精度的双轴平移台,通过机械式拼接的方法,就能以较低的成本解决这一问题。此外,为了拾取结合面所在的区域,也需要通过图像处理的办法选择感兴趣区域(ROI),以提高测量效率,提升测量自动化程度。
对于平移台与结构光扫描仪的联合测量,最常见的做法是结构光扫描仪不动,平移台携带被测物体移动,例如,山东大学的张鹏于2017年发表的《基于快速扫描相移算法的三维形貌测试仪开发》。由于扫描仪景深有限,工件尺寸各异,这样每次测量时扫描仪的位姿就不得不改变以适应工件形状。这种方法明显不适于尺寸重量较大的螺栓连接部件,而且需要经常对设备进行标定。或者结构光扫描仪平直安装于单轴平移台上,默认成像平面与运动轴平行。这种方法无法保证高精度测量。而对于螺栓连接结合面ROI的分割提取,目前的主流商业软件中尚无专门的算法,需要使用者手动选取或者自行开发。显然,基于结构光扫描仪的螺栓连接结合面共面度测量方法的研究是很有必要的。
发明内容
本发明为克服现有技术的缺陷,发明了双轴平移台和结构光扫描仪的栓接结合面共面度测量方法,该方法通过两个平移台的运动模块以及标准陶瓷球完成了结构光扫描仪的位置标定,并且通过图像处理算法提取螺栓连接结合面的特征区域几何形貌,最终完成共面度的计算。该方法具有设备简单、易于操作、自动化程度高、精度较高等特点,有效扩展了结构光扫描仪的测量应用范围,可以应用于大尺寸部件螺栓连接结合面共面度的测量,是一种具有广泛应用前景的测量方法。
本发明采用的技术方案是一种双轴平移台和结构光扫描仪的栓接结合面共面度测量方法。该方法先分别操作平移台的两个轴,每次使扫描仪沿平移台的单一轴平动;通过拍摄标准陶瓷球靶标,采集靶球球心坐标与扫描仪平移距离并与靶标检定数据进行比对,分别求解平移台两个轴的标定参数;然后换下靶标安装待测工件,操作平移台平动扫描仪以拍摄待测工件的螺栓连接结合面,通过图像处理算法提取特征区域形貌点云,根据平动距离进行点云拼接,筛去粗差点,最后计算共面度。测量方法的具体步骤如下:
第一步、安置固定结构光扫描仪、双轴平移台以及支架。
将X轴平移台1与Y轴平移台2组成双轴平移台,并使双轴平移台固定,要求Y轴平移台2应尽量关于X轴平移台1对称垂直放置,以确保双轴平移台运行稳定并减轻振动。之后将结构光扫描仪3通过夹具与Y轴平移台2刚性连接。同样应确保结构光扫描仪3关于Y轴平移台2的载物台对称放置,使结构光扫描仪3随载物台移动时稳定平动而不发生振动。然后,安放固定支架4。固定支架4应不与X轴平移台1或Y轴平移台2接触,且不妨碍后者的运动。最后,通电预热,并检查设备。
第二步、标定平移台的X轴。
将标准陶瓷球靶标5沿X轴平移台1通过热熔胶固定于升降架6上。再松开升降架6,沿支架4上的轨道调整升降架6高度,并通过拧紧螺丝将升降架6与支架4固定。要求确保标准陶瓷球靶标5位于结构光扫描仪3的可测量范围内,调节结构光扫描仪的参数,最终保证结构光扫描仪成像清晰。标准陶瓷球靶标5上共有6个半径为R的哑光陶瓷球。
操控X轴平移台,将结构光扫描仪移动到适当位置,拍摄第1号靶球的表面点云,并记录当前结构光扫描仪在此平移台上的位置CI1,其中,角标中罗马字母I表示组号,数字1表示球号。将深度图的信息转换为点云的坐标,通过最小二乘拟合,得出位置CI1下靶球在扫描仪测量坐标系的球心坐标(XI1,YI1,ZI1)。球心坐标求解具体步骤为:
基于球面方程,可得出:
将球面点云坐标(x1,y1,z1)、(x2,y2,z2)……(xn,yn,zn)带入其中,得到最小二乘方程组:
最终求出球心坐标(XI1,YI1,ZI1)与半径RI1的最优解。半径R作为一个检验值判断(XI1,YI1,ZI1)是否有误。多次拍摄,参照拉依达准则,过滤粗差点,得到平均值。
之后重复上述操作5次,操控X轴平移台,分别移动结构光扫描仪到其他位置,拍摄标准陶瓷靶标上余下5个靶球,求出不同位置下球心相对于测量坐标系的坐标(XI2,YI2,ZI2)……(XI6,YI6,ZI6),并分别记录每次结构光扫描仪在此平移台上移动到的位置CI2……CI6。
对于高定位精度、高运动精度的单轴平移台,其运动过程可视为平移,旋转可以忽略不计。X轴平移台的运动机构只沿着一个轴L1移动,Y轴平移台的运动机构只沿着一个轴L2移动。由于结构光扫描仪内部的结构以及安装过程中装配工艺、夹具制造问题,轴L1与轴L2无法确保垂直;轴L1必然与扫描仪测量坐标系OMXMYMZM中的轴XM、YM、ZM存在夹角α1、θ1、β1,而轴L2必然与扫描仪测量坐标系OMXMYMZM中的轴XM、YM、ZM存在夹角α2、θ2、β2,这6个夹角余弦值也就是需求解的标定参数。
而另一方面,由于初等几何学的知识,
cos2α1+cos2θ1+cos2β1=1,cos2α2+cos2θ2+cos2β2=1(几何约束)。 (5)
将点系平移变换后的球心距和标定参数的几何约束,带入点系的距离约束中,可得:
(XIi-XIj)·cosα1+(YIi-YIj)·cosθ1+(ZIi-ZIj)·cosβ1=((XIi-XIj)2+(YIi-YIj)2+(ZIi-ZIj)2+(CIi-CIj)2-Lij 2)/(2(CIj-CIi))。 (6)
为了方便表述,设有函数φ(PIi,PIj,CIi,CIj,Lij)即:
φ(PIi,PIj,CIi,CIj,Lij)=((XIi-XIj)2+(YIi-YIj)2+(ZIi-ZIj)2+(CIi-CIj)2-Lij 2)/(2(CIj-CIi)) (7)
这样就把复杂的三元二次方程转化为较为简单的三元一次方程,降低计算复杂程度。每两个点组成一个方程,一共有n(n-1)/2个方程。将不同位置下的球心坐标带入其中,得到最小二乘方程组如下:
由此解出X轴的标定参数cosα1、cosθ1、cosβ1。
第三步、标定平移台的Y轴。
将标准陶瓷球靶标5卸下后,沿Y轴平移台2方向通过热熔胶固定于升降架6上。再松开升降架6,沿支架4上的轨道调整升降架6高度,并通过拧紧螺丝将升降架6与支架4固定。要求确保标准陶瓷球靶标5位于结构光扫描仪3的可测量范围内。然后操作Y轴平移台2,分别测量结构光扫描仪3于Y轴平移台2位置CII1……CII6时,靶球分别在其测量坐标系下的球心坐标(XII1,YII1,ZII1)……(XII6,YII6,ZII6)。参考第二步中的公式与方法,解出Y轴标定参数cosα2、cosθ2、cosβ2。至此,测量系统的6个标定参数求解完毕。
第四步、拍摄工件螺栓连接结合面,提取特征区域。
卸下标准陶瓷球靶标5,将待测工件7通过螺栓连接的方式安装于升降架6上;松开升降架6,调整其高度后拧紧螺丝将升降架6与支架4固定;确保待测工件7的待测特征区域位于结构光扫描仪3的测量范围内。要求待测区域面与铅锤方向的倾斜角度在10度以内。然后操控X轴平移台1与Y轴平移台2,平移结构光扫描仪3至合适位置,拍摄相应区域。记录每次拍摄的深度图Image1、Image2……Imagem(m≥2),记录结构光扫描仪3在X轴平移台1与Y轴平移台2上的位置坐标(AI1,BII1)、(AI2,BII2)……(AIm,BIIm)。
对于深度图Imagei(1≤i≤m),先将含有点云信息的灰度图进行二值化得到图像BWi。采用canny算子求BWi的边缘得到BWCi,再将BWCi中的孔洞进行填充得到BWCFi。设特征区域的估算面积为Si/mm2,扫描仪在XY方向上的分辨率为dmm/pixel,那么将BWCFi中像素面积大于1.5·Si/d2小于Si/(4d2)的区域过滤即可得到初始特征区域PROIi;对BWi进行分水岭算法,得到图像WSi,再将WSi中像素面积小于10000的区域过滤得到WSFi;采用大小为30*30像素的矩形结构元素对WSFi进行膨胀运算,得到孔的特征区域WSFDi;PROIi与WSFDi相减,最终得到结合面的特征区域ROIi;根据特征区域,将深度图中灰度值还原为特征几何形貌的点云矩阵PCi。对所有图片均以此操作。
第五步、图像拼接,求解共面度。
根据位置坐标与标定参数进行图像拼接。以获取PC1所在的(AI1,BII1)为拼接的基准位置。PCi的拍摄位置为(AIi,BIIi),将PCi向(AI1,BII1)位置转化。设PCi内某点坐标(xr,yr,zr)i,则转化后该点坐标为:
将(AI2,BII2)到(AIm,BIIm)位置下拍摄的点云矩阵PC2……PCm向位置(AI1,BII1)转化,最终得到拼接后的点云矩阵PC。
多个平面的共面度,可视为这些平面最小包络面的间距。应先拟合出其最小二乘平面,再将包络面上下平移求出最小包络面。设拟合出的最小二乘平面应满足方程:
式中a、b、d为待求系数,ε为残差。那么将点云矩阵PC坐标(x1,y1,z1)、(x2,y2,z2)……(xn,yn,zn)带入其中,可得最小而成矩阵:
最终,求出平面系数a、b、d。将每个点的坐标带入式(10),求得每点的残差ε1、ε2……εn。求出ε的平均值X与标注差σ为:
本发明的有益效果是采用基于结构光扫描仪的螺栓连接结合面共面度测量方法,该方法通过两个平移台的运动模块以及标准陶瓷球完成了结构光扫描仪的位置标定,并且通过图像处理算法提取螺栓连接结合面的特征区域几何形貌,最终完成共面度的计算。设备结构简单,操作流程容易,自动化程度较高,而且适应性强。与传统方法相比,该方法不需要重复标定,降低了综合成本,提高了测量精度,是一种具有广泛应用前景的标定方法。
附图说明
图1是测试系统的标定示意图,图2是待测工件测量示意图。其中,1-X轴平移台,2-Y轴平移台,3-结构光扫描仪,4-支架,5-标准陶瓷球靶标,6-升降架,7-待测工件。
图3是坐标系平移变换示意图。图中L1、L2表示单轴平移台的运动方向,OMXMYMZM表示测量坐标系的初始位置;OMIXMIYMIZMI与OMIIXMIIYMIIZMII分别表示结构光扫描仪沿X轴平移台和和轴平移台平移后的测量坐标系。
图4是图像处理流程示意图。
图5是单双轴平移台与结构光扫描仪的螺栓连接结合面共面度测量的流程图。
图6是最终的拼接效果示意图。
具体实施方式
下面结合附图和技术方案详细说明本发明的具体实施。
本实施例选用的结构光扫描仪3为LMI公司的3506型号测量仪,X轴平移台1与Y轴平移台2均为PI公司M-531平移台。选用的标准陶瓷球靶标5由6个标准半径5mm的哑光陶瓷球和长度250mm的碳纤维板组成。支架4与升降架6均由铝型材组装而成。
附图5是单双轴平移台与结构光扫描仪的螺栓连接结合面共面度测量的流程图。具体步骤如下:
第一步、安置固定结构光扫描仪、双轴平移台以及待测物支架。
将X轴平移台1与Y轴平移台2组成双轴平移台并使双轴平移台固定,之后将结构光扫描仪3通过夹具与Y轴平移台2刚性连接。同样应确保结构光扫描仪3关于Y轴平移台的载物台对称放置。然后安放固定支架4,固定支架4应不与X轴平移台1或Y轴平移台2接触,且不妨碍后者的运动,如图1、2所示。最后,通电预热并检查设备。
第二步、标定平移台的X轴。
先将标准陶瓷球靶标5沿X轴平移台固定于升降架6上,再调整升降架6的高度,并固定。然后操作X轴平移台1,分别测量结构光扫描仪于X轴平移台位置CI1……CI6时,靶球分别在其测量坐标系下的球心坐标(XI1,YI1,ZI1)……(XI6,YI6,ZI6),图3是坐标系平移变换示意图。所测数据如表1所示。
表1
经三坐标机标定的部分球心距如表2所示:(单位:mm)
表2
参考式(8),将6个位置数据、6个球心坐标以及经三坐标机标定的球心距分别代入其中。每两个点可构成一个方程,这样可得最多15个方程。这里仅用其中6个,在确保精度的前提下提高效率。
最终求得,cosα1=0.2762,cosθ1=0.9534,cosβ1=0.1213。从而完成了X轴平移台1的标定。
第三步、标定平移台的Y轴。
将标准陶瓷球靶标5,沿Y轴平移台2方向固定于升降架6上,再将升降架6调整高度,并固定。要求确保标准陶瓷球靶标5位于结构光扫描仪3的可测量范围内。然后操作Y轴平移台2,分别测量结构光扫描仪3于Y轴平移台位置CII1……CII6时,靶球分别在其测量坐标系下的球心坐标(XII1,YII1,ZII1)……(XII6,YII6,ZII6)。图3是坐标系平移变换示意图,球心坐标如表3所示。
表3
参考第二步中的公式与方法,解出Y轴标定参数:
cosα2=0.9546、cosθ2=0.2618、cosβ2=0.1421。
第四步、拍摄待测工件7螺栓连接结合面,提取特征区域。
卸下标准陶瓷球靶标5,将待测工件7安装于升降架6上,调整高度,确保工件的待测特征区域位于结构光扫描仪3的测量范围内。然后操控X轴平移台1与Y轴平移台2,平移结构光扫描仪3至合适位置,拍摄相应区域。记录每次拍摄的深度图Image1、Image2、Image3、Image4,记录结构光扫描仪在X轴平移台与Y轴平移台上的位置坐标。图4是图像处理流程示意图,得到的数据如表4所示。
表4
将四张深度图Image1、Image2、Image3、Image4经过图4所示的图像处理流程,选择特征区域点云。根据位置坐标(AI1,BII1)、(AI2,BII2)、(AI3,BII3)、(AI4,BII4)与标定参数cosα1、cosθ1、cosβ1、cosα2、cosθ2、cosβ2,将四个点云进行拼接,得到拼接后的点云如图6所示。
而残差ε的标准差为:δ=5.937×10-3。根据p=6δ,最终求得共面度p=35.62μm。
本发明是通过双轴平移台与结构光扫描仪的螺栓连接结合面共面度测量方法,该方法通过两个平移台的运动模块以及标准陶瓷球完成了结构光扫描仪的位置标定,并且通过图像处理算法提取螺栓连接结合面的特征区域几何形貌,最终完成共面度的计算。该方法具有设备简单、易于操作、自动化程度高、精度度测量较高等特点,有效扩展了结构光扫描仪的测量应用范围,该方法可以应用于复杂结构工件螺栓组的结合面共面度测量,是一种具有广泛应用前景的标定方法。
Claims (1)
1.一种双轴平移台和结构光扫描仪的栓接结合面共面度测量方法;其特征是,该方法先分别操作平移台的两个轴,每次使扫描仪沿平移台的单一轴平动;通过拍摄标准陶瓷球靶标,采集靶球球心坐标与扫描仪平移距离并与靶标检定数据进行比对,分别求解平移台两个轴的标定参数;然后换下靶标安装待测工件,操作平移台平动扫描仪以拍摄工件的螺栓连接结合面,通过图像处理算法提取特征区域形貌点云,根据平动距离进行点云拼接,筛去粗差点,最终计算共面度;测量方法的具体步骤如下:
第一步、安置固定结构光扫描仪、双轴平移台以及待测物支架;
将X轴平移台(1)与Y轴平移台(2)组成双轴平移台并使双轴平移台固定,要求Y轴平移台(2)应尽量关于X轴平移台(1)对称垂直放置以确保双轴平移台运行稳定并减轻振动;之后将结构光扫描仪(3)通过夹具与Y轴平移台(2)刚性连接;同样应确保结构光扫描仪(3)关于Y轴平移台(2)的载物台对称放置,使结构光扫描仪(3)随载物台移动时稳定平动而不发生振动,安放固定支架(4);然后,通电预热并检查设备;
第二步、标定平移台的X轴;
将标准陶瓷球靶标(5)沿X轴平移台(1)通过热熔胶固定于升降架(6)上;再松开升降架(6),沿支架(4)上的轨道调整升降架(6)高度,并通过拧紧螺丝将升降架(6)与支架(4)固定;要求确保标准陶瓷球靶标(5)位于结构光扫描仪(3)的可测量范围内,调节结构光扫描仪的参数,最终保证结构光扫描仪成像清晰;标准陶瓷球靶标(5)上共有6个半径为R的哑光陶瓷球;
操控X轴平移台(1),将结构光扫描仪移动到适当位置,拍摄第一个靶球的表面点云,并记录当前结构光扫描仪在此平移台上的位置CI1,其中,角标中罗马字母I表示组号,数字1表示球号;将深度图的信息转换为点云的坐标,通过最小二乘拟合,得出位置CI1下靶球在扫描仪测量坐标系的球心坐标(XI1,YI1,ZI1);球心坐标求解步骤为:
基于球面方程,得出:
将球面点云坐标(x1,y1,z1)、(x2,y2,z2)……(xn,yn,zn)带入其中,得到最小二乘方程组:
最终,求出球心坐标(XI1,YI1,ZI1)与半径RI1的最优解;半径R作为一个检验值判断(XI1,YI1,ZI1)是否有误;多次拍摄,参照拉依达准则,过滤粗差点,得到平均值;
之后重复上述操作5次,操控X轴平移台,分别移动结构光扫描仪到其他位置拍摄标准陶瓷靶标上余下5个靶球,求出不同位置下球心相对于测量坐标系的坐标(XI2,YI2,ZI2)……(XI6,YI6,ZI6),并分别记录每次结构光扫描仪在此平移台上移动到的位置CI2……CI6;
X轴平移台的运动机构只沿着一个轴L1移动,Y轴平移台的运动机构只沿着一个轴L2移动;由于结构光扫描仪内部的结构以及安装过程中装配工艺、夹具制造问题,轴L1与轴L2无法确保垂直;轴L1必然与扫描仪测量坐标系OMXMYMZM中的轴XM、YM、ZM存在夹角α1、θ1、β1,而轴L2必然与扫描仪测量坐标系OMXMYMZM中的轴XM、YM、ZM存在夹角α2、θ2、β2,这6个夹角余弦值也就是需求解的标定参数;
由初等几何学的知识,几何约束为:
cos2α1+cos2θ1+cos2β1=1,cos2α2+cos2θ2+cos2β2=1 (5)
将点系平移变换后的球心距和标定参数的几何约束,带入点系的距离约束中,得到:
(XIi-XIj)·cosα1+(YIi-YIj)·cosθ1+(ZIi-ZIj)·cosβ1=((XIi-XIj)2+(YIi-YIj)2+(ZIi-ZIj)2+(CIi-CIj)2-Lij 2)/(2(CIj-CIi)) (6)
为了方便表述,设有函数φ(PIi,PIj,CIi,CIj,Lij)即:
把复杂的三元二次方程转化为较为简单的三元一次方程,降低计算复杂程度;每两个点组成一个方程,一共有15个方程;将不同位置下的球心坐标带入其中,得到最小二乘方程组如下:
由此,解出X轴的标定参数cosα1、cosθ1、cosβ1;
第三步、标定平移台的Y轴;
将标准陶瓷球靶标(5)沿Y轴平移台(2)方向通过热熔胶固定于升降架(6)上,再松开升降架(6),沿支架(4)上的轨道调整升降架(6)高度,并通过拧紧螺丝将升降架(6)与支架(4)固定;要求确保标准陶瓷球靶标(5)位于结构光扫描仪(3)的可测量范围内;然后操作Y轴平移台(2),分别测量结构光扫描仪(3)于Y轴平移台(2)位置CII1……CII6时,靶球分别在其测量坐标系下的球心坐标(XII1,YII1,ZII1)……(XII6,YII6,ZII6);参考第二步中的公式与方法,解出Y轴标定参数cosα2、cosθ2、cosβ2;至此,测量系统的6个标定参数求解完毕;
第四步、拍摄工件螺栓连接结合面,提取特征区域;
卸下标准陶瓷球靶标(5),将待测工件(7)通过螺栓连接的方式安装于升降架(6)上,松开升降架(6),调整其高度后拧紧螺丝将升降架(6)与支架(4)固定;确保待测工件(7)的待测特征区域位于结构光扫描仪(3)的测量范围内;要求待测区域面与铅锤方向的倾斜角度在10度以内;然后操控X轴平移台(1)与Y轴平移台(2),平移结构光扫描仪(3)至合适位置,拍摄相应区域;记录每次拍摄的深度图Image1、Image2……Imagem(m≥2),记录结构光扫描仪(3)在X轴平移台(1)与Y轴平移台(2)上的位置坐标(AI1,BII1)、(AI2,BII2)……(AIm,BIIm);
对于深度图Imagei(1≤i≤m),先将含有点云信息的灰度图进行二值化得到图像BWi;采用canny算子求BWi的边缘得到BWCi,再将BWCi中的孔洞进行填充得到BWCFi;设特征区域的估算面积为Si/mm2,扫描仪在XY方向上的分辨率为dmm/pixel,那么将BWCFi中像素面积大于1.5·Si/d2小于Si/(4d2)的区域过滤即可得到初始特征区域PROIi;对BWi进行分水岭算法,得到图像WSi,再将WSi中像素面积小于10000的区域过滤得到WSFi,采用大小为30*30像素的矩形结构元素对WSFi进行膨胀运算,得到孔的特征区域WSFDi;PROIi与WSFDi相减,最终得到结合面的特征区域ROIi;根据特征区域,将深度图中灰度值还原为特征几何形貌的点云矩阵PCi;对所有图片均以此操作;
第五步、图像拼接,求解共面度;
根据位置坐标与标定参数进行图像拼接;以获取PC1所在的(AI1,BII1)为拼接的基准位置;PCi的拍摄位置为(AIi,BIIi),将PCi向(AI1,BII1)位置转化;设PCi内某点坐标(xr,yr,zr)i,则转化后该点坐标为:
将(AI2,BII2)到(AIm,BIIm)位置下拍摄的点云矩阵PC2……PCm向位置(AI1,BII1)转化,最终得到拼接后的点云矩阵PC;
多个平面的共面度,可视为这些平面最小包络面的间距;应先拟合出其最小二乘平面,再将包络面上下平移求出最小包络面;设拟合出的最小二乘平面应满足方程:
式中,a、b、d为待求系数,ε为残差;将点云矩阵PC坐标(x1,y1,z1)、(x2,y2,z2)……(xn,yn,zn)带入其中,得到最小而成矩阵:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910819433.2A CN110440726B (zh) | 2019-08-31 | 2019-08-31 | 双轴平移台与结构光扫描仪的栓接结合面共面度测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910819433.2A CN110440726B (zh) | 2019-08-31 | 2019-08-31 | 双轴平移台与结构光扫描仪的栓接结合面共面度测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110440726A CN110440726A (zh) | 2019-11-12 |
CN110440726B true CN110440726B (zh) | 2020-08-14 |
Family
ID=68438661
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910819433.2A Active CN110440726B (zh) | 2019-08-31 | 2019-08-31 | 双轴平移台与结构光扫描仪的栓接结合面共面度测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110440726B (zh) |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090056424A1 (en) * | 2007-08-29 | 2009-03-05 | Cornell Research Foundation, Inc. | Microscope Rheometer for Measuring Shear and Compression Properties of Biological Samples |
JP5570902B2 (ja) * | 2010-07-28 | 2014-08-13 | 本田技研工業株式会社 | 車両ピッチ角の推定装置 |
CN102654387B (zh) * | 2012-05-25 | 2014-07-02 | 南京理工大学 | 一种基于空间曲面约束的工业机器人在线标定装置及其实现方法 |
EP2847542B1 (de) * | 2012-08-07 | 2016-03-23 | Carl Zeiss Industrielle Messtechnik GmbH | Vorrichtung mit interner optischer referenz |
CN103954241B (zh) * | 2014-04-22 | 2017-01-25 | 华南理工大学 | 一种基于结构光的ic引脚共面度测量系统及其测量方法 |
CN104374338B (zh) * | 2014-09-28 | 2017-04-26 | 北京航空航天大学 | 一种基于固定相机和单靶标的单轴旋转角的视觉测量方法 |
CN105252341B (zh) * | 2015-09-02 | 2017-06-23 | 大连理工大学 | 五轴数控机床动态误差视觉测量方法 |
CN105423917B (zh) * | 2015-12-01 | 2018-05-29 | 中国科学院西安光学精密机械研究所 | 位置敏感探测器定位误差的标定方法 |
CN106354094B (zh) * | 2016-09-14 | 2018-10-16 | 大连理工大学 | 基于空间标准球的机床随动激光扫描坐标标定方法 |
CN106354092B (zh) * | 2016-09-24 | 2018-08-07 | 大连理工大学 | 一种随动与轮廓误差自适应实时补偿方法 |
CN206803969U (zh) * | 2017-06-06 | 2017-12-26 | 富加宜连接器(东莞)有限公司 | 一种点激光共面度测试装置 |
-
2019
- 2019-08-31 CN CN201910819433.2A patent/CN110440726B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110440726A (zh) | 2019-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Burner et al. | Videogrammetric model deformation measurement technique | |
CN111486802B (zh) | 基于自适应距离加权的旋转轴标定方法 | |
CN110146038B (zh) | 筒形件装配转角的分布式单目相机激光测量装置及方法 | |
CN109029299B (zh) | 舱段销孔对接转角的双相机测量装置及测量方法 | |
CN109357631B (zh) | 一种基于激光位移传感器的测量系统中心标定方法 | |
CN111076880B (zh) | 一种考虑相机姿态变化的大跨桥梁多点挠度测量方法 | |
CN111982019B (zh) | 基于线结构光传感器的叶片截面轮廓高精度检测方法 | |
CN111536902A (zh) | 一种基于双棋盘格的振镜扫描系统标定方法 | |
CN110926364B (zh) | 基于线结构光的叶片检测方法 | |
CN108562233B (zh) | 利用二次曲线不变量的轴类零件直径尺寸在线测量方法 | |
CN110926365B (zh) | 一种基于线结构光检测物标定方法 | |
CN106892133B (zh) | 负载摄像装置的飞行器云台稳定性测试方法 | |
CN111366908A (zh) | 一种激光雷达转台及其测量装置和测量方法 | |
CN109447971A (zh) | 一种摆线齿轮型线误差精确提取及关键形位精度快速分析方法 | |
CN113155047A (zh) | 长距离孔距测量装置、方法、存储介质、设备及轨道车辆 | |
CN111024047A (zh) | 基于正交双目视觉的六自由度位姿测量装置及方法 | |
CN107202553B (zh) | 全视角扫描测量系统及其靶标测量方法 | |
CN113702384A (zh) | 一种回转构件表面缺陷检测装置、检测方法及标定方法 | |
Deng et al. | Active 3-D thermography based on feature-free registration of thermogram sequence and 3-D shape via a single thermal camera | |
CN112697074B (zh) | 动态待测物角度测量仪及测量方法 | |
CN110440726B (zh) | 双轴平移台与结构光扫描仪的栓接结合面共面度测量方法 | |
CN113916128A (zh) | 一种基于光笔式视觉测量系统的提高精度的方法 | |
Wang et al. | Robust vision-based method for wing deflection angle measurement with defocus images | |
CN102519671B (zh) | 一种用于测量陀螺仪静平衡的基于双目视觉的空间位姿测量装置 | |
Jing et al. | The technologies of close-range photogrammetry and application in manufacture |
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 |