莫仲念，周輝林，段榮行，王玉皞

(南昌大學(xué)信息工程學(xué)院電子信息工程系，江西南昌 330031)

?

基于DSM-CSI的非線(xiàn)性逆散射算法研究

莫仲念，周輝林，段榮行，王玉皞

(南昌大學(xué)信息工程學(xué)院電子信息工程系，江西南昌 330031)

由于反演問(wèn)題中散射體所處區(qū)域和散射體個(gè)數(shù)信息的缺失，進(jìn)而導(dǎo)致非線(xiàn)性逆散射方法的待重構(gòu)空間維數(shù)高和運(yùn)算代價(jià)高.基于線(xiàn)性逆散射方法在重構(gòu)目標(biāo)區(qū)域和目標(biāo)個(gè)數(shù)信息的計(jì)算代價(jià)低，本文提出一種DSM(Direct Sampling Method)線(xiàn)性方法和CSI(Contrast Source Inversion)非線(xiàn)性方法相結(jié)合的非線(xiàn)性逆散射方法.該方法利用DSM方法的重構(gòu)結(jié)果為CSI方法提供散射體所處區(qū)域和散射體個(gè)數(shù)的先驗(yàn)信息，進(jìn)而降低待重構(gòu)空間的維數(shù)以及非線(xiàn)性逆散射問(wèn)題的計(jì)算量.實(shí)驗(yàn)計(jì)算結(jié)果證明了該方法的有效性.

逆散射；直接抽樣方法；對(duì)比源反演

1 引言

電磁場(chǎng)逆散射方法是利用測(cè)量散射場(chǎng)數(shù)據(jù)和電磁場(chǎng)前向模型，重構(gòu)目標(biāo)的幾何形狀或電性能參數(shù)，近年來(lái)已經(jīng)廣泛應(yīng)用于目標(biāo)識(shí)別、生物醫(yī)學(xué)成像、遙感、地球物理成像、無(wú)損測(cè)試和評(píng)價(jià)等領(lǐng)域[1～3].大多數(shù)電磁場(chǎng)逆散射算法在某種形式上都依賴(lài)于散射體內(nèi)外場(chǎng)的Lippmann-Schwinger方程或者域積分方程[4].由于這兩個(gè)方程所具有的非線(xiàn)性和病態(tài)性，使其求解和計(jì)算過(guò)程極具挑戰(zhàn).

針對(duì)非線(xiàn)性問(wèn)題導(dǎo)致計(jì)算復(fù)雜度高的問(wèn)題，近年來(lái)提出一種基于玻恩近似[5,6]的線(xiàn)性化方法，該方法假設(shè)散射體為弱散射體，其隱含兩個(gè)條件：(1)目標(biāo)與背景介質(zhì)的對(duì)比度小；(2)目標(biāo)尺寸小于波長(zhǎng).在該假定條件下，目標(biāo)區(qū)域的總場(chǎng)近似等于入射場(chǎng).在此基礎(chǔ)上運(yùn)用TSVD和Tikhonov等正則化方法求解該線(xiàn)性病態(tài)方程.但線(xiàn)性電磁場(chǎng)逆散射成像方法只能定性實(shí)現(xiàn)散射體目標(biāo)的檢測(cè)定位和粗略地估計(jì)其幾何形狀，不能定量重構(gòu)散射體的電性能參數(shù).而在某些應(yīng)用場(chǎng)景中除了需重構(gòu)目標(biāo)幾何形狀外，還需重構(gòu)目標(biāo)的電性能參數(shù)[7].

