一种磁-热-声耦合成像的电导率图像重建方法
【专利摘要】一种磁-热-声耦合成像的电导率图像重建方法，首先根据接收到的电磁超声信号利用时间反演法获取导电物体的热声源分布，然后根据热声源分布再重建得到导电物体的电导率分布。具体过程为，首先利用时间反演法获取某一段层面上的热声源，然后利用插值法获取导电物体的热声源，接着根据激励电流，利用毕奥-萨伐尔定律计算得到一次磁矢位的空间分量A1，对导电物体进行空间离散，给出电导率的初值[σ]0，然后根据电流连续性定理，已知[σ]0和A1求解标量电位空间分量[φ]1，再接着将标量电位的空间分量[φ]1代入热声源与电导率满足的关系式，得到更新后的电导率[σ]1；利用[σ]1替换[σ]0，重复以上描述的过程，直到电导率相对误差满足ε＝||([σ]1-[σ]0)/[σ]0||2≤ε0，即可停止迭代过程。
【专利说明】一种磁-热-声耦合成像的电导率图像重建方法
【技术领域】
[0001] 本发明涉及一种电导率图像的重建方法，特别涉及一种磁-热-声耦合成像的电 导率图像重建方法。
【背景技术】
[0002] 目前传统电阻抗成像技术的灵敏度和空间分辨率不高，主要因为电阻抗成像通常 采用频率较低的电磁波作为激励。由于波长远远大于成像体，导致电磁场探测对比度高，但 分辨率低。毋庸置疑，单一场都有其物理局限性。因此多物理场成像技术受到越来越多的 关注，即将一种物理场作用于生物组织，转换为另一种物理场进行检测，由一种物理场提供 分辨率，另一种物理场提供对比度，实现对比度和分辨率的同时提高。电磁场和超声相结合 的多物理场成像技术正是考虑到电磁场对人体组织电导率的高对比度和超声波探测的高 分辨率特性，成为人们的研宄热点，磁热声成像作为一种新型的多物理场成像技术在最近 一年受到重视。
[0003] 磁热声成像是由新加坡南洋理工大学在2013年首次提出的新型的电阻抗成像方 法，通过对导电物体施加低于20MHz的交变磁场，在导电物体内部产生感应电场，进而产生 焦耳热，激发热弹性的声信号，检测声信号进行成像。该方法是一种以交变磁场作为激励 源，基于生物组织内部焦耳热吸收率的差异，以超声作为信息载体的无损生物医学影像技 术。与微波热声成像技术相比，激励源的频率降低，可以深入到导电体的更深处，使磁热声 图像扩展到人体组织的深层。由测量的超声信号到电导率的重建分为两个过程，首先由测 量的超声信号重建热声源分布，然后利用热声源分布重建电导率分布，目前的相关文献和 专利只重建了热声源（S = 〇 E2 ( 〇 )，这里E为电场强度的空间分量），而没有提及电导率 〇的重建。显然，电场强度E与电导率σ的分布有关，从热声源S中重建出电导率σ是非 常困难的。
【发明内容】
[0004] 本发明的目的是克服现有的磁热声成像方法无法给出电导率分布的不足，提出一 种基于线性泊松方程直接迭代的磁-热-声耦合成像的电导率图像重建方法。本发明可以 精确的重建导电物体的电导率。
[0005] 本发明基于磁热声成像的原理：利用激励线圈对导电物体施加 MHz电流激励，导 电物体在电磁激励下产生焦耳热，进而产生超声信号，利用超声换能器接收超声信号，并对 接收到的超声信号利用电导率图像重建算法进行图像重建，即可获取导电物体的电导率图 像分布。
[0006] 本发明磁-热-声耦合成像的电导率图像重建方法包括五个步骤：第一步首先获 取导电物体的磁热声信号；第二步根据声压波动方程获取导电物体的热声源分布；第三步 求解标量电位空间分量；第四步利用求解的标量电位空间分量计算新的电导率；第五步设 定的电导率与求解的电导率进行相对误差计算，如果满足要求则新的电导率即为求解的电 导率，否则把新的电导率设定为电导率初值，重新进行步骤三和步骤四，直到设定电导率与 求解电导率之间的相对误差满足要求，则停止迭代，得到最终的导电物体的电导率。
[0007] 具体步骤如下：
[0008] 第一步：获取导电物体的磁热声信号
[0009] 为了获取高分辨率和高信噪比的磁热声信号，本发明采用准连续脉冲激励方法对 导电物体施加电磁激励。所述的准连续脉冲激励的脉冲宽度在2ns-0. Ius可调，连续脉冲 个数可调。电磁激励源产生的准连续脉冲激励信号通过激励线圈作用到导电物体上，激励 线圈与导电物体通过空气隔离，导电物体在感应电流作用下产生焦耳热，进而产生电磁热 声信号，电磁热声信号通过耦合剂耦合到超声换能器，超声换能器接收到电磁超声信号后 通过数据处理采集电路进行放大、滤波等处理，并存储。
[0010] 第二步：获取热声源分布
[0011] 已知磁热声成像的声压波动方程：
【权利要求】
1. 一种磁-热-声耦合成像的电导率图像重建方法，其特征在于：所述的磁-热-声 耦合成像的电导率图像重建方法包括以下步骤： 第一步：获取导电物体磁热声信号 采用准连续脉冲激励方法对导电物体施加电磁激励；所述的准连续脉冲激励的脉冲宽 度在2ns-0.Ius可调，连续脉冲个数可调，电磁激励源产生的准连续脉冲激励信号通过激 励线圈（1)作用到导电物体（3)上，激励线圈（1)与导电物体（3)通过空气隔离，也即，激励 线圈⑴位于水槽⑵的耦合剂之外，导电物体⑶和超声换能器⑷位于水槽的耦合剂之 内，导电物体（3)位于激励线圈（1)产生的激励磁场之中；导电物体（3)在感应电流作用下 产生焦耳热，进而产生电磁热声信号，电磁热声信号通过耦合剂耦合到超声换能器（4)，超 声换能器（4)接收到电磁超声信号后通过数据处理采集电路进行放大、滤波处理，并存储； 第二步：获取热声源分布 已知磁热声成像的声压波动方程：
其中r为热声源位置坐标，p(r，t)是声压，Cs为热声源在介质中的传播声速，Cp为导电 物体（3)的比热容，β为导电物体（3)的热膨胀系数，S(r)为热声源分布，δ(t)是狄拉克 函数，t是时间项；
选取导电物体的某一断层面z=Ztl，根据声压波动方程求解该断层面上的热声源分布S(X，y，Ztl)，并求解导电物体上的所有断层面上的热声源分布S(X，y，z); 热声源同时是电导率和电场强度的函数，因此也表示为： S= 〇E2= 〇E·E (2) 其中，σ为导电物体的电导率，E为激励电流在导电物体内产生的电场强度； 第三步：求解标量电位空间分量 根据电流连续性定理： ▽·σ(▽沴+ 4) =O (3) 其中A1为一次磁矢位的空间分量，一次磁矢位与导电物体无关，是激励源在真空中产 生的磁矢位，Φ是标量电位的空间分量，▽为哈密顿算符； 对导电物体进行空间离散，给出电导率的初值[σ]°，同时根据激励电流大小，利用毕 奥-萨伐尔定律计算得到一次磁矢位的空间分量A1，将[。]°和A入公式（3)，结合电绝 缘边界条件，进行线性有限元求解，求解得到标量电位空间分量[Φ]S 第四步：求解电导率 有公式（2)可以得出：
将热声源S(x，y，z)和标量电位的空间分量[Φ]1代入公式（4)，得到更新后的电导率 [σ]1; 第五步：迭代求解电导率 根据第四步的求解的电导率，计算是否满足ε=II([(^^[0]°)/^]°! |2彡ε。 其中^为给定的最大相对误差，ε为相对误差，[。]°给定的电导率初值；如果满足 则电导率[ο]1即为所求解的电导率；如果不满足要求，利用[〇] 1替换[〇] °，重复进行第 三步和第四步，直到满足要求为止。
【文档编号】A61B5/05GK104434094SQ201410767707
【公开日】2015年3月25日 申请日期:2014年12月14日 优先权日:2014年12月14日
【发明者】夏慧, 刘国强, 李士强, 夏正武, 杨延菊, 刘宇 申请人:中国科学院电工研究所