CN107049315A - 基于最优化迭代方法的注入电流式热声电阻率图像重建方法 - Google Patents
基于最优化迭代方法的注入电流式热声电阻率图像重建方法 Download PDFInfo
- Publication number
- CN107049315A CN107049315A CN201710089515.7A CN201710089515A CN107049315A CN 107049315 A CN107049315 A CN 107049315A CN 201710089515 A CN201710089515 A CN 201710089515A CN 107049315 A CN107049315 A CN 107049315A
- Authority
- CN
- China
- Prior art keywords
- mrow
- thermoacoustic
- resistivity
- formula
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/053—Measuring electrical impedance or conductance of a portion of the body
-
- 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
- Health & Medical Sciences (AREA)
- Life Sciences & Earth 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)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
一种基于最优化迭代方法的注入电流式热声电阻率图像重建方法,通过电极向成像目标体注入电流,在成像目标体中产生焦耳热,引起热膨胀,产生注入电流式热声信号:利用超声换能器获取目标体各断层热声信号,根据注入电流式热声声压波动方程,获取目标体每个断层上的热声源分布,利用插值法获取目标体整体的热声源分布;然后对导电物体的电阻率进行空间离散,设定电阻率的初值,根据欧姆定律,利用线性有限元方程重建目标体的矢量电位;再根据热声源和电流密度关系,结合电流连续性定理,获取电阻率所满足的关于热声源与矢量电位的方程,定义目标函数;最后将求解的矢量电位代入目标函数,利用目标函数最小化原则,求解导电物体的电阻率。
Description
技术领域
本发明涉及一种电阻率图像重建方法,特别涉及一种注入电流式电阻率图像重建方法。
背景技术
由于激励频率的限制,传统电阻抗成像技术的灵敏度和空间分辨率不高。单一场都有其物理局限性,多物理场成像由一种物理场提供分辨率,另一种物理场提供对比度,实现对比度和分辨率的同时提高。电磁场和超声相结合的多物理场成像技术正是考虑到电磁场对人体组织电阻率的高对比度和超声波探测的高分辨率特性,成为人们研究的热点,磁热声成像正是一种新兴的多物理场成像技术。
磁热声成像是由新加坡南洋理工大学的Feng在2013年首次提出的新型的电阻抗成像方法,其原理为:通过对导电成像体施加MHz量级的交变磁场,在目标体内部产生感应电场,进而产生焦耳热,激发热弹性的超声信号,检测超声信号进行成像。与微波热声成像相比,允许更低的功率进行高效的成像,并且具有便携式成像的潜力,同时,激励源的频率降低,使得磁场穿透组织更深,也避免的辐射。
作为一种新型的多物理场成像方法,2013年Feng利用金属铜仿体,检测到磁热声信号,并得到铜仿体的热声图像,并未进行电阻率图像重建,且生物组织不同于金属铜仿体,磁场和感应电流作用产生的洛伦兹力较弱。专利“一种磁热声成像的电阻率重建方法”(201410771496.2)公布了一种基于磁热声效应的电阻率重建方法,其在热函数基础上提出了电阻率重建方法,依然为采用线圈激励方式。同时,外加激励线圈产生时变磁场在目标体中产生二次磁场和感应电流,磁场和感应电流作用同时产生洛伦兹力和焦耳热,即磁声效应和磁热声效应共存,如何区分磁声效应和磁热声效应是仍需解决的问题。
基于此,采用注入电流式热声成像方法,可以避开磁声效应和磁热声效应共存问题,注入电流式热声成像与磁热声成像在激励方式、从热函数到电阻率的重建方法上均不相同。
发明内容
本发明的目的是克服现有磁热声成像存在的问题,提出一种基于最优化迭代方法的注入电流式热声电阻率图像重建方法。本发明避免了热声成像过程中磁声效应的干扰,同时采用注入电流式激励,可增强热声效应,也实现了目标体电阻率的图像重建。
注入电流式热声成像原理为:通过注入电极向成像目标体注入电流,在成像目标体中产生焦耳热,引起热膨胀,产生超声信号,采用超声换能器进行检测,根据检测的超声信号,重建热声源分布和电阻率。
本发明基于最优化迭代方法的注入电流式热声电阻率图像重建方法具体步骤如下:
第一步为获取注入电流式热声信号:利用超声换能器获取目标体各断层热声信号,每个断层扫描一周进行热声信号的检测,即检测超声信号;第二步根据注入电流式热声声压波动方程,获取目标体每个断层上的热声源分布,利用插值法获取目标体整体的热声源分布;第三步对导电物体的电阻率进行空间离散,设定电阻率的初值,根据欧姆定律,利用线性有限元方程重建的到目标体的矢量电位;第四步根据热声源和电流密度关系,结合电流连续性定理,获取电阻率所满足的关于热声源与矢量电位的方程,定义目标函数;第五步将求解的矢量电位代入目标函数,利用目标函数最小化原则,求解导电物体的电阻率。具体分述如下:
第一步:获取导电物体热声信号
激励源通过注入电极A和注入电极B,对目标体注入脉冲电流,目标体在电流作用下产生焦耳热,进而产生热膨胀,激发超声信号,超声信号通过耦合剂耦合到超声换能器,超声换能器接收到信号后,通过检测系统进行信号的放大、滤波、采集和存储,超声换能器在控制器控制下对目标体进行扫描检测;
第二步:求解目标体热声源分布
已知热声成像的声压波动方程:
其中r为超声换能器位置,p(r,t)是声压,cs为介质中的声速,CP为目标体(3)的比热容,β为目标体(3)的热膨胀系数,δ(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)。移动超声换能器在不同断层面上进行扫描检测,求解不同断层面上的热声源分布,目标体(3)上整体热声源分布可由分层检测计算或z方向上插值得到;
第三步:获取矢量电位
热声源S同时是电阻率和电流密度的函数,热声源可以表示为:
S=ρJ2=ρJ·J (3)
其中,ρ为目标体(3)的电阻率,J为目标体(3)内电流密度分布;
根据电流连续性定理,电流的散度为零,引入矢量电位,则有:
其中,T为矢量电位,为矢量电位的旋度,为哈密顿算符,J为目标体(3)内电流密度分布;
利用欧姆定律,J=σE=E/ρ,则有:
满足的边界条件为:
其中,ΓA,B为注入电极位置,Γg为除去电极之外的目标体边界,A0为电极与目标体接触面积,I为注入电流,n表示区域的外法向单位矢量;
对导电物体进行空间离散,并给出电阻率的初值[ρ]0,将电阻率的初值[ρ]0代入公式(5)中,考虑在成像体的边界施加边界条件(6),采用有限元方法求解矢量电位[T]1;
第四步:定义目标函数
由式(3)和式(4)可得:
根据公式(7)定义目标函数f(ρ):
其中,S(r)为通过测量声压由时间反演法计算得到的热声源分布,S为可以通过公式(5)、(6)和(7),第三步和第四步,由电阻率和矢量电位计算得到的热声源。
第五步:目标函数最优化迭代求解
将矢量电位[T]1代入目标函数,通过最优化算法更新电阻率进行迭代,直到误差满足最优终止条件ε为所要求的精度,最终求解目标体电阻率。
第五步的迭代方法为:
(1)若误差满足最优终止条件则计算过程结束,得到目标体电阻率;
(2)若不满足,采用高斯牛顿法或非线性阻尼最小二乘法,计算雅克比矩阵,求解电阻率参数增量△ρ;
(3)将电阻率参数增量△ρ代入上一个电阻率值[ρ]i,得到新的电阻率[ρ]i+1;
(4)计算新的矢量电位,反复迭代下去,直至误差满足最优终止条件
附图说明
图1本发明重建方法图像重建过程示意图;
图2本发明重建方法热声信号检测示意图;
图中:1激励源、2注入电极A、3目标体、4注入电极B、5超声换能器。
具体实施方式
以下结合附图和具体实施方式对本发明做进一步说明。
如图1所示,本发明基于最优化迭代方法的注入电流式热声电阻率图像重建方法步骤如下:
第一步为获取注入电流式热声信号:利用超声换能器获取目标体各断层热声信号,每个断层扫描一周进行热声信号的检测,即检测超声信号;第二步根据注入电流式热声声压波动方程,获取目标体每个断层上的热声源分布,利用插值法获取目标体整体的热声源分布;第三步对导电物体的电阻率进行空间离散,设定电阻率的初值,根据欧姆定律,利用线性有限元方程重建的到目标体的矢量电位;第四步根据热声源和电流密度关系,结合电流连续性定理,获取电阻率所满足的关于热声源与矢量电位的方程,定义目标函数;第五步将求解的矢量电位代入目标函数,利用目标函数最小化原则,求解导电物体的电阻率。具体分述如下:
第一步:获取导电物体热声信号
如图2所示,激励源1通过注入电极A2、注入电极B 4对目标体3注入脉冲电流,目标体3在电流作用下产生焦耳热,进而产生热膨胀,激发超声信号,超声信号通过耦合剂耦合到超声换能器5,超声换能器5接收到信号后通过检测系统进行信号的放大、滤波、采集 和存储,超声换能器5在控制器控制下对目标体进行扫描检测;
第二步:求解目标体热声源分布
已知热声成像的声压波动方程:
其中r为超声换能器位置,p(r,t)是声压,cs为介质中的声速,CP为目标体3的比热容,β为目标体3的热膨胀系数,δ(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)。移动超声换能器在不同断层面上进行扫描检测,求解不同断层面上的热声源分布,目标体3上整体热声源分布可由分层检测计算或z方向上插值得到;
第三步:获取矢量电位
热声源S同时是电阻率和电流密度的函数,可以表示为:
S=ρJ2=ρJ·J (3)
其中,ρ为目标体3的电阻率,J为目标体3内电流密度分布;
根据电流连续性定理,电流的散度为零,引入矢量电位,则有:
其中,T为矢量电位,为矢量电位的旋度,为哈密顿算符;
利用欧姆定律,J=σE=E/ρ,则有:
满足的边界条件为:
其中,ΓA,B为注入电极位置,Γg为除去电极之外的目标体边界,A0为电极与目标体接触面积,I为注入电流,n表示区域的外法向单位矢量;
对导电物体进行空间离散,并给出电阻率的初值[ρ]0,将电阻率的初值[ρ]0代入公式(5)中,考虑在成像体的边界施加边界条件(6),采用有限元方法求解矢量电位[T]1。
第四步:定义目标函数
由式(3)和式(4)可得:
根据公式(7)定义目标函数f(ρ):
其中,S(r)为通过测量声压由时间反演法计算得到的热声源分布,S为可以通过公式(5)、(6)和(7),以及第三步和第四步,由电阻率和矢量电位计算得到的热声源。
第五步:目标函数最优化迭代求解
将矢量电位[T]1代入目标函数,通过最优化算法更新电阻率进行迭代,直到误差满足最优终止条件ε为所要求的精度,最终求解目标体电阻率。
第五步的迭代方法为:
(1)若误差满足最优终止条件则计算过程结束,得到目标体电阻率;ε为所要求的精度;
(2)若不满足,采用高斯牛顿法或非线性阻尼最小二乘法,计算雅克比矩阵,求解电阻率参数增量△ρ;
(3)将电阻率参数增量代入上一个电阻率值[ρ]i,得到新的电阻率[ρ]i+1;
(4)计算新的矢量电位,反复迭代,直至误差满足最优终止条件
Claims (1)
1.一种基于最优化迭代方法的注入电流式热声电阻率图像重建方法,所述的注入电流式热声电导率图像重建方法基于注入电流式热声成像原理,所述的注入电流式热声成像方法通过电极向成像目标体注入电流,在成像目标体中产生焦耳热,引起热膨胀,产生超声信号,根据检测的超声信号,重建热声源分布和电阻率,其特征在于:所述的基于最优化迭代方法的注入电流式热声电阻率图像重建方法为:
第一步为获取注入电流式热声信号:利用超声换能器获取目标体各断层热声信号,每个断层扫描一周进行热声信号的检测;第二步根据注入电流式热声声压波动方程,获取目标体每个断层上的热声源分布,利用插值法获取目标体整体的热声源分布;第三步对导电物体的电阻率进行空间离散,设定电阻率的初值,根据欧姆定律,利用线性有限元方程重建目标体的矢量电位;第四步根据热声源和电流密度关系,结合电流连续性定理,获取电阻率所满足的关于热声源与矢量电位的方程,定义目标函数;第五步将求解的矢量电位代入目标函数,利用目标函数最小化原则,求解导电物体的电阻率;
所述的基于最优化迭代方法的注入电流式热声电阻率图像重建方法具体步骤如下:
第一步:获取导电物体热声信号
激励源(1)通过注入电极A(2)和注入电极B(4)对目标体(3)注入脉冲电流,目标体(3)在电流作用下产生焦耳热,进而产生热膨胀,激发超声信号,超声信号通过耦合剂耦合到超声换能器(5),超声换能器(5)接收到信号后通过检测系统进行信号的放大、滤波、采集和存储,超声换能器(5)在控制器控制下对目标体进行扫描检测;
第二步:求解目标体热声源分布
第一步采用注入电极A(2)和注入电极B(4),向目标体注入电流后,引发热膨胀,激发声信号,检测到声压信号,采用声压信号重建目标体上整体热声源分布;
第三步:获取矢量电位
根据电流连续性定理,电流的散度为零,引入矢量电位,则有:
<mrow>
<mi>J</mi>
<mo>=</mo>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,T为矢量电位,为矢量电位的旋度,为哈密顿算符,J为目标体(3)内电流密度分布;
利用欧姆定律,J=σE=E/ρ,则有:
<mrow>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>T</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>
满足的边界条件为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>T</mi>
<mo>&CenterDot;</mo>
<mi>n</mi>
<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>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>T</mi>
<mo>&CenterDot;</mo>
<mi>n</mi>
<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>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,ΓA,B为注入电极位置,Γg为除去电极之外的目标体边界,A0为电极与目标体接触面积,I为注入电流,n表示区域的外法向单位矢量;
对导电物体进行空间离散,并给出电阻率的初值[ρ]0,将电阻率的初值[ρ]0代入公式(5)中,考虑在成像体的边界施加边界条件(6),采用有限元方法求解矢量电位[T]1;
第四步:定义目标函数
由式S=ρJ2=ρJ·J和式(4)可得:
<mrow>
<mi>&rho;</mi>
<mo>=</mo>
<mfrac>
<mi>S</mi>
<mrow>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>T</mi>
<mo>&CenterDot;</mo>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>T</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
根据公式(7)定义目标函数f(ρ):
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,S(r)为通过测量声压由时间反演法计算得到的热声源分布,S为通过公式(5)、(6)和(7),以及第三步和第四步,由电阻率和矢量电位计算得到的热声源,ρ为目标体(3)的电阻率,J为目标体(3)内电流密度分布;
第五步:目标函数最优化迭代求解
将矢量电位[T]1代入目标函数,通过最优化算法更新电阻率进行迭代,直到误差满足最优终止条件ε为所要求的精度,最终求解目标体电阻率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710089515.7A CN107049315B (zh) | 2017-02-20 | 2017-02-20 | 基于最优化迭代方法的注入电流式热声电阻率图像重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710089515.7A CN107049315B (zh) | 2017-02-20 | 2017-02-20 | 基于最优化迭代方法的注入电流式热声电阻率图像重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107049315A true CN107049315A (zh) | 2017-08-18 |
CN107049315B CN107049315B (zh) | 2020-02-11 |
Family
ID=59621834
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710089515.7A Active CN107049315B (zh) | 2017-02-20 | 2017-02-20 | 基于最优化迭代方法的注入电流式热声电阻率图像重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107049315B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111887807A (zh) * | 2020-03-13 | 2020-11-06 | 中国石油大学(华东) | 一种纳秒脉冲电场诱导热声信号处理方法及系统 |
CN115372416A (zh) * | 2022-02-23 | 2022-11-22 | 重庆文理学院 | 一种感应式脉冲压缩磁声检测方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5465730A (en) * | 1991-07-26 | 1995-11-14 | British Technology Group Ltd. | Electrical impedance tomography |
CN104434099A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 利用线性双旋度方程进行磁热声成像的电阻率重建方法 |
CN104458818A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 一种基于线性泊松方程的磁热声成像电导率重建方法 |
CN104434100A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 一种磁热声成像的电阻率重建方法 |
CN104473639A (zh) * | 2014-12-14 | 2015-04-01 | 中国科学院电工研究所 | 一种基于最优化迭代算法的磁热声成像电阻率重建方法 |
-
2017
- 2017-02-20 CN CN201710089515.7A patent/CN107049315B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5465730A (en) * | 1991-07-26 | 1995-11-14 | British Technology Group Ltd. | Electrical impedance tomography |
CN104434099A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 利用线性双旋度方程进行磁热声成像的电阻率重建方法 |
CN104458818A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 一种基于线性泊松方程的磁热声成像电导率重建方法 |
CN104434100A (zh) * | 2014-12-14 | 2015-03-25 | 中国科学院电工研究所 | 一种磁热声成像的电阻率重建方法 |
CN104473639A (zh) * | 2014-12-14 | 2015-04-01 | 中国科学院电工研究所 | 一种基于最优化迭代算法的磁热声成像电阻率重建方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111887807A (zh) * | 2020-03-13 | 2020-11-06 | 中国石油大学(华东) | 一种纳秒脉冲电场诱导热声信号处理方法及系统 |
CN115372416A (zh) * | 2022-02-23 | 2022-11-22 | 重庆文理学院 | 一种感应式脉冲压缩磁声检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107049315B (zh) | 2020-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107064302B (zh) | 一种注入电流式热声成像电导率重建方法 | |
Adam et al. | Survey on medical imaging of electrical impedance tomography (EIT) by variable current pattern methods | |
Mariappan et al. | Magnetoacoustic tomography with magnetic induction: bioimepedance reconstruction through vector source imaging | |
Teniou et al. | A multimodal image reconstruction method using ultrasonic waves and electrical resistance tomography | |
Li et al. | Multi-excitation magnetoacoustic tomography with magnetic induction for bioimpedance imaging | |
CN105816156A (zh) | 一种结合热声效应信号的电导率磁声成像装置及成像方法 | |
CN104146708B (zh) | 一种基于声换能器接收特性的磁声图像重建方法及装置 | |
CN104473639B (zh) | 一种基于最优化迭代算法的磁热声成像电阻率重建方法 | |
CN104688224B (zh) | 一种应用于声学非均匀媒介磁声耦合成像重建方法 | |
Liang et al. | Nonstationary image reconstruction in ultrasonic transmission tomography using Kalman filter and dimension reduction | |
Mariappan et al. | B-scan based acoustic source reconstruction for magnetoacoustic tomography with magnetic induction (MAT-MI) | |
Liu et al. | Magnetoacoustic tomography with current injection | |
CN106885842B (zh) | 一种注入电流式热声成像电阻率重建方法 | |
CN104434099B (zh) | 利用线性双旋度方程进行磁热声成像的电阻率重建方法 | |
CN104434094A (zh) | 一种磁-热-声耦合成像的电导率图像重建方法 | |
CN107049315A (zh) | 基于最优化迭代方法的注入电流式热声电阻率图像重建方法 | |
CN104473640B (zh) | 一种磁热声成像的电导率重建方法 | |
CN106037638A (zh) | 一种去除热声效应影响的电导率磁声成像装置及成像方法 | |
CN104458818B (zh) | 一种基于线性泊松方程的磁热声成像电导率重建方法 | |
CN106580249B (zh) | 一种注入电流式热声成像方法 | |
CN106859600A (zh) | 基于建立标准值的声电成像方法 | |
CN104644170B (zh) | 基于声电效应的三维功率密度成像方法 | |
CN104434100B (zh) | 一种磁热声成像的电阻率重建方法 | |
Guillermin et al. | Quantitative non-linear ultrasonic imaging of targets with significant acoustic impedance contrast—An experimental study | |
CN107536608A (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 |