技术领域
[0001] 本
发明涉及一种医学成像方法及装置,特别涉及一种磁声
电阻抗成像方法及装置。
背景技术
[0002] 传统的电阻抗成像的灵敏度和空间
分辨率不高,为解决此问题,美国明尼苏达大学(YuanXu and Bin He 2005 Magnetoacoustic tomography with magnetic induction(MAT-MI)Phys.Med.Biol.50 5175-5187)提出了感应式磁声电阻抗成像方法。该方法的主要原理为:成像目标置于静
磁场B0中,用脉冲磁场B1
电流激励成像目标后在其中产生感应电流J,电流在静磁场作用下产生洛仑兹
力F=J×B0,从而激发出声压
信号p。这种成像方法依据理论
基础为声压
波动方程 (这里cs为声速)。提取声压信号,通过洛
伦兹力散度重建公式 得到洛伦
兹力的散度 由于 则可获得电流
密度的旋度 再利用电流密
度的无散性 以及电流密度的法向分量为零 的边界条件,重建出电流密度分布,进一步利用 重建得到成像目标的电导率σ分布。该方法将
电磁感应技术与超声技术相结合,其分辨率得到很好的改善。然而,该方法依据的理论基础仍存在
缺陷,引起重建方法存在奇异性问题。具体为:(1)在洛伦兹力的散度 重建过程中,由于在成像体边界处,电流密度的法向分量不连续,导致洛伦兹力F不连续,进一步导致洛伦兹力的散度 不存在,或者说,在边界处 存在奇异性,则边界处重建的 不合理;(2)在电流密度的重建过程中,由于边界处电流密度的不连续性,导致电流密度的旋度 不存在,或者说,在边界处 存在奇异性,因此,重建的电流密度J不合理;(3)在电导率的重建过程中,公式中既包含 又包含重建出的J,不合理再次被放大。
发明内容
[0003] 本发明的目的是克服
现有技术的缺点,提出一种磁声电阻抗成像方法及装置。本发明可用于测量
生物体电阻抗分布。
[0004] 本发明利用连接激励线圈的脉冲激励器,使激励线圈产生瞬变电流JS,在成像体中产生感应电流,感应电流在静磁场B0作用下产生洛仑兹力,洛伦兹力引起成像体质点振动,激发超声信号,超声信号传播到成像体外,利用成像体外放置的
加速度
传感器测量加速度传感器所在
位置的位移
波形u。根据位移波形u与电导率的非线性关系,借助于如下理论基础和重建方法得到成像体的电导率图像:
[0005]
生物组织的特性与
流体类似,在洛伦兹力作用下,生物组织中的动力学方程为[0006]
[0007] 这里ρ0为生物组织密度,cs为声速。
[0008] 与明尼苏达大学依据的理论不同,本发明采用的方程中不存在洛伦兹力的散度项,在边界处洛伦兹力不连续并不影响数学物理模型的合理性。本发明的主要创新在于通过测量的位移波形信号重建洛伦兹力,而不是采用声压信号重建洛伦兹力的散度,因为洛伦兹力的散度有奇异性,导致重建时图像分辨率低。
[0009] 本发明的图像重建过程包括3个步骤,具体重建过程如下:
[0010] 步骤1,根据测得的位移波形信号,重建成像体
断层面的洛伦兹力F。
[0011] 本发明采用加速度传感器测量位移波形信号,根据测量的位移波形信号,采用时间翻转方法重建成像体断层面的洛伦兹力F,具体重建公式如下:
[0012]
[0013] 其中,∑是环绕成像体的检测面,加速度传感器布置在检测面上,rd加速度传感器的位置,n是单位法向矢量,r是成像体内断层面上点的位置,dSd是检测面上的面元,u″(rd,|rd-r|/cs)是位移波形u的二阶时间导数,B0同前文所述,为静磁场的磁通密度,J为断层面上的电流密度。
[0014] 步骤2,根据步骤1得出的洛仑兹力F重建出断层面的电流密度J。
[0015] 将F=J×B0展开,并考虑到静磁场B0=B0ez,有
[0016] F=JyB0ex-JxB0ey
[0017] 这里ex、ey和ez分别为x,y,z三个方向的单位矢量,Jx和Jy分别是电流密度的x,y分量。
[0018] 由于静磁场B0已知,则断层面上的电流密度为
[0019] J=Fy/B0ex-Fx/B0ey (2)
[0020] 这里Fx和Fy分别是洛伦兹力的x,y分量。
[0021] 步骤3,根据断层面上的电流密度J重建电导率分布σ。
[0022] 引入矢量磁位A和标量电位Ф,将重建区域Ω剖分成n个单元,第i个单元区域为Ωi,给定单元区域Ωi的电导率初值σi,利用
电磁场领域对于
涡流问题常用的A-Ф有限元数值计算方法,计算出成像体的矢量磁位A和标量电位Ф,进一步计算出电流密度Jc,使之与步骤2重建出的电流密度J之差在最小二乘的意义下最小,即对如下目标函数进行最小化:
[0023]
[0024] 其中‖Jc-J‖2表示Jc和J之差的二范数,则有更新后的电导率为,[0025] i=1,2,…,n (3)
[0026] 通过
迭代,最终获得断层面上的电导率分布。
[0027] 本发明方法具有以下特征:
[0028] (1)本发明方法采用加速度传感器
测量传感器所在位置的位移波形;
[0029] (2)依据测量的位移波形量重建洛伦兹力,进一步重建电导率分布。
[0030] 应用本发明方法的磁声电阻抗成像系统包括脉冲激励器、激励线圈、静磁体、加速度传感器阵列、
数据采集器、计算机。其中激励线圈和静磁体置于成像体两侧。加速度传感器阵列环绕成像体分布。脉冲激励器通过
电缆连接激励线圈。加速度传感器阵列、数据采集器、计算机依次连接。
附图说明
[0031] 图1为电导率成像系统的原理图。
具体实施方式
[0032] 以下结合附图和具体实施方式进一步说明本发明。
[0033] 本发明磁声电阻抗成像方法是利用连接激励线圈的脉冲激励器使激励线圈产生瞬变电流JS,在成像体中产生感应电流,感应电流在静磁场B0作用下产生洛仑兹力F,洛伦兹力F引起成像体质点振动,激发超声信号,超声信号传播到成像体外。本发明利用加速度传感器测量加速度所在位置的位移波形u。根据位移波形u与电导率的非线性关系,重建得到成像体的电导率图像。
[0034] 如附图所示,应用本发明方法的磁声电阻抗成像系统包括脉冲激励器、激励线圈、静磁体、加速度传感器阵列、数据采集器、计算机。其中激励线圈和静磁体置于成像体两侧。加速度传感器阵列环绕成像体分布。脉冲激励器通过电缆连接激励线圈。加速度传感器阵列、数据采集器、计算机依次连接。
[0035] 本发明装置的工作过程为:启动脉冲激励器,连接在脉冲激励器上的激励线圈产生瞬变电流,在成像体中产生感应电流,感应电流在静磁场作用下产生洛仑兹力,洛伦兹力引起成像体质点振动,激发超声信号,超声信号传播到成像体外,加速度传感器测量到位移波形信号,通过数据采集器采集位移波形信号,将此信号存储在计算机中。
[0036] 根据存储的位移波形信号,利用公式(1)~(3)所示的位移波形信号和电导率的非线性关系实现图像重建。具体步骤为:
[0037] 步骤1,根据测得的位移波形信号,重建成像体断层面的洛伦兹力F。
[0038] 本发明采用加速度传感器测量位移波形信号,根据测量的位移波形信号,采用时间翻转方法重建成像体断层面的洛伦兹力F,具体重建公式如下:
[0039]
[0040] 其中,∑是环绕成像体的检测面,加速度传感器布置在检测面上,rd加速度传感器的位置,n是单位法向矢量,r是成像体内断层面上点的位置,dSd是检测面上的面元,u″(rd,|rd-r|/cs)是位移波形u的二阶时间导数,B0同前文所述,为静磁场的磁通密度,J为断层面上的电流密度。
[0041] 步骤2,根据步骤1得出的洛仑兹力F重建出断层面的电流密度J。
[0042] 将F=J×B0展开,并考虑到静磁场B0=B0ez,有
[0043] F=JyB0ex-JxB0ey
[0044] 这里ex、ey和ez分别为x,y,z三个方向的单位矢量,Jx和Jy分别是电流密度的x,y分量。
[0045] 由于静磁场B0已知,则断层面上的电流密度为
[0046] J=Fy/B0ex-Fx/B0ey
[0047] 这里Fx和Fy分别是洛伦兹力的x,y分量。
[0048] 步骤3,根据断层面上的电流密度J重建电导率分布σ。
[0049] 引入矢量磁位A和标量电位Ф,将重建区域Ω剖分成n个单元,第i个单元区域为Ωi,给定单元区域Ωi的电导率初值σi,利用电磁场领域对于涡流问题常用的A-Ф有限元数值计算方法,计算出成像体的矢量磁位A和标量电位Ф,进一步计算出电流密度Jc,使之与步骤2重建出的电流密度J之差在最小二乘的意义下最小,即对如下目标函数进行最小化:
[0050]
[0051] 其中‖Jc-J‖2表示Jc和J之差的二范数,则有更新后的电导率为,[0052] i=1,2,…,n
[0053] 通过迭代,最终获得断层面上的电导率分布,在显示器上输出重建的电导率。
[0054] 在本
实施例中,激励线圈的平均半径为5cm,电感为17μH。瞬变电流的波形的脉冲宽度为1μs,峰值为200A。
[0055] 在本实施例中,成像体为一半径为50mm,厚度为4mm的薄圆盘状琼脂,将琼脂置于流体中。薄圆盘状琼脂和流体的电导率分别为1s/m和0,二者声速cs均为1.5mm/μs。加速度传感器阵列包含20个传感器,置于与薄圆盘状琼脂同轴的半径为100mm的圆周上,各自之间在圆周上间隔18度。