鄧 娟 陳素華 沙 洪 趙 舒 任超世*

1（南昌大學(xué)信息工程學(xué)院電子信息工程系，南昌 330031）

2（中國(guó)醫(yī)學(xué)科學(xué)院＆北京協(xié)和醫(yī)學(xué)院生物醫(yī)學(xué)工程研究所，天津 300192）

引言

電阻抗斷層成像以人體組織、器官的阻抗分布或阻抗變化為依據(jù)，通過配置于人體體表的電極系統(tǒng)，提取與人體生理、病理狀態(tài)相關(guān)的組織或器官的電特性信息，給出反映組織、器官功能狀態(tài)及其變化規(guī)律的功能性圖像結(jié)果。國(guó)外的EIT研究一直很活躍，一些小組已由實(shí)驗(yàn)室研究向臨床應(yīng)用研究過渡［1－3］。國(guó)內(nèi)20世紀(jì)90年代末陸續(xù)進(jìn)行EIT系統(tǒng)研究，已建立多個(gè)EIT實(shí)驗(yàn)裝置和研究平臺(tái)，并在乳腺腫瘤檢測(cè)和腦、腹部滲血監(jiān)護(hù)和肺功能測(cè)量等方面向臨床應(yīng)用研究發(fā)展［4－7］。我國(guó)的阻抗成像研究起步較晚，與國(guó)外的差距依然存在，但現(xiàn)已具備相當(dāng)?shù)幕A(chǔ)，從方法學(xué)研究向應(yīng)用基礎(chǔ)和臨床應(yīng)用研究過渡的時(shí)機(jī)已經(jīng)成熟［8］。

EIT系統(tǒng)主要包括數(shù)據(jù)檢測(cè)系統(tǒng)和圖像重構(gòu)算法兩個(gè)部分，兩者相互依賴，互相補(bǔ)充。以往的EIT方法學(xué)研究旨在探索EIT成像原理，希望影響因素最少，通常假設(shè)EIT檢測(cè)系統(tǒng)是理想的，具有最高的檢測(cè)靈敏度和最佳的信噪比。然而，實(shí)用化EIT測(cè)量與理想仿真時(shí)的條件差別很大。比如，由于檢測(cè)系統(tǒng)有限的測(cè)量精度和動(dòng)態(tài)范圍，測(cè)量數(shù)據(jù)精度和噪聲的關(guān)系十分重要；對(duì)于不同位置目標(biāo)的成像，檢測(cè)系統(tǒng)信噪比對(duì)不同激勵(lì)-測(cè)量模式和圖像重建算法應(yīng)用效果的發(fā)揮和影響也有差異。

Rafiei-Naeini和McCann認(rèn)為，提高EIT檢測(cè)系統(tǒng)的信噪比，可以有效改善系統(tǒng)的檢測(cè)靈敏度和空間分辨率，而其中電流激勵(lì)子系統(tǒng)信噪比的提高對(duì)整個(gè)系統(tǒng)信噪比的提高將有很大貢獻(xiàn)。他們采用模擬和數(shù)字電子技術(shù)，設(shè)計(jì)了一種信噪比為80 dB低噪聲的電流激勵(lì)源子系統(tǒng)。該系統(tǒng)輸出阻抗達(dá)到10 MΩ，激勵(lì)頻率在10 kHz～4 MHz之間變化時(shí)產(chǎn)生的漂移較小［9］。Alder等認(rèn)為，一個(gè)優(yōu)化的EIT檢測(cè)系統(tǒng)應(yīng)對(duì)場(chǎng)域電導(dǎo)率變化有較好的分辨率，而這個(gè)分辨率可以表示為EIT系統(tǒng)的電極數(shù)和激勵(lì)電流幅度、興趣區(qū)電導(dǎo)率變化量、目標(biāo)位置和大小、激勵(lì)測(cè)量模式和測(cè)量噪聲的函數(shù)。該分辨率正比于與阻抗變化幅度，反比于噪聲幅度［10］。鑒于EIT對(duì)測(cè)量邊界的噪聲十分敏感，F(xiàn)ranciane通過在仿真數(shù)據(jù)中加入高斯白噪聲，研究了相鄰激勵(lì)和相對(duì)激勵(lì)的抗噪性能，發(fā)現(xiàn)相鄰激勵(lì)無(wú)法較好地解決含噪信號(hào)圖像重建的問題，而相對(duì)激勵(lì)則能較好地重建且成像誤差較小［11］。Wang等提出，在小電壓信號(hào)測(cè)量上，采取比大信號(hào)花費(fèi)更長(zhǎng)的采集時(shí)間來(lái)進(jìn)行噪聲補(bǔ)償?shù)奶幚矸椒ǎ?2］。羅辭勇等提出，采用一種交叉測(cè)量協(xié)議，不單純采用相鄰測(cè)量模式，通過增大測(cè)量電極間隔，獲取測(cè)量數(shù)據(jù)中較大的電壓進(jìn)行圖像重建，并從提高EIT檢測(cè)信號(hào)的抗噪性能、改善靈敏度和降低病態(tài)性等方面表明了交叉測(cè)量模式的實(shí)用化效果［13］。在算法方面，Kao等提出Tikhonov和Noser類組合正則化算法，降低噪聲影響的同時(shí)準(zhǔn)確定位目標(biāo)［14］。

