一种基于局部细剖分的节点电阻抗成像方法技术

技术编号:20913642 阅读:34 留言:0更新日期:2019-04-20 09:08
本发明专利技术公开了一种基于局部细剖分的节点电阻抗成像方法,方法包括步骤:构建具有扰动物的目标模式,扰动物和背景区域构成待求解的目标区域,采用指定幅值的电流输入目标区域进行激励,基于激励获取边界电压矩阵U;对目标区域以三角形为单元通过有限元法剖分,建立FEM模型;基于FEM模型采用一步牛顿法获取目标区域的第一阻抗信息,根据第一阻抗信息对目标区域做局部细剖分;根据FEM模型中的节点与单元的关系,求出经局部细剖分的目标区域的节点雅克比矩阵,并采用Newton‑Raphson法求解逆问题,获取经局部细剖分的目标区域节点的第二阻抗信息;将第二阻抗信息转换成单元阻抗,基于单元阻抗进行图像重构成像;本发明专利技术能在减少计算量的同时提升成像的精度。

A Node Electrical Impedance Imaging Method Based on Local Subdivision

The invention discloses a node impedance imaging method based on local fine partition, which includes steps: constructing a target mode with perturbation, constituting a target area to be solved by perturbation and background area, exciting the target area by current input with specified amplitude, obtaining boundary voltage matrix U based on excitation, and adopting finite element method for the target area with triangle as unit. Based on FEM model, the first impedance information of the target area is obtained by one-step Newton method, and the target area is subdivided locally according to the first impedance information. According to the relationship between nodes and elements in FEM model, the node Jacobian matrix of the target area subdivided locally is obtained, and the inverse problem is solved by Newton Raphson method to obtain the local subdivision. The second impedance information of the node in the target area is transformed into the unit impedance, and the image reconstruction image is performed based on the unit impedance. The present invention can reduce the calculation amount and improve the imaging accuracy.

