Electrical Impedance Tomography,EIT

电阻抗成像(Electrical Impedance Tomography,EIT)是一种无创的以人体内部的电阻率分布为目标的重建体内组织图像的技术。人体是一个大的生物电导体,各组织、器官均有一定的阻抗,当人体的局部器官发生病变时,局部部位的阻抗必然与其他部位不同,因而可以通过阻抗的测量来对人体器官的病变进行诊断。

历史

1978年 1983 年 商用
John G. Webster发明 EIT 可视化人类前臂的横截面(David C. Barber and Brian H. Brown.) https://web.archive.org/web/20101202035314/http://www.maltronint.com/imaging.htm

2DEIT和3D EIT

从电阻抗成像的维度可以将胸腔阻抗成像的维度分为2DEIT和3D EIT。2D EIT即用一层电极带重建出电极层位置的阻抗分布图像,其优点是数据量少,电极相对好固定,鲁棒性好,目前有较为成熟的2DEIT模型和重建算法,Drager公司发布的产品也是2D EIT成像;其缺点是能提供的信息有限,在2DEIT中考虑的是电极层位置的阻抗变化,但是实际的测量过程中,在电极层以外的阻抗变化也会对电流的分布产生影响,所以在理论上2DEIT并不能重建出真实的阻抗分布情况。3D EIT即通过多个电极层实现对成像目标的三维重建,其优点是可以提供更多的信息,目前国内有河北工业大学的李颖研究小组,天津大学的王稼祥研究小组等研究小组在进行一些3DEIT肺部成像的尝试;但是3D EIT的缺点也很明显,它的数据量大,系统复杂,鲁棒性差,受模型影响较大,电极固定困难等。以目前的研究成果而言,3DEIT系统在肺部成像中是否真的有意义或者有必要研究还需要进一步的评估,因为3D EIT 系统在多个方面需要更多的代价。第一是硬件成本,3DEIT系统需要更多的采集通道,相应的硬件成本会增加,系统的设计复杂度会大大增加。第二,3DEIT系统需要付出更多的时间成本,对于一个16电极的系统的每次采集数据量是256,对于一个32个电极的系统的每次数据量是1024,EIT每扫描一帧图像的数据量会随着电极数呈平方倍增长,所以在系统的数据采集上就会消耗大量的时间成本;而图像的重建时间也会大大的增加,基于此,如果用普通的串行采集系统,很难实现3DEIT的实时性,如果采用全并行的设计,系统的成本会大幅度增加,并且通道的一致性也会受到影响。第三,系统的可靠性会受到影响,3DEIT意味着需要安装更多的电极,需要更复杂的系统,引入更多的通道一致性噪声。其中电极系统的不稳定,通道的一致性噪声都会直接影响重建图像的质量,如果不能保证电极系统的可靠和稳定,3DEIT系统会失去其价值。第四,3D EIT系统虽然能提供比2D EIT更多的信息,但是这些信息是否对医学诊断有实际的指导意义还需要进一步的评估。
————————————————
版权声明:本文为CSDN博主「zyk1120102464」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/zyk1120102464/article/details/78143944

原理

  • 科罗拉多州立大学ACE 1 EIT系统

电阻抗成像技术实际上就是:通过人体特定部位注入已知电压来测量在体表所引起的电流,或者注入一已知电流来测量在体表所引起的电压,利用所测量的电流电压值,依照一定的重建算法,计算出人体内部各组织、器官在电场作用下所呈现的阻抗分布,利用计算机产生断层成像。

通过粘性电极、电极带或导电电极背心围绕感兴趣的身体部位定位表面电极后,将在两个或多个驱动电极上施加通常为几毫安、频率为 10–100 kHz 的交流电。剩余的电极将用于测量产生的电压。然后将针对许多“刺激模式”重复该过程,例如连续的相邻电极对,直到完成整个圆圈,并且图像重建可以由结合了复杂数学算法和先验数据的数字工作站进行和显示。算法有线性反投影法(linear back projection,LBP) 法,奇异值分解,以Newton-Raphsom 算法、Tikhonov迭代算法和神经网络深度学习等。

符号定义

admittivity 复电导率形成原因假若物质处于随时间而变的电场,则电导率为复数(或对于异向性物质,电导率为复矩阵),称为导纳率(admittivity)。

