发明名称 一种焊接缺陷超声相控阵扇形扫描图像特征提取方法
摘要 一种焊接缺陷超声相控阵扇形扫描图像特征提取方法,属于焊接缺陷研究技术领域。本发明解决了焊接缺陷分类识别时直接采用缺陷超声相控阵扇形扫描图像进行缺陷识别,图像数据维数大,分类模型复杂,学时间长的难题。主要步骤:取焊接缺陷超声相控阵扇形扫描图像,构建缺陷图像数据矩阵;计算缺陷图像数据矩阵的协方差矩阵;求协方差矩阵的特征值与特征向量;确定缺陷图像主成分个数,建立焊接缺陷特征表达函数。本发明的焊接缺陷超声相控阵扇形扫描图像特征提取方法既保留了缺陷图像决大部分信息,表征了缺陷的类型,保证了分类模型的正确识别率,又降低了缺陷图像数据的维数,极大地提高了分类模型的学速度。
申请公布号 CN103268491A 申请公布日期 2013.08.28
申请号 CN201310209299.7 申请日期 2013.05.30
申请人 哈尔滨工业大学 发明人 栾亦琳;刚铁;冯吉才;张秉刚;孙涛
分类号 G06K9/46(2006.01)I 主分类号 G06K9/46(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 1.一种焊接缺陷超声相控阵扇形扫描图像特征提取方法,其特征在于:所述方法的具体过程为:步骤一、取n幅缺陷图像,每幅缺陷图像含有p数据,其中p=k×l×t,k为每幅缺陷图像中超声波A信号数目,l为每个超声波A信号中含有的数据,t为数据中包含的回波幅度和缺陷三维坐标信息,这样共得到n×p个数据,缺陷图像数据矩阵X如下:<maths num="0001"><![CDATA[<math><mrow><mi>X</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mn>11</mn></msub></mtd><mtd><msub><mi>x</mi><mn>12</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>x</mi><mrow><mn>1</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>21</mn></msub></mtd><mtd><msub><mi>x</mi><mn>22</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>x</mi><mrow><mn>2</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>x</mi><mi>np</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤二、计算缺陷图像数据矩阵X的协方差矩阵∑:∑=(σ<sub>ij</sub>)<sub>p×p</sub>    (1-9)其中σ<sub>ij</sub>=Cov(X<sub>i</sub>,X<sub>j</sub>)(i,j=1,2,…,p),Cov(X<sub>i</sub>,X<sub>j</sub>)是求X<sub>i</sub>和X<sub>j</sub>的协方差;步骤三、求协方差矩阵∑的特征值与特征向量;由公式(1-10)求得∑的p个的特征值λ<sub>i</sub>(i=,2,…,p),按由大到小的顺序排列为λ<sub>1</sub>≥λ<sub>2</sub>…≥λ<sub>p</sub>;根据公式(1-11)计算λ<sub>i</sub>对应的标准正交特征向量α<sub>i</sub>(i=1,2,…,p):|λE-∑|=0    (1-10)(λE-∑)α=0    (1-11)步骤四、确定缺陷图像主成分个数m(m<p),m确定过程为:利用公式(1-6)求得第m个缺陷图像主成分的方差贡献率,利用公式(1-7)求得前m个缺陷图像主成分累积贡献率:<maths num="0002"><![CDATA[<math><mrow><msub><mi>&lambda;</mi><mi>m</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>&lambda;</mi><mi>j</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>&lambda;</mi><mi>j</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>当前m个缺陷图像主成分累积贡献率满足大于等于85%的条件:<img file="FDA00003276676600014.GIF" wi="425" he="142" />即可确定m的值;步骤五、在得到缺陷图像主成分个数m后,构建焊接缺陷特征表达函数:F<sub>i</sub>=α<sub>i1</sub>X<sub>1</sub>+α<sub>i2</sub>X<sub>2</sub>+…+α<sub>ip</sub>X<sub>p</sub> i=1,2,…,m  (1-12)其中F<sub>1</sub>称为特征1,F<sub>2</sub>称为特征2,以此类推,F<sub>m</sub>称为特征m。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号