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

news2024/9/25 7:25:03

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)。

在这里插入图片描述

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

γ = σ + i ∗ ω ∗ ϵ γ 复电导率 ( a d m i t t i v i t y ) σ 电导率 ( c o n d u c t i v i t y ) ω 交流频率 ( A C f r e q u e n c y ,实验用的 125 k i l o h e r t z ) ϵ 介电常数 ( p e r m i t t i v i t y ) u 电势 ( e l e c t r i c p o t e n t i a l , 或者说内部的电压 ) \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公式的推导

d i v g e r e n c e ( γ ) ∗ g r a d ( 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 的电阻抗断层扫描工具包,用于设计和构建健康和运动传感设备。

硬件

芯片用的 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的直接重建方法应用于功能性肺成像

  • 电阻抗成像百度百科
  • https://en.wikipedia.org/wiki/Electrical_impedance_tomography
  • What is the complex conductivity of a metal?

电阻抗成像技术(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

  • Deep D-Bar:使用深度神经网络进行实时电阻抗层析成像
    在这里插入图片描述
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好像还是一个留学项目相关的东西

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/176002.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

SLF4J基本使用

文章目录1. 日志门面2. SLF4J的使用2.1 SLF4J入门2.2 为什么要使用SLF4J作为日志门面?2.3 绑定日志的实现(Binding)2.4 logback 整合 slf4j2.5 log4j 整合 slf4j2.6 jul 整合 slf4j3. 桥接旧的日志框架(Bridging)3.1 S…

完全二叉树的权值

题目描述 给定一棵包含 N 个节点的完全二叉树,树上每个节点都有一个权值,按从 上到下、从左到右的顺序依次是 A1​,A2​,⋅⋅⋅AN​,如下图所示: 现在小明要把相同深度的节点的权值加在一起,他想知道哪个深度的节点 权值之和最大?如果有多个深度的权值和同为最大,请你输…

react用高阶组件优化文件结构 帮助建立高阶组件应用思路

其实高阶组件是一个将组件写的更灵活的方式,他的应用场景在业务开发中会非常多样 这里 我们演示一种 主要还是解决问题的思想最重要 或者是 这个不叫解决问题 而是设计组件结构的思路 我们来模拟一个场景 在src下有一个 components 文件夹目录 在 components 下有…

3.开发社区首页

数据库设计如下: 任务:开发社区首页,显示前十个帖子;开发分页组件,分页显示所有帖子。 开发流程:entity->dao->mapper.xml->service->controller 由于关注后端,所以未使用前后分离…

【Potplayer】如何用Potplayer实现对视频的局部缩放播放?如何用potplayer更舒心地播放竖屏视频?

一、问题背景 有时候我们从互联网上下载得到一个竖屏视频,用电脑播放时,左右两边都是黑的,电脑屏幕的高宽度没有得到良好应用。而如果恰好这个竖屏视频其实只有一小部分才是有效区域/重点内容,比如中央区域,上下区域都…

win10怎么玩经典扫雷?繁星软件园推荐大家试试Windows7Games吧,亲测好用的扫雷Win10电脑版下载!!!

扫雷Win10电脑版(Windows7Games)是专门为微软最新版本的Windows 10/11操作系统打造的一款非常不错的消除类PC休闲小游戏安装包, Windows7Games游戏安装包为中文版,包含经典扫雷、空当接龙、红心大战、国际象棋等微软低版本操作系统Windows自带的多款超级…

我的硕士前半生

本篇文章属于随笔类,它可能无法对你起到直接帮助,它只是我这个普通学生一年半以来的一些足迹与思考。本文首发于我的个人博客 Forever Young 我的本科像硕士,有实验室有工位,桌子超级大,有老师带有学长可以问。 我的硕…

分享122个ASP源码,总有一款适合您

ASP源码 分享122个ASP源码,总有一款适合您 下面是文件的名字,我放了一些图片,文章里不是所有的图主要是放不下..., 122个ASP源码下载链接:https://pan.baidu.com/s/1-hwBdIxuv1NvpapRJsGAcQ?pwdgs1l 提取码&#x…

设计模式_行为型模式 -《迭代器模式》

设计模式_行为型模式 -《迭代器模式》 笔记整理自 黑马程序员Java设计模式详解, 23种Java设计模式(图解框架源码分析实战) 概述 定义 提供一个对象来顺序访问聚合对象中的一系列数据,而不暴露聚合对象的内部表示。 结构 迭代器…

手撕Pytorch源码#3.Dataset类 part3

写在前面手撕Pytorch源码系列目的:通过手撕源码复习了解高级python语法熟悉对pytorch框架的掌握在每一类完成源码分析后,会与常规深度学习训练脚本进行对照本系列预计先手撕python层源码,再进一步手撕c源码版本信息python:3.6.13p…

Couplet | 用Python写一副对联送给诸位科研汪!~

1写在前面 小伙伴们大家新年好啊!!!🥰 又是一年新春到,玉兔祝福要记牢:蹦蹦跳跳身体棒,平平淡淡精神爽,红红火火财运旺,和和气气朋友广,简简单单幸福长。&…

Windows SDK编程 初学笔记

#include "windows.h"int WINAPI WinMain(HINSTANCE hinstance, HINSTANCE hPreInstance, PSTR szCmdLine, int iCmdShow) {MessageBox(NULL, TEXT("来见见世面"), TEXT("Say Hi"), MB_OK);return 0; } MessageBox第一个参数为句柄,第…

Socket通信

什么是Socket?

Java基础——运算符与表达式

目录 Eclipse下载 安装 使用 运算符 键盘录入 Eclipse下载 安装 使用 Eclipse的概述(磨刀不误砍柴工)——是一个IDE(集成开发环境)Eclipse的特点描述(1)免费 (2)纯Java语言编写 (3)免安装 &#xff08…

【手把手教你学51单片机】

注:本文章转载自《手把手教你学习51单片机》!因转载需要原文链接,故无法选择转载! 如若侵权,请联系我进行删除!上传至网络博客目的为了记录自己学习的过程的同时,同时能够帮助其他一同学习的小伙…

AJAX Axios 总结

AJAX & Axios1. AJAX1.1 作用①与服务器进行数据交换②异步交互异步和同步1.2 基本使用1.3 案例SelectUserServlet&#xff1a;register.html&#xff1a;register.html中的<script2. Axios异步框架2.1 基本使用2.2 案例axiosServlet&#xff1a;axios-demo.html&#x…

Elasticsearch7.8.0版本高级查询—— 聚合查询文档

目录一、初始化文档数据二、聚合查询文档2.1、概述2.2、对某个字段取最大值 max 示例2.3、对某个字段取最小值 min 示例2.4、对某个字段求和sum 示例2.5、对某个字段取平均值 avg 示例2.6、对某个字段的值进行去重之后再取总数 示例三、State 聚合查询文档3.1、概述3.2、示例一…

目标检测论文解读复现【NO.24】改进 YOLOv5s 的轨道障碍物检测模型轻量化研究

前言此前出了目标改进算法专栏&#xff0c;但是对于应用于什么场景&#xff0c;需要什么改进方法对应与自己的应用场景有效果&#xff0c;并且多少改进点能发什么水平的文章&#xff0c;为解决大家的困惑&#xff0c;此系列文章旨在给大家解读最新目标检测算法论文&#xff0c;…

Cadence PCB仿真使用Allegro PCB SI生成反射仿真报告及报告导读图文教程

🏡《Cadence 开发合集目录》   🏡《Cadence PCB 仿真宝典目录》 目录 1,概述2,生成报告3,报告导读4,总结1,概述 本文简单介绍使用Allegro PCB SI生成网络的反射性能评估的报告的方法,及反射报告要点导读。 2,生成报告 第1步,选择需要生成报告的网络,然后单击右…

(侯捷C++)1.2面向对象高级编程(上)

1.整体结构 2.三大函数&#xff1a;拷贝构造&#xff0c;拷贝赋值&#xff0c;析构 拷贝构造&#xff1a;第一次出现对象&#xff0c;使用拷贝构造进行创建&#xff0c;例如&#xff1a;String s3(s1)。拷贝赋值&#xff1a;对象已经构造&#xff0c;重新赋值&#xff0c;例如…