鑒于國(guó)內(nèi)外EIT研究正從方法學(xué)研究向應(yīng)用基礎(chǔ)和臨床應(yīng)用研究過渡，為針對(duì)具體應(yīng)用目標(biāo)構(gòu)建EIT檢測(cè)系統(tǒng)，以及評(píng)價(jià)實(shí)用化EIT系統(tǒng)的性能，在不同位置單目標(biāo)和多目標(biāo)EIT成像中，探索不同信噪比的實(shí)用化EIT系統(tǒng)對(duì)不同的重建算法的影響以及相應(yīng)圖像重建效果的定量評(píng)價(jià)十分必要。

1 方法

筆者通過仿真的方法模擬不同信噪比的實(shí)用化EIT檢測(cè)系統(tǒng)和不同位置的成像目標(biāo)，比較不同算法的成像結(jié)果和性能，因此構(gòu)建了基于Matlab的EIT-GUI軟件研究平臺(tái)。應(yīng)用這一研究平臺(tái)，可進(jìn)行自適應(yīng)剖分，EIT問題正向計(jì)算和多種算法的圖像重建。圖像重建過程中，剖分層數(shù)、電極數(shù)目和激勵(lì)-測(cè)量模式可選。假設(shè)EIT圖像重建的是圓形場(chǎng)域，應(yīng)用平臺(tái)進(jìn)行7層三角剖分，并在場(chǎng)域表面安置16電極，如圖1所示。其中，場(chǎng)域表面編號(hào)1～16位置處用于放置電極，圖1中表示采用電極對(duì)（1，2）進(jìn)行相鄰激勵(lì)，采用電極對(duì)（3，4）進(jìn)行相鄰電壓測(cè)量，這里假設(shè)電極為點(diǎn)電極。

圖1 EIT場(chǎng)域電極位置及激勵(lì)－測(cè)量方式Fig.1 The electrodes position and excitementmeasurement pattern of EIT field

1.1 不同信噪比實(shí)用化EIT系統(tǒng)的模擬

針對(duì)16電極EIT系統(tǒng)，通過仿真獲取EIT場(chǎng)域邊界電壓作為理想EIT系統(tǒng)的采集信號(hào)。由于測(cè)量數(shù)據(jù)眾多，考慮所有檢測(cè)數(shù)據(jù)的影響，將所有數(shù)據(jù)檢測(cè)通道視為一個(gè)整體，對(duì)理想系統(tǒng)所有通道采集的數(shù)據(jù)整體添加高斯白噪聲，模擬一定信噪比的EIT檢測(cè)系統(tǒng)的測(cè)量電壓。SNR定義為

式中，i=1，2，…，n，n表示一次圖像重建所用的測(cè)量數(shù)據(jù)。voi是理想仿真條件下獲得的邊界電壓信號(hào)，vfi是仿真條件含有高斯分布噪聲的仿真檢測(cè)信號(hào)。模擬SNR分別為40、60、80 dB的實(shí)用化EIT系統(tǒng)。

1.2 EIT圖像重建算法