實(shí)際應(yīng)用中絕大多數(shù)散射體屬于強(qiáng)散射體，使得線(xiàn)性算法的假設(shè)不成立.隨著計(jì)算機(jī)性能和存儲(chǔ)性能的不斷提升，近年來(lái)非線(xiàn)性逆散射算法成為研究的熱點(diǎn).其求解思路將逆散射問(wèn)題轉(zhuǎn)化為優(yōu)化問(wèn)題，通過(guò)構(gòu)建目標(biāo)函數(shù)和迭代優(yōu)化方法實(shí)現(xiàn)目標(biāo)的幾何形狀和電性能參數(shù)非線(xiàn)性重構(gòu).非線(xiàn)性算法則有對(duì)比源反演(Contrast Source Inversion,CSI)[8]、子空間優(yōu)化方法(Subspace Optimization Method,SOM)[9]、玻恩迭代法(Born Iterative Method,BIM)[10]及其變形(Distorted Born Iterative Method,DBIM)[11]等.非線(xiàn)性逆散射算法需已知目標(biāo)所處區(qū)域和目標(biāo)個(gè)數(shù)的先驗(yàn)信息，但該先驗(yàn)信息通常是未知的，且重構(gòu)結(jié)果極大地受初始值的影響[7].有鑒于此，本文提出一種DSM[12]線(xiàn)性算法和CSI非線(xiàn)性算法相結(jié)合的混合成像方法.該方法利用DSM的快速運(yùn)算優(yōu)勢(shì)，實(shí)現(xiàn)快速目標(biāo)定性成像，重構(gòu)結(jié)果為CSI方法提供目標(biāo)個(gè)數(shù)及所處區(qū)域的先驗(yàn)信息，并將其作為CSI方法的初始假設(shè)，從而克服傳統(tǒng)CSI方法對(duì)先驗(yàn)信息的要求，及降低待重構(gòu)空間的維數(shù)，最終降低非線(xiàn)性逆散射問(wèn)題的復(fù)雜度和提升優(yōu)化迭代速度.

2 電磁場(chǎng)非線(xiàn)性方程

考慮一個(gè)二維成像幾何模型如圖1所示，假定背景介質(zhì)均勻無(wú)界，其中存在一個(gè)有界的非均勻目標(biāo)區(qū)域D，區(qū)域D中包含一個(gè)或多個(gè)散射目標(biāo)，散射目標(biāo)的位置和對(duì)比度函數(shù)χ(r)未知.在TM極化條件下，同一頻率發(fā)射天線(xiàn)和接收天線(xiàn)分別按圓形軌跡T和R等間隔分布，其位置矢量分別為rj(j=1,2,…,Nj)和rq(q=1,2,…,NS).在目標(biāo)區(qū)域D內(nèi)的總場(chǎng)滿(mǎn)足積分方程

r∈D

(1)

(2)

其中ε(r′)和σ(r′)分別為散射目標(biāo)的介電常數(shù)和導(dǎo)電率.

在接收天線(xiàn)處接收到的散射場(chǎng)滿(mǎn)足積分方程

(3)

將目標(biāo)區(qū)域D劃分為M個(gè)面積相等的方形網(wǎng)格，網(wǎng)格中心的位置矢量為rm(m=1,2,…,M)，每個(gè)網(wǎng)格面積用ΔS表示.則積分方程(1)和(3)可用離散形式表示

m=1,2,…,M

(4)

rq=1,2,…，NS

(5)

將方程(4)用矩陣形式表示

同理將方程(5)寫(xiě)成矩陣形式

(7)

3 基于DSM-CSI的混合成像方法

基于線(xiàn)性算法的快速成像結(jié)果能夠滿(mǎn)足非線(xiàn)性算法對(duì)目標(biāo)區(qū)域和目標(biāo)個(gè)數(shù)先驗(yàn)信息的要求，本章將提出一種線(xiàn)性和非線(xiàn)性算法相結(jié)合的混合非線(xiàn)性成像算法.此算法能夠降低非線(xiàn)性算法求解逆散射問(wèn)題中未知網(wǎng)格信息個(gè)數(shù)，減少計(jì)算復(fù)雜度和降低計(jì)算代價(jià).

3.1 DSM方法

本文采用一種DSM線(xiàn)性算法直接估計(jì)散射體目標(biāo)形狀和位置，相比MUSIC(Multiple Signal Classification)[13,14]和LSM(Linear Sampling Method)[15]線(xiàn)性電磁場(chǎng)逆散射成像算法，該算法并不涉及任何的矩陣運(yùn)算，所以DSM算法運(yùn)算代價(jià)低.在DSM算法中，在第m個(gè)目標(biāo)網(wǎng)格處定義一個(gè)定性評(píng)價(jià)函數(shù)

