CN107677976B - 一种自适应核磁共振梯度预加重波形产生装置及方法 - Google Patents
一种自适应核磁共振梯度预加重波形产生装置及方法 Download PDFInfo
- Publication number
- CN107677976B CN107677976B CN201710884240.6A CN201710884240A CN107677976B CN 107677976 B CN107677976 B CN 107677976B CN 201710884240 A CN201710884240 A CN 201710884240A CN 107677976 B CN107677976 B CN 107677976B
- Authority
- CN
- China
- Prior art keywords
- gradient
- data
- preemphasis
- field waveform
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种自适应核磁共振梯度预加重波形产生装置,包括用户计算机、采样计算机、梯度控制模块、梯度功率放大器模块和梯度信号采集模块;还公开了一种自适应核磁共振梯度预加重波形产生方法,本发明基于梯度电流预加重补偿原理,首先在不进行预加重处理时,通过梯度功率放大器模块的电流监控端口采集梯度线圈上失真的实际梯度磁场波形数据,上传到用户计算机后计算出预设梯度磁场波形数据与实际梯度磁场波形数据的偏差,根据此偏差获得涡流曲线,然后通过最小二乘法拟合出一组幅度常数和时间常数,再通过拟合出来的幅度常数和时间常数对梯度预加重参数进行预加重处理,以此来抵消涡流对梯度磁场的影响。
Description
技术领域
本发明涉及核磁共振仪器技术领域,具体涉及一种自适应核磁共振梯度预加重波形产生装置,还涉及一种自适应核磁共振梯度预加重波形产生方法,用于减小核磁共振实验中由于梯度涡流引起的梯度磁场畸变和失真。
背景技术
在核磁共振成像系统中,梯度系统负责产生线性变化的梯度磁场,用于对被测体进行空间编码,梯度系统通过在梯度线圈上施加梯度电流来产生所需要的梯度磁场。由于梯度线圈处于磁体内部,而磁体内包含很多金属部件,如磁屏蔽外壳、超导线圈、匀场线圈和射频线圈等,那么根据电磁感应定律,在梯度电流进行快速开关的突变过程中,会在这些金属部件内部感应出阻碍磁场变化的电流,这些感应出来的暂态电流称之为涡流(EddyCurrent)。涡流的存在会严重影响梯度磁场的变化,使磁场产生畸变和失真,导致梯度编码出现错误和偏差,如果不进行补偿,则会在磁共振图像中形成伪影,严重影响成像质量。因此,磁共振成像中必须寻找有效的方法来消除涡流对梯度磁场的影响。
为了消除梯度涡流的影响,目前采用的最普遍方法是对梯度电流进行预加重处理,即对梯度电流进行预补偿来改善梯度电流质量,从而产生期望的梯度磁场输出效果。这种方法主要是根据滤波器的瞬态响应进行设计,在梯度波形进入梯度功率放大器之前,由预加重滤波器给原始的梯度电流叠加上补偿电流,使梯度波形预失真,然后将预失真的波形输出到梯度功率放大器,最后在梯度线圈上得到满足要求的梯度波形。
补偿电流的大小满足多组e指数函数的叠加形式,方向与涡流方向相反。其大小可以表示为:其中Ie(t)为叠加后的补偿电流,N为补偿滤波器的级数,Ai和τi分别为补偿电流曲线的幅度常数和时间常数。因此,使用梯度预加重的方法进行涡流补偿即变为寻找一组或多组合适的幅度常数Ai和时间常数τi,然后用这些参数产生补偿电流叠加到原始梯度上,使得补偿电流刚好与涡流的变化趋势相反,抵消涡流对梯度磁场的影响,从而达到补偿涡流的目的。
目前常用的寻找补偿电流幅度常数Ai和时间常数τi的方法是基于核磁共振信号强度的梯度预加重调节方法,使用的脉冲序列示意图如图1所示。
该方法通常先对被测梯度线圈施加一段时间的驱动电流后关闭,此后的一段时间内信号都会受涡流磁场的影响。依次设置四组射频激发(RF)和采样(ACQ)得到四组自由感应衰减(FID)信号,其中最后一组激发采样设置的延时比较长,确保其不受涡流的影响,作为补偿调节的参照。根据梯度涡流的e指数衰减规律,启动扫描后可以观察到越靠近梯度脉冲的信号离散越快,信号幅度越小,后面采集的信号由于涡流影响逐渐减小,幅度逐渐增大。接着进行涡流补偿调节,假设有三组不同时间常数和幅度常数的涡流存在系统中,先对时间最长的一组涡流进行补偿,设置一组时间常数较长的补偿电流参数,将第三个FID信号波形恢复到与第四个FID相同。然后再设置时间常数中等和较小的两组补偿电流参数,依次补偿第二个FID和第一个FID,直到4个FID信号形状与幅度都基本一致,最后得到的三组时间常数和幅度常数则作为该梯度的涡流补偿参数。
上述基于核磁共振信号强度的梯度预加重调节方法虽然可以获得涡流补偿参数,但是依然存在以下不足:1.需要调试人员手动反复设置梯度预加重参数才能获得较好补偿效果,调试过程费时费力;2.时间常数和幅度常数的设置依赖于专业调试人员的经验,普通用户难以完成调试;3.该方法只能获取梯度下降沿产生的补偿电流常数,而梯度上升沿的补偿参数则无法获取,不能精确补偿涡流的影响。4.由于需要通过观察核磁共振信号强度来调节预加重参数,因此调试结果极易受到磁场不均匀性、样品弛豫时间、脉冲序列参数等因素的影响。
发明内容
针对现有技术中存在的上述不足之处,本发明提出了一种自适应核磁共振梯度预加重波形产生装置,还提供一种自适应核磁共振梯度预加重波形产生方法。
本发明的上述目的通过以下技术方案实现:
一种自适应核磁共振梯度预加重波形产生装置,包括:
用户计算机,用于将梯度预加重参数、原始梯度数据和脉冲序列控制参数发送到采样计算机;还用于根据接收到实际梯度磁场波形数据生成涡流补偿参数,其中,涡流补偿参数包括幅度常数Ai和时间常数τi;还用于将涡流补偿参数设置到梯度预加重参数;
采样计算机,用于将梯度预加重参数、原始梯度数据和脉冲序列控制参数分别对应转换为梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据,并将梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据分别发送到梯度控制模块和梯度信号采集模块;还用于接收梯度信号采集模块发送的实际梯度磁场波形数据,并将实际梯度磁场波形数据发送到用户计算机;
梯度控制模块,用于将梯度预加重二进制数据叠加到原始梯度二进制数据上,得到预加重之后的梯度数据,然后将预加重之后的梯度数据转换成模拟梯度小信号,并将模拟梯度小信号输出到梯度功率放大器模块;
梯度功率放大器模块,用于将模拟梯度小信号进行放大之后驱动梯度线圈产生梯度磁场,梯度功率放大器模块包括用于监控梯度线圈上实际梯度磁场波形的电流监控接口;
梯度信号采集模块,用于根据接收到的脉冲序列二进制数据采集梯度功率放大器模块的电流监控端口上的实际梯度磁场波形数据,并将实际梯度磁场波形数据上传到采样计算机。
如上所述的幅度常数Ai和时间常数τi使得涡流曲线ΔG(n)与理论上的涡流曲线表达式之间的方差E达到最小:涡流曲线ΔG(n)为经过离散化处理的预设梯度磁场波形数据与实际梯度磁场波形数据的幅度差,n∈{1~M},M为取样点数,N为预加重补偿滤波器的级数。
一种自适应核磁共振梯度预加重波形产生方法,包括以下步骤:
步骤1、用户计算机首先将梯度预加重参数设置为初始值,用户计算机将初始的幅度常数Ai和时间常数τi设置到梯度预加重参数中;
步骤2、用户计算机将梯度预加重参数、原始梯度数据和脉冲序列控制参数发送到采样计算机;
步骤3、采样计算机接收到用户计算机发送的梯度预加重参数、原始梯度数据和脉冲序列控制参数后,分别对应转换为梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据,并将梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据分别发送到梯度控制模块和梯度信号采集模块;
步骤4、梯度控制模块将梯度预加重二进制数据叠加到原始梯度二进制数据上,得到预加重之后的梯度数据,然后将预加重之后的梯度数据转换成模拟梯度小信号,模拟梯度小信号作为梯度控制模块的输出信号,输出到梯度功率放大器模块;
步骤5、模拟梯度小信号经过梯度功率放大器模块放大之后驱动梯度线圈产生梯度磁场,梯度信号采集模块根据接收到的脉冲序列二进制数据开始采集梯度功率放大器模块的电流监控端口上的实际梯度磁场波形数据,并将实际梯度磁场波形数据上传到采样计算机;
步骤6、采样计算机接收到实际梯度磁场波形数据后,将实际梯度磁场波形数据上传到用户计算机,用户计算机接收到实际梯度磁场波形数据后,求出涡流曲线ΔG(n),
如果涡流曲线ΔG(n)小于规定的阈值,将幅度常数Ai和时间常数τi作为涡流补偿参数进行保存,结束调试过程;
如果涡流曲线ΔG(n)大于或等于规定的阈值,则求取幅度常数Ai和时间常数τi,使得涡流曲线ΔG(n)与理论上的涡流曲线表达式之间的方差E达到最小,并将幅度常数Ai和时间常数τi设置到梯度预加重参数中返回步骤2,直到涡流曲线ΔG(n)的值小于规定的阈值,其中,n∈{1~M},M为取样点数,N为预加重补偿滤波器的级数。
如上所述的步骤6包括以下步骤:
步骤6.1、用户计算机首先计算经过离散化处理的预设梯度磁场波形数据与实际梯度磁场波形数据的幅度差作为涡流曲线ΔG(n),计算公式为ΔG(n)=G0(n)-G(n),其中G0(n)为预设梯度磁场波形数据,G(n)为实际梯度磁场波形数据,n∈{1~M},M为取样点数;
如果涡流曲线ΔG(n)小于规定的阈值,将幅度常数Ai和时间常数τi作为涡流补偿参数进行保存,结束调试过程;
如果涡流曲线ΔG(n)大于或等于规定的阈值,则进入步骤6.2,直到涡流曲线ΔG(n)的值小于规定的阈值。
步骤6.2、通过最小二乘法准则进行拟合,求出一组幅度常数Ai和时间常数τi,使涡流曲线ΔG(n)与理论上的涡流曲线表达式之间的方差E达到最小,拟合计算的公式为:
上式中M为取样点数,N为预加重补偿滤波器的级数,T1和T2分别为设定时间范围的最小值和最大值,
步骤6.3、将步骤6.2获得的幅度常数Ai和时间常数τi通过用户计算机1设置到梯度预加重参数中,返回步骤2。
如上所述的幅度常数Ai的取值范围为[-1,1],时间常数τi的单位为ms,最小值T1和最大值T2分别为1ms和1000ms。
本发明基于梯度电流预加重补偿原理,首先在不进行预加重处理时,通过梯度功率放大器的电流监控端口采集梯度线圈上失真的实际梯度磁场波形数据,上传到用户计算机后,计算出预设梯度波形与实际梯度磁场波形的幅度差,根据幅度差获得涡流曲线数据,然后用曲线拟合的方法求出预加重补偿电流对应的多组幅度常数Ai和时间常数τi,再将Ai和τi发送到梯度控制模块对梯度波形进行预加重处理,以此来抵消涡流对梯度磁场的影响。梯度涡流的补偿效果通过再次采集预加重之后的实际梯度磁场波形数据来衡量,如果预加重之后采集的实际梯度磁场波形数据与预设梯度波形之间的幅度差小于预先设定的阈值,则补偿结束,否则继续使用上述方法求出预加重参数进一步对涡流进行补偿,直到幅度差小于设定的阈值后结束。本发明可以实现自动测试和拟合预加重参数的功能,能够快速完成梯度涡流的补偿,并能分别对梯度上升沿和下降沿计算补偿参数,实现更精确的涡流补偿,且不受磁场不均匀性、样品弛豫时间和脉冲序列参数等因素的影响,具有很强的抗干扰性,有效克服了现有涡流补偿方法的缺点和不足。
本发明与现有技术相比,具有以下优点:
1、可以实现计算机自动测试和拟合预加重参数的功能,快速完成梯度涡流的补偿,补偿过程不依赖于专业人员的调试经验;
2、能够分别对梯度上升沿和下降沿引起的涡流进行补偿;
3、可以同时对X、Y、Z三个方向上的梯度涡流进行测量和计算补偿参数,极大提高了涡流补偿的效率;
4、本发明采用直接采集梯度磁场波形的方法获得涡流曲线,调试结果与磁场不均匀性、样品弛豫时间和脉冲序列参数等因素无关,具有很强的抗干扰性;
附图说明
图1为基于核磁共振信号强度的梯度预加重调节方法的脉冲序列示意图;
图2为本发明的原理示意图;
图3为本发明工作流程图;
图4为根据预设梯度波形和实际梯度波形计算涡流曲线的示意图,其中(a)为G0(n)波形示意图,(b)为G(n)波形示意图,(c)为ΔG(n)波形示意图。
具体实施方式
以下结合附图对本发明的技术方案作进一步详细描述。
实施例1
如图2所示,一种自适应核磁共振梯度预加重波形产生装置,包括:
用户计算机1,用于将梯度预加重参数、原始梯度数据和脉冲序列控制参数发送到采样计算机2;还用于根据接收到实际梯度磁场波形数据生成涡流补偿参数,其中,涡流补偿参数包括幅度常数Ai和时间常数τi;还用于将涡流补偿参数设置到梯度预加重参数。
幅度常数Ai和时间常数τi由如下方式生成:
用户计算机1首先计算经过离散化处理的预设梯度磁场波形数据与实际梯度磁场波形数据的幅度差ΔG(n),计算公式为ΔG(n)=G0(n)-G(n),其中G0(n)为预设梯度磁场波形数据,G(n)为实际梯度磁场波形数据,计算得到的ΔG(n)即为涡流曲线,n∈{1~M},M为取样点数;
求取幅度常数Ai和时间常数τi,使得涡流曲线ΔG(n)与理论上的涡流曲线表达式之间的方差E达到最小。
拟合计算的公式为:
上式中M为取样点数,N为预加重补偿滤波器的级数,Ai和τi分别为补偿电流曲线的幅度常数和时间常数,同时对所求幅度常数Ai和时间常数τi的取值范围进行限定,其中幅度常数Ai的取值范围为[-1,1],时间常数τi的单位一般为ms,最小值T1和最大值T2分别为1ms和1000ms,幅度常数Ai和时间常数τi的求解变成一个有约束条件的多变量函数最优化问题,可以利用多变量函数寻优的方法求得数值解,从而得到涡流补偿参数Ai和τi;
用户计算机1可以是通用台式计算机、工作站或移动笔记本电脑等提供图形用户界面的设备,用户计算机1与采样计算机2之间通过以太网进行数据通讯。
采样计算机2,用于将梯度预加重参数、原始梯度数据和脉冲序列控制参数分别对应转换为梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据,并将梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据分别发送到梯度控制模块3和梯度信号采集模块5;还用于接收梯度信号采集模块5发送的实际梯度磁场波形数据,并将实际梯度磁场波形数据发送到用户计算机1。
所述的梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据位宽可以为48bits,每2个48bits的二进制数据为一组,每一组数据控制一个脉冲序列事件。采样计算机2是配备PCIe高速数据传输接口的通用台式计算机,采样计算机2与梯度控制模块3和梯度信号采集模块5之间通过PCIe总线进行高速数据通讯。
梯度控制模块3,用于将梯度预加重二进制数据叠加到原始梯度二进制数据上,得到预加重之后的梯度数据,然后将预加重之后的梯度数据转换成模拟梯度小信号,并将模拟梯度小信号输出到梯度功率放大器模块4。
梯度控制模块3可以由FPGA控制模块和DAC数据转换模块实现,FPGA控制模块的核心器件为Xilinx xc7a100tfgg676pkg,DAC数据转换模块的核心器件为AD5764CSU。
梯度功率放大器模块4,用于将模拟梯度小信号进行放大之后驱动梯度线圈产生梯度磁场。
梯度功率放大器模块4上还提供了用于监控梯度线圈上实际梯度磁场波形的电流监控接口401,该接口用于采集实际梯度磁场波形数据。梯度功率放大器模块4的核心部件可以选用CopleyMODEL 265P。
梯度信号采集模块5,用于根据接收到的脉冲序列二进制数据采集梯度功率放大器模块4的电流监控端口401上的实际梯度磁场波形数据,并将实际梯度磁场波形数据上传到采样计算机2。
梯度信号采集模块5接收并解析从采样计算机2发送过来的脉冲序列二进制数据,用于控制实际梯度磁场波形数据的采样时序。梯度信号采集模块5可以由FPGA控制模块和ADC数据转换模块实现,FPGA控制模块的核心器件为Xilinx xc7a100tfgg676pkg,ADC数据转换模块的核心器件为AD9244BST-65。
一种自适应核磁共振梯度预加重波形产生方法,其工作流程如图3所示,具体工作流程如下:
步骤1、用户计算机1首先将梯度预加重参数设置为初始值,用户计算机1初始的幅度常数Ai和时间常数τi设置到梯度预加重参数中;
步骤2、用户计算机1将梯度预加重参数、原始梯度数据和脉冲序列控制参数发送到采样计算机2;
步骤3、采样计算机2接收到用户计算机1发送的梯度预加重参数、原始梯度数据和脉冲序列控制参数后,分别对应转换为梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据,并将梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据分别发送到梯度控制模块3和梯度信号采集模块5;
步骤4、梯度控制模块3将梯度预加重二进制数据叠加到原始梯度二进制数据上,得到预加重之后的梯度数据,然后将预加重之后的梯度数据转换成模拟梯度小信号,模拟梯度小信号作为梯度控制模块3的输出信号,输出到梯度功率放大器模块4;
步骤5、模拟梯度小信号经过梯度功率放大器模块4放大之后驱动梯度线圈产生梯度磁场,梯度信号采集模块5根据接收到的脉冲序列二进制数据开始采集梯度功率放大器模块4的电流监控端口401上的实际梯度磁场波形数据,并将实际梯度磁场波形数据上传到采样计算机2;
步骤6、采样计算机2接收到实际梯度磁场波形数据后,将实际梯度磁场波形数据上传到用户计算机1,用户计算机1接收到实际梯度磁场波形数据后,通过如下所述方法求出涡流补偿参数:
步骤6.1、用户计算机1首先计算经过离散化处理的预设梯度磁场波形数据与实际梯度磁场波形数据的幅度差作为涡流曲线ΔG(n),计算公式为ΔG(n)=G0(n)-G(n),其中G0(n)为预设梯度磁场波形数据,G(n)为实际梯度磁场波形数据,其中,n∈{1~M},M为取样点数;
如果涡流曲线ΔG(n)小于规定的阈值,便可认为补偿了涡流的影响,将幅度常数Ai和时间常数τi作为涡流补偿参数进行保存,结束调试过程;
如果涡流曲线ΔG(n)大于或等于规定的阈值,则进入步骤6.2,直到涡流曲线ΔG(n)的值小于规定的阈值。
步骤6.2、涡流曲线的示意图如图4所示,在本实施例中,预设梯度磁场波形数据G0(n)的波形如图4(a)所示,实际梯度磁场波形数据G(n)的波形如图4(b)所示,则计算二者的幅度差得到的涡流曲线ΔG(n)如图4(c)所示,涡流曲线ΔG(n)由I1(t)和I2(t)两段曲线组成,其中t1t2段为梯度上升沿的涡流曲线I1(t),t3t4段为梯度下降沿的涡流曲线I2(t)。通过最小二乘法准则进行拟合,求出一组幅度常数Ai和时间常数τi,使实际测量得到的涡流曲线ΔG(n)与理论上的涡流曲线表达式之间的方差E达到最小,拟合计算的公式为:
上式中M为取样点数,N为预加重补偿滤波器的级数,Ai和τi分别为补偿电流曲线的幅度常数和时间常数,同时对所求幅度常数Ai和时间常数τi的取值范围进行限定,其中幅度常数Ai的取值范围为[-1,1],时间常数τi的单位一般为ms,T1和T2分别为设定时间范围的最小值和最大值,优选的最小值T1和最大值T2分别为1ms和1000ms,幅度常数Ai和时间常数τi的求解变成一个有约束条件的多变量函数最优化问题,可以利用多变量函数寻优的方法求得数值解,从而得到幅度常数Ai和时间常数τi;
步骤6.3、根据理论上的涡流曲线表达式幅度常数Ai和时间常数τi就是获得补偿电流曲线所需要设置的涡流补偿参数。将步骤6.2获得的幅度常数Ai和时间常数τi通过用户计算机1设置到梯度预加重参数中,返回步骤2。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (4)
1.一种自适应核磁共振梯度预加重波形产生装置,其特征在于,包括:
用户计算机(1),用于将梯度预加重参数、原始梯度数据和脉冲序列控制参数发送到采样计算机(2);还用于根据接收到实际梯度磁场波形数据生成涡流补偿参数,其中,涡流补偿参数包括幅度常数Ai和时间常数τi;还用于将涡流补偿参数设置到梯度预加重参数;
采样计算机(2),用于将梯度预加重参数、原始梯度数据和脉冲序列控制参数分别对应转换为梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据,并将梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据分别发送到梯度控制模块(3)和梯度信号采集模块(5);还用于接收梯度信号采集模块(5)发送的实际梯度磁场波形数据,并将实际梯度磁场波形数据发送到用户计算机(1);
梯度控制模块(3),用于将梯度预加重二进制数据叠加到原始梯度二进制数据上,得到预加重之后的梯度数据,然后将预加重之后的梯度数据转换成模拟梯度小信号,并将模拟梯度小信号输出到梯度功率放大器模块(4);
梯度功率放大器模块(4),用于将模拟梯度小信号进行放大之后驱动梯度线圈产生梯度磁场,梯度功率放大器模块(4)包括用于监控梯度线圈上实际梯度磁场波形的电流监控接口(401);
梯度信号采集模块(5),用于根据接收到的脉冲序列二进制数据采集梯度功率放大器模块(4)的电流监控端口(401)上的实际梯度磁场波形数据,并将实际梯度磁场波形数据上传到采样计算机(2),
所述的幅度常数Ai和时间常数τi使得涡流曲线ΔG(n)与理论上的涡流曲线表达式之间的方差E达到最小:涡流曲线ΔG(n)为经过离散化处理的预设梯度磁场波形数据与实际梯度磁场波形数据的幅度差,n∈{1~M},M为取样点数,N为预加重补偿滤波器的级数。
2.一种自适应核磁共振梯度预加重波形产生方法,其特征在于,包括以下步骤:
步骤1、用户计算机(1)首先将梯度预加重参数设置为初始值,用户计算机(1)将初始的幅度常数Ai和时间常数τi设置到梯度预加重参数中;
步骤2、用户计算机(1)将梯度预加重参数、原始梯度数据和脉冲序列控制参数发送到采样计算机(2);
步骤3、采样计算机(2)接收到用户计算机(1)发送的梯度预加重参数、原始梯度数据和脉冲序列控制参数后,分别对应转换为梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据,并将梯度预加重二进制数据、原始梯度二进制数据和脉冲序列二进制数据分别发送到梯度控制模块(3)和梯度信号采集模块(5);
步骤4、梯度控制模块(3)将梯度预加重二进制数据叠加到原始梯度二进制数据上,得到预加重之后的梯度数据,然后将预加重之后的梯度数据转换成模拟梯度小信号,模拟梯度小信号作为梯度控制模块(3)的输出信号,输出到梯度功率放大器模块(4);
步骤5、模拟梯度小信号经过梯度功率放大器模块(4)放大之后驱动梯度线圈产生梯度磁场,梯度信号采集模块(5)根据接收到的脉冲序列二进制数据开始采集梯度功率放大器模块(4)的电流监控端口(401)上的实际梯度磁场波形数据,并将实际梯度磁场波形数据上传到采样计算机(2);
步骤6、采样计算机(2)接收到实际梯度磁场波形数据后,将实际梯度磁场波形数据上传到用户计算机(1),用户计算机(1)接收到实际梯度磁场波形数据后,求出涡流曲线ΔG(n),
如果涡流曲线ΔG(n)小于规定的阈值,将幅度常数Ai和时间常数τi作为涡流补偿参数进行保存,结束调试过程;
如果涡流曲线ΔG(n)大于或等于规定的阈值,则求取幅度常数Ai和时间常数τi,使得涡流曲线ΔG(n)与理论上的涡流曲线表达式之间的方差E达到最小,并将幅度常数Ai和时间常数τi设置到梯度预加重参数中返回步骤2,直到涡流曲线ΔG(n)的值小于规定的阈值,其中,n∈{1~M},M为取样点数,N为预加重补偿滤波器的级数。
3.根据权利要求2所述的一种自适应核磁共振梯度预加重波形产生方法,其特征在于,所述的步骤6包括以下步骤:
步骤6.1、用户计算机(1)首先计算经过离散化处理的预设梯度磁场波形数据与实际梯度磁场波形数据的幅度差作为涡流曲线ΔG(n),计算公式为ΔG(n)=G0(n)-G(n),其中G0(n)为预设梯度磁场波形数据,G(n)为实际梯度磁场波形数据,n∈{1~M},M为取样点数;
如果涡流曲线ΔG(n)小于规定的阈值,将幅度常数Ai和时间常数τi作为涡流补偿参数进行保存,结束调试过程;
如果涡流曲线ΔG(n)大于或等于规定的阈值,则进入步骤6.2,直到涡流曲线ΔG(n)的值小于规定的阈值,
步骤6.2、通过最小二乘法准则进行拟合,求出一组幅度常数Ai和时间常数τi,使涡流曲线ΔG(n)与理论上的涡流曲线表达式之间的方差E达到最小,拟合计算的公式为:
上式中M为取样点数,N为预加重补偿滤波器的级数,T1和T2分别为设定时间范围的最小值和最大值,
步骤6.3、将步骤6.2获得的幅度常数Ai和时间常数τi通过用户计算机1设置到梯度预加重参数中,返回步骤2。
4.根据权利要求3所述的一种自适应核磁共振梯度预加重波形产生方法,其特征在于,所述的幅度常数Ai的取值范围为[-1,1],时间常数τi的单位为ms,最小值T1和最大值T2分别为1ms和1000ms。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710884240.6A CN107677976B (zh) | 2017-09-26 | 2017-09-26 | 一种自适应核磁共振梯度预加重波形产生装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710884240.6A CN107677976B (zh) | 2017-09-26 | 2017-09-26 | 一种自适应核磁共振梯度预加重波形产生装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107677976A CN107677976A (zh) | 2018-02-09 |
CN107677976B true CN107677976B (zh) | 2019-07-16 |
Family
ID=61138119
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710884240.6A Active CN107677976B (zh) | 2017-09-26 | 2017-09-26 | 一种自适应核磁共振梯度预加重波形产生装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107677976B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108303665B (zh) * | 2018-02-27 | 2019-12-13 | 奥泰医疗系统有限责任公司 | 一种磁共振成像系统涡流补偿方法 |
CN108872900A (zh) * | 2018-07-02 | 2018-11-23 | 中国科学院武汉物理与数学研究所 | 一种时分复用核磁共振梯度数字预加重装置 |
EP3617731A1 (en) * | 2018-08-29 | 2020-03-04 | Erasmus University Medical Center Rotterdam | Magnetic resonance imaging based on transient response signals |
CN109765512B (zh) * | 2019-01-18 | 2021-09-17 | 上海联影医疗科技股份有限公司 | 磁共振梯度系统及其涡流补偿方法及装置 |
DE102019215046B4 (de) | 2019-09-30 | 2023-09-21 | Siemens Healthcare Gmbh | Verfahren zur Kompensation von Wirbelströmen bei einer Erfassung von Messdaten mittels Magnetresonanz |
CN110954853B (zh) * | 2019-12-03 | 2021-10-12 | 华东师范大学 | 一种梯度预补偿波形数字分辨率提高方法及装置 |
CN111521963B (zh) * | 2020-05-06 | 2021-02-19 | 电子科技大学 | 一种梯度波形速率可自动适配输出的方法及磁共振系统 |
CN112881959B (zh) * | 2020-06-10 | 2021-12-21 | 苏州纽迈分析仪器股份有限公司 | 一种用于磁共振成像的梯度涡流补偿方法及系统 |
CN115276676B (zh) * | 2021-04-30 | 2024-06-04 | 瑞昱半导体股份有限公司 | 发送器电路、补偿值校正装置与同相与正交不平衡补偿值校正方法 |
CN113296037A (zh) * | 2021-05-21 | 2021-08-24 | 电子科技大学 | 一种高场磁共振梯度控制器 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6025715A (en) * | 1997-12-15 | 2000-02-15 | King; Kevin F. | Method for compensating an MRI system for time dependent gradient field distortion |
CN102103196A (zh) * | 2009-12-18 | 2011-06-22 | 东软飞利浦医疗设备系统有限责任公司 | 具有可调数字预加重的磁共振梯度波形产生系统及方法 |
CN102298129A (zh) * | 2011-05-30 | 2011-12-28 | 苏州安科医疗系统有限公司 | 一种用于核磁共振成像系统的涡流测量及补偿方法 |
CN104181479A (zh) * | 2013-05-23 | 2014-12-03 | 上海联影医疗科技有限公司 | 用于磁共振成像系统的涡流补偿方法 |
CN104614694A (zh) * | 2015-01-27 | 2015-05-13 | 华东师范大学 | 一种磁共振梯度涡流补偿方法 |
CN106291423A (zh) * | 2016-09-07 | 2017-01-04 | 厦门大学 | 核磁共振仪梯度预加重调节装置 |
CN106772134A (zh) * | 2017-03-01 | 2017-05-31 | 中国科学院武汉物理与数学研究所 | 一种自动磁场补偿的装置与方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7783441B2 (en) * | 2003-04-17 | 2010-08-24 | Northern Digital Inc. | Eddy current detection and compensation |
JP5416960B2 (ja) * | 2008-12-17 | 2014-02-12 | 株式会社東芝 | 磁気共鳴イメージング装置 |
-
2017
- 2017-09-26 CN CN201710884240.6A patent/CN107677976B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6025715A (en) * | 1997-12-15 | 2000-02-15 | King; Kevin F. | Method for compensating an MRI system for time dependent gradient field distortion |
CN102103196A (zh) * | 2009-12-18 | 2011-06-22 | 东软飞利浦医疗设备系统有限责任公司 | 具有可调数字预加重的磁共振梯度波形产生系统及方法 |
CN102298129A (zh) * | 2011-05-30 | 2011-12-28 | 苏州安科医疗系统有限公司 | 一种用于核磁共振成像系统的涡流测量及补偿方法 |
CN104181479A (zh) * | 2013-05-23 | 2014-12-03 | 上海联影医疗科技有限公司 | 用于磁共振成像系统的涡流补偿方法 |
CN104614694A (zh) * | 2015-01-27 | 2015-05-13 | 华东师范大学 | 一种磁共振梯度涡流补偿方法 |
CN106291423A (zh) * | 2016-09-07 | 2017-01-04 | 厦门大学 | 核磁共振仪梯度预加重调节装置 |
CN106772134A (zh) * | 2017-03-01 | 2017-05-31 | 中国科学院武汉物理与数学研究所 | 一种自动磁场补偿的装置与方法 |
Non-Patent Citations (1)
Title |
---|
一种具有数字预加重的磁共振成像梯度波形发生器;徐勤 等;《波谱学杂志》;20060331;第23卷(第1期);第11-16页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107677976A (zh) | 2018-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107677976B (zh) | 一种自适应核磁共振梯度预加重波形产生装置及方法 | |
US11698430B2 (en) | Eddy current mitigation systems and methods | |
CN108303665B (zh) | 一种磁共振成像系统涡流补偿方法 | |
CN102298129B (zh) | 一种用于核磁共振成像系统的涡流测量及补偿方法 | |
Faifer et al. | Characterization of voltage instrument transformers under nonsinusoidal conditions based on the best linear approximation | |
CN108196215B (zh) | 一种直流电子式互感器延时时间测试方法及装置 | |
CN104375105B (zh) | 用于线圈的自动高频匀场配置 | |
Serrano et al. | Why magnet: Quantifying the complexity of modeling power magnetic material characteristics | |
CN104808251A (zh) | 一种提高Overhauser磁力仪拉莫尔信号测频精度的方法及其电路 | |
JPH11235325A (ja) | Mrシステムの磁界補償方法および装置 | |
DE102010061568A1 (de) | System und Verfahren zur Modellierung der durch den Gradientenspulenbetrieb hervorgerufenen Magnetfelddrift | |
CN109765512A (zh) | 磁共振梯度系统及其涡流补偿方法及装置 | |
CN104614659A (zh) | 自动化测试系统和方法 | |
CN100533168C (zh) | 确定对涡流场进行补偿设置的方法 | |
CN110515020A (zh) | 一种噪声系数分析仪接收通道最佳线性增益校准补偿方法 | |
CN104614694A (zh) | 一种磁共振梯度涡流补偿方法 | |
CN106526520B (zh) | 一种光学电流互感器光源波动干扰试验系统 | |
CN104635187B (zh) | 磁共振系统的自适应有源自动匀场方法和系统 | |
CN105301545B (zh) | 电磁兼容性测试系统的一种校准方法 | |
CN203849435U (zh) | 一种检测核磁共振找水仪器的微弱信号检测平台 | |
CN106680756A (zh) | 直流电压互感器准确度等级测试方法、系统及装置 | |
CN103885015A (zh) | 对磁共振扫描序列进行梯度延时补偿的方法、装置和系统 | |
Meng et al. | NMR log data de-noising method based on a variable order wavelet packet domain adaptive filtering | |
CN104572388B (zh) | 数据滤波系统、数据滤波装置及方法 | |
CN102004854A (zh) | 实时获取配电网非线性负荷等值阻抗参数的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |