首页 / 专利库 / 诊断设备和程序 / 医学影像学 / 层析成像 / 电阻抗成像 / 一种电阻抗成像系统频差FNOSER成像方法

一种电阻抗成像系统频差FNOSER成像方法

阅读:1004发布:2020-05-17

专利汇可以提供一种电阻抗成像系统频差FNOSER成像方法专利检索,专利查询,专利分析的服务。并且本 发明 涉及一种基于高 精度 信号 测量的 电阻 抗成像 系统的频差FNOSER成像方法,包括高精度信号测量的 电阻抗成像 系统和成像方法,电阻抗成像系统包括壳体、 电极 阵列、激励模 块 、测量模块、中央控制处理 电路 、USB通讯模块、上位机;成像方法特征是:(1)采用电阻抗成像测量系统采集激励 电流 频率 为fr时的边界 电压 ,并将其作为参考电压Vp;(2)采用电阻抗成像测量系统采集激励电流频率为fm时的边界电压,并将其作为测量电压Vu;(3)计算被测模型的雅克比矩阵J;(4)将测量值及计算值代入FNOSER 算法 计算 电阻率 分布;(5)绘制电阻率分布图。本发明可以从算法 角 度上抑制测量系统的系统误差,提出一种确定初始电阻率的实用可行的方法,并可快速、准确的重构被测模型的电阻率分布,对电阻抗动态成像效果起到明显的改进作用。,下面是一种电阻抗成像系统频差FNOSER成像方法专利的具体信息内容。

1.一种电阻抗成像系统频差FNOSER成像方法,其特征在于:
1)、构建高精度信号测量的电阻抗成像系统
所述电阻抗成像系统包括壳体、电极阵列、激励模、测量模块、中央控制处理电路、USB通讯模块、上位机;其中电极阵列和测量模块安装在同一个探头盒中,让电极阵列的引出线能就近接入测量模块;
2)、采集激励电流频率为参考频率fr时的边界电压
采用前述电阻抗成像系统采集激励电流频率为参考频率fr时的边界电压,将其作为参考电压Vp,边界电压包括每一个电极注入时其它电极测量得到的电压;
3)、采集激励电流频率为测量频率fm时的边界电压
采用上述电阻抗成像系统采集激励电流频率为测量频率fm时的边界电压,将其作为测量电压Vu,边界电压包括每一个电极注入时其它电极测量得到的电压;
4)、计算被测模型的雅克比矩阵J
a)、被测模型的建立及剖分:通过人体躯干参数建立一个统一的模型;剖分采用开源剖分工具DistMesh软件包,得到被测模型的节点数N和单元数M;
b )、 被 测 模 型 离 散 成 M 个 单 元 , 电 阻 率
具有M个元素,每个元素为各个单元上的电阻率
为电阻率 分布下,各次电流注入下,使用有限元法计算各次电流
注入下相邻电极上的电压;雅克比矩阵J由下式确定:
(1)
5)应用FNOSER算法计算电导率分布
假设已知某一模型的真实电阻率 分布下,依次施加电流激励,测量的相邻电极上的电压为 ,电阻抗成像就是要根据该测量的数据 来计算场域内电阻率的
分布,作为计算机仿真,能够用有限元法求解正问题来获得该数据;
针对同一模型,但电阻率分布未知,如果计算得到其电阻率为 ,则在该电阻率 分布下,依次从相邻电极上注入电流,并利用有限方法计算各次电流注入下相邻电极上的电压,如果计算所得模型电阻率为 即为已知电阻率 ,则应有:
(2)
由于各种误差的存在,上式的相等不可能实现,只可能使 与 近似相等;并通过寻找某一电阻率 分布,使得该电阻率 分布下计算的电压 与实际测量的电压之间的差最小,即 最小,用两组电压各元素差值的平方和来评
估 与 之间的差别,即构造误差函数
(3)
上式中 是电阻率 的函数,电阻抗成像逆问题就是寻找电阻率 分布使得最小;
求解使 最小的 值,采用通常的求极值方法:
令 对 的各元素的偏导数为0,即:
(4)
用 来表示上式右端,即:
(5)
故式(4)可简化表达为:
(6)
上式共有M个非线性方程,构成一非线性方程组,以矩阵形式表示为:
(7)
其中:
(8)
(9)
以NR法求解方程组(7),其迭代公式为:
(10)
FNOSER计算:根据式(10)只迭代一步就完成图像重建,即
(11)
将步骤2、3得到的边界电压以及步骤4计算的雅克比矩阵代入步骤5的公式(5)就可以计算 ,然后应用步骤2、3得到的边界电压、步骤4计算的雅克比矩阵以及公式(5)算出的 就可以根据式(11)计算电阻率分布;
在激励电流频率分别为参考频率fr和测量频率fm时,采集参考电压数据集Vp和测量电压数据集Vu,并用参考电压Vp代替公式(2)~(5)中的 ,测量电压Vu代替公式(2)~(5)中的 ;这样相当于用参考频率fr的激励电流确定初始电阻率,并且其与真实电阻率很接近,然后用参考频率fr的激励电流测量,再按照FNOSER算法求解电阻率分布;
6)、绘制电阻率分布图
根据步骤5得到的最终电阻率 分布,绘制电导率分布图。