(8)

式中，||·||表示取L2范數(shù)，〈·〉表示取內(nèi)積，假如Φj(m)≈1，則此網(wǎng)格在散射體內(nèi)；相反若Φj(m)≈0，此網(wǎng)格在散射體外.

3.2 CSI反演方法

(9)

同理方程(7)寫(xiě)為

(10)

目標(biāo)函數(shù)的構(gòu)建是逆散射算法的重要組成部分，定義目標(biāo)函數(shù)為

(11)

上式右端第一項(xiàng)為歸一化的場(chǎng)數(shù)據(jù)方程誤差，第二項(xiàng)為歸一化的場(chǎng)狀態(tài)方程誤差.該目標(biāo)函數(shù)是關(guān)于ωj的二次函數(shù)，但關(guān)于χ的非線(xiàn)性函數(shù).實(shí)際應(yīng)用表明，CSI算法是最有效的非線(xiàn)性逆散射算法之一.

3.3 DSM-CSI方法

本文結(jié)合以上兩種方法提出一種DSM-CSI方法，該方法利用DSM的重構(gòu)結(jié)果為CSI方法提供散射目標(biāo)個(gè)數(shù)及其所處區(qū)域信息，進(jìn)而減少待重構(gòu)空間χ維數(shù)，從而降低整個(gè)非線(xiàn)性算法的計(jì)算量.具體結(jié)合DSM-CSI算法流程如算法1所示.

算法1 DSM-CSI算法流程

1.利用DSM方法為CSI方法提供散射體先驗(yàn)信息

10.當(dāng)式(11)小于某個(gè)設(shè)定值時(shí)則終止迭代輸出χn，否則n=n+1，返回4.

4 實(shí)驗(yàn)結(jié)果分析

在本實(shí)驗(yàn)中散射體對(duì)比度實(shí)部即相對(duì)介電常數(shù)大致范圍已知，可由相對(duì)介電常數(shù)誤差來(lái)衡量反演重構(gòu)算法誤差.定義相對(duì)介電常數(shù)誤差

(12)

式中χmax為反演重構(gòu)對(duì)比度中最大值，χr為散射體真實(shí)介電常數(shù)最大值.當(dāng)相對(duì)介電常數(shù)誤差小于終止條件時(shí)，反演迭代終止.

在本節(jié)中，分別運(yùn)用DSM、CSI和DSM-CSI方法對(duì)散射體逆散射反演成像，并對(duì)重構(gòu)結(jié)果進(jìn)行對(duì)比分析.本實(shí)驗(yàn)中散射體的導(dǎo)電率未知，為驗(yàn)證算法的成像效果，本文中CSI和DSM-CSI算法的成像結(jié)果只給出相對(duì)介電常數(shù)分布圖并與散射體目標(biāo)的真實(shí)相對(duì)介電常數(shù)對(duì)比.

利用來(lái)自菲涅爾實(shí)驗(yàn)室的數(shù)據(jù)(twodiel-TM 8f.exp)[16]重構(gòu)散射體.圖2(a)為實(shí)驗(yàn)中散射體橫截面，圖中描述了散射體的實(shí)際個(gè)數(shù)、位置和形狀大小，兩個(gè)圓柱形散射體分別位于原點(diǎn)的兩側(cè)，半徑為15mm，圓心到原點(diǎn)距離為45mm，背景介質(zhì)為自由空間，散射體的相對(duì)介電常數(shù)在3±0.3之間.36根發(fā)射天線(xiàn)均勻分布于圓心在實(shí)驗(yàn)裝置中心的圓上，半徑為0.72m,49根接收天線(xiàn)分布于半徑為0.76m的圓上，在反演中選取頻率f=4GHz的數(shù)據(jù)重構(gòu)散射體參數(shù).

