Recognition algorithm of small-diameter tube X-ray welding defect image

XIAO Yang1 ，GAO Weixin1，2，DENG Guohao2 （1.School of Electronic Engineering，Xi'an Shiyou University，Xi'an 71OO65，China; 2.Key Laboratory of Gas-Oil Logging Technology，Xi'an Shiyou University，Xi’an 710065，China）

Abstract： To address the current situation of low accuracy rate of small-diameter tube welding image defect detection，by combining image feature analysis and sparse dictionary learning，a small-diameter tube welding defect detection algorithm based on image segmentation is proposed.Firstly，using two-step image segmentation way acquires the region of interest which is in small-diameter tube welding image.Secondly，the suspected defectregion is obtained by extracting welding defect.Finally，we propose a mathematical model of the dictionary matrix of small-diameter tube welding defects with the objective of minimizingcorelations between diferent types of atoms and solve it by using K-SVD algorithm. After that，the dictionary matrix is used to clasify circular defects，strip defects and noise. To improve the real-time performance of the system，we use parallel programming to accelerate the image segmentation algorithm.The results show that the recognition rate of the proposed method is O.974 for circular defects and O.967 for strip defects，and the recognition speed is fast，which enables the effctive recognition of defects in small-diameter tube welding image.

Key words： smal-diameter tube;welding defect; image segmentation; sparse dictionary learning

0 前言

小徑管結(jié)構在油氣地面裝置中的應用十分廣泛。在焊接過程中，受外部環(huán)境和焊接條件的影響，其焊縫區(qū)域不可避免的產(chǎn)生諸如：氣孔，裂紋，夾渣和焊瘤等缺陷[1]。這些缺陷不僅影響工件的外觀，還會影響工件的結(jié)構強度和性能。更嚴重的是，缺陷會對工件使用者造成一些潛在的安全隱患[2-4]。因此，焊縫質(zhì)量檢測是焊接過程結(jié)束后的一道關鍵環(huán)節(jié)。無損檢測技術可以在不影響被檢測對象性能的基礎上探測物體在制造與使用過程中的缺陷問題[5-6]。在施工現(xiàn)場的檢測流程中，主要采用無損檢測方式作為焊縫缺陷的檢測方法。