复电导率将电流密度的大小和相位与施加电场的大小和相位相关联,类似于交流电路中的电流和电压与复电导相关。当您记得只使用了电流和电场的复值的实部并且您用复指数表示复电导率时,最容易看出这一点σ=∣σ∣exp(−iϕ)σ= |σ| exp( − i ϕ )σ=∣σ∣exp(−iϕ)。 |σ|是比例常数的大小和φ是电流场相位( current-field phase shift)。

γ=σ+i∗ω∗ϵγ复电导率(admittivity)σ电导率(conductivity)ω交流频率(ACfrequency,实验用的125kilohertz)ϵ介电常数(permittivity)u电势(electricpotential,或者说内部的电压)\gamma = \sigma+ i* \omega* \epsilon\\ \gamma 复电导率(admittivity)\\ \sigma 电导率(conductivity)\\ \omega 交流频率(AC frequency,实验用的125 kilohertz)\\ \epsilon 介电常数(permittivity)\\ u 电势(electric potential,或者说内部的电压)γ=σ+i∗ω∗ϵγ复电导率(admittivity)σ电导率(conductivity)ω交流频率(ACfrequency,实验用的125kilohertz)ϵ介电常数(permittivity)u电势(electricpotential,或者说内部的电压)

公式推导视频

生物化学 电阻抗成像OpenEIT 番外篇 EIT公式的推导

divgerence(γ)∗grad(u)=0∈Ωu=f∈∂Ωdivgerence(\gamma) * grad(u) = 0 \in \Omega \\ u=f \in \partial \Omega divgerence(γ)∗grad(u)=0∈Ωu=f∈∂Ω









MATLAB实现

文献: Jennifer L Mueller 和 Samuli Siltanen(SIAM 2012)所著的《线性和非线性逆问题及其实际应用 》一书的第 12-15 章 。下面我们参考[MS2012]这本书。
         D-bar 方法绕过了这些问题:它是一种基于从 EIT 数据计算电导率的非线性傅里叶变换然后反转变换的直接方法。无需解决任何直接问题。正则化由非线性频域中的低通滤波提供。由于低通滤波器,重建是模糊的,但这可以通过各向异性扩散、频域填充或深度学习在后处理中得到帮助。

更多有关EIT的论文

  • 以下不需要access的多篇数学描述pdf

  • http://www.siltanen-research.net/publ/Dbar_demystified.pdf

  • http://www.siltanen-research.net/publ/KnudsenLassasMuellerSiltanen2007.pdf

  • https://www.sce.carleton.ca/faculty/adler/publications/2019/hamilton-2019-compare-d-bar.pdf

  • https://www.researchgate.net/publication/259367835_A_Data-Driven_Edge-Preserving_D-bar_Method_for_Electrical_Impedance_Tomography

  • 以下的需要access

  • 通过 Schur 补码特性将基于样本的先验引入 D-Bar 方法(IEEE access)

  • Effect of Domain Shape Modeling and Measurement Errors on the 2-D D-Bar Method for EIT(IEEE access)

  • 用于电阻抗断层扫描的 D-bar 方法——揭秘(access)

  • Dbar方法在多分量非线性方程中的应用(知网access)

  • 用于EIT图像重建的深度神经网络模型构建方法(知网access)

EIT-kit

EIT-kit:MIT 设计的 ESP32 电阻抗断层扫描工具包

麻省理工学院 (MIT) 计算机科学与人工智能实验室 (CSAIL) 的一组工程师在上个月的ACM 用户界面软件与技术 (UIST) 研讨会上发表了一篇关于人机交互 (HCI) 的论文。在他们的论文中,Junyi Zhu、Jackson C. Snowden、Joshua Verdejo、Emily Chen、Paul Zhang、Hamid Ghaednia、Joseph H. Schwab 和 Stefanie Mueller 介绍了 EIT-kit,这是一种基于 ESP32 的电阻抗断层扫描工具包,用于设计和构建健康和运动传感设备。

硬件

嘉立创EDA支持 Eagle 文件导入(但是仅支持 Eagle v6.0 及以上)。

芯片用的 ESP32

openeit

https://openeit.github.io/docs/html/index.html

https://www.crowdsupply.com/mindseye-biomedical/spectra

pyeit : A python based framework for Electrical Impedance Tomography

pyeit paper

https://github.com/eitcom/pyEIT

该项目中包含三种算法 - Graz Consensus、Gauss-Newton 或 Jacobian 方法以及 Back Projection。每个都有自己的优点和缺点,可以调整。Spectra 套件使用 EIDORS 的 python 端口,称为 pyEIT。这意味着您可以在 python 中使用大部分 EIDORS 功能。

