首页 / 专利库 / 物理 / 斯托克定律 / 一种基于磁声电原理的电导率成像系统

一种基于磁声电原理的电导率成像系统

阅读:838发布:2020-07-01

专利汇可以提供一种基于磁声电原理的电导率成像系统专利检索,专利查询,专利分析的服务。并且一种基于磁声电原理的电导率成像系统,其成像平台将采集到的磁声电 信号 传输到图像重建模 块 。所述的成像平台包括声场驱动激励模块、 磁场 激励模块和检测模块;声场驱动激励模块产生声场激励源;磁场激励模块为开放式磁体,放置于 生物 体附近,产生的非均匀静磁场作用于生物体,生物体在非均匀静磁场 和声 场的作用下产生动生源 电流 ,检测模块采集动生源电流,将其转化为磁声电 电压 信号。本 发明 利用 电极 检测磁声 电信号 ,通过图像重建模块得到成像体的电导率分布图像,实现 生物组织 目标区域内的电导率分布的检测。,下面是一种基于磁声电原理的电导率成像系统专利的具体信息内容。

1.一种基于磁声电原理的电导率成像系统,其特征在于,所述的成像系统包括成像平台和图像重建模;成像平台与图像重建模块连接,将采集到的磁声电信号通过传输线传输到图像重建模块;所述的成像平台包括声场驱动激励模块、磁场激励模块和检测模块;声场驱动激励模块产生声场激励源;磁场激励模块为开放式磁体,放置于生物体附近,产生的非均匀静磁场作用于生物体,生物体在非均匀静磁场和声场的作用下产生动生源电流,检测模块采集动生源电流将其转化为磁声电电压信号;
所述的声场驱动激励模块由超声驱动激励源、超声阵列及耦合囊组成;超声阵列的一端连接超声驱动激励,超声阵列的另一端和耦合水囊接触;超声驱动激励源激励超声阵列产生声波;耦合水囊中充满水,耦合水囊填充于超声阵列和生物体之间的空间,使超声阵列产生的超声波能够传播到生物体中;
检测模块由电极、滤波电路、放大电路及信号采集装置组成;电极和生物体接触,检测生物体表面的电压信号;电极的一端和生物体连接,电极的另一端和滤波电路的输入端连接,滤波电路的输出端接放大电路的输入端,放大电路的输出端连接信号采集装置的输入端,信号采集装置的输出端和图像重建模块相连;
图像重建模块根据信号采集装置输出的生物体的磁声电电压信号重建电导率分布。
2.按照权利要求1所述的基于磁声电原理的电导率成像系统,其特征在于,所述的图像重建模块利用所述生物体的磁声电电压信号,以及已知的非均匀静磁分布信息,采用图像重建算法实现电导率分布的重建;所述的图像重建模块采用两种电导率重建算法实现电导率分布的重建,重建算法一为直接代数迭代电导率重建算法、重建算法二为等效均匀场代数迭代电导率重建算法。
3.按照权利要求2所述的基于磁声电原理的电导率成像系统,其特征在于,所述的算法一为直接代数迭代电导率重建算法,包括以下三个步骤:
(1)利用互易定理建立实际测量过程与互易过程物理量的对应关系;
所述的实际测量过程为:在超声驱动激励模块产生的超声波和磁场激励模块产生的非均匀静磁场作用下,质点的振速为v,非均匀静磁场为B0(r),所述的磁场由磁场激励模块的开放式磁体产生,所述的非均匀静磁场B0(r)强度的分布是已知的;
所述的互易过程为:关闭声场驱动激励模块,并使磁场激励模块不工作,向检测模块中的电极通入I安培的直流电,此电流在生物体中产生的电流密度为Jr(r);
在磁场激励模块产生的非均匀静磁场B0(r)的激励下,基于互易定理得到磁声电电压分布u(r,t)和振动速度势 互易过程电流密度Jr(r),以及磁场激励模块产生的非均匀静磁场B0(r)之间的关系式为:
式(1)中,t表示超声波的传播时间,Ω表示生物体所在的区域,r表示场点,即生物体所在区域Ω中的点, 表示振动速度势,ρ0为生物体的密度,Jr(r)为互易过程中生物体的电流密度,B0(r)为磁场激励模块产生的非均匀静磁场分布, 为散度算符;
(2)根据实际测量过程中测得的磁声电电压u(r,t),依据公式(1)重建互易过程电流密度Jr(r)和静磁场B0(r)之间的关系
振动速度势 满足的格林函数为 将其带入公式(1),根据格林函
数的对称性,可知磁声电电压u(r,t)满足:
式(2)中,r'为源点,表示超声阵列所在区域的点;
利用公式(2)得到互易过程的电流密度Jr(r)和开放式磁体产生的非均匀静磁场B0(r)之间的关系式,同时,由于ρ0为常数,因此得到:
式(3)中,c0表示声场驱动激励模块产生超声波的速度,trd=2T0-t+|r-r'|/c0表示反转场的时间,t表示超声波的传播时间,T0表示反转u(r,t)的时刻,r'表示超声阵列所在区域的点,r表示生物体所在区域Ω中的点,S表示生物体所在区域Ω的表面,n是生物体所在区域Ω外边界处的单位法向矢量,u'(r',trd)表示u(r',trd)的一阶导数,u”(r',trd)表示u(r',trd)二阶导数;
利用式(3)实现变量 的重建;
(3)根据变量H(r)重建电导率分布;
利用公式(3)实现变量 的分布重建之后,在生物体组织每一个断
层面z0上变量H(r)为 表示为 即:
其中各个变量括号中(x,y,z0)表示相对应的变量在生物体内部每一个断层面z0上坐标;
Jr(x,y,z0)为z0断层面上互易过程的电流密度,利用欧姆定律,该电流密度Jr(x,y,z0)表示为互易过程电位的梯度和电导率之间的乘积,即:
因此f(x,y,z0)表示为
其中ur(x,y,z0)是互易过程的电位在z0断层面的分布,σ(x,y,z0)表示生物体的电导率在z0断层面的分布;
互易过程的电位在z0断层面的分布满足如下关系式:
其中,I为互易过程中注入I安培的直流电,rA和rB表示检测模块中电极对的位置,Γ表示目标体所在区域Ω的表面,n表示表面Γ的外法线方向的单位矢量;
每一个断层面z0上的 的值表示为f(x,y,z0),由f(x,y,z0)重建
电导率σ的方法步骤如下:
1)将生物体划分为一系列的子块,认为这些子块内部的电导率均匀,给出生物体的电导率分布矩阵[σ]的初值,通常选择电导率的初值为0.1S/m,给定误差精度ε;
2)根据公式(5)计算得到互易过程的电标位ur(x,y,z0);
3)利用公式(4)重建z=z0的断层面的电导率分布;
4)利用步骤3)得到的各个断层的电导率分布,通过线性差值得到各个子块的电导率分布,如此得到第k次迭代的生物体整个三维区域的电导率分布σk;
k k+1
5)计算第k次迭代得到的生物体电导率分布σ和第k+1次生物体电导率分布σ 之间的相对误差,比对该相对误差是否满足给定的误差精度ε,如果满足给定的误差精度ε,停止迭代;否则把第k次得到的生物体的电导率分布σk作为初始电导率分布,转到步骤2),以上过程依次迭代下去,直到相邻两次计算得到的生物体电导率分布的相对误差满足精度要求;
利用公式(1)和(2)能够实现 的重建,再利用以上步骤1)~步骤5)实现生物体电导率图像的重建。
4.按照权利要求2所述的基于磁声电原理的电导率成像系统,其特征在于,所述的算法二等效均匀场代数迭代电导率重建算法基于磁声电电压信号和非均匀静磁场B0(r)之间的关系,以非均匀静磁场B0(r)作为磁场激励源的磁声电电压信号等效为均匀静磁场激励的磁声电电压信号的原理如下公式(6)-(14)所示:
给定三维模型,在非均匀静磁场B0(r)的激励下,生物体中等效电流源对应的电流I(t)和声场驱动激励模块产生的超声波在生物体中的振速v之间的关系为:
I(t)=∫sσv×B0(r)·dS            (6)
其中S表示等效电流源流过的面的面元;
根据声学原理,超声动冲量M和振速v需满足:
其中ρ0为生物体的密度, 为梯度算符;
将式(7)代入式(6),可得:
其中n表示等效电流源流过的面的面元的法向方向的单位矢量;
进一步利用斯托克斯公式和矢量恒等式,式(8)化简为:
其中,l表示面元的外边缘的线,且l的正向与S的外法线方向符合右手定理,S的法向方向为n;
考虑实际应用时,超声阵列发出的能量频率范围包含极少的直流频率,表示超声阵列产生的波包的净动量为零,因此式(9)等号右边的第一项为零;
同时在实际检测中,电极只能检测到生物体的一部分电流,因此信号采集装置采集到的电流只是电流I(t)的一部分,定义采集电压和电流之间的比例为α,检测到的磁声电电压U(t)表示为:
式(11)即检测到的磁声电电压U(t)、非均匀静磁场B0(r)、电导率σ以及密度ρ0之间的关系式;
实际应用时, 为一个很小的量,因此式(11)等号右边第二项
为一个小量,忽视此小量,故式(11)简化为:
公式(12)中,当等效电流源流过的面的面元的法向方向的单位矢量n的方向与的方向相同时,式(12)表示为:
其中 取此变量的ex方向分量,则公式(13)表示为:
其中,Boy和Boz分别为B0(r)沿y方向和z方向的磁场分布,θ表示 和B0(r)之间的夹
式(14)为非均匀静磁场激励下,磁声电电压信号等效为均匀静磁场激励下的磁声电电压信号的理论关系式,非均匀静磁场激励下,在电导率发生变化的位置磁声电信号会被放大或者缩小 倍;
所述算法二等效均匀场代数迭代电导率重建算法步骤为:
首先利用式(14)将非均匀静磁场激励下采集的磁声电电压信号uin(r,t)调整为等效的均匀场激励下采集的磁声电电压信号uho(r,t);
然后将公式(1)-(5)中所有的非均匀静磁场B0(r)替换为均匀静磁场B0;
最后利用算法一直接代数迭代电导率重建算法的步骤1)-5),实现电导率的重建。

