CN106580249A - 一种注入电流式热声成像方法 - Google Patents
一种注入电流式热声成像方法 Download PDFInfo
- Publication number
- CN106580249A CN106580249A CN201610991573.4A CN201610991573A CN106580249A CN 106580249 A CN106580249 A CN 106580249A CN 201610991573 A CN201610991573 A CN 201610991573A CN 106580249 A CN106580249 A CN 106580249A
- Authority
- CN
- China
- Prior art keywords
- injection current
- formula
- objective body
- distribution
- imaging method
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
一种注入电流式热声成像方法,通过电极向成像目标体注入电流,在成像目标体中产生焦耳热,引起热膨胀,产生超声信号,利用超声探头检测超声信号。根据检测的超声信号,重建声源和目标体的电导率分布。本发明可用于高分辨率的电阻抗成像,用于疾病的早期诊断。
Description
技术领域
本发明涉及一种热声成像方法,特别涉及一种医学注入电流式热声成像方法。
背景技术
热声成像是目标体吸收能量,产生热膨胀而产生超声信号,利用超声进行对目标体结构进行重建。现有的热声成像根据激励源可分为光声成像、微波热声成像和磁热声成像等。光声成像是利用激光对目标体进行照射,目标体吸收能量而产生膨胀,进而产生超声信号,光声成像深度可到毫米级别。微波热声是利用微波代替激光,对目标体进行辐照。微波热声成像深度可到厘米级别。磁热声利用低频磁场(20MHz左右)激励目标体,在目标体内产生感应电场,产生焦耳热,进而产生热膨胀而产生超声信号。磁热声成像具有更深的成像深度,可达10厘米以上。
对于生物组织成像,微波热声和磁热声的成像深度虽然能到厘米级别,但微波热声和激励源功率需要数十千瓦,磁热声的激励磁场由大脉冲功率源驱动线圈产生,大功率的激励源笨重、功耗大、体积大、成本高,不方便仪器的设计和使用。
医学中人体组织发生病变,其功能病变要早于形态病变,在形态结构尚未变化之前,病变组织的电特性已经发生变化。因此,通过对组织电特性的检测,可实现疾病的早期诊断和早期治疗。电阻抗成像是一种对人体内部组织电阻抗进行成像的技术,具有无损害、低成本等特点。三十多年来,已经发展了多种电阻抗成像技术。传统电阻抗成像包括注入式电阻抗成像、感应式电阻抗成像、磁感应成像等,但传统电阻抗成像的测量信息对物体内部的电导率反映灵敏度低,加之测量数据有限以及测量噪声的影响,重建电导率的逆问题存在严重的病态性,尽管可以利用各种数学技巧来处理逆问题,但难以从本质上解决算法的稳定性、收敛性等问题,因此成像分辨率低。近些年提出的核磁共振电阻抗成像和磁声成像。核磁共振电阻抗成像设备昂贵;磁声成像也需要大功率磁激励源,也存在激励源笨重、功耗大、体积大的问题。
发明内容
为了克服上述现有方法的不足,本发明提出一种注入电流式热声成像方法。
本发明采用电极注入电流作为激励。通过电极向成像目标体注入电流,在成像目标体中 产生焦耳热,引起热膨胀,产生超声信号,根据检测的超声信号,重建导电物体的电导率分布。
本发明的目标体产生的超声信号是由注入的电流产生,注入电流I为脉冲电流,脉冲电流的频谱在微波以下,主要分布在100MHz以内,采用较低的频率具有更深的成像深度。成像目标体的电导率分布σ,注入电流I在成像目标体中产生的电场强度E,电流密度J=σE,成像目标体中产生的热量主要是目标体电导率引起的电导损耗,即焦耳热,其单位体积的功率密度为σE2,目标体中产生的超声声压满足方程:
式中,p为声压,c为声速,β为等压膨胀系数,Cp为物体的比热容,t为时间,▽为拉普拉斯算子。
本发明的成像方法步骤如下:
步骤1:通过电极向成像目标体注入电流,通过布置在目标体周边的超声探头接收产生的超声信号;
首先通过电极对目标体施加脉冲电流激励,目标体内由于电流感应电场的作用产生焦耳热,进而产生热声信号;所产生的热声信号通过耦合剂传播到超声探头;超声探头接收到超声信号后通过对超声信号处理、采集子系统进行前置放大、滤波、二级放大处理后,进行存储;
步骤2:利用接收的超声信号,重建声源分布;
已知注入电流式热声成像的声压波动方程:
式中,p(r,t)是在时刻t和位置r的声压,cs是热声源在耦合剂与目标体中的传播声速,CP和β分别表示目标体的比热容和体积膨胀系数,Q(r′)是r′位置的热声源分布;
对方程(1)利用时间反演法进行求解,得到热声源的重建公式:
式中,r为超声换能器位置、r′为热声源位置、S为超声换能器的检测面,n为检测面上超声换能器位置r的单位矢量。
选取目标体的某一断层面z=z0,断层面z=z0上利用时间反演法对方程(1)求解,获取热声源Q(x,y,z0)。目标体上的热声源Q(x,y,z)通过断层数据在z方向的插值得到,或者通过 移动超声换能器分层计算得到;
步骤3:获取目标体标量电位空间分量
热声源分布同时是电导率和电场强度的函数,因此热声源表示为:
Q=σE2 (3)
式(3)中,Q为热声源分布,σ为目标体的电导率,E为目标体内电场强度分布;
本发明所述的注入电流式热声成像方法主要用于生物组织的成像,考虑到生物组织的电导率较低,电场强度的空间分布表示为:
E=-▽φ (4)
式(4)中,φ是标量电势,▽是哈密顿算符;
根据电流连续性定理,有:
▽·(σ▽φ)=0 (5)
对目标体进行空间离散,给出电导率的初值[σ]0,将[σ]0代入公式(5),利用线性有限元方法求解得到标量电位的空间分量[φ]0。
步骤4:最优化迭代重建成像目标体电导率
由公式(3)和公式(4)得到:
定义优化目标函数:
f(σ)=[Q-σ(▽φ)·(▽φ)]2 (7)
将标量电位的空间分量[φ]0代入优化目标函数中,通过最优化算法进行迭代,直到误差满足终止条件,所得的电导率即为最终目标体的电导率。
本发明中,所述的激励是注入脉冲电流,检测的是超声信号,重建的是目标体的电阻率,是一种高分辨的电阻抗成像方法。
可根据实际情况,选择注入电流的脉冲宽度和脉冲电流类型。对于浅层成像,对分辨率的要求较高时,采用较窄脉冲;对于较深目标体成像,由于超声的传播衰减,采用较宽脉冲;对于信噪比较低时,采用编码脉冲、脉冲压缩技术提高信噪比。
本发明采用电极注入电流激励目标体产生超声信号,激励源的功率小、轻便、成本低、使用方便、灵活、可产生复杂波形,信噪比好、成像质量高。可实现低成本、高质量的高分辨率的电阻抗成像,可用于疾病的早期诊断。
附图说明
图1本发明原理示意图;其中,1激励源,2导线,3注入电流的电极,4成像目标体,5成像目标体中的不同电导率物体,6超声探头阵列;
图2本发明应用实例图;其中,7激励源,8注入电流的电极,9生物组织体,10耦合剂,11超声探头,12放大采集器,13成像计算机。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步说明。
本发明原理如图1所示,激励源1通过电极3向目标体4注入电流作为激励,在目标体4中产生焦耳热,引起热膨胀,产生超声信号,利用超声探头阵列6检测超声信号,根据检测的超声信号重建声源和目标体电导率分布。
激励源1通过电极3向目标体4注入电流I为脉冲电流,脉冲电流的频谱在微波以下,主要分布在100MHz以内,采用较低的频率将具有更深的成像深度,对于医学成像应用,根据生物组织的电阻抗特性,选择脉冲宽度为0.1μs~10μs,目标体4的电导率分布σ,注入电流I在目标体4中产生的电场强度E,电流密度J=σE,目标体3中产生的热量主要是目标体4的电导率引起的电导损耗,即焦耳热,其单位体积的功率密度为σE2,产生的超声声压满足方程:其中,p为声压,c为声速,β为等压膨胀系数,Cp为物体的比热容,t为时间,▽为拉普拉斯算子。
利用注入电流式热声实现电阻抗成像步骤如下:
步骤1:通过电极向成像目标体注入电流,通过布置在目标体周边的超声探头接收产生的超声信号;
首先通过电极对目标体施加脉冲电流激励,目标体内由于电流感应电场的作用产生焦耳热,进而产生热声信号;所产生的热声信号通过耦合剂传播到超声探头;超声探头接收到超声信号后通过对超声信号处理、采集子系统进行前置放大、滤波、二级放大处理后,进行存储;
步骤2:利用接收的超声信号,重建声源分布;
已知注入电流式热声成像的声压波动方程:
其中p(r,t)是在时刻t和位置r的声压,cs是热声源在耦合剂与目标体中的传播声速,CP和β分别表示目标体的比热容和体积膨胀系数,Q(r′)是r′位置的热声源分布;
对方程(1)利用时间反演法进行求解,得到热声源的重建公式:
其中r为超声换能器位置、r′为热声源位置、S为超声换能器的检测面和n为检测面上超声换能器位置r的单位矢量。
选取目标体的某一断层面z=z0,断层面z=z0上利用时间反演法对方程(1)求解,获取热声源Q(x,y,z0)。目标体上的热声源目标体上的热声源Q(x,y,z)通过断层数据在z方向的插值得到,或者通过移动超声换能器分层计算得到;
步骤3:获取目标体标量电位空间分量;
热声源分布同时是电导率和电场强度的函数,因此热声源表示为:
Q=σE2 (3)
其中,Q为热声源分布,σ为目标体的电导率,E为目标体内电场强度分布;
本发明所述的注入电流式热声成像方法主要用于生物组织的成像,考虑到生物组织的电导率较低,电场强度的空间分布表示为:
E=-▽φ (4)
其中φ是标量电势,▽是哈密顿算符;
根据电流连续性定理,有:
▽·(σ▽φ)=0 (5)
对目标体进行空间离散,给出电导率的初值[σ]0,将[σ]0代入公式(5),利用线性有限元方法求解得到标量电位的空间分量[φ]0;
步骤4:最优化迭代重建成像目标体电导率;
由公式(2)和公式(3)得到:
定义优化目标函数:
f(σ)=[Q-σ(▽φ)·(▽φ)]2 (7)
将标量电位的空间分量[φ]0代入优化目标函数中,通过最优化算法进行迭代,直到误差满足终止条件,所得的电导率即为最终目标体的电导率。
可根据实际情况,对注入电流的脉冲宽度、脉冲电流类型进行选择。对于浅层成像,对分辨率的要求较高时,采用较窄脉冲;对于较深目标体成像,由于超声的传播衰减,采用较宽脉冲;对于信噪比较低时,采用编码脉冲、脉冲压缩技术提高信噪比。
电极2可根据目标体3的结构、外形、电阻率进行选择,选用接触面较大的电极,减小 接触电阻,同时使得电流在目标体3的分布较为均匀。
如图2所示,在医学中,人体组织发生病变,其功能病变要早于形态病变,在形态结构尚未变化之前,病变组织的电特性已经发生变化。因此,通过对组织电特性的检测,可实现疾病的早期诊断和早期治疗,实现“治未病”的目的。利用注入电流式热声成像疾病诊断的应用实例如图2所示,激励源7的输出电压范围为0-1000V,输出电流范围为0-10A,调节输出电压和输出电流大小,激励电流脉冲宽度为0.5μs,电极8贴在生物组织体9的表面,激励源7通过电极8向生物组织体9注入电流作为激励,在生物组织体9中产生焦耳热,引起热膨胀,产生超声信号,将中心频率为1MHz超声相控阵探头11通过耦合剂10与生物组织体9接触,利用超声探头11检测超声信号,超声信号经过放大采集器12进行信号放大和采集,放大采集器12的放大倍数为60dB,采集速率为50Msps,将采集的数据通过以太网传输到计算机13,根据检测的超声信号进行超声成像以及重建电阻抗图像,在计算机13中显示成像结果,根据超声图像和电阻抗图像进行疾病诊断。
Claims (4)
1.一种注入电流式热声成像方法,其特征在于,所述的注入电流式热声成像方法通过电极向成像目标体注入电流,在成像目标体中产生焦耳热,引起热膨胀,产生超声信号,根据检测的超声信号,重建声源分布和电导率分布。
2.根据权利要求1所述的注入电流式热声成像方法,其特征在于,所述的超声信号由注入的电流产生,成像目标体的电导率分布σ,注入电流I在成像目标体中产生的电场强度E,电流密度J=σE,产生的超声声压满足方程:
其中,p为声压,c为声速,β为等压膨胀系数,Cp为物体的比热容,t为时间,为拉普拉斯算子。
3.根据权利要求1所述的注入电流式热声成像方法,其特征在于,所述的注入电流式热声成像方法步骤如下:
步骤1:通过电极向成像目标体注入电流,通过布置在目标体周边的超声探头接收产生的超声信号;
步骤2:利用接收的超声信号重建声源分布;
已知注入电流式热声成像的声压波动方程:
式中,p(r,t)是在时刻t和位置r的声压,cs是热声源在耦合剂与目标体中的传播声速,CP和β分别表示目标体的比热容和体积膨胀系数,Q(r′)是r′位置的热声源分布;
对方程(1)利用时间反演法进行求解,得到热声源的重建公式:
式中,r为超声换能器位置、r′为热声源位置、S为超声换能器的检测面,n为检测面上超声换能器位置r的单位矢量;
选取目标体的某一断层面z=z0,断层面z=z0上利用时间反演法对方程(1)求解,获取热声源Q(x,y,z0);目标体上的热声源Q(x,y,z)通过断层数据在z方向的插值得到,或者通过移动超声换能器分层计算得到;
步骤3:获取目标体标量电位空间分量;
热声源分布同时是电导率和电场强度的函数,因此热声源表示为:
Q=σE2 (3)
式(3)中,Q为热声源分布,σ为目标体的电导率,E为目标体内电场强度分布;
所述的注入电流式热声成像方法主要用于生物组织的成像,考虑到生物组织的电导率较低,电场强度的空间分布表示为:
式(4)中,φ是标量电势,是哈密顿算符;
根据电流连续性定理,有:
对目标体进行空间离散,给出电导率的初值[σ]0,将[σ]0代入公式(5),利用线性有限元方法求解得到标量电位的空间分量[φ]0;
步骤4:最优化迭代重建成像目标体电导率;
由公式(3)和公式(4)得到:
定义优化目标函数:
将标量电位的空间分量[φ]0代入优化目标函数中,通过最优化算法进行迭代,直到误差满足终止条件,所得的电导率即为最终目标体的电导率。
4.根据权利要求1所述的注入电流式热声成像方法,其特征在于,所述的注入电流为脉冲电流,脉冲电流的频谱分布在100MHz以下。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610991573.4A CN106580249B (zh) | 2016-11-08 | 2016-11-08 | 一种注入电流式热声成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610991573.4A CN106580249B (zh) | 2016-11-08 | 2016-11-08 | 一种注入电流式热声成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106580249A true CN106580249A (zh) | 2017-04-26 |
CN106580249B CN106580249B (zh) | 2019-07-05 |
Family
ID=58591109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610991573.4A Active CN106580249B (zh) | 2016-11-08 | 2016-11-08 | 一种注入电流式热声成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106580249B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108169332A (zh) * | 2017-12-14 | 2018-06-15 | 太原理工大学 | 基于热声成像的钢丝绳芯检测方法和装置 |
CN111887807A (zh) * | 2020-03-13 | 2020-11-06 | 中国石油大学(华东) | 一种纳秒脉冲电场诱导热声信号处理方法及系统 |
CN112443315A (zh) * | 2020-11-23 | 2021-03-05 | 中国科学院电工研究所 | 一种磁声电成像测井方法及其装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040127783A1 (en) * | 2002-09-26 | 2004-07-01 | Kruger Robert A. | Tissue scanner |
US20100041987A1 (en) * | 2006-12-21 | 2010-02-18 | Universiteit Twente | Imaging Apparatus And Method |
CN104434099A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 利用线性双旋度方程进行磁热声成像的电阻率重建方法 |
CN104434100A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 一种磁热声成像的电阻率重建方法 |
CN104473640A (zh) * | 2014-12-14 | 2015-04-01 | 中国科学院电工研究所 | 一种磁热声成像的电导率重建方法 |
-
2016
- 2016-11-08 CN CN201610991573.4A patent/CN106580249B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040127783A1 (en) * | 2002-09-26 | 2004-07-01 | Kruger Robert A. | Tissue scanner |
US20100041987A1 (en) * | 2006-12-21 | 2010-02-18 | Universiteit Twente | Imaging Apparatus And Method |
CN104434099A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 利用线性双旋度方程进行磁热声成像的电阻率重建方法 |
CN104434100A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 一种磁热声成像的电阻率重建方法 |
CN104473640A (zh) * | 2014-12-14 | 2015-04-01 | 中国科学院电工研究所 | 一种磁热声成像的电导率重建方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108169332A (zh) * | 2017-12-14 | 2018-06-15 | 太原理工大学 | 基于热声成像的钢丝绳芯检测方法和装置 |
CN111887807A (zh) * | 2020-03-13 | 2020-11-06 | 中国石油大学(华东) | 一种纳秒脉冲电场诱导热声信号处理方法及系统 |
CN111887807B (zh) * | 2020-03-13 | 2021-08-17 | 中国石油大学(华东) | 一种纳秒脉冲电场诱导热声信号处理方法及系统 |
CN112443315A (zh) * | 2020-11-23 | 2021-03-05 | 中国科学院电工研究所 | 一种磁声电成像测井方法及其装置 |
CN112443315B (zh) * | 2020-11-23 | 2023-09-26 | 中国科学院电工研究所 | 一种磁声电成像测井方法及其装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106580249B (zh) | 2019-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Imaging electrical impedance from acoustic measurements by means of magnetoacoustic tomography with magnetic induction (MAT-MI) | |
CN102640014B (zh) | 图像生成装置和图像生成方法 | |
Li et al. | Multi-excitation magnetoacoustic tomography with magnetic induction for bioimpedance imaging | |
Mariappan et al. | Magnetoacoustic tomography with magnetic induction: bioimepedance reconstruction through vector source imaging | |
Omar et al. | Near‐field thermoacoustic imaging with transmission line pulsers | |
Mariappan et al. | Magnetoacoustic tomography with magnetic induction for high‐resolution bioimepedance imaging through vector source reconstruction under the static field of MRI magnet | |
JP2015507947A5 (zh) | ||
Teniou et al. | A multimodal image reconstruction method using ultrasonic waves and electrical resistance tomography | |
CN106580249B (zh) | 一种注入电流式热声成像方法 | |
CN104473639B (zh) | 一种基于最优化迭代算法的磁热声成像电阻率重建方法 | |
Mariappan et al. | B-scan based acoustic source reconstruction for magnetoacoustic tomography with magnetic induction (MAT-MI) | |
CN100493442C (zh) | 基于声透镜与偏振检测的实时光声成像方法 | |
Huang et al. | Blood vessel imaging using radiofrequency-induced second harmonic acoustic response | |
Sun et al. | Acoustic dipole radiation based electrical impedance contrast imaging approach of magnetoacoustic tomography with magnetic induction | |
Song et al. | An integrated simulation approach and experimental research on microwave induced thermo-acoustic tomography system | |
CN104473640B (zh) | 一种磁热声成像的电导率重建方法 | |
CN104434099B (zh) | 利用线性双旋度方程进行磁热声成像的电阻率重建方法 | |
CN104586363A (zh) | 基于图像块稀疏系数的快速光声成像图像重建方法 | |
CN111948291A (zh) | 一种纳秒脉冲电场诱导热声成像的电导率重建方法 | |
Deng et al. | Image quality improvement of magneto-acousto-electrical tomography with Barker coded excitation | |
CN106885842B (zh) | 一种注入电流式热声成像电阻率重建方法 | |
Hakakzadeh et al. | Unipolar back-projection algorithm for photoacoustic tomography | |
Bitton et al. | A 3-D high-frequency array based 16 channel photoacoustic microscopy system for in vivo micro-vascular imaging | |
Zheng et al. | A computer-based simulator for intravascular photoacoustic images | |
CN104434100B (zh) | 一种磁热声成像的电阻率重建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |