1、一种生物电阻抗断层参数成像方法的建立作者：李镐炜，董秀珍，王跃科，史学涛，付峰 【关键词】 图像处理 【Abstract】 AIM: To reconstruct the electrical impedance parameter tomogram of tissues. METHODS: Backprojection was utilized to reconstruct the electrical impedance tomogram at every driving frequency before the ColeCole model of bioelectrical spectr
2、oscopy was analyzed and computed in complex plane. Then, image reconstruction was performed. RESULTS: The image of the dispersion parameter of ColeCole model was reconstructed using the data acquired from a 32electrode phantom. CONCLUSION: The reconstructed image of the scheme is acceptable and the
3、method used in this study can yield better images after some optimization. 【Keywords】 backprojection; ColeCole model; electric impedance; tomography, X ray computed; image processing, computerassisted 【摘要】 目的： 对组织进行生物电阻抗断层参数成像方法： 用反投影法重建出各激励频率下阻抗断层图像，在复平面上对生物阻抗谱的 ColeCole模型进行分析计算.随后进行图像重建 结果： 对 32电极
4、物理模型上采集的实验数据用此法分析,得到了 ColeCole模型散射系数的重建图像结论： 所述方法成像结果较准确，是一种可望经改进获得优良成像质量的参数成像方法 【关键词】 反投影;ColeCole 模型;电阻抗；体层摄影术，X 线计算机；图像处理，计算机辅助 0 引言 电阻抗断层成像（EIT）是在生物医学工程领域发展起来的一种无损检测与图像重建技术因其成本低廉、无辐射、可动态显示组织生理活动等优点，有广阔的应用前景人体不同组织和体液有不同的电参数，组织病变时，电阻率和介电常数均会变化无论健康普查或临床辅助诊断，生物电阻抗断层参数成像（EIPT）相对通常的 EIT技术，对生物体内阻抗信息挖掘的
5、更充分. 结合生物电阻抗的 ColeCole模型1求取组织内部参数时，可得到更为丰富的内部组织信息本文在 EIT原理和生物电阻抗 ColeCole模型的基础上，建立了一种生物电阻抗断层参数成像的方法. 1 原理 1.1EIT 原理生物电阻抗断层参数成像首先面临的是“逆”问题，同时需考虑组织阻抗的 ColeCole模型，通过边界测量数据求解生物组织内部的 ColeCole模型的关键参数生物组织的阻抗分布复杂，正问题与逆问题均超出仅用欧姆定律既能求解的范围，结合电磁场理论，可由偏微分方程中的椭圆方程边值问题描述2:=0(1)=(2)n=v(3)其中， 为复介电常数, 为未知电势分布函数,和 v分别
6、为边界区域的电压和边界电流密度， 为物体所占空间区域， 为其边界公式(1)、(2)构成已知全部边界场值的第一边界条件问题（Dirichlet 问题）,公式(1)、(3)构成已知全部边界的法向导数的第二边界条件问题（Neumann 问题） 反投影法3是一种针对 EIT逆问题较成熟的算法，并由此衍生出多种变体其基本原理是先用有限元法将被测组织离散化为多个单元，对每一单元进行分析；然后在不同电极上施加激励信号并测量其他电极电位，计算出不同电力线方向的阻抗数据；最后将与分析单元相关的阻抗数据按照线性叠加原理求和，总和即与该单元的阻抗值成正比这是分析生物阻抗的特征参数的基础. 1.2ColeCole 模
7、型根据 ColeCole理论1 ，生物组织的阻抗特性取决于细胞膜的低漏电特性所导致的容抗 Cm和细胞质与细胞间质所导致的电阻 Ri、Re,组织阻抗的经验公式是4： Z(f)=R+（R0-R）/1+(jf/fc)a (4)其中 R0Re 为低频时电阻，RRi/Re 为高频时电阻,fc1/2RiCm 为组织的特征频率, 为散射系数从复变函数理论可知,在复平面,Z(f)是处于第 4像限的一段圆弧，如 Fig 1所示. 2 参数成像方法 已知 3个频率下测量出生物组织的阻抗，结合组织电阻抗 ColeCole图可解出其相应的电阻抗参数 、R、R0、fc 等. 不妨记 ColeCole圆周上三点为： (f
8、i)Re(fi)jIm(fi) (i=1，2，3)(5)则 Cole圆的圆心 Y对应的复数 a+jb的实部与虚部经推导，可用(fi)的实部与虚部表示为：aIm(f1)-Im(f2)-Re2(f1)-Re2(f3) Im(f3)-Im(f1)-1-Re2(f2)-Re2(f3) Im(f2)-Im(f3)-1/2Re(f3)-Re(f2)Im(f2)-Im(f3)-1-Re(f1)-Re(f3) Im(f3)-Im(f1)-1 (6) bRe(f3)-Re(f2) Im(f2)-Im(f3)-1 Im(f1)Im(f3)-Re2(f1)-Re2(f3) Im(f3)-Im(f1)-1-Re(f1
9、)-Re(f3) Im(f3)-Im(f1)-1 Im(f2)Im(f3)+Re2(f2)-Re2(f3) Im(f2)-Im(f3)-1 /2Re(f3)-Re(f2) Im(f2)-Im(f3)-1-Re(f1)-Re(f3)Im(f3)-Im(f1)-1(7) 根据圆心 ajb 和(fi)即可进一步得出其他参数从 EIT图像重构的角度考虑，若要准确确定待测物体内部的复阻抗分布，需对测量数据进行静态成像重构，这对测量系统硬件和重构算法的要求都很高，因此需做变通处理. 如 Fig 1所示，圆周上四点 O，A， B， C，如已知以其中某一点 O为起始点，其余三点为终点的向量，同样可根据公式(6
10、)、(7)确定圆心坐标，进而在复平面唯一确定三向量终点所在的圆周方程故可以用 4个向量的 3个独立差向量求解圆周方程. 设待测组织在含有 4个频率的安全电流激励下的边界电压是(fi)(i=1,2,3,4),其中 f1f2f3f4对于 n个电极的临近驱动临近测量系统，(fi)i1，i2，in(n-1)/2 (8)若设(fj)(fj)-(f1)(j=2,3,4)(9) 则有Re(fj)Re(fj)-Re(f1)(10)Im(fj)Im(fj)-Im(f1) (11)同样对于(fj，k)=(fj，k)-(f1，k) (j=2,3,4) (12)有Re(fj，k)=Re(fj，k)-Re(f1，k)
11、(13)Im(fj，k)=Im(fj，k)-Im(f1，k) (14)对于反投影图像重构3 ，已知边界测量数据Boundaryb1，b2，bn(n-1)/2 ，若经过反投影重构计算的内部阻抗分布数据为 Distrd1，d2，dfem （下标 fem为采用有限元法对组织进行剖分的单元总数目） ，将反投影重构运算用映射 T表示，则可有 DistrTBoundary (15)因为反投影法在微小扰动时，线性较好，故可认为 DistrTBoundary (16)故此，可得到(fj，k)k1 T(fj) (17)即Re(fj，k)k1TRe(fj) (18)Im(fj，k) k1TIm(fj) (19)利
12、用(18)、(19)根据不同频率边界数据采用反投影图像重构得到Re(fi，k)和Im(fi，k)后，对于每个单元 k，由公式(6)，(7)确定圆心 Y的坐标，进而在复平面唯一确定三向量(OA), (OB), (OC)终点 A， B， C 所在的圆周方程 以散射系数为例，此时的可以近似地表示为 =2arctg(a/b)/ (20)其中，a, b为 ColeCole圆的圆心 Y对应的复数 a+jb在以 O为原点的复平面的数值根据所有单元的数值即可构成组织阻抗的散射系数断层图像. 3 结果 用在物理模型上采集的多频边界数据5进行了初步的验证实验成像目标是实验人员探入盐水槽中的手指.盐水槽直径 D=3
13、0 cm，高H=10 cm.进行重构运算时,边界数据采用(fj)/|(f1)|的形式. 重构图像如 Fig 4所示Fig 2标示了成像目标在物理模型中的方位, Fig 3为重构数据采用|(f2)|(f1)|/|(f1)|，以反投影法得到的准静态断层图像作为对照，Fig 4 为散射系数的断层图像 4 讨论 从 Fig 4看，有一定的图像噪声.这是参数图像重 Fig 4Parameter tomogram of the object(Dispersion parameter)构运算过程中，不同频率下的边界数据经过反投影、坐标变换、参数计算等多重运算过程，测量系统固有的噪声与数据处理系统的逐次计算误
14、差积累的结果同时，因为数据处理中用到了反投影算法，反投影的伪影数据参与了后续数据处理，使得目标边缘有所畸变对于成像中的噪声，可设法优化计算处理流程并研究适当的图像后处理方法加以改进对于重构图像边缘的伪影，可以设计适合 EIT的滤波反投影算法，以消减参与参数图像重构运算中的伪影数据的影响. 需要明确的是,在复平面内,由向量(OA), (OB), (OC)确定的圆的方程与(OA), (OB), (OC)所在的 ColeCole圆弧存在一个坐标平移的换算关系这种情况下，只有在准确的(f1,k) (即向量 OO)的模值和相角的基础上，R、R0、fc 的断层图像才有真实意义或参考价值对于在复平面中作为角
15、度量的散射系数 ，当(f1,k)在低频段时，其相角较小， 受到坐标变换甚至尺度伸缩的影响很小,重构图像可直接反映散射系数 的分布. 从成像整体结果来看，ColeCole 模型参数中的散射系数的断层图像在形态上与成像目标在物理模型中的位置和形状较吻合，定位准确，是一种可行和有效的参数成像方法相信经过改进，成像质量会有进一步改观. 【参考文献】 1 Cole KS, Cole RH. Dispersion and absorption in dielectrics J. J Chem Phys, 1941;9:341. 2 Saacson D, Isaacson EL. Comment on Ca
16、lderons paper: “On a increase boundary value problem”J. Math Comput, 1989;52:553. 3 Barber DC, Seagar AD. Fast reconstruction of resistance images J. Clin Phys Physiol Meas,1987;8 (Suppl.2A):47. 4 Tang M, Peng CL. Bioimpedance measurement theory and technology J. J Biomed Eng, 1997;14(2):152. 5史学涛. 用于电阻抗多频及参数成像的多频数据采集技术研究.D.西安:第四军医大学生物医学工程系，2001.
