一种基于线性泊松方程的磁热声成像电导率重建方法
【专利摘要】一种基于线性泊松方程的磁热声成像电导率重建方法，激励线圈在导电物体上产生电磁热声信号，超声换能器接收电磁热声信号，超声信号处理、采集子系统采集和处理此信号，控制电路控制电流激励源、超声换能器和超声信号处理、采集子系统的同步。超声换能器电磁热声信号进行断层圆周扫描，获取每个断层圆周上的电磁超声信号，最后结合图像重建算法实现电导率图像的重建。电导率图像重建方法首先定义一个满足热声源、电导率、一次磁失位空间分量和标量电位空间分量的目标函数，并给定电导率的初值，在热声源分布已知的情况下，根据电流连续性定理，求解标量电位的空间分量，然后将标量电位的空间分量和磁失位空间分量代入目标函数，重建电导率的分布。
【专利说明】一种基于线性泊松方程的磁热声成像电导率重建方法
【技术领域】
[0001] 本发明涉及一种电导率图像的重建方法，特别涉及一种基于线性泊松方程的磁热 声成像电导率重建方法。
【背景技术】
[0002] 目前传统电阻抗成像技术的灵敏度和空间分辨率不高，主要因为电阻抗成像通常 采用频率较低的电磁波作为激励，由于波长远远大于成像体，导致电磁场探测对比度高，但 分辨率低。毋庸置疑，单一场都有其物理局限性。因此多物理场成像技术受到越来越多的 关注，即将一种物理场作用于生物组织，转换为另一种物理场进行检测，由一种物理场提供 分辨率，另一种物理场提供对比度，实现对比度和分辨率的同时提高。电磁场和超声相结合 的多物理场成像技术正是考虑到电磁场对人体组织电导率的高对比度和超声波探测的高 分辨率特性，成为人们的研宄热点，磁热声成像作为一种新型的多物理场成像技术在最近 一年受到重视。
[0003] 磁热声成像是由新加坡南洋理工大学在2013年首次提出的新型的电阻抗成像方 法，通过对导电物体施加低于20MHz的交变磁场，在导电物体内部产生感应电场，进而产生 焦耳热，激发热弹性的声信号，检测声信号进行成像。该方法是一种以交变磁场作为激励 源，基于生物组织内部焦耳热吸收率的差异，以超声作为信息载体的无损生物医学影像技 术。与微波热声成像技术相比，激励源的频率降低，可以深入到导电体的更深处，使磁热声 图像扩展到人体组织的深层。由测量的超声信号到电导率的重建分为两个过程，首先由测 量的超声信号重建热声源分布，然后利用热声源分布重建电导率分布，目前的相关文献和 专利只重建了热声源（S= 〇E2 ( 〇 )，这里E为电场强度的空间分量），而没有提及电导率 〇的重建。显然，电场强度E与电导率〇的分布有关，从热声源S中重建出电导率〇是非 常困难的。
【发明内容】
[0004] 本发明的目的是克服现有的磁热声成像方法无法给出电导率分布的不足，提出一 种基于线性泊松方程最优化迭代的磁热声成像电导率重建方法，本发明可以精确的重建导 电物体的电导率。
[0005] 本发明基于磁热声成像原理：利用激励线圈对导电物体施加MHz电流激励，在导 电物体内产生焦耳热，进而产生超声信号，利用超声换能器接收超声信号，对接收到的超声 信号借助电导率图像重建子系统重建导电物体的电导率图像。
[0006] 本发明磁热声成像导电率的重建方法主要包括五个步骤：第一步获取有效的磁热 声信号；第二步利用磁热声信号，采用时间反演法获取热声源分布；第三步对导电物体进 行空间离散，给出电导率的初值，根据电流连续性定理和磁热声分布，求解线性泊松方程获 取导电物体标量电位的空间分量；第四步设定泊松方程满足的目标函数，；第五步将标量电 位的空间分量和一次磁失位空间分量代入目标函数，通过最优化算法进行迭代，直到误差 满足终止条件，最终获取导电物体的电导率分布。
[0007] 具体步骤如下：
[0008] 第一步：获取磁热声信号
[0009] 首先由具有线性调频电流激励源产生线性调频电流激励信号，产生的电流激励 信号经激励线圈作用在导电物体上，导电物体由于感应电流的作用产生焦耳热，进而产生 电磁热声信号。本发明线性调频信号的调制频率范围在10MHz-30MHz，载波频率范围在 0. 5MHz-5MHz范围可选。利用超声换能器接收产生的电磁热声信号，然后通过超声信号处 理、采集子系统对信号进行采集并处理，控制电路实现对电流激励源、超声换能器和超声信 号处理、采集子系统的控制和同步，利用超声换能器对产生的电磁热声信号进行断层圆周 扫描，获取每个断层圆周上的电磁热声信号，将获得电磁热声信号进行存储。
[0010] 第二步：获取热声源分布
[0011] 根据电磁热声的产生原理，得出电磁热声成像的声压波动方程：
【权利要求】
1. 一种基于线性泊松方程的磁热声成像电导率重建方法，其特征在于：所述的磁热声 成像电导率重建方法包括以下步骤： 第一步：获取磁热声信号 由具有线性调频电流激励源（1)产生线性调频电流激励信号经激励线圈作用在导电 物体（4)上，导电物体（4)由于感应电流的作用产生焦耳热，进而产生电磁热声信号；利用 超声换能器（6)接收产生的电磁热声信号，然后通过超声信号处理、采集子系统（7)对信号 进行采集和处理；控制电路实现对电流激励源、超声换能器（6)和超声信号处理、采集子系 统（7)的控制和同步；利用超声换能器（6)对产生的电磁热声信号进行断层圆周扫描，获取 每个断层圆周上的电磁热声信号，将获得电磁热声信号进行存储； 第二步：获取热声源分布 根据电磁热声的产生原理，得出电磁热声成像的声压波动方程：
其中r为超声换能器与热声源之间的距离，p(r,t)是声压，cx为热声源在介质中的传 播声速，Cp为导电物体的比热容，β为导电物体的热膨胀系数，δ(t)是狄拉克函数，S(r) 是热声源分布，t是时间项； 根据声压波动方程，利用时间反演法得到热函数的重建公式：
其中R= |r' -r|，?=|，r'是超声换能器的位置，Sd是超声换能器所在的平面， K 是声压对时间的一阶导数，η是r'位置Sd的法线矢量； 选取导电物体的某一断层面z=Ztl，利用公式（2)获取断层面上的热声源S(x，y，Ztl); 导电物体上所有热声源分布通过断层数据S(X，y，Ztl)在z方向的插值得到，或者通过移动 超声换能器通过测量的方式获取； 第三步：获取导电物体标量电位空间分量 导电物体的热声源直接与导电物体的电场强度和电导率相关： S= 〇E2= 〇E·E(3) 其中E是磁场强度； 因生物组织的电导率较低，电场强度的空间分量E表示为： E ^-Vφ-A1 (4) 其中A1为一次磁矢位的空间分量，一次磁矢位与导电物体无关，是电流激励源在真空 中产生的磁矢位，Φ是标量电位的空间分量，▽是哈密顿算符； 根据电流连续性定理，有： ya{y￠+A1)= 〇 (5) 其中，σ为导电物体的电导率，对导电物体进行空间离散，给出导电物体电导率的初值[σ]°，将[〇 ]°和一次磁矢位的空间分量A1R入公式（5)，利用线性有限元方法求解得到标 量电位的空间分量[φ]1; 第四步：定义目标函数 由公式⑶和公式⑷得到：
定义目标函数：
第五步：最优化迭代求解导电物体电导率： 将标量电位的空间分量[Φ]1和一次磁矢位的空间分量Ai代入目标函数，通过最优化 算法进行迭代，直到误差满足终止条件，所得的电导率即为最终导电物体的最终电导率。
【文档编号】G01R27/02GK104458818SQ201410773989
【公开日】2015年3月25日 申请日期:2014年12月14日 优先权日:2014年12月14日
【发明者】刘国强, 夏慧, 夏正武, 李士强, 杨延菊, 刘宇 申请人:中国科学院电工研究所