说明书全文

一种电阻抗成像系统频差FNOSER成像方法

技术领域

[0001] 本发明针对电阻抗成像系统(Electrical Impedance Tomography, EIT)中逆问题重构过程中,系统误差不能消除、初始电阻率不好选择的瓶颈问题,发明了一种基于频差FNOSER成像算法

背景技术

[0002] 电阻抗成像技术(Electrical Impedance Tomography, EIT)是一种特点鲜明的成像技术,主要应用于临床医学诊断。EIT技术通过配置于人体体表的电极阵列对人体施加电激励并测量其电响应信号,从中提取与人体生理、病理状态相关的组织与器官的电特性信息,不但可以反映解剖学结构,更能给出功能性图像结果。EIT技术不使用核素或射线,对人体无创、无电离辐射损伤,可以多次测量、重复使用。其设备成本低廉,不要求特殊的工作环境,便于大规模推广。
[0003] 电阻抗成像方法可分为静态和差分成像,后者有时也称为动态成像。在差分成像中,为获得一幅图像需要两组独立的电压测量值,而重构图像则显示了被测对象内部电阻抗分布的差异或是变化,这种变化的来源可能是时间、频率或其它因素。差分成像通常是求解线性化的重构问题(逆问题),采用诸如滤波反投影法(Filtered Back-Projection Method, FBP)、快速一步顿法(Fast Newton’s One-Step Error Reconstructor, FNOSER)等方法,成像速度能得到保证。更为重要之处在于,差分成像在很大程度上能消除很多未知模型参数的影响,换言之,差分成像对系统误差不敏感,其原因在于两次电压测量在基本相同的条件下进行。这一点对实际的临床应用特别重要。
[0004] 反观静态成像(Static Imaging),虽然能够重构出绝对电阻抗分布,理论上更有价值,但由于众所周知EIT逆问题的病态性,少许的重构输入误差则会引起很大的成像误差甚至导致成像失败。所以静态重构本质上比差分重构更为困难。
[0005] 如前所述,差分重构需要两组输入数据——参考电压数据集和测量电压数据集,前者通常是采用均匀场域下的测量数据,这在仿真和物理模型实验中容易得到,但实际临床应用中,几乎无法找到电阻率均匀分布的情况。因此,差分成像虽然具有对系统误差不敏感的优点,但其应用局限性较大。所以在临床中,只好退而求其次,在人体电阻抗分布存在变化的条件下应用动态EIT成像得到差分图像。

发明内容