说明书全文

一种基于磁声电原理的电导率成像系统

技术领域

[0001] 本发明涉及一种开放式磁声电电导率成像系统。

背景技术

[0002] 实现疾病的早期诊断意义重大,一方面可以减少患者治疗过程中的痛苦,以及高昂的医疗费,更重要的是可以提高存活率,研究表明组织的电特性的变化早于组织结构的病变,因此以电参数为成像目标参数的成像方法可以实现疾病的早期诊断。磁声电成像以电参数为成像目标参数,具有良好应用前景的医学成像方法兼具高对比度和高分辨率的优势。
[0003] 1998年,Han wen等人提出霍尔效应成像,并且在4.0T磁共振磁场静磁场的激励下,利用电极检测到培根肉的实验信号,但是实验系统的配置并没有过多提及。2003年,Hanwen在专利US6520911B1中详细论述了此成像系统,在系统的配置过程中,利用均匀静磁场作为磁声电成像的磁场激励源。2007年Y.Xu、S Haider等人在霍尔效应的基础上提出了基于互易定理的磁声电成像,利用实验的方法得到了互易电流密度的分布图,同时对互易过程的电流密度做了仿真分析,然而并没有给出电导率的重建方法。针对均匀静磁场作为磁声电成像的磁场激励源,电导率重建算法有很多学者做了研究工作,2014年,法国里昂一大的Ammari H,Grasland-Mongrain P.等人报道了磁声电成像的理论分析和不同信噪比条件下的仿真研究(Ammari H等.,2014)。同年,中国科学院电工研究所的郭亮利用时间发转法、压缩感知和拟迭代算法,通过仿真分析实现了电极检测式磁声电成像电导率的重建(Guo L等.,2014)。Kunyansky L 2017年报道了电导率边界的重建算法,设计了一个3D扫描平台,利用电极检测到了牛肉组织的磁声电信号,并且进行了牛肉组织分界面信息的重建,并没有通过实验信号重建得到电导率的图像。以上的报道的实验平台要么是在均匀静磁场的激励下,要么是针对的目标体是一个小型的仿体结构,因此假定静磁场均匀是合理的。
[0004] 综上,目前还没有实际医用系统的报道。对于实际的医用系统有两种解决方案实现静磁场,(1)利用均匀的静磁场作为磁场激励;(2)利用非均匀静磁场作为磁场激励源。对于方案(1)均匀静磁场的实现可以借鉴磁共振成像中静磁场的设计理念,采用永磁体、电磁体或者超导磁体实现,电导率的重建可以借鉴前期的研究成果,但是均匀度高的静磁场的实现会极大的提高静磁场产生装置的成本,进一步大大提高磁声电成像临床应用设备的造价。而且封闭的环境减小了成像的区域,也不适用于幽闭症病人。对于方案(2)不需要完全均匀的静磁场产生装置,极大降低了医用系统的造价,但是,目前研究的重建算法并不适用。
[0005] 目前存在的磁声电成像的平台难以转化成医用的成像系统,存在以下缺点:(1)目前理论模型针对均匀的静磁场,实现完全均匀的静磁场大大提高了设备的造价;(2)非均匀静磁场作为磁场激励源可降低设备的造价,但是目前无相应的重建算法。