采用5種圖像重建算法，分別為反投影算法（lines back projection，LBP）、基于牛頓－拉夫遜（Newton-Raphson）算法的Tikhonov正則化算法（NR Tikhonov）和Tikhonov-Noser組合正則化算法（NR Tikhonov-Noser），以及Landweber二階迭代算法（Landweber）和Newton-Raphson算法。

LBP是基于反投影理論的算法。Landweber二階迭代法基于靈敏度系數(shù)理論，是最速下降法的一種變形，以迭代的方法逼近靈敏度系數(shù)矩陣的逆陣［15］。牛頓-拉夫遜算法是無(wú)約束極小化的一種迭代算法，它使用二次函數(shù)近似代替目標(biāo)函數(shù)［16］。為解決EIT圖像重建的不適定性問題，可以進(jìn)行正則化。Tikhonov正則化算法采用單位矩陣，作為正則化矩陣，該方法可以很好地去除噪聲，但是由于給系數(shù)矩陣加入了單位矩陣容易產(chǎn)生偽像［17］。Cheney等采用對(duì)角矩陣進(jìn)行正則化，稱為NOSER類正則化，該方法可提供正確的目標(biāo)位置，但重建時(shí)噪聲過大［18］。Kao等人于2006年提出Tikhonov和Noser類組合正則化算法，組合正則化方程為

g=（STS+εT×I+εN×diag（STS））－1STz（2）式中，εN代表Noser類正則化參數(shù)，而εT代表Tikhonov正則化參數(shù)；S代表靈敏度系數(shù)矩陣，而I代表單位矩陣；diag（STS）表示由正定陣STS生成的對(duì)角矩陣。

1.3 EIT圖像質(zhì)量評(píng)價(jià)方法

研究表明，同一算法用于不同的成像目標(biāo)（大小或位置不同），不同算法用于同一成像目標(biāo)，其重建結(jié)果在圖像質(zhì)量、信息量上有較大差異［19］。醫(yī)學(xué)圖像質(zhì)量的評(píng)價(jià)一般分為主觀評(píng)價(jià)方法和客觀評(píng)價(jià)方法。主觀評(píng)價(jià)方法常根據(jù)事先規(guī)定的評(píng)價(jià)尺度或自己的經(jīng)驗(yàn)，按視覺效果對(duì)待測(cè)圖像質(zhì)量進(jìn)行判斷，易受到觀察者的知識(shí)背景、觀測(cè)目的和環(huán)境等影響，難以用數(shù)學(xué)模型表達(dá)，因此無(wú)法廣泛推廣應(yīng)用?？陀^質(zhì)量評(píng)價(jià)方法依據(jù)數(shù)學(xué)模型給出的量化指標(biāo)衡量醫(yī)學(xué)圖像質(zhì)量，彌補(bǔ)了主觀評(píng)價(jià)方法的不足，已經(jīng)成為醫(yī)學(xué)圖像質(zhì)量評(píng)價(jià)研究的重點(diǎn)［20］。建立客觀、定量的EIT圖像評(píng)價(jià)方法，也是目前EIT研究和應(yīng)用中急需解決的問題之一。EIT通過有限元剖分方法求解，各剖分單元的電導(dǎo)率仍然是通過數(shù)值填充來(lái)進(jìn)行顯示的。對(duì)于這種填充的圖像，引入圖像重建誤差函數(shù)ER定量評(píng)價(jià)圖像重建誤差，用圖像的結(jié)構(gòu)相似度函數(shù)SSIM，定量比較和評(píng)價(jià)重建圖像對(duì)目標(biāo)分辨和定位的能力［21－22］。

1.3.1 圖像重建均方誤差ER

定義EIT圖像重建誤差函數(shù)ER為式中，M為有限元剖分的總單元數(shù)，GX（p）為成像模型各單元電導(dǎo)率，GY（p）為重建圖像各單元電導(dǎo)率，ER表示圖像重建前后每個(gè)單元的重建圖像電導(dǎo)率和模型電導(dǎo)率的誤差。

通過ER可獲得重建圖像和成像模型之間的差別，可靈敏地反映重建圖像的質(zhì)量。ER值越小，表示成像結(jié)果與模型的差別越小，圖像重建質(zhì)量越高。

1.3.2 結(jié)構(gòu)相似度函數(shù)SSIM

