薛倩 馬敏 王化祥 詹湘琳

(中國民航大學(xué)航空自動化學(xué)院1,天津 300300;天津大學(xué)電氣工程與自動化學(xué)院2,天津 300072)

用于復(fù)合材料缺陷可視化檢測的l1正則化算法

薛倩1馬敏1王化祥2詹湘琳1

(中國民航大學(xué)航空自動化學(xué)院1,天津 300300;天津大學(xué)電氣工程與自動化學(xué)院2,天津 300072)

隨著復(fù)合材料在民用飛機(jī)中越來越廣泛的應(yīng)用,復(fù)合材料缺陷的在線檢測對無損檢測技術(shù)提出了更高要求。為提高平面電容層析成像(ECT)系統(tǒng)的成像質(zhì)量,從而提高其對復(fù)合材料缺陷的檢測精度,對基于l1正則化的成像算法進(jìn)行了研究。采用交替最小化方法進(jìn)行數(shù)值計(jì)算,并結(jié)合迭代閾值法,對重建圖像進(jìn)行降噪。試驗(yàn)結(jié)果表明,與ECT中常用的Landweber、Tikhonov正則化等算法相比,該方法能夠有效分辨物場中的不同介質(zhì),提高了成像質(zhì)量和被測物體邊緣分辨率。

層析成像 無損檢測 圖像重建 復(fù)合材料 正則化

0 引言

先進(jìn)復(fù)合材料(acomposite materials,ACM)因其質(zhì)量輕,隔熱、隔音、減震效果好,比強(qiáng)度高、抗腐蝕、耐高(低)溫等特點(diǎn),在現(xiàn)代民用飛機(jī)中獲得了廣泛應(yīng)用[1]。新研制的波音系列飛機(jī)的機(jī)翼、機(jī)身等主承力機(jī)構(gòu)均由ACM制成,空客A380也大量采用了ACM[2]。由于ACM具有各向異性等特點(diǎn),其損傷機(jī)理、失效形式都與金屬材料顯著不同,這給無損檢測技術(shù)帶來了新的挑戰(zhàn)。當(dāng)非金屬復(fù)合材料發(fā)生疲勞或損傷時(shí),其介電常數(shù)也發(fā)生相應(yīng)變化。因此,理論上可利用基于電容敏感機(jī)理的電容層析成像(electrical capacitance tomography,ECT)技術(shù)實(shí)現(xiàn)在線監(jiān)測[3-4]。平面ECT技術(shù)用于材料無損檢測的研究起步較晚,基本還處于仿真試驗(yàn)階段[5-7],主要的困難在于該技術(shù)空間分辨率較低。因此,本文通過圖像重建算法的研究,提高系統(tǒng)的空間分辨率。

1 系統(tǒng)結(jié)構(gòu)與數(shù)學(xué)模型

平面ECT系統(tǒng)主要由三部分組成:平面陣列電容傳感器、電容測量與數(shù)據(jù)采集電路和成像計(jì)算機(jī)。平面ECT系統(tǒng)結(jié)構(gòu)框圖如圖1所示。

圖1 平面ECT系統(tǒng)結(jié)構(gòu)框圖Fig.1 Structure of the planar ECT system

ECT系統(tǒng)通過傳感器的陣列電極獲得敏感場內(nèi)介質(zhì)分布狀況的邊界信息。電容測量及數(shù)據(jù)采集電路對傳感器的輸出電容值進(jìn)行濾波、變換、放大,并將數(shù)據(jù)傳給成像計(jì)算機(jī)。成像計(jì)算機(jī)利用圖像重建算法重構(gòu)被測介質(zhì)的分布情況,同時(shí)向數(shù)據(jù)采集系統(tǒng)發(fā)送控制信號。

根據(jù)似穩(wěn)場理論的Maxwell方程,對于均勻、線性、各向同性的介質(zhì),ECT敏感場內(nèi)任意一點(diǎn)滿足拉普拉斯方程,即:

由電場的唯一性定理可知,滿足拉普拉斯方程和給定邊界條件的φ(x,y)是給定靜電場的唯一解。本文中ECT系統(tǒng)采用電壓激勵(lì),激勵(lì)電極的電勢為U,其余電極電勢為0,符合Direchlet邊界條件,即:

式中:?Ω為場域邊界。

電勢分布φ(x,y)與介電常數(shù)分布e(x,y)通過積分方程相聯(lián)系。基于擾動理論和有限元法,得到線性化和離散化的圖像重建模型為:

將式(3)歸一化處理后,表示為:

式中:C為M×1維的測量數(shù)據(jù)向量,M代表獨(dú)立測量數(shù);g為N×1維的介電常數(shù)分布圖像向量,N代表被測場域的剖分單元數(shù);S為M×N維的靈敏度矩陣。