【技术实现步骤摘要】
一种基于局部细剖分的节点电阻抗成像方法
本专利技术属于电阻抗成像
,具体涉及一种基于局部细剖分的节点电阻抗成像方法。
技术介绍
传统的医学影像是借助于某种介质与人体的相互作用,把人体内部组织器官的形态结构、密度、功能等以影像的方式表现出来。而EIT则是根据人体内不同器官组织在不同的生理、病理状态下具有不同的电阻抗特性,采用各种激励模式,给人体施加安全驱动电流/电压,通过驱动电流或电压在人体的测量响应信息,重建人体内部的电阻抗分布或其变化的图像。在EIT技术中,根据成像目标的不同,可以分为动态成像和静态成像两类。静态成像对被测区域的实际阻抗值进行成像,能反映阻抗分布的实际情况,与动态成像相比,所包含的信息较大,对理论和临床诊断具有重要的意义,因此静态成像算法成为当前EIT技术研究的重点。然而,实际测量时存在模型误差和测量噪声,成像效果迅速变坏,产生含有噪声的图像,使其成像效果严重失真。动态成像利用两个时刻的电压差值,经过相应的算法重建两时刻阻抗差值的分布情况。其优点是通过相减运算可以消除数据中的干扰和和噪声信号,对系统测量数据的精度要求不会太高,所以可用于对某些疾病的发展过程进行实时图像监测。然而,在临床应用中前一时刻的数据有时获取不到,或者两个时刻被测物体的电阻抗值一样则无法重构出阻抗差值的分布情况。EIT成像是一个非线性逆问题,在成像理论和具体实现等方面都存在一些问题有待解决,这些问题在一定程度上限制了EIT技术的应用和发展。如何实现系统的高精度、高分辨率和算法的快速收敛是目前EIT技术的主要问题。对电阻抗重构算法的研究显示:修正的牛顿-拉夫逊法具有重构误差小、收敛性好和对初值要求不高的优点,其缺点在于单步迭代时间长,存储空间要求很大,从而限制了其进入实际阻抗图像重构;等位线法是目前最快的一种阻抗重构算法,但由于其在理论推倒上的不严格性所得结果具有较大误差。
技术实现思路
本专利技术目的是针对上述现有技术中存在的问题,提供一种基于局部细剖分的节点电阻抗成像方法,该方法能够对扰动目标进行有效成像,具体技术方案如下:一种基于局部细剖分的节点电阻抗成像方法,所述方法包括步骤:S1、构建具有扰动物的放置目标模式,采用指定幅值的电流输入目标区域进行激励,基于所述激励获取边界电压矩阵U;S2、对求解的目标区域以三角形为剖分单元通过有限元法剖分,建立FEM模型;S3、基于所述FEM模型采用一步牛顿法求解获取放置目标在所述激励下的第一阻抗信息,根据所述第一阻抗信息确定所述目标区域的阻抗分布图像,并对目标区域做局部细剖分;S4、根据FEM模型中的节点与所述剖分单元的关系,求出经局部细剖分的所述目标区域的节点雅克比矩阵,并采用Newton-Raphson法求解逆问题,获取经局部细剖分的所述目标区域节点的第二阻抗信息;S5、根据FEM模型中节点与所述剖分单元的关系,计算每个所述剖分单元的阻抗,根据所述阻抗将所述第二阻抗信息转换成单元阻抗,基于所述单元阻抗进行图像重构,得到阻抗成像图。进一步的,所述基于所述激励获取边界电压矩阵U包括:S11、在所述目标区域周围设定指定数量的电极,选取指定所述电极,依次通过所述电极的两个相邻电极注入指定幅值的电流,分别获取所述电极与所述两个相邻电极之间的电压;S12、依次选取所述目标区域周围所有的所述电极,执行步骤S11,获取构成所述便捷电压矩阵U的电压值。进一步的,所述步骤S2包括:基于所述目标区域,采用有限元法剖分将所述区域分成指定数量的同心圆,并设定同心圆中所有圆的半径,基于所述半径计算对应的每一所述同心圆的上的节点数目,采用makemeshdata.m函数编号所有所述节点,并按照设定规则连接所有所述节点构成所述FEM模型。进一步的,所述步骤S3包括:基于所述第一阻抗信息确定所述扰动物在所述FEM模型中对应的剖分单元,采用连接所述剖分单元各边中点的方法对所述目标区域做局部细剖分。进一步的,所述步骤S4包括:S41、基于局部细剖分后的所述成像区域确定单元雅克比矩阵,并由所述单元雅克比矩阵确定所述节点雅克比矩阵;S42、Newton-Raphson法中相应的正则化矩阵R转换成与与局部细剖分后的FEM模型相对应的节点正则化矩阵RN。进一步的,所述采用Newton-Raphson法求解目标区域阻抗分布的过程如下:a.设定放置目标的均匀阻抗为ρ0=1进行初始化,并基于所述均匀阻抗计算放置目标相应的电极电位分布U0;b.设定迭代条件||U-U0||<ε,若满足||U-U0||<ε,则终止迭代,并进行图像重构;否则转至步骤c;c.基于所述阻抗信息计算单元雅克比矩阵J(ρ);d.根据所述单元雅克比矩阵计算所述节点雅克比矩阵JN(ρ)和步长h;e.根据ρk+1=ρk+h,计算节点阻抗ρk+1;再将节点电导率转换为对应的单元电导率计算相应的电极电位分布Uk+1,转步骤b,b中的U0要更改为Uk+1,进行条件判断;进一步的,所述步骤S5中,所述FEM模型中每个单元的单元阻抗由节点阻抗转换构成;其中,所述单元阻抗由三个节点阻抗构成,所述单元阻抗的阻抗值为三个所述节点阻抗的阻抗值的平均值。本专利技术的基于局部细剖分的节点电阻抗成像方法,通过激励电流来获取边界电压矩阵,并通过有限元法剖分建立FEM模型,对FEM模型采用一步牛顿法获得相对应的第一阻抗信息,并基于第一阻抗信息对成像区域做细剖分处理,同时,采用Newton-Raphson法进行逆求解,得到局部细剖分后成像区域的第二阻抗信息,基于第二阻抗信息转换成FEM模型中的单元阻抗,根据单元阻抗进行图像重构实现成像操作;与现有技术相比,本专利技术的有益效果为:采用局部细剖分,减少计算量的同时提高了成像的精度;EIT采用Newton-Raphson类算法,Jacobian和Hessian矩阵的计算量很大,再加上各种矩阵的乘法、求逆等计算,计算量会更大,十分耗时,通过节点雅克比矩阵进行计算,在保证分辨率的前提下缩短运算时间,提高实时性;采用部细剖分的节点电阻抗重构成像方法,获得图像精度较高,图像定位准确。附图说明图1为本专利技术实施例中所述基于局部细剖分的节点电阻抗成像方法的流程图示意;图2为本专利技术实施例中电阻抗成像系统的电压测量方式;图3为本专利技术实施例中所述对成像区域的细剖分图示意。具体实施方式为了使本
的人员更好地理解本专利技术方案,下面将结合本专利技术实施例中的附图,对本专利技术实施例中的技术方案进行清楚、完整地描述。在本专利技术实施例中,提供了一种基于局部细剖分的节点电阻抗成像方法,具体参阅图1,所述方法包括步骤:S1、构建具有扰动物的放置目标模式,采用指定幅值的电流输入目标区域进行激励,基于激励获取边界电压矩阵U;优选的,本专利技术将水槽作为测试对象,并通过放置异物在水槽中作为扰动物;具体可参阅图2,边界电压矩阵U的获取采用相邻测量方式:通过相邻电极注入电流,轮流测量其他相邻电极对之间的电压;然后切换到下一个相临电极对其进行激励,再测量其他相邻,非激励电极对上的电压;最后,重复上述过程直到所有的相临电极对轮换一周,优选的,本实施例采用16电极模式,每一个激励下获取13个电压值,获取包含13*16个电压值的边界电压矩阵U;在获取边界电压过程中采用幅值为1mA的电流进行激励进行实验本文档来自技高网
...

