CN107633502A - 一种轴孔装配自动对中的靶心识别方法 - Google Patents
一种轴孔装配自动对中的靶心识别方法 Download PDFInfo
- Publication number
- CN107633502A CN107633502A CN201710623740.4A CN201710623740A CN107633502A CN 107633502 A CN107633502 A CN 107633502A CN 201710623740 A CN201710623740 A CN 201710623740A CN 107633502 A CN107633502 A CN 107633502A
- Authority
- CN
- China
- Prior art keywords
- point set
- point
- cylinder
- ladder
- face
- 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
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种轴孔装配自动对中的靶心识别方法,该方法包括以下步骤:a.在待装配孔孔内的后端安装阶梯靶标,梯靶标的轴线与待装配孔处于共轴状态,采集阶梯靶标图像;b.步骤a中获得到阶梯靶标图像进行处理,得到边缘平滑过渡的阶梯靶标图像;c.对步骤b中边缘平滑过渡的阶梯靶标图像处理,得到一组拟合椭圆中心点集合;d.将步骤c中得到的拟合椭圆中心点集合分为第一圆柱体前端面点集和第二圆柱体前端面点集。步骤e.将步骤d中的所述第一圆柱体前端面点集和第二圆柱体前端面点集,分别运用最小二乘法进行椭圆拟合,的阶梯靶标的轴线。该轴孔装配自动对中的靶心识别方法,计算复杂度低,图像处理效率高且稳定。
Description
技术领域
本发明属于图像特征参数提取技术领域,具体涉及一种轴孔装配自动对中的靶心识别方法。
背景技术
随着国际制造技术取得了突飞猛进的发展,机械装配的质量占据越来越重要的地位,是研究制造高性能装备的关键。旋转机构在目前汽车、飞机设计过程中应用非常普遍,如何保证该机构孔轴的装配精准性是一直困扰当前一线工程技术人员的难题。装配之前,孔轴处于分离状态,装配时要求孔轴对中,任何装配过程的误差都有可能降低产品使用寿命,加剧磨损,对工件的使用造成很大的危害。
目前旋转机械孔轴的对中装配过程,传统的方法是由工人们肉眼观察孔轴的配合间隙,为了满足同轴度、轴间隙等硬性指标的要求,有时必须要反复调整,耗费大量的工时成本,并且过多人为因素的干预导致装配精度难以保证。为了降低工人的劳动强度,减少人为因素对装配精度的影响,在装配线规划中,针对孔轴对中装配过程应用机器视觉技术将会显著提高装配效果。
发明内容
本发明所要解决的技术问题在于针对上述现有技术的不足,提供一种计算复杂度低,图像处理效率高且稳定的轴孔装配自动对中的靶心识别方法。
为解决上述技术问题,本发明采用的技术方案是,一种轴孔装配自动对中的靶心识别方法,该方法包括以下步骤:
步骤a.使用阶梯靶标,阶梯靶标包括同轴设置且位于前部的第一圆柱体和位于后部的第二圆柱体,且第一圆柱体的直径小于第二圆柱体的直径,第一圆柱体和第二圆柱体的前端面上、且绕其边缘一周均设置有4A个圆形图案,且对应端面上的圆形图案的圆心分别位于对应的圆的圆周上,所形成的对应的两个圆为同心圆,其中:A为大于1的自然数,且两个圆柱体前端面的A的取值不要求相同;
在待装配孔内的安装阶梯靶标,该阶梯靶标上的图案朝向待装配孔前端,且阶梯靶标轴线与待装配孔处于共轴状态,在待装配孔前端朝向孔内采集阶梯靶标,得阶梯靶标图像;
步骤b.对步骤a中获得到阶梯靶标图像进行处理,得到边缘平滑过渡的阶梯靶标图像;
步骤c.对步骤b中边缘平滑过渡的阶梯靶标图像处理,获得阶梯靶标图像边缘轮廓集合,对边缘轮廓集合中的每个子集分别进行最小二乘椭圆拟合,得到一组拟合椭圆中心点集合;
步骤d.将步骤c中得到的拟合椭圆中心点集合分为第一圆柱体前端面点集和第二圆柱体前端面点集;
步骤e.将步骤d中的所述第一圆柱体前端面点集和第二圆柱体前端面点集,分别运用最小二乘法进行椭圆拟合,即得到对应的第一圆柱体前端面中心点和第二圆柱体前端面中心点;连接两个中心点,即为阶梯靶标的轴线,靶标轴线的位姿即为待装配孔的位姿。
进一步地,该步骤d的具体过程如下:
步骤d1.设处理后的阶梯靶标图像左上角第一个像素点为图像坐标系原点,水平向右为X轴正方向,竖直向下为Y轴正方向;在图像坐标系中,将所述拟合椭圆中心点按照距离图像坐标系原点由小到大的原则重新排序,获得有序点集Si(i=1,2,3,4,5,6…);
步骤d2.遍历有序点集Si,搜索满足正方形分布的坐标点,首先设S1为正方形的一个顶点Pt1,以顶点Pt1为对角线的一个点,遍历有序点集Si,在该对角线上寻找正方形的第二个顶点Pt2,该对角线的斜率k满足如下要求:0.5<k<3;
步骤d3.搜索正方形的第三个顶点,设Pt1和Pt2两点间对角线对应的一次函数为:y1=ax+b,该对角线中点为其中:Pt1(m1,n1)和Pt2(m2,n2);Pt1到M点的距离为L,过M点作Pt1与Pt2连线的中垂线,得到该中垂线的一次函数为:设中垂线上一点为:遍历有序点集Si,若存在Si满足则将Si设为Pt3点;
若Si不满足重新执行步骤d2,搜寻其他满足条件的Pt1、Pt2再执行步骤d3,直至存在Si满足则将Si设为Pt3点;
步骤d4.搜寻正方形的第四个顶点,以点Pt1,Pt3确定向量VecX,点Pt2,Pt3确定向量VecY,设Pt1为坐标原点,向量VecX所在方向为X轴,向量VecY所在方向为Y轴,构建新坐标系;由正方形顶点排列特性得到参考点Pt4’在新构建的坐标系中可表示为:Pt4’=Pt1+VecY,遍历有序点集Si,保存满足的Si点集,设为点集Qi,点集Qi中满足||Pt4'-Qi||取得最小值的点即为Pt4点,若存在Pt4点则执行步骤d6,若不存在Pt4点,则执行步骤d5;
步骤d5.如果以S1点作为第一个顶点Pt1,不存在另外3个点满足步骤d2、d3、d4,则遍历点集Si,依次令Si(i=2,3,4…)为Pt1,再继续步骤d2、d3、d4;
步骤d6.将搜寻到的满足正方形分布的4个坐标点作为采集的阶梯靶标图像的参照点,计算椭圆圆心Si(i=1,2,3…)到参照点的距离Gi(i=1,2,3…),并求得所有椭圆圆心到参照点的距离和的均值Ver;
若Gi>Ver,且则该点属于外圆;
若Gi<Ver且则该点属于内圆;
其中比例值P为阶梯靶标的内外圆半径之比;
得到内圆点集和外圆点集,内圆点集为第一圆柱体前端面点集,外圆点集为第二圆柱体前端面点集;
其中:内圆为第一圆柱体前端面上的圆形图案的圆心形成的圆,外圆为第二圆柱体前端面上的前端面上的圆形图案的圆心形成的圆;
步骤d7.若阶梯靶标倾斜角度过大,前端面对后端面的圆形图案遮挡,经过分类处理后只能得到第一圆柱体前端面点集,则再次遍历点集Si剔除得到的第一圆柱体前端面点集,即可得到第二圆柱体前端面点集。
进一步地,该步骤b的具体过程如下:
步骤b1.对步骤a中获得到阶梯靶标图像进行二值化处理,得到二值化阶梯靶标图像;
步骤b2.对b1中获得的二值化阶梯靶标图像进行滤波处理,得到边缘平滑过渡的二值化阶梯靶标图像。
进一步地,该步骤b1中采用最大类间方差法计算阈值,以得到二值化阶梯靶标图像。
进一步地,该步骤b2中的具体处理过程如下:使用n*n的像素矩形沿着二值化阶梯靶标图像水平向右和竖直向下两个方向移动,将中心像素灰度值替换为相邻n*n像素矩形框内像素灰度的平均值,得到边缘平滑过渡的二值化阶梯靶标图像;其中n为大于0的自然数。
进一步地,该步骤c中,将边缘平滑过渡的二值化阶梯靶标图像采用Canny算子获得阶梯靶标图像边缘轮廓集合。
进一步地,每一个圆形图案内均设置有光源,所述第一圆柱体和第二圆柱体的前端面上均贴合安装有挡光板,每一个挡光板上均对应开设有圆孔,圆孔的位置与对应的圆形图案相一致,且圆孔的直径小于圆形图案的直径。
本发明一种轴孔装配自动对中的靶心识别方法具有如下优点:数学推导简单,计算复杂度低,图像处理效率高且稳定,可以有效的完成对中过程靶心特征的识别,反馈的位姿状态信息可用于自动装配过程的实时调整,极大提高了孔轴装配的效率和精度。
附图说明
图1是本发明一种轴孔装配自动对中的靶心识别方法中待装配孔和待装配孔使用的阶梯靶标的结构示意图;
图2是本发明中工业相机采集的原始靶标图像;
图3是本发明中经过二值化处理后的靶标图像;
图4是本发明中经过边缘提取后的轮廓图像;
图5是本发明中基于正方形分布阵列的分类原理图;
图6是本发明中阶梯靶标两端面中心识别效果图。
其中:1.阶梯靶标;1-1.第一圆柱体;1-2.第二圆柱体;1-3圆形图案;1-4.挡光板。
具体实施方式
本发明一种轴孔装配自动对中的靶心识别方法,该方法包括以下步骤:
步骤a.使用阶梯靶标1,如图1所示,所述阶梯靶标1包括同轴设置且位于前部的第一圆柱体1-1和位于后部的第二圆柱体1-2,且第一圆柱体1-1的直径小于第二圆柱体1-2的直径,所述第一圆柱体1-1和第二圆柱体1-2的前端面上、且绕其边缘一周均设置有4A个圆形图案1-3,且对应端面上的圆形图案1-3的圆心分别位于对应的圆的圆周上,所形成的对应的两个圆为同心圆,其中:A为大于1的自然数,且两个圆柱体前端面的A的取值不要求相同;每一个圆形图案内均设置有光源,所述第一圆柱体1-1和第二圆柱体1-2的前端面上均贴合安装有挡光板1-4,每一个挡光板1-4上均对应开设有圆孔,圆孔的位置与对应的圆形图案相一致,且圆孔的直径小于圆形图案的直径。挡光板对光源进行遮挡,防止光线衍射带来误差。
在待装配孔内的安装阶梯靶标,所述阶梯靶标1上的图案朝向待装配孔前端,且阶梯靶标轴线与待装配孔处于共轴状态,在待装配孔前端朝向孔内采集阶梯靶标,得阶梯靶标图像。
步骤b.对步骤a中获得到阶梯靶标图像进行处理,得到边缘平滑过渡的阶梯靶标图像;具体包括如下步骤:
步骤b1.对步骤a中获得到阶梯靶标图像进行二值化处理,得到二值化阶梯靶标图像;
步骤b2.对b1中获得的二值化阶梯靶标图像进行滤波处理,得到边缘平滑过渡的二值化阶梯靶标图像。步骤b2中的具体处理过程如下:使用n*n的像素矩形沿着二值化阶梯靶标图像水平向右和竖直向下两个方向移动,将中心像素灰度值替换为相邻n*n像素矩形框内像素灰度的平均值,得到边缘平滑过渡的二值化阶梯靶标图像;其中n为大于0的自然数。
步骤c.对步骤b中边缘平滑过渡的二值化阶梯靶标图像采用Canny算子,获得阶梯靶标图像边缘轮廓集合,对边缘轮廓集合中的每个子集分别进行最小二乘椭圆拟合,得到一组拟合椭圆中心点集合;
步骤d.将步骤c中得到的拟合椭圆中心点集合分为第一圆柱体前端面点集和第二圆柱体前端面点集;具体过程如下:
步骤d1.设处理后的阶梯靶标图像左上角第一个像素点为图像坐标系原点,水平向右为X轴正方向,竖直向下为Y轴正方向;在图像坐标系中,将所述拟合椭圆中心点按照距离图像坐标系原点由小到大的原则重新排序,获得有序点集Si(i=1,2,3,4,5,6…);
步骤d2.遍历有序点集Si,搜索满足正方形分布的坐标点,首先设S1为正方形的一个顶点Pt1,以顶点Pt1为对角线的一个点,遍历有序点集Si,在该对角线上寻找正方形的第二个顶点Pt2,该对角线的斜率k满足如下要求:0.5<k<3;
步骤d3.搜索正方形的第三个顶点,设Pt1和Pt2两点间对角线对应的一次函数为:y1=ax+b,该对角线中点为其中:Pt1(m1,n1)和Pt2(m2,n2);Pt1到M点的距离为L,过M点作Pt1与Pt2连线的中垂线,得到该中垂线的一次函数为:设中垂线上一点为:遍历有序点集Si,若存在Si满足则将Si设为Pt3点;
若Si不满足重新执行步骤d2,搜寻其他满足条件的Pt1、Pt2再执行步骤d3,直至存在Si满足则将Si设为Pt3点;
步骤d4.搜寻正方形的第四个顶点,以点Pt1,Pt3确定向量VecX,点Pt2,Pt3确定向量VecY,设Pt1为坐标原点,向量VecX所在方向为X轴,向量VecY所在方向为Y轴,构建新坐标系;由正方形顶点排列特性得到参考点Pt4’在新构建的坐标系中可表示为:Pt4’=Pt1+VecY,遍历有序点集Si,保存满足的Si点集,设为点集Qi,点集Qi中满足||Pt4'-Qi||取得最小值的点即为Pt4点,若存在Pt4点则执行步骤d6,若不存在Pt4点,则执行步骤d5;
步骤d5.如果以S1点作为第一个顶点Pt1,不存在另外3个点满足步骤d2、d3、d4,则遍历点集Si,依次令Si(i=2,3,4…)为Pt1,再继续步骤d2、d3、d4;
步骤d6.将搜寻到的满足正方形分布的4个坐标点作为采集的阶梯靶标图像的参照点,计算椭圆圆心Si(i=1,2,3…)到参照点的距离Gi(i=1,2,3…),并求得所有椭圆圆心到参照点的距离和的均值Ver;
若Gi>Ver,且则该点属于外圆;
若Gi<Ver且则该点属于内圆;
其中比例值P为阶梯靶标的内外圆半径之比;
得到内圆点集和外圆点集,内圆点集为第一圆柱体1-1前端面点集内圆点集为第二圆柱体1-2前端面点集;
其中:内圆为第一圆柱体1-1前端面上的圆形图案的圆心形成的圆,外圆为第二圆柱体1-2前端面上的前端面上的圆形图案的圆心形成的圆;
步骤d7.若阶梯靶标倾斜角度过大,前端面对后端面的圆形图案遮挡,经过分类处理后只能得到第一圆柱体1-1前端面点集,则再次遍历点集Si剔除得到法人第一圆柱体1-1前端面点集,即可得到第二圆柱体1-2前端面点集。
步骤e.将步骤d中的所述第一圆柱体前端面点集和第二圆柱体前端面点集,分别运用最小二乘法进行椭圆拟合,即得到对应的第一圆柱体前端面中心点和第二圆柱体前端面中心点;连接两个中心点,即为阶梯靶标1的轴线,靶标轴线的位姿即为待装配孔的位姿。
将安装于待装配孔内的阶梯靶标拆除,选择合适的另一个阶梯靶标,装配到将与待装配孔配装的待装配轴上,另一个阶梯靶标的轴线与待装配轴的轴线相重合,确定另一个阶梯靶标的轴线,该另一个阶梯靶标的轴线的位姿即为待装配轴的位姿。本发明由“西北工业大学研究生创意创新种子基金资助。”
实施例
步骤a.用工业相机采集阶梯靶标1,得阶梯靶标图像,如图2所示;
步骤b.对阶梯靶标图像进行二值化处理,运用最大类间方差法计算得到最佳分割阈值μ,μ=124,得二值化阶梯靶标图像,如图3所示;
对二值化阶梯靶标图像进行处理,采用均值滤波器对以当前像素点(x,y)为中心的的5×5邻域内的所有像素灰度求均值,将均值作为当前像素的新灰度值,得到边缘平滑过渡的二值化阶梯靶标图像;
步骤c.对步骤b中获得的边缘平滑过渡的二值化阶梯靶标图像处理,采用Canny算子获得阶梯靶标图像边缘轮廓集合,如图4所示,对边缘轮廓集合中的每个子集分别进行最小二乘椭圆拟合,得到一组拟合椭圆中心点集合,椭圆中心点坐标为(单位:像素):
(2845.97,3109.94)、(3553.48,2817.48)、(2138.67,2816.52)、(2845.99,2554.95)
(3161.02,2425.02)、(2531.10,2424.04)、(3846.93,2110.11)、(3291.93,2110.11)
(2401.15,2109.01)、(1846.15,2109.01)、(3162.02,1795.06)、(2532.03,1794.1)
(2847.08,1664.18)、(3554.47,1402.63)、(2139.59,1401.66)、(2847.09,1109.18);
步骤d.根据所使用阶梯靶标的特点,使用一种基于正方形分布阵列的方法对对所述拟合椭圆中心点集合分类为第一圆柱体端面点集和二圆柱体端面点集,如图5所示,具体步骤如下:
步骤d1.设边缘平滑过渡的二值化阶梯靶标图像左上角第一个像素点为图像坐标系原点,水平向右为X轴正方向,竖直向下为Y轴正方向;设图像坐标原点O为图像左上角顶点,将所述拟合椭圆中心点坐标按照距离图像坐标系原点O由小到大的原则重新排序,获得有序点集Si(i=1,2,3,4,5,6…);
步骤d2.遍历点集Si,搜索满足正方形分布的坐标点,首先设S1为正方形的一个顶点Pt1,以顶点Pt1为对角线的一个点,遍历点集Si,在该对角线上寻找正方形的第二个顶点Pt2,考虑到靶标位姿变化范围较小的因素,该对角线的斜率k满足如下要求:0.8<k<1.25;遍历点集Si如果存在满足条件的中心点则设为Pt2,否则重新选择Pt1点;
步骤d3.搜寻正方形第三个顶点,设Pt1和Pt2两点间对角线对应的一次函数:y1=ax+b,该对角线中点其中:Pt1(m1,n1)和Pt2(m2,n2);Pt1到M点的距离为L,L=1000.46,过M点作Pt1与Pt2连线的中垂线,得到该重垂线的一次函数为设中垂线上一点遍历点集Si,若存在Si满足则将Si设为Pt3点;
若Si不满足重新执行步骤e2,搜寻其他满足条件的Pt1、Pt2再执行步骤e3,直至存在Si满足则将Si设为Pt3点;
步骤d4.搜寻正方形的第四个顶点,以点Pt1,Pt3确定向量VecX,点Pt2,Pt3确定向量VecY,设Pt1为坐标原点,向量VecX所在方向为X轴,向量VecY所在方向为Y轴,建立坐标系;由正方形顶点排列特性得到参考点Pt4’在新构建的坐标系中可表示为:Pt4’=Pt1+VecY,遍历点集Si,保存满足的Si点集,设为点集Qi,点集Qi中满足||Pt4'-Qi||取得最小值的点即为Pt4点,若存在Pt4点则执行步骤e6,若不存在Pt4点,则执行步骤e5;
步骤d5.如果以S1点作为第一个顶点Pt1,不存在另外3个点满足步骤e2、e3、e4,则遍历点集Si,依次令Si(i=2,3,4…)为Pt1,再继续步骤e2、e3、e4;如果没有满足要求的4个点,则识别失败退出;
步骤d6.将搜寻到满足正方形分布的4个坐标点的中心作为整个模式分类中心的初步位置,然后计算椭圆圆心Si(i=1,2,3…)到模式分类中心的距离Gi(i=1,2,3…),并求得所有椭圆圆心到模式分类中心的距离和的均值Ver=722.9;
若Gi>Ver且则该点属于外圆;
若Gi<Ver且则该点属于内圆;
其中比例值P为阶梯靶标的内外圆半径之比;
比例值为阶梯靶标实际的内外圆半径之比;
得到内圆点集和外圆点集,内圆点集为第一圆柱体1-1前端面点集内圆点集为第二圆柱体1-2前端面点集;
其中:内圆为第一圆柱体1-1上的圆形图案的圆心形成的圆,外圆为第二圆柱体1-2上的圆形图案的圆心形成的圆。
步骤e.将内外圆的两类点集再次运用最小二乘法进行椭圆拟合即可得到第一圆柱体前端面靶标和第二圆柱体前端面靶标的中心,如图6所示。第一圆柱体前端面和第二圆柱体前端面靶标中心点坐标分别为:(2847.1,1109.18)、(2846.54,2109.56)。连接两个中心点,即为阶梯靶标1的轴线,靶标轴线的位姿即为待装配孔的位姿。
Claims (7)
1.一种轴孔装配自动对中的靶心识别方法,其特征在于,该方法包括以下步骤:
步骤a.使用阶梯靶标(1),所述阶梯靶标(1)包括同轴设置且位于前部的第一圆柱体(1-1)和位于后部的第二圆柱体(1-2),且第一圆柱体(1-1)的直径小于第二圆柱体(1-2)的直径,所述第一圆柱体(1-1)和第二圆柱体(1-2)的前端面上、且绕其边缘一周均设置有4A个圆形图案(1-3),且对应端面上的圆形图案(1-3)的圆心分别位于对应的圆的圆周上,所形成的对应的两个圆为同心圆,其中:A为大于1的自然数,且两个圆柱体前端面的A的取值不要求相同;
在待装配孔内的安装阶梯靶标,所述阶梯靶标(1)上的图案朝向待装配孔前端,且阶梯靶标轴线与待装配孔处于共轴状态,在待装配孔前端朝向孔内采集阶梯靶标,得阶梯靶标图像;
步骤b.对步骤a中获得到阶梯靶标图像进行处理,得到边缘平滑过渡的阶梯靶标图像;
步骤c.对步骤b中边缘平滑过渡的阶梯靶标图像处理,获得阶梯靶标图像边缘轮廓集合,对边缘轮廓集合中的每个子集分别进行最小二乘椭圆拟合,得到一组拟合椭圆中心点集合;
步骤d.将步骤c中得到的拟合椭圆中心点集合分为第一圆柱体前端面点集和第二圆柱体前端面点集;
步骤e.将步骤d中的所述第一圆柱体前端面点集和第二圆柱体前端面点集,分别运用最小二乘法进行椭圆拟合,即得到对应的第一圆柱体前端面中心点和第二圆柱体前端面中心点;连接两个中心点,即为阶梯靶标(1)的轴线,靶标轴线的位姿即为待装配孔的位姿。
2.按照权利要求1所述的一种轴孔装配自动对中的靶心识别方法,其特征在于,所述步骤d的具体过程如下:
步骤d1.设处理后的阶梯靶标图像左上角第一个像素点为图像坐标系原点,水平向右为X轴正方向,竖直向下为Y轴正方向;在图像坐标系中,将所述拟合椭圆中心点按照距离图像坐标系原点由小到大的原则重新排序,获得有序点集Si(i=1,2,3,4,5,6…);
步骤d2.遍历有序点集Si,搜索满足正方形分布的坐标点,首先设S1为正方形的一个顶点Pt1,以顶点Pt1为对角线的一个点,遍历有序点集Si,在该对角线上寻找正方形的第二个顶点Pt2,该对角线的斜率k满足如下要求:0.5<k<3;
步骤d3.搜索正方形的第三个顶点,设Pt1和Pt2两点间对角线对应的一次函数为:y1=ax+b,该对角线中点为其中:Pt1(m1,n1)和Pt2(m2,n2);Pt1到M点的距离为L,过M点作Pt1与Pt2连线的中垂线,得到该中垂线的一次函数为:设中垂线上一点为:遍历有序点集Si,若存在Si满足则将Si设为Pt3点;
若Si不满足重新执行步骤d2,搜寻其他满足条件的Pt1、Pt2再执行步骤d3,直至存在Si满足则将Si设为Pt3点;
步骤d4.搜寻正方形的第四个顶点,以点Pt1,Pt3确定向量VecX,点Pt2,Pt3确定向量VecY,设Pt1为坐标原点,向量VecX所在方向为X轴,向量VecY所在方向为Y轴,构建新坐标系;由正方形顶点排列特性得到参考点Pt4’在新构建的坐标系中可表示为:Pt4’=Pt1+VecY,遍历有序点集Si,保存满足的Si点集,设为点集Qi,点集Qi中满足||Pt4'-Qi||取得最小值的点即为Pt4点,若存在Pt4点则执行步骤d6,若不存在Pt4点,则执行步骤d5;
步骤d5.如果以S1点作为第一个顶点Pt1,不存在另外3个点满足步骤d2、d3、d4,则遍历点集Si,依次令Si(i=2,3,4…)为Pt1,再继续步骤d2、d3、d4;
步骤d6.将搜寻到的满足正方形分布的4个坐标点作为采集的阶梯靶标图像的参照点,计算椭圆圆心Si(i=1,2,3…)到参照点的距离Gi(i=1,2,3…),并求得所有椭圆圆心到参照点的距离和的均值Ver;
若Gi>Ver,且则该点属于外圆;
若Gi<Ver且则该点属于内圆;
其中比例值P为阶梯靶标的内外圆半径之比;
得到内圆点集和外圆点集,内圆点集为第一圆柱体(1-1)前端面点集内圆点集为第二圆柱体(1-2)前端面点集;
其中:内圆为第一圆柱体(1-1)前端面上的圆形图案的圆心形成的圆,外圆为第二圆柱体(1-2)前端面上的前端面上的圆形图案的圆心形成的圆;
步骤d7.若阶梯靶标倾斜角度过大,前端面对后端面的圆形图案遮挡,经过分类处理后只能得到第一圆柱体(1-1)前端面点集,则再次遍历点集Si剔除得到法人第一圆柱体(1-1)前端面点集,即可得到第二圆柱体(1-2)前端面点集。
3.按照权利要求1或2所述的一种轴孔装配自动对中的靶心识别方法,其特征在于,所述步骤b的具体过程如下:
步骤b1.对步骤a中获得到阶梯靶标图像进行二值化处理,得到二值化阶梯靶标图像;
步骤b2.对b1中获得的二值化阶梯靶标图像进行滤波处理,得到边缘平滑过渡的二值化阶梯靶标图像。
4.按照权利要求3所述的一种轴孔装配自动对中的靶心识别方法,其特征在于,步骤b1中采用最大类间方差法计算阈值,以得到二值化阶梯靶标图像。
5.按照权利要求3所述的一种轴孔装配自动对中的靶心识别方法,其特征在于,所述步骤b2中的具体处理过程如下:使用n*n的像素矩形沿着二值化阶梯靶标图像水平向右和竖直向下两个方向移动,将中心像素灰度值替换为相邻n*n像素矩形框内像素灰度的平均值,得到边缘平滑过渡的二值化阶梯靶标图像;其中n为大于0的自然数。
6.按照权利要求5所述的一种轴孔装配自动对中的靶心识别方法,其特征在于,所述步骤c中,将边缘平滑过渡的二值化阶梯靶标图像采用Canny算子获得阶梯靶标图像边缘轮廓集合。
7.按照权利要求1-6中任一项所述的一种轴孔装配自动对中的靶心识别方法,其特征在于,每一个圆形图案内均设置有光源,所述第一圆柱体(1-1)和第二圆柱体(1-2)的前端面上均贴合安装有挡光板(1-4),每一个挡光板(1-4)上均对应开设有圆孔,圆孔的位置与对应的圆形图案相一致,且圆孔的直径小于圆形图案的直径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710623740.4A CN107633502B (zh) | 2017-07-27 | 2017-07-27 | 一种轴孔装配自动对中的靶心识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710623740.4A CN107633502B (zh) | 2017-07-27 | 2017-07-27 | 一种轴孔装配自动对中的靶心识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107633502A true CN107633502A (zh) | 2018-01-26 |
CN107633502B CN107633502B (zh) | 2020-09-29 |
Family
ID=61099595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710623740.4A Active CN107633502B (zh) | 2017-07-27 | 2017-07-27 | 一种轴孔装配自动对中的靶心识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107633502B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109048903A (zh) * | 2018-08-17 | 2018-12-21 | 上海非夕机器人科技有限公司 | 机械臂插轴入孔的方法及设备 |
CN111085837A (zh) * | 2019-12-28 | 2020-05-01 | 成都行必果光电科技有限公司 | 一种自动化装配位姿态测量方法 |
CN111322967A (zh) * | 2020-03-04 | 2020-06-23 | 西北工业大学 | 一种面向阶梯轴与孔装配过程的对中方法 |
CN112634345A (zh) * | 2020-12-18 | 2021-04-09 | 福州华虹智能科技股份有限公司 | 基于最小二乘法的图像特征信息提取矿井壁图像的方法 |
CN113111687A (zh) * | 2020-01-13 | 2021-07-13 | 阿里巴巴集团控股有限公司 | 数据处理方法、系统及电子设备 |
CN116930192A (zh) * | 2023-07-26 | 2023-10-24 | 常州润来科技有限公司 | 一种高精度铜管缺陷检测方法及系统 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101699217A (zh) * | 2009-11-03 | 2010-04-28 | 武汉大学 | 一种用于工业零件的同心圆检测方法 |
CN102062577A (zh) * | 2010-12-06 | 2011-05-18 | 柳州科路测量仪器有限责任公司 | 激光轮对自动检测机及其用于轮对自动检测的方法 |
CN102218652A (zh) * | 2011-05-09 | 2011-10-19 | 无锡引域智能机器人有限公司 | 一种利用机器人的柔顺性实现轴孔装配的装置和方法 |
US20130163849A1 (en) * | 2010-09-14 | 2013-06-27 | Ronny Jahnke | Apparatus and method for automatic inspection of through-holes of a component |
CN203981114U (zh) * | 2014-08-05 | 2014-12-03 | 吉林大学 | 基于球面坐标的汽车视觉检测系统的摄像机标定靶标 |
CN104424601A (zh) * | 2013-08-28 | 2015-03-18 | 深圳市智信精密仪器有限公司 | 异形体装配件的对中装配方法及装置 |
CN104568295A (zh) * | 2014-12-08 | 2015-04-29 | 太原航空仪表有限公司 | 嵌入式大气数据系统压力故障监控处理方法 |
CN105320940A (zh) * | 2015-10-08 | 2016-02-10 | 江苏科技大学 | 一种基于中心-边界连接模型的交通标志检测方法 |
US9367908B1 (en) * | 2011-04-21 | 2016-06-14 | Ivan Bekey | Overcoming inappropriate apparent rotation of vehicle wheels in videos by image modification |
CN106204610A (zh) * | 2016-07-20 | 2016-12-07 | 张家港长安大学汽车工程研究院 | 基于图像处理的轴孔装配自动对中系统 |
CN106651989A (zh) * | 2016-12-21 | 2017-05-10 | 北京信息科技大学 | 用于微小孔对中的嵌套圆拟合方法及装置 |
-
2017
- 2017-07-27 CN CN201710623740.4A patent/CN107633502B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101699217A (zh) * | 2009-11-03 | 2010-04-28 | 武汉大学 | 一种用于工业零件的同心圆检测方法 |
US20130163849A1 (en) * | 2010-09-14 | 2013-06-27 | Ronny Jahnke | Apparatus and method for automatic inspection of through-holes of a component |
CN102062577A (zh) * | 2010-12-06 | 2011-05-18 | 柳州科路测量仪器有限责任公司 | 激光轮对自动检测机及其用于轮对自动检测的方法 |
US9367908B1 (en) * | 2011-04-21 | 2016-06-14 | Ivan Bekey | Overcoming inappropriate apparent rotation of vehicle wheels in videos by image modification |
CN102218652A (zh) * | 2011-05-09 | 2011-10-19 | 无锡引域智能机器人有限公司 | 一种利用机器人的柔顺性实现轴孔装配的装置和方法 |
CN104424601A (zh) * | 2013-08-28 | 2015-03-18 | 深圳市智信精密仪器有限公司 | 异形体装配件的对中装配方法及装置 |
CN203981114U (zh) * | 2014-08-05 | 2014-12-03 | 吉林大学 | 基于球面坐标的汽车视觉检测系统的摄像机标定靶标 |
CN104568295A (zh) * | 2014-12-08 | 2015-04-29 | 太原航空仪表有限公司 | 嵌入式大气数据系统压力故障监控处理方法 |
CN105320940A (zh) * | 2015-10-08 | 2016-02-10 | 江苏科技大学 | 一种基于中心-边界连接模型的交通标志检测方法 |
CN106204610A (zh) * | 2016-07-20 | 2016-12-07 | 张家港长安大学汽车工程研究院 | 基于图像处理的轴孔装配自动对中系统 |
CN106651989A (zh) * | 2016-12-21 | 2017-05-10 | 北京信息科技大学 | 用于微小孔对中的嵌套圆拟合方法及装置 |
Non-Patent Citations (3)
Title |
---|
WEN-HAN QIAN 等: "On the Optimization of Automobile Panel Fitting", 《PROCEEDINGS OF IEEE INTERNATIONAL CONFERENCE ON ROBOTICS AND AUTOMATION》 * |
朱宁宁 等: "一种平面靶心的提取算法", 《大地测量与地球动力学》 * |
程云勇 等: "基于误差控制的薄壁叶片椭圆弧形前后缘建模方法", 《计算机辅助设计与图形学学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109048903A (zh) * | 2018-08-17 | 2018-12-21 | 上海非夕机器人科技有限公司 | 机械臂插轴入孔的方法及设备 |
CN109048903B (zh) * | 2018-08-17 | 2021-08-27 | 上海非夕机器人科技有限公司 | 机械臂插轴入孔的方法及设备 |
CN111085837A (zh) * | 2019-12-28 | 2020-05-01 | 成都行必果光电科技有限公司 | 一种自动化装配位姿态测量方法 |
CN113111687A (zh) * | 2020-01-13 | 2021-07-13 | 阿里巴巴集团控股有限公司 | 数据处理方法、系统及电子设备 |
CN111322967A (zh) * | 2020-03-04 | 2020-06-23 | 西北工业大学 | 一种面向阶梯轴与孔装配过程的对中方法 |
CN111322967B (zh) * | 2020-03-04 | 2021-05-07 | 西北工业大学 | 一种面向阶梯轴与孔装配过程的对中方法 |
CN112634345A (zh) * | 2020-12-18 | 2021-04-09 | 福州华虹智能科技股份有限公司 | 基于最小二乘法的图像特征信息提取矿井壁图像的方法 |
CN112634345B (zh) * | 2020-12-18 | 2023-03-31 | 福州华虹智能科技股份有限公司 | 基于最小二乘法的图像特征信息提取矿井壁图像的方法 |
CN116930192A (zh) * | 2023-07-26 | 2023-10-24 | 常州润来科技有限公司 | 一种高精度铜管缺陷检测方法及系统 |
CN116930192B (zh) * | 2023-07-26 | 2024-02-09 | 常州润来科技有限公司 | 一种高精度铜管缺陷检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107633502B (zh) | 2020-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107633502A (zh) | 一种轴孔装配自动对中的靶心识别方法 | |
CN106251353A (zh) | 弱纹理工件及其三维位姿的识别检测方法及系统 | |
CN101464999B (zh) | 三维色阶比对系统及方法 | |
CN105354880B (zh) | 一种基于线激光扫描的喷砂机器人自动路径生成方法 | |
CN107516098B (zh) | 一种基于边缘曲率角的目标轮廓三维信息提取方法 | |
CN107617876B (zh) | 一种基于机器视觉监测的装配对中方法 | |
CN106778541B (zh) | 一种基于视觉的多层梁中外梁孔的识别与定位方法 | |
CN107622499A (zh) | 一种基于目标二维轮廓模型的识别与空间定位方法 | |
CN106204610A (zh) | 基于图像处理的轴孔装配自动对中系统 | |
CN109345554A (zh) | 一种基于rgb-d相机的粘连蘑菇视觉原位测量方法 | |
CN103530889A (zh) | 一种基于骨架提取的测量棒编码标记点解码方法 | |
CN109064481B (zh) | 一种机器视觉定位方法 | |
CN110223345A (zh) | 基于点云的配电线路作业对象位姿估计方法 | |
CN106023186B (zh) | 一种基于机器视觉的毛坯轮毂钻孔中心精确定位方法 | |
CN107578430A (zh) | 一种基于自适应权值和局部熵的立体匹配方法 | |
CN108596925A (zh) | 基于视觉的太阳花模组表面螺钉孔位置图像处理方法 | |
CN109544513A (zh) | 一种钢管端面缺陷提取识别的方法 | |
CN105825515B (zh) | 一种用于自主空中加油的加油锥套图像检测定位方法 | |
CN108088444A (zh) | 基于三维激光的室内点云地图生成系统及方法 | |
CN107346437A (zh) | 车身侧视图参数化模型的自动提取方法 | |
CN104992448B (zh) | 机器人防损伤采摘葡萄的自动定位方法 | |
CN109359604A (zh) | 面向巡检机器人的阴影干扰下仪表识别方法 | |
CN104599267A (zh) | 切割轨迹生成方法及装置 | |
CN110288619A (zh) | 基于视觉的太阳花模组表面螺钉孔位置检测方法 | |
CN111815575B (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 |