[0006] 本发明的目的就是提出一种基于高精度信号测量的电阻抗成像系统的频差FNOSER成像方法。针对FNOSER算法,参考电压数据集Vp和测量电压数据集Vu分别来自两个不同频率下的测量结果,将这两个频率命名为参考频率fr和测量频率fm,则构成了FNOSER频差重构算法。FNOSER频差重构能得到清晰稳定的阻抗分布图像,而且对误差有较好的容受度,解决了系统误差不能消除、初始电阻率不好选择的瓶颈问题。
[0007] 本发明提出一种基于高精度信号测量的电阻抗成像系统的频差FNOSER成像方法,其特征在于:
[0008] 1、构建高精度信号测量的电阻抗成像系统
[0009] 所述电阻抗成像系统包括壳体、电极阵列、激励模、测量模块、中央控制处理电路、USB通讯模块、上位机;其中测量电极和测量电路安装在同一探头盒中,大大缩短了测量电路与测量电极之间的距离,让测量电极的引出线能就近接入测试系统,从而形成一种新颖的用于电阻抗成像的高精度测量系统。
[0010] 2、采集激励电流频率为fr时的边界电压
[0011] 采用具有USB通讯接口的电阻抗成像测量系统采集激励电流频率为fr时的边界电压,将其作为参考电压Vp,边界电压包括每一个电极注入时其它电极测量得到的电压。
[0012] 3、采集激励电流频率为fm时的边界电压
[0013] 采用具有上述电阻抗成像测量系统采集激励电流频率为fm时的边界电压,将其作为测量电压Vu,边界电压包括每一个电极注入时其它电极测量得到的电压。
[0014] 4、计算被测模型的雅克比矩阵J
[0015] 1)、被测场域模型的建立及剖分:由于人体的个体差异,电阻抗成像系统不可能建立一个与实际被测场域完全一致的模型,只能通过人体躯干参数建立一个统一的模型;剖分采用美国麻省理工大学的开源剖分工具DistMesh软件包,得到被测模型的节点数N和单元数M;
[0016] 2)、 被 测 模 型 离 散 成 M 个 单 元, 电 阻 率具有M个元素,每个元素为各个单元上的电阻率;
为某一设定电阻率 分布下,各次电流注入下,使用有限元法计算
各次电流注入下相邻电极上的电压;雅克比矩阵J由下式确定:
[0017] (1)。
[0018] 5、应用频差FNOSER算法计算电阻率分布
[0019] 在某一设定电阻率 分布下,依次施加电流激励,测量的相邻电极上的电压为,电阻抗成像就是要根据该测量数据 来计算场域内电阻率的分布。作为计算机仿真,可以用有限元法求解正问题来获得该数据。
[0020] 如果计算得到的电阻率为 ,则在该电阻率分布下,依次从相邻电极上注入电流,并利用有限方法计算各次电流注入下相邻电极上的电压 ,如果电阻率为 即为设定电阻率 ,则应有:
[0021] (2)
[0022] 由于各种误差的存在,上式的相等不可能实现,只可能使 与 近似相等。因此一般不去寻找真正的电阻率分布 以满足等式(2),而是寻找某一电阻率分布 ,使得该电阻率分布下计算的电压 与实际测量的电压 之间的差最小,即最小,通常用两组电压各元素差值的平方和来评估 与之间的差别,即构造误差函数
[0023] (3)
[0024] 上式中 是电阻率 的函数,电阻抗成像逆问题就是寻找电阻率分布 使得最小,这在数学上是一个最小二乘问题。
[0025] 将场域离散成M个单元。则电阻率向量 具有M个元素,每个元素为各个单元上的电阻率。为求解使 最小的 值,可采用通常的求极值方法。
[0026] 令 对 的各元素的偏导数为0,即:
[0027] (4)
[0028] 用 来表示上式右端,即:
[0029] (5)
[0030] 故式(4)可简化表达为:
[0031] (6)
[0032] 上式共有M个非线性方程,构成一非线性方程组,以矩阵形式表示为:
[0033] (7)
[0034] 其中:
[0035] (8)
[0036] (9)
[0037] 以NR法求解方程组(7),其迭代公式为:
[0038] (10)
[0039] 由于牛顿法进行电阻抗成像重构时,初始电阻率的选取至关重要,选择不当有可能使计算过程无法收敛,而且需要反复迭代来修正电阻率的分布,直到满足精度要求为止,计算量通常非常大,重构一幅图像一般需数小时以上,所以并不能用在实际应用中。
[0040] 如果牛顿法中设定的初始电阻率足够接近于真实电阻率分布,则牛顿法中的迭代次数会大大减小,甚至只需要在设定的初始电阻率的基础上迭代一次,就可以得到场域电阻率分布的近似解,即使该近似解并不精确,它在数值上与真实解之间仍存在较大误差,但如果从重建的电阻率分布图像上仍可以看出场域内目标的位置、大小,这样重建算法的速度、收敛性均得到了解决,这就是FNOSER算法的基本思想。
[0041] FNOSER算法最关键的问题是如何设定初始电阻率 。只要 选择适当,就可以根据式(10)只迭代一步就完成图像重建,即
[0042] (11)
[0043] 将步骤1、2得到的边界电压以及步骤3计算的雅克比矩阵代入步骤4的公式(5)就可以计算 ,然后应用步骤1、2得到的边界电压、步骤3计算的雅克比矩阵以及公式(5)算出的 就可以根据式(11)计算最终电阻率分布。
[0044] 虽然FNOSER的思想可以解决重建算法速度及收敛性等问题,但是在实际应用中使设定的初始电阻率接近真实电阻率非常困难。为解决这个问题,在激励电流频率分别为参考频率fr和测量频率fm时,采集参考电压数据集Vp和测量电压数据集Vu,并用参考电压Vp代替公式(2)~(5)中的 ,测量电压Vu代替公式(2)~(5)中的 。这样相当于用参考频率fr的激励电流确定初始电阻率,并且其与真实电阻率很接近,然后用参考频率fr的激励电流测量,再按照FNOSER算法求解电阻率分布。
[0045] 6、绘制电阻率分布图
[0046] 根据步骤4得到的最终电导率分布 ,绘制电导率分布图。
[0047] 本发明能够应用于:
[0048] 1. 电阻抗成像系统临床应用;
[0049] 2. 电阻抗成像系统物理实验;
[0050] 3. 电阻抗成像系统软件仿真。
[0051] 本发明可达到电阻抗成像逆问题重构时间小于2s性能指标。
[0052] 因此,本发明带来的积极效果是显而易见的,即应用本发明,能够在短时间之内计算被测场域的电阻率分布,并且绘出分布图。
[0053] 不仅提高电阻抗成像的精度、速度,而且消除了系统误差对成像的影响,特别是提出一种确定初始电阻率的实用可行的方法。
[0054] 基于上述因素,通过频差FNOSER算法成像,可实现电阻抗实时动态成像,节约时间,对电阻抗成像系统临床应用起到推进作用。附图说明
[0055] 图1是本发明电阻抗成像系统的结构图;
[0056] 图2是本发明的算法流程图
[0057] 图3是本发明的实施例的实验槽;
[0058] 图4是本发明的实施例的计算场域图;
[0059] 图5是本发明的实施例的电阻率分布图。

