CN102894974B - 一种磁声电成像系统及成像方法 - Google Patents

Abstract

一种磁声电成像系统，主要包括超声驱动激励源、超声探头阵列、控制系统、磁体系统、电极对或线圈、信号检测处理系统。通过控制系统使超声探头阵列处于发射或测量两种模式。磁声电成像系统的成像方法，其特征在于采用声辐射力在生物组织内部产生聚焦区，并采用超声回波技术测量聚焦区位置和质点振动速度。由于聚焦区面积小，使得聚焦区内静磁场、质点振动速度和假想过程的电流密度可以近似为均匀值，从而建立起测量电压与假想过程的电流密度的线性关系。

Description

一种磁声电成像系统及成像方法

技术领域

本发明涉及一种医学成像方法及装置，特别涉及一种磁声电成像方法及装置。

背景技术

传统的电阻抗成像的灵敏度和空间分辨率不高，为解决此问题，人们提出了各种新型成像方法。磁声电成像就是一种新型的具有良好应用前景的医学成像方法。其成像原理是将一束超声波注入成像体，成像体中的局部离子随超声波的传播而振动，振动的离子在静磁场的作用下受到洛伦兹力作用而引起电荷分离，进而在成像体内形成局部的电场，通过贴放在成像体上的接收电极或与成像体非接触的接收线圈检测电信号重构电阻抗图像。

1997年，Han Wen等人提出霍尔效应成像的概念，并给出一维模型，利用普通超声探头和薄铜片搭建的简单实验设备。在2007年，Y.Xu、S Haider等人，在Han Wen概念提出的基础上，提出基于互易定理的磁声电成像，实验仍旧采用一维铜片样本，利用电极对进行测量，并根据互易定理推导出测量电压与声场和电磁场之间简单关系式，在系统的配置方面没有提及。

上述成像方法在提高灵敏度和空间分辨率方面具有潜在的优势，但是目前尚停留在理论研究和离体实验系统的研究阶段。采用超声激励生物组织，成像的横向分辨率取决于焦斑的大小，而轴向分辨率取决于超声脉冲的脉宽大小。生物组织内部的质点振动速度无法作为已知量，只能假定生物组织声学特性均匀，并采用声场正问题数值模拟的方法获得生物体内不同位置的质点振动速度，这与实际情况有所偏离。

发明内容

本发明的目的是克服现有技术的如上所述的缺点，提出的一种新的磁声电成像装置和方法。

本发明的主要原理为：将两束频率接近的超声波聚焦注入生物组织，在生物组织一定深度内产生声辐射力，形成焦斑大小的聚焦区，聚焦区中局部离子随声波的传播而振动，振动频率为两束超声波的差频。振动的离子在有静磁场的情况下，会受到洛仑兹力的作用，引起电荷分离，产生库仑电场，利用贴放在生物体表面的电极对或非接触的线圈测量电压信号，并通过超声检测装置检测聚焦区位置和质点振动速度。这样，聚焦区位置和质点振动速度以及静磁场均为已知参数，避免了现有只能假定生物组织声学特性均匀的数值模拟方法的不足。进而利用互易定理建立测量电压信号与聚焦区电导率的非线性关系，重构电阻抗图像。本发明所述的成像方法，成像的横向分辨率和轴向分辨率均取决于焦斑的大小。

本发明磁声电成像系统主要包括超声驱动激励源、超声探头阵列、控制系统、磁体系统、电极对或线圈，以及信号检测处理系统。超声驱动激励源和超声探头连接，控制系统和超声探头连接；电极对或线圈、信号检测处理系统、图像重建系统依次连接。通过控制系统使超声探头阵列处于发射或测量两种模式。

本发明方法根据以下建立的电压信号与电导率的非线性关系，重建断层面的电导率图像。

重建过程包括四个步骤：

1、利用互易定理建立实际测量过程与假想过程物理量的对应关系