2 圖像重建算法研究

2.1 數(shù)學(xué)模型變換

在實(shí)際ECT系統(tǒng)中,可測量的邊界值數(shù)目往往遠(yuǎn)小于剖分單元數(shù),即線性方程組(4)是欠定的,因此解不存在或不唯一。通常將式(4)的求解轉(zhuǎn)化為求最小二范數(shù)解的凸優(yōu)化問題[8-9],即:

式中:A為式(4)中的靈敏度矩陣S;b為測量數(shù)據(jù)向量C;x為要求解的圖像向量g。

求解此問題常用的算法有共軛梯度法和Landweber迭代法等。由于測量值C中不可避免地包含噪聲,而靈敏度矩陣S的條件數(shù)通常很大,為了較準(zhǔn)確和穩(wěn)定地重建g,一般采用Tikhonov正則化方法將圖像重建問題轉(zhuǎn)化為:

由于泛函數(shù)Ja(x)是光滑的,可直接令其導(dǎo)數(shù)為零,求得Ja(x)的極小點(diǎn)x=(ATA+αI)-1ATb作為Ax= b的正則化近似解。由此可證明ATA+αI的逆總是存在的,故解存在。當(dāng)A是線性算子時(shí),Ja(x)是嚴(yán)格凸的,故存在唯一的極小點(diǎn)。

基于L2范數(shù)的Tikhonov正則化方法,容易使解過度光滑化。復(fù)合材料缺陷檢測中被測區(qū)域的介質(zhì)往往存在不連續(xù)性,采用該方法不可避免地會降低空間分辨率;且如果A是非線性的,Ja(x)不再是嚴(yán)格凸的,可能有多個(gè)局部極小點(diǎn)。不同的迭代初值可能導(dǎo)致不同的極值點(diǎn),可能最后得到不合理的解。為削弱圖像的平滑性,可考慮基于稀疏性的l1正則化。1992年, Rudin等人提出總變差(total variation,TV)正則化算法[10]并將其應(yīng)用于圖像去噪。后來,該算法分別由Fang(2004)、Soleimani(2005)和Wang(2007)等人引入ECT圖像重建[11-13]?；赥V正則化的圖像重建算法以成像區(qū)域內(nèi)像素值的總變差作為懲罰項(xiàng),由于同樣的TV值允許函數(shù)有不同的表達(dá)形式,包括不連續(xù)的函數(shù)形式,這就使得介質(zhì)分布不連續(xù)的區(qū)域得到較好的重建。該算法在提高逆問題求解穩(wěn)定性的同時(shí),對介質(zhì)分布不連續(xù)區(qū)域的分辨力也明顯提高,具有良好的保邊緣性。ECT中圖像總變差的離散化形式為:

可見,TV正則化相當(dāng)于利用梯度算子將圖像稀疏化,然后求最小l1范數(shù)解。實(shí)際上,TV正則化可通過引入分離變量轉(zhuǎn)化為l1正則化形式[14]。由于增加了l1范數(shù)的懲罰項(xiàng),優(yōu)化問題不再是光滑的,求解難度增加。一種魯棒的算法是原始-對偶內(nèi)點(diǎn)法。由于該算法效率很低,Wang等利用總變差表達(dá)式的一種二次可微的近似形式,采用牛頓-拉夫遜迭代算法求解[13]。

除TV正則化外,由于大部分圖像都可經(jīng)離散余弦變換(discrete cosine transform,DCT)、離散傅里葉變換(discrete Fourier transform,DFT)或小波變換等獲得稀疏表達(dá)形式,在圖像重建中也可將稀疏化的圖像向量的l1范數(shù)作為懲罰項(xiàng),以提高空間分辨率,因此,圖像重建問題轉(zhuǎn)化為基于l1范數(shù)的凸優(yōu)化問題,即:

式中:μ＞0為權(quán)重參數(shù)。

基于凸優(yōu)化理論,稀疏向量x∈iN可由少量線性測量值b=Ax∈iM重建。由此,壓縮傳感領(lǐng)域求解l1問題的各種算法都可用于測量數(shù)據(jù)有限的CT/ECT圖像重建,唯一的問題在于算法效率。2011年,Yang J F等將優(yōu)化領(lǐng)域的經(jīng)典算法——交替方向法(alternating direction method,ADM)擴(kuò)展到8種不同形式的l1范數(shù)最小化問題[15]?；谠鰪VLagrangian函數(shù)和交替最小化思想的交替方向算法,在測量數(shù)據(jù)受噪聲污染的情況下能獲得較準(zhǔn)確的解。與現(xiàn)存的許多一階算法相比,該算法效率較高且通用性強(qiáng)。

