CN107064302A - 一种注入电流式热声成像电导率重建方法 - Google Patents
一种注入电流式热声成像电导率重建方法 Download PDFInfo
- Publication number
- CN107064302A CN107064302A CN201710089729.4A CN201710089729A CN107064302A CN 107064302 A CN107064302 A CN 107064302A CN 201710089729 A CN201710089729 A CN 201710089729A CN 107064302 A CN107064302 A CN 107064302A
- Authority
- CN
- China
- Prior art keywords
- mrow
- formula
- signal
- mfrac
- injection current
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/24—Probes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02483—Other human or animal parts, e.g. bones
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02863—Electric or magnetic parameters
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
一种注入电流式热声成像电导率重建方法,基于电流注入方式的热声成像原理。通过电极向成像目标体注入脉冲电流,在成像目标体中产生焦耳热,引起热膨胀,产生超声信号,用超声换能器接收超声信号,对接收到的超声信号进行处理和采集,采用电导率图像重建算法获取目标体的电导率图像。具体步骤为:1、首先获取注入电流式热声信号;2、利用获取的注入电流式热声信号重建目标体热声源;3、利用热声源,采用非线性有限元求解方法重建标量电位;4、利用重建的标量电位重建电导率。
Description
技术领域
本发明涉及一种电导率图像重建方法,特别涉及一种注入电流式热声成像的电导率重建方法。
背景技术
由于激励频率的限制,传统电阻抗成像技术的灵敏度和空间分辨率不高。单一场都有其物理局限性,多物理场成像由一种物理场提供分辨率,另一种物理场提供对比度,实现对比度和分辨率的同时提高。电磁场和超声相结合的多物理场成像技术正是考虑到电磁场对人体组织电导率的高对比度和超声波探测的高分辨率特性,成为人们研究的热点,磁热声成像正是一种新兴的多物理场成像技术。
磁热声成像是由新加坡南洋理工大学的Feng在2013年首次提出的新型的电阻抗成像方法,其原理为:通过对导电成像体施加MHz量级的交变磁场,在目标体内部产生感应电场,进而产生焦耳热,激发热弹性的超声信号,检测超声信号进行成像。与微波热声成像相比,允许更低的功率进行高效的成像,并且具有便携式成像的潜力,同时,激励源的频率降低,使得磁场穿透组织更深,也避免的辐射。
作为一种新型的多物理场成像方法,2013年Feng利用金属铜仿体,检测到磁热声信号,并得到铜仿体的热声图像,并未进行电导率图像重建,且生物组织不同于金属铜仿体,磁场和感应电流作用产生的洛伦兹力较弱。专利“一种磁热声成像的电导率重建方法”(201410773988.5)公布了一种基于磁热声效应的电导率重建方法,其在热函数基础上提出了电导率重建方法,依然为采用线圈激励方式。同时,外加激励线圈产生时变磁场在目标体中产生二次磁场和感应电流,磁场和感应电流作用同时产生洛伦兹力和焦耳热,即磁声效应和磁热声效应共存,如何区分磁声效应和磁热声效应是仍需解决的问题。
基于此,采用注入电流式热声成像方法,可以避开磁声效应和磁热声效应共存问题,注入电流式热声成像与磁热声成像在激励方式、从热函数到电导率的重建方法上均不相同。
发明内容
本发明的目的是克服现有磁热声成像存在的问题,提出一种基于注入电流式热声成像的电导率重建方法,利用热声源重建电导率。本发明避免了热声成像过程中磁声效应的干扰,同时采用注入电流式激励,可增强热声效应,可实现目标体电导率图像的精确重建。
注入电流式热声成像原理为:通过注入电极向成像目标体注入电流,在成像目标体中产生焦耳热,引起热膨胀,产生超声信号。采用超声换能器进行检测,根据检测的超声信号,重建热声源和电导率。
本发明注入电流式热声成像的电导率图像重建包括四个步骤:1、首先获取注入电流式热声信号,即检测超声信号;2、利用获取的注入电流式热声信号重建目标体热声源;3、利用热声源,采用非线性有限元求解方法重建标量电位;4、利用重建的标量电位重建电导率。
具体描述如下:
第一步:获取注入电流式热声信号
信号发生器和功率放大器组成脉冲激励源,通过注入电极A和注入电极B对目标体注入电流,目标体在电流作用下产生焦耳热,进而产生热膨胀,激发超声信号,超声信号通过耦合剂耦合到超声换能器,超声换能器接收到信号后通过信号处理器进行信号的放大和滤波,经数据采集***后进行存储,超声换能器在控制器控制下对目标体进行扫描检测;
第二步:获取目标体热声源
已知热声成像的声压波动方程:
其中r为超声换能器位置,p(r,t)是声压,cs为介质中的声速,CP为目标体(5)的比热容,β为目标体的热膨胀系数,δ(t)是狄拉克函数,S(r)为热声源分布,t为时间,为拉普拉斯算符;
热声源分布的时间反演法重建公式为:
其中,R为标量,R=|r′-r|,R为矢量,eR为单位矢量,r'是超声探头的位置,r为热声源位置,Sd是超声探头所在的平面,p′是声压对时间的一阶导数,n是r′位置Sd的法线矢量,β为目标体的热膨胀系数;
选取目标体的某一断层面z=z1,超声换能器在此断层面上进行圆周扫描,采集超声信号,利用方程(1)和公式(2),求解z=z1断层面上的热声源分布S(x,y,z1),移动超声换能器在不同断层面上进行扫描检测,求解不同断层面上的热声源分布,目标体上整体热声源S可由分层检测计算或z方向上插值得到;
热声源S同时是电导率和电场强度的函数,可以表示为:
S=σE2=σE·E (3)
其中,σ为目标体的电导率,E为目标体内电场强度分布;
第三步:求解标量电位
所述的注入电流式热声成像方法,电场强度的空间分布表示为:
式(4)中,φ是标量电位,是哈密顿算符;
针对生物组织,采用电准静态近似,根据电流连续性定理,有:
其中,为散度符号,是标量电位φ的梯度;
由公式(3)和(4),可得到
将公式(6)代入公式(5)中
满足的边界条件为:
其中,ΓA,B为注入电极位置,Γg为除去电极之外的目标体边界,A0为电极与目标体接触面积,I为注入电流,为φ的法向导数;
将热声源S代入公式(7),结合边界条件公式(8),进行有限元法求解,即可重建得到标量电位φ;
第四步:求解电导率
将标量电位φ代入公式(6),即可重建电导率σ。
附图说明
图1本发明重建方法所涉及的注入电流式热声信号获取装置示意图;
图中:1信号发生器、2功率放大器、3水槽、4注入电极A、5目标体、6注入电极B、7超声换能器、8信号处理器、9数据采集***、10图像重建模块、11控制器。
具体实施方式
以下结合附图和具体实施方式对本发明做进一步说明。
本发明重建方法所涉及的注入电流式热声信号获取装置主要包括激励***、检测***、控制器和目标体四个部分。如图1所示,所述的激励***包括信号发生器1、功率放大器2、注入电极A4和注入电极B6。所述的检测***包括超声换能器7、信号处理器8、数据采集***9和图像重建模块10。信号发生器1驱动功率放大器2,通过注入电极A4和注入电极B6向目标体5注入电流信号。超声换能器7与目标体5之间通过耦合剂耦合,超声换能器7的输出端连接信号处理器8的输入端,信号处理器8的输出端连接数据采集***9的输入端,数据采集***9连接图像重建模块10。控制器11实现对超声换能器7的旋转扫描运动控制。
注入电流式热声成像原理为:通过注入电极A4和注入电极B6向成像目标体5注入电流,在成像目标体5中产生焦耳热,引起热膨胀,产生超声信号,采用超声换能器7进行检测超声信号,根据检测的超声信号,由图像重建模块10重建热声源和电导率。
本发明注入电流式热声成像的电导率图像重建包括四个步骤:1、首先获取注入电流式热声信号,即检测超声信号;2、利用获取的超声信号重建目标体热声源;3、利用热声源,采用非线性有限元求解方法重建标量电位;4、利用重建的标量电位重建电导率。
图像重建的具体过程描述如下:
第一步:获取注入电流式热声信号
信号发生器1和功率放大器2组成脉冲激励源,脉冲激励源通过注入电极A4和注入电极B6对目标体5注入电流,目标体5在电流作用下产生焦耳热,进而产生热膨胀,激发超声信号,超声信号通过耦合剂耦合到超声换能器7,超声换能器7接收到信号后,通过信号处理器8进行信号的放大和滤波,经数据采集***9后进行存储,超声换能器7在控制器11控制下对目标体5进行扫描检测;
第二步:获取目标体热声源
已知热声成像的声压波动方程:
其中r为超声换能器位置,p(r,t)是声压,cs为介质中的声速,CP为目标体(5)的比热容,β为目标体(5)的热膨胀系数,δ(t)是狄拉克函数,S(r)为热声源分布,t为时间,为拉普拉斯算符;
热声源分布的时间反演法重建公式为:
其中,R为标量,R=|r′-r|,R为矢量,eR为单位矢量,r'是超声探头的位置,r为热声源位置,Sd是超声探头所在的平面,p′是声压对时间的一阶导数,n是r′位置Sd的法线矢量,β为目标体的热膨胀系数;
选取目标体的某一断层面z=z1,超声换能器在此断层面上进行圆周扫描,采集超声信号,利用方程(1)和公式(2),求解z=z1断层面上的热声源分布S(x,y,z1),移动超声换能器在不同断层面上进行扫描检测,求解不同断层面上的热声源分布,目标体上整体热声源S可由分层检测计算或z方向上插值得到;
热声源S同时是电导率和电场强度的函数,可以表示为:
S=σE2=σE·E (3)
其中,σ为目标体的电导率,E为目标体内电场强度分布;
第三步:求解标量电位
所述的注入电流式热声成像方法,电场强度的空间分布表示为:
式(4)中,φ是标量电位,是哈密顿算符;
针对生物组织,采用电准静态近似,根据电流连续性定理,有:
其中,为散度符号,是标量电位φ的梯度;
由公式(3)和(4),可得到
将公式(6)代入公式(5)中
满足的边界条件为:
其中,ΓA,B为注入电极位置,Γg为除去电极之外的目标体边界,A0为电极与目标体接触面积,I为注入电流,为φ的法向导数;
将热声源S代入公式(7),结合边界条件(8),进行有限元法求解,即可重建得到标量电位φ;
第四步:求解电导率
将标量电位φ代入公式(6),即可重建电导率σ。
Claims (1)
1.一种注入电流式热声成像电导率重建方法,所述的注入电流式热声电导率图像重建方法基于注入电流式热声成像原理,通过电极向成像目标体注入脉冲电流,在成像目标体中产生焦耳热,引起热膨胀,产生超声信号,用超声换能器接收超声信号,对接收到的超声信号进行处理和采集,采用电导率图像重建算法获取目标体的电导率图像,其特征在于:所述的注入电流式热声成像电导率重建方法包括以下步骤:
第一步:获取注入电流式热声信号
信号发生器(1)和功率放大器(2)组成脉冲激励源,通过注入电极A(4)和注入电极B(6)对目标体(5)注入电流,目标体(5)在电流作用下产生焦耳热,进而产生热膨胀,激发超声信号,超声信号通过耦合剂耦合到超声换能器(7),超声换能器(7)接收到信号后通过信号处理器(8)进行信号的放大和滤波,经数据采集***(9)后进行存储,由图像重建模块(10)重建热声源和电导率,超声换能器(7)在控制器(11)控制下对目标体(5)进行扫描检测;
第二步:获取目标体热声源
根据第一步采用一对注入电极A、注入电极B,向目标体注入电流后,引发热膨胀,激发声信号,检测到声压信号,采用声压信号重建目标体上整体热声源;
第三步:求解标量电位
所述的注入电流式热声成像方法,电场强度的空间分布表示为:
<mrow>
<mi>E</mi>
<mo>=</mo>
<mo>-</mo>
<mo>&dtri;</mo>
<mi>&phi;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式(4)中,φ是标量电位,是哈密顿算符;
针对生物组织,采用电准静态近似,根据电流连续性定理,有:
<mrow>
<mo>&dtri;</mo>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>&sigma;</mi>
<mo>&dtri;</mo>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为散度符号,是标量电位φ的梯度;
由公式S=σE2=σE·E和公式(4),得到
<mrow>
<mi>&sigma;</mi>
<mo>=</mo>
<mfrac>
<mi>S</mi>
<mrow>
<mi>E</mi>
<mo>&CenterDot;</mo>
<mi>E</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mi>S</mi>
<mrow>
<mo>&dtri;</mo>
<mi>&phi;</mi>
<mo>&CenterDot;</mo>
<mo>&dtri;</mo>
<mi>&phi;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
将公式(6)代入公式(5)中
<mrow>
<mo>&dtri;</mo>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mi>S</mi>
<mrow>
<mo>&dtri;</mo>
<mi>&phi;</mi>
<mo>&CenterDot;</mo>
<mo>&dtri;</mo>
<mi>&phi;</mi>
</mrow>
</mfrac>
<mo>&dtri;</mo>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
满足的边界条件为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mi>S</mi>
<mrow>
<mo>&dtri;</mo>
<mi>&phi;</mi>
<mo>&CenterDot;</mo>
<mo>&dtri;</mo>
<mi>&phi;</mi>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&phi;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>A</mi>
<mo>,</mo>
<mi>B</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<mo>&PlusMinus;</mo>
<mfrac>
<mi>I</mi>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&phi;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<msub>
<mi>&Gamma;</mi>
<mi>g</mi>
</msub>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,ΓA,B为注入电极位置,Γg为除去电极之外的目标体边界,A0为电极与目标体接触面积,I为注入电流,为φ的法向导数;
将热声源S代入公式(7),结合边界条件公式(8),进行有限元法求解,即可重建得到标量电位φ;
第四步:求解电导率
将标量电位φ代入公式(6),即可重建电导率σ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710089729.4A CN107064302B (zh) | 2017-02-20 | 2017-02-20 | 一种注入电流式热声成像电导率重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710089729.4A CN107064302B (zh) | 2017-02-20 | 2017-02-20 | 一种注入电流式热声成像电导率重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107064302A true CN107064302A (zh) | 2017-08-18 |
CN107064302B CN107064302B (zh) | 2019-10-11 |
Family
ID=59621760
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710089729.4A Active CN107064302B (zh) | 2017-02-20 | 2017-02-20 | 一种注入电流式热声成像电导率重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107064302B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111458566A (zh) * | 2020-05-20 | 2020-07-28 | 重庆文理学院 | 一种储能材料电导率的非接触式检测方法及*** |
CN111887807A (zh) * | 2020-03-13 | 2020-11-06 | 中国石油大学(华东) | 一种纳秒脉冲电场诱导热声信号处理方法及*** |
CN111948291A (zh) * | 2020-03-13 | 2020-11-17 | 中国石油大学(华东) | 一种纳秒脉冲电场诱导热声成像的电导率重建方法 |
CN112694974A (zh) * | 2020-11-26 | 2021-04-23 | 中国石油大学(华东) | 针对纳秒脉冲电场消融动态监测***构建及监测方法 |
CN112914539A (zh) * | 2021-03-12 | 2021-06-08 | 中国科学院电工研究所 | 一种磁热声温度成像方法与装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040267346A1 (en) * | 2003-06-30 | 2004-12-30 | Shelso Susan I. | Stent grip and system for use therewith |
CN104458818A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 一种基于线性泊松方程的磁热声成像电导率重建方法 |
CN104473640A (zh) * | 2014-12-14 | 2015-04-01 | 中国科学院电工研究所 | 一种磁热声成像的电导率重建方法 |
-
2017
- 2017-02-20 CN CN201710089729.4A patent/CN107064302B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040267346A1 (en) * | 2003-06-30 | 2004-12-30 | Shelso Susan I. | Stent grip and system for use therewith |
CN104458818A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 一种基于线性泊松方程的磁热声成像电导率重建方法 |
CN104473640A (zh) * | 2014-12-14 | 2015-04-01 | 中国科学院电工研究所 | 一种磁热声成像的电导率重建方法 |
Non-Patent Citations (4)
Title |
---|
HUANG L ET AL.: "Quantitative thermoacoustic tomography:Recovery of conductivity maps of heterogeneous media", 《APPLIED PHYSICS LETTERS》 * |
XU MINGHUA ET AL: "Analytic explanation of spatial resolution related to bandwidth and detector aperture size in thermoacoustic or photoacoustic reconstruction", 《PHYSICAL REVIEW》 * |
张顺起 等: "注入电流式磁声成像的电导率模型构建和实验研究", 《中国生物医学工程学报》 * |
李俊霖 等: "注入式磁声成像中热声效应的初探", 《生物医学工程研究》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111887807A (zh) * | 2020-03-13 | 2020-11-06 | 中国石油大学(华东) | 一种纳秒脉冲电场诱导热声信号处理方法及*** |
CN111948291A (zh) * | 2020-03-13 | 2020-11-17 | 中国石油大学(华东) | 一种纳秒脉冲电场诱导热声成像的电导率重建方法 |
CN111458566A (zh) * | 2020-05-20 | 2020-07-28 | 重庆文理学院 | 一种储能材料电导率的非接触式检测方法及*** |
CN112694974A (zh) * | 2020-11-26 | 2021-04-23 | 中国石油大学(华东) | 针对纳秒脉冲电场消融动态监测***构建及监测方法 |
CN112914539A (zh) * | 2021-03-12 | 2021-06-08 | 中国科学院电工研究所 | 一种磁热声温度成像方法与装置 |
CN112914539B (zh) * | 2021-03-12 | 2022-08-05 | 中国科学院电工研究所 | 一种磁热声温度成像方法与装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107064302B (zh) | 2019-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107064302A (zh) | 一种注入电流式热声成像电导率重建方法 | |
CN102788836B (zh) | 一种磁声显微成像方法及成像*** | |
CN105816156B (zh) | 一种结合热声效应信号的电导率磁声成像装置及成像方法 | |
CN102085096B (zh) | 注入电流式磁声耦合成像装置 | |
CN102805621B (zh) | 一种磁声电成像***及成像方法 | |
CN101791219A (zh) | 一种磁声电阻抗成像方法及装置 | |
CN104146708B (zh) | 一种基于声换能器接收特性的磁声图像重建方法及装置 | |
CN102894974A (zh) | 一种磁声电成像***及成像方法 | |
Liu et al. | Magnetoacoustic tomography with current injection | |
CN104434094B (zh) | 一种磁‑热‑声耦合成像的电导率图像重建方法 | |
CN104013388B (zh) | 基于低频连续波的磁声耦合成像激励与检测方法及装置 | |
CN106037638A (zh) | 一种去除热声效应影响的电导率磁声成像装置及成像方法 | |
CN106885842B (zh) | 一种注入电流式热声成像电阻率重建方法 | |
CN104434099B (zh) | 利用线性双旋度方程进行磁热声成像的电阻率重建方法 | |
CN102590625A (zh) | 磁声耦合成像的微弱声信号频域检测处理方法 | |
CN111948291B (zh) | 一种纳秒脉冲电场诱导热声成像的电导率重建方法 | |
CN104473640B (zh) | 一种磁热声成像的电导率重建方法 | |
CN104434101B (zh) | 一种磁热声成像方法及其成像*** | |
CN104473639A (zh) | 一种基于最优化迭代算法的磁热声成像电阻率重建方法 | |
CN107049315A (zh) | 基于最优化迭代方法的注入电流式热声电阻率图像重建方法 | |
CN104458818B (zh) | 一种基于线性泊松方程的磁热声成像电导率重建方法 | |
CN106580249B (zh) | 一种注入电流式热声成像方法 | |
CN104434100B (zh) | 一种磁热声成像的电阻率重建方法 | |
CN103837580B (zh) | 一种基于超声和电磁超声相结合的双模无损检测方法 | |
CN103142216B (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 |