CN103156605B - 基于声电效应的生物组织电流密度成像方法 - Google Patents
基于声电效应的生物组织电流密度成像方法 Download PDFInfo
- Publication number
- CN103156605B CN103156605B CN201310104468.0A CN201310104468A CN103156605B CN 103156605 B CN103156605 B CN 103156605B CN 201310104468 A CN201310104468 A CN 201310104468A CN 103156605 B CN103156605 B CN 103156605B
- Authority
- CN
- China
- Prior art keywords
- biological tissue
- current density
- signal
- acoustoelectric
- ultrasonic
- 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.)
- Expired - Fee Related
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种基于声电效应的生物组织电流密度成像方法,包括以下步骤:(1)采用超声激励源产生超声波并经超声驱动、超声聚焦于生物组织内部;(2)对聚焦于生物组织内部的超声波进行连续定位,从而完成对生物组织整体或局部的扫描;(3)信号采集:超声波聚焦定位于某一空间位置,则会产生一个与该空间位置电流密度对应的高频电信号,得到生物组织整体或局部的一组信号;(4)根据采集到对生物组织整体或局部扫描后得到信号,计算和重建生物组织内部的电流密度。本发明的电流密度成像方法是一种动态成像,能够直观地反映生物组织内部活动情况,能都得到高精度、高分辨率的电流密度成像,便于医学诊断。
Description
技术领域
本发明涉及生物组织成像领域,具体地讲是一种基于声电效应的生物组织电流密度成像方法。
背景技术
医学图像具有直观、形象和信息量大的特点,在临床诊断中具有重要地位。随着医学成像技术的进步,医学影像已经从解剖结构成像发展到功能成像,即医学图像不仅可以显示疾病解剖结构的形态变化,而且可以反映机体器官的功能及相关的生理、生化变化,能够检测到其生理活动的异常,从而为疾病诊断提供了重要的信息,推动了医学的进步与发展。以心脏为例,心律失常是危害人类健康的重大临床病症,在杂乱的颤动波中确定病灶位置是心律失常诊断、治疗和生理机制解释的关键。
电阻抗成像技术是根据物体内部不同组织以及组织在不同功能状态下具有不同电阻抗这一原理,通过布放于体表的一系列电极作为传感器,然后给物体注入一定的电流(对人体而言需注入小的安全电流),再测量体表电位来重建物体内部的阻抗分布图像。电阻抗成像是一种非入侵性的功能成像技术,但是由于其被理论和实验证实的固有的空间分辨率限制,目前尚未在心脏电兴奋测量中得到广泛的应用。
在心电测量中,一对电极称为一个导联,心脏电兴奋测量实质是从体表导联测量电压获取心脏内部电流密度分布的过程。心电正问题是已知信号源求解体表电位,而逆问题是根据测量到的心脏表面的电位信号来推算心脏中产生这些电活动的源的位置和大小,从而进行病理、生理和认知方面的功能分析。如同其他领域的逆问题一样,心电逆问题中已知的只是某一方程组的解,而对产生这个解函数的源函数分布及方程的算子并不了解,它不可避免地遇到理论上的困难:在数学上逆问题不具备唯一解;无论选择哪种数值计算方法求解心电产生源,均不能改变其解的病态特性,即当问题涉及的参数与数据即使存在很小的误差或干扰时,也会引起解的不确定性,甚至会出现大幅度的剧烈振荡;测量数据的不完备性也将导致心电源求解问题的困难。长期以来,人们利用各种可能的数学理论研究心电逆问题的数值求解方法,以期望活得高精度、稳定、合理的解。其研究大致可以分为两类方法,一是基于等效偶极子的参数定位法,只有在心电活动局限在很小的区域范围内时才是对真实情况的近似,实际活动区域可能向任意方向扩展,呈现不同的空间几何形状,这时等效偶极子的逆映射结果可能不理想;一是基于电流分布模型的图像重建法。
基于传统的心电测量的电流密度成像需要尽可能多的导联场来求解病态的逆问题,由于电极深入体表在组织损伤的代价下可以改善电流源位置的确定,但是这一逆问题受限于组织入侵程度和成像空间分辨率之间的折衷。对于心脏这样的生物组织,侵入组织来布置电极会造成组织器官的严重损害,难以在该领域推广应用。同时,采用布置导联对生物组织进行电激励,并采集产生的低频电信号,由于产生的电信号强度较弱,难以得到清晰的成像。
因此,基于以上这些技术问题,可见单一的电阻抗成像技术难以实现有效的高精度、高分辨率成像,造成EIT技术在临床应用的进展一直较为缓慢。
发明内容
本发明要解决的技术问题是,提供一种采用超声激励,并与生物组织自身的电信号协同作用,并能实现高精度、高分辨率成像的基于声电效应的生物组织电流密度成像方法。
本发明的技术解决方案是,提供以下步骤的基于声电效应的生物组织电流密度成像方法,包括以下各步骤:
(1)采用超声激励源产生超声波并经超声驱动、超声聚焦于生物组织内部;
(2)对聚焦于生物组织内部的超声波进行连续定位,从而完成对生物组织整体或局部的扫描;
(3)信号采集:超声波聚焦定位于某一空间位置,则会产生一个与该空间位置电流密度相对应的高频信号,由多路探测导联采集信号并经过放大和滤波后输出,在完成超声波聚焦扫描后,便能得到生物组织整体或局部的一组信号;
(4)将采集后的电压信号导入计算机,根据采集到对生物组织整体或局部扫描后得到信号,计算和重建生物组织内部的电流密度,从而实现生物组织的电流密度成像。
采用本发明的方法,与现有技术相比,本发明具有以下优点:本发明无需采用外部电激励,将超声波聚焦于生物组织上,使得超声波与生物组织内部电信号之间产生协同作用,在超声波的定位扫描下,能够得到相应点的高频电信号,能采集到大量均匀分布的信息,并根据这些点反映的数据信息计算和重建生物组织的电流密度分布,从而得到生物组织整体或局部的电流密度成像;由于生物组织不同时间的电活动是不同的,因此本发明的电流密度成像方法是一种动态成像,能够直观地反映生物组织内部活动情况,能都得到高精度、高分辨率的电流密度成像,便于医学诊断。
作为改进,步骤4中的电流密度计算和重建是指,生物组织内部电流在聚焦超声束作用下,生物组织会发生声电效应现象,在生物组织外部可以无创地测量到这个高频电信号的影响;根据声电效应理论,在生物组织内部的某作用点(x′,y′,z′),其探测导联声电效应信号V为:
其中Td为超声脉冲持续时长,K是声电相互作用的常数,ρ0为生物组织的导电阻抗,J为内部电流密度,P0是超声声压,L是电极导联,F是聚焦焦距,b(x,y,z)是超声束,a(t)是超声波波形;
通过逐点扫描的方式完成整个目标区域的测量,点的扫描构成时间序列;声电信号经过传导,可以经测量电极测量到:
M(x′,y,z′)=ρ0(x′,y′,z′)L(x′,y,z′)·JI(x′,y′,z′)
由于生物电电流密度是个低频信号,其分布具有稀疏性特征,可以通过一下优化得到电流密度分布:
min‖(M(x′,y′,z′)-ρ0(x′,y′,z′)L(x′,y′,z′)·JI(x′,y′,z′)‖2
作为改进,步骤3中的信号放大和滤波是通过以下方法实现的:多路探测导联通过一个截断频率为480kHz的模拟高通滤波器连接到一个增益为20dB、带宽为2~5MHz的差分放大器上,输出的信号再通过放大器再放大29dB,信号的采样频率为5MHz。该信号处理方法能够将声电信号(超声与电产生的协同效应的信号一般为高频信号)从采集到的信号中分离和放大,进一步提高了成像的清晰度。
作为改进,步骤1和2中,超声波的聚焦范围为短轴为1mm,长轴为3-4mm的椭圆形,以1mm为步长在生物组织表面进行生物组织内部聚焦定位扫描;该聚焦范围相比导联场达厘米级的积分范围小的多,超声束起到类似狄拉克(Dirac)函数的筛选作用。
附图说明
图1为本发明基于声电效应的生物组织电流密度成像方法的流程示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
如图1所示,本发明的基于声电效应的生物组织电流密度成像方法,包括以下步骤:
(1)采用超声激励源产生超声波并经超声驱动、超声聚焦于生物组织内部;
(2)对聚焦于生物组织内部的超声波进行连续定位,从而完成对生物组织整体或局部的扫描;
(3)信号采集:超声波聚焦定位于某一空间位置,则会产生一个与该空间位置电流密度相对应的高频信号,由多路探测导联采集信号并经过放大和滤波后输出,在完成超声波聚焦扫描后,便能得到生物组织整体或局部的一组信号;
(4)将采集后的电压信号导入计算机,根据采集到对生物组织整体或局部扫描后得到信号,计算和重建生物组织内部的电流密度,从而实现生物组织的电流密度成像。
步骤4中的电流密度计算和重建是指,生物组织内部电流在聚焦超声束作用下,生物组织会发生声电效应现象,在生物组织外部可以无创地测量到这个高频电信号的影响;根据声电效应理论,在生物组织内部的某作用点(x′,y′,z′),其探测导联声电效应信号V为:
其中Td为超声脉冲持续时长,K是声电相互作用的常数,ρ0为生物组织的导电阻抗,J为内部电流密度,P0是超声声压,L是电极导联,F是聚焦焦距,b(x,y,z)是超声束,a(t)是超声波波形;
通过逐点扫描的方式完成整个目标区域的测量,点的扫描构成时间序列;声电信号经过传导,可以经测量电极测量到:
M(x′,y′,z′)=ρ0(x′,y′,z′)L(x′,y′,z′)·JI(x′,y,z′)
由于生物电电流密度是个低频信号,其分布具有稀疏性特征,可以通过一下优化得到电流密度分布:
min‖(M(x′,y′,z′)-ρ0(x′,y′,z′)L(x′,y′,z′)·JI(x′,y′,z′)‖2。
步骤3中的信号放大和滤波是通过以下方法实现的:多路探测导联通过一个截断频率为480kHz的模拟高通滤波器连接到一个增益为20dB、带宽为2~5MHz的差分放大器上,输出的信号再通过放大器再放大29dB,信号的采样频率为5MHz。
步骤1和2中,超声波的聚焦范围为短轴为1mm,长轴为3-4mm的椭圆形,以1mm为步长在生物表面进行生物组织内部聚焦定位扫描;该聚焦范围相比导联场达厘米级的积分范围小的多,超声束起到类似狄拉克(Dirac)函数一样的筛选作用。
以上仅就本发明较佳的实施例作了说明,但不能理解为是对权利要求的限制。本发明不仅局限于以上实施例,其具体结构允许有变化。总之,凡在本发明独立权利要求的保护范围内所作的各种变化均在本发明的保护范围内。
Claims (3)
1.一种基于声电效应的生物组织电流密度成像方法,包括以下步骤:
(1)采用超声激励源产生超声波并经超声驱动、超声聚焦于生物组织内部;
(2)对聚焦于生物组织内部的超声波进行连续定位,从而完成对生物组织整体或局部的扫描;
(3)信号采集:超声波聚焦定位于某一空间位置,则会产生一个与该空间位置电流密度对应的高频电信号,由多路探测导联采集信号并经过放大和滤波后输出,在完成超声波聚焦扫描后,便能得到生物组织整体或局部的一组信号;
其特征在于:还包括如下步骤:
(4)根据采集到的生物组织整体或局部扫描后得到信号,计算和重建生物组织内部的电流密度,从而实现生物组织的电流密度成像;
步骤4中的电流密度计算和重建是指,生物组织内部电流在聚焦超声束作用下,生物组织会发生声电效应现象,在生物组织外部可以无创地测量到这个高频电信号的影响;根据声电效应理论,在生物组织内部的某作用点(x′,y′,z′),其探测导联声电效应信号V为:
其中K是声电相互作用的常数,ρ0为生物组织的导电阻抗,P0是超声声压,L是电极导联,F是聚焦焦距,b(x,y,z)是超声束,a(t)是超声波波形;
通过逐点扫描的方式完成整个目标区域的测量,点的扫描构成时间序列;声电信号经过传导,可以经测量电极测量到:
M(x′,y′,z′)=ρ0(x′,y′,z′)L(x′,y′,z′)·JI(x′,y′,z′)
由于生物电电流密度的分布具有稀疏性特征,可以通过以下优化来获得电流密度分布:
min||(M(x′,y′,z′)-ρ0(x′,y′,z′)L(x′,y′,z′)·JI(x′,y′,z′)||2。
2.根据权利要求1所述的基于声电效应的生物组织电流密度成像方法,其特征在于:步骤3中的信号放大和滤波是通过以下方法实现的:多路探测导联通过一个截断频率为480kHz的模拟高通滤波器连接到一个增益为20dB、带宽为2~5MHz的差分放大器上,输出的信号再通过放大器再放大29dB,信号的采样频率为5MHz。
3.根据权利要求1所述的基于声电效应的生物组织电流密度成像方法,其特征在于:步骤1和2中,超声波的聚焦范围为短轴为1mm,长轴为3-4mm的椭圆形,以1mm为步长在生物组织表面进行生物组织内部聚焦定位扫描。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310104468.0A CN103156605B (zh) | 2013-03-26 | 2013-03-26 | 基于声电效应的生物组织电流密度成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310104468.0A CN103156605B (zh) | 2013-03-26 | 2013-03-26 | 基于声电效应的生物组织电流密度成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103156605A CN103156605A (zh) | 2013-06-19 |
CN103156605B true CN103156605B (zh) | 2016-04-13 |
Family
ID=48580505
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310104468.0A Expired - Fee Related CN103156605B (zh) | 2013-03-26 | 2013-03-26 | 基于声电效应的生物组织电流密度成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103156605B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104644219B (zh) * | 2015-02-16 | 2017-02-22 | 天津大学 | 基于声电效应的功率密度成像方法 |
CN107536608A (zh) * | 2017-09-19 | 2018-01-05 | 天津大学 | 基于灵敏度场理论的非侵入式声电成像方法 |
CN108829238B (zh) * | 2018-05-04 | 2021-01-08 | 天津大学 | 一种基于超声调制的脑-机接口通讯方法 |
CN109645999B (zh) * | 2018-11-29 | 2022-04-26 | 天津大学 | 一种基于声电效应的4d经颅聚焦超声神经成像方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102894974A (zh) * | 2012-10-16 | 2013-01-30 | 中国科学院电工研究所 | 一种磁声电成像系统及成像方法 |
-
2013
- 2013-03-26 CN CN201310104468.0A patent/CN103156605B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102894974A (zh) * | 2012-10-16 | 2013-01-30 | 中国科学院电工研究所 | 一种磁声电成像系统及成像方法 |
Non-Patent Citations (1)
Title |
---|
measurement of a 2D electric dipole filed using the acousto-electric effect;Ragnar olafsson et.al.;《Medical Imaging 2007: Ultrasonic Imaging and Signal Processing》;20071231;第6531卷;65310S-1-65310S-10 * |
Also Published As
Publication number | Publication date |
---|---|
CN103156605A (zh) | 2013-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6345806B2 (ja) | 心臓組織のマッピングのための医療装置 | |
US8057390B2 (en) | High-resolution mapping of bio-electric fields | |
US20150065845A1 (en) | Measuring apparatus and its method | |
Olafsson et al. | Ultrasound current source density imaging | |
Qin et al. | Ultrasound current source density imaging of the cardiac activation wave using a clinical cardiac catheter | |
Deán-Ben et al. | Optoacoustic image formation approaches—a clinical perspective | |
CN103156605B (zh) | 基于声电效应的生物组织电流密度成像方法 | |
WO2019082892A1 (ja) | 超音波診断システム及び超音波診断方法 | |
CN110575627B (zh) | 一种快速获取目标神经治疗能量投送位点的物理标测装置 | |
JP2017529938A (ja) | 心臓組織のマッピング用医療装置 | |
Qin et al. | 4D acoustoelectric imaging of current sources in a human head phantom | |
Bera et al. | Electrical impedance tomography (EIT): a harmless medical imaging modality | |
CN105249933A (zh) | 光声分子三维成像仪 | |
CN103156604B (zh) | 超声协同的生物组织电阻抗成像方法 | |
JP6614910B2 (ja) | 光音響装置 | |
CN115708933A (zh) | 基于超声成像的聚焦超声治疗系统 | |
Kumar et al. | An experimental validation of bio-impedance technique for medical & non-medical application | |
Witte et al. | 1A-4 Acoustoelectric Detection of Current Flow in a Neural Recording Chamber | |
US20210059535A1 (en) | Handheld ultrasound transducer array for 3d transcranial and transthoracic ultrasound and acoustoelectric imaging and related modalities | |
Provost et al. | Ultrafast acoustoelectric imaging | |
KR20140094955A (ko) | 온도 영상 생성 방법 및 장치, 이를 포함한 초음파 시스템 | |
CN114209355B (zh) | 深部神经超声自动定位和标测方法、装置、设备和介质 | |
KR101956460B1 (ko) | 초음파 의료영상 장치를 이용한 미세석회화 조직 검출 방법 및 그 초음파 의료영상 장치 | |
Zheng et al. | A targeting method based on acoustic backscatter for treatment planning in tissue ablation using focused ultrasound | |
JP2019518544A (ja) | 電気音響画像誘導治療 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Xu Wenlong Inventor after: Xu Bingqiao Inventor after: Zhang Ning Inventor after: Song Xinghui Inventor before: Xu Wenlong Inventor before: Xu Bingqiao Inventor before: Zhang Ning |
|
COR | Change of bibliographic data | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160413 Termination date: 20170326 |
|
CF01 | Termination of patent right due to non-payment of annual fee |