首页 / 专利库 / 测量仪器和方法 / 电阻抗成像 / 一种磁声电阻抗成像方法及装置

一种磁声电阻抗成像方法及装置

阅读:816发布:2020-07-29

专利汇可以提供一种磁声电阻抗成像方法及装置专利检索,专利查询,专利分析的服务。并且一种磁声 电阻 抗成像 方法,其特征在于利用 加速 度 传感器 测量由于脉冲 磁场 激励的感应 涡流 在静态磁场作用下产生的洛伦兹 力 振动位移 波形 信号 ,根据位移波形信号与电导率的非线性关系,通过图像重建得到成像体的电导率图像。应用本 发明 方法的装置,包括脉冲激励器、激励线圈、静磁体、加速度传感器阵列、 数据采集 器和计算机。其中激励线圈和静磁体置于成像体两侧。加速度传感器阵列环绕成像体分布。脉冲激励器通过 电缆 连接激励线圈。加速度传感器阵列、数据采集器、计算机依次连接。,下面是一种磁声电阻抗成像方法及装置专利的具体信息内容。

1.一种磁声电阻抗成像方法,利用连接激励线圈的脉冲激励器使激励线圈产生瞬变电流JS,在成像体中产生感应电流,感应电流在静磁场B0作用下产生洛仑兹F,洛伦兹力F引起成像体质点振动,激发超声信号,超声信号传播到成像体外,其特征在于利用加速传感器测量所在位置的位移波形u,根据位移波形u与电导率的非线性关系,重建得到成像体的电导率图像。
2.根据权利要求1所述的磁声电阻抗成像方法,其特征在于图像重建过程包括三个步骤:
步骤1、根据测得的位移波形u重建断层面的洛仑兹力F=J×B0;
步骤2、根据洛仑兹力F重建断层面的电流密度J;
步骤3、根据断层面上的电流密度J重建电导率σ分布;
其中:B0为静磁场磁场强度。
3.根据权利要求2所述的磁声电阻抗成像方法,其特征在于所述的步骤2重建断层面的电流密度J的具体过程为:
根据公式F=JyB0ex-JxB0ey
计算断层面上的电流密度J为:
J=Fy/B0ex-Fx/B0ey
其中:Fx和Fy分别是洛伦兹力的x,y分量,Jx和Jy分别是电流密度的x,y分量,ex和ey分别为x,y两个方向的单位矢量。
4.根据权利要求2所述的磁声电阻抗成像方法,其特征在于所述的步骤3根据断层面上的电流密度J重建电导率σ分布的具体过程为:
引入矢量磁位A和标量电位Ф,将重建区域Ω剖分成n个单元,第i个单元区域为Ωi,给定单元区域Ωi的电导率初值σi,利用A-Ф有限元数值计算方法,计算出成像体的矢量磁位A和标量电位Ф,进一步计算出电流密度Jc,使之与步骤2重建出的电流密度J的差在最小二乘的意义下最小,即对如下目标函数进行最小化:
其中‖Jc-J‖2表示Jc和J之差的二范数,则有更新后的电导率为,
i=1,2,…,n
通过迭代,最终获得断层面上的电导率分布。
5.应用权利要求1所述的磁声电阻抗成像方法的电阻抗成像系统,包括脉冲激励器、激励线圈、静磁体、数据采集器和计算机,激励线圈和静磁体置于成像体两侧,脉冲激励器通过电缆连接激励线圈;加速度传感器阵列、数据采集器、计算机依次连接,其特征在于所述的电阻抗成像系统包括加速度传感器阵列,加速度传感器阵列环绕成像体分布。

说明书全文

一种磁声电阻抗成像方法及装置

技术领域

[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度。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