定義圖像X，Y的結(jié)構(gòu)相似度為［23］

式中，uX，uY，σXY分別表示已有電導(dǎo)率分布的模型圖像X和重建的電導(dǎo)率分布圖像Y的均值、方差和協(xié)方差，圖像的均值、方差（標(biāo)準(zhǔn)差）可以看作是對(duì)圖像亮度、對(duì)比度的粗略估計(jì)，而協(xié)方差則是重建圖像相對(duì)于成像模型的非線性改變。因此，l（X，Y），c（X，Y），s（X，Y）可以分別看作亮度比較、對(duì)比度比較和圖像結(jié)構(gòu)比較［23］。在圖像中，這3項(xiàng)是獨(dú)立的，單純的亮度、對(duì)比度改變并不會(huì)影響第3項(xiàng)。SSIM為

SSIM（X，Y）∈［0，1］，當(dāng) EIT重建圖像與檢測(cè)模型圖像完全相同時(shí)，SSIM（X，Y）=1，即結(jié)構(gòu)相似度SSIM的值越接近于1，成像質(zhì)量越好。

2 不同信噪比條件下不同算法的仿真成像

仿真實(shí)驗(yàn)中，采用16電極EIT系統(tǒng)，激勵(lì)電流為1mA，均勻場(chǎng)電導(dǎo)率設(shè)置為0.1 Ω－1m－1，成像目標(biāo)電導(dǎo)率設(shè)置為0.01 Ω－1m－1。由于場(chǎng)域中不同位置的成像目標(biāo)成像結(jié)果的差異較大，設(shè)放置在成像場(chǎng)域中不同位置的單個(gè)目標(biāo)模型為A、B、C，兩個(gè)目標(biāo)模型為D和4個(gè)目標(biāo)模型為E。圖2～圖4分別為在40、60、和80 dB信噪比EIT系統(tǒng)條件下，采用不同算法對(duì)不同目標(biāo)的仿真成像結(jié)果。表1～表3為采用圖像重建誤差函數(shù)ER和結(jié)構(gòu)相似度函數(shù)SSIM對(duì)圖2～圖4的定量評(píng)價(jià)結(jié)果。

圖2 40 dB EIT系統(tǒng)下不同算法的圖像重建結(jié)果Fig.2 The reconstruction results of different algorithms when the SNR is 40dB

表1 40 dB EIT系統(tǒng)下用不同算法重建圖像的定量評(píng)價(jià)Tab.1 The quantitative evaluation of reconstruction image of different algorithms when the SNR is 40 dB

圖2的成像結(jié)果表明，在信噪比為40 dB EIT測(cè)量系統(tǒng)條件下，對(duì)于場(chǎng)域中心的目標(biāo)A，僅LBP和NR Tikhonov-Noser可以成像。對(duì)所設(shè)其他4個(gè)目標(biāo)模型成像，5種算法成像效果的差異較大。采用ER和SSIM定量評(píng)價(jià)的結(jié)果表明，單目標(biāo)成像中，當(dāng)目標(biāo)由場(chǎng)域中心向測(cè)量淺表移動(dòng)（依次是目標(biāo)A、B、C），5種算法的SSIM都增大，這提示算法對(duì)目標(biāo)的定位能力增大。當(dāng)目標(biāo)位置都靠近淺表但個(gè)數(shù)增加時(shí)（依次為單目標(biāo)C、雙目標(biāo)D和4個(gè)目標(biāo)E），5種算法的ER都呈現(xiàn)出增大的趨勢(shì)，這表明圖像重建的誤差隨成像目標(biāo)個(gè)數(shù)的增加而增大。

圖3 60 dB EIT系統(tǒng)下不同算法的圖像重建結(jié)果Fig.3 The reconstruction results of different algorithms when the SNR is 60 dB

表2 60 dB EIT系統(tǒng)下不同算法重建圖像的定量評(píng)價(jià)Tab.2 The quantitative evaluation of reconstruction image of different algorithms when the SNR is 60 dB

圖4 80 dB EIT系統(tǒng)下不同算法的圖像重建Fig.4 The reconstruction results of different algorithms when the SNR is 80 dB

