专利摘要

一种注入电流式热声成像电阻率重建方法，基于电流注入方式的热声成像原理。通过电极向成像目标体注入脉冲电流，在成像目标体中产生焦耳热，引起热膨胀，产生超声信号，用超声换能器接收超声信号，对接收到的超声信号进行处理和采集，采用电阻率图像重建算法获取目标体的电阻率图像。具体步骤为：1、首先获取注入电流式热声信号；2、利用获取的注入电流式热声信号重建目标体热声源；3、利用热声源，采用非线性有限元求解方法重建矢量电位；4、利用重建的矢量电位重建电阻率。

权利要求

1.一种注入电流式热声成像电阻率重建方法，其特征在于，所述的注入电流式热声电导率图像重建方法基于注入电流式的热声成像原理，通过电极向成像目标体注入脉冲电流，在成像目标体中产生焦耳热，引起热膨胀，产生超声信号，用超声换能器接收超声信号，对接收到的超声信号进行处理和采集，采用电阻率图像重建算法获取目标体的电阻率图像；

热声源的获取方法为：已知热声成像的声压波动方程：

${▿}^{2}p(r,t)-\frac{1}{{c_s}^2}\frac{{\partial }^2}{\partial t^2}p(r,t)=-\frac{\beta }{C_P}S\left(r\right){\delta }^{\prime }\left(t\right)---\left(1\right)$

其中r为超声换能器位置，p(r,t)是声压，c_s为介质中的声速，C_P为目标体(3)的比热容，β为目标体(3)的热膨胀系数，δ(t)是狄拉克函数，S(r)为热声源分布，t为时间， 为拉普拉斯算符；

选取目标体的某一断层面z＝z₁，超声换能器在此断层面上进行圆周扫描，采集超声信号，利用方程(1)，求解z＝z₁断层面上的热声源分布S(x,y,z₁)，移动超声换能器在不同断层面上进行扫描检测，求解不同断层面上的热声源分布，目标体(3)上整体热声源S可由分层检测计算或z方向上插值得到；

热声源S同时是电阻率和电流密度的函数表示为：

S＝ρJ²＝ρJ·J (3)

其中，ρ为目标体(3)的电阻率，J为目标体(3)内电流密度分布；

所述的注入电流式热声成像电阻率重建方法包括以下步骤：

第一步：获取注入电流式热声信号

激励源(1)通过注入电极A(2)和注入电极B(4)对目标体(3)注入脉冲电流，目标体(3)在电流作用下产生焦耳热，进而产生热膨胀，激发超声信号，超声信号通过耦合剂耦合到超声换能器(5)，超声换能器(5)接收到信号后通过检测系统(6)进行信号的放大、滤波、采集和存储，在上位机(7)中实现图像重建，超声换能器(5)在控制器(8)控制下对目标体进行扫描检测；

第二步：获取目标体热声源

根据第一步采用一对注入电极A(2)、注入电极B(4)，向目标体注入电流后，引发热膨胀，激发声信号，检测到声压信号，采用声压信号重建目标体上整体热声源；

第三步：求解矢量电位

根据电流连续性定理，电流的散度为零， 引入矢量电位，则有：

$J=▿\times T---\left(4\right)$

其中，T为矢量电位， 为矢量电位的旋度， 为哈密顿算符；

利用欧姆定律，J＝σE＝E/ρ，则有：

$▿\times (\rho ▿\times T)=0---\left(5\right)$

由式(3)和式(4)可得：

$\rho =\frac{S}{▿\times T·▿\times T}---\left(6\right)$

将式(6)代入式(5)中，得到：

$▿\times (\frac{S}{▿\times T·▿\times T}▿\times T)=0---\left(7\right)$

满足的边界条件为：