硬件

  • 芯片用的ADUCM350(精密微控制器 | 亚德诺(ADI)半导体)

  • 他们使用的标准硬件好像是这个众筹项目。,正在众筹的是Altium版本的,eagle是现在开源了的,eagle学习推荐《高质量PCB设计入门》

电极材料

py-dbar

D-bar 方法可用于二维框架中的电阻抗层析成像。该方法基于 Nachman 的理论重建程序和 Siltanen 建立的正则化策略,并为数值实现提供了一个自然框架。

https://pypi.org/project/py-dbar/

实验数据

  • https://blog.fips.fi/tomography/eit/the-d-bar-method-for-electrical-impedance-tomography-experimental-data/

参考与更多

  • 生物电阻抗分析在疾病临床监测和诊断中的应用

  • 开源生物医学成像系统,为每个人提供安全、可访问的成像

  • 电阻抗成像硬件系统中的常用概念

Jennifer Mueller - Direct reconstruction methods in EIT applied to functional pulmonary imaging EIT的直接重建方法应用于功能性肺成像

电阻抗成像技术(EIT)与人工智能相结合未来发展如何?

链接:https://www.zhihu.com/question/409512600/answer/2089095978

电阻抗断层成像技术(Electrical impedance tomography, EIT), 是一种以重构电阻抗为目标的成像技术,目前发展的较好的是动态成像,主要用于肺功能监测,通气监测等。如下图[1]:动态成像,也称为时差成像,即重构的时不同时刻阻抗的变化值。对于肺部癌症诊断来说,我们很难获取癌症前的正常组织的电阻抗测量数据,在段时间内癌变组织的阻抗变化也较小,因此动态成像无法进行癌变组织的检测。对于癌症检测来说需要利用准静态或静态EIT,而该项技术截止目前还没有取得较好的研究成果。虽然成像很难,但与电阻抗相关的生物电阻抗分析(Bioelectrical impedance analysis,BIA)技术在癌症检测领域也取的了很好的进展[2]。[1]ZONG Z, WANG Y, WEI Z. A Review of Algorithms and Hardware Implementations in Electrical Impedance Tomography (Invited) [J]. Progress In Electromagnetics Research, 2020, 169(59-71).[2] BAYFORD R, BERTEMES-FILHO P, FRERICHS I. Topical issues in electrical impedance tomography and bioimpedance application research [J]. Physiol Meas, 2021, 41(12): 120301.

深度方法 Deep D-bar: Real time Electrical ImpedanceTomography Imaging with Deep Neural Networks

https://github.com/asHauptmann/DeepDbar
机器学习增强的二维材料电阻抗层析成像
[1]Hong S, Lee J, Bae J, et al. A 10.4 mWelectrical impedance tomography SoC for portable real-time lung ventilationmonitoring system[C]//Solid-State Circuits Conference (A-SSCC), 2014 IEEEAsian. IEEE, 2014: 193-196.[2] Hong S, Lee K, Ha U, et al. A 4.9 mΩ-Sensitivity MobileElectrical Impedance Tomography IC for Early Breast-Cancer Detection System[J].Solid-State Circuits, IEEE Journal of, 2015, 50:245-257.



一些演示视频

德尔格学院:全身麻醉中肺保护性通气及EIT的应用-马加海教授 1:30:14

醉德我心-胸阻抗断层成像技术(EIT)的临床应用-张军教授 1:11:40

EIT在围术期临床应用进展-潘晓蕊 40:13

生物电阻抗成像仪(EIT)成像演示 01:08

基于电阻抗成像技术的手势识别设备原型机00:35

PS: EIT好像还是一个留学项目相关的东西