2.2 數(shù)值求解算法

對式(8)中的變量進(jìn)行下列代換:

通過Legendre變換求得f和g的共軛函數(shù):

則式(8)的對偶問題可表示為:

令λ=-p,則式(10)等價(jià)于:

引入向量z∈iN,消除δz(q)項(xiàng),將式(11)轉(zhuǎn)化為帶約束的最大化問題,即:

利用增廣Lagrangian函數(shù),將式(12)轉(zhuǎn)化為無約束最小化問題,即:

式中:x∈iN為Lagrangian乘子(即原問題中的圖像向量x);β＞0為懲罰系數(shù)。

相對普通的Lagrangian函數(shù),增廣Lagrangian函數(shù)多一個(gè)關(guān)于約束的懲罰項(xiàng)。采用交替方向(ADM)算法,在每次迭代中交替更新原始變量x和對偶變量λ、z,先給x和z(或λ)賦初值,通過最小化增廣拉格朗日函數(shù),得到新的λ(或z),并更新z(或l),由新的l和z更新乘子x,依此迭代,最終x將收斂到原問題的最優(yōu)解。本文采用線性反投影(linear back projection,LBP)法重建結(jié)果作為x的初值。

式(13)關(guān)于λ的最小化問題是一個(gè)最小二乘問題。通過對λ求導(dǎo),得:

給定xk和zk,采用最速下降法求解式(14),迭代步驟如下。

由λ=λk+1和給定的x=xk,利用式(13)對z的導(dǎo)數(shù)更新向量z:

最后,利用λ=λk+1和z=zk+1更新x:

φ(x)與φd(λ)間的對偶間隙為:

由此,迭代終止條件設(shè)為:對偶間隙小于某一正數(shù)或解已經(jīng)穩(wěn)定,即Δ＜ε或‖xk+1-xk‖2/‖xk+1‖2＜θε,其中參數(shù)ε＞0,θ∈(0,1]。

2.3 迭代方法選取閾值

為提高重建圖像的空間分辨率,利用迭代方法自適應(yīng)設(shè)定閾值,將像素值小于閾值的像素置零,以消除重建圖像中的偽影。

①迭代次數(shù)k=0,設(shè)置初始閾值T0為圖像全部像素的平均值。

②將重建圖像像素分為兩部分:Xa和Xb。其中, Xa代表值大于Tk的像素,Xb代表其他像素。

③分別計(jì)算Xa和Xb中像素的平均值:Ta=mean(Xa)、Tb=mean(Xb)。

3 試驗(yàn)驗(yàn)證

3.1 試驗(yàn)條件

試驗(yàn)用平面ECT系統(tǒng)如圖2所示。系統(tǒng)采用4× 4的矩形電極傳感器,電極尺寸為2 cm×2 cm,電極間距為0.5 cm。檢測對象為長、寬、高均為2 cm的正方體有機(jī)玻璃,介電常數(shù)約為3,分別放置在傳感器的不同位置,以測試ECT敏感場不同位置處的檢測精度。

3.2 試驗(yàn)結(jié)果

分別用Landweber算法、Tikhonov正則化算法和l1正則化算法進(jìn)行圖像重建,結(jié)果如圖2、圖3和圖4所示。

圖2 Landweber算法對不同模型的重建圖像Fig.2 The reconstructed images using Landweber algorithm for different models

圖3 Tikhonov正則化算法對不同模型的重建圖像Fig.3 The reconstructed images using Tikhonov regularization algorithm for different models

圖4 l1正則化算法對不同模型的重建圖像Fig.4 The reconstructed images using l1regularization algorithm for different models

由圖2可知,Landweber算法對敏感場中心和偏上側(cè)位置的物體均可較好地重建,但物體形狀有失真,邊緣處模糊情況較嚴(yán)重,重建物體的尺寸隨著物體放置位置的不同也發(fā)生變化;對兩個(gè)物體的重建圖像,可分辨出兩個(gè)物體,但形狀和尺寸的失真更為嚴(yán)重。由圖3可知,Tikhonov正則化算法對重建圖像質(zhì)量的改善不明顯,可探測出物體的存在和大致位置,但重建的形狀和尺寸存在比較嚴(yán)重的失真。由圖4可見,l1正則化算法有效改進(jìn)了成像質(zhì)量,提高了物體邊緣的空間分辨率,從而改善了平面ECT系統(tǒng)的檢測性能。

4 結(jié)束語