$\left\{\begin{array}{c}▿\times T·n{|}_{{\Gamma }_{A,B}}=\pm \frac{I}{A_0}\\ ▿\times T·n{|}_{{\Gamma }_g}=0\end{array}\right.---\left(8\right)$

其中，Γ_A,B为注入电极位置，Γ_g为除去电极之外的目标体边界，A₀为电极与目标体接触面积，I为注入电流，n表示区域的外法向单位矢量；

将热声源S代入公式(7)，结合边界条件(8)，进行有限元法求解，即可重建得到矢量电位T；

第四步：求解电阻率

将矢量电位T代入公式(6)，即可重建电阻率ρ。

说明书

技术领域

本发明涉及一种电阻率图像重建方法，特别涉及一种注入电流式电阻率图像重建方法。

背景技术

由于激励频率的限制，传统电阻抗成像技术的灵敏度和空间分辨率不高。单一场都有其物理局限性，多物理场成像由一种物理场提供分辨率，另一种物理场提供对比度，实现对比度和分辨率的同时提高。电磁场和超声相结合的多物理场成像技术正是考虑到电磁场对人体组织电阻率的高对比度和超声波探测的高分辨率特性，成为人们研究的热点，磁热声成像正是一种新兴的多物理场成像技术。

磁热声成像是由新加坡南洋理工大学的Feng在2013年首次提出的新型的电阻抗成像方法，其原理为：通过对导电成像体施加MHz量级的交变磁场，在目标体内部产生感应电场，进而产生焦耳热，激发热弹性的超声信号，检测超声信号进行成像。与微波热声成像相比，允许更低的功率进行高效的成像，并且具有便携式成像的潜力，同时，激励源的频率降低，使得磁场穿透组织更深，也避免的辐射。

作为一种新型的多物理场成像方法，2013年Feng利用金属铜仿体，检测到磁热声信号，并得到铜仿体的热声图像，并未进行电阻率图像重建，且生物组织不同于金属铜仿体，磁场和感应电流作用产生的洛伦兹力较弱。专利“一种磁热声成像的电阻率重建方法”(201410771496.2)公布了一种基于磁热声效应的电阻率重建方法，其在热函数基础上提出了电阻率重建方法，依然为采用线圈激励方式。同时，外加激励线圈产生时变磁场在目标体中产生二次磁场和感应电流，磁场和感应电流作用同时产生洛伦兹力和焦耳热，即磁声效应和磁热声效应共存，如何区分磁声效应和磁热声效应是仍需解决的问题。

基于此，采用注入电流式热声成像方法，可以避开磁声效应和磁热声效应共存问题，注入电流式热声成像与磁热声成像在激励方式、从热函数到电阻率的重建方法上均不相同。

发明内容

本发明的目的是克服现有磁热声成像存在的缺点，提出一种基于注入电流式热声成像的电阻率重建方法。本发明利用热声源重建电阻率，避免了热声成像过程中磁声效应的干扰，同时采用注入电流式激励，可增强热声效应，可实现目标体电阻率图像的精确重建。

注入电流式热声成像原理为：通过注入电极向成像目标体注入电流，在成像目标体中产生焦耳热，引起热膨胀，产生超声信号，采用超声换能器检测超声信号，根据检测的超声信号重建热声源和电阻率。

本发明注入电流式热声成像的电阻率图像重建包括四个步骤：1、首先获取注入电流式热声信号，即检测超声信号；2、利用获取的注入电流式热声信号重建目标体热声源；3、利用热声源，采用非线性有限元求解方法重建矢量电位；4、利用重建的矢量电位重建电阻率。

电阻率图像重建过程具体描述如下：

第一步：获取注入电流式热声信号

激励源通过注入电极A和注入电极B对目标体注入脉冲电流，目标体在电流作用下产生焦耳热，进而产生热膨胀，激发超声信号，超声信号通过耦合剂耦合到超声换能器，超声换能器接收到信号后通过检测系统进行信号的放大、滤波、采集和存储，超声换能器在控制器控制下对目标体进行扫描检测；

第二步：获取目标体热声源

已知热声成像的声压波动方程：

其中r为超声换能器位置，p(r,t)是声压，c_s为介质中的声速，C_P为目标体(3)的比热容，β为目标体(3)的热膨胀系数，δ(t)是狄拉克函数，S(r)为热声源分布，t为时间， 为拉普拉斯算符；

热声源分布的时间反演法重建公式为：

其中,R为标量，R＝|r′-r|，R为矢量， e_R为单位矢量，r'是超声探头的位置，r为热声源位置，S_d是超声探头所在的平面，p′是声压对时间的一阶导数，n是r′位置S_d的法线矢量，β为目标体的热膨胀系数；

选取目标体的某一断层面z＝z₁，超声换能器在此断层面上进行圆周扫描，采集超声信号，利用方程(1)和公式(2)求解z＝z₁断层面上的热声源分布S(x,y,z₁)，移动超声换能器在不同断层面上进行扫描检测，求解不同断层面上的热声源分布，目标体(3)上整体热声源S可由分层检测计算或z方向上插值得到；

热声源S同时是电阻率和电流密度的函数，可以表示为：

S＝ρJ²＝ρJ·J (3)

其中，ρ为目标体的电阻率，J为目标体内电流密度分布；

第三步：求解矢量电位

根据电流连续性定理，电流的散度为零， 引入矢量电位，则有：

其中，T为矢量电位， 为矢量电位的旋度， 为哈密顿算符；

利用欧姆定律，J＝σE＝E/ρ，则有：

其中，ρ为目标体的电阻率， 为矢量电位的旋度。

由式(3)和式(4)可得：

将式(6)代入式(5)中，得到：

满足的边界条件为：

其中，Γ_A,B为注入电极位置，Γ_g为除去电极之外的目标体边界，A₀为电极与目标体接触面积，I为注入电流，n表示区域的外法向单位矢量；

将热声源S代入公式(7)，结合边界条件(8)，进行有限元法求解，即可重建得到矢量电位T；

第四步：求解电阻率

将矢量电位T代入公式(6)，即可重建电阻率ρ。

附图说明

图1本发明重建方法所涉及的注入电流式热声信号获取装置示意图；

图中：1激励源、2注入电极A、3目标体、4注入电极B、5超声换能器、6检测系统、7上位机、8控制器、9水槽。

具体实施方式

以下结合附图和具体实施方式对本发明做进一步说明。

激励源1连接注入电极A2和注入电极B4，向目标体3中注入电流信号。超声换能器5与目标体3之间通过耦合剂耦合，超声换能器5的输出端连接检测系统6的输入端，检测系统6的输出端连接上位机7，在上位机7中实现图像重建。控制器8实现对超声换能器5的旋转扫描运动控制。目标体3、注入电极A2、注入电极B4和超声换能器5置于水槽9中。

注入电流式热声成像原理为：通过注入电极向成像目标体注入电流，在成像目标体中产生焦耳热，引起热膨胀，产生超声信号，采用超声换能器进行检测，根据检测的超声信号，重建热声源和电阻率。

本发明注入电流式热声成像的电阻率图像重建包括四个步骤：1、首先获取注入电流式热声信号，即检测超声信号；2、利用获取的注入电流式热声信号重建目标体热声源；3、利用热声源，采用非线性有限元求解方法重建矢量电位；4、利用重建的矢量电位重建电阻率。

图像重建过程具体描述如下：

第一步：获取注入电流式热声信号

激励源1通过注入电极A2和注入电极B4对目标体3注入脉冲电流，目标体3在电流作用下产生焦耳热，进而产生热膨胀，激发超声信号，信号通过耦合剂耦合到超声换能器5，超声换能器5接收到信号后通过检测系统6进行信号的放大、滤波、采集和存储，超声换能器5在控制器8的控制下对目标体进行扫描检测；

第二步：获取目标体热声源

已知热声成像的声压波动方程：

其中r为超声换能器位置，p(r,t)是声压，c_s为介质中的声速，C_P为目标体(3)的比热容，β为目标体(3)的热膨胀系数，δ(t)是狄拉克函数，S(r)为热声源分布，t为时间， 为拉普拉斯算符；

热声源分布的时间反演法重建公式为：

其中,R为标量，R＝|r′-r|，R为矢量， e_R为单位矢量，r'是超声探头的位置，r为热声源位置，S_d是超声探头所在的平面，p′是声压对时间的一阶导数，n是r′位置S_d的法线矢量，β为目标体的热膨胀系数；

选取目标体的某一断层面z＝z₁，超声换能器在此断层面上进行圆周扫描，采集超声信号，利用方程(1)、(2)，求解z＝z₁断层面上的热声源分布S(x,y,z₁)，移动超声换能器在不同断层面上进行扫描检测，求解不同断层面上的热声源分布，目标体3上整体热声源S可由分层检测计算或z方向上插值得到；

热声源S同时是电阻率和电流密度的函数，可以表示为：

S＝ρJ²＝ρJ·J (3)

其中，ρ为目标体3的电阻率，J为目标体3内电流密度分布；

第三步：求解矢量电位

根据电流连续性定理，电流的散度为零， 引入矢量电位，则有：

其中，T为矢量电位， 为矢量电位的旋度， 为哈密顿算符；

利用欧姆定律，J＝σE＝E/ρ，则有：

其中，ρ为目标体的电阻率， 为矢量电位的旋度。

由式(3)和式(4)可得：

将式(6)代入式(5)中，得到：

满足的边界条件为：

其中，Γ_A,B为注入电极位置，Γ_g为除去电极之外的目标体边界，A₀为电极与目标体接触面积，I为注入电流，n表示区域的外法向单位矢量；

将热声源S代入公式(7)，结合边界条件(8)，进行有限元法求解，即可重建得到矢量电位T；

第四步：求解电阻率

将矢量电位T代入公式(6)，即可重建电阻率ρ。

一种注入电流式热声成像电阻率重建方法专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。