（1）实际测量过程为：在超声激励和静磁场共同作用聚焦区中，质点振动速度为Ve_x，静磁场为Be_z，所述的静磁场由磁体系统产生，其对应的磁感应强度B是已知的，质点振动速度V是通过超声回波测定的。电极对或线圈测量到的电压信号为U，e_x和e_z分别为直角坐标系下x和z方向的单位矢量。

（2）假想过程为：去除超声激励和静磁场，将实际测量过程中的电极对或线圈，改为注入电极或激励线圈，并通入单位电流，设此时在聚焦区的电流密度的y分量为J，聚焦区很小，则根据互易定理有：

U＝∫_SJe_y·(Ve_x×Be_z)ds （1）

上式中，e_y为直角坐标系下y方向的单位矢量。

2、根据测得的电压U重建超声激励和静磁场共同作用聚焦区内假想过程的电流密度的y分量；

实际上超声激励和静磁场共同作用聚焦区面积S很小，可以假定所述的聚焦区内质点振动速度、静磁场、电流密度均匀，则根据公式（1），有：

$J=\left|\frac{U}{\mathrm{VBS}}\right|---\left(2\right)$

公式（2）中||表示绝对值，J表示电流密度y分量分布，V为质点振动速度。

3、通过控制系统控制超声探头阵列依次扫描移动所述的聚焦区的位置，重复步骤1和步骤2，直到所有聚焦区覆盖整个断层面，从而获得断层面内假想过程的电流密度y分量分布J。

4、根据断层面内假想过程的电流密度y分量分布J重建电导率分布。

（1）给定断层面内所有所述聚集区内的电导率初值σ，通常选择电导率初值为1s/m，求解电位u满足的方程获得断层面内电位分布，进一步，利用计算电流密度y分量。这里δ为狄拉克函数，r为断层面内任意位置坐标，r_a和r_b分别为电极对的位置坐标。

（2）对目标函数||J-J_c||₂极小化，迭代获得电导率分布。这里|| ||₂为二范数。

本发明磁声电成像系统的工作过程如下：

所述的超声驱动激励源发射脉冲超声激励信号，所发射的激励信号通过超声驱动电缆发送至超声探头阵列。超声探头阵列通过耦合剂与生物体接触。通过控制系统控制超声探头阵列发射聚焦超声，在生物体的组织内部一定深度的局部区域中激发声辐射力，引起生物体组织局部区域断层面内的局部聚焦区质点振动。磁体系统在所述的生物体组织局部区域断层面产生均匀静磁场，这里静磁场的大小是确定的也是已知的。聚焦区质点振动速度和均匀静磁场方向垂直，质点振动速度和均匀静磁场二者共同作用引起聚焦区内产生洛伦兹力。由于洛伦兹力的作用，引起聚焦区电荷分离，产生库仑电场，在所述的断层面内产生电流分布。利用贴放在所述的断层面上的电极对或线圈测量电压信号，通过信号检测处理系统记录电压信号。超声探头阵列通过控制系统切换到测量模式，通过超声回波测量获得聚焦区的位置坐标以及质点振动速度，并通过信号检测处理系统记录。利用测量到的电压信号、聚焦区位置信息、聚焦区质点振动速度，以及已知的静磁场信息，采用图像重建系统重建获得电导率图像。

本发明的主要激励装置为超声探头阵列和静磁场，而检测装置为电极对或线圈、超声探头阵列。静磁场技术和超声波激励检测技术皆为医学领域传统技术，易于实现。这种磁声电成像方法综合了超声成像和电成像技术的优点，有广泛的应用前景和潜在的应用价值。

附图说明

图1本发明装置结构示意图；

图2本发明磁声电成像系统装置在采用电极对检测过程中与生物体相对位置关系图；

图3本发明磁声电成像系统装置在采用线圈检测过程中与生物体相对位置关系图。

