CN102008303A - 一种电阻抗成像系统频差fnoser成像方法 - Google Patents
一种电阻抗成像系统频差fnoser成像方法 Download PDFInfo
- Publication number
- CN102008303A CN102008303A CN 201010525155 CN201010525155A CN102008303A CN 102008303 A CN102008303 A CN 102008303A CN 201010525155 CN201010525155 CN 201010525155 CN 201010525155 A CN201010525155 A CN 201010525155A CN 102008303 A CN102008303 A CN 102008303A
- Authority
- CN
- China
- Prior art keywords
- voltage
- resistivity
- electrical impedance
- distribution
- fnoser
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明涉及一种基于高精度信号测量的电阻抗成像系统的频差FNOSER成像方法，包括高精度信号测量的电阻抗成像系统和成像方法，电阻抗成像系统包括壳体、电极阵列、激励模块、测量模块、中央控制处理电路、USB通讯模块、上位机；成像方法特征是：(1)采用电阻抗成像测量系统采集激励电流频率为fr时的边界电压，并将其作为参考电压Vp；(2)采用电阻抗成像测量系统采集激励电流频率为fm时的边界电压，并将其作为测量电压Vu；(3)计算被测模型的雅克比矩阵J；(4)将测量值及计算值代入FNOSER算法计算电阻率分布；(5)绘制电阻率分布图。本发明可以从算法角度上抑制测量系统的系统误差，提出一种确定初始电阻率的实用可行的方法，并可快速、准确的重构被测模型的电阻率分布，对电阻抗动态成像效果起到明显的改进作用。
Description
技术领域
本发明针对电阻抗成像系统（Electrical Impedance Tomography, EIT）中逆问题重构过程中，系统误差不能消除、初始电阻率不好选择的瓶颈问题，发明了一种基于频差FNOSER成像算法。
背景技术
电阻抗成像技术（Electrical Impedance Tomography, EIT）是一种特点鲜明的成像技术，主要应用于临床医学诊断。EIT技术通过配置于人体体表的电极阵列对人体施加电激励并测量其电响应信号，从中提取与人体生理、病理状态相关的组织与器官的电特性信息，不但可以反映解剖学结构，更能给出功能性图像结果。EIT技术不使用核素或射线，对人体无创、无电离辐射损伤，可以多次测量、重复使用。其设备成本低廉，不要求特殊的工作环境，便于大规模推广。
电阻抗成像方法可分为静态和差分成像，后者有时也称为动态成像。在差分成像中，为获得一幅图像需要两组独立的电压测量值，而重构图像则显示了被测对象内部电阻抗分布的差异或是变化，这种变化的来源可能是时间、频率或其它因素。差分成像通常是求解线性化的重构问题（逆问题），采用诸如滤波反投影法(Filtered Back-Projection Method, FBP)、快速一步牛顿法(Fast Newton’s One-Step Error Reconstructor, FNOSER)等方法，成像速度能得到保证。更为重要之处在于，差分成像在很大程度上能消除很多未知模型参数的影响，换言之，差分成像对系统误差不敏感，其原因在于两次电压测量在基本相同的条件下进行。这一点对实际的临床应用特别重要。
反观静态成像(Static Imaging)，虽然能够重构出绝对电阻抗分布，理论上更有价值，但由于众所周知EIT逆问题的病态性，少许的重构输入误差则会引起很大的成像误差甚至导致成像失败。所以静态重构本质上比差分重构更为困难。
如前所述，差分重构需要两组输入数据——参考电压数据集和测量电压数据集，前者通常是采用均匀场域下的测量数据，这在仿真和物理模型实验中容易得到，但实际临床应用中，几乎无法找到电阻率均匀分布的情况。因此，差分成像虽然具有对系统误差不敏感的优点，但其应用局限性较大。所以在临床中，只好退而求其次，在人体电阻抗分布存在变化的条件下应用动态EIT成像得到差分图像。
发明内容
本发明的目的就是提出一种基于高精度信号测量的电阻抗成像系统的频差FNOSER成像方法。针对FNOSER算法，参考电压数据集V p 和测量电压数据集V u 分别来自两个不同频率下的测量结果，将这两个频率命名为参考频率f r 和测量频率f m ，则构成了FNOSER频差重构算法。FNOSER频差重构能得到清晰稳定的阻抗分布图像，而且对误差有较好的容受度,解决了系统误差不能消除、初始电阻率不好选择的瓶颈问题。
本发明提出一种基于高精度信号测量的电阻抗成像系统的频差FNOSER成像方法，其特征在于：
1、构建高精度信号测量的电阻抗成像系统
所述电阻抗成像系统包括壳体、电极阵列、激励模块、测量模块、中央控制处理电路、USB通讯模块、上位机；其中测量电极和测量电路安装在同一探头盒中,大大缩短了测量电路与测量电极之间的距离,让测量电极的引出线能就近接入测试系统,从而形成一种新颖的用于电阻抗成像的高精度测量系统。
2、采集激励电流频率为f r 时的边界电压
采用具有USB通讯接口的电阻抗成像测量系统采集激励电流频率为f r 时的边界电压，将其作为参考电压V p ，边界电压包括每一个电极注入时其它电极测量得到的电压。
3、采集激励电流频率为f m 时的边界电压
采用具有上述电阻抗成像测量系统采集激励电流频率为f m 时的边界电压，将其作为测量电压V u ，边界电压包括每一个电极注入时其它电极测量得到的电压。
4、计算被测模型的雅克比矩阵J
1）、被测场域模型的建立及剖分：由于人体的个体差异，电阻抗成像系统不可能建立一个与实际被测场域完全一致的模型，只能通过人体躯干参数建立一个统一的模型；剖分采用美国麻省理工大学的开源剖分工具DistMesh软件包，得到被测模型的节点数N和单元数M；
2）、被测模型离散成M个单元，电阻率 
具有M个元素，每个元素为各个单元上的电阻率；为某一设定电阻率
分布下，各次电流注入下，使用有限元法计算各次电流注入下相邻电极上的电压；雅克比矩阵J由下式确定：
5、应用频差FNOSER算法计算电阻率分布
在某一设定电阻率
分布下，依次施加电流激励，测量的相邻电极上的电压为
，电阻抗成像就是要根据该测量数据
来计算场域内电阻率的分布。作为计算机仿真，可以用有限元法求解正问题来获得该数据。
如果计算得到的电阻率为
，则在该电阻率分布下，依次从相邻电极上注入电流，并利用有限方法计算各次电流注入下相邻电极上的电压
，如果电阻率为即为设定电阻率
，则应有：
由于各种误差的存在，上式的相等不可能实现，只可能使
与
近似相等。因此一般不去寻找真正的电阻率分布
以满足等式（2），而是寻找某一电阻率分布，使得该电阻率分布下计算的电压
与实际测量的电压
之间的差最小，即最小，通常用两组电压各元素差值的平方和来评估
与
之间的差别，即构造误差函数
上式中
是电阻率
的函数，电阻抗成像逆问题就是寻找电阻率分布
使得
最小，这在数学上是一个最小二乘问题。
将场域离散成M个单元。则电阻率向量
具有M个元素，每个元素为各个单元上的电阻率。为求解使
最小的值，可采用通常的求极值方法。
令
对
的各元素的偏导数为0，即：
用
来表示上式右端，即：
故式（4）可简化表达为：
上式共有M个非线性方程，构成一非线性方程组，以矩阵形式表示为：
（7）
其中：
（8）
以NR法求解方程组（7），其迭代公式为：
由于牛顿法进行电阻抗成像重构时，初始电阻率的选取至关重要，选择不当有可能使计算过程无法收敛，而且需要反复迭代来修正电阻率的分布，直到满足精度要求为止，计算量通常非常大，重构一幅图像一般需数小时以上，所以并不能用在实际应用中。
如果牛顿法中设定的初始电阻率足够接近于真实电阻率分布，则牛顿法中的迭代次数会大大减小，甚至只需要在设定的初始电阻率的基础上迭代一次，就可以得到场域电阻率分布的近似解，即使该近似解并不精确，它在数值上与真实解之间仍存在较大误差，但如果从重建的电阻率分布图像上仍可以看出场域内目标的位置、大小，这样重建算法的速度、收敛性均得到了解决，这就是FNOSER算法的基本思想。
FNOSER算法最关键的问题是如何设定初始电阻率
。只要
选择适当，就可以根据式（10）只迭代一步就完成图像重建，即
将步骤1、2得到的边界电压以及步骤3计算的雅克比矩阵代入步骤4的公式（5）就可以计算
，然后应用步骤1、2得到的边界电压、步骤3计算的雅克比矩阵以及公式（5）算出的
就可以根据式（11）计算最终电阻率分布。
虽然FNOSER的思想可以解决重建算法速度及收敛性等问题，但是在实际应用中使设定的初始电阻率接近真实电阻率非常困难。为解决这个问题，在激励电流频率分别为参考频率f r 和测量频率f m 时，采集参考电压数据集V p 和测量电压数据集V u ，并用参考电压V p 代替公式（2）～（5）中的
，测量电压V u 代替公式（2）～（5）中的
。这样相当于用参考频率f r 的激励电流确定初始电阻率，并且其与真实电阻率很接近，然后用参考频率f r 的激励电流测量，再按照FNOSER算法求解电阻率分布。
6、绘制电阻率分布图
根据步骤4得到的最终电导率分布
，绘制电导率分布图。
本发明能够应用于：
1. 电阻抗成像系统临床应用；
2. 电阻抗成像系统物理实验；
3. 电阻抗成像系统软件仿真。
本发明可达到电阻抗成像逆问题重构时间小于2s性能指标。
因此，本发明带来的积极效果是显而易见的，即应用本发明，能够在短时间之内计算被测场域的电阻率分布，并且绘出分布图。
不仅提高电阻抗成像的精度、速度，而且消除了系统误差对成像的影响，特别是提出一种确定初始电阻率的实用可行的方法。
基于上述因素，通过频差FNOSER算法成像，可实现电阻抗实时动态成像，节约时间，对电阻抗成像系统临床应用起到推进作用。
附图说明
图1是本发明电阻抗成像系统的结构图；
图2是本发明的算法流程图；
图3是本发明的实施例的实验水槽；
图4是本发明的实施例的计算场域图；
图5是本发明的实施例的电阻率分布图。
具体实施方式
下面结合附图给出的一个非限定性的实施例，对本发明作进一步的说明。
参见附图1，图中的高精度信号测量的电阻抗成像系统，包括壳体、电极阵列1、激励模块2、测量模块3、中央控制处理电路4、USB通讯模块5、上位机6；其特征是：
上位机6通过USB通讯模块5发送测量指令到中央控制处理电路4，中央控制处理电路4控制激励模块2输出激励电流，电流通过电极阵列1注入到被测对象，测量模块3采集电极阵列上的电压，并由中央控制处理电路4经USB通讯模块5送到上位机。
本实施例采用高精度信号测量的电阻抗成像系统测量实验水槽点边界电压数据，计算场域设定为16cm×8cm的长方形，如图3电极阵列位于上表面正中，以电极阵列的正中为坐标原点，电极数为8个，以1cm的间距排列成线阵，单个电极尺寸为5mm×15mm。图3中采用直径2cm，高2cm的香蕉块作为成像目标。
1、采集激励电流频率为f r 时的边界电压
本实施例中采用20kHz的激励电流测量参考电压V p ，为抑制接触阻抗的影响，在测量的时候不测激励电极上的电位差。所以为一V p 5乘8矩阵：
2、采集激励电流频率为f m 时的边界电压
本实施例中采用1MHz的激励电流测量参考电压V m ，为抑制接触阻抗的影响，在测量的时候不测激励电极上的电位差。所以为一V m 5乘8矩阵：
3、计算被测模型的雅克比矩阵J
计算场域设定为16cm×8cm的长方形，成像目标为5mm×5mm的方形区域。用DistMesh软件包剖分，单元数为3200，节点数为4653。并应用步骤2中公式（1）计算雅克比矩阵J，由于雅克比矩阵的物理意义是每一次电流激励时，每一个单元电导率变化在每一个边界电压上的反映，所以雅克比矩阵是一个3200×5×8的三维矩阵，为了计算存储方便，在这里转化为3200×40的二维矩阵：
4、应用频差FNOSER算法计算电阻率分布
由前三个步骤已得到的边界电压以及雅克比矩阵，现代入步骤4的公式（5）及（11）就可以最终电阻率分布。电导率分布为：
5、绘制电阻率分布图
根据步骤4得到的最终电导率分布，绘制电导率分布图，如说明书附图4。
Claims (1)
1.一种电阻抗成像系统频差FNOSER成像方法，其特征在于：
1）、构建高精度信号测量的电阻抗成像系统
所述电阻抗成像系统包括壳体、电极阵列、激励模块、测量模块、中央控制处理电路、USB通讯模块、上位机；其中测量电极和测量电路安装在同一个探头盒中,让测量电极的引出线能就近接入测试系统；
2）、采集激励电流频率为f r 时的边界电压
采用前述电阻抗成像测量系统采集激励电流频率为f r 时的边界电压，将其作为参考电压V p ，边界电压包括每一个电极注入时其它电极测量得到的电压；
3）、采集激励电流频率为f m 时的边界电压
采用上述电阻抗成像测量系统采集激励电流频率为f m 时的边界电压，将其作为测量电压V u ，边界电压包括每一个电极注入时其它电极测量得到的电压；
4）、计算被测模型的雅克比矩阵J
a）、被测场域模型的建立及剖分：通过人体躯干参数建立一个统一的模型；剖分采用开源剖分工具DistMesh软件包，得到被测模型的节点数N和单元数M；
b）、被测模型离散成M个单元，电阻率 
具有M个元素，每个元素为各个单元上的电阻率；
为某一设定电阻率
分布下，各次电流注入下，使用有限元法计算各次电流注入下相邻电极上的电压；雅克比矩阵J由下式确定：
4）应用FNOSER算法计算电导率分布
在某一设定电阻率
分布下，依次施加电流激励，测量的相邻电极上的电压为
，电阻抗成像就是要根据该测量数据
来计算场域内电阻率的分布。作为计算机仿真，可以用有限元法求解正问题来获得该数据；
如果计算得到的电阻率为
，则在该电阻率分布下，依次从相邻电极上注入电流，并利用有限方法计算各次电流注入下相邻电极上的电压
，如果电阻率为
即为设定电阻率
，则应有：
由于各种误差的存在，上式的相等不可能实现，只可能使
与
近似相等；并通过寻找某一电阻率分布
，使得该电阻率分布下计算的电压与实际测量的电压
之间的差最小，即
最小，用两组电压各元素差值的平方和来评估
与
之间的差别，即构造误差函数
（3）
上式中
是电阻率的函数，电阻抗成像逆问题就是寻找电阻率分布
使得
最小；
求解使
最小的
值，采用通常的求极值方法。
令
对的各元素的偏导数为0，即：
用来表示上式右端，即：
故式（4）可简化表达为：
上式共有M个非线性方程，构成一非线性方程组，以矩阵形式表示为：
其中：
（9）
以NR法求解方程组（7），其迭代公式为：
FNOSER计算：根据式（10）只迭代一步就完成图像重建，即
（11）
将步骤1、2得到的边界电压以及步骤3计算的雅克比矩阵代入步骤4的公式（5）就可以计算，然后应用步骤1、2得到的边界电压、步骤3计算的雅克比矩阵以及公式（5）算出的
就可以根据式（11）计算最终电阻率分布。
在激励电流频率分别为参考频率f r 和测量频率f m 时，采集参考电压数据集V p 和测量电压数据集V u ，并用参考电压V p 代替公式（2）～（5）中的
，测量电压V u 代替公式（2）～（5）中的
；这样相当于用参考频率f r 的激励电流确定初始电阻率，并且其与真实电阻率很接近，然后用参考频率f r 的激励电流测量，再按照FNOSER算法求解电阻率分布；
5）、绘制电阻率分布图
根据步骤4得到的最终电导率分布，绘制电导率分布图。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010105251559A CN102008303B (zh) | 2010-10-29 | 2010-10-29 | 一种电阻抗成像系统频差fnoser成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010105251559A CN102008303B (zh) | 2010-10-29 | 2010-10-29 | 一种电阻抗成像系统频差fnoser成像方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102008303A true CN102008303A (zh) | 2011-04-13 |
| CN102008303B CN102008303B (zh) | 2012-05-23 |
Family
ID=43838795
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2010105251559A Expired - Fee Related CN102008303B (zh) | 2010-10-29 | 2010-10-29 | 一种电阻抗成像系统频差fnoser成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102008303B (zh) |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102599907A (zh) * | 2012-04-09 | 2012-07-25 | 成都晨德科技有限公司 | 基于网格位移模型的电阻抗断层成像方法 |
| CN102599908A (zh) * | 2012-04-09 | 2012-07-25 | 成都晨德科技有限公司 | 基于带平衡因子的网格位移模型的电阻抗断层成像方法 |
| CN104321011A (zh) * | 2012-05-21 | 2015-01-28 | 通用电气公司 | 用于断层摄影成像的方法和系统 |
| CN108175408A (zh) * | 2018-01-15 | 2018-06-19 | 西安交通大学 | 一种生物电阻抗信号产生方法及模拟装置 |
| CN108652623A (zh) * | 2018-04-18 | 2018-10-16 | 天津大学 | 电学谱表征的自适应宽频电阻抗层析成像方法 |
| CN109101742A (zh) * | 2018-08-28 | 2018-12-28 | 南京航空航天大学 | 基于电阻抗成像损伤监测的复合材料强度预测方法 |
| CN109662712A (zh) * | 2018-12-25 | 2019-04-23 | 中国人民解放军第四军医大学 | 一种基于阻抗频谱约束的动态多频电阻抗断层成像方法 |
| CN110186583A (zh) * | 2019-06-13 | 2019-08-30 | 南京航空航天大学 | 基于电阻抗成像的陶瓷基复合材料高温部件温度测量方法 |
| CN111035387A (zh) * | 2019-12-27 | 2020-04-21 | 上海交通大学 | 提升电阻率成像获取图片质量的方法 |
| CN111543993A (zh) * | 2020-04-29 | 2020-08-18 | 北京航空航天大学 | 一种动态细胞阻抗成像装置 |
| CN112051309A (zh) * | 2020-08-06 | 2020-12-08 | 南方电网科学研究院有限责任公司 | 一种基于电阻抗成像的直流电缆附件无损检测方法及装置 |
| CN112147061A (zh) * | 2020-09-25 | 2020-12-29 | 辽宁工程技术大学 | 一种接地网状态电极聚焦式电阻抗成像系统及方法 |
| CN112401865A (zh) * | 2020-11-11 | 2021-02-26 | 中国科学技术大学 | 基于超形状的电阻抗成像方法 |
| CN113116333A (zh) * | 2021-04-13 | 2021-07-16 | 杭州永川科技有限公司 | 脑部阻抗分布差异数据的检测方法及系统 |
| CN113533439A (zh) * | 2021-06-25 | 2021-10-22 | 中国农业大学 | 玉米穗水分分布的电阻抗成像方法、系统及电子设备 |
| CN114224313A (zh) * | 2021-12-27 | 2022-03-25 | 深圳融昕医疗科技有限公司 | 电阻抗成像方法及计算机可读存储介质 |
| CN115245322A (zh) * | 2022-07-29 | 2022-10-28 | 南京航空航天大学 | 三维电极传感器、检测系统及方法、基于块稀疏贝叶斯学习数据处理方法 |
| CN115670422A (zh) * | 2023-01-04 | 2023-02-03 | 北京华睿博视医学影像技术有限公司 | 电阻抗成像系统 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101194834A (zh) * | 2006-12-05 | 2008-06-11 | 重庆博恩富克医疗设备有限公司 | 一种生物电阻抗测量方法及测量装置 |
-
2010
- 2010-10-29 CN CN2010105251559A patent/CN102008303B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101194834A (zh) * | 2006-12-05 | 2008-06-11 | 重庆博恩富克医疗设备有限公司 | 一种生物电阻抗测量方法及测量装置 |
Non-Patent Citations (2)
| Title |
|---|
| 《Automation Congress, 2008. WAC 2008. World》 20081002 Ciyong Luo,Minyou Chen,Ping Wang,Wei He A novel voltage measurement approach for image reconstruction of Electrical impedance tomography 1-4 1 , 2 * |
| 《计算机仿真》 20040930 黄嵩等 电阻抗成像中Jacobi矩阵的一种快速仿真算法 60-62,70 1 第21卷, 第9期 2 * |
Cited By (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102599907A (zh) * | 2012-04-09 | 2012-07-25 | 成都晨德科技有限公司 | 基于网格位移模型的电阻抗断层成像方法 |
| CN102599908A (zh) * | 2012-04-09 | 2012-07-25 | 成都晨德科技有限公司 | 基于带平衡因子的网格位移模型的电阻抗断层成像方法 |
| CN102599908B (zh) * | 2012-04-09 | 2014-07-23 | 成都晨德科技有限公司 | 基于带平衡因子的网格位移模型的电阻抗断层成像方法 |
| CN102599907B (zh) * | 2012-04-09 | 2015-03-04 | 思澜科技（成都）有限公司 | 基于网格位移模型的电阻抗断层成像方法 |
| CN104321011A (zh) * | 2012-05-21 | 2015-01-28 | 通用电气公司 | 用于断层摄影成像的方法和系统 |
| US9581627B2 (en) | 2012-05-21 | 2017-02-28 | General Electric Company | Method and system for tomographic imaging |
| CN108175408A (zh) * | 2018-01-15 | 2018-06-19 | 西安交通大学 | 一种生物电阻抗信号产生方法及模拟装置 |
| CN108652623A (zh) * | 2018-04-18 | 2018-10-16 | 天津大学 | 电学谱表征的自适应宽频电阻抗层析成像方法 |
| CN109101742A (zh) * | 2018-08-28 | 2018-12-28 | 南京航空航天大学 | 基于电阻抗成像损伤监测的复合材料强度预测方法 |
| CN109662712A (zh) * | 2018-12-25 | 2019-04-23 | 中国人民解放军第四军医大学 | 一种基于阻抗频谱约束的动态多频电阻抗断层成像方法 |
| CN109662712B (zh) * | 2018-12-25 | 2021-09-03 | 中国人民解放军第四军医大学 | 一种基于阻抗频谱约束的动态多频电阻抗断层成像方法 |
| WO2020248854A1 (zh) * | 2019-06-13 | 2020-12-17 | 南京航空航天大学 | 基于电阻抗成像的陶瓷基复合材料高温部件温度测量方法 |
| CN110186583A (zh) * | 2019-06-13 | 2019-08-30 | 南京航空航天大学 | 基于电阻抗成像的陶瓷基复合材料高温部件温度测量方法 |
| CN111035387A (zh) * | 2019-12-27 | 2020-04-21 | 上海交通大学 | 提升电阻率成像获取图片质量的方法 |
| CN111543993B (zh) * | 2020-04-29 | 2022-01-28 | 北京航空航天大学 | 一种动态细胞阻抗成像装置 |
| CN111543993A (zh) * | 2020-04-29 | 2020-08-18 | 北京航空航天大学 | 一种动态细胞阻抗成像装置 |
| CN112051309A (zh) * | 2020-08-06 | 2020-12-08 | 南方电网科学研究院有限责任公司 | 一种基于电阻抗成像的直流电缆附件无损检测方法及装置 |
| CN112051309B (zh) * | 2020-08-06 | 2022-07-22 | 南方电网科学研究院有限责任公司 | 一种基于电阻抗成像的直流电缆附件无损检测方法及装置 |
| CN112147061A (zh) * | 2020-09-25 | 2020-12-29 | 辽宁工程技术大学 | 一种接地网状态电极聚焦式电阻抗成像系统及方法 |
| CN112147061B (zh) * | 2020-09-25 | 2023-03-10 | 辽宁工程技术大学 | 一种接地网状态电极聚焦式电阻抗成像系统及方法 |
| CN112401865A (zh) * | 2020-11-11 | 2021-02-26 | 中国科学技术大学 | 基于超形状的电阻抗成像方法 |
| CN112401865B (zh) * | 2020-11-11 | 2024-05-17 | 中国科学技术大学 | 基于超形状的电阻抗成像方法 |
| CN113116333A (zh) * | 2021-04-13 | 2021-07-16 | 杭州永川科技有限公司 | 脑部阻抗分布差异数据的检测方法及系统 |
| CN113533439A (zh) * | 2021-06-25 | 2021-10-22 | 中国农业大学 | 玉米穗水分分布的电阻抗成像方法、系统及电子设备 |
| CN114224313A (zh) * | 2021-12-27 | 2022-03-25 | 深圳融昕医疗科技有限公司 | 电阻抗成像方法及计算机可读存储介质 |
| CN115245322A (zh) * | 2022-07-29 | 2022-10-28 | 南京航空航天大学 | 三维电极传感器、检测系统及方法、基于块稀疏贝叶斯学习数据处理方法 |
| CN115245322B (zh) * | 2022-07-29 | 2024-02-02 | 南京航空航天大学 | 基于块稀疏贝叶斯学习数据处理方法 |
| CN115670422A (zh) * | 2023-01-04 | 2023-02-03 | 北京华睿博视医学影像技术有限公司 | 电阻抗成像系统 |
| CN115670422B (zh) * | 2023-01-04 | 2023-03-17 | 北京华睿博视医学影像技术有限公司 | 电阻抗成像系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102008303B (zh) | 2012-05-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102008303B (zh) | 一种电阻抗成像系统频差fnoser成像方法 | |
| Lines et al. | Mathematical models and numerical methods for the forward problem in cardiac electrophysiology | |
| Woo et al. | A robust image reconstruction algorithm and its parallel implementation in electrical impedance tomography | |
| Niederer et al. | Verification of cardiac tissue electrophysiology simulators using an N-version benchmark | |
| Jehl et al. | A fast parallel solver for the forward problem in electrical impedance tomography | |
| Buist et al. | A deformable finite element derived finite difference method for cardiac activation problems | |
| CN106959325A (zh) | 一种平面阵列电极电容传感器、成像系统及成像方法 | |
| CN109646000B (zh) | 一种基于局部细剖分的节点电阻抗成像方法 | |
| CN108376124B (zh) | 一种用于电学成像的多导体体系导纳矩阵快速计算方法 | |
| CN102688041A (zh) | 基于电极采用交错位方式排列的三维电阻抗断层成像方法 | |
| CN109662712A (zh) | 一种基于阻抗频谱约束的动态多频电阻抗断层成像方法 | |
| Rocha et al. | A macro finite-element formulation for cardiac electrophysiology simulations using hybrid unstructured grids | |
| CN104851080A (zh) | 一种基于tv的三维pet图像重建方法 | |
| Zhang et al. | A numerical computation forward problem model of electrical impedance tomography based on generalized finite element method | |
| CN109758149B (zh) | 一种结合时差成像和频差成像的电阻抗成像方法 | |
| CN115808650A (zh) | 基于瞬时线化的电特性断层成像方法、系统、设备及介质 | |
| CN102915385B (zh) | 一种基于时域梯形法差分的互连线模型降阶方法 | |
| CN102599907A (zh) | 基于网格位移模型的电阻抗断层成像方法 | |
| Stenroos et al. | On the accuracy of collocation and Galerkin BEM in the EEG/MEG forward problem | |
| Zhang et al. | A meshfree method for simulating myocardial electrical activity | |
| CN110501587B (zh) | 一种基于电力线分布的Radon逆变换电容层析成像方法 | |
| RU2690107C2 (ru) | Устройство и способ формирования изображений распределения электрического импеданса с высоким разрешением | |
| Dumont et al. | The best of two worlds: The expedite boundary element method | |
| Buist et al. | Torso coupling techniques for the forward problem of electrocardiography | |
| Cho et al. | A fast EIT image reconstruction method for the two-phase flow visualization |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| ASS | Succession or assignment of patent right |
Owner name: STATE GRID SHANDONG ELECTRIC POWER CORPORATION ZIB Free format text: FORMER OWNER: CHONGQING UNIVERSITY Effective date: 20150803 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| C53 | Correction of patent for invention or patent application | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wang Tao Inventor after: Kong Gang Inventor after: Zhang Shidong Inventor after: Yang Lei Inventor after: Zhang Ning Inventor after: Yu Yang Inventor after: Gong Fangwei Inventor after: Yang Chao Inventor after: He Wei Inventor after: Li Bing Inventor after: He Chuanhong Inventor after: Liu Xinghua Inventor after: Zhang Li Inventor after: Luo Haijun Inventor after: Xu Zheng Inventor after: Sun Xuefeng Inventor after: Liu Hongzheng Inventor after: Guo Zhihong Inventor after: Chen Yufeng Inventor after: Luo Liyong Inventor after: Lv Juntao Inventor after: Xu Tianxi Inventor before: He Wei Inventor before: Li Bing Inventor before: He Chuanhong Inventor before: Zhang Li Inventor before: Luo Haijun Inventor before: Xu Zheng |
|
| COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: HE WEI LI BING HE CHUANHONG ZHANG LI LUO HAIJUN XU ZHENG TO: WANG TAO LIU XINGHUA SUN XUEFENG LIU HONGZHENG GUO ZHIHONG CHEN YUFENG LUO LIYONG LV JUNTAO XU TIANXI KONG GANG ZHANG SHIDONG YANG LEI ZHANG NING YU YANG GONG FANGWEI YANG CHAO HE WEI LI BING HE CHUANHONG ZHANG LI LUO HAIJUN XU ZHENG |
|
| TR01 | Transfer of patent right |
Effective date of registration: 20150803 Address after: 255000 Zhangdian District, Shandong, Beijing North Road, No. 67, No. Patentee after: Zibo Power Supply Company of State Grid Shandong Province Electric Power Company Address before: 400044 Shapingba District Sha Street, No. 174, Chongqing Patentee before: Chongqing University |
|
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120523 Termination date: 20151029 |
|
| EXPY | Termination of patent right or utility model |