通過定量評(píng)價(jià)結(jié)果對(duì)5種成像算法的性能進(jìn)行分析發(fā)現(xiàn)，無(wú)論目標(biāo)位于何處，NR Tikhonov-Noser算法都可以有效成像，且成像結(jié)果都具有最大的SSIM，除在對(duì)目標(biāo)B的成像中其ER略大于LBP外，還具有最小的ER，其成像效果最好。LBP算法的ER僅次于NR Tikhonov-Noser算法，提示其成像誤差較?。坏瑫r(shí)其SSIM值也是最小，表明其重建圖像對(duì)目標(biāo)的分辨和定位能力較差，雖然能識(shí)別目標(biāo)，但往往有干擾目標(biāo)識(shí)別的成像噪聲存在。盡管如此，LBP能對(duì)場(chǎng)域中心的目標(biāo)A成像，這一優(yōu)勢(shì)也是除NR Tikhonov-Noser算法外其他3種算法所不能及的。NR Tikhonov、Landweber和Newton-Raphson算法的成像效果接近，其ER和SSIM值互相交叉，它們的成像效果明顯不如NR Tikhonov-Noser正則化算法的成像效果。對(duì)于它們可成像的目標(biāo)（除目標(biāo)A外），其成像誤差ER比LBP大，但結(jié)構(gòu)相似度SSIM也大，各有可取之處。

EIT測(cè)量系統(tǒng)信噪比達(dá)到60 dB時(shí)，圖3的結(jié)果表明僅NR Tikhonov和Landweber二階迭代算法不能對(duì)場(chǎng)域中心的單目標(biāo)A成像，5種算法對(duì)各個(gè)目標(biāo)的成像效果相對(duì)于40 dB時(shí)也有較大改善。表3采用函數(shù)ER和SSIM定量評(píng)價(jià)，發(fā)現(xiàn)目標(biāo)由場(chǎng)域中心向測(cè)量淺表移動(dòng)時(shí)5種算法的SSIM都增大。當(dāng)目標(biāo)個(gè)數(shù)增加，5種算法的ER都呈現(xiàn)出增大的趨勢(shì)，這與40 dB時(shí)類似。根據(jù)表3的評(píng)價(jià)結(jié)果來(lái)比較5種算法的性能，發(fā)現(xiàn)NR Tikhonov-Noser算法的優(yōu)勢(shì)仍然較為明顯，對(duì)所有成像目標(biāo)的圖像重建結(jié)果都有最小的ER，除目標(biāo)B的SSIM略小于Landweber算法外還具有最大的SSIM，成像效果是最好的，尤其是靠近淺表的目標(biāo)C、D、E。NR Tikhonov、Landweber和Newton-Raphson算法的成像效果相近，成像和定量評(píng)價(jià)的結(jié)果表明，它們的成像質(zhì)量相對(duì)于40 dB時(shí)改善較大，提示EIT系統(tǒng)信噪比的提高對(duì)這3種正則化算法成像效果的改善具有較好的作用。在40 dB信噪比測(cè)量系統(tǒng)下，LBP算法可對(duì)所有目標(biāo)成像且其ER較小，但在60 dB下優(yōu)勢(shì)已不明顯。

表3 80 dB EIT系統(tǒng)下用不同算法重建圖像的定量評(píng)價(jià)Tab.3 The quantitative evaluation of reconstruction image of different algorithms when the SNR is 80dB

圖4和表3的結(jié)果表明，當(dāng)系統(tǒng)信噪比達(dá)到80 dB時(shí)，5種算法均可對(duì)設(shè)定的5個(gè)目標(biāo)有效成像。其中，NR Tikhonov和Landweber算法獲得了最為明顯的改善，對(duì)中心目標(biāo)A實(shí)現(xiàn)了有效成像。相對(duì)于60 dB下的圖像重建，NR Tikhonov-Noser算法的重建效果略有提高，改善不大。對(duì)于同一目標(biāo)的成像，相對(duì)于60 dB時(shí)LBP的ER明顯減小，SSIM增大，重建效果獲得了較大的改善。Newton-Raphson的圖像重建也獲得了更優(yōu)的結(jié)果。提示具有80 dB信噪比這一參數(shù)的EIT系統(tǒng)足以保證本研究采用的5種算法對(duì)不同成像目標(biāo)有好的成像效果。對(duì)于NR Tikhonov-Noser算法，再繼續(xù)提高EIT系統(tǒng)的信噪比性能，通過算法來(lái)獲得成像效果改善的可能性不大。