图中：10超声驱动激励源，20超声探头阵列，30耦合剂，40磁体系统，50电极对或线圈，60信号检测处理系统，70图像重建系统，80控制系统，90生物体，501和502一对电极，503线圈，110声辐射力聚焦区，110断层面，130聚焦区质点振动速度，140磁体系统40产生的静磁场。

具体实施方式

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

如图1所示，本发明磁声电成像系统主要包括超声驱动激励源10、超声探头阵列20，耦合剂30，磁体系统40，电极对或线圈50，信号检测处理系统60，图像重建系统70，控制系统80。另外还包括生物体90。

超声驱动激励源10与超声探头阵列20连接，控制系统80与超声探头阵列20连接，电极对或线圈50、信号检测处理系统60、图像重建系统70依次连接。磁体系统40可以采用永磁体、电磁体或超导磁体实现。

本发明磁声电成像系统装置采用电极对检测过程中与生物体相对位置关系如图2所示。超声探头阵列20通过耦合剂30和生物体90的某一部位接触，120为所述生物体90某一部位的断层面，110为断层面内的声辐射力聚焦区，130为聚焦区质点振动速度。磁体系统40在断层面120处产生均匀静磁场140。电极对501和502贴放在生物体的断层面120的外边界上，离开耦合剂30的任意位置，两个电极501和502之间的距离任意。

本发明磁声电成像系统装置采用线圈检测过程中与生物体相对位置关系如图3所示。除了不采用电极对501和502外，其他部件之间的位置关系与图2所示的采用电极对检测过程中与生物体相对位置关系的本发明装置相同。线圈非接触地位于生物体的断层面120的外边界的外侧，离开耦合剂30的任意位置。

本发明工作过程如下：

所述的超声驱动激励源10发射脉冲超声激励信号，所述的激励信号通过超声驱动电缆发送至超声探头阵列20。超声探头阵列20通过耦合剂30与生物体90接触，控制系统80控制超声阵列探头发射聚焦超声，在生物组织内部一定深度的局部区域中激发声辐射力，引起断层面120内的局部聚焦区110质点振动，质点振动速度为130。磁体系统40在断层面120处产生均匀静磁场140，质点振动速度130和均匀静磁场140方向垂直。由质点振动速度130和均匀静磁场140共同作用引起聚焦区110内产生洛伦兹力，由于洛伦兹力的作用，引起聚焦区电荷分离，产生库仑电场，在断层面内产生电流分布。如图2所示利用贴放在断层面120上的电极对501和502测量电压信号，或者如图3所示利用断层面120外的非接触线圈503测量电压信号，通过信号检测处理系统60记录电压信号。通过控制系统80，超声探头阵列20切换到测量模式，通过回波测量获得聚焦区110的位置坐标以及质点振动速度130，并通过信号检测处理系统60记录。利用测量到的电压信号、聚焦区位置信息、聚焦区质点振动速度130，采用图像重建系统70重建获得电导率图像。

重建断层面的电导率图像的步骤如下：

重建过程包括四个步骤：

1、利用互易定理建立实际测量过程与假想过程物理量的对应关系；

（1）实际测量过程为：在超声激励和静磁场共同作用的聚焦区110中，质点振动速度130为Ve_x，静磁场为Be_z，所述的均匀静磁场由磁体系统40产生，其对应的磁感应强度B是已知的，质点振动速度V是通过超声回波测定的。电极对501和502或线圈503测量到的电压信号为U，e_x和e_z分别为直角坐标系下x和z方向的单位矢量。

（2）假想过程为：去除超声激励和静磁场，将实际测量过程中的电极对或线圈改为注入电极或激励线圈，并通入单位电流，设此时在聚焦区的电流密度的y分量为J，聚焦区很小，则根据互易定理有：

U＝∫_SJe_y·(Ve_x×Be_z)ds （1）

上式中，e_y为直角坐标系下y方向的单位矢量。

