基于同步測量的電阻抗層析成像系統(tǒng)研究*

黃華芳，楊 濤，陳曉艷，楊永政

（天津科技大學電子信息與自動化學院，天津300222）

基于同步測量的電阻抗層析成像系統(tǒng)研究*

黃華芳，楊 濤，陳曉艷*，楊永政

（天津科技大學電子信息與自動化學院，天津300222）

針對APT模式EIT系統(tǒng)存在的異步測量、結構復雜、調試困難等問題，研究設計了一種基于同步測量的電阻抗層析成像系統(tǒng)，從硬件和軟件兩個方面闡述了如何實現同步測量。系統(tǒng)硬件由PXI集成控制器、開關模塊、數據采集模塊和自制電流源組成，簡化了系統(tǒng)結構，提高了系統(tǒng)的穩(wěn)定性和精確度。本系統(tǒng)采用LabVIEW開發(fā)了專用軟件包，實現了信號的同步采集，并用改進正則化算法重建了電導率分布圖像。對比實驗顯示了同步測量和異步測量的差異，結果表明同步測量EIT系統(tǒng)圖像更清晰，偽影更少，平均相對誤差降低了48.99%。本文提供了一種可行的EIT系統(tǒng)同步測量方法，該方法適用于多電極（電極數≥16）場域的無損傷成像應用。

電阻抗層析成像；同步測量；正則化算法；PXI；LabVIEW

電阻抗層析成像技術（Electrical Impedance Tomography，EIT）是繼形態(tài)、結構成像之后，出現的一種新興的功能性成像技術，電阻抗層析成像系統(tǒng)通常由數據采集系統(tǒng)和成像系統(tǒng)組成［1］。其基本原理是根據不同場域具有不同電阻抗這一物理原理，通過給場域注入交流激勵電流/電壓，測量邊界電壓信息，利用圖形重建算法來重構場域內部的電導率分布或其變化的圖像［2］。EIT的最大優(yōu)勢是能夠持續(xù)地反映場域內部結構及狀態(tài)的變化，是一種正處于迅速發(fā)展中并有廣泛應用前景的圖像重建技術［3］，具有無損傷、功能性成像、設備簡單、使用方便等特點［4-5］。

目前，大多數電阻抗成像系統(tǒng)都采用APT工作模式，邊界測量信號是按照順/逆時針順序依次采集及處理［6-9］，這種異步測量會導致測量數據的解調結果不準確，成像效果不理想。針對上述問題，本文旨在研究能實現同步測量、結構簡單、配置靈活的新型EIT系統(tǒng)。擬采用PXI（PCI Extensions for In?strumentation）平臺構建EIT系統(tǒng)，開發(fā)LabVIEW專用軟件，實現同步測量，并驗證該方法的可行性與優(yōu)越性。

1 理論基礎

1.1 正交序列解調算法

對EIT系統(tǒng)而言，當正弦電流注入激勵電極時，測量電極接收到的是交變的電壓信號，必須通過解調才能得到阻抗的實部和虛部信息，一般采用正交序列解調算法［10］進行解調。設

式中，p（n）是激勵電流信號，作為正交參考信號，u（n）是測量電壓信號，q（n）是同相參考信號，A為測量電壓幅值，φ為測量電壓相角，N為一個信號周期的采樣點數，n為第幾個采樣點（0≤n≤N-1）。根據正交序列解調算法，可計算測量電壓的實部分量R和虛部分量I。

進而得到成像系統(tǒng)所需的幅值信息

以上算法中的測量電壓相角φ表征的是由介質引起的測量信號與激勵信號之間的相位差。在異步測量時，除第一個測量電極的電壓信號能較準確表征該相位差外，隨后測量的電壓信號中都會因時間滯后引入相角誤差，這就導致正交序列解調算法得到的實部分量和虛部分量不準確，最終影響成像質量［11］。電極數越多，異步測量引入的相角誤差越大。針對多電極EIT系統(tǒng)（電極數≥16），采用同步測量可解決上述問題。

1.2 改進正則化算法

EIT系統(tǒng)完成邊界信息采集后，需要對被測對象的電導率分布進行圖像重建，本系統(tǒng)采用改進正則化算法［12］進行逆問題的求解。

電阻抗層析成像逆問題的數學模型可以表示為

式中，U為被測對象的邊界電壓測量值；J為靈敏度矩陣；σ為被測對象的電導率分布矩陣。

采用最小二乘法進行逆問題求解，

取JTJ的每一行向量或列向量的最大元素構成對角陣Q作為阻尼項，對JTJ施加阻尼作用，改進正則化算法如式（9）所示。