3 討論與結(jié)論

筆者立足于實(shí)用化EIT測(cè)量與理想的仿真成像的差異，仿真研究了不同位置單個(gè)目標(biāo)和多個(gè)目標(biāo)的成像中信噪比對(duì)5種圖像重建算法性能的影響，發(fā)現(xiàn)EIT系統(tǒng)的信噪比不足將影響成像算法的應(yīng)用效果。在40、60、80 dB信噪比下，對(duì)不同算法的圖像重建結(jié)果和評(píng)價(jià)結(jié)果綜合分析如下：

1）無(wú)論系統(tǒng)信噪比是40、60還是80 dB，在單目標(biāo)成像中，當(dāng)目標(biāo)由場(chǎng)域中心向測(cè)量淺表移動(dòng)（依次是目標(biāo)A、B、C），5種算法的仿真成像圖與目標(biāo)模型圖的結(jié)構(gòu)相似度都增大，提示算法對(duì)目標(biāo)的定位能力增大；當(dāng)目標(biāo)位置都靠近淺表但個(gè)數(shù)增加時(shí)（依次為單目標(biāo)C、雙目標(biāo)D和4個(gè)目標(biāo)E），5種算法的ER都呈現(xiàn)出增大的趨勢(shì)，表明圖像重建的誤差隨成像目標(biāo)個(gè)數(shù)的增加而增大。

2）在5種算法中，NR Tikhonov-Noser算法表現(xiàn)出較大的優(yōu)勢(shì)。一方面，體現(xiàn)在該算法可在信噪比較差的EIT系統(tǒng)下有效成像；另一方面，在同樣的信噪比EIT系統(tǒng)和相同成像目標(biāo)條件下，該算法的成像質(zhì)量更高。究其原因，主要是算法加入了單位矩陣提高了抗噪性，加入了對(duì)角矩陣能減少單位陣帶來(lái)的偽像，提高了算法對(duì)目標(biāo)的分辨和定位能力，從而提高了成像質(zhì)量。

3）對(duì)于場(chǎng)域中心的目標(biāo)A，Newton-Raphson算法在40 dB信噪比時(shí)無(wú)法成像，NR Tikhonov，Landweber在40和60 dB時(shí)均不能獲得有效圖像。究其原因，一方面，是目標(biāo)靠近場(chǎng)域中心，形狀也不是完全的中心對(duì)稱，而這些算法都采用了靈敏度系數(shù)矩陣，由靈敏度系數(shù)理論可知，場(chǎng)域邊界靈敏度較高且電極附近靈敏場(chǎng)變化劇烈，而中心區(qū)域靈敏度低且變化緩慢；另一方面，NR Tikhonov-Noser算法也基于靈敏度系數(shù)矩陣，而40 dB即可對(duì)A成像，這三種算法抗噪性不及NR Tikhonov-Noser算法；再者，實(shí)驗(yàn)中將所有通道的測(cè)量數(shù)據(jù)視為一個(gè)整體加入高斯白噪聲，也可能是導(dǎo)致無(wú)法成像的一個(gè)因素。事實(shí)上，由于方法和實(shí)驗(yàn)實(shí)際情況的差異，在研究過程中研究者們加入和定義信噪比的方法不盡相同。如羅辭勇等人在研究不同測(cè)量方法的抗噪性的過程中，就定義信噪比為最大測(cè)量電壓的信噪比，他們?cè)?0 dB時(shí)即獲得了對(duì)場(chǎng)域中心的圓形規(guī)則目標(biāo)的有效成像結(jié)果［13］。這些因素和影響有待于實(shí)際EIT檢測(cè)系統(tǒng)實(shí)驗(yàn)的進(jìn)一步驗(yàn)證。

