一种实用的动态电阻抗断层成象算法.pdfVIP

维普资讯 3 一 ； 第 19卷 第 3期 中 国 生 物 医 学 工 程 学 报 Vo1．19 3 2o00年 9月 CHINESE JOURNALOFBIOMEDICALENGINEERING Septentzr200O 一 种实用的动态 电阻抗断层成象算法 董秀珍 秦明 汤孟兴 付 峰 史学涛 尤富 ～ 一 ～ 第四晕 ji强物医学工程系医电工程教研室，西安710032) ， 在电阻抗断层成象技术(EIT)的研究中提 出并推导 了一种较为严格的图象重构算法一广义逆 法 ．在此基础上针对 EIT的病态性问题 ，对该算法进行了改进 ，并在一种新 的计算机有限元模型上 实现 了该算法 ，重构 出了反映阻抗变化量的图象 (即动态成象)。改进的广义逆法突破了原广义逆 往在剖分规模 上的限制 ，能在较大 的剖分规模上进行 ．更有利于 图象 的重构 。将该算法用 于计 算 机仿真研究 ，结果表 明：改进的广义逆法与现有算法一等位线法相 比，成象定位更准确、重构 图象 的 误差小，在信号信噪 比较高时能较好地重构 出阻抗分布复杂 的模型 最后 ．对改进该算法 的噪声 性能进行 了讨 论 关键 词： 电阻抗 ；断层成象；算法 ；广义逆矩阵 ～ 分类 号： P,318．19 0 引 言 电阻抗断层成象 (EIT)技术是医学成像技术的一个新方 向。它可分为静态 (static)和动态 (dynamic)两种。动态 EIT技术以被测物体 (人体)电阻抗的变化量的相对值为成象 目标 ．通 过 向该物体注入驱动 电流 (电压)．测量边界电压 (电流 )，利用测量到的边界数据重构出反映物 体阻抗变化相对值的图象(目前的研究基本限于二维)…。而静态 EIT与此不 同之处在于其 要得到的是阻抗分布的绝对值 。动态 EIT技术是 目前在进行临床应用研究并取得了较好结 果的 EIT技术 J，但其成象有一定的误差 ，造成误差 的原因之一是现有的算法和重构模型都 有不足之处。这里提出了一种较为严格的动态 图象重构算法一广义逆法 ，并对算法进行 了改 进 ．在一种新的有限元模型上实现 了该算法 ，用此算法在仿真模型上重构出了较为理想的图 象 。 1 广义逆重构算法的推导 当给被测物体 (人体 )施加驱动 电流后 ．便在物体 内形成一个 电场 ，场 电位 函数 由 Laplace方程给 出： = 一 (1／c) · c · 口 其 中 c是场域 电导率分布(人体的阻抗这里可认为伫似是实的 ，是纯 电阻[2])。重建阻抗 断层 图像就是 巳知 求 c。而动态成像则是 已知 的变化量求 c的变化量 。由于 是 c的函 数 ，这是一个非线性 问题 ，直接求解 f很困难 。 Laplace方程可写成 ： 口： Rr · (1) Rr=一lnc (c为区域 电导率 ) 设 是方程 (1)的解 ，而当 Rr变化至Rr+Rp时，方程变为： 口： ( Rr+VRp)-