選取0.15m×0.15m的目標(biāo)區(qū)域，將其劃分為120×120個(gè)網(wǎng)格，每個(gè)網(wǎng)格的大小等于1.25mm×1.25mm，DSM方法重構(gòu)結(jié)果如圖2(b)所示，從圖中可以確定散射目標(biāo)所處區(qū)域和目標(biāo)個(gè)數(shù)并跟實(shí)際場(chǎng)景相符，但DSM方法并不能重構(gòu)散射體的對(duì)比度.利用CSI算法，將目標(biāo)區(qū)域剖分為40×40個(gè)網(wǎng)格，成像結(jié)果如圖2(c)所示，所重構(gòu)相對(duì)介電常數(shù)值跟實(shí)際介電常數(shù)值非常相近.運(yùn)用DSM-CSI方法，根據(jù)圖2(b)的成像結(jié)果，可將目標(biāo)區(qū)域縮小為兩個(gè)0.06m×0.06m的目標(biāo)區(qū)域，為保持與CSI方法網(wǎng)格尺寸相同，將每個(gè)目標(biāo)區(qū)域剖分為16×16個(gè)網(wǎng)格，反演結(jié)果如圖2(d)所示，散射體相對(duì)介電常數(shù)反演值也非常接近實(shí)際值.對(duì)比圖2(c)和圖2(d)，CSI算法和DSM-CSI算法的反演質(zhì)量相差不大.DSM-CSI方法和CSI方法的相對(duì)介電常數(shù)誤差與迭代次數(shù)如圖3所示，根據(jù)終止條件式(12)，可知DSM-CSI迭代次數(shù)比CSI迭代次數(shù)少得多.通過(guò)對(duì)比分析，運(yùn)用CSI方法需要重構(gòu)出40×40個(gè)網(wǎng)格信息，而DSM-CSI方法僅僅需要重構(gòu)16×32個(gè)網(wǎng)格信息，有效的減少了逆散射問(wèn)題的未知量個(gè)數(shù).DSM-CSI的核心算法是CSI算法，由于DSM-CSI方法減少了待重構(gòu)空間的維數(shù)，所以DSM-CSI方法的每次迭代時(shí)間比CSI方法少得多，因此DSM-CSI方法總的運(yùn)算時(shí)間大大少于CSI方法.

采用來(lái)自菲涅爾實(shí)驗(yàn)數(shù)據(jù)(FoamTwinDielTM)[17]，運(yùn)用三種不同的反演方法逆散射成像得出的結(jié)果如圖4所示.實(shí)驗(yàn)中散射體的橫截面如圖4(a)所示，一個(gè)半徑為15.5mm的圓柱體包含于一個(gè)半徑為40mm的圓柱體內(nèi)，而一個(gè)半徑同樣為15.5mm的圓柱體與大圓柱體相切.大圓柱體和小圓柱體的相對(duì)介電常數(shù)分別為1.45±0.15和3±0.3.在實(shí)驗(yàn)中，18根發(fā)射天線(xiàn)均勻分布于圓心在實(shí)驗(yàn)裝置中心的圓上，半徑為1.67m，241根接收天線(xiàn)分布于半徑為1.67m的圓上用于接收散射場(chǎng)數(shù)據(jù)，利用頻率f=5GHz的實(shí)驗(yàn)室數(shù)據(jù)重構(gòu)反射體目標(biāo)形狀、位置或?qū)Ρ榷确植?

