一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法
技术领域
本发明涉及智能装配技术领域,尤其涉及一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法。
背景技术
对缝点云间隙阶差提取技术应用领域广泛,包括机械工程、航空航天等,目前在飞机装配领域中存在许多对缝间隙阶差结构,但是其间隙阶差结构多样、尺寸不均匀,导致临界点难以准确提取,尽管3D扫描技术发展较快,扫描设备获取的点云数据依然存在较多噪声与密度缺陷,这给间隙阶差提取增加了难度。
除此之外,飞行器结构对缝的点云数据十分庞大,灵活的提取所需的间隙阶差也变得具有挑战性,针对此问题,本发明提出一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法。
发明内容
本发明的目的是为了解决现有技术中存在的缺点,而提出的一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法。
为了实现上述目的,本发明采用了如下技术方案:
一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法,利用数模缝隙边界离散点作为测量位置点,结合PCA及包围盒方法提取局部点云;利用PCA-SVM算法,使得超平面能够分割对缝两侧点云;利用对点云三角划分建立拓扑关系的方法,提取点云的边界点;根据点云的边界点与超平面的向量距离提取出边缘点,最后提取临界点,根据建立的数学模型提取间隙阶差,具体步骤如下:
S1:采用等弦长离散法获取模型缝隙边界离散点;
S2:利用SVM算法训练出分割超平面;
S3:分别对分割后两侧点云进行三角网格划分;
S4:采用弦高差法对点云进行去噪,排除跳点对后续曲率估计的影响;
S5:对所有扫描线建立点序关系,筛选出平面拟合点;
S6:计算间隙阶差值。
优选地:所述S1中,以离散点作为测量位置点,构建测量位置点的k邻域,利用PCA方法估计对缝点云的主法向,并以此为包围盒子坐标系的z轴,以扫描线为x轴,建立包围盒,包围盒的大小由T-Scan扫描线的宽度、线间距、点间距与缝隙结构曲率确定,具体步骤如下:
A2:对曲线进行等弦长离散,以
为圆心,做半径为
的圆,圆与曲线的交点即为
。对等弦长离散法进行编程实现,各离散点作为对缝测量位置点p;
A3:构建测量位置点的k邻域,对于点p的k邻域的拟合局部平面问题可以转化为局部邻域协方差矩阵M进行特征值分界,M最小特征值对应的特征向量即为点p的法向量,p点的局部邻域协方差矩阵如式
所示,式中:
为局部邻域重心;
A4:以测量位置点P为球心,包围盒宽度W为直径,取P点的局部邻域G
1。对G
1的所有点构建k邻域,半径为0.8mm,根据PCA算法,估计G
1所有点的法向量,视点方向统一,对所有法向量进行最小二乘滤波去除杂点,并取中值,获取局部邻域G
1的单位法向量
,根据PCA估计出的法向,以单位法向量
作为以P为原点的坐标系的z轴,以扫描线方向为x轴,其单位向量为
,则y轴单位向量为
,根据右手定则建立以点P为原点的坐标系S
1。在S
1坐标系下包围盒两个顶点坐标为P
min和P
max,根据公式
对S
1下的P
min和P
max进行坐标变换,得到在原世界坐标系S下的两顶点的坐标;
A5:通过世界坐标系下包围盒两个顶点,建立包围盒分割点云,其中,包围盒的内点满足条件:xmin≤x≤xmax,ymin≤y≤ymax,zmin≤z≤zmax,满足上述条件的点即为子点云。
优选地:所述S2中,利用PCA-SVM算法训练出分割超平面的方法为:
B1:利用PCA估计的法向量与扫描线方向构建对齐数据坐标系与训练数据坐标系;
B2:将两坐标系对齐变换,调整超平面法向,使得点云法向对齐,
B3:对点到超平面的向量距离排序,根据有无阶跃判断是否需要调整,并使得阶跃的两点的中点在超平面上,以调整超平面位移项。
优选地:所述S3中,三角网格划分方法为:采用邻接三角形数量判断法提取两侧三角网格的边界点,定义邻接三角形的个数为1的边为边界边,个数为2的边为内部边,对每一条边做边界边判断,然后从边界边中提取边界点。
优选地:所述S4中,其具体方法为:
C1:利用公式将两侧边界点投影到坐标系S1的xoy面;
C2:设其法向量为(a,b,c)。点(x
i,y
i,z
i)在xoy面ax+by+cz+d=0上的投影点为
,然后利用曲率筛选出对缝侧边缘角点;
C3:采用弦高差法对点云进行去噪,排除跳点对后续曲率估计的影响;
C4:分别遍历两侧点云,对所有点的k近邻点估计曲率,曲率较小的点视为角点,提取出对缝侧边缘点上下角点,即为边缘点;
C5:根据边缘点,提取点云临界点。
优选地:所述C4中,k的取值为:k=3。
优选地:所述C5中,根据所述C4中得到的边缘点,提取云临界点,由于扫描线点云是有序序列,其具体方法为:
C11:将得到的边缘点与点云所有点进行逐一比较,得出边缘点在所有点云中的序号,依照此序号即能获取单条扫描线边缘点之前的相邻若干点序号;
C12:根据第i个边缘点pi0确定的扫面线各点位置关系,pij为第i条扫描线边缘点pi0向前第j个点。
优选地:所述S5中,其具体方法如下:
D1:对单条扫描线点的k=3邻域求取曲率半径;
D2:设定阈值λ2,根据曲率半径判断是否存在大曲率边界,若存在大曲率,则获取最左侧对应点序号L,筛选出临界点区域,在临界点区域的左侧提取出平面拟合点pij(L+5≤J≤20+L),若不存在大曲率,则说明临界点区域不明显,直接提取左侧平面拟合点即可,并采用稳健特征值法获取拟合平面。
优选地:所述S6中,设定阈值d,设dij是第i条扫描线第j个点到平面α的距离,当dij>d,表明该点pij属于临界点区域,将属于临界点区域且dij最小的点作为临界点qij,表示第i条扫描线的第j(1≤j≤2)个临界点,随后即可根据数学模型,计算得出间隙阶差值。
本发明的有益效果为:
1.本发明对缝点云间隙阶差提取的精度、效率和稳定性,很好地解决了点云质量缺陷对特征点提取的影响。
附图说明
图1为本发明提出的一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法中测量位置点离散示意图;
图2为本发明提出的一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法中包围盒分割点云示意图;
图3为本发明提出的一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法中的PCA-SVM分割点云流程示意图;
图4为本发明提出的一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法中的超平面调整示意图;
图5为本发明提出的一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法中的边界点提取示意图;
图6为本发明提出的一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法中的边缘点提取示意图;
图7为本发明提出的一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法中的临界点提取示意图;
图8为本发明提出的一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法中的数学模型结构示意图。
具体实施方式
下面结合具体实施方式对本专利的技术方案作进一步详细地说明。
在本专利的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“设置”应做广义理解,例如,可以是固定相连、设置,也可以是可拆卸连接、设置,或一体地连接、设置。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本专利中的具体含义。
一种面向飞机蒙皮对缝的三维扫描线点云间隙阶差提取方法,利用数模缝隙边界离散点作为测量位置点,结合PCA及包围盒方法提取局部点云;利用PCA-SVM算法,使得超平面能够分割对缝两侧点云;利用对点云三角划分建立拓扑关系的方法,提取点云的边界点;根据点云的边界点与超平面的向量距离提取出边缘点,最后提取临界点,根据建立的数学模型提取间隙阶差,具体步骤如下:
S1:采用等弦长离散法获取模型缝隙边界离散点;
S2:利用SVM算法训练出分割超平面;
S3:分别对分割后两侧点云进行三角网格划分;
S4:采用弦高差法对点云进行去噪,排除跳点对后续曲率估计的影响;
S5:对所有扫描线建立点序关系,筛选出平面拟合点;
S6:计算间隙阶差值。
所述S1中,以离散点作为测量位置点,构建测量位置点的k邻域,利用PCA方法估计对缝点云的主法向,并以此为包围盒子坐标系的z轴,以扫描线为x轴,建立包围盒,包围盒的大小由T-Scan扫描线的宽度、线间距、点间距与缝隙结构曲率确定,具体步骤如下:
A2:对曲线进行等弦长离散,以
为圆心,做半径为
的圆,圆与曲线的交点即为
。对等弦长离散法进行编程实现,各离散点作为对缝测量位置点p;
A3:构建测量位置点的k邻域,对于点p的k邻域的拟合局部平面问题可以转化为局部邻域协方差矩阵M进行特征值分界,M最小特征值对应的特征向量即为点p的法向量,p点的局部邻域协方差矩阵如式
所示,式中:
为局部邻域重心;
A4:以测量位置点P为球心,包围盒宽度W为直径,取P点的局部邻域G
1。对G
1的所有点构建k邻域,半径为0.8mm,根据PCA算法,估计G
1所有点的法向量,视点方向统一,对所有法向量进行最小二乘滤波去除杂点,并取中值,获取局部邻域G
1的单位法向量
,根据PCA估计出的法向,以单位法向量
作为以P为原点的坐标系的z轴,以扫描线方向为x轴,其单位向量为
,则y轴单位向量为
,根据右手定则建立以点P为原点的坐标系S
1。在S
1坐标系下包围盒两个顶点坐标为P
min和P
max,根据公式
对S
1下的P
min和P
max进行坐标变换,得到在原世界坐标系S下的两顶点的坐标;
A5:通过世界坐标系下包围盒两个顶点,建立包围盒分割点云,其中,包围盒的内点满足条件:xmin≤x≤xmax,ymin≤y≤ymax,zmin≤z≤zmax,满足上述条件的点即为子点云。
所述S2中,利用PCA-SVM算法训练出分割超平面的方法为:
B1:利用PCA估计的法向量与扫描线方向构建对齐数据坐标系与训练数据坐标系;
B2:将两坐标系对齐变换,调整超平面法向,使得点云法向对齐,
B3:对点到超平面的向量距离排序,根据有无阶跃判断是否需要调整,并使得阶跃的两点的中点在超平面上,以调整超平面位移项。
所述S3中,三角网格划分方法为:采用邻接三角形数量判断法提取两侧三角网格的边界点,定义邻接三角形的个数为1的边为边界边,个数为2的边为内部边,对每一条边做边界边判断,然后从边界边中提取边界点。
所述S4中,其具体方法为:
C1:利用公式将两侧边界点投影到坐标系S1的xoy面;
C2:设其法向量为(a,b,c)。点(x
i,y
i,z
i)在xoy面ax+by+cz+d=0上的投影点为
,然后利用曲率筛选出对缝侧边缘角点;
C3:采用弦高差法对点云进行去噪,排除跳点对后续曲率估计的影响;
C4:分别遍历两侧点云,对所有点的k近邻点估计曲率,曲率较小的点视为角点,提取出对缝侧边缘点上下角点,即为边缘点;
C5:根据边缘点,提取点云临界点。
所述C4中,k的取值为:k=3。
所述C5中,根据所述C4中得到的边缘点,提取云临界点,由于扫描线点云是有序序列,其具体方法为:
C11:将得到的边缘点与点云所有点进行逐一比较,得出边缘点在所有点云中的序号,依照此序号即能获取单条扫描线边缘点之前的相邻若干点序号;
C12:根据第i个边缘点pi0确定的扫面线各点位置关系,pij为第i条扫描线边缘点pi0向前第j个点。
所述S5中,其具体方法如下:
D1:对单条扫描线点的k=3邻域求取曲率半径;
D2:设定阈值λ2,根据曲率半径判断是否存在大曲率边界,若存在大曲率,则获取最左侧对应点序号L,筛选出临界点区域,在临界点区域的左侧提取出平面拟合点pij(L+5≤J≤20+L),若不存在大曲率,则说明临界点区域不明显,直接提取左侧平面拟合点即可,并采用稳健特征值法获取拟合平面。
所述S6中,设定阈值d,设dij是第i条扫描线第j个点到平面α的距离,当dij>d,表明该点pij属于临界点区域,将属于临界点区域且dij最小的点作为临界点qij,表示第i条扫描线的第j(1≤j≤2)个临界点,随后即可根据数学模型,计算得出间隙阶差值。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。