生物化学 电阻抗成像OpenEIT,Dbar算法,数据集等(暂记)相关推荐

  1. 生物化学 电阻抗成像OpenEIT 番外篇 EIT公式

    EIT简介 摘要电阻抗断层扫描(EIT)是一种成像方式,使用无害的电流探测患者或物体.电流通过放置在靶表面上的电极馈送,数据由在电极处测量的电压组成,这些电压由一组线性独立的电流注入模式产生.EIT旨 ...

  2. 电阻抗断层成像开源软件eidors3.8的入门使用方法

    声明 1)本文仅供学术交流,非商用.具体引用的资料请看参考文献.如果某部分不小心侵犯了大家的利益,请联系博主删除. 2)本人才疏学浅,整理总结的时候难免出错,还望各位前辈不吝指正,谢谢. 联系方式:k ...

  3. comsol如何定义狄利克雷边界_COMSOL与Visual C++三维电阻抗有限元联合建模与仿真研究...

    COMSOL 与 Visual C++ 三维电阻抗有限元联合建模与仿 真研究 张 涛,章伟睿,徐灿华,代 萌,杨 滨,史学涛,董秀珍,付 峰 [摘 要] [ 摘要 ] 目的:充分利用 COMSOL 软 ...

  4. 《MATLAB智能算法30个案例》:第28章 支持向量机的分类——基于乳腺组织电阻抗特性的乳腺癌诊断

    <MATLAB智能算法30个案例>:第28章 支持向量机的分类--基于乳腺组织电阻抗特性的乳腺癌诊断 1. 前言 2. MATLAB 仿真示例 3. 小结 1. 前言 <MATLAB ...

  5. 基于deap脑电数据集的脑电情绪识别二分类算法(附代码)

    想尝试一下脑电情绪识别的各个二分类算法. 代码主要分为三部分:快速傅里叶变换处理(fft).数据预处理.以及各个模型处理. 采用的模型包括:决策树.SVM.KNN三个模型(模型采用的比较简单,可以直接 ...

  6. 机器学习算法 01 —— K-近邻算法(数据集划分、归一化、标准化)

    文章目录 系列文章 K-近邻算法 1 什么是K-近邻算法 1.1 K-近邻算法(KNN)概念 1.2 举例说明K-近邻算法 1.3 K-近邻算法流程总结 2 K-近邻算法API初步使用 2.1 Sci ...

  7. 前端电商 sku 的全排列算法

    前端电商 sku 的全排列算法 什么是sku 聊聊常见的需求 解决思路 思路分解 上代码 什么是sku 针对电商而言: 1.SKU是指一款商品,每款都有出现一个SKU,便于电商品牌识别商品. 2.一款 ...

  8. 数据结构与算法随学随记

    1.数据结构是一门研究非数值计算的程序设计问题中的操作对象,以及它们之间的关系和操作等相关问题的学科. 2.传统上,我们把数据结构分为逻辑结构和物理结构. 逻辑结构:是指数据对象中数据元素之间的相互关 ...

  9. 《安富莱嵌入式周报》第297期:开源生物医学成像系统,可肺部成像,C算法合集500例,突出极致运算速度,数值方法书籍,芯片级激光隔离,3D打印机固件Marlin

    往期周报汇总地址:嵌入式周报 - uCOS & uCGUI & emWin & embOS & TouchGFX & ThreadX - 硬汉嵌入式论坛 - P ...

最新文章

  1. 关于Sqlserver中 BCP命令,实现数据导出功能。
  2. mysql什么格式转换_MySQL日期格式转换
  3. 微型计算机1tb等于,计算机应用基础第一套试卷100分
  4. 小米手机60帧录屏_手机录屏怎样只录手机内部声音不录入外部声音?教你三种方法,一定能帮到你...
  5. Ubuntu20.04下面运行applet与freemind部署到web上(调研+找到替换方案)
  6. cx oracle dpi 1002,cx_Oracle连接数据库总结
  7. kangle虚拟主机空间销售网站源码
  8. 【SpringCloud】Spring cloud Alibaba Sentinel 服务熔断
  9. 黑客利用未修复的 Atlassian 服务器攻击电信运营商和 ISP
  10. 数据结构与算法(Python)第三天
  11. 建议能在园子里面发布Silverlight2.0应用
  12. Java和jsp编程中应该注意的几个常见问题
  13. 3dmax用vr渲染还是cr渲染器?
  14. 实验九 哈希表的查找操作
  15. Siemens NX 1930中文版
  16. 实时系统性能测试指标及方法
  17. mac10.15切换大小写指示灯不亮
  18. vi 和 vim 编辑器的使用
  19. 苹果手机的计算机怎么设置快捷键大全,教程方法;苹果电脑快捷键大全最常用的都在这里了电脑技巧-琪琪词资源网...
  20. PlantCV中文文档

热门文章

  1. Python基本图形绘制
  2. golang自动构建脚本
  3. 衬线字体和无衬线字体
  4. 原来将Excel表格转换成应用程序如此简单
  5. Redis学习之incr命令
  6. 联想高管加盟乐视,手机市场风云突变
  7. 如何备份整个硬盘为镜像文件?
  8. python基本数据类型之字典详解
  9. 知网搜论文没有服务器响应,CNKI知网论文检测使用中的十四个常见问题
  10. 调研字节码插桩技术,用于互联网分布式系统监控设计和实现!