发明内容

[0006] 本发明的目的是克服现有技术的缺点,提出一种基于磁声电原理的电导率成像系统。
[0007] 本发明基于磁声电原理的电导率成像系统包括成像平台和图像重建模。成像平台与图像重建模块连接,将采集到的磁声电信号通过传输线传输到图像重建模块。
[0008] 所述的成像平台包括声场驱动激励模块、磁场激励模块和检测模块。声场驱动激励模块产生声场激励源,即声波超声波在空气中衰减很快,为了能更好的传播于生物体中,需要耦合囊,耦合水囊和声场驱动激励模块完全接触,以减小声波的衰减。磁场激励模块产生的非均匀静磁场作用于生物体,生物体在非均匀静磁场和声场的作用下会产生动生源电流,检测模块采集动生源电流,将其转化为磁声电电压信号。
[0009] 所述的声场驱动激励模块由超声驱动激励源、超声阵列及耦合水囊组成。超声阵列的一端连接超声驱动激励,超声阵列的另一端和耦合水囊接触。超声驱动激励源激励超声阵列产生超声波。耦合水囊中充满水,耦合水囊填充于超声阵列和生物体之间的空间,使超声阵列产生的超声波能够传播到生物体中。
[0010] 磁场激励模块为开放式磁体,放置于生物体附近。
[0011] 检测模块由电极、滤波电路、放大电路及信号采集装置组成。电极和生物体接触,检测生物体表面的电压信号,滤波和放大电路实现检测信号的滤波和放大,最终信号采集装置实现信号的采集。所述的电极的一端和生物体连接,电极的另一端和滤波电路的输入端连接,滤波电路的输出端接放大电路的输入端,放大电路的输出端连接信号采集装置的输入端,信号采集装置的输出端和图像重建模块相连。
[0012] 图像重建模块根据信号采集装置输出的生物体的磁声电电压信号重建电导率分布。
[0013] 本发明基于磁声电原理的电导率成像系统的工作过程如下:
[0014] 所述声场驱动激励模块的超声驱动激励源产生脉冲激励信号,作用于超声阵列,超声阵列通过耦合水囊与生物体耦合。超声阵列发射超声波,在生物体组织内产生超声振动,引起生物体组织的局部质点振动。磁场激励模块在生物体组织振动区域产生非均匀静磁场,生物体组织振动的离子在非均匀静磁场作用下受到洛伦兹的作用而产生电荷分离,进而在生物体中形成局部电场,产生电流分布。贴放在生物体的电极测量此电信号,经检测模块的滤波电路、放大电路的滤波和放大,该电压信号由信号采集装置输出至图像重建模块。图像重建模块利用该生物体的磁声电电压信号,以及已知的非均匀静磁分布信息,采用图像重建算法实现电导率分布的重建。
[0015] 所述的图像重建模块采用两种电导率重建算法重建电导率分布,重建算法一为直接代数迭代电导率重建算法、重建算法二为等效均匀场代数迭代电导率重建算法。
[0016] 算法一,直接代数迭代电导率重建算法,包括以下三个步骤:
[0017] 1、利用互易定理建立实际测量过程与互易过程物理量的对应关系
[0018] (1)实际测量过程为:在超声驱动激励模块产生的超声波和磁场激励模块产生的非均匀静磁场作用下,质点的振速为v,非均匀静磁场为B0(r),所述的磁场由磁场激励模块的开放式磁体产生,所述的非均匀静磁场B0(r)强度的分布是已知的。
[0019] (2)互易过程为:关闭声场驱动激励模块,并使磁场激励模块不工作,向检测模块中的电极通入I安培的直流电,此电流在生物体中产生的电流密度为Jr(r)。
[0020] 在磁场激励模块产生的非均匀静磁场B0(r)的激励下,基于互易定理可以得到磁声电电压分布u(r,t)和振动速度势 互易过程电流密度Jr(r),以及磁场激励模块产生的非均匀静磁场B0(r)之间的关系式:
[0021]
[0022] 式(1)中,t表示超声波的传播时间,Ω表示生物体所在的区域,r表示场点,即生物体所在区域Ω中的点, 表示振动速度势,ρ0为生物体的密度,Jr(r)为互易过程中生物体的电流密度,B0(r)为磁场激励模块产生的非均匀静磁场分布,▽·为散度算符。
[0023] 2、根据实际测量过程中测得的磁声电电压u(r,t),依据公式(1)重建互易过程电流密度Jr(r)和静磁场B0(r)之间的关系▽·(Jr(r)×B0(r))
[0024] 式(1)中振动速度势 满足的格林函数为
[0025] 根据公式(1)和格林函数的对称性,可知磁声电电压u(r,t)满足:
[0026]
[0027] 式(2)中,r'为源点,表示超声阵列所在区域的点。
[0028] 利用公式(2)可以得到互易过程的电流密度Jr(r)和开放式磁体产生的非均匀静磁场B0(r)之间的关系式,同时,由于ρ0为常数,因此得到:
[0029]
[0030] 式(3)中,c0表示声场驱动激励模块产生超声波的速度,trd=2T0-t+|r-r'|/c0表示反转场的时间,t表示超声波的传播时间,T0表示反转u(r,t)的时刻,r'表示超声阵列所在区域的点,r表示生物体所在区域Ω中的点,S表示生物体所在区域Ω的表面,n是生物体所在区域Ω边界处的单位法向矢量,u'(r',trd)表示u(r',trd)的一阶导数,u”(r',trd)表示u(r',trd)二阶导数。
[0031] 利用式(3)可以实现变量H(r)=▽·(Jr(r)×B0(r))的重建。
[0032] 3、根据变量H(r)重建电导率分布
[0033] 利用公式(3)实现变量H(r)=▽·(Jr(r)×B0(r))的分布之后,在生物体组织每一个断层面z0上变量H(r)为 可表示为 即:
[0034]
[0035] 其中各个变量括号中(x,y,z0)表示相对应的变量在生物体内部每一个断层面z0上坐标。
[0036] Jr(x,y,z0)为z0断层面上互易过程的电流密度,利用欧姆定律,该电流密度Jr(x,y,z0)可以表示为互易过程电位的梯度和电导率之间的乘积,即:
[0037] Jr(x,y,z)=-σ(x,y,z)▽ur(x,y,z)
[0038] 因此f(x,y,z0)可以表示为:
[0039]
[0040] 其中ur(x,y,z0)是互易过程的电位在z0断层面的分布,σ(x,y,z0)表示生物体的电导率在z0断层面的分布。
[0041] 互易过程的电位在z0断层面的分布满足如下关系式:
[0042]
[0043] 其中的I为互易过程中注入I安培的直流电,rA和rB表示检测模块中电极对的位置,Γ表示目标体所在区域Ω的表面,n表示表面Γ的外法线方向的单位矢量。
[0044] 每一个断层面z0上的 的值表示为f(x,y,z0),由f(x,y,z0)重建电导率σ的方法步骤如下:
[0045] 1)将生物体划分为一系列的子块,认为这些子块内部的电导率均匀,给出生物体的电导率分布矩阵[σ]的初值,通常选择电导率的初值为0.1S/m,给定误差精度ε;
[0046] 2)根据公式(5)计算得到互易过程的电标位ur(x,y,z0);
[0047] 3)利用公式(4)重建z=z0的断层面的电导率分布;
[0048] 4)利用步骤3)得到的各个断层的电导率分布,通过线性差值得到各个子块的电导率分布,如此得到第k次迭代的生物体整个三维区域的电导率分布σk。
[0049] 5)计算第k次迭代得到的生物体电导率分布σk和第k+1次生物体电导率分布σk+1之间的相对误差,比对该相对误差是否满足给定的误差精度ε,如果满足给定的误差精度ε,停止迭代。否则把第k次得到的生物体的电导率分布σk作为初始电导率分布,转到步骤2),以上过程依次迭代下去,直到相邻两次计算得到的生物体电导率分布的相对误差满足精度要求。
[0050] 利用公式(1)和(2)可实现建▽·(Jr(r)×B0(r))的重建,再利用以上步骤1)~步骤5)实现生物体电导率图像的重建,此电导率重建算法称为直接代数迭代电导率重建算法。
[0051] 算法二,等效均匀场代数迭代电导率重建算法如下:
[0052] 等效均匀场代数迭代电导率重建算法基于磁声电电压信号和非均匀静磁场B0(r)之间的关系,以非均匀静磁场B0(r)作为磁场激励源的磁声电电压信号等效为均匀静磁场激励的磁声电电压信号的原理如下公式(6)-(14)。
[0053] 给定三维模型,在非均匀静磁场B0(r)的激励下,生物体中等效电流源对应的电流I(t)和声场驱动激励模块产生的超声波在生物体中的振速v之间的关系为:
[0054] I(t)=∫sσv×B0(r)·dS      (6)
[0055] 其中S表示等效电流源流过的面的面元。
[0056] 根据声学原理,超声动冲量M和振速v需满足:
[0057]
[0058] 其中ρ0为生物体的密度,▽为梯度算符。
[0059] 将式(7)代入式(6),可得:
[0060]
[0061] 其中n表示等效电流源流过的面的面元的法向方向的单位矢量。
[0062] 进一步利用斯托克斯公式和矢量恒等式,式(8)化简为:
[0063]
[0064] 其中l表示面元的外边缘的线,且l的正向与S的外法线方向符合右手定理,S的法向方向为n。
[0065] 考虑实际应用时,超声阵列发出的能量频率范围包含极少的直流频率,表示超声阵列产生的波包的净动量为零,因此式(9)等号右边的第一项为零。
[0066]
[0067] 同时在实际检测中,电极只能检测到生物体的一部分电流,因此信号采集装置采集到的电流只是电流I(t)的一部分,定义采集电压和电流之间的比例为α,检测到的磁声电电压U(t)可表示为:
[0068]
[0069] 式(11)即检测到的磁声电电压U(t)、非均匀静磁场B0(r)、电导率σ以及密度ρ0之间的关系式。
[0070] 实际应用时,▽×B0(r)为一个很小的量,因此式(11)等号右边第二项为一个小量,忽视此小量,故式(11)可简化为:
[0071]
[0072] 公式(12)中,当等效电流源流过的面的面元的法向方向的单位矢量n的方向与的方向相同时,式(12)可表示为:
[0073]
[0074] 其中 取此变量的ex方向分量,则公式(13)可表示为:
[0075]
[0076] 其中Boy和Boz分别为B0(r)沿y方向和z方向的磁场分布θ表示 和B0(r)之间的夹
[0077] 式(14)为非均匀静磁场激励下,磁声电电压信号等效为均匀静磁场激励下的磁声电电压信号的理论关系式,非均匀静磁场激励下,在电导率发生变化的位置磁声电信号会被放大或者缩小 倍。
[0078] 算法二,等效均匀场代数迭代电导率重建算法,流程为:
[0079] 首先利用式(14)将非均匀静磁场激励下采集的磁声电电压信号uin(r,t)调整为等效的均匀场激励下采集的磁声电电压信号uho(r,t);
[0080] 然后将公式(1)-(5)中所有的非均匀静磁场B0(r)替换为均匀静磁场B0;
[0081] 最后利用算法一直接代数迭代电导率重建算法的步骤1)-5),实现电导率的重建。
[0082] 在等效均匀场代数迭代电导率重建算法中,将非均匀静磁场激励下采集的磁声电电压信号uin(r,t)调整为均匀静磁场激励下采集的磁声电电压信号uho(r,t)之后,不仅可以利用直接代数迭代电导率重建算法实现电导率的重建,也可以利用其他均匀静磁场激励下电导率的重建算法获得电导率的图像。附图说明
[0083] 图1本发明电导率成像系统组成示意图;
[0084] 图2本发明电导率成像系统实施例的结构示意图;
[0085] 图3生物体和耦合水囊以及超声阵列之间的位置关系;
[0086] 图中A1超声驱动激励源,A2超声阵列,A3耦合水囊,A4磁场激励模块,A5生物体,A6电极,A7放大电路,A8滤波电路,A9信号采集装置,A10图像重建模块,A11单个超声探头,A12磁场激励模块产生的磁场,A13质点振动速度。
[0087] 图4重建算法一流程图
[0088] 图5重建算法二流程图。