2、根据测得的电压U重建超声激励和静磁场共同作用聚焦区内假想过程的电流密度的y分量；

实际上超声激励和静磁场共同作用聚焦区面积S很小，可以假定所述的聚焦区内质点振动速度、静磁场、电流密度均匀，则根据公式（1），有：

$J=\left|\frac{U}{\mathrm{VBS}}\right|---\left(2\right)$

公式（2）中||表示绝对值。

3、通过控制系统控制超声探头阵列依次扫描移动所述的聚焦区的位置，重复步骤1和步骤2，直到所有聚焦区覆盖整个断层面，从而获得断层面内假想过程的电流密度y分量分布J。

4、根据断层面内假想过程的电流密度y分量分布J重建电导率分布；

（1）给定断层面内所有所述聚集区内的电导率初值σ，通常选择电导率初值为1s/m，求解电位u满足的方程获得断层面内电位分布，进一步，利用计算电流密度y分量。这里δ为狄拉克函数，r为断层面内任意位置坐标，r_a和r_b分别为电极对的位置坐标。

（2）对目标函数||J-J_c||₂极小化，迭代获得电导率分布。这里|| ||₂为二范数。

Claims (2)

1.一种磁声电成像系统，其特征在于，所述的磁声电成像系统包括超声驱动激励源(10)、超声探头阵列(20)、耦合剂(30)、磁体系统(40)、电极对或线圈(50)、信号检测处理系统(60)、图像重建系统(70)和控制系统(80)；所述的超声驱动激励源(10)与超声探头阵列(20)连接，控制系统(80)与超声探头阵列(20)连接；电极对或线圈(50)、信号检测处理系统(60)、图像重建系统(70)依次顺序连接；所述的磁体系统(40)为永磁体或电磁体或超导磁体；超声驱动激励源(10)提供两种频率接近的超声波；所述的控制系统(80)控制超声探头阵列处于发射模式或测量模式；贴放在生物体组织局部区域断层面(120)上的电极对或线圈(50)测量电压信号，信号检测处理系统(60)记录所述的电压信号；超声探头阵列(20)通过控制系统(80)切换到测量模式，通过超声回波测量获得聚焦区的位置坐标以及质点振动速度，由信号检测处理系统(60)记录。

2.一种应用权利要求1所述的磁声电成像系统的成像方法，其特征在于采用声辐射力在生物组织内部产生聚焦区，采用超声回波技术测量聚焦区位置和质点振动速度；

所述的聚焦区面积小，使得聚焦区内静磁场、质点振动速度和假想过程的电流密度能够近似为均匀值，从而建立起测量电压与假想过程的电流密度的线性关系；

所述的超声驱动激励源(10)发射脉冲超声激励信号，所发射的激励信号发送至超声探头阵列(20)；超声探头阵列(20)通过耦合剂(30)与生物体(90)接触；控制系统(80)控制超声探头阵列(20)发射聚焦超声，在生物体的组织内部的局部区域中激发声辐射力，引起所述生物体的组织内部局部区域断层面内的局部聚焦区质点振动；磁体系统(40)在所述的生物体组织局部区域断层面产生均匀静磁场；所述的聚焦区质点振动速度和均匀静磁场方向垂直，质点振动速度和均匀静磁场二者共同作用引起聚焦区内产生洛伦兹力；由于洛伦兹力的作用，引起聚焦区电荷分离，产生库仑电场，在所述的断层面内产生电流分布；利用贴放在所述的断层面上的电极对或线圈(50)测量电压信号，通过信号检测处理系统(60)记录电压信号；超声探头阵列(20)通过控制系统(80)切换到测量模式，通过超声回波测量获得聚焦区的位置坐标以及质点振动速度，并通过信号检测处理系统(60)记录；利用测量到的电压信号、聚焦区位置信息、聚焦区质点振动速度，以及已知的静磁场信息，采用图像重建系统重建获得电导率图像。