【技术保护点】
1.一种基于局部细剖分的节点电阻抗成像方法,其特征在于,所述方法包括步骤:S1、构建具有扰动物的放置目标模式,采用指定幅值的电流输入目标区域进行激励,基于所述激励获取边界电压矩阵U;S2、对求解的目标区域以三角形为剖分单元通过有限元法剖分,建立FEM模型;S3、基于所述FEM模型采用一步牛顿法求解获取放置目标在所述激励下的第一阻抗信息,根据所述第一阻抗信息确定所述目标区域的阻抗分布图像,并对目标区域做局部细剖分;S4、根据FEM模型中的节点与所述剖分单元的关系,求出经局部细剖分的所述目标区域的节点雅克比矩阵,并采用Newton‑Raphson法求解逆问题,获取经局部细剖分的所述目标区域节点的第二阻抗信息;S5、根据FEM模型中节点与所述剖分单元的关系,计算每个所述剖分单元的阻抗,根据所述阻抗将所述第二阻抗信息转换成单元阻抗,基于所述单元阻抗进行图像重构,得到阻抗成像图。

【技术特征摘要】
1.一种基于局部细剖分的节点电阻抗成像方法,其特征在于,所述方法包括步骤:S1、构建具有扰动物的放置目标模式,采用指定幅值的电流输入目标区域进行激励,基于所述激励获取边界电压矩阵U;S2、对求解的目标区域以三角形为剖分单元通过有限元法剖分,建立FEM模型;S3、基于所述FEM模型采用一步牛顿法求解获取放置目标在所述激励下的第一阻抗信息,根据所述第一阻抗信息确定所述目标区域的阻抗分布图像,并对目标区域做局部细剖分;S4、根据FEM模型中的节点与所述剖分单元的关系,求出经局部细剖分的所述目标区域的节点雅克比矩阵,并采用Newton-Raphson法求解逆问题,获取经局部细剖分的所述目标区域节点的第二阻抗信息;S5、根据FEM模型中节点与所述剖分单元的关系,计算每个所述剖分单元的阻抗,根据所述阻抗将所述第二阻抗信息转换成单元阻抗,基于所述单元阻抗进行图像重构,得到阻抗成像图。2.如权利要求1所述的基于局部细剖分的节点电阻抗成像方法,其特征在于,所述基于所述激励获取边界电压矩阵U包括:S11、在所述目标区域周围设定指定数量的电极,选取指定所述电极,依次通过所述电极的两个相邻电极注入指定幅值的电流,分别获取所述电极与所述两个相邻电极之间的电压;S12、依次选取所述目标区域周围所有的所述电极,执行步骤S11,获取构成所述便捷电压矩阵U的电压值。3.如权利要求1所述的基于局部细剖分的节点电阻抗成像方法,其特征在于,所述步骤S2包括:基于所述目标区域,采用有限元法剖分将所述区域分成指定数量的同心圆,并设定同心圆中所有圆的半径,基于所述半径计算对应的每一所述同心圆的上的节点数目,采用makemeshdata.m函数编号所有所述节点,并按照...

【专利技术属性】
技术研发人员:戎舟成民民程大元李若愚
申请(专利权)人:南京邮电大学
类型:发明
国别省市:江苏,32

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1