式中，α為正則系數，取值為0～1之間，具體值可通過建模、多元擬合和求最優(yōu)值的方法獲得。

2 系統(tǒng)硬件設計

本系統(tǒng)由被測對象（含16個電極）、激勵源模塊、切換開關模塊、數據采集模塊、控制器模塊和上位計算機組成，系統(tǒng)結構如圖1所示。上位機通過PCI總線與PXI-1033內置控制器連接，對每塊PXI板卡進行控制。首先PXI-5404產生正弦電壓信號，經VCCS（Voltage Control Current Source）電路轉化為電流信號后作為系統(tǒng)激勵源；該信號經由PXI-2593所控激勵電極注入被測對象，再由PXI-2520所控測量電極同時將14路電壓信號送至NI 5751和PXI-7952R進行數據采集；采集到的電壓數據通過PXI總線經由PXI-1033傳至上位計算機；最后，在上位機軟件LabVIEW中進行數據處理，并調用圖像重建程序對被測對象的電導率分布進行重構。

圖1 電阻抗層析成像系統(tǒng)結構圖

2.1 激勵源模塊

本系統(tǒng)采用電流激勵-電壓測量工作方式。激勵源由PXI-5404板卡和VCCS電路組成。信號發(fā)生器PXI-5404具有1.07 μHz分辨率，可編程輸出頻率9 kHz～100 MHz、幅值0.5 V～1 V的正弦波電壓信號，經16位運算放大器AD8021構建的VCCS電路可實現壓流轉換。VCCS電路原理如圖2所示，電壓信號由R1左端輸入，電流信號由R5右端輸出。本系統(tǒng)的激勵電流頻率為100 kHz，幅值為2 mA。

圖2 VCCS電路原理圖

2.2 三模態(tài)電極切換模塊

在EIT系統(tǒng)中，每個電極有三種工作模態(tài)，即電流激勵、電壓測量、接地，各電極均需在這三種模態(tài)下進行切換，選用PXI-2593板卡和PXI-2520板卡實現切換功能。為實現同步測量，本系統(tǒng)中每個電極由三條線纜連接，分別對應電極的激勵、測量和接地模態(tài)。PXI-2593多路復用器用于選通激勵電極流入端；PXI-2520通用開關用于選通激勵電極的流出端和測量端，其通道均由單刀單擲電樞繼電器構成，每個繼電器為非鎖存繼電器，具有極低的導通電阻，能有效隔離通道間干擾，切換時間為1 ms，滿足系統(tǒng)對測量通道一致性及同步測量的要求。三模態(tài)電極切換電路如圖3所示。

圖3 三模態(tài)電極切換電路

2.3 同步采集模塊

鑒于系統(tǒng)整體性能與數據采集速度（精度）要求，選取NI 5751板卡作為適配器、PXI-7952R板卡作為FPGA模塊，共同構成數據采集模塊。

NI 5751是模數轉換適配器，可實現16通道同步采樣，采樣速率50 MS/s，輸入電壓范圍2 VPP，分辨率14位。經LabVIEW編程可實現FPGA模塊PXI-7952R對適配器NI 5751的控制，實現多路信號同步采集，克服常規(guī)依次采集所造成的測量精度下降問題，保證了系統(tǒng)的成像質量。

3 系統(tǒng)軟件設計

本系統(tǒng)采用LabVIEW 2013和LabVIEW 2013 FPGA完成軟件開發(fā)工作，主要包括三部分：用戶界面GUI設計、數據采集和圖像重構。

3.1 GUI設計

在LabVIEW平臺上設計了用戶界面，實現參數設置、數據采集、曲線繪制和系統(tǒng)成像的功能。在參數設置界面，可設置PXI-5404的地址、激勵源幅值及頻率，PXI-2593和PXI-2520的地址、工作模式，PXI-7952R FPGA的地址等參數。曲線繪制界面可實現對實時采集數據的曲線繪制和歷史數據的曲線繪制，界面如圖4所示。系統(tǒng)成像界面可選擇圖像重建算法及對應算法參數。

圖4 LabVIEW曲線繪制界面

3.2 數據采集程序

數據采集程序流程如圖5所示。首先進行各功能板卡的初始化，設好地址及工作模式；接著進入選通開關程序，LabVIEW控制PXI-2593選通激勵電極，控制PXI-2520選通測量和接地電極；然后對14路測量電壓信號進行同步采集，數據處理后存在上位機中；之后選通下一組激勵、測量和接地電極，繼續(xù)測量下一組數據，直到16組激勵電極均已選通，便完成了一幅圖像重構所需的所有數據采集。重復以上測量過程，即可得到下一幅圖像數據。