具体实施方式

[0060] 下面结合附图给出的一个非限定性的实施例,对本发明作进一步的说明。
[0061] 参见附图1,图中的高精度信号测量的电阻抗成像系统,包括壳体、电极阵列1、激励模块2、测量模块3、中央控制处理电路4、USB通讯模块5、上位机6;其特征是:
[0062] 上位机6通过USB通讯模块5发送测量指令到中央控制处理电路4,中央控制处理电路4控制激励模块2输出激励电流,电流通过电极阵列1注入到被测对象,测量模块3采集电极阵列上的电压,并由中央控制处理电路4经USB通讯模块5送到上位机。
[0063] 本实施例采用高精度信号测量的电阻抗成像系统测量实验水槽点边界电压数据,计算场域设定为16cm×8cm的长方形,如图3电极阵列位于上表面正中,以电极阵列的正中为坐标原点,电极数为8个,以1cm的间距排列成线阵,单个电极尺寸为5mm×15mm。图3中采用直径2cm,高2cm的香蕉块作为成像目标。
[0064] 1、采集激励电流频率为fr时的边界电压
[0065] 本实施例中采用20kHz的激励电流测量参考电压Vp,为抑制接触阻抗的影响,在测量的时候不测激励电极上的电位差。所以为一Vp5乘8矩阵:
[0066] 。
[0067] 2、采集激励电流频率为fm时的边界电压
[0068] 本实施例中采用1MHz的激励电流测量参考电压Vm,为抑制接触阻抗的影响,在测量的时候不测激励电极上的电位差。所以为一Vm5乘8矩阵:
[0069] 。
[0070] 3、计算被测模型的雅克比矩阵J
[0071] 计算场域设定为16cm×8cm的长方形,成像目标为5mm×5mm的方形区域。用DistMesh软件包剖分,单元数为3200,节点数为4653。并应用步骤2中公式(1)计算雅克比矩阵J,由于雅克比矩阵的物理意义是每一次电流激励时,每一个单元电导率变化在每一个边界电压上的反映,所以雅克比矩阵是一个3200×5×8的三维矩阵,为了计算存储方便,在这里转化为3200×40的二维矩阵:
[0072] 。
[0073] 4、应用频差FNOSER算法计算电阻率分布
[0074] 由前三个步骤已得到的边界电压以及雅克比矩阵,现代入步骤4的公式(5)及(11)就可以最终电阻率分布。电导率分布为:
[0075] 。
[0076] 5、绘制电阻率分布图
[0077] 根据步骤4得到的最终电导率分布 ,绘制电导率分布图,如说明书附图4。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