基于视觉的薄板对接焊缝实时检测方法技术

本发明专利技术属于检测技术领域，特别涉及一种基于视觉的薄板对接焊缝实时检测方法。首先对焊缝图像进行列扫描和横扫描，得到焊缝的稀疏特征点列，再对特征点列进行基于粒子滤波的二次曲线拟合，实现对薄板对接焊缝的特征曲线参数的检测。由于采用稀疏点列，大大提高了检测速度和跟踪性能，由于采用的粒子滤波算法，不仅可以减小计算量，还可以有效地避免点固点等干扰的影响并提高了检测精度，其实时性能和抗干扰性能都达到了实际焊接的要求，检测实时性能达到了３ｍｓ／帧，检测精度达到了０．０５ｍｍ。

【技术实现步骤摘要】

本专利技术属于检测
，特别涉及一种。
技术介绍
薄板对接焊缝具有较窄的焊缝(1mm以下)，较高的焊接速度(1m/min以上)，采用常规的结构光焊缝跟踪难以检测如此小的焊缝，采用常规的视觉检测技术如边缘检测算子，小波变换等，因为消耗较多的计算用时，难以达到薄板对接焊缝的实时速度，并且点固点等干扰对这些算法的影响较大，检测正确性和精度都难以保障。
技术实现思路
本专利技术提供了一种，其特征在于，包括以下步骤 (1)使用工业微型CCD对薄板对接焊缝进行检测，扫描得到焊缝特征点列根据检测得到焊缝图像，进行横扫描，得到焊缝横扫描的灰度值，根据灰度突变的最大的地方，作为焊缝特征点，最后得到整幅图像的焊缝特征点列，该特征点列的数目远小于整幅图像的数据，因此可以大大减少运算量，提高运算速度； (2)设图像数据点的采样间隔为T，焊缝点列特征点沿垂直方向的坐标为tk，k为计数标志，焊缝特征点列的沿着垂直轴线的位置、速度、加速度和加加速度分别为sk，vk，ak，jk，则 其中jk取值为白噪声随机变量，且相邻两个时间间隔的焊缝图像位移 (3)设特征点在tk位置的状态变量为Xk＝T，系统状态转移矩阵Φk，k-1，系统噪声矩阵Γk-1和量测矩阵Hk可阵描述为 在此设定下，根据粒子数确定下一个时刻Xk的粒子，特征点列的参数估计系统可描述为 i＝1，2…N Zk＝HkXk+Vk 其中，Wk-1i为均方差为1的高斯白噪声，N为粒子数，R为粒子半径，Vk为量测误差且Vk＝δ·sk＝-KΔs，其中[]代表图像象素取整； (4)设定各个粒子的初始权重为i＝1，2…N； (5)确定观测量及权重更新 进行归一化处理其中，σ为常数， 为k-1时刻的焊缝模板在(HXkHi，HXkVi)的位置上与k时刻焊缝图像相关值； 粒子权重更新得到 (6)对系统进行状态估计 (7)粒子重采样 重采样是对粒子权重为最小的粒子，从权值大的粒子上衍生出一些粒子来代替该粒子，其权值重新设置为1/N； (8)i＝i+1返回步骤(6)重新确定状态Xki，直到Xki和上次状态估计的距离时，获得第次的状态估计Xk，采用其状态估计对焊缝采用非遍历递归最小二乘算法进行二次曲线拟合，得到由焊缝特征曲线参数即可实现对薄板对接焊缝进行跟踪； (9)k＝k+1返回步骤(5)进行下一时刻的状态估计和跟踪。 所述焊缝特征点列通过对焊缝图像进行一维横扫描或列扫描得到。 本专利技术的有益效果为采用了焊缝稀疏特征点列的非遍历递归最小二乘算法，大大提高了检测速度和跟踪性能；采用的粒子滤波算法，不仅可以减小计算量，还可以有效地避免点固点等干扰的影响并提高了检测精度，其实时性能和抗干扰性能都达到了实际焊接的要求，检测实时性能达到了3ms/帧，检测精度达到了0.05mm。 附图说明 图1为集中箱的焊缝照片； 图2为图像横扫描和特征检测结果； 图3为整幅图像的特征点列检测结果； 图4为焊缝特征参数检测流程； 图5为焊缝曲线检测结果。 具体实施例方式 本专利技术提供了一种，下面通过附图说明和具体实施方式对本专利技术做进一步说明。 本专利技术使用中国专利CN200920105035.6所述的装置对薄板对接焊缝进行检测。 根据检测得到焊缝图像，进行横扫描，得到焊缝横扫描的灰度值，根据灰度突变的最大的地方，作为焊缝特征点，如图2所示，最后得到整幅图像的焊缝特征点列，如图3所示，该特征点列的数目远小于整幅图像的数据，因此可以大大减少运算量，提高运算速度。 粒子滤波(particle filter)是求解贝叶斯概率的一种实用算法，它通过非参数化的蒙特卡罗模拟方法来实现递推贝叶斯滤波，适用于任何能用状态空间模型表示的非线性系统，以及传统卡尔曼滤波无法表示的非线性系统，其方法灵活，计算量少，速度快，容易实现等优点。 设图像数据点的采样间隔为T，焊缝点列特征点沿垂直方向的坐标为tk，焊缝特征点列的沿着垂直轴线的位置、速度、加速度和加加速度分别为sk，vk，ak，jk。 设jk为白噪声随机变量，且相邻两个时间间隔的焊缝图像位移为 设特征点在tk位置的状态变量为Xk＝T 系统状态转移矩阵Φk，k-1，系统噪声矩阵Γk-1和量测矩阵Hk可阵描述为 在此设定下，特征点列的参数估计系统可描述为 Xk＝Φk，k-1Xk-1+Γk-1Wk-1(4) Zk＝HkXk+Vk 其量测误差为Vk＝δsk＝-KΔs (5) 其中 代表图像象素取整。 采用粒子滤波进行焊缝跟踪流程如下 1)系统初始化 设定计数标志k＝0，设定粒子数N＝512，和粒子半径R＝10，确定系统初始状态变量X0，其中s0，H，s0，V为初始位置，v0，H，v0，V，a0，H，a0，V为零； 2)状态转移 根据粒子数确定下一个时刻Xk的粒子 i＝1，2…N(6) 其中Wk-1i为均方差为1的高斯白噪声； 设定各个粒子的权重为i＝1，2…N； 3)确定观测量及权重更新 归一化处理得到 σ为常数，选取10； 为k-1时刻的焊缝模板在(HXkHi，HXkVi)的位置上与k时刻焊缝图像相关值； 粒子权重更新 4)状态估计 基于粒子滤波的焊缝跟踪的算法流程如图4所示； 5)粒子重采样 重采样即是在出现某些粒子权值太小时，从权值大的粒子上衍生出一些粒子来代替，其权值重新设置为1/N； 6)i＝i+1，返回步骤(4)重新确定状态Xki直到Xki和上次状态估计的距离时，获得第次的状态估计Xk，采用其状态估计对焊缝采用非遍历递归最小二乘算法进行二次曲线拟合，得到 焊缝特征曲线检测结果如图5所示，由焊缝特征曲线参数即可实现对薄板对接焊缝进行跟踪。 7)k＝k+1，返回步骤(3)进行下一时刻的状态估计和跟踪。权利要求1.，其特征在于，包括以下步骤(1)使用工业微型CCD对薄板对接焊缝进行检测，扫描得到焊缝特征点列根据检测得到焊缝图像，进行横扫描，得到焊缝横扫描的灰度值，根据灰度突变的最大的地方，作为焊缝特征点，最后得到整幅图像的焊缝特征点列，该特征点列的数目远小于整幅图像的数据，因此可以大大减少运算量，提高运算速度；(2)设图像数据点的采样间隔为T，焊缝点列特征点沿垂直方向的坐标为tk，k为计数标志，焊缝特征点列的沿着垂直轴线的位置、速度、加速度和加加速度分别为sk，vk，ak，jk，则其中jk取值为白噪声随机变量，且相邻两个时间间隔的焊缝图像位移(3)设特征点在tk位置的状态变量为Xk＝T，系统状态转移矩阵Φk，k-1，系统噪声矩阵Γk-1和量测矩阵Hk可阵描述为在此设定下，根据粒子数确定下一个时刻Xk的粒子，特征点列的参数估计系统可描述为i＝1，2...NZk＝HkXk+Vk其中，Wk-1i为均方差为1的高斯白噪声，N为粒子数，R为粒子半径，Vk为量测误差且Vk＝δ·sk＝-KΔs，其中[]代表图像象素取整；(4)设定各个粒子的初始权重为i＝1，2…N；(5)确定观测量及权重更新进行归一化处理其中，σ为常数本文档来自技高网...