圖5 軟件控制流程圖

3.3 圖像重構程序

本系統(tǒng)采用LabVIEW調用MATLAB程序的方法，進行被測對象的電導率分布圖像重構。MATLAB程序流程如圖6所示，在載入靈敏度矩陣和空場、滿場數據并導入正則系數之后，先對采集的數據進行預處理，用空場數據減去滿場數據得到相對數據，再用改進正則化算法進行逆問題的求解，接著對求解數據進行歸一化處理，最后繪制重建圖像。

在LabVIEW中進入系統(tǒng)成像界面，選擇使用正則化算法，設定正則系數α的值后即可得到重構的電導率分布圖像。

圖6 圖像重構程序流程圖

4 實驗及結果

系統(tǒng)實物如圖7所示。被測對象選用直徑20 cm、高25 cm的圓柱水槽，水槽內盛有電導率0.332 S/m、溫度15℃的生理鹽水，電極選用直徑4 mm的金屬鈦圓形點電極。

圖7 基于同步測量的電阻抗層析成像系統(tǒng)

4.1 系統(tǒng)性能測試

本系統(tǒng)電極數16個，激勵源頻率100 kHz，激勵源幅值2 mA，采樣頻率5 MHz，測量電壓范圍-1 V～+1 V，根據數據采集模塊分辨率可知能分辨的最小測量電壓為［1-（-1）］/214=0.122 mV，測量精度達±0.2‰。

對系統(tǒng)信噪比SNR（Signal to Noise Ratio）進行測試，按式（10）計算

4.2 靜態(tài)實驗

在水槽中依次放入直徑5 cm、4 cm、3 cm三種尺寸的有機玻璃棒，分別標記為粉色、藍色和黑色以示區(qū)分。在LabVIEW中選擇正則化算法進行成像，正則系數α=0.7，與實驗室前期研制的基于FPGA的異步測量EIT系統(tǒng)［13］進行對比實驗，成像效果如表1所示。

表1 兩套系統(tǒng)成像效果對比

從表1可以看出，本系統(tǒng)更好地反映了有機玻璃棒的形狀、大小及位置，聚焦性增強了，偽影減少了，成像質量明顯提高。

選取圖像相對誤差RE（Relative Error）作為重建圖像質量的衡量指標，按式（11）計算

式中，M為重建圖像的總單元數；σ（j）為重建電導率分布；σ*（j）為仿真電導率分布［12］。計算結果如表2所示。

表2 兩套系統(tǒng)圖像相對誤差對比

從表2可以看出，三種情況下同步測量系統(tǒng)的圖像相對誤差都比異步測量系統(tǒng)小，平均相對誤差降低了48.99%。

4.3 動態(tài)實驗

經實際測試，本系統(tǒng)采集速度為7 frame/s，將直徑5 cm的有機玻璃棒放在電極8和電極9之間，在距離水槽圓心5 cm處順時針勻速旋轉一周，用時約14 s，采集到100幅圖像數據，選取其中8幅，動態(tài)成像效果如圖8（a）-（h）所示。

圖8顯示的是第5幅、第18幅、第31幅、第44幅、第57幅、第70幅、第83幅、第96幅的重建圖像，這些圖像如實地反映了有機玻璃棒的運動過程，說明本系統(tǒng)的成像結果能夠很好地識別有機玻璃棒的運動軌跡。

圖8 動態(tài)成像圖

5 結論

本文設計了基于同步測量的電阻抗層析成像系統(tǒng)，實現了對被測對象的電導率分布圖像重構。PXI平臺極大地簡化了硬件系統(tǒng)，縮短了硬件平臺的搭建時間；三模態(tài)電極和同步采集模塊的使用實現了同步測量，提高了系統(tǒng)的采集和成像精度；在軟件方面，LabVIEW極大地提高了軟件的開發(fā)效率，縮短了程序開發(fā)周期。由于PXI平臺的可擴展性，方便今后搭建更加復雜的系統(tǒng)。本系統(tǒng)集成度高，穩(wěn)定性好，測量準確，抗干擾能力強，擴展方便，為無創(chuàng)醫(yī)學成像和工業(yè)成像領域提供了理想的解決方案，尤其適合于多電極（電極數≥16）場域的無損傷成像應用。

［1］Barber D C，Brown B H.Applied Potential Tomography［J］.Jour?nal of Physics E：Scientific Instruments，1984，17（9）：723-733.