具体实施方式

[0089] 以下结合附图和具体实施方式进一步说明本发明。
[0090] 本发明基于磁声电原理的电导率成像系统包括成像平台和图像重建模块。成像平台与图像重建模块A10连接,将成像平台的采集到的磁声电信号传输到重建算法模块A10。
[0091] 如图1、图2所示,所述的成像平台包括声场驱动激励模块、磁场激励模块A4和检测模块。声场激励源模块生成超声波。磁场激励模块A4产生磁场激励A12作用于生物体,生物体在磁场和声场的作用下产生动生源电流,检测模块检测到此电流信号,图像重建模块A10根据此电压信号重建电导率分布。
[0092] 所述的声场驱动激励模块由超声驱动激励源A1、超声阵列A2及耦合水囊A3组成。超声阵列A2的一端接超声驱动激励源A1,超声阵列A2的另一端和耦合水囊A3接触,耦合水囊A3填充于超声阵列A2和生物体A5之间,与超声阵列A2和生物体A5接触良好,如图3所示,超声驱动激励源A1激励超声阵列A2产生的超声波能够传播到生物体A5中。
[0093] 磁场激励模块A4产生的非均匀静磁场A12由开放式磁体产生。
[0094] 检测模块由电极A6、滤波电路A8、放大电路A7及信号采集装置A9组成。电极A6和生物体A5接触,检测生物体表面的电压信号,检测的电压信号经滤波电路A8和放大电路A7滤波和放大,由信号采集装置A9采集。所述的电极A6的一端和生物体A5连接,电极A6的另一端和滤波电路A8的输入端连接,滤波电路A8的输出端连接放大电路A7的输入端,放大电路A7的输出端连接信号采集装置A9的输入端,信号采集装置A9的输出端和图像重建模块A10相连。
[0095] 图像重建模块A10根据信号采集装置A9输出的生物体A5的电压信号重建电导率分布。
[0096] 本发明的基于磁声电原理的电导率成像系统的工作过程如下:
[0097] 所述声场驱动激励模块的超声驱动激励源A1产生脉冲激励信号,作用于超声阵列A2,超声阵列A2通过耦合水囊A3与生物体A5耦合。超声阵列A2发射超声波,在生物体A5组织内产生超声振动,引起生物体组织的局部质点振动。磁场激励模块A4在生物体A5组织振动区域产生非均匀静磁场A12,生物体A5组织振动的离子在非均匀静磁场A12作用下受到洛伦兹力的作用而产生电荷分离,进而在生物体A5中形成局部电场,产生电流分布。贴放在生物体A5的电极A6测量此电信号,经检测模块的滤波电路A8、放大电路A7的滤波和放大,该电压信号由信号采集装置A9输出至图像重建模块A10。图像重建模块A10利用该生物体A5电压信号,以及已知的非均匀静磁A12分布信息,采用图像重建算法实现电导率分布重建。所述的图像重建模块A10采用两种电导率重建算法重建电导率分布,重建算法一为直接代数迭代电导率重建算法、重建算法二为等效均匀场代数迭代电导率重建算法。
[0098] 重建算法一的流程如图4所示,步骤如下:
[0099] 1)利用u(r,t)和公式(2)重建得到变量H(r)的分布;
[0100] 2)将生物体划分为一系列的子块,认为这些子块内部的电导率均匀,给出生物体的电导率的分布矩阵[σ]的初值,通常选择电导率的初值为0.1S/m,给定误差精度ε;
[0101] 3)根据公式(5)计算得到互易过程的电标位ur(x,y,z0);
[0102] 4)利用公式(4)重建z=z0的断层面的电导率分布;
[0103] 5)利用步骤3)得到的各个断层的电导率分布,通过线性差值得到各个子块的电导率分布,如此得到第k次迭代的生物体整个三维区域的电导率分布σk;
[0104] 6)计算第k次迭代得到的生物体电导率分布σk和第k+1次生物体电导率分布σk+1之间的相对误差,比对该相对误差是否满足给定的误差精度ε,如果满足给定的误差精度ε,停止迭代。否则把第k次得到的生物体的电导率分布σk作为初始电导率分布,转到步骤2),以上过程依次迭代下去,直到相邻两次计算得到的生物体电导率分布的相对误差满足精度要求。
[0105] 重建算法二的流程如图5所示,步骤如下:
[0106] 1)利用式(14)首先将非均匀静磁场B0(r)激励下采集的磁声电电压信号uin(r,t)调整为等效的均匀场B0激励下采集的磁声电电压信号uho(r,t);
[0107] 2)然后将公式(1)-(5)中所有的非均匀静磁场B0(r)替换为均匀静磁场B0;
[0108] 3)最后利用算法一的步骤1)-6),实现电导率的重建。
[0109] 在等效均匀场代数迭代电导率重建算法中,将非均匀静磁场B0(r)激励下采集的磁声电电压信号uin(r,t)调整为均匀静磁场B0激励下采集的磁声电电压信号uho(r,t)之后,不仅可以利用直接代数迭代电导率重建算法实现电导率的重建,也可以利用其他均匀静磁场激励下电导率的重建算法获得电导率的图像。
高效检索全球专利

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

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

申请试用

分析报告

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

申请试用

QQ群二维码
意见反馈