【技术保护点】
基于视觉的薄板对接焊缝实时检测方法，其特征在于，包括以下步骤：　（１）使用工业微型ＣＣＤ对薄板对接焊缝进行检测，扫描得到焊缝特征点列根据检测得到焊缝图像，进行横扫描，得到焊缝横扫描的灰度值，根据灰度突变的最大的地方，作为焊缝特征点，最 后得到整幅图像的焊缝特征点列，该特征点列的数目远小于整幅图像的数据，因此可以大大减少运算量，提高运算速度；　（２）设图像数据点的采样间隔为Ｔ，焊缝点列特征点沿垂直方向的坐标为ｔ↓［ｋ］，ｋ为计数标志，焊缝特征点列的沿着垂直轴线的位置、 速度、加速度和加加速度分别为：ｓ↓［ｋ］，ｖ↓［ｋ］，ａ↓［ｋ］，ｊ↓［ｋ］，则　＊＊＊　其中：ｊ↓［ｋ］取值为白噪声随机变量，且相邻两个时间间隔的焊缝图像位移Ｚ↓［ｋ］＝ｓ↓［ｋ，Ｈ］＝Ｋ（ｖ↓［ｋ－１］Ｔ＋ａ↓［ｋ－１］Ｔ↑ ［２］／２）；　（３）设特征点在ｔ↓［ｋ］位置的状态变量为Ｘ↓［ｋ］＝［ｓ↓［ｋ］，ｖ↓［ｋ］，ａ↓［ｋ］］↑［Ｔ］，系统状态转移矩阵Φ↓［ｋ，ｋ－１］，系统噪声矩阵Γ↓［ｋ－１］和量测矩阵Ｈ↓［ｋ］可阵描述为：　＊＊＊　 在此设定下，根据粒子数确定下一个时刻Ｘ↓［ｋ］的粒子，特征点列的参数估计系统可描述为：　Ｘ↓［ｋ］↑［ｉ］＝Φ↓［ｋ，ｋ－１］Ｘ↓［ｋ－１］↑［ｉ］＋Γ↓［ｋ－１］ＲＷ↓［ｋ－１］↑［ｉ］　ｉ＝１，２…Ｎ　Ｚ↓［ｋ］＝Ｈ↓［ｋ ］Ｘ↓［ｋ］＋Ｖ↓［ｋ］　其中，Ｗ↓［ｋ－１］↑［ｉ］为均方差为１的高斯白噪声，Ｎ为粒子数，Ｒ为粒子半径，Ｖ↓［ｋ］为量测误差且Ｖ↓［ｋ］＝δ.ｓ↓［ｋ］＝［ＫΔｓ］－ＫΔｓ，其中［］代表图像象素取整；　（４）设定各个粒子的初始 权重为：ω↓［ｋ］↑［ｉ］＝１／Ｎ　ｉ＝１，２…Ｎ；　（５）确定观测量及权重更新：　Ｐ（Ｚ↓［ｋ］↑［ｉ］｜Ｘ↓［ｋ］↑［ｉ］）＝ｅｘｐ（－１／２σ↑［２］ｒ↓［ｌ↓［ｋ，ｋ－１］］（ＨＸ↓［ｋＨ］↑［ｉ］，ＨＸ↓［ｋＶ］↑［ ｉ］）），进行归一化处理：Ｐ（Ｚ↓［ｋ］｜Ｘ↓［ｋ］↑［ｉ］）＝Ｐ（Ｚ↓［ｋ］↑［ｉ］｜Ｘ↓［ｋ］↑［ｉ］）／∑Ｐ（Ｚ↓［ｋ］↑［ｉ］｜Ｘ↓［ｋ］↑［ｉ］），其中，σ为常数，ｒ↓［ｌ↓［ｋ，ｋ－１］］（ＨＸ↓［ｋＨ］↑［ｉ］，ＨＸ↓［ｋＶ］↑［ｉ］）为ｋ－１时刻的焊缝模板在（ＨＸ↓［ｋＨ］↑［ｉ］，ＨＸ↓［ｋＶ］↑［ｉ］）的位置上与ｋ时刻焊缝图像相关值；　粒子权重更...

【技术特征摘要】

【专利技术属性】
技术研发人员：郑军，潘际銮，
申请(专利权)人：清华大学，
类型：发明
国别省市：11[中国|北京]