［2］董秀珍.生物電阻抗成像研究的現狀與挑戰(zhàn)［J］.中國生物醫(yī)學工程學報，2008，27（5）：641-643.

［3］孫進平，程吉寬，杜巖.電阻抗成像的修正阻尼最小二乘重建算法［J］.電路與系統(tǒng)學報，2000，5（2）：47-51.

［4］Brown B H.Electrical Impedance Tomography（EIT）：A Review［J］.Journal of Medical Engineering and Technology，2003，27（3）：97-108.

［5］任超世，李章勇，王妍，等.我國電阻抗層析成像實用化應用研究展望［J］.航天醫(yī)學與醫(yī)學工程，2010，23（4）：309-312.

［6］Goharian M，Soleimani M，Jegatheesan A，et al.A DSP Based Multi-Frequency 3D Electrical Impedance Tomography System［J］.Annals of Biomedical Engineering，2008，36（9）：1594-1603.

［7］張凌烽，王化祥.基于FPGA和DSP的電阻層析成像數據采集系統(tǒng)［J］.傳感技術學報，2011，24（7）：1011-1016.

［8］Gaggero P O，Adler A，Brunner J，et al.Electrical Impedance Tomography System Based on Active Electrodes［J］.Physiological Measurement，2012，33（5）：831-847.

［9］侯海嶺，王化祥，陳曉艷.基于多周期欠采樣技術的生物電阻抗成像硬件系統(tǒng)［J］.傳感技術學報，2015，28（2）：211-216.

［10］趙德春，任超春，沙洪，等.EIT高精度數字解調方法誤差分析［J］.儀器儀表學報，2010，31（9）：1933-1938.

［11］陳曉艷，楊永政，杜萌.無解調電阻抗成像方法研究［J］.傳感技術學報，2015，28（6）：870-875.

［12］陳曉艷，房曉東.一種新的正則化圖像重建算法及參數優(yōu)化［J］.天津科技大學學報，2014，29（6）：74-77.

［13］Chen X Y，Yang T，Yang Y Z.A PXI-Based Biomedical Electrical Impedance Tomography System［C］//Applied Mechanics and Materials，2014，670-671：1205-1209.

黃華芳（1979-），女，副教授，碩士生導師，天津科技大學電子信息與自動化學院，主要研究方向為測控技術與電力電子技術，huanghf@tust.edu.cn；

楊 濤（1990-），男，碩士研究生，現就讀于天津科技大學電子信息與自動化學院，主要研究方向為電阻抗成像技術，13820772636@163.com；

陳曉艷（1973-），女，教授，博士，天津科技大學電子信息與自動化學院，主要研究方向為電學參數檢測與處理，cxywxr@tust.edu.cn；

楊永政（1988-），男，碩士研究生，現就讀于天津科技大學電子信息與自動化學院，主要研究方向為電阻抗成像技術，yyz1988@163.com。

Study on Electrical Impedance Tomography System Based on Synchronous Measurement*

HUAN Guafang，YANG Tao，CHEN Xiaoyan*，YANG Yongzheng
（College of Information and Automation，Tianjin University of Science and Technology，Tianjin 300222，China）

Avoiding the disadvantages of traditional APT-pattern electrical impedance tomography（EIT）systems，such as asynchronous measurement，complex structure and debugging difficulty，a synchronous measuring EIT sys?tem is proposed and designed.This paper introduces an approach of realizing synchronous measuring in terms of hardware and software.The hardware platform consists of PXI-based integrated controller module，switch module，data acquisition module and lab-made current generator.Thishardware scheme simplifiesthe structure and improves the stability and accuracy of the system.A specified software kit is developed by LabVIEW to carry out data syn?chronous collecting and images reconstructing.An improved regularization algorithm is applied to reconstruct the in?ternal conductivity distribution.The comparative experiments are taken to illustrate the differences between syn?chronous and asynchronous measuring.The results show that the images are clearer and the artifacts are less.The average relative error of EIT images based on synchronous measuring is decreased by 48.99%.This paper provides a feasible and available method of synchronous measuring for EIT system，which is suitable for noninvasive imaging applications in numerous electrodes（electrode number≥16）fields.

electrical impedance tomography；synchronous measurement；regularization algorithm；PXI；LabVIEW

TP391

A

1004-1699(2015)10-1575-06

??7230

10.3969/j.issn.1004-1699.2015.10.027

項目來源：國家自然科學基金項目（61301246）；天津市高校優(yōu)秀青年教師資助計劃

2015-05-22 修改日期：2015-07-21