CN114894090A - 一种基于线激光的回转内腔尺寸检测方法 - Google Patents
一种基于线激光的回转内腔尺寸检测方法 Download PDFInfo
- Publication number
- CN114894090A CN114894090A CN202210501827.5A CN202210501827A CN114894090A CN 114894090 A CN114894090 A CN 114894090A CN 202210501827 A CN202210501827 A CN 202210501827A CN 114894090 A CN114894090 A CN 114894090A
- Authority
- CN
- China
- Prior art keywords
- point
- image
- coordinate
- points
- camera
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 23
- 238000000605 extraction Methods 0.000 claims abstract description 8
- 230000009466 transformation Effects 0.000 claims description 22
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000012216 screening Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000003287 optical effect Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 238000013519 translation Methods 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 18
- 238000011156 evaluation Methods 0.000 abstract description 2
- 230000005540 biological transmission Effects 0.000 abstract 1
- 238000000691 measurement method Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009776 industrial production Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 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/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
-
- 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/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于线激光的回转内腔尺寸检测方法,测量工具对回转内腔的尺寸进行测量;测量工具包括测量平台及安装于平台表面的传送带、线激光发射器以及摄像机;其中待测物放置在传送带表面,穿过线激光发发射器的发射区域,同时摄像机拍摄通过过程中的线激光轮廓图片。经过光条中心点提取算法获得图像坐标系的二维坐标,拟合计算出线激光图像角点的二维坐标,将角点图像坐标转换到空间坐标,从而计算距离实现对回转内腔尺寸的测量。本发明能够取代传统的三坐标测量仪和游标卡尺,对于各种结构匹配更灵活,检测状态评估更加精确。
Description
技术领域
本发明涉及机器视觉测量、精密测量,利用单目线结构光测量系统以及利用机器视觉测量技术领域,具体涉及一种基于线激光的回转内腔尺寸检测方法。
背景技术
如今,回转内腔广泛应用于医学用品行业以及航空航天用品等行业。长期以来,对回转内腔尺寸的检测主要依赖于游标卡尺两点测量以及三坐标测量机床等,但游标卡尺两点测量存在回转内腔直径磨损不易取点、取点不精确的问题,且无法实现在线动态测量,三坐标测量机床效率低,只能应用于精度要求相对较低,检测速度要求不高的工业生产,不能满足部分工业生产中快速检测的目的。
现有技术公开号为:CN102601681B公开了一种回转腔体的立式在线装卸及精密测量平台及测量方法,由易受冲击振动的装卸子平台和精密的测量子平台两部分组成,其中的测量子平台是通过地基隔振措施,保证了基准面精密测量环境与传送带复杂扰动环境的完全分离并没有具体涉及到测量方法的改进。
目前,随着机器视觉测量技术发展与日臻成熟,以线激光、结构光为代表的三维视觉测量技术广泛应用于螺距测量、航空发动机叶片检测、焊缝位姿识别研究以及钢轨磨耗量研究等方面。
因此,针对这类问题,将激光扫描用于检测回转内腔尺寸的成为可能。
发明内容
1.所要解决的技术问题:
针对上述技术问题,本发明提供一种基于线激光的回转内腔尺寸检测方法,采用线激光的三维视觉测量方法,其能够有效克服现有基于游标卡尺测量辅助检测的缺点。
2.技术方案:
一种基于线激光的回转内腔尺寸检测方法,其特征在于:用测量工具对回转内腔的尺寸进行测量;所述测量工具包括测量平台及安装于平台表面的传送带、线激光发射器、以及摄像机;所述待测物放置在传送带表面,穿过线激光发发射器的发射区域,同时摄像机拍摄通过过程中的激光图片,经过光条中心点提取算法获得图像坐标系的二维坐标,拟合计算出线激光图像角点的二维坐标,将角点图像坐标转换到空间坐标,从而计算距离;所述经过光条中心点提取算法获得图像坐标系的二维坐标,拟合计算出线激光图像角点的二维坐标,将角点图像坐标转换到空间坐标,从而计算距离具体包括以下步骤:
步骤一:获得回转内腔的线激光轮廓图片,所述回转内腔的线激光轮廓图片包括一条椭圆线条与一条直线线条;使用光条中心点提取方法获得轮廓形状图片的二维图像位置信息,并且用坐标点集表示;
步骤二:筛选出符合条件的的二维坐标点,作为所需拟合椭圆与直线的轮廓坐标点;
步骤三:进行直线与椭圆的拟合,求出两个交点的二维图像坐标;
步骤四:将步骤三生成的交点坐标转换到空间坐标;
步骤五:计算两个交点之间的距离。
进一步地,步骤一中:所述轮廓二维坐标点集表示为:[(X1,Y1),(X2,Y2),(X3,Y3)........(Xn,Yn)],其中n为整数,(Xn,Yn)表示第n个散点对应的坐标点。
进一步地,步骤二具体包括:
S21:在步骤一生成的轮廓二维坐标点集中找出第r个点为图像最低点,其坐标表示为(Xmin,Ymin),;所述图像最低点为轮廓形状图片中y坐标数值最小的点;计算每个散点和最低点连线的斜率,其计算公式如下式(1):
则所有散点与最低点连线的斜率集合表示为:k=[k1 … kn] (2);
在(1)(2)式中,r、m均为小于等于n的整数;
忽略由于散点集中产生的极大值,在散点与最低点连线的斜率集中,判断出两个极大值的点,则该两个点即为相邻的直线与椭圆之间的拟合点的分界点;
S22:直线拟合点筛选;从图像最低点(Xmin,Ymin)的左半边、右半边分别筛选拟合点;过图像最低点(Xmin,Ymin)做垂直与X轴的直线,将图像分割为左半边与右半边;
当从图像最低点左半边计算时,此时计算包括图像最低点,斜率k取到极大值时对应的点为第i个点,i∈(1,r),则左半边直线拟合点的序号为第1个点到第i*α个点,其中α为拟合系数,且0.95<α<1;
当从图像最低点右半边计算时,斜率取到极大值时对应的点为第j个点,j∈(r,n),则右半边直线拟合点的序号为第j*β个点到第n个点,其中β为拟合系数,且1<β<1.05;
S23:椭圆拟合点筛选;将步骤S22生成的序号为第i*β到j*α个对应的散点作为椭圆拟合点。
进一步地,步骤三具体包括以下步骤:
S31:将步骤S22筛选出的直线拟合点,采用matlab中的polyfit指令进行表述,其表述的形式为直线方程y=kx+b;
S32:对步骤S23筛选出的椭圆拟合点采用最小二乘法拟合为如下的椭圆方程进行表述:
f(X,Y)=AX2+BXY+CY2+DX+EY+F (3);
其中椭圆方程进行表述包括以下步骤:
S321:将平面内某点(X0,Y0)到椭圆方程所代表的曲线的距离表示为f(X0,Y0),对步骤S23筛选出的椭圆拟合点进行最小二乘处理,即优化目标为式(4):
S322:根据最小值来确定椭圆方程的系数A、B、C、D、E、F;
根据极值原理,对系数A、B、C、D、E、F求偏导,从而求得系数A,B,C,D,E,F的值;
S323:利用计算出的系数A,B,C,D,E,F的值计算出椭圆图像的五个参数:所述五个参数包括未知参数(θ,XC,YC)和形状参数(A,B);其中(XC,YC)为椭圆的几何中心;A、B分别为该椭圆的长半轴、短半轴;θ为椭圆的长短轴倾角;具体的计算过程如下:
几何中心:
长半轴短半轴:
长短轴倾角:
根据求出的五个参数,使用matlab函数ellipse即可画出椭圆效果图;
S33:根据步骤S31、S32拟合成的直线方程与椭圆方程,得到椭圆与直线交点分别为(X1’,Y1’),(X2’,Y2’)。
进一步地,步骤四具体包括以下步骤:
S41:将拟合算法中的图像坐标系转换成图像的像素坐标系;
将交点的坐标在X轴坐标不变的情况下,Y轴坐标减去图像的宽度h,即为:
X1’p=X1’
X2’p=X2’
Y1’P=|h-Y1’|
Y2’P=|h-Y2’|
则两个交点的图像坐标变成量纲为像素的像素坐标(X1’,Y1’P),(X2’,Y2’P);
S42:分别对两个交点的坐标进行三维坐标转换;每个交点的变换均包括以下步骤:
S421:设交点p的像素坐标为(xp,yp),不考虑镜头畸变,则其对应的图像坐标为(xn,yn),像素坐标系和图像坐标系原点均在图片左上角,则该两个坐标的变换关系为:
(9)式中,dx、dy分别表示单位像素沿x轴、y轴方向的物理距离;dx、dy均为摄像机自身属性,这些参数被称为摄像机的内参;
S422:以摄像机的光心Oc为坐标原点,建立摄像机坐标系;根据世界坐标系与摄像机坐标系的位置关系,通过坐标系刚性变换建立P点的摄像机坐标系与世界坐标系的变换关系,该变换关系如下:
(10)式中,P点的世界坐标为(Xw,Yw,Zw),P点的摄像机坐标为(Xc,Yc,Zc);R为两个坐标之间的旋转变换矩阵,该矩阵为单位正交的3×3矩阵;T矩阵为平移变换矩阵,该矩阵为3×1的列向量;
S423:根据摄像机的线性成像模型,图像坐标系与摄像机坐标系的关系为:
(11)式中,f为镜头的焦距;
S424:联立式(9)(10)(11)得到如下(12)式,即为摄像机光心Oc至点P的连线的直线方程;
(12)式中,A’表示内参:t为平移向量:r1 r2分别为旋转矩阵R的第一列和第二列;
S43:采用光平面标定技术在世界坐标系(Ow-XWYWZW)中标定出光平面方程;
设光平面方程为:
awxw+bwyw+cwzw+dw=0 (13);
(13)式中,(Xw,Yw,Zw)为点P在三维世界坐标方程中的表示;其中aw,bw,cw,dw为光平面方程的系数;
根据式(12)、(13)即能够计算出P点的三维世界坐标。
进一步地,步骤五具体为:两个支点的三维世界坐标分别为(Xw1,Yw1,Zw1)、(Xw2,Yw2,Zw2);则两个支点之间的距离d为:
有益效果:
(1)本发明的一种基于线激光的回转内腔尺寸检测方法,通过线激光发射器发送激光信号至回转内腔,通过摄像机采集摄像机拍摄线激光图片,进行图片分析,能够取代传统的三坐标测量仪和游标卡尺,对于各种结构匹配更灵活,检测状态评估更加精确。
(2)本发明的一种基于线激光的回转内腔尺寸检测方法,相对三坐标测量仪和游标卡尺检测成本大大降低,不用在三维空间拟合计算交点,而是在二维求出交点转换成三维再计算距离,效率显著提高。
(3)本发明提供的一种基于线激光的回转内腔尺寸检测方法,其实用性强,使用方便,适用面广。
附图说明
图1为本发明的流程图;
图2为本发明中实现线激光的回转内腔尺寸检测的装置的示意图;
图3为本发明中步骤一中采集的回转内腔的线激光轮廓灰度图片;
图4为经过步骤一的光条中心点提取方法获得的线激光轮廓图片;
图5为经过步骤二中所有散点与最低点连线的斜率散点图;
图6为经过直线拟合后的直线拟合效果图像;
图7为经过椭圆拟合后的椭圆拟合效果图像;
图8为本发明中涉及的测量坐标系示意图。
具体实施方式
下面结合附图对本发明进行具体的说明。
一种基于线激光的回转内腔尺寸检测方法,其特征在于:用测量工具对回转内腔的尺寸进行测量;所述测量工具如附图2所示,包括测量平台及安装于平台表面的传送带、线激光发射器以及摄像机;其中待测物放置在传送带表面,穿过线激光发发射器的发射区域,同时摄像机拍摄通过过程中的线激光轮廓图片,通过对线激光轮廓图片进行处理得到待测物即回转内腔的尺寸;所述通过对线激光轮廓图片进行处理得到待测物即回转内腔的尺寸具体如附图1所示,包括以下步骤:
步骤一:获得回转内腔的线激光轮廓图片,如附图3所示,所述回转内腔的线激光轮廓图片包括一条椭圆线条与一条直线线条;使用光条中心点提取方法获得轮廓形状图片的二维图像位置信息,如附图4所示,并且用坐标点集表示。
步骤二:筛选出符合条件的的二维坐标点,作为所需拟合椭圆与直线的轮廓坐标点;
步骤三:进行直线与椭圆的拟合,求出两个交点的二维图像坐标;
步骤四:将步骤三生成的交点坐标转换到空间坐标;
步骤五:计算两个交点之间的距离。
进一步地,步骤一中:所述轮廓二维坐标点集表示为:[(X1,Y1),(X2,Y2),(X3,Y3)........(Xn,Yn)],其中n为整数,(Xn,Yn)表示第n个散点对应的坐标点。
进一步地,步骤二具体包括:
S21:在步骤一生成的轮廓二维坐标点集中找出第r个点为图像最低点,其坐标表示为(Xmin,Ymin),;所述图像最低点为轮廓形状图片中y坐标数值最小的点;计算每个散点和最低点连线的斜率,其计算公式如下式(1):
则所有散点与最低点连线的斜率集合表示为:k=[k1 … kn] (2);
在(1)(2)式中,r、m均为小于等于n的整数。
上述过程中的图像最低点如附图4中的三个点中的中间点即为该图像的图像最低点。
忽略由于散点集中产生的极大值,在散点与最低点连线的斜率集中,判断出两个极大值的点,则该两个点即为相邻的直线与椭圆之间的拟合点的分界点;分界点如附图4中的三个点中的左右两个点,以及附图5中的两个点
由于光条中心点提取算法的特性,在提取的过程中,当计算点序号接近r时候会产生一个极大值点,如附图5中的横坐标为1400附近的点;因此需要忽略该极大值点。
S22:直线拟合点筛选;从图像最低点(Xmin,Ymin)的左半边、右半边分别筛选拟合点;过图像最低点(Xmin,Ymin)做垂直与X轴的直线,将图像分割为左半边与右半边;
当从图像最低点左半边计算时,此时计算包括图像最低点,斜率k取到极大值时对应的点为第i个点,i∈(1,r),则左半边直线拟合点的序号为第1个点到第i*α个点,其中α为拟合系数,且0.95<α<1;
当从图像最低点右半边计算时,斜率取到极大值时对应的点为第j个点,j∈(r,n),则右半边直线拟合点的序号为第j*β个点到第n个点,其中β为拟合系数,且1<β<1.05;
S23:椭圆拟合点筛选;将步骤S22生成的序号为第i*β到j*α个对应的散点作为椭圆拟合点。
进一步地,具体包括以下步骤:
S31:将步骤S22筛选出的直线拟合点,采用matlab中的polyfit指令进行表述,其表述的形式为直线方程y=kx+b;其生成的效果图如附图6所示;
S32:对步骤S23筛选出的椭圆拟合点采用最小二乘法拟合为如下的椭圆方程进行表述:
f(X,Y)=AX2+BXY+CY2+DX+EY+F (3);
其中椭圆方程进行表述包括以下步骤:
S32:对步骤S23筛选出的椭圆拟合点采用最小二乘法拟合为如下的椭圆方程进行表述:
f(X,Y)=AX2+BXY+CY2+DX+EY+F (3);
其中椭圆方程进行表述包括以下步骤:
S321:将平面内某点(X0,Y0)到椭圆方程所代表的曲线的距离表示为f(X0,Y0),对步骤S23筛选出的椭圆拟合点进行最小二乘处理,即优化目标为式(4):
S322:根据最小值来确定椭圆方程的系数A、B、C、D、E、F;
根据极值原理,对系数A、B、C、D、E、F求偏导,从而求得系数A,B,C,D,E,F的值;
S323:利用计算出的系数A,B,C,D,E,F的值计算出椭圆图像的五个参数:所述五个参数包括未知参数(θ,XC,YC)和形状参数(A,B);其中(XC,YC)为椭圆的几何中心;A、B分别为该椭圆的长半轴、短半轴;θ为椭圆的长短轴倾角;具体的计算过程如下:
几何中心:
长半轴短半轴:
长短轴倾角:
根据求出的五个参数,使用matlab函数ellipse即可画出椭圆效果图,如附图7所示。
S33:根据步骤S31、S32拟合成的直线方程与椭圆方程,得到椭圆与直线交点分别为(X1’,Y1’),(X2’,Y2’),如附图7中所描黑的两个点。
进一步地,步骤四具体包括以下步骤:
S41:将拟合算法中的图像坐标系转换成图像的像素坐标系;如附图8所示,为本发明中涉及的测量坐标系示意图。
将交点的坐标在X轴坐标不变的情况下,Y轴坐标减去图像的宽度h,即为:
X1’p=X1’
X2’p=X2’
Y1’P=|h-Y1’|
Y2’P=|h-Y2’|
则两个交点的图像坐标变成量纲为像素的像素坐标(X1’,Y1’P),(X2’,Y2’P);
S42:分别对两个交点的坐标进行三维坐标转换;每个交点的变换均包括以下步骤:
S421:设交点p的像素坐标为(xp,yp),不考虑镜头畸变,则其对应的图像坐标为(xn,yn),像素坐标系和图像坐标系原点均在图片左上角,则该两个坐标的变换关系为:
(9)式中,dx、dy分别表示单位像素沿x轴、y轴方向的物理距离;dx、dy均为摄像机自身属性,这些参数被称为摄像机的内参;
S422:以摄像机的光心Oc为坐标原点,建立摄像机坐标系;根据世界坐标系与摄像机坐标系的位置关系,通过坐标系刚性变换建立P点的摄像机坐标系与世界坐标系的变换关系,该变换关系如下:
(10)式中,P点的世界坐标为(Xw,Yw,Zw),P点的摄像机坐标为(Xc,Yc,Zc);R为两个坐标之间的旋转变换矩阵,该矩阵为单位正交的3×3矩阵;T矩阵为平移变换矩阵,该矩阵为3×1的列向量;
S423:根据摄像机的线性成像模型,图像坐标系与摄像机坐标系的关系为:
(11)式中,f为镜头的焦距;
S424:联立式(9)(10)(11)得到如下(12)式,即为摄像机光心Oc至点P的连线的直线方程;
(12)式中,A’表示内参:t为平移向量:r1 r2分别为旋转矩阵R的第一列和第二列;
S43:采用光平面标定技术在世界坐标系(Ow-XWYWZW)中标定出光平面方程;
设光平面方程为:
awxw+bwyw+cwzw+dw=0 (13);
(13)式中,(Xw,Yw,Zw)为点P在三维世界坐标方程中的表示;其中aw,bw,cw,dw为光平面方程的系数。
根据式(12)、(13)即能够计算出P点的三维世界坐标。
进一步地,步骤五具体为:两个支点的三维世界坐标分别为(Xw1,Yw1,Zw1)、(Xw2,Yw2,Zw2);则两个支点之间的距离d为:
虽然本发明已以较佳实施例公开如上,但它们并不是用来限定本发明的,任何熟习此技艺者,在不脱离本发明之精神和范围内,自当可作各种变化或润饰,因此本发明的保护范围应当以本申请的权利要求保护范围所界定的为准。
Claims (6)
1.一种基于线激光的回转内腔尺寸检测方法,其特征在于:用测量工具对回转内腔的尺寸进行测量;所述测量工具包括测量平台及安装于平台表面的传送带、线激光发射器以及摄像机;其中待测物放置在传送带表面,穿过线激光发发射器的发射区域,同时摄像机拍摄通过过程中的线激光轮廓图片,通过对线激光轮廓图片进行处理得到待测物即回转内腔的尺寸;所述通过对线激光轮廓图片进行处理得到待测物即回转内腔的尺寸具体包括以下步骤:
步骤一:获得回转内腔的线激光轮廓图片,所述回转内腔的线激光轮廓图片包括一条椭圆线条与一条直线线条;使用光条中心点提取方法获得轮廓形状图片的二维图像位置信息,并且用坐标点集表示;
步骤二:筛选出符合条件的的二维坐标点,作为所需拟合椭圆与直线的轮廓坐标点;
步骤三:进行直线与椭圆的拟合,求出两个交点的二维图像坐标;
步骤四:将步骤三生成的交点坐标转换到空间坐标;
步骤五:计算两个交点之间的距离。
2.根据权利要求1所述的一种基于线激光的回转内腔尺寸检测方法,其特征在于:步骤一中:所述轮廓二维坐标点集表示为:[(X1,Y1),(X2,Y2),(X3,Y3)........(Xn,Yn)],其中n为整数,(Xn,Yn)表示第n个散点对应的坐标点。
3.根据权利要求2所述的一种基于线激光的回转内腔尺寸检测方法,其特征在于:步骤二具体包括:
S21:在步骤一生成的轮廓二维坐标点集中找出第r个点为图像最低点,其坐标表示为(Xmin,Ymin),;所述图像最低点为轮廓形状图片中y坐标数值最小的点;计算每个散点和最低点连线的斜率,其计算公式如下式(1):
则所有散点与最低点连线的斜率集合表示为:k=[k1…kn] (2);
在(1)(2)式中,r、m均为小于等于n的整数;
忽略由于散点集中产生的极大值,在散点与最低点连线的斜率集中,判断出两个极大值的点,则该两个点即为相邻的直线与椭圆之间的拟合点的分界点;
S22:直线拟合点筛选;从图像最低点(Xmin,Ymin)的左半边、右半边分别筛选拟合点;过图像最低点(Xmin,Ymin)做垂直与X轴的直线,将图像分割为左半边与右半边;
当从图像最低点左半边计算时,此时计算包括图像最低点,斜率k取到极大值时对应的点为第i个点,i∈(1,r),则左半边直线拟合点的序号为第1个点到第i*α个点,其中α为拟合系数,且0.95<α<1;
当从图像最低点右半边计算时,斜率取到极大值时对应的点为第j个点,j∈(r,n),则右半边直线拟合点的序号为第j*β个点到第n个点,其中β为拟合系数,且1<β<1.05;
S23:椭圆拟合点筛选;将步骤S22生成的序号为第i*β到j*α个对应的散点作为椭圆拟合点。
4.根据权利要求3所述的一种基于线激光的回转内腔尺寸检测方法,其特征在于:步骤三具体包括以下步骤:
S31:将步骤S22筛选出的直线拟合点,采用matlab中的polyfit指令进行表述,其表述的形式为直线方程y=kx+b;
S32:对步骤S23筛选出的椭圆拟合点采用最小二乘法拟合为如下的椭圆方程进行表述:
f(X,Y)=AX2+BXY+CY2+DX+EY+F (3);
其中椭圆方程进行表述包括以下步骤:
S321:将平面内某点(X0,Y0)到椭圆方程所代表的曲线的距离表示为f(X0,Y0),对步骤S23筛选出的椭圆拟合点进行最小二乘处理,即优化目标为式(4):
S322:根据最小值来确定椭圆方程的系数A、B、C、D、E、F;
根据极值原理,对系数A、B、C、D、E、F求偏导,从而求得系数A,B,C,D,E,F的值;
S323:利用计算出的系数A,B,C,D,E,F的值计算出椭圆图像的五个参数:所述五个参数包括未知参数(θ,XC,YC)和形状参数(A,B);其中(XC,YC)为椭圆的几何中心;A、B分别为该椭圆的长半轴、短半轴;θ为椭圆的长短轴倾角;具体的计算过程如下:
几何中心:
长半轴短半轴:
长短轴倾角:
根据求出的五个参数,使用matlab函数ellipse即可画出椭圆效果图;
S33:根据步骤S31、S32拟合成的直线方程与椭圆方程,得到椭圆与直线交点分别为(X1’,Y1’),(X2’,Y2’)。
5.根据权利要求4所述的一种基于线激光的回转内腔尺寸检测方法,其特征在于:步骤四具体包括以下步骤:
S41:将拟合算法中的图像坐标系转换成图像的像素坐标系;
将交点的坐标在X轴坐标不变的情况下,Y轴坐标减去图像的宽度h,即为:
X1’p=X1’
X2’p=X2’
Y1’P=|h-Y1’|
Y2’P=|h-Y2’|
则两个交点的图像坐标变成量纲为像素的像素坐标(X1’,Y1’P),(X2’,Y2’P);
S42:分别对两个交点的坐标进行三维坐标转换;每个交点的变换均包括以下步骤:
S421:设交点p的像素坐标为(xp,yp),不考虑镜头畸变,则其对应的图像坐标为(xn,yn),像素坐标系和图像坐标系原点均在图片左上角,则该两个坐标的变换关系为:
(9)式中,dx、dy分别表示单位像素沿x轴、y轴方向的物理距离;dx、dy均为摄像机自身属性,这些参数被称为摄像机的内参;
S422:以摄像机的光心Oc为坐标原点,建立摄像机坐标系;根据世界坐标系与摄像机坐标系的位置关系,通过坐标系刚性变换建立P点的摄像机坐标系与世界坐标系的变换关系,该变换关系如下:
(10)式中,P点的世界坐标为(Xw,Yw,Zw),P点的摄像机坐标为(Xc,Yc,Zc);R为两个坐标之间的旋转变换矩阵,该矩阵为单位正交的3×3矩阵;T矩阵为平移变换矩阵,该矩阵为3×1的列向量;
S423:根据摄像机的线性成像模型,图像坐标系与摄像机坐标系的关系为:
(11)式中,f为镜头的焦距;
S424:联立式(9)(10)(11)得到如下(12)式,即为摄像机光心Oc至点P的连线的直线方程;
(12)式中,A’表示内参:t为平移向量:r1 r2分别为旋转矩阵R的第一列和第二列;
S43:采用光平面标定技术在世界坐标系(Ow-XWYWZW)中标定出光平面方程;
设光平面方程为:
awxw+bwyw+cwzw+dw=0 (13);
(13)式中,(Xw,Yw,Zw)为点P在三维世界坐标方程中的表示;其中aw,bw,cw,dw为光平面方程的系数;
根据式(12)、(13)即能够计算出P点的三维世界坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210501827.5A CN114894090B (zh) | 2022-05-09 | 2022-05-09 | 一种基于线激光的回转内腔尺寸检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210501827.5A CN114894090B (zh) | 2022-05-09 | 2022-05-09 | 一种基于线激光的回转内腔尺寸检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114894090A true CN114894090A (zh) | 2022-08-12 |
CN114894090B CN114894090B (zh) | 2023-06-23 |
Family
ID=82721407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210501827.5A Active CN114894090B (zh) | 2022-05-09 | 2022-05-09 | 一种基于线激光的回转内腔尺寸检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114894090B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10221074A (ja) * | 1997-01-31 | 1998-08-21 | Nippon Steel Corp | 三次元座標測定方法 |
CN103499302A (zh) * | 2013-09-27 | 2014-01-08 | 吉林大学 | 基于结构光视觉成像系统的凸轮轴直径尺寸在线测量方法 |
CN107014312A (zh) * | 2017-04-25 | 2017-08-04 | 西安交通大学 | 一种振镜式线激光扫描三维测量系统的整体标定方法 |
CN108562233A (zh) * | 2018-01-15 | 2018-09-21 | 吉林大学 | 利用二次曲线不变量的轴类零件直径尺寸在线测量方法 |
US10107594B1 (en) * | 2017-11-24 | 2018-10-23 | Huntercraft Limited | Analysis method of electronic spotting scope for automatically analyzing shooting accuracy |
CN110033491A (zh) * | 2019-04-15 | 2019-07-19 | 南京工程学院 | 一种相机标定方法 |
CN111121655A (zh) * | 2019-12-18 | 2020-05-08 | 浙江大学 | 一种共面等大多孔型工件位姿与孔径视觉检测方法 |
DE102019201526A1 (de) * | 2019-02-06 | 2020-08-06 | Ford Global Technologies, Llc | Verfahren und System zum Erfassen und Messen der Position eines Bauteils gegenüber einer Referenzposition sowie der Verschiebung und der Verdrehung eines sich relativ zu einem Bezugssystem bewegenden Bauteils |
CN113702384A (zh) * | 2021-07-19 | 2021-11-26 | 南京工程学院 | 一种回转构件表面缺陷检测装置、检测方法及标定方法 |
-
2022
- 2022-05-09 CN CN202210501827.5A patent/CN114894090B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10221074A (ja) * | 1997-01-31 | 1998-08-21 | Nippon Steel Corp | 三次元座標測定方法 |
CN103499302A (zh) * | 2013-09-27 | 2014-01-08 | 吉林大学 | 基于结构光视觉成像系统的凸轮轴直径尺寸在线测量方法 |
CN107014312A (zh) * | 2017-04-25 | 2017-08-04 | 西安交通大学 | 一种振镜式线激光扫描三维测量系统的整体标定方法 |
US10107594B1 (en) * | 2017-11-24 | 2018-10-23 | Huntercraft Limited | Analysis method of electronic spotting scope for automatically analyzing shooting accuracy |
CN108562233A (zh) * | 2018-01-15 | 2018-09-21 | 吉林大学 | 利用二次曲线不变量的轴类零件直径尺寸在线测量方法 |
DE102019201526A1 (de) * | 2019-02-06 | 2020-08-06 | Ford Global Technologies, Llc | Verfahren und System zum Erfassen und Messen der Position eines Bauteils gegenüber einer Referenzposition sowie der Verschiebung und der Verdrehung eines sich relativ zu einem Bezugssystem bewegenden Bauteils |
CN110033491A (zh) * | 2019-04-15 | 2019-07-19 | 南京工程学院 | 一种相机标定方法 |
CN111121655A (zh) * | 2019-12-18 | 2020-05-08 | 浙江大学 | 一种共面等大多孔型工件位姿与孔径视觉检测方法 |
CN113702384A (zh) * | 2021-07-19 | 2021-11-26 | 南京工程学院 | 一种回转构件表面缺陷检测装置、检测方法及标定方法 |
Non-Patent Citations (2)
Title |
---|
ZHU YONGGUO 等: "Dynamic Synthesis Correction of Deviation for Aircraft Wing-fuselage Docking Assembly Based on Laser Tracker and Machine Vision", JOURNAL OF MECHANICAL ENGINEERING, vol. 55, no. 24, pages 187 - 96 * |
黄鹏 等: "激光跟踪仪三维坐标转换综合优化方法", 计算机集成制造系统, vol. 21, no. 11, pages 2912 - 2920 * |
Also Published As
Publication number | Publication date |
---|---|
CN114894090B (zh) | 2023-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108986070B (zh) | 一种基于高速视频测量的岩石裂缝扩展实验监测方法 | |
CN109579695B (zh) | 一种基于异构立体视觉的零件测量方法 | |
CN112419429B (zh) | 一种基于多视角的大型工件表面缺陷检测标定方法 | |
CN110766669B (zh) | 一种基于多目视觉的管线测量方法 | |
CN109751964B (zh) | 一种高精度非接触式管径测量方法及装置 | |
CN106996748A (zh) | 一种基于双目视觉的轮径测量方法 | |
CN113137920A (zh) | 一种水下测量的设备与水下测量方法 | |
CN110260818B (zh) | 一种基于双目视觉的电子连接器鲁棒检测方法 | |
CN110966956A (zh) | 一种基于双目视觉的三维检测装置和方法 | |
CN110223355B (zh) | 一种基于双重极线约束的特征标志点匹配方法 | |
CN109974618B (zh) | 多传感器视觉测量系统的全局标定方法 | |
CN113375584B (zh) | 单目三维采集结构及单目三维多标距视觉系统、引伸计 | |
CN114627177A (zh) | 一种基于图像分割的飞机蒙皮间隙与阶差测量方法 | |
CN106918306A (zh) | 基于光场单相机的工业产品三维形貌实时检测系统 | |
Zhong et al. | Stereo-rectification and homography-transform-based stereo matching methods for stereo digital image correlation | |
JPH06137840A (ja) | 視覚センサの自動キャリブレーション装置 | |
CN113008158A (zh) | 多线激光轮胎花纹深度测量方法 | |
CN113702384A (zh) | 一种回转构件表面缺陷检测装置、检测方法及标定方法 | |
CN116402792A (zh) | 一种基于三维点云的空间孔位对接方法 | |
CN115540775A (zh) | 一种ccd单相机的3d视频引伸计 | |
CN111738971B (zh) | 一种基于线激光双目立体视觉的电路板立体扫描检测方法 | |
Su et al. | Measuring wear of the grinding wheel using machine vision | |
CN113916128A (zh) | 一种基于光笔式视觉测量系统的提高精度的方法 | |
CN110966937B (zh) | 一种基于激光视觉传感的大型构件三维构形拼接方法 | |
CN114894090B (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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20220812 Assignee: Hanzhong (Nanjing) Technology Co.,Ltd. Assignor: NANJING INSTITUTE OF TECHNOLOGY Contract record no.: X2024980000528 Denomination of invention: A method for measuring the size of a rotating inner cavity based on line laser Granted publication date: 20230623 License type: Common License Record date: 20240112 |