CN108074277A - 一种面向数字化制造的弯管测量和余量定位方法 - Google Patents
一种面向数字化制造的弯管测量和余量定位方法 Download PDFInfo
- Publication number
- CN108074277A CN108074277A CN201611015476.8A CN201611015476A CN108074277A CN 108074277 A CN108074277 A CN 108074277A CN 201611015476 A CN201611015476 A CN 201611015476A CN 108074277 A CN108074277 A CN 108074277A
- Authority
- CN
- China
- Prior art keywords
- mtd
- point
- mtr
- bend pipe
- line
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- 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/10028—Range image; Depth image; 3D point clouds
-
- 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/20036—Morphological image processing
-
- 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/30172—Centreline of tubular or elongated structure
Abstract
本发明涉及一种面向数字化制造的弯管测量和余量定位方法。其步骤为:借助于激光跟踪仪和靶标进行多传感器的全局标定,将多传感器坐标系统一到全局坐标系;多个传感器严格同步采集弯管图像,边缘检测以及形态学细化得到弯管的中心线,借助于双向极线匹配算法进行中心线三维重建;将中心线离散为空间点云,计算控制点、直线段端点以及弧线段圆弧中心点,利用中心线作为引导线进行三维造型;最后利用最近点迭代法匹配设计模型和测量模型,计算余量位置,进而实现弯管的自寻位切割和焊接,加快弯管数字化制造进程。
Description
技术领域
本发明涉及计算机视觉技术领域,具体的说是一种面向数字化制造的弯管测量和余量定位方法。
背景技术
管路系统是航空、航天、汽车等工业领域重要的生命线,如人体血管一样,负责液压、冷气、燃气等介质的传输,目前对成型的弯管测量方式主要是接触式的,例如靠模法、三坐标测量机、关节式坐标测量机等,这些方法存在对比装置制造困难,容易造成损伤或变形,效率低,过分依赖操作者主观判断,无法重建三维模型以及无定位基准等缺陷。
先进制造技术是改造传统制造业的有效手段,现如今数字化制造技术是先进制造技术的核心技术,如何缩短研制时间、降低研制成本、提高产品质量是数字化制造的主要目的,然而弯管数字化制造是飞机快速制造的瓶颈,其关键问题在于弯管精确的三维重建以及自寻位切割。
唐正宗等人在专利“一种三维弯管多相机视觉检测方法及其系统”中公开一种弯管检测算法,该算法简单概述了弯管三维检测算法,并未对中心线的直线段和弧线段的分割方法进行分析,其次公开的算法中计算了弯管参数,仅仅利用中垂线中点得到了控制点,并没有计算直线段端点以及弧线段中心点。
发明内容
针对现有技术中存在的上述不足之处,本发明要解决的技术问题是提供一种面向数字化制造的弯管测量和余量定位方法,实现弯管的自寻位切割和焊接,加快弯管的快速数字化制造。
本发明为实现上述目的所采用的技术方案是:一种面向数字化制造的弯管测量和余量定位方法,包括以下步骤:
S1:利用单个激光跟踪仪和一个平面靶标对多个视觉传感器进行全局标定,将不同位置的视觉传感器的坐标系均统一到激光跟踪仪的全局坐标系下;
S2:同步触发多个视觉传感器采集弯管图像,采用边缘检测算法对弯管图像进行边缘检测,采用形态学细化方法获取弯管图像的中心线,剔除中心线上的分支,采用双向极线匹配算法进行中心线三维重建,获取弯管空间中心线;
S3:将弯管空间中心线离散为中心线三维点云,计算空间中心线的空间三维控制点;
S4:弯管由若干直线段和圆弧段组成,计算弯管空间中心线的每条直线段端点和每个圆弧段圆心位置;
S5:利用得到的直线段端点、空间圆弧圆心、圆弧半径以及管件管径分段进行扫掠得到复杂管件三维模型;
S6:采用最近点迭代ICP匹配算法,建立测量坐标系和设计坐标系之间的转换关系Γ和误差测度d(Γ),通过对误差测度d(Γ)进行最小化迭代,得到实际转换关系Γ,从而确定与设计坐标系的控制点相匹配的测量坐标系下的实际切割位置;
S7:按照给定的切割量对弯管两端的切割位置进行三维定位,对相机和机器人,机器人和切割机之间的位置关系进行标定,利用机器人协助将弯管放到切割机上进行自动切割。
所述剔除中心线上的分支,包括以下步骤:
形态学细化获得弯管中心线骨架,骨架为只包含一个像素宽的轮廓区域,将轮廓区域转换为8领域;用以下三种8领域形式进行说明:
其中1为前景,0为背景;
在8领域内,将存在三个或三个以上连接点的点进行标记:
其中,1为前景,0为背景,2为标记的连接点;
利用连接点将中心线轮廓区域分割成多个只包含一个或两个连接点的轮廓区域,统计每个轮廓区域的点数,剔除小于L个点的轮廓区域,保留的最终轮廓为管件的中心线轮廓,其中L=80。
所述计算空间三维控制点,包括以下步骤:
将空间中心线离散为中心线三维点云,对连续m个点(Pi,…Pi+m-1)作出拟合直线Li(ai,bi,ci,xi,yi,zi),再顺序对m个点(Pi+m,…Pi+2m-1)作出拟合直线Li+1(ai+1,bi+1,ci+1,xi+1,yi+1,zi+1);其中,ai,bi,ci分别为空间直线Li的参数,(xi,yi,zi)表示点Pi的坐标;ai+1,bi+1,ci+1分别为空间直线Li的参数,(xi+1,yi+1,zi+1)表示点Pi+m的坐标;
那么,两条空间直线Li与Li+1的空间夹角为若θi小于设定的阈值,则定义点云集PP1={Pi,…Pi+m-1,Pi+m,…Pi+2m-1},依次类推,直到θi大于设定的阈值时,点云集PP1={Pi,……Pi+lm-1}为同一条直线上的点,其中点云集PP1中点的个数为l×m,其中,选取m=10;
用上述方法将所有点云分割为直线段点云PP={PP1,PP2,…,PPn}和弧线段的点云;
由于弯管加工后会产生回弹力,相邻的两条直线并不共面,因此为了求取控制点,以点集PP1和点集PP2为例,利用最小二乘法将点集PP1、点集PP2的点拟合为一平面π:a0x+b0y+c0z=d0;其中,a0,b0,c0,d0为平面π的参数;
将第一组点云PP1拟合为一条直线L1,其参数方程为:在直线L1上取点A(x1,y1,z1),过该点且垂直于平面的直线方程为 该直线方程和平面π方程联立可得点A在平面π上的投影点B(xb,yb,zb),即 其中,a1,b1,c1为直线L1的参数;
同理,可以求得直线方程L1和平面方程π的交点为C(xc,yc,zc),因此过点B、C的直线方程为即为直线L1的投影方程;同理,得到相邻直线段L2的投影方程;
计算L1和L2的两条投影直线的交点,得到控制点D2;同理,可得控制点(D2,D3…,Dn),结合弯管两端的两个端点D1和Dn+1,可得整个空间中心线的所有控制点(D1,D2,…,Dn+1);其中,控制点的总数是n+1个,n为直线段PP={PP1,PP2,…,PPn}的个数。
所述计算空间中心线的所有直线段的端点和空间圆弧圆心位置,包括以下步骤:
由于相邻两条直线L1和L2的夹角有两个,因此利用空间方向向量D2D1和D2D3,计算夹角其中,D1,D2,D3为空间中心线上三个相邻的控制点;
对于非变弯曲半径的管路分支管件的加工,各个管件的弯曲半径是个常数r,因此,有由于直线段端点B1在直线L1上,点B1的坐标值为(a1T12+x1,b1T12+y1,c1T12+z1);其中比例参数同理,可以得到所有直线段的端点B(B1,B2,…,B2n);其中,直线段的端点总数为2n,n为直线段的个数;
求取三个平面:垂直于L1且过B1的平面、垂直于L2且过B2的平面、以及前面所述P1,P2点云拟合的平面π的交点,即得到圆弧的圆心Oc1;同理得到所有空间圆弧段的圆心位置O(Oc1,Oc2,…,Oc(n-1)),其中,圆弧的圆心总数为n-1。
所述实际转换关系Γ的计算包括以下步骤:
定义测量模型控制点集P={pj},j=1,…,N;定义设计模型控制点集X={xj},j=1,…,N,对点集P,X做如下处理:
其中
计算3×3的控制点的协方差矩阵其中,t表示矩阵的转置;
对加权协方差矩阵H进行SVD分解H=UΛVt;其中,旋转矩阵R和平移矩阵T分别为:
R=VUt;
通过不断迭代协方差矩阵H得到p′j=Rpj+T来最小化代价函数d(Г),得到最佳的刚体转换关系Г:
Г(pj)=Rpj+T
其中,Г包含了3×3的旋转矩阵R和3×1的平移矩阵T。
本发明具有以下有益效果及优点:
1.本发明所述的弯管自寻位切割方法实现弯管的三维重建。
2.本发明所述的弯管自寻位切割方法可以借助机器人实现弯管的自寻位切割,给出弯管的切割位置以及切割余量。
3.本发明所述的弯管自寻位切割方法可以有效地缩短弯管研制时间、降低研制成本、提高产品质量。
附图说明
图1为本发明方法流程图;
图2为本发明方法中弯管直线段端点和圆弧中心点位置定位图。
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
如图1所示,本发明所述一种面向数字化制造的弯管测量和余量定位方法的具体步骤如下:
步骤1:利用单个激光跟踪仪和一个平面靶标对多个视觉传感器进行全局标定,将不同位置的视觉传感器的坐标系均统一到激光跟踪仪的全局坐标系下;
步骤2:同步触发多个视觉传感器采集弯管图像,采用边缘检测算法对弯管图像进行边缘检测,采用形态学细化方法获取弯管图像的中心线,剔除中心线上的分支,采用双向极线匹配算法进行中心线三维重建,获取弯管空间中心线;
形态学细化获得弯管中心线骨架,骨架为只包含一个像素宽的轮廓区域,将轮廓区域转换为8领域。以下列三种8领域形式为例进行说明:
其中1为前景,0为背景;
在8领域内,将存在三个或三个以上连接点的点进行标记:
其中1为前景,0为背景,2为标记的连接点;
利用连接点将中心线轮廓区域分割成多个只包含一个或两个连接点的轮廓区域,统计每个轮廓区域的点数,剔除小于L个点的轮廓区域,保留的最终轮廓为管件的中心线轮廓。其中L=80;
步骤3:将弯管空间中心线离散为中心线三维点云,计算空间中心线的空间三维控制点;
将空间中心线离散为中心线三维点云,对连续m个点(Pi,…Pi+m-1)作出拟合直线Li(ai,bi,ci,xi,yi,zi),再顺序对m个点(Pi+m,…Pi+2m-1)作出拟合直线Li+1(ai+1,bi+1,ci+1,xi+1,yi+1,zi+1);其中,ai,bi,ci分别为空间直线Li的参数,(xi,yi,zi)表示点Pi的坐标;ai+1,bi+1,ci+1分别为空间直线Li的参数,(xi+1,yi+1,zi+1)表示点Pi+m的坐标;
那么,两条空间直线Li与Li+1的空间夹角为若θi小于设定的阈值,则定义点云集PP1={Pi,…Pi+m-1,Pi+m,…Pi+2m-1},依次类推,直到θi大于设定的阈值时,点云集PP1={Pi,……Pi+lm-1}为同一条直线上的点,其中点云集PP1中点的个数为l×m,其中,选取m=10;
用上述方法将所有点云分割为直线段点云PP={PP1,PP2,…,PPn}和弧线段的点云;
由于弯管加工后会产生回弹力,相邻的两条直线并不共面,因此为了求取控制点,以点集PP1和点集PP2为例,利用最小二乘法将点集PP1、点集PP2的点拟合为一平面π:a0x+b0y+c0z=d0;其中,a0,b0,c0,d0为平面π的参数;
将第一组点云PP1拟合为一条直线L1,其参数方程为:在直线L1上取点A(x1,y1,z1),过该点且垂直于平面的直线方程为 该直线方程和平面π方程联立可得点A在平面π上的投影点B(xb,yb,zb),即 其中,a1,b1,c1为直线L1的参数;
同理,可以求得直线方程L1和平面方程π的交点为C(xc,yc,zc),因此过点B、C的直线方程为即为直线L1的投影方程;同理,得到相邻直线段L2的投影方程;
如图2所示,计算L1和L2的两条投影直线的交点,得到控制点D2;同理,可得控制点(D2,D3…,Dn),结合弯管两端的两个端点D1和Dn+1,可得整个空间中心线的所有控制点(D1,D2,…,Dn+1);其中,控制点的总数是n+1个,n为直线段PP={PP1,PP2,…,PPn}的个数。
步骤4:弯管由若干直线段和圆弧段组成,计算弯管空间中心线的每条直线段端点和每个圆弧段圆心位置;
由于相邻两条直线L1和L2的夹角有两个,因此利用空间方向向量D2D1和D2D3,计算夹角其中,D1,D2,D3为空间中心线上三个相邻的控制点;
对于非变弯曲半径的管路分支管件的加工,各个管件的弯曲半径是个常数r,因此,有由于直线段端点B1在直线L1上,点B1的坐标值为(a1T12+x1,b1T12+y1,c1T12+z1);其中比例参数同理,可以得到所有直线段的端点B(B1,B2,…,B2n);其中,直线段的端点总数为2n,n为直线段的个数;求取三个平面:垂直于L1且过B1的平面、垂直于L2且过B2的平面、以及前面所述P1,P2点云拟合的平面π的交点,即得到圆弧的圆心Oc1;同理得到所有空间圆弧段的圆心位置O(Oc1,Oc2,…,Oc(n-1)),其中,圆弧的圆心总数为n-1。
步骤5:利用得到的直线段端点、空间圆弧圆心、圆弧半径以及管件管径分段进行扫掠得到复杂管件三维模型;
步骤6:采用最近点迭代ICP匹配算法,建立测量坐标系和设计坐标系之间的转换关系Γ和误差测度d(Γ),通过对误差测度d(Γ)进行最小化迭代,得到实际转换关系Γ,从而确定与设计坐标系的控制点相匹配的测量坐标系下的实际切割位置;
定义测量模型控制点集P={pj},j=1,…,N;定义设计模型控制点集X={xj},j=1,…,N,对点集P,X做如下处理:
其中
计算3×3的控制点的协方差矩阵其中,t表示矩阵的转置;
对加权协方差矩阵H进行SVD分解H=UΛVt;其中,旋转矩阵R和平移矩阵T分别为:
R=VUt;
通过不断迭代协方差矩阵H得到p′j=Rpj+T来最小化代价函数d(Γ),得到最佳的刚体转换关系Γ,:
Γ(pj)=Rpj+T
其中,Γ包含了3×3的旋转矩阵R和3×1的平移矩阵T;
步骤7:按照设计模型给定的切割量对弯管两端的切割位置进行三维定位,对相机和机器人,机器人和切割机之间的位置关系进行标定,利用机器人协助将弯管放到切割机上进行自动切割。
为验证本发明方法的有效性和正确性,采用AVT公司的Manta系列相机(分辨率:2048×1088),相机配置了Computar公司的8mm精工镜头,所有的仿真实验均在Windows 7操作系统下采用Visual Studio 2008软件实现。
Claims (5)
1.一种面向数字化制造的弯管测量和余量定位方法,其特征在于,包括以下步骤:
S1:利用单个激光跟踪仪和一个平面靶标对多个视觉传感器进行全局标定,将不同位置的视觉传感器的坐标系均统一到激光跟踪仪的全局坐标系下;
S2:同步触发多个视觉传感器采集弯管图像,采用边缘检测算法对弯管图像进行边缘检测,采用形态学细化方法获取弯管图像的中心线,剔除中心线上的分支,采用双向极线匹配算法进行中心线三维重建,获取弯管空间中心线;
S3:将弯管空间中心线离散为中心线三维点云,计算空间中心线的空间三维控制点;
S4:弯管由若干直线段和圆弧段组成,计算弯管空间中心线的每条直线段端点和每个圆弧段圆心位置;
S5:利用得到的直线段端点、空间圆弧圆心、圆弧半径以及管件管径分段进行扫掠得到复杂管件三维模型;
S6:采用最近点迭代ICP匹配算法,建立测量坐标系和设计坐标系之间的转换关系Г和误差测度d(Γ),通过对误差测度d(Γ)进行最小化迭代,得到实际转换关系Γ,从而确定与设计坐标系的控制点相匹配的测量坐标系下的实际切割位置;
S7:按照给定的切割量对弯管两端的切割位置进行三维定位,对相机和机器人,机器人和切割机之间的位置关系进行标定,利用机器人协助将弯管放到切割机上进行自动切割。
2.根据权利要求1所述的一种面向数字化制造的弯管测量和余量定位方法,其特征在于,所述剔除中心线上的分支,包括以下步骤:
形态学细化获得弯管中心线骨架,骨架为只包含一个像素宽的轮廓区域,将轮廓区域转换为8领域;用以下三种8领域形式进行说明:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
</mtable>
</mfenced>
其中1为前景,0为背景;
在8领域内,将存在三个或三个以上连接点的点进行标记:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,1为前景,0为背景,2为标记的连接点;
利用连接点将中心线轮廓区域分割成多个只包含一个或两个连接点的轮廓区域,统计每个轮廓区域的点数,剔除小于L个点的轮廓区域,保留的最终轮廓为管件的中心线轮廓,其中L=80。
3.根据权利要求1所述的一种面向数字化制造的弯管测量和余量定位方法,其特征在于,所述计算空间三维控制点,包括以下步骤:
将空间中心线离散为中心线三维点云,对连续m个点(Pi,…Pi+m-1)作出拟合直线Li(ai,bi,ci,xi,yi,zi),再顺序对m个点(Pi+m,…Pi+2m-1)作出拟合直线Li+1(ai+1,bi+1,ci+1,xi+1,yi+1,zi+1);其中,ai,bi,ci分别为空间直线Li的参数,(xi,yi,zi)表示点Pi的坐标;ai+1,bi+1,ci+1分别为空间直线Li的参数,(xi+1,yi+1,zi+1)表示点Pi+m的坐标;
那么,两条空间直线Li与Li+1的空间夹角为若θi小于设定的阈值,则定义点云集PP1={Pi,…Pi+m-1,Pi+m,…Pi+2m-1},依次类推,直到θi大于设定的阈值时,点云集PP1={Pi,……Pi+lm-1}为同一条直线上的点,其中点云集PP1中点的个数为l×m,其中,选取m=10;
用上述方法将所有点云分割为直线段点云PP={PP1,PP2,…,PPn}和弧线段的点云;
由于弯管加工后会产生回弹力,相邻的两条直线并不共面,因此为了求取控制点,以点集PP1和点集PP2为例,利用最小二乘法将点集PP1、点集PP2的点拟合为一平面π:a0x+b0y+c0z=d0;其中,a0,b0,c0,d0为平面π的参数;
将第一组点云PP1拟合为一条直线L1,其参数方程为:在直线L1上取点A(x1,y1,z1),过该点且垂直于平面的直线方程为L2:该直线方程和平面π方程联立可得点A在平面π上的投影点B(xb,yb,zb),即 其中,a1,b1,c1为直线L1的参数;
同理,可以求得直线方程L1和平面方程π的交点为C(xc,yc,zc),因此过点B、C的直线方程为即为直线L1的投影方程;同理,得到相邻直线段L2的投影方程;
计算L1和L2的两条投影直线的交点,得到控制点D2;同理,可得控制点(D2,D3…,Dn),结合弯管两端的两个端点D1和Dn+1,可得整个空间中心线的所有控制点(D1,D2,…,Dn+1);其中,控制点的总数是n+1个,n为直线段PP={PP1,PP2,…,PPn}的个数。
4.根据权利要求1所述的一种面向数字化制造的弯管测量和余量定位方法,其特征在于,所述计算空间中心线的所有直线段的端点和空间圆弧圆心位置,包括以下步骤:
由于相邻两条直线L1和L2的夹角有两个,因此利用空间方向向量D2D1和D2D3,计算夹角其中,D1,D2,D3为空间中心线上三个相邻的控制点;
对于非变弯曲半径的管路分支管件的加工,各个管件的弯曲半径是个常数r,因此,有由于直线段端点B1在直线L1上,点B1的坐标值为(a1T12+x1,b1T12+y1,c1T12+z1);其中比例参数同理,可以得到所有直线段的端点B(B1,B2,…,B2n);其中,直线段的端点总数为2n,n为直线段的个数;
求取三个平面:垂直于L1且过B1的平面、垂直于L2且过B2的平面、以及前面所述P1,P2点云拟合的平面π的交点,即得到圆弧的圆心Oc1;同理得到所有空间圆弧段的圆心位置O(Oc1,Oc2,…,Oc(n-1)),其中,圆弧的圆心总数为n-1。
5.根据权利要求1所述的一种面向数字化制造的弯管测量和余量定位方法,其特征在于,所述实际转换关系Γ的计算包括以下步骤:
定义测量模型控制点集P={pj},j=1,…,N;定义设计模型控制点集X={xj},j=1,…,N,对点集P,X做如下处理:
<mrow>
<msubsup>
<mi>p</mi>
<mi>j</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>p</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mo>,</mo>
<msubsup>
<mi>x</mi>
<mi>j</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
其中
计算3×3的控制点的协方差矩阵其中,t表示矩阵的转置;
对加权协方差矩阵H进行SVD分解H=UΛVt;其中,旋转矩阵R和平移矩阵T分别为:
R=VUt;
<mrow>
<mi>T</mi>
<mo>=</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<mi>R</mi>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mo>;</mo>
</mrow>
通过不断迭代协方差矩阵H得到p′j=Rpj+T来最小化代价函数d(Г),得到最佳的刚体转换关系Г:
<mrow>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>&Gamma;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
Г(pj)=Rpj+T
其中,Γ包含了3×3的旋转矩阵R和3×1的平移矩阵T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611015476.8A CN108074277B (zh) | 2016-11-18 | 2016-11-18 | 一种面向数字化制造的弯管测量和余量定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611015476.8A CN108074277B (zh) | 2016-11-18 | 2016-11-18 | 一种面向数字化制造的弯管测量和余量定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108074277A true CN108074277A (zh) | 2018-05-25 |
CN108074277B CN108074277B (zh) | 2021-09-21 |
Family
ID=62160119
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611015476.8A Active CN108074277B (zh) | 2016-11-18 | 2016-11-18 | 一种面向数字化制造的弯管测量和余量定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108074277B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110849291A (zh) * | 2019-11-26 | 2020-02-28 | 二重(德阳)重型装备有限公司 | 大型弯管弯曲半径检测方法 |
CN111047702A (zh) * | 2019-12-18 | 2020-04-21 | 成都飞机工业(集团)有限责任公司 | 一种基于双目视觉的法兰弯管自动焊接方法 |
CN111489432A (zh) * | 2019-12-16 | 2020-08-04 | 西安航天发动机有限公司 | 一种基于点云数据的弯管重建及余量计算方法 |
CN111754462A (zh) * | 2020-05-28 | 2020-10-09 | 新拓三维技术(深圳)有限公司 | 一种三维弯管的视觉检测方法及系统 |
CN112729150A (zh) * | 2019-10-14 | 2021-04-30 | 中国科学院沈阳自动化研究所 | 机器人夹持下的导管测量方法 |
CN112950627A (zh) * | 2021-04-01 | 2021-06-11 | 上海柏楚电子科技股份有限公司 | 用于激光切割的检测及控制方法和系统 |
CN113205086A (zh) * | 2021-07-05 | 2021-08-03 | 武汉瀚迈科技有限公司 | 一种基于椭圆拟合的圆截面弯管类零件特征参数识别算法 |
CN114677374A (zh) * | 2022-05-27 | 2022-06-28 | 杭州键嘉机器人有限公司 | 一种下颌神经管中心线提取与半径计算的方法 |
CN114850691A (zh) * | 2022-04-12 | 2022-08-05 | 西安航天发动机有限公司 | 一种定制化导管余量自动去除工艺方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2295931A1 (de) * | 2009-09-02 | 2011-03-16 | MSG Maschinenbau GmbH | Vorrichtung und Verfahren zum Vermessen der Form eines Gegenstands |
CN103862330A (zh) * | 2012-12-14 | 2014-06-18 | 中国科学院沈阳自动化研究所 | 基于机器视觉的弯管磁研磨自动导航方法 |
CN104036513A (zh) * | 2014-06-26 | 2014-09-10 | 西安新拓三维光测科技有限公司 | 一种三维弯管多相机视觉检测方法及其系统 |
CN105809736A (zh) * | 2016-03-14 | 2016-07-27 | 北京理工大学 | 一种管路的三维重建方法及装置 |
CN105913483A (zh) * | 2016-03-31 | 2016-08-31 | 百度在线网络技术(北京)有限公司 | 一种三维交叉道路模型生成的方法及装置 |
-
2016
- 2016-11-18 CN CN201611015476.8A patent/CN108074277B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2295931A1 (de) * | 2009-09-02 | 2011-03-16 | MSG Maschinenbau GmbH | Vorrichtung und Verfahren zum Vermessen der Form eines Gegenstands |
CN103862330A (zh) * | 2012-12-14 | 2014-06-18 | 中国科学院沈阳自动化研究所 | 基于机器视觉的弯管磁研磨自动导航方法 |
CN104036513A (zh) * | 2014-06-26 | 2014-09-10 | 西安新拓三维光测科技有限公司 | 一种三维弯管多相机视觉检测方法及其系统 |
CN105809736A (zh) * | 2016-03-14 | 2016-07-27 | 北京理工大学 | 一种管路的三维重建方法及装置 |
CN105913483A (zh) * | 2016-03-31 | 2016-08-31 | 百度在线网络技术(北京)有限公司 | 一种三维交叉道路模型生成的方法及装置 |
Non-Patent Citations (1)
Title |
---|
张天: "基于多目视觉的管路数字化测量方法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112729150A (zh) * | 2019-10-14 | 2021-04-30 | 中国科学院沈阳自动化研究所 | 机器人夹持下的导管测量方法 |
CN112729150B (zh) * | 2019-10-14 | 2021-10-22 | 中国科学院沈阳自动化研究所 | 机器人夹持下的导管测量方法 |
CN110849291A (zh) * | 2019-11-26 | 2020-02-28 | 二重(德阳)重型装备有限公司 | 大型弯管弯曲半径检测方法 |
CN110849291B (zh) * | 2019-11-26 | 2021-05-11 | 二重(德阳)重型装备有限公司 | 大型弯管弯曲半径检测方法 |
CN111489432A (zh) * | 2019-12-16 | 2020-08-04 | 西安航天发动机有限公司 | 一种基于点云数据的弯管重建及余量计算方法 |
CN111047702A (zh) * | 2019-12-18 | 2020-04-21 | 成都飞机工业(集团)有限责任公司 | 一种基于双目视觉的法兰弯管自动焊接方法 |
CN111754462A (zh) * | 2020-05-28 | 2020-10-09 | 新拓三维技术(深圳)有限公司 | 一种三维弯管的视觉检测方法及系统 |
CN112950627A (zh) * | 2021-04-01 | 2021-06-11 | 上海柏楚电子科技股份有限公司 | 用于激光切割的检测及控制方法和系统 |
CN112950627B (zh) * | 2021-04-01 | 2023-01-20 | 上海柏楚电子科技股份有限公司 | 用于激光切割的检测及控制方法和系统 |
CN113205086A (zh) * | 2021-07-05 | 2021-08-03 | 武汉瀚迈科技有限公司 | 一种基于椭圆拟合的圆截面弯管类零件特征参数识别算法 |
CN114850691A (zh) * | 2022-04-12 | 2022-08-05 | 西安航天发动机有限公司 | 一种定制化导管余量自动去除工艺方法 |
CN114677374A (zh) * | 2022-05-27 | 2022-06-28 | 杭州键嘉机器人有限公司 | 一种下颌神经管中心线提取与半径计算的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108074277B (zh) | 2021-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108074277A (zh) | 一种面向数字化制造的弯管测量和余量定位方法 | |
CN104700451B (zh) | 基于迭代就近点算法的点云配准方法 | |
CN106023298B (zh) | 基于局部泊松曲面重建的点云刚性配准方法 | |
CN104748683B (zh) | 一种数控机床工件在线自动测量装置及测量方法 | |
CN111080627A (zh) | 一种基于深度学习的2d+3d大飞机外形缺陷检测与分析方法 | |
CN110070567B (zh) | 一种地面激光点云配准方法 | |
JP2014500484A (ja) | 連続b‐スプライン変形を利用してタイヤ表面の三次元画像を前処理する方法 | |
CN107330901B (zh) | 一种基于骨架的物体构件分解方法 | |
CN103278138A (zh) | 一种复杂结构薄部件三维位置及姿态的测量方法 | |
CN107680137A (zh) | 一种基于球体投影公切线的多相机标定及参数优化方法 | |
CN104778716B (zh) | 基于单幅图像的卡车车厢体积测量方法 | |
JPH08136220A (ja) | 物品の位置検出方法およびその装置 | |
CN115063465B (zh) | 一种基于激光雷达的无人车行驶路况建模方法 | |
CN111489432A (zh) | 一种基于点云数据的弯管重建及余量计算方法 | |
CN104615880B (zh) | 一种三维激光雷达点云匹配的快速icp方法 | |
CN112508895A (zh) | 一种基于曲面配准的螺旋桨叶片质量评估方法 | |
CN112465966A (zh) | 一种集倾斜摄影测量与三维激光扫描的陡崖三维建模方法 | |
CN114995398A (zh) | 路径生成方法、装置、存储介质、处理器及电子装置 | |
CN106813569A (zh) | 一种基于线结构光的汽车轮胎三维定位方法 | |
CN113689394A (zh) | 一种基于三维点云的隧道初支平整度检测方法 | |
Xu et al. | A new welding path planning method based on point cloud and deep learning | |
CN109191583A (zh) | 一种基于各向异性mls曲面精确对齐方法 | |
CN112037282A (zh) | 一种基于关键点以及骨架的飞机姿态估计方法及系统 | |
CN115578534A (zh) | 用于焊缝的3d模型重构方法 | |
CN108536932A (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 |