摘要
电阻抗断层成像技术(Electrical Impedance Tomography,EIT)是当今生物医学工程学重大研究课题之一,具有功能成像,无损伤和医学图像监护三大突出优势,是一种理想的、具有诱人应用前景的无损伤医学成像技术。
EIT电极系统直接与人接触,位于系统的最前端,是整个系统的关键部分之一。目前电阻抗断层成像技术现正处于一个从实验室研究阶段向临床应用过渡的时期,临床EIT直接与人体接触,电极系统面临一系列的问题需要研究和改进。而电极系统上发生的事件,包括有用信息、噪声、伪差、接触阻抗、极化电压等,都会进入后续电路被放大、传输,参与信号处理,影响图像重建结果,无用信息有时甚至会淹没有用信息,对系统检测灵敏度和图像重建质量的影响都非常大。场域边界移动形变造成的电极系统不确定性和电极排布方式也是影响成像结果的因素。根据实验室研究和临床研究的不同需求,研究电极特性,优化电极系统性能,具有十分重要的意义。
论文回顾了EIT研究历史,分析了EIT电极系统目前亟待解决的问题。在介绍EIT成像基本原理、硬件系统、成像算法的基础上,提出EIT图像重建质量的评价方法,建立了EIT成像软硬件平台和测试系统,比较了几种重建算法的成像效果及其局限。
基于强制等势点有限元模型,研究了电极类型、电极结构参数、电极数和激励测量模式对成像结果的影响,提出了电极系统优化设计方案。设计了一种用于EIT实验、参数可变的盐水槽电极系统,电极片拆卸方便,结构简单。
进行了EIT电极阻抗特性研究。提出了EIT电极阻抗特性的意义和重要性。对目前EIT人体实验通常采用的一次性心电电极进行了交直流阻抗特性测试和研究,提出了使用注意事项。
针对肺功能EIT检测时,由于电极运动和形变导致的图像伪差问题,采用改进的正则化算法和基于椭圆模型成像方法,改善图像质量。
设计了手持式平面电极系统,研究了平面电极阵列EIT成像方法。建立了3D有限元模型,应用基于等势面反投影法的改进的加权重建算法,仿真重建出了多深度层的断层图像,为EIT乳腺测量新方法奠定了基础。
Electrical impedance tomography (EIT) is one of the most important research topics in today's biomedical engineering. It has the following outstanding advantages: noninvasive, harmless, functional imaging and medical image monitoring. EIT is an attractive medical imaging technology with potential clinical application.
Electrode system lies in the front of the whole EIT system and contacts with human body directly. It is one of the most pivotal parts of EIT system. There are a lot of problems should be studied and improved in the electrode system of EIT. The useless information from electrode system, such as noise, contact impedance and polarization voltage etc, goes into the following circuits as useful signals. They are amplified, transferred and participate in signal processing process, and therefore influence the result of reconstruction finally. The uncertainty made by the change of boundary shape and electrode movement, as well as the arrange function of electrode system, are also affect the quality of reconstruction image greatly. These problems are even serious, as EIT is stepping from experimental study into clinical study currently. The study of both the property of electrode system and optimizing the electrode system is very important.
In this paper, firstly, the research history of EIT has been reviewed. Then the problems of current electrode systems have been presented. The principle of EIT, hardware system and image reconstruction method have been introduced. Several criteria functions used to evaluate the quality of reconstructed image have been proposed. The software, hardware and test platform have developed. Saveral reconstruction algorithms have been compared.
We studied the influence of electrode type, electrode structure parameter, electrode count and the exitation function on the quality of image reconstruction, based on the coercive equipotential node finite element model. Some optimizing design suggestion has been proposed. A new electrode install method with simple structure has been designed, which is convenient to remove and install the electrode slices. And the electrode parameter can be changed freely.
The impedance characteristic of EIT electrode, especially the amplitude and the stability of the alternating current impedance, is very important. Since the one-off ECG electrode has been widely used in clinical EIT measurement, we studied and compared more than ten ECG electrodes. Some suggestions need attention have been proposed.
To minish the image error made by electrode movement and boundary deform in EIT pulmonary function measurement, we used a modified regularization algorithm and an elliptical model based algorithm to improve the image quality.
An EIT method with planar electrode array has been proposed and a hand-held planar electrode system has been designed. Based on a 3-dimentional finite element model, a modified weighted backprojection algorithm, which is back-projected along equipotential surfaces of the 3D electric filed, is used to reconstruct slice images of 3D conductivity distribution, which is parallel to the electrode array. The method using planar electrode array establishs the foundation for EIT breast mesurement.
引文
[1]Barber D.C.and Brown B.H.,1984,Applied petential tomography,J.Phsy.E:Sci.instrum.,17:723-733
[2]任超世,生物电阻抗测量技术,中同医疗器械信息,2004,10(1):21-25。
[3]Henderson R P and Webster J G 1978 An impedance camera for spatially specific measurements of the thorax IEEE Trans.Biomd.Eng 25 250-254.
[4]A.M.Dijkstra,B.H.Brown,A.D.Leathard,et al,Review Clinical Application of Electrical Impedance Tomography,Journal of Medical Engineering & Technology,1993,17(3):89-98
[5]Kerner T E 2001 Electrical impedance tomography for breast imaging,Phd thesis,Dartmouth College.
[6]Frerichs I,Electrical impedance tomography(EIT) in applications related to lung and ventilation:a review of experimental and clinical activities,Physiol.Meas.200021:R1-21
[7]Yang W Q and York T,A new ac-based capacitance tomography system,Proc.IEE,1999,146:47-53.
[8]Jin Keun Seo,Ohin Kwon,et al,Megnetic resonance Electrical Impedance Tomography(MREIT),Ⅻ International Conference on Electrical Bio-lmpedance & Ⅴ Electrical Impedance Tomography.Gdansk,Poland,June 20-24,2004:19-22.
[9]David Holder.Electrical Impedance Tomography of Brain Function,Ⅻ International Conference on Electrical bio-lmpedance & Ⅴ Electrical Impedance Tomography,Gdansk,Poland,June 20-24 2004:7-12
[10]http://www.rpi.edu.
[11]史学涛 霍旭阳 尤富生 付峰 刘锐岗 徐灿华 董秀珍。颅内出血电阻抗成像系统及初步动物实验。航天医学与医学工程。2007,20(1):24-27
[12]廖琪梅,董秀珍,付峰,陈广生,阎艳,张聚良,凌瑞,漆家学,郭青霞。乳腺癌浸润部位对其电阻抗频谱特性影响的研究。航天医学与医学工程。2007,20(4):269-273
[13]王化祥 李璐 王春艳。用于肺功能监测的电阻抗成像系统。医疗卫生装备。2007(7):1-3
[14]盛剑宏,电磁场数值分析,北京:科学出版社,1984
[15]毕德显,电磁场理论,北京:电子工业出版社,1985
[16]李忠元,电磁场边界元素法,北京:北京工业出版社,1987
[17]倪光正,钱秀英,电磁场数值计算,北京:高等教育出版社,1992
[18]李开泰,有限元方法及应用,西安:西安交通大学出版社,1992
[19]曾余庚,电磁场有限单元法,北京:科学出版社,1982
[20]David S Holder,Electrical Impedance Tomography,Methods,History and Applications.IOP Publishing Ltd.2005.
[21]董秀珍,秦明新,汤孟兴,付峰,史学涛,尤富生,“一种实用的动态电阻抗断层成像算法”,中国生物医学工程学报,2000,19(3):pp.348-352
[22]王超,医学电阻抗断层成像的研究,大津大学博士论文,2001
[23]秦明新、李世俊、董秀珍、刘锐岗、付峰、尤富生、史学涛,“基于敏感性定理的EIT重构算法仿真研究”,生物医学工程学杂志,2002:19(2):pp.236-238
[24]刘锐岗,董秀珍,秦明新,尤富生,史学涛,付峰,汤孟兴,“基于物理模型的静态电阻抗断层成像初步结果”,第四军医大学学报,2001,22(4):pp.297-300
[25]田海燕,何为,杨浩,“电阻抗断层成像技术重建算法的研究”,2002年生物电磁学学术研讨会论文,pp.32-36
[26]李文波,王博亮,孙桂菊,向飞,刘希顺,“一种用于电阻抗成像的拟Newton 算法-基于修IEBk的开关算法”,北京生物医学工程,2000,19(1):PP:21-27
[27]侯卫东,莫玉龙,“电阻抗成象中求解反问题的新方法”,计算机工程,2001,27(5),pp.71-78
[28]魏颖,赵进创,土师,陆增喜,“基于RBF神经网络的电阻层析成像算法的研究”,仪器仪表学报,2001,22(4):pp.369-371
[29]A.1.Nachman,"Global uniqueness for a two-dimensional inverse boundary value problem," Ann.Math,vol.143,pp.71-96,1996
[30]Kolehmainen V,Lassas M and Ola P 2006 Inverse conductivity problem with an imperfectly known boundary SIAM J.Appl.Math.66(2) 365-383
[31]董秀珍,秦明新,刘锐岗等.影响生物电阻抗断层成像质量的因素[J].中国医学物理学杂忐,2001,18(4):209-210.
[32]Barber D C.EIT:the view from Sheffield.Electrical impdeance tomography methods,history and applications.UK:Bristol,lnst Phys Pub,2005:348-373.
[33]Avis N J,Barber DC.Adjacent or polar drive:image reconstruction implications in electrical impedance tomography systems employing filtered back projection[C].Paris:IEEE Press,1992:1689-1690.
[34]王磊,医用电极对心脑电图机测量的影响,实用医技杂志2007年7月第14卷第20期:2625-2626。
[35]郑万松,董秀珍,尤富生,李正,张国鹏,赵大庆.单个电极接触不良对腹部电阻抗断层成像的影响[J],第四军医大学学报,2006年02期,86-88。
[36]Assenheimer M et al.2001.The T-scan technology:electrical impedance as a diagnostic tool for breast cancer detection.Physiol.Meas.22,1-8。
[37]Kerner T E 2001 Electrical impedance tomography for breast imaging,Phd thesis,Dartmouth College.
[38]Wei Wang,M Tang,J Yang et al,Electrical impedance mammography a new imaging modal ity for breast cancer screening.Proceedings of the 9th Annual Conference of LSSCB Dec 2000 Leicester UK 58-59.
[39]Wang Huaxiang,Wang Chao,Yin Wuliang,Optimum design of the structure of the electrode for a medical EIT system,Meas.Sci.Technol,2001,12:1020-1023。
[40]尤富生,董秀珍,史学涛,付峰,刘锐岗,李江,电阻抗断层成像中临近和对向驱动模式的研究,第四军医大学学报,2004,25(1),pp88-91
[41]Tin-Su Pan,Andrew E.Yagle,Numerical study of multigraid implementation of some interative image reconstruction algorithms,IEEE Tansactions on medical imaging,1991,10(4):843-852
[42]Tin-Su Pan,Andrew E.Yagle,acceleration of Landweber-type algorithms by suppression of projection on the maximum singular vector,IEEE Transactions on Medical imaging,1992,11(4):479-487
[43]Polydorides N,Image reconstruction algorithms for soft-field tomography,Ph.D Thesis,UMIST,England,2002.
[44]Tikhonov A N and Arsenin V Y,Solutions of ill-posed problems,WH Winston,1977.
[45]Claude Cohen-Bacrie,Yves Goussard,Regularized Reconstruction in Electrical Impedance Tomography Using a Variance Uniformization Constraint,IEEE Transactions on Medical Imaging,1997,16(5):562-571.
[46]Yang Wuqiang,Peng Lihui,Image reconstruction algorithms for electrical capacitance tomography,Measurement Science and Technology,2003,14:1-13
[47]唐磊,电学成像系统图像重建算法研究及可视化软件设计,博士学位论文,天津大学,2008.
[48]丁永祥,夏巨愖,有限元分析过程中的等值线图生成算法,模具技术,1994,3:47-52
[49]Frerichs I,Electrical impedance tomography(EIT) in applications related to lung and ventilation:a review of experimental and clinical activities,Physiol.Meas.2000 21:R1-21
[50]Kolehmainen V,Vauhkonen M,Karjalainen P A and Kaipio J P,Assessment of errors in static electrical impedance tomography with adjacent and trigonometric current patterns Physiol.Meas.1997(18):289-303
[51]Harris N D,Suggett A J,Barber D C and etal.Applied potential tomography:a new technique for monitoring pulmonary function,Clin.Phys.Physiol.Meas,1988(9):79-85.
[52]Lionheart W R B,Boundary shape and electrical impedance tomography Inverse Problems 1998 14:139-147
[53]Gersing E,Hofman B,Osypka M.Influence of changing peripheral geometry on electrical impedance tomography measurements.Med.Biol.Eng.Comput,199634:359-361.
[54]Kolehmainen V,Voutilainen A and Kaipio J P Estimation of non-stationary region boundaries in EIT—state estimation approach Inverse Problems 2001 171937-1956
[55]Kolehmainen V,Lassas M and Ola P Inverse conductivity problem with an imperfectly known boundary SIAMJ.Appl.Math.2006 66(2):365-383
[56]Soleimani M,Gomez—Laberge C,Adler A.Imaging of conductivity changes and electrode movement in EIT,Physiol.Meas,2006 27:103-113.
[57]Michael J.Young著,邱仲潘等译,Visual C++6从入门到精通,北京:电子工业出版社,1999.
[58]马艺馨,电阻层析成象技术及其在气/液两相泡状流检测中的应用,1999年。
[59]Zhang J and Patterson R P.EIT images of ventilation:what contributes to the resistivity changes? Physiol.Meas.2005,26:81-92
[60]Harris N D,Suggett A J,Barber D C and etal.Applied potential tomography:a new technique for monitoring pulmonary function,Clin.Phys.Physiol.Meas,1988(9):79-85.
[61]李现勇等,Visual C++串口通信技术与工程实践,北京人民邮电出版社,2004.[62]Wang Yan,Sha Hong and Ren Chaoshi,Optimum design of electrode structure and parameters in electrical impedance tomography.Physiol.Meas.27(2006)291-306.
[63]王妍,沙洪,任超世。电阻抗断层成像线电极的仿真研究,生物医学工程学杂志,2007,V01.24(5):986-989。
[64]Barber D C and Brown B H.,Errors in reconstruction of resistivity images using a linear reconstruction technique.Clin.Phys.Physiol.Meas.1988,9:101-4
[65]Adler A,Guardo R and Berthiaume Y 1996 Impedance imaging of lung ventilation:do we need to account for chest expansion? IEEE Trans.Biomed.Eng.43 414-20
[66]Katusyuki Sakamoto,Hiroshi Kanai,The measurement of lung impedance by impedance plethysmography,Proceedings of the Vth ICEBI Aug,1981,Tokyo,pp 361-364
[67]V Cherepenin,A Karpov,A Korjenevsky,V Kornienko,A Mazaletskaya,D Mazourov and D Meister.A 3D electrical impedance tomography(EIT) system for breast cancer detection.Physio.Meas.22(2001)9-18.
[68]何世钧,电容层析成像系统的研究与应用,博士学位论文,天津大学,2005。
[69]Ping Hua et al,Using compound electrodes in electrical impedance tomography,IEEE Trans.Biome.Engn.,Vol.40,No.1,1993,pp 29-34
[70]Blott B H,Daniell G J and Meeson S 1998 Electrical impedance tomography with compensation for electrode positioning variations Phys.Med.Biol.43 1731-9
[71]一次性使用心电电极,中华人民共和国医药行业标准,YY/T 0196-2005。
[72]Wang Zhou,Bovik A C,Sheikh H R,et al.Image Quality Assessment:From Error Measurement to Structural Similarity.IEEE Transactions on Image Processing,2004,13(4):600-613.
[73]Wang Zhou,Bovik A C.A universal image quality.IEEE Signal Processing Letters,2002,9(3):81-84.