一种应用于声学非均匀媒介磁声耦合成像重建方法
【技术领域】
[0001] 本发明涉及一种磁声耦合成像重建算法。特别是涉及一种利用换能器所检测的磁 声信号点的相关性求解声速分布，进而重建原始声源分布的应用于声学非均匀媒介磁声耦 合成像重建方法。
【背景技术】
[0002] 磁声耦合成像是一种结合电阻抗成像和超声成像的技术的新型生物组织电特性 成像技术。其原理是：将被测生物组织置于稳恒磁场中，对生物组织给予稳恒磁场方向平行 的脉冲磁场，进而在其内部产生感应电流（或者通过电极直接往生物组织内部注入脉冲电 流）；被测组织内部的注入电流在稳恒磁场的作用下受洛伦兹力影响产生与激励信号同频 率的位移，进而产生高频振动，向外界辐射超声波；包含物体内部的声信号在生物组织内部 通过一系列复杂的反射、折射等变化传递到体表，被换能器所接收并用于重建内部电导率 信息。
[0003] 由磁声耦合成像的原理可知，一个电流分布为? (变量上端的波浪线指代变量包 含时间量）的媒介置于一个静态磁场Btl中，可以得到声压的波动方程如下：
【主权项】
1. 一种应用于声学非均匀媒介磁声耦合成像重建方法，其特征在于，是在信号采集过 程中使用线阵换能器阵列采集相互对称的声压数据，利用对称位置换能器采集数据之间的 相关性求解声穿越时间，结合代数迭代重建算法重建样本声速分布，最后利用所求声速分 布求解样本初始声源分布，具体包括如下阶段： 1) 在扫描过程中采用平面线阵换能器阵列采集样本内部发出的磁声信号； 2) 求解换能器之间的声穿越时间，包括： (1)设定换能器检测的声压满足的条件；（2)进行声压补偿；（3)求解对称位置换能器 所检测声压的相关系数；(4)求解两个轴线相互重叠的换能器聚焦区轴线部位的声穿越时 间；（5)求解得到整个扫描过程的声穿越时间； 3) 求解被测样本的声速分布图，包括： (1)建立对称位置换能器连线的声速平均矩阵；（2)修正平均声速矩阵；（3)划分各单 元网格；(4)求解声穿越时间；（5)对结果再次修正之后得到每个网格实际声速； 4) 求解原始声源分布，包括： (1)声速重排及声速空间扩展；（2)反转声源求解及声压转置；（3)声压数据重采样； (4)求解原始声源；（5)提取样本所在区域声源分布。
2. 根据权利要求1所述的一种应用于声学非均匀媒介磁声耦合成像重建方法，其特征 在于，阶段1)在扫描过程中每一个换能器都要满足如下条件： (1) 换能器扫描过中完成一次扫描后进行旋转，旋转次数要满足2M-1次； (2) 每一个角度的聚焦区域移动次数为N，移动间距Ad保持不变，且满足NXAd多R， 其中R为能容纳样本的最小正方形区域的边长，M和N是大于等于1的整数； (3) 扫描过程中，换能器旋转半径r要满足： r > sfl.R； (4) 处于对称位置的换能器，在聚焦区轴线重叠的区域内每一个点满足如下条件： Q+L产2r 其中U为所述点到左边换能器的距离，L2为所述点到右边换能器的距离，r为换能器 扫描过程中的旋转半径。
3. 根据权利要求1所述的一种应用于声学非均匀媒介磁声耦合成像重建方法，其特征 在于，阶段2)所述的求解换能器之间的声穿越时间，包括： (1) 设定换能器检测的声压满足如下条件：
其中A为换能器聚焦区轴线的延长线与换能器的交点，t为换能器接收的时间轴，； 为换能器所接收的声压，J为感应电流密度，B为稳恒磁场强度，r'为样本上某一点的位置，T(r0, -r0)为重叠的聚焦区轴线间的声传播时间。 (2) 进行声压补偿 定义换能器接收信号的修正值S(r(l，t)满足：
(3) 求解对称位置换能器所检测的声压的相关系数 相对于位置A的换能器接收的声压修正后为S(r&t)，设定与r(l位置对称的信号采集 点位置为-IV故所述信号采集点-r(l接收的声压修正后为S(_r&t)，考虑声的传播特性，位 置A和位置-r啲声压信
号修正值在时间和幅度上均具有高度相关，定义心。（〇满足： 为聚焦区域轴线相互重叠的一组换能器所接收信号的相关系数； (4) 求解两个换能器聚焦区轴线相互重叠部位的声穿越时间 由相关性原理及卷积定理可以得到，当t=TCt^，-ij时，式TCt^，-ij可以得到最大 值，故由下式得到两个换能器聚焦区轴线相互重叠部位的声穿越时间
(5) 求解得到整个扫描过程的声穿越时间 对阶段1)扫描所得到的MXN对重叠的聚焦区轴线位置，依次重复步骤（1)~步骤 (4)，得到MXN组穿越时间。
4.根据权利要求1所述的一种应用于声学非均匀媒介磁声耦合成像重建方法，其特征 在于，阶段3)所述的求解被测样本的声速分布图，具体包括： (1) 建立对称位置换能器连线的传播时间矩阵 对于扫描过程中换能器动态聚焦区的位置分布，定义传播时间矩阵TNM的任何一个元 素要满足如下公式： Ti=Krj,-rj) ；1 ^i^NXM 其中b-ri代表第i对聚焦区轴线延长线与换能器交点所处的位置； (2) 修正传播时间矩阵TNM，得到声穿越时间T'm 设定耦合剂的所处区域作为背景，背景声速为C(l，对于扫描过程中换能器的动态聚焦 区位置分布，得到对称位置换能器之间的距离等同于旋转直径，故排除背景声速C(l影响后 的声穿越时间T'm的元素满足：
Am为元素全为1的MXN矩阵； 由上式得到声穿越时间T'M，T'_是被测样本被换能器聚焦区轴线穿过部位的声穿越 时间再减去相同距离耦合剂的声穿越时间； (3) 划分各单元网格 将一个感兴趣的包含被测样本的扫描空间的区域分成nXn的网格，网格大小为AS， 每个网格的声穿越时间为tj;l彡j彡n2,网格大小As应当满足如下： 上式中
为磁声耦合成像自由空间中的最大声波频率，满足： fmax^ 2f in fin为磁声耦合成像的激励信号频率，cmin为磁声耦合系统中的最小声速； (4) 求解声穿越时间 根据步骤（3)的网格，按如下公式进行迭代重建，求解每一个网格的声穿越时间值tj:
k为迭代次数，X为松弛因子，a为权重因子，当网格m对声穿越时间h有贡献时取值 为1，否则为〇; (5) 对结果再次修正之后得到每个网格实际声速 步骤（4)求解得到的声穿越时间为去除背景声速之后的声穿越时间，对结果再次修正 之后得到每个网格实际声速Cj;l<j<n2满足： Av, . 9 &=-1 <./ <tr〇 h
5.根据权利要求4所述的一种应用于声学非均匀媒介磁声耦合成像重建方法，其特征 在于，步骤（4)所述的迭代重建步骤包括： a) 给第j个网格的声穿越时间武初值： =t(p<j<n2 式中，初值选取依据人体组织的平均声速和耦合剂的背景声速，满足： (0)=A5_A5Cpc0 其中匕为人体组织的平均声速； b) 计算第i对换能器间的声穿越时间修正后的估计值：
c) 计算误差： ^=T；~T； d) 计算第j个未知量的修正值：
e) 对tj的值进行修正： T(k+l) = T(k) +^Cl f) k=k+1 ;重复步骤b)到f)，直到完成所有声穿越时间方程，即完成一轮迭代； g)以上一轮迭代作为初值，重复b)到f)进行新一轮迭代，直到得到符合收敛要求的结 果。
6.根据权利要求1所述的一种应用于声学非均匀媒介磁声耦合成像重建方法，其特征 在于，阶段4)所述的求解原始声源分布，具体包括： (1) 声速重排及声速空间扩展 利用阶段3)所求得到的每个网格实际声速Cfl<j<n2,进行声速重排，得到nXn的声速分布矩阵;将所述声速分布矩阵向外扩充为一个mXm;m>n的自由空间，声 速分布为CmXm;扩充部分的网格的声速为c^ (2) 反转声源求解及声压转置 对扫描过程中换能器接收到的声压信号为P(r，t)，其中r表示空间位置；依据时间反 演理论，将采集的声压信号反转作为声源，定义用于反演的声源项为q(r，t)满足： q(r,t) =p(r,Td-t) ； 1 ^i^NXM Td表示检测过程中，空间每一个点的检测时间，故得基于时间反演的波动方程如下：
c(r)为声速分布CmXm圆柱坐标形式，P(l(r，Td)为需要求解的原始声源分布； (3) 声压数据重采样 定义时间间隔
其中，为磁声耦合成像实验中的最大声速， 依据时间间隔At对反演的声源项q(r，t)进行重采样，采样结果用反演声源Q(t)表 示：
(4) 时间和空间交替迭代求解原始声源 依据声波波动方程，结合磁声耦合成像声源产生机制，得到不均匀媒介中的声压和振 动速度的表达式，如下：
其中，vx，vy表示质点二维平面x，y方向的振动速度，表示声压，P表示物质密度，c表示声速，t表示波动持续时间。采用时域有限差分的方法，在时间和空间上交替离散上述 表达式，利用反演声源Q求解下一个时刻的质点振动速度，并利用质点振动速度求解声压， 依此交替； (5)提取样本所在区域声源分布 依据自由空间大小，将原始声源分布Ptl(r，Td)，以自由空间中心为基准点，半径为r之 外的所有区域声压置零，所得到的声压值P〇即为初始声源分布。
7.根据权利要求6所述的一种应用于声学非均匀媒介磁声耦合成像重建方法，其特征 在于，步骤（4)所述的不均匀媒介中的声压和振动速度的表达式具体求解过程如下： a) 确定所述表达式求解的初始条件，所述的初始条件包括：磁通密度B，物质密度P， 声速媒介声速分布CmXm和采集的声压信号Q(1At); b) 离散化不均匀媒介中的声压和振动速度的表达式
将过程a)中确定的初始条件代入过程b)中所给出的声压方程中，求解得到初始时刻ta= 0的声压p d) 将过程c)所求ta时刻的声压P(l代入过程b)中所给出的振动速度方程中，求解 时刻的振速Vx,vy; e) 依据过程d)所求tb时刻的振速vx，vy，将换能器检测的反演声源Q代入过程b)中 所给出的声压方程中，求解^ =6 11时刻的声压P，存储数据； f) 如果ta=Td，转到过程d)，否则结束； 当时ta=Td时，依据上述步骤所求解得到的声压分布p(r，Td)，即包含原始声源分布。
【专利摘要】一种应用于声学非均匀媒介磁声耦合成像重建方法，是在信号采集过程中使用线阵换能器阵列采集相互对称的声压数据，利用对称位置换能器采集数据之间的相关性求解声穿越时间，结合代数迭代重建算法重建样本声速分布，最后利用所求声速分布求解样本初始声源分布。本发明的一种应用于声学非均匀媒介磁声耦合成像重建方法，可以有效消除样本的声学特性不均匀带来的重建误差，同时也能够通过信号的相关性得到样本的声速分布，为声源的精确重建提供了必要的条件。
【IPC分类】A61B5-05
【公开号】CN104688224
【申请号】CN201510149821
【发明人】刘志朋, 张伟, 殷涛, 马任
【申请人】中国医学科学院生物医学工程研究所
【公开日】2015年6月10日
【申请日】2015年3月31日