| 一种实用的动态电阻抗断层成象算法 在电阻抗断层成象技术(EIT)的研究中提出并推导了一种较为严格的图象重构算法-广义逆法，在此基础上针对EIT的病态性问题，对该算法进行了改进，并在一种新的计算机有限元模型上实现了该算法，重构出了反映阻抗变化量的图象(即动态成象)。改进的广义逆法突破了原广义逆法在剖分规模上的限制，能在较大的剖分规模上进行，更有利于图象的重构。将该算法用于计算机仿真研究，结果表明：改进的广义逆法与现有算法-等位线法相比，成象定位更准确、重构图象的误差小，在信号信噪比较高时能较好地重构出阻抗分布复杂的模型。最后，对改进该算法的噪声性能进行了讨论。 分类号： R318.19 A PRACTICAL RECONSTRUCTION ALGORITHM OF DYNAMIC ELECTRICAL IMPEDANCE TOMOGRAPHY Dong Xiuzhen, Qin Mingxin, Tang Mengxing, Fu Feng, Shi Xuetao, You Fusheng (Bio Medical Engineering Department of the Forth Military Medical University, Xian 710032) ABSTRACT A rigid electrical impedance tomography(EIT) reconstruction algorithm-general inversion algorithm was presented. To improve the ill-posed condition, we modified the algorithm and implemented it on a new FEM model, and get the image of impedance change(dynamic image). This modified general inversion algorithm(MGIA) could be implemented on a bigger dimension FEM model, so could do good to reconstruction. When implementing this algorithm on computer simulative condition, we found that the MGIA had a more accurate location ability and a smaller reconstruction error than the now-generally-used algorithm-equi-potential-back-projection algorithm. When having a relatively high SNR data, it could reconstruct images of complicated model. Finally, noise performance and improving method were discussed. Key words: Electrical impedance; Tomography; General inversion matrix 0 引 言 电阻抗断层成象(EIT)技术是医学成像技术的一个新方向。它可分为静态(static)和动态(dynamic)两种。动态EIT技术以被测物体(人体)电阻抗的变化量的相对值为成象目标，通过向该物体注入驱动电流(电压)，测量边界电压(电流)，利用测量到的边界数据重构出反映物体阻抗变化相对值的图象(目前的研究基本限于二维)［1］。而静态EIT与此不同之处在于其要得到的是阻抗分布的绝对值。动态EIT技术是目前在进行临床应用研究并取得了较好结果的EIT技术［2］，但其成象有一定的误差，造成误差的原因之一是现有的算法和重构模型都有不足之处。这里提出了一种较为严格的动态图象重构算法-广义逆法，并对算法进行了改进，在一种新的有限元模型上实现了该算法，用此算法在仿真模型上重构出了较为理想的图象。 1 广义逆重构算法的推导 当给被测物体(人体)施加驱动电流后，便在物体内形成一个电场，场电位函数v由Laplace方程给出： 其中c是场域电导率分布(人体的阻抗这里可认为伫似是实的，是纯电阻［2］)。重建阻抗断层图像就是已知v求c。而动态成像则是已知v的变化量求c的变化量。由于v是c的函数，这是一个非线性问题，直接求解c很困难。 Laplace方程可写成： (1) 设vr是方程(1)的解，而当Rr变化至Rr+Rp时，方程变为： (2) 方程(2)的解设为v=vr+vp,则： 这里引入一个假设：电导率变化量Rp较小，则vp与vr相比较起来也较小，方程可近似为［3］： (3) 这表明在区域电导率变化较小的情况下，电导率对数的变化与电压的变化呈一种近似的线性关系。因此可以用一个线性的变换来表示这个关系： Δv=F*Δr(4) Δv表示电压变化量，Δr表示电导率对数的变化，F称为正向矩阵。一般情况下，F不是方阵，可设F＋是F的穆尔-彭罗斯广义逆阵，F阵为mn的矩阵，其中n＜m，则根据穆尔-彭罗斯广义逆阵性质定理有： F＋*F=E (其中E为单位阵)(7) 由式(4)可得： F＋*Δv=F＋*F*Δr(8) 将(7)式代入得： F＋*Δv=Δr(9) 以上推导表明：只要得到正向矩阵F的广义逆阵F＋，即可求得阻抗变化的分布。 2 算法以及模型的改进 广义逆法在理论上是可行的，但在具体实现过程中，由于F矩阵的病态性，使其对误差非常敏感。而计算机舍入误差、数值计算误差、测量误差又不可避免，所以，在有限元剖分(算法在计算机上的实现需要建立重构模型并将其离散，现在广泛采用成熟的有限元方法建模并将区域离散进行求解)规模较大时，该算法的结果是发散的。针对这一问题，对算法进行了改进，改进方法如下： 在求广义逆阵的时候需要先进行F矩阵的奇异值分解： F=U∑VT U和V代表矩阵的特征向量，∑是一个矩阵，对角线元素代表矩阵的奇异值。这里奇异值表明数据对重构图象的影响情况，小的奇异值表示图象受噪声影响大［4］。一个矩阵的病态性由矩阵的最大和最小奇异值的比决定。随着最小奇异值趋于零，病态性将越来越严重。而小的奇异值除了与EIT问题本身有关外，还与各种误差有关。这里，将分解得到的奇异值进行了修正，取一定小的值α为门限，小于α的奇异值均置为α，然后再求逆。这个过程相当于用一个病态性较小的矩阵近似代替原有的病态性较大的矩阵，虽然损失了一定的有用信息，但换取了病态性的改善，使广义逆法能在大的剖分规模上进行。 现在EIT技术研究中常用的有限元模型剖分方案不同，但单元大小基本一致(如图1)［5］，若按这种模型计算其正向矩阵F，则会发现其在一定规模剖分时就会病态严重。主要原因在于阻抗成象技术存在一个特点：边界测量对被测区域中心部分(即远离电极的部分)的阻抗变化不敏感，使算得的正向矩阵的条件数较大。若根据此正向矩阵进行图象重构，则区域中心阻抗的重构对噪声和误差非常敏感，一点误差将导致阻抗的很大变化，而且随着剖分元素的增多而尤为突出。对于 12
| 
