发明名称 一种电子束焊接熔池形状参数的视觉检测方法
摘要 一种电子束焊接熔池形状参数的视觉检测方法,它涉及电子束焊接熔池形状参数的视觉检测方法,本发明要解决现有熔池图像处理中存在噪声干扰大,无法获得连续,光滑的电子束焊熔池边缘的问题,本发明对电子束焊彩色熔池图像视觉传感系统进行标定,然后启动系统采集熔池图像,采用二值形态学图像处理算法提取熔池图像边缘,最后利用熔池形状参数提取算法对熔池形状参数进行提取;本发明提取的熔池边缘图像完全消除了噪声的影响,边缘连续光滑,熔池形状参数提取精确度高,整个处理过程所需时间大约为30ms,完全可以满足实时检测的需要,本发明应用于焊接自动化领域。
申请公布号 CN102519387B 申请公布日期 2014.03.05
申请号 CN201110331374.8 申请日期 2011.10.27
申请人 哈尔滨工业大学 发明人 张秉刚;石铭霄;陈国庆;冯吉才
分类号 G01B11/24(2006.01)I;G01B11/03(2006.01)I;G06T5/00(2006.01)I 主分类号 G01B11/24(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 金永焕
主权项 1.一种电子束焊接熔池形状参数的视觉检测方法,其特征在于所述的电子束焊接熔池形状参数的视觉检测方法具体过程如下:步骤A、将待焊工件放入真空室内固定;步骤B、启动真空室内摄像机,对电子束焊彩色熔池图像视觉传感系统进行标定:一、以焊接工件表面的焊接起点为坐标原点o,焊接方向为x轴正方向,垂直于焊接方向向下为y轴正方向,建立坐标系oxy,即为工件坐标系;二、然后以工件坐标系为参照系画一组同心圆,其中最小圆的半径为1毫米,每个同心圆的半径以1毫米为单位递增,最小的圆与oy轴的交点为A点,最小的圆与ox轴的交点为B点;三、对步骤二所画的坐标系进行拍照;四、以照片平面的左上角为原点O,原点O正下方为Y轴正方向,与Y轴垂直向右为X轴正方向,建立坐标系OXY,即为图像坐标系;五、然后通过下列公式计算出图像坐标系中单位像素在工件坐标系中xy方向上的实际尺寸:<maths num="0001"><![CDATA[<math><mrow><msub><mi>w</mi><mi>x</mi></msub><mo>=</mo><mfrac><mn>1</mn><msqrt><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>B</mi></msub><mo>-</mo><msub><mi>X</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Y</mi><mi>B</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>w</mi><mi>y</mi></msub><mo>=</mo><mfrac><mn>1</mn><msqrt><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>A</mi></msub><mo>-</mo><msub><mi>X</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Y</mi><mi>A</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac></mrow></math>]]></maths>W<sub>x</sub>为熔池宽度;W<sub>y</sub>为熔池长度;X<sub>o</sub>为图像坐标系坐标原点;X<sub>B</sub>为图像坐标系坐标B点;X<sub>A</sub>为图像坐标系坐标A点;步骤C、提取熔池图像边缘:(1)采用二值形态学提取熔池图像边缘,具体过程如下:设熔池图像灰度值为0~255,设熔池图像灰度阈值为k,k分别取值1,2…,m,将比k大的熔池图像灰度值设为w<sub>1</sub>(k),平均灰度值为M<sub>1</sub>(k),方差为σ<sub>1</sub>(k);将比k小的熔池图像灰度值设为w<sub>2</sub>(k),平均灰度值为M<sub>2</sub>(k),方差为σ<sub>2</sub>(k);设熔池图像灰度值的总平均值为M<sub>T</sub>,则类内方差<img file="FDA0000416693880000013.GIF" wi="94" he="85" />类间方差<img file="FDA0000416693880000014.GIF" wi="56" he="69" />和灰度值的总平均值M<sub>T</sub>的数学表达式分别为:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>W</mi><mn>2</mn></msubsup><mo>=</mo><msub><mi>w</mi><mn>1</mn></msub><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>B</mi><mn>2</mn></msubsup><mo>=</mo><msub><mi>w</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>M</mi><mn>1</mn></msub><mo>-</mo><msub><mi>M</mi><mi>T</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>M</mi><mn>1</mn></msub><mo>-</mo><msub><mi>M</mi><mi>T</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>w</mi><mn>1</mn></msub><msub><mi>w</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><msub><mi>M</mi><mn>1</mn></msub><mo>-</mo><msub><mi>M</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mo>(</mo><msub><mi>w</mi><mn>1</mn></msub><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><mo>)</mo></mrow></mfrac></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>M</mi><mi>T</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>M</mi><mn>1</mn></msub><msub><mi>w</mi><mn>1</mn></msub><mo>+</mo><msub><mi>M</mi><mn>2</mn></msub><msub><mi>w</mi><mn>2</mn></msub></mrow><mrow><mo>(</mo><msub><mi>w</mi><mn>1</mn></msub><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><mo>)</mo></mrow></mfrac></mrow></math>]]></maths>当<img file="FDA0000416693880000022.GIF" wi="136" he="70" />取得最大值,即<img file="FDA0000416693880000023.GIF" wi="54" he="70" />为最大值,<img file="FDA0000416693880000024.GIF" wi="65" he="66" />为最小值时,阈值k为最佳阈值,通过最佳阈值k对图像进行阈值分割,将图像转化为目标和背景分离的二值图像;(2)熔池图像腐蚀处理:通过以下数学表达式获得熔池图像腐蚀E(X)的处理结果:<maths num="0006"><![CDATA[<math><mrow><mi>E</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mo>{</mo><mi>a</mi><mo>/</mo><msub><mi>B</mi><mi>a</mi></msub><mo>&Subset;</mo><mi>X</mi><mo>}</mo><mo>=</mo><mi>X&Theta;B</mi></mrow></math>]]></maths>其中,X为被处理图像;B为结构元素;a为像素点;(3)熔池图像膨胀处理:通过以下数学表达式获得熔池图像膨胀D(X)的处理结果:<maths num="0007"><![CDATA[<math><mrow><mi>D</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mo>{</mo><mi>a</mi><mo>/</mo><msub><mi>B</mi><mi>a</mi></msub><mo>&cap;</mo><mi>X</mi><mo>}</mo><mo>=</mo><mi>X</mi><mo>&CirclePlus;</mo><mi>B</mi></mrow></math>]]></maths>其中,X为被处理图像;B为结构元素;a为像素点;(4)熔池图像边缘检测,具体过程如下:熔池图像边缘F(X)是通过熔池图像膨胀D(X)与熔池图像腐蚀E(X)相减获得的,数学表达式如下式所示:F(X)=D(X)-E(X)步骤D、对熔池形状参数进行提取,具体过程如下:(1)熔宽的提取:在图像坐标系OXY中以熔池中心点o(X<sub>0</sub>,Y<sub>0</sub>)为起点沿oy方向搜索边缘点,当遇到边缘点A时停止搜索,记下A点的坐标值A(X<sub>A</sub>,Y<sub>A</sub>),则熔池宽度w<sub>a</sub>的计算公式为:<maths num="0008"><![CDATA[<math><mrow><msub><mi>w</mi><mi>a</mi></msub><mo>=</mo><mn>2</mn><msqrt><msup><mrow><mo>(</mo><msub><mi>X</mi><mn>0</mn></msub><mo>-</mo><msub><mi>Y</mi><mi>A</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Y</mi><mn>0</mn></msub><mo>-</mo><msub><mi>Y</mi><mi>A</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths>通过系统标定后得到的图像像素与实际尺寸的对应关系即可得到熔宽的实际尺寸w<sub>b</sub>,具体计算公式如下:<maths num="0009"><![CDATA[<math><mrow><msub><mi>w</mi><mi>b</mi></msub><mo>=</mo><mn>2</mn><msub><mi>w</mi><mi>y</mi></msub><msqrt><msup><mrow><mo>(</mo><msub><mi>X</mi><mn>0</mn></msub><mo>-</mo><msub><mi>X</mi><mi>A</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Y</mi><mn>0</mn></msub><mo>-</mo><msub><mi>Y</mi><mi>A</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths>上式所述的w<sub>y</sub>为熔池长度;(2)熔池长度的提取:在图像坐标系OXY中以熔池中心点o(X<sub>0</sub>,Y<sub>0</sub>)为起点沿ox方向搜索边缘点,当遇到边缘点B时停止搜索,记下B点的坐标值B(X<sub>B</sub>,Y<sub>B</sub>),则熔池长度l的计算公式为:<maths num="0010"><![CDATA[<math><mrow><mi>l</mi><mo>=</mo><mn>2</mn><msqrt><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>o</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>B</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Y</mi><mi>o</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>B</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths>通过系统标定后得到的图像像素与实际尺寸的对应关系即可得到熔池长度的实际尺寸,即:<maths num="0011"><![CDATA[<math><mrow><mi>l</mi><mo>=</mo><mn>2</mn><msub><mi>w</mi><mi>x</mi></msub><msqrt><msup><mrow><mo>(</mo><msub><mi>X</mi><mn>0</mn></msub><mo>-</mo><msub><mi>X</mi><mi>B</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Y</mi><mn>0</mn></msub><mo>-</mo><msub><mi>Y</mi><mi>B</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths>其中,在焊接过程中采用导电布缠绕在摄像机上。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号