4）不同的圖像重建算法和不同成像目標(biāo)對(duì)成像系統(tǒng)的信噪比要求不一致。本仿真結(jié)果可指導(dǎo)實(shí)用化EIT系統(tǒng)的構(gòu)建過程中系統(tǒng)信噪比指標(biāo)的設(shè)計(jì)，應(yīng)根據(jù)其應(yīng)用對(duì)象（如檢測(cè)目標(biāo)位置）和擬采用的算法，避免造成信噪比不足無(wú)法有效成像，或避免過高的信噪比給硬件設(shè)計(jì)帶來(lái)難度。對(duì)于NR Tikhonov-Noser算法，40 dB的信噪比即可對(duì)場(chǎng)域內(nèi)目標(biāo)成像，60 dB的信噪比足夠獲得場(chǎng)域中心處較高質(zhì)量的圖像。Newton-Raphson算法在60 dB信噪比的系統(tǒng)中時(shí)可對(duì)場(chǎng)域中心目標(biāo)成像，而對(duì)于NR Tikhonov和Landweber算法80 dB才能實(shí)現(xiàn)這一目標(biāo)。對(duì)于已經(jīng)構(gòu)建好的系統(tǒng)，可根據(jù)實(shí)際成像結(jié)果估計(jì)其信噪比，從而對(duì)一定的成像目標(biāo)選擇更適合的圖像重建算法。NR Tikhonov-Noser算法可作為一般情況下的首選算法。除NR Tikhonov-Noser算法外，系統(tǒng)信噪比在40 dB左右或者40～60 dB時(shí)可選擇LBP。低于60 dB的系統(tǒng)，對(duì)于場(chǎng)域中心目標(biāo)的成像，不建議選擇NR Tikhonov、Landweber和Newton-Raphson算法。在80 dB時(shí)，雖然LBP算法相對(duì)于40和60 dB時(shí)自身性能有所改善，但相對(duì)于其他算法成像效果不夠理想，所以不建議選擇，除此之外的其余4種算法均可選。

EIT方法學(xué)研究采用理想化的測(cè)量系統(tǒng)，然而實(shí)用化EIT系統(tǒng)的信噪比總是有限的。在EIT檢測(cè)系統(tǒng)信噪比一定的情況下，系統(tǒng)往往缺乏必要的數(shù)據(jù)分辨能力，影響算法應(yīng)用效果的發(fā)揮。不同的算法由于各自的原理和抗噪性能不同，在不同的信噪比條件下表現(xiàn)出不同的性能，得到不同的圖像重建結(jié)果。在構(gòu)建EIT檢測(cè)系統(tǒng)時(shí)，信噪比指標(biāo)應(yīng)根據(jù)其應(yīng)用目標(biāo)和擬采用的算法進(jìn)行設(shè)計(jì)。已經(jīng)構(gòu)建好一定信噪比的系統(tǒng)，應(yīng)根據(jù)具體成像目標(biāo)，選擇更適合的圖像重建算法。

本研究?jī)H探索了不同的EIT檢測(cè)系統(tǒng)信噪比對(duì)不同圖像重建算法應(yīng)用效果的影響。在實(shí)際情況下，除了噪聲引起的誤差，電極間不匹配也會(huì)引起的誤差；實(shí)際電極也不是點(diǎn)電極，而是具有一定大小的面，電極與皮膚的接觸阻抗也將引起誤差；電極位置移動(dòng)，也會(huì)給成像引入偽差；檢測(cè)系統(tǒng)靈敏度的不足，也會(huì)影響到圖像重建的效果，等等。這些都是EIT走向?qū)嵱没倪^程中必須解決的關(guān)鍵問題，有待于深入研究。

［1］Trokhanova OV，Okhapkin MB，Korjenevsky AV.Dualfrequency electrical impedance mammography for the diagnosis of non-malignant breast disease［J］.Physiol Meas，2008，29（6）：S331－S344.

［2］Demidenko E.An analytic solution to the homogeneous EIT problem on the 2D disk and its application to estimation of electrode contact impedances［J］.Physiol Meas，2011，32（9）：1453－1471.

［3］Chakraborti KL，Selvamurth W.Clinical application of Electrical Impedance Tomography in the Present Health Scenario of India［J］.Journal of Physics：Conference Series，2010，224（1）：69－72.

［4］徐桂芝，王明時(shí)，李有余，等.醫(yī)用電阻抗成像系統(tǒng)的模塊化設(shè)計(jì)［J］.天津大學(xué)學(xué)報(bào)，2006，39（S1）：133－137.