假定已知散射體位于0.3m×0.3m的目標(biāo)區(qū)域內(nèi)，實(shí)際位置以及形狀未知.運(yùn)用DSM算法，成像結(jié)果如圖4(b)所示，可以看出DSM方法雖不能成功重構(gòu)散射體的形狀，但能精確重構(gòu)散射體所處區(qū)域.將0.3m×0.3m的目標(biāo)區(qū)域劃分為50×50個(gè)網(wǎng)格，每個(gè)網(wǎng)格為6mm×6mm，利用CSI方法成像如圖4(c)所示，可以看出三個(gè)散射體的位置及形狀大小，且成功反演三個(gè)散射體目標(biāo)的相對(duì)介電常數(shù).根據(jù)圖4(b)，DSM重構(gòu)結(jié)果為CSI方法重新選取目標(biāo)區(qū)域大小為0.12m×0.12m，將目標(biāo)區(qū)域剖分為20×20個(gè)網(wǎng)格，網(wǎng)格尺寸與圖4(c)相同，DSM-CSI反演成像如圖4(d)所示，DSM-CSI成像效果跟CSI方法一致.DSM-CSI方法和CSI方法的相對(duì)介電常數(shù)誤差與迭代次數(shù)如圖5所示，從圖中可知DSM-CSI方法的相對(duì)介電常數(shù)誤差隨著迭代次數(shù)增加而減少速度比CSI方法略快，在式(12)的終止條件下，兩種算法的迭代次數(shù)相差不大.相比CSI算法，DSM-CSI算法極大減少了逆散射問(wèn)題中的未知量個(gè)數(shù)，這極大降低了重構(gòu)問(wèn)題的運(yùn)算代價(jià)，從而降低運(yùn)算時(shí)間.在逆散射反演問(wèn)題中運(yùn)算時(shí)間是衡量一個(gè)算法優(yōu)越性的重要指標(biāo)，相比單一的CSI方法，DSM-CSI算法具有極大的速度優(yōu)勢(shì).

5 結(jié)論

本文所提出DSM-CSI方法利用DSM為CSI方法提供先驗(yàn)信息，在此基礎(chǔ)上重新選取目標(biāo)區(qū)域D，然后運(yùn)用CSI算法反演出散射體的相關(guān)參數(shù).數(shù)值實(shí)驗(yàn)結(jié)果證明，DSM-CSI方法和CSI方法的反演質(zhì)量一致，但DSM-CSI方法的計(jì)算代價(jià)大大低于CSI方法.

[1]Tuchin V V,Tuchin V.Tissue Optics:Light Scattering Methods and Instruments for Medical Diagnosis[M].Bellingham:SPIE Press,2007.

[2]Randazzo A,Estatico C.A regularisation scheme for electromagnetic inverse problems:application to crack detection in civil structures[J].Nondestructive Testing and Evaluation,2012,27(3):189-197.

[3]Colton D,Kress R.Inverse Acoustic and Electromagnetic Scattering Theory[M].Springer Science & Business Media,2012.

[4]Van Den Berg P M,Kleinman R E.A contrast source inversion method[J].Inverse Problems,1997,13(6):1607.

[5]Poli L,Oliveri G,Massa A.Microwave imaging within the first-order Born approximation by means of the contrast-field Bayesian compressive sensing[J].IEEE Transactions on Antennas and Propagation,2012,60(6):2865-2879.

[6]鄭虎,李廉林,李芳.玻恩近似下的無(wú)相位檢測(cè)電磁成像方法[J].電波科學(xué)學(xué)報(bào),2008,23(4):669-673.

Zheng Hu,Li Lianlin,Li Fang.EM inverse scattering from phase less data of total field under Born approximation[J].Chinese Journal of Radio Science,2008,23(4):669-673.(in Chinese)

[7]Catapano I,Randazzo A,Slob E,et al.GPR imaging via qualitative and quantitative approaches[A].Civil Engineering Applications of Ground Penetrating Radar[M/OL].Springer International Publishing,2015.239-280.

[8]Jeffrey I,Zakaria A,LoVetri J.Microwave imaging by mixed-order discontinuous Galerkin contrast source inversion[A].General Assembly and Scientific Symposium (URSI GASS)[C].USA:IEEE,2014.1-4.

[9]Oliveri G,Zhong Y,Chen X,et al.Multiresolution subspace-based optimization method for inverse scattering problems[J].JOSA A,2011,28(10):2057-2069.

[10]趙延文,聶在平.多重網(wǎng)格技術(shù)與波恩迭代法相結(jié)合的反演新方法[J].電子學(xué)報(bào),2003,31(2):171-174.

Zhao Yanwen,Nie Zaiping.A novel inversion method combining born iterative method with multi-grid technique[J].Acta Electronica Sinica,2003,31(2):171-174.(in Chinese)