本文以平面ECT技術(shù)在航空復(fù)合材料缺陷檢測中的應(yīng)用為背景,利用l1正則化方法改進(jìn)ECT圖像重建算法。16電極平面ECT樣機(jī)試驗(yàn)結(jié)果表明,該方法明顯提高了重建圖像質(zhì)量,特別是敏感場內(nèi)物質(zhì)界面處的空間分辨率,提高了平面ECT系統(tǒng)的檢測精度,為平面ECT系統(tǒng)用于航空復(fù)合材料檢測提供了一定的試驗(yàn)論證。

[1] 楊乃賓.新一代大型客機(jī)復(fù)合材料結(jié)構(gòu)[J].航空學(xué)報(bào),2008,29(3): 596-603.

[2] 周正干,肖鵬,劉航航.航空復(fù)合材料先進(jìn)超聲無損檢測技術(shù)[J].航空制造技術(shù),2013(4):40-43.

[3] Xue Qian,Wang Huaxiang,Yang Chengyi,et al.Dynamical lag correlation exponent based method for gas-solid flow velocity measurement using twin-plane electrical capacitance tomography[J].Measurement Science and Technology,2012,23(8):085301-085307.

[4] 王挺.基于電容敏感機(jī)理的飛機(jī)復(fù)合材料檢測方法研究[D].天津:中國民航大學(xué),2013.

[5] Thiele S,Silva M J D,Uwe H.Capacitance planar array sensor for fast multiphase flow imaging[J].IEEE Sensors Journal,2009,9(5): 533-540.

[6] Chen D X,Hu X H,Yang W Q.Design of a security screening system with a capacitance sensor matrix operating in single-electrode mode[J].Measurement Science and Technology,2011,22(11): 114026-114032.

[7] Sun J T,Yang W Q.Fringe effect of electrical capacitance and resistance tomography sensors[J].Measurement Science and Technology,2013,24(7):074002.

[8] 張玲玲,王化祥,范文茹.基于1范數(shù)的電阻層析成像圖像重建算法[J].天津大學(xué)學(xué)報(bào),2011,44(9):786-790.

[9] Zhang Lingling,Xu Genqi,Xue Qian,et al.An iterative thresholding algorithm for the inverse problem of electrical resistance tomography[J]. Flow Measurement and Instrumentation,2013,33:244-250.

[10] Rudin L I,Osher S,Fatemi E.Nonlinear total variation based noise removal algorithms[J].Physica D,1992,60(1-4):259-268.

[11] Fang W F.A nonlinear image reconstruction algorithm for electrical capacitance tomography[J].Measurement Science and Technology, 2004,15(10):2124-2132.

[12] Soleimani M,Lionheart W R B.Nonlinear image reconstruction for electrical capacitance tomography using experimental data[J]. Measurement Science and Technology,2005,16(10):1987-1996.

[13] Wang H X,Tang L,Cao Z.An image reconstruction algorithm based on total variation with adaptive mesh refinement for ECT[J].Flow Measurement and Instrumentation,2007,18(5-6):262-267.

[14] Wang Y,Yang J,Yin W,et al.A new alternating minimization algorithm for total variation image reconstruction[J].SIAM Journal on Image Sciense,2008(1):948-951.

[15] Yang J F,Zhang Y.Alternating direction algorithms for l1-problems in compressive sensing[J].SIAM Journal on Scientific Computing, 2011,33(1):250-278.

l1Regularization Algorithm Used in Visualization Defect Inspection for Composite Materials

Along with the composite materials have been used more widely in civil aircrafts,the online inspection of the defects of composite materials requests higher level non-destructive inspecting technology.In order to improve the imaging quality of planar electrical capacitance tomography(ECT)system for enhancing the inspection accuracy for defects of composite materials,the imaging algorithm based on l1-regularization is researched,by applying alternating minimization method for numerical calculation,and combining with the iterative threshold method,the noise reduction is conducted for the reconstructed images.The experimental results show that comparing with the commonly used Landweber or Tikhonov regularization algorithms,this method can effectively distinguish different media in object field,and improve the imaging quality and the resolution for the edge of the measured object.

Tomography Non-destructive testing Image reconstruction Composite materials Regularization

TP207

A

國家自然科學(xué)基金青年科學(xué)基金資助項(xiàng)目(編號:61102096);

天津市科技支撐計(jì)劃基金資助項(xiàng)目(編號:11ZCGHHZ00700);

中央高?；究蒲袠I(yè)務(wù)費(fèi)中國民航大學(xué)專項(xiàng)基金資助項(xiàng)目(編號: 3122013C007);

中國民航大學(xué)科研啟動基金資助項(xiàng)目(編號:2013QD01S)。

修改稿收到日期:2014-04-02。

薛倩(1987-),女,2013年畢業(yè)于天津大學(xué)控制科學(xué)與工程專業(yè),獲博士學(xué)位,講師;主要從事無損檢測技術(shù)的研究。