小徑管焊縫缺陷的無損檢測方法主要包括超聲檢測和射線檢測。其中，超聲相控陣檢測技術可以檢測缺陷的高度和埋深，且具有高精度和環(huán)保的特點，成為當前小徑管缺陷檢測領域的研究熱點?；诔暀z測的方法，江野等人將超聲相控陣檢測技術應用于奧氏體異種鋼薄壁小口徑管道對接焊縫的裂紋檢測。陳小明等人[8探究了將超聲相控陣檢測技術應用于蛇形管高壓加熱器超小徑管焊縫檢測的可行性。

近年來，伴隨射線檢測技術逐漸成熟，部分學者也試想將射線檢測技術應用于小徑管焊縫缺陷檢測中。Xie等人將超聲相控陣檢測和射線檢測對小徑管焊縫缺陷的檢測效果進行對比，結(jié)果表明射線檢測對體積型缺陷的檢測更為精準，檢測結(jié)果更直觀且可永久保存。賈瀟等人[1°的研究也表明在對缺陷類型的判別上，射線檢測更具有優(yōu)勢。然而，目前未見基于射線檢測技術實現(xiàn)小徑管焊縫缺陷檢測的相關算法和文獻。

綜上所述，該文從射線檢測的角度出發(fā)，針對小徑管結(jié)構，通過兩步分割獲取感興趣區(qū)域（theregionofinterest，ROI），提取缺陷疑似局部圖像（suspecteddefectregion，SDR），構造字典矩陣數(shù)學模型并分析了字典矩陣原子數(shù)對識別精度的影響，最終確定了具體的原子數(shù)，使用字典矩陣對數(shù)據(jù)庫中的圖像進行檢測，驗證了所提方法的可行性，在驗證過程結(jié)束后，對部分缺陷識別錯誤的原因進行可能性分析。

1 試驗方法

1.1 圖像分析

實際工業(yè)生產(chǎn)中，小徑管的X射線檢測圖像如圖1所示，圖1（b）為小徑管在雙壁雙影橢圓透照下的圖像，圖1（c）為小徑管在雙臂雙影垂直透照下的圖像。

垂直透照下，小徑管焊縫圖像對于面積型缺陷的檢出率較高，且面積型缺陷如裂紋等對焊管的危害性相比于氣孔等缺陷更大，故選取垂直透照下的小徑管圖像作為研究對象。

由圖1可知，射線圖像包含背景區(qū)域和小徑管區(qū)域，缺陷存在于小徑管中的焊縫區(qū)域。由于焊縫圖像在采集時會受到光照和表面反射的影響，所以背景存在一些較暗的陰影區(qū)域[\"]，陰影區(qū)域像素值范圍與小徑管區(qū)域相差較大，可以通過圖像分割的方式提取小徑管區(qū)域。部分小徑管射線圖像表現(xiàn)為低對比度，分割時要求算法具有較強的魯棒性，且實際的小徑管射線圖像都具有高噪聲的特點，分割算法應具有良好的降噪效果。

1.2 ROI提取

X射線焊縫圖像背景起伏，背景區(qū)域的存在既降低算法的運行速度，又對缺陷識別產(chǎn)生干擾。因此，提取ROI區(qū)域是進一步實現(xiàn)缺陷識別的重要中間環(huán)節(jié)。從原始圖像可知，焊縫區(qū)域處于小徑管區(qū)域內(nèi)，難以通過分割算法直接提取焊縫區(qū)域。故采用兩步分割的方法，即先提取小徑管區(qū)域，再從小徑管區(qū)域圖像中提取焊縫區(qū)域。具體的圖像分割算法流程如下。

1.2.1 小徑管區(qū)域提取

（1）圖像增強與降噪：利用Gamma變換增強小徑管X射線圖像S得到圖像 S₁ ，使用中值濾波和均值濾波算法對圖像 S₁ 進行降噪處理得到圖像 S₂ 。

（2）二值化：使用Ostu分割 S₂ 得到二值圖像 B₁ 。

（3）獲取小徑管區(qū)域連通域：使用種子填充法提取B₁ 中所有連通域區(qū)域，并用模板匹配算法對每個連通域進行處理，含有模板T的連通域即為小徑管區(qū)域L。

（4）獲取小徑管區(qū)域圖像：利用小徑管區(qū)域連通域

L的外接矩形獲得對應原圖 S₂ 的小徑管區(qū)域，記為 M₁ 。

1.2.2 焊縫區(qū)域提取

（1）二值化：使用 Ostu 分割圖像 M₁ 得到二值圖像 B₂ 。（2）獲取最大連通域：使用種子填充法提取圖像B₂ 連通域，對比連通域面積，將小于最大面積的連通域設定為背景區(qū)域，處理后的圖像記為 B₃ 。（3）統(tǒng)計圖像 B₃ 每列像素的平均值，若平均值超過一定閾值，記此列為焊縫子列，所有焊縫子列的組合區(qū)域即為感興趣焊縫區(qū)域 ROI 。

為驗證小徑管焊縫分割算法的合理性，使用某石油公司提供的焊接數(shù)據(jù)庫對小徑管X射線圖像進行ROI提取試驗。數(shù)據(jù)庫中共有3535張現(xiàn)場實拍圖像，通過上述方法對3535張圖像進行分割試驗，結(jié)果表明有3487張圖像可以準確提取ROI，有效分割率為98.6% ，剩余48張圖像識別錯誤經(jīng)檢查原因為質(zhì)量很差不具有分割意義或中心標記與管道邊緣嚴重重合。分割試驗過程中，所提方法表現(xiàn)出良好的魯棒性，對對接接頭焊口圖像和馬鞍形焊口圖像的ROI分割都表現(xiàn)出良好的效果分別如圖2和圖3所示。

分割過程中，即使部分圖像質(zhì)量較差或者圖像中管道區(qū)域有彎曲如圖4和圖5所示，分割算法也可以準確分割。

1.3字典矩陣模型及求解

小徑管焊縫區(qū)域內(nèi)缺陷類型較多，不同類型的缺陷具有不同的特征，對焊管的危害性也不同。因此，在缺陷的識別過程中，不僅要確定是否為缺陷，還要精準判定缺陷的類型。圓形缺陷和線形缺陷是小徑管焊縫區(qū)域內(nèi)最為常見和危害性較大的兩類缺陷，噪聲是缺陷識別過程中的干擾因素，將三者作為識別對象可有效覆蓋大部分缺陷類型。利用文獻[12]中的灰度密度聚類算法可獲得焊縫區(qū)域內(nèi)圓形缺陷，線形缺陷和噪聲所對應的SDR圖像。圖6為SDR示意圖，圖中的外切矩形區(qū)域為疑似缺陷的外切矩形，擴展5個像素點后的區(qū)域為SDR。按其方法提取出的部分SDR如圖7所示。SDR面積較小，傳統(tǒng)基于特征提取的方法難以獲得圖像的顯著特征，識別過程中容易誤判一些細小缺陷。稀疏字典學習可通過稀疏系數(shù)直接識別SDR類型，避免了特征提取過程，在小徑管焊縫缺陷識別過程中具有獨特優(yōu)勢。

稀疏字典學習是從龐大數(shù)量的樣本中找出最具代表性的原子組成字典，并用字典中的原子對樣本進行稀疏描述[13]。設 為待檢測樣本，稀疏字典學習的表達式為

y=d₁x₁+d₂x₂+…+d_nx_n=DX

式中 D=[d₁，d₂，…，d_k]∈R^m×k ，表示字典矩陣中的原子數(shù)為 k ，表示待檢測樣本 y 可由 n 個 k 維列向量線性表示。

小徑管焊縫缺陷字典矩陣中，原子是由SDR的像素點按照列優(yōu)先的原則逐列排列而成，表示為

D={d₁，d₂，…，d_cd，d_cd+1，d_cd+2，…，d_cd+Ld，d_cd+Ld+1，

d_Cd+Ld+2，…，d_Cd+Ld+Nd}

式中： C_d 表示圓形缺陷原子數(shù)； L_d 表示線形缺陷原子數(shù)； N_d 表示噪聲原子數(shù)。稀疏系數(shù) X 中的非0系數(shù)與待檢測圖像類別有關，與類別無關的系數(shù)表示為0。對X 的求解可等效為求2范數(shù)最小化問題，即

式（3）可等效為

可使用正交匹配追蹤（orthogonalmatchingpursuit，OMP）算法求解式（4），得到矩陣 X 的1個廣義解。

1.3.1 構建字典矩陣模型

在小徑管焊縫缺陷識別過程中，可通過學習的方式獲得字典矩陣 ，學習過程等效為求解下式

式中： Y 為訓練樣本矩陣； ε 為矩陣 X 的稀疏性。

小徑管 ΔX 射線焊縫圖像SDR可分為圓形SDR，線形SDR和噪聲SDR共3類。相同類型SDR具有一定的相似性，不同類型SDR表現(xiàn)出不同的特征，可以認為特定種類的SDR可由特定字典進行稀疏描述。因此，只需分別構建3種SDR的特定字典并進行組合，即可實現(xiàn)字典矩陣 的構建， D 的形式如式（2）所示。若類別數(shù)為 C ，可通過式（6）求解特定種類SDR對應的字典矩陣 D_κ 。

式中： Y_K 為第 K 類訓練樣本； X_K 為第 K 類的稀疏系數(shù)。

字典矩陣 可表示為

D=[D₁，D₂，…，D_K]

式（6）雖可求解矩陣 D_K ，但并未對 D_κ 中的原子數(shù)進行限制。實際構成 時，訓練樣本可能較多，此時k 值較大，較大的 k 值雖有利于對更加多元化的數(shù)據(jù)進行表示，但也會增加不同缺陷類型原子間的相關性。若 D 中的原子完全正交，則僅需少量原子就可覆蓋缺陷的全部特征，極大地減少可行解的數(shù)量。鑒于此，該文從小徑管X射線焊縫圖像缺陷識別的唯一性和實時性角度出發(fā)，提出一種以不同類型缺陷原子相關性最小為目標的小徑管焊縫缺陷字典矩陣學習模型，即

s.t.?q，∣∣x_q∣∣?ε

式中： Θ 為所有樣本的集合， P∈{ 圓形SDR，線形SDR，噪聲SDR ，其中 P_i 表示對應類型原子區(qū)域的第i 列元素，如當 P 為圓形SDR原子時， P_i 表示圓形SDR原子區(qū)域的第 i 列元素，其余變量與式（6）相同。

式（8）中，字典矩陣 和稀疏編碼矩陣 X 均為未知量，可以使用K均值奇異值分解（K-SVD）算法，先確定 后計算 X ，再固定 X 更新 的方式交替更新求解。

1.3.2 字典矩陣求解

利用K-SVD算法求解字典矩陣的算法流程如圖8所示，算法運行前需預先設定字典矩陣的原子個數(shù)k 。圖8中，mink取值為20，maxk取值為80。若 k 值過小則無法完整描述SDR的特征， k 值過大則使類間相似性原子增多，字典冗余性增大。再考慮到要將3種缺陷類型的字典矩陣通過式（7）進行組合，故將各缺陷類型字典矩陣原子數(shù)的最大值設為80。

確定3種缺陷的字典后，按照圓形缺陷、線形缺陷、噪聲的順序?qū)⒃舆M行排序，記錄此時字典矩陣中圓形缺陷原子數(shù) C_d ，線形缺陷原子數(shù) L_d 和噪聲原子數(shù)N_d 。對任意測試樣本 y ，在 D 確定后可根據(jù)式（1）求解稀疏系數(shù) X ，根據(jù) X 中非0系數(shù)的分布可以判定測試樣本的缺陷類型，判別方法如式（9）式（10）和式（11）所示。若滿足式（9）判定為圓形缺陷，滿足式（10）判定為線形缺陷，滿足式（11）判定為噪聲。

學習過程中，可能會出現(xiàn)因訓練樣本不同導致3種缺陷字典中的原子數(shù)互不相同的情況。為探究這種因素的影響，設計了1個 k 值為80的字典，字典中圓形缺陷SDR設定為50，線形缺陷SDR設定為25，噪聲SDR設定為5，利用此字典進行識別試驗。結(jié)果表明， X 中的圓形缺陷段系數(shù)具有明顯優(yōu)勢，識別結(jié)果傾向于判定為圓形缺陷。將3種缺陷SDR數(shù)量對調(diào)，識別效果沒有明顯變化。當3種缺陷SDR分別設定為27，27和26時，從識別結(jié)果看3種缺陷的識別率沒有明顯差別。因此，為獲得更好的識別效果，應保證組合后字典矩陣中各缺陷SDR大致持平，故令 C_d=L_d =N_d=k 。

1.4數(shù)據(jù)準備及評價指標

該文試驗條件為InterCorei7-9700、NIDIAGe-ForceGTX1660、16GB內(nèi)存、64位Windows操作系統(tǒng)、VisualStudio2019，PyCharm2022.1.2。為驗證所提算法的有效性，對某石油鉆采公司提供的小徑管X射線圖像進行識別試驗，數(shù)據(jù)集為3535幅，大小為900像素 ×1 200 像素的焊縫圖像。在獲取數(shù)據(jù)集后，利用該文提出的圖像分割算法分割數(shù)據(jù)集圖像，提取缺陷SDR并將SDR歸一化為 15×15 的模板，保存在SDR數(shù)據(jù)庫中。隨機提取SDR數(shù)據(jù)庫中的 40% 構建字典矩陣，其余 60% 用以驗證字典矩陣分類的有效性。在對分類模型的性能進行描述時，選取混淆矩陣中的準確率 A ，召回率 R ，精準率 P 和F1-score共4項指標，計算方式依照文獻12」，選取準確率 A 作為衡量分類模性能的主要指標。

2 試驗結(jié)果與分析

根據(jù)圖8的K-SVD算法求解字典矩陣的算法分別求解圓形缺陷字典矩陣 D_c 、線形缺陷字典矩陣 D_L 和噪聲字典矩陣 δ?D_N ，并依據(jù)式（7）構成小徑管X射線焊縫缺陷字典矩陣 D 。不同原子數(shù)下，字典矩陣 的識別率和3種缺陷原子間的互相關性如圖9所示。

從圖9可知，不同類型原子間的互相關性 μ 隨 k 值增大而增大。當 k=57 時，字典矩陣 的缺陷識別率達到最大值為 97.2% ，當 k 在一定范圍內(nèi)繼續(xù)增大時，缺陷識別率在 97% 上下小幅震蕩，在 k 值為58，59，62和64時識別率也可達到最大值，若 k 值增大到80附近，此時識別率因字典冗余性過大有所降低。綜上所述，從識別實時性以及不同類原子間互相關性最小兩個角度考慮，應在保證識別率的前提下使 k 值盡可能小。

選取 k=57 作為最佳原子數(shù)，此時字典矩陣 D 的原子數(shù)為171。使用字典矩陣 D 分別對圓形缺陷和線形缺陷的測試樣本進行識別，圓形缺陷識別率為97.4% ，線形缺陷識別率為 96.7% ，分類結(jié)果如表1所示，表1中混淆矩陣的4項指標為所有缺陷相對噪聲。實際試驗中還發(fā)現(xiàn)一些未能正確識別的缺陷SDR，這些SDR或缺陷區(qū)域面積極小，或缺陷顏色較淺。該文統(tǒng)計了未正確識別SDR的缺陷區(qū)域面積 d_s 歸一化后和灰度跨度 P_r ，正確識別SDR的缺陷區(qū)域面積 d_s 的統(tǒng)計結(jié)果與未正確識別SDR的缺陷區(qū)域面積 d_s 的統(tǒng)計結(jié)果對比如圖10所示.正確識別的SDR的灰度跨度 P_r 的統(tǒng)計結(jié)果與未正確識別SDR的灰度跨度 P_r 的統(tǒng)計結(jié)果對比如圖11所示。由于數(shù)據(jù)量過小的原故，目前試驗結(jié)果無法證明缺陷面積對識別準確率的影響，但從圖11可知，正確識別的SDR與未正確識別的SDR在灰度跨度上具有明顯差異，未正確識別的SDR灰度跨度普遍低于43。究其原因，該文認為當圖像灰度跨度低于43時，缺陷與背景間對比度較低，稀疏系數(shù) X 中的非0值分布分散，因此識別時易產(chǎn)生誤判。數(shù)據(jù)庫中也存在部分SDR表現(xiàn)出具有其它缺陷的圖像特征，這種類型的缺陷在分類時更容易出錯，對這種難以分類的SDR還需進行更加深入的研究。所提算法包括ROI區(qū)域分割，SDR提取和缺陷識別共3個步驟，在算法運行過程中，步驟一耗時遠大于步驟2和步驟3。從系統(tǒng)實時性考慮，應使用GPU對步驟1進行加速，以處理1幅圖像為例，平均加速效果如表2所示。

隨著計算機和人工智能的發(fā)展，許多學者從傳統(tǒng)的主觀提取X射線焊縫圖像特征，轉(zhuǎn)向用自動學習的方式提取X射線焊縫圖像特征[14-15]。近年來，深度學習的發(fā)展也使越來越多的學者將各種不同的深度學習模型應用到自己的研究當中并取得優(yōu)異的成果。為了證明所提算法的優(yōu)越性，將相同數(shù)據(jù)集在VGG16，GoogleNet，ResNet18和YoLoV3共4個網(wǎng)絡模型上進行測試對比，結(jié)果如表3所示。

從表3可以看出所提算法相比于一些深度學習模型對焊縫缺陷的檢測精度更高，這可能是由于缺陷面積相對于圖像過小導致深度學習模型無法提取有效特征的緣故，而稀疏字典學習由于不需要提取特征，因此相比于深度學習模型而言在解決小面積對象分類問題上更有優(yōu)勢。從檢測速度的角度出發(fā)，所提算法檢測速度較慢，在之后的研究中可嘗試通過深化程序的并行程度或者結(jié)合基于深度模型的圖像分割方法加快檢測速度。

3結(jié)論

（1）采用兩步圖像分割法可準確提取出小徑管焊縫圖像感興趣ROI，極大降低了提取SDR和缺陷分類的成本。經(jīng)過并行加速后，該分割方法在實際應用中是快速有效的，ROI分割時間降低為1.12s，有效分割率達到 98.6% 。（2）與深度學習方法相比，通過稀疏字典學習的方式對焊縫缺陷進行識別避免了特征提取的難題，學習成本更低，計算速度更快，準確率也更高，在實踐中展現(xiàn)出了巨大的優(yōu)越性。（3）字典矩陣僅需較小的原子數(shù)就可以達到很高的缺陷識別準確率，即使學習過程完成，也可隨時將典型缺陷原子加入字典，最終形成過完備字典。

參考文獻

[1] 周兆逸，張亞南，王肖鋒，等．基于改正二維主成分分析的焊縫表面缺陷檢測[J]．焊接學報，2021，42（11）：70 -76.

[2] YangL，WangHX，HuoBY，etal.Anautomaticweldingdefectlocation algorithmbased on deep learning[J].NDTamp;EInternational：Independent Nondestructive Testing andEvaluation，2021，120：102435.

[3] SunJ，LiC，WuXJ，etal.AnEffective methodofwelddefect detection and classification based on machine vision[J].Ieee Transactions onIndustrial Informatics，2O19，15（12）：6322 -6333.

[4]ZhangY，You D，Gao X，etal.Welding defects detection-based on deep learning with multiple optical sensors duringdisklaser welding of thick plates[J]. Journal ofManufactur-ingSystems，2019，51：87-94.

［5］魏博文，高煒欣，劉夢溪．激光超聲結(jié)合CNN的鑄件缺陷檢測方法[J].激光與紅外，2022，52（9）：1327-1334.

[6]ZhangL，ZhangYJ，DaiBC，et al.Welding defect detec-tionbased on local image enhancement[J]. IET Image Pro-cessing，2019，13（13）：2647-2658.

[7] 江野，王東，秦小陽，等．奧氏體異種鋼薄壁小徑管焊縫裂紋超聲相控陣檢測［J]．焊接學報，2018，39（8）：119 -122.

[8] 陳小明，譚云華，賴傳理，等．相控陣技術在超小徑管焊縫檢測應用研究[J].鍋爐技術，2021，52（S1）：55-58

[9] XieL，Wang T，ZhangY.Comparative experimental studyon phased array and X-ray detection of small diameterpipeweld[J]. Journal ofPhysics，2021，1885（3）：032023.

［10］賈瀟，趙?。嗫仃嚦暸c射線檢測對石油化工管道缺陷檢測的對比分析[J]．中國石油和化工標準與質(zhì)量，2022，42（15）：174-175.

［11］張小琳，劉祎，白賚諷，等．基于背景估計的焊縫缺陷檢測[J]．電子測量技術，2022，45（14）：116-122.

[12] 高煒欣，胡玉衡，穆向陽，等．基于聚類的埋弧焊X射線焊縫圖像缺陷分割算法及缺陷模型［J]．焊接學報，2012，33（4）：37-41.

[13］孫春紅，丁廣太，方坤．基于稀疏字典學習的羊絨與羊毛分類[J]．紡織學報，2022，43（4）：28-32.

[14] Hou W H， Zhang D S，Wei Y，et al.Review on computeraided weld defect detection from radiography images[J].AppliedSciences，2020，10（5）：1878-1893.

[15］修延飛，李海超，胡廣澤，等．一種用于穿孔塞焊焊縫特征提取的視覺識別算法［J]．焊接學報，2020，41（2）： 75-79.肖揚簡介：碩士；主要研究方向為焊縫缺陷無損檢測；Email：xymoh981017@163.com.

通信

作者簡介：高煒欣，博士，教授，碩士研究生導師；Email：wxgao@ xsyu. edu. cn.