［5］Sha Hong，Wang Yan，Zhao Shu，et al.Integrated data collection in electrical impedance tomography［C］.//Scharfetter H，Merva R，eds.IFMBE Proceedings.Berlin：Springer，2007：348－351.

［6］史學(xué)濤，霍旭陽(yáng)，尤富生，等.顱內(nèi)出血電阻抗成像系統(tǒng)及初步動(dòng)物實(shí)驗(yàn)［J］.航天醫(yī)學(xué)與醫(yī)學(xué)工程，2007，20（1）：24－27

［7］劉丹丹，烏日?qǐng)D，王超，等.基于醫(yī)學(xué)阻抗技術(shù)的乳腺癌檢測(cè)方法［J］.電子測(cè)量技術(shù)，2008，31（4）：60－64.

［8］任超世，李章勇，王妍，等.電阻抗斷層成像應(yīng)用基礎(chǔ)與臨床應(yīng)用的一些研究進(jìn)展［J］.中國(guó)生物醫(yī)學(xué)工程學(xué)報(bào)，2010，29（2）：300－304.

［9］Rafiei-Naeini M，McCann H.Low-noise current excitation subsystem for medical EIT ［J］.Physiol Meas，2008，29（6S），173－184.

［10］Adler A，Gaggero P，Maimaitijiang YS.Distinguishability in EIT using a hypothesis-testing model［J］. Journal of Physics：Conference Series，2010，224（1）：56 － 59.

［11］Franciane CP，Luis PSB.Somenumerical results on the inuence of measurement strategies and load patterns in the EIT inverse problem ［J］.Journal of Physics：Conference Series，2010，224（1）：145－148.

［12］Wang W，Brown BH，Leatilard AD，et al，Noise equalisation within EIT images［J］.Physiol Meas，1995，15（S2A）：211 －216.

［13］羅辭勇，陳民鈾，王平.電阻抗成像交叉測(cè)量模式的抗噪性能研究［J］.儀器儀表學(xué)報(bào)，2009，1（1）：1－6.

［14］Kao TJ，Isaacson D，Newell JC，et al.A 3D reconstruction algorithm for EIT using a handheld probe for breast cancer detection［J］.Physiol Meas，2006，27（5）：1 － 11.

［15］Yang WQ，Spink DM，York TA，et al.An image reconstruction algorithm based on Landweber iteration method for electrical capacitance tomography［J］.Meas Sci Technol，1999，10（11）：1065－1069.

［16］Yorkey TJ，Webster JG，Tompkins WJ.An improved peurtbration technique for electrical impedance imaging with some criticisms［J］.IEEE Trans Biomed Eng，1987，BME － 34（11）：898 －901.

［17］Vauhkonen M，Vadasz D，Kaejalainen PA，et al.Tikhonov regularization andpriorinform ationin electricalimpedance tomography［J］.IEEE Trans on Medical Imaging，1998，17（2）：285－293.

［18］Cheney M，Isaacson D，Newell JC，et al.NOSER：an algorithm for solving the inverse conductivity problem.International Journal of Imaging Systems and Technology［J］.1990，2（2）：66 － 75.

［19］董秀珍，秦明新，付峰，等.生物電阻抗斷層成像技術(shù)的研究進(jìn)展［J］.第四軍醫(yī)大學(xué)學(xué)報(bào)，1999，20（3）：252－255.

［20］段影影，馬建華，陳武凡，等.改進(jìn)的結(jié)構(gòu)相似醫(yī)學(xué)圖像質(zhì)量評(píng)價(jià)方法［J］.計(jì)算機(jī)工程與應(yīng)用，2010，46（2）：145－149.

［21］呂婧華，宮偉彥，王妍，等.基于梯度誤差電阻抗斷層成像技術(shù)的圖像評(píng)價(jià)［J］.中國(guó)組織工程研究與臨床康復(fù)，Vol 12，2008（52）：10277－10280.

［22］Wang Yan， ShaHong， RenChaoshi. Optimum designof electrode structure and parameters in electricalimpedance tomography［J］.Physiol Meas，2006，27（3）：291 – 306.

［23］Wang Z， Bovik AC， Sheikh HR， etal. Image quality assessment：from error measurement to structural similarity［J］.IEEE Trans on Image Processing，2004，13（4）：600 －612.