[11]Ran-Duc T,Linh-Trung N,Do M N.Modified distorted born iterative method for ultrasound tomography by random sampling[A].International Symposium on Communications and Information Technologies (ISCIT)[C].USA:IEEE,2012.1065-1068.

[12]Ito K,Jin B,Zou J.A direct sampling method to an inverse medium scattering problem[J].Inverse Problems,2012,28(2):025003.

[13]樊祥,程正東,馬東輝,等.兩信號(hào)源的相關(guān)性對(duì)MUSIC算法分辨性能的影響[J].電子學(xué)報(bào),2008,36(12):2315-2318.

Fan Xiang,Cheng Zhengdong,Ma Donghui,et al.Effect of correlation of two signal sources on resolution performance of MUSIC algorithm[J].Acta Electronica Sinica,2008,36(12):2315-2318.(in Chinese)

[14]郭躍,王宏遠(yuǎn),周陬.陣元間距對(duì)MUSIC算法的影響[J].電子學(xué)報(bào)，2007，35(9):1675-1679.

Guo Yue,Wang Hongyuan,Zhou Zou.Effect of array elements spacing on MUSIC algorithm[J].Acta Electronica Sinica,2007,35(9):1675-1679.(in Chinese)

[15]Haddar H,Lechleiter A,Marmorat S.An improved time domain linear sampling method for Robin and Neumann obstacles[J].Applicable Analysis,2014,93(2):369-390.

[16]Belkebir K,Saillard M.Special section:Testing inversion algorithms against experimental data[J].Inverse Problems,2001,17(6):1565-1571.

[17]Geffrin J M,Sabouroux P,Eyraud C.Free space experimental scattering database continuation:experimental set-up and measurement precision[J].Inverse Problems,2005,21(6):S117.

莫仲念 男，1989年生，廣西人，現(xiàn)為南昌大學(xué)在讀碩士研究生，主要從事超寬帶穿墻雷達(dá)成像，逆散射成像方法研究.

周輝林 男，1979年生，江西人，現(xiàn)為南昌大學(xué)教授、理學(xué)博士，從事超寬帶雷達(dá)成像、雷達(dá)信號(hào)處理和目標(biāo)識(shí)別等領(lǐng)域的研究.

E-mail:zhouhuilin@ncu.edu.cn

段榮行 男，1977年生，江西人，現(xiàn)為南昌大學(xué)副教授，主要研究方向?yàn)閺?fù)雜系統(tǒng)可靠性分析、故障診斷及信息融合，主持和參與多項(xiàng)國(guó)家自然科學(xué)基金和省部級(jí)項(xiàng)目.

DSM-CSI Based Nonlinear Inverse Scattering Method forComplex Objects Reconstruction

MO Zhong-nian,ZHOU Hui-lin,DUAN Rong-xing,WANG Yu-hao

(DepartmentofElectronic&InformationEngineering,NanchangUniversity,Nanchang,Jiangxi330031,China)

Due to the absence of information concerning the location and number of scattering objects in inverse scattering problem,it results in nonlinear inverse scattering method calling for high refactoring space dimension and high computational cost.Based on the low computational burden with the linear inverse scattering methods in reconstructing the location and number of the targets,we propose a nonlinear inverse scattering method in this paper,which combined with DSM linear and CSI nonlinear inverse scattering methods.The reconstruction results of DSM are used as the prior information to provide the location and the number of the targets for the CSI method,which reduce the dimension of refactoring space and the amount of computation.Experimental results show the effectiveness of the proposed method.

inverse scattering problem;direct sampling method (DSM);contrast source inversion method (CSI)

2015-02-06；

2015-05-18；責(zé)任編輯：孫瑤

國(guó)家自然科學(xué)基金(No.61062009，No.71461021，No.61261010)；江西省自然科學(xué)基金(No.2015BAB207001)；江西省科技支撐計(jì)劃(No.20151BBE50090)

O451

A

0372-2112 (2016)10-2501-06

??學(xué)報(bào)URL:http://www.ejournal.org.cn

10.3969/j.issn.0372-2112.2016.10.031