CN112881959B - 一种用于磁共振成像的梯度涡流补偿方法及系统 - Google Patents
一种用于磁共振成像的梯度涡流补偿方法及系统 Download PDFInfo
- Publication number
- CN112881959B CN112881959B CN202010524731.1A CN202010524731A CN112881959B CN 112881959 B CN112881959 B CN 112881959B CN 202010524731 A CN202010524731 A CN 202010524731A CN 112881959 B CN112881959 B CN 112881959B
- Authority
- CN
- China
- Prior art keywords
- gradient
- eddy current
- time
- echo
- pulse sequence
- 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
Images
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
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
- G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (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
技术领域
本发明涉及核磁共振成像技术领域,具体涉及一种用于磁共振成像的梯度涡流补偿方法及系统。
背景技术
梯度脉冲,又称为脉冲场梯度(pulsed field gradient),可使磁场随空间发生线性变化,进而使磁场中不同位置的原子核处在不同共振频率下,实现检测样品中不同位置原子核的空间编码,是磁共振成像的关键环节。然而涡流所引起的梯度实际输出波形的改变,破坏了梯度空间编码的准确性,导致图像产生畸变和伪影等失真现象,严重影响磁共振成像质量。因此,消除梯度涡流带来的影响对提升磁共振图像质量至关重要。目前梯度涡流的测量与补偿方法,依赖外部设备,精度较低,并且存在操作复杂,耗时长的缺点。
发明内容
有鉴于此,本发明实施例提供了一种用于磁共振成像的梯度涡流补偿方法及系统,解决现有技术中梯度涡流的测量与补偿方法,依赖外部设备,精度较低,并且存在操作复杂,耗时长的问题。
本发明实施例提供了一种用于磁共振成像的梯度涡流补偿方法,包括:采集预设采样脉冲序列的回波信号,所述预设采样脉冲序列用于测量涡流数据;利用预设算法对所述回波信号进行数据处理,计算得到涡流成分参数;根据所述涡流成分参数生成预加重波形;通过所述预加重波形对磁共振成像梯度脉冲进行涡流补偿,得到预期的目标梯度波形。
可选地,所述采集预设采样脉冲序列的回波信号的步骤,包括:步骤S10:将测试样品放置于预设测试区域,完成预扫描,得到所述测试区域的校正参数;步骤S20:预设置观测时间,利用所述观测时间、测试区域的校正参数及预设采样脉冲序列参数,运行预设采样脉冲序列,得到回波信号,所述预设采样脉冲序列参数包括:测试梯度数据、散相梯度数据及聚相梯度数据,其中所述测试梯度数据包括:测试梯度强度、测试梯度持续时间及测试梯度方向;步骤S30:改变所述测试梯度方向,将预设采样脉冲序列参数替换为改变后的预设采样脉冲序列参数,重复步骤S20,直到完成所有梯度方向的回波信号的采集。
可选地,所述回波信号包括:基础回波信号和参考回波信号,
所述预设置观测时间,利用所述观测时间、测试区域的校正参数及预设采样脉冲序列参数,运行预设采样脉冲序列,得到回波信号的步骤,包括:步骤201:预设置第一观测时间;步骤202:利用第一观测时间、测试区域的校正参数及第一预设采样脉冲序列参数,运行预设采样脉冲序列,得到基础回波信号;步骤203:将所述第一预设采样脉冲序列参数中的测试梯度强度及测试梯度持续时间置为零,得到第二预设采样脉冲序列参数,根据第一观测时间运行预设采样脉冲序列,得到参考回波信号;步骤204:将第二观测时间预设置为第一观测时间,并返回执行步骤S202,直到完成预设采样脉冲序列的所有回波信号的采集。
可选地,所述利用预设算法对所述回波信号进行数据处理,计算得到涡流成分参数的步骤,包括:对所述回波信号进行取模值操作,得到模值回波信号;获取所述模值回波信号的主峰,对所述模值回波信号的主峰进行拟合和插值回归,得到回波信号的聚相时间;通过采样脉冲中心到采样开始时刻的时间加上对应的聚相时间得到回波时间,所述聚相时间包括:对应基础回波信号的聚相时间及对应参考回波信号的参考聚相时间;利用所述回波时间、所述聚相时间和所述聚相梯度数据计算涡流成分参数。
可选地,所述利用所述回波时间、所述聚相时间和所述聚相梯度数据计算涡流成分参数的步骤,包括:将回波时间内所述聚相时间与所述参考聚相时间之差确定为延迟时间;利用所述聚相梯度数据中的聚相梯度强度及所述延迟时间,计算得到涡流积分值;利用所述回波时间、涡流积分值、涡流积分函数模型,基于预设约束模型计算得到涡流成分参数,所述涡流成分参数包括:指数函数的幅度及指数函数的衰减时间常数。
可选地,通过以下公式计算涡流积分值:
Ieddy[m]=Grph(δrph[m]-δrph[ref]),
其中,Ieddy[m]表示第m组涡流积分值,Grph表示聚相梯度强度,δrph[m]表示第m组聚相时间,δrph[ref]表示参考聚相时间。
可选地,所述利用所述回波时间、涡流积分值及涡流积分函数模型,计算得到涡流成分参数的步骤,包括:
利用所述回波时间及涡流积分值,通过以下公式计算涡流函数模型:
其中,Geddy(t)表示涡流函数模型,n表示指数成分的序号,n=1,2,…,N,A[n]表示指数函数的幅度,T[n]表示指数函数的衰减时间常数;
利用所述观测时间、所述回波时间、所述涡流函数模型及所述涡流积分值,通过以下公式计算涡流成分参数:
其中,Ieddy[m]表示第m组涡流积分值,τobs[m]表示第m组观测时间,τecho[m]表示第m组回波时间,Geddy(t)表示涡流函数模型。
本发明实施例还提供了一种用于磁共振成像的梯度涡流补偿系统,包括:采集模块,用于采集预设采样脉冲序列的回波信号;计算模块,用于利用预设算法对所述回波信号进行数据处理,计算得到涡流成分参数;波形生成模块,用于根据所述涡流成分参数生成预加重波形;校正模块,用于通过所述预加重波形对磁共振成像梯度脉冲进行涡流补偿,得到预期的目标梯度波形。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储计算机指令,从而执行本发明实施例提供的用于磁共振成像的梯度涡流补偿方法。
本发明实施例还提供了一种电子设备,包括:存储器和处理器,所述存储器和所述处理器之间互相通信连接,所述存储器中存储有计算机指令,所述处理器通过执行所述计算机指令,用于所述计算机指令执行本发明实施例提供的用于磁共振成像的梯度涡流补偿方法。
本发明技术方案,具有如下优点:
1.本发明提供的用于磁共振成像的梯度涡流补偿方法及系统,通过采集预设采样脉冲序列的回波信号,利用预设算法对回波信号进行数据处理,计算得到涡流成分参数,对磁共振成像梯度脉冲进行梯度涡流补偿,得到预期的目标梯度波形;测试与补偿的过程中无需额外硬件设备支持,噪声对确定回波位置产生的影响小,通过回波位置分析数据的方法抗噪能力较强,且不易受主磁场及射频场不稳定性等问题的影响,稳健性高;并且操作便捷,信号采集速度快,数据处理结果可靠,能够高效、准确地实现梯度涡流的测量与补偿。
2.本发明提供的用于磁共振成像的梯度涡流补偿方法及系统,对测试样品的形状、材质没有严格要求,只需单组分、弛豫时间较长的液体样品;无需精确测量样品位置,样品一经放置,直至完成所有梯度方向的涡流测量前,无需再次移动其位置;整个样品所有质子均参与信号贡献,信号量大,信噪比高,有效解决了传统方法数据采集过程操作繁琐、耗时长等问题,克服了传统方法数据分析结果易受噪声等不确定性因素影响的问题。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中提供的传统用于磁共振成像的梯度涡流补偿方法的流程图;
图2为本发明实施例提供的信号相位法的示意图;
图3为本发明实施例提供的梯度回波相位法的示意图;
图4为本发明实施例提供的采样脉冲序列的示意图;
图5为本发明实施例提供的运行预设采样脉冲序列的一个具体示例的流程图;
图6为本发明实施例提供的未受涡流影响的参考回波信号与受涡流影响的基础回波信号原理示意图;
图7为本发明实施例提供的模值回波信号的示意图;
图8为本发明实施例提供的主峰拟合和插值回归的示意图;
图9为本发明实施例提供的涡流积分的拟合效果图;
图10为本发明实施例提供的用于磁共振成像的梯度涡流补偿系统的模块组成示意图;
图11为本发明实施例提供的计算机设备一个具体示例的组成图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
此外,下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
核磁共振能够做到无损(非侵入)检测、定量检测,以及多模态检测,基于这些优势,核磁共振技术及其应用受到广泛关注,并将在今后很长一段时间内持续发展。并且核磁共振的信号采集通过脉冲序列来实现,这里的脉冲主要指射频脉冲和梯度脉冲。其中梯度脉冲,又称为脉冲场梯度(pulsed field gradient),可使磁场随空间发生线性变化,进而使磁场中不同位置的原子核处在不同共振频率下,是磁共振成像的关键环节,梯度脉冲的输入波形一般是矩形或梯形,然而实际输出却无法满足。因为梯度磁场在经历瞬态变化过程时,梯度线圈周围的金属部件内部会产生阻碍磁场变化的感应电流,并形成涡旋形的闭合回路,即涡电流,简称涡流(eddy current)。由涡流所引起的梯度实际输出波形的改变,破坏了梯度空间编码的准确性,导致图像产生畸变和伪影等失真现象,涡流随时间的变化表现为多指数衰减的形式,通过设置磁共振设备的预加重参数,可产生相应的过电流成分,灵活地实现涡流补偿。
如图1所示,该用于磁共振成像的梯度涡流补偿方法具体包括:
步骤S1:采集预设采样脉冲序列的回波信号,其中预设采样脉冲序列用于测量涡流数据。
本发明实施例中,采集预设采样脉冲序列的回波信号,其中预设采样脉冲序列用于测量涡流数据,设置预设采样脉冲序列的各项参数,包括设置观测时间τobs后采集一组回波信号,改变τobs再次采集,直到完成所有信号采集。例如共采集了8组不同τobs条件下的回波信号,其中τobs从0.1到1000毫秒对数等间隔取值,当τobs取1000毫秒时,可认为测试梯度的涡流已完全衰减,故该组回波信号可同时作为参考回波信号使用。
目前梯度涡流的测量方法,如图2所示,信号相位法通过关闭梯度开关,等待一定时间(τ)采集磁共振信号,则能够计算出信号的相位变化率,上述信号相位法存在操作繁琐(测试过程需要改变样品位置并对样品精确定位)、信号量小、准确度低的问题;如图3所示,梯度回波相位法通过关闭梯度开关,等待一定时间(τ)后利用射频脉冲激发信号,再施加一对极性相反的梯度脉冲,采集回波信号,对采集到的回波信号进行傅里叶变换,则涡流的大小可以计算得到,但是梯度回波相位法存在理论漏洞:(1)由于涡流的存在,回波峰比预期提早出现,导致真实回波时间比序列设置回波时间小,且涡流越大回波峰出现越早,真实回波时间越小,故用该方法测得的涡流小于真值,涡流越大该方法的偏差越大;(2)该方法测得的涡流大小为整个回波时间里涡流的几何平均值,并非观测时刻处的涡流大小,真实涡流水平被低估。本发明实施例提供的用于磁共振成像的梯度涡流补偿方法无需额外硬件设备,如线圈、电路、数字信号接收机等的支持;抗噪能力较强,且不易受主磁场及射频场不稳定性等问题的影响。
需要说明的是,本发明实施中的预设采样脉冲序列是根据实际核磁共振成像系统进行设定的,在实际应用中,也可以根据实际图像需要对预设采样脉冲序列进行修改,并且本发明实施例仅举例说明采集的8组信号以及参考回波信号的确定方法,本发明并不以此为限。
步骤S2:利用预设算法对回波信号进行数据处理,计算得到涡流成分参数。
本发明实施例中,信号采集完毕后,对采集到的回波信号进行数据处理,其中回波信号包括:基础回波信号和参考回波信号,利用预设算法对回波信号进行数据处理,计算得到涡流成分参数的过程为:首先对采集到的回波信号进行取模值操作,然后对取模值操作后得到的模值回波信号的主峰进行拟合和插值回归,得到高分辨的聚相时间;之后通过聚相时间计算出回波时间τecho并计算回波时间内测试梯度涡流的积分值Ieddy;以上述8组回波信号为例,参考信号的聚相时间等于第8组信号的聚相时间,涡流积分函数模型中,除待求解的指数函数的幅度A和衰减时间常数T外,其余变量全部确定。
需要说明的是,本发明实施例仅仅举例说明,采集的回波信号为8组数据,并且计算聚相时间、回波时间以及求解涡流积分函数的方法都可以选择多种方式进行求解,本发明并不以此为限。
步骤S3:根据涡流成分参数生成预加重波形。
本发明实施例中,计算出涡流成分参数后,可准确拟合实际梯度涡流情况,将成分参数输入预加重模块,生成预加重波形实现涡流补偿。
步骤S4:通过预加重波形对磁共振成像梯度脉冲进行涡流补偿,得到预期的目标梯度波形。
本发明实施例中,通过预加重波形对磁共振成像梯度脉冲进行涡流补偿,得到预期的目标梯度波形,其中,预加重(pre-emphasis)是从梯度脉冲输入端解决涡流的技术,是在输入梯度电流时,通过梯度线圈的预加重单元额外输入一个“过电流”,形成失真的梯度波形,失真部分若与涡流幅度相同、方向相反,则两者可相互抵消,最后得到预期的梯度波形。即设置磁共振设备的预加重参数,可产生相应的过电流成分,灵活地实现涡流补偿,最终得到预期的目标梯度波形,避免了由涡流所引起的梯度实际输出波形的改变,破坏了梯度空间编码的准确性,导致图像产生畸变和伪影等失真现象,严重影响磁共振成像质量的问题。
本发明提供的用于磁共振成像的梯度涡流补偿方法,通过采集预设采样脉冲序列的回波信号,利用预设算法对回波信号进行数据处理,计算得到涡流成分参数,对磁共振成像梯度脉冲进行梯度涡流补偿,得到预期的目标梯度波形;测试与补偿的过程中无需额外硬件设备支持,噪声对确定回波位置产生的影响小,通过回波位置分析数据的方法抗噪能力较强,且不易受主磁场及射频场不稳定性等问题的影响,稳健性高;并且操作便捷,信号采集速度快,数据处理结果可靠,能够高效、准确地实现梯度涡流的测量与补偿。
在一具体实施例中,执行步骤S1的过程可以具体包括如下步骤:
步骤S10:将测试样品放置于预设测试区域,完成预扫描,得到测试区域的校正参数。
本发明实施例中,优选纯水作为测试样品,放置样品后进行包括寻找中心频率、匀场、寻找射频脉冲功率与90°脉宽在内的常规预扫描操作,完成预扫描,得到测试区域的校正参数,为后续测试更好的排除噪声干扰,进而达到优化预设采样脉冲序列的目的;其中,优选纯水作为测试样品是因为本发明实施例中对测试样品的形状、材质没有严格要求,只需单组分、弛豫时间较长的液体样品,因此选用生活中常见的纯水作为测试样品,并且整个样品所有质子均参与信号贡献,信号量大,信噪比高,无需精确测量样品位置,即预设测试区域,样品一经放置,直至完成所有梯度方向的涡流测量前,无需再次移动其位置。
需要说明的是,本发明实施例中,仅仅举例说明完成预扫描后得到的测试区域的校正参数,实际应用还可以包括其他参数,可根据实际精度需求进行计算,本发明并不以此为限。
步骤S20:预设置观测时间,利用观测时间、测试区域的校正参数及预设采样脉冲序列参数,运行预设采样脉冲序列,得到回波信号,预设采样脉冲序列参数包括:测试梯度数据、散相梯度数据及聚相梯度数据,其中测试梯度数据包括:测试梯度强度、测试梯度持续时间及测试梯度方向。
本发明实施例中,利用如图4所示的采样脉冲序列采集磁共振信号,采集前,首先设置预设采样脉冲序列参数,例如设置测试梯度方向,如x、y或z;设置采样间隔或谱宽(两者互为倒数);设置用于信号采集的射频脉冲宽度;设置测试梯度强度Gtest,其持续时间δtest可设置为一较大值,例如20毫秒;设置散相梯度强度Gdph及散相时间δdph;设置聚相梯度强度Grph;并且信号采集从散相梯度脉冲关闭,即聚相梯度脉冲开启处开始,一直持续到保证能够观测到完整的回波信号,聚相梯度脉冲在信号采集期间始终开启,直到信号采集结束时关闭,观察采集到的磁共振信号的模值,其最大值所处位置即回波出现时刻,从信号采集开始时刻到回波出现时刻的这段时间即聚相时间δrph。预设置观测时间τobs,即从测试梯度脉冲关闭到射频脉冲中心的时间,运行脉冲序列采集磁共振信号,得到回波信号。
需要说明的是,本发明实施例仅仅举例说明设置的预设采样脉冲序列参数,其中预设采样脉冲序列参数包括:测试梯度数据、散相梯度数据及聚相梯度数据,并且测试梯度数据包括:测试梯度强度、测试梯度持续时间及测试梯度方向;散相梯度数据包括:散相梯度强度、散相梯度持续时间及散相梯度方向;聚相梯度数据包括:聚相梯度强度、聚相梯度持续时间及聚相梯度方向;还可以根据实际需要对预设采样脉冲序列参数进行调整与增加,本发明并不以此为限。
步骤S30:改变测试梯度方向,将预设采样脉冲序列参数替换为改变后的预设采样脉冲序列参数,重复步骤S20,直到完成所有梯度方向的回波信号的采集。
在一具体实施例中,如图5所示,执行步骤S20的过程可以具体包括如下步骤:
步骤S201:预设置第一观测时间。
步骤S202:利用第一观测时间、测试区域的校正参数及第一预设采样脉冲序列参数,运行预设采样脉冲序列,得到基础回波信号。
本发明实施例中,利用第一观测时间、测试区域的校正参数及第一预设采样脉冲序列参数,运行预设采样脉冲序列,得到基础回波信号。需要说明的是,本发明实施例中的基础回波信号就是在第一观测时间内,设置测试区域的校正参数及第一预设采样脉冲序列参数,然后运行预设采样脉冲序列,得到的回波信号作为基础回波信号。
步骤S203:将第一预设采样脉冲序列参数中的测试梯度强度及测试梯度持续时间置为零,得到第二预设采样脉冲序列参数,根据第一观测时间运行预设采样脉冲序列,得到参考回波信号。
本发明实施例中,采集一组未受测试梯度涡流影响的信号作为参考回波信号,其聚相时间为δrph[ref],如图6所示,其中可以通过将第一预设采样脉冲序列参数中的测试梯度强度及测试梯度持续时间置为零,得到第二预设采样脉冲序列参数,根据第一观测时间运行预设采样脉冲序列,得到参考回波信号,需要说明的是,本发明实施例中,还可以通过不改变测试梯度Gtest的强度和持续时间,设置观测时间τobs为一较大值,如1000毫秒,然后采集信号,可认为此时涡流已完全衰减,在实际应用中,还可以根据实际系统功能要求,利用其他方法获取参考回波信号即未受涡流影响的回波信号,本发明并不以此为限。
步骤S204:将第二观测时间预设置为第一观测时间,并返回执行步骤S202,直到完成预设采样脉冲序列的所有回波信号的采集。
本发明实施例中,改变第一观测时间,用第二观测时间代替第一观测时间,即改变τobs的大小并重复上述步骤S202操作,总共进行M次信号采集,得到M组回波信号,以m代表采集的序号,则δrph[m]代表第m次信号采集的聚相时间,如图6所示,通过本发明实施例的方法采集到的磁共振信号,其回波出现位置受测试梯度脉冲产生的涡流影响发生了偏移。
在一具体实施例中,执行步骤S2的过程可以具体包括如下步骤:
步骤S21:对回波信号进行取模值操作,得到模值回波信号。
本发明实施例中,将采集到的预设采样脉冲序列的回波信号进行取模值操作,得到模值回波信号,其中主要是要将采集到的磁共振信号的回波信号中的复数信号进行取模。
步骤S22:获取模值回波信号的主峰,对模值回波信号的主峰进行拟合和插值回归,得到回波信号的聚相时间。
本发明实施例中,获取模值回波信号的主峰,对模值回波信号的主峰进行拟合和插值回归,得到回波信号的聚相时间。以任意梯度方向采集的第m组回波信号为例,对复数信号数据取模值,并通过回波出现时刻,即模值回波信号最大值所处的位置来确定聚相时间δrph[m],由于存在采样时间间隔、数据噪声等影响因素,对主回波峰进行拟合并插值回归可以得到更精确的聚相时间,其中采用多项式拟合主回波峰,当多项式的次数不小于3时,得到的结果比较准确、可靠,拟合得到多项式的各项系数后,可基于该多项式函数模型进行插值回归,得到高精度的聚相时间。
需要说明的是,本发明实施例中仅举例说明可以通过对主回波峰进行拟合并插值回归得到更精确的聚相时间,在实际应用中还可以根据实际经验及系统需求进行设定,并且对主回波峰进行拟合所采用的函数模型并不唯一,还可以采用如洛伦兹线型函数、高斯线型函数、广延指数分布函数等进行拟合,本发明并不以此为限。
步骤S23:通过采样脉冲中心到采样开始时刻的时间加上对应的聚相时间得到回波时间,其中聚相时间包括:对应基础回波信号的聚相时间及对应参考回波信号的参考聚相时间。
步骤S24:利用回波时间、聚相时间和聚相梯度数据计算涡流成分参数。
在一具体实施例中,执行步骤S24的过程可以具体包括如下步骤:
步骤S241:将回波时间内聚相时间与参考聚相时间之差确定为延迟时间。
本发明实施例中,以第m组回波信号为例,该组信号的回波时间τecho[m]为采样射频脉冲中心到采样开始时刻的时间加上其对应的聚相时间δrph[m],在这段回波时间内,所有梯度及其涡流随时间的积分值总和为零。对于参考信号,其聚相时间为δrph[ref],整个回波时间内不存在测试梯度脉冲产生的涡流,由于散相梯度与测试梯度方向相同,而聚相梯度与之方向相反,故对于测试信号,其回波时刻因涡流的叠加而产生延迟,延迟的这段时间内,聚相梯度脉冲起到了弥补涡流的作用。
步骤S242:利用聚相梯度数据中的聚相梯度强度及延迟时间,计算得到涡流积分值。
本发明实施例中,利用聚相梯度数据中的聚相梯度强度及延迟时间,通过以下公式计算涡流积分值:
Ieddy[m]=Grph(δrph[m]-δrph[ref]) (1)
其中,Ieddy[m]表示第m组涡流积分值,Grph表示聚相梯度强度,δrph[m]表示第m组聚相时间,δrph[ref]表示参考聚相时间。
步骤S243:利用回波时间、涡流积分值、涡流积分函数模型,基于预设约束模型计算得到涡流成分参数,涡流成分参数包括:指数函数的幅度及指数函数的衰减时间常数。
本发明实施例中,利用回波时间及涡流积分值,通过以下公式计算涡流函数模型:
其中,Geddy(t)表示涡流函数模型,n表示指数成分的序号,n=1,2,…,N,A[n]表示指数函数的幅度,T[n]表示指数函数的衰减时间常数;
利用观测时间、回波时间、涡流函数模型及涡流积分值,通过以下公式计算涡流成分参数:
其中,Ieddy[m]表示第m组涡流积分值,τobs[m]表示第m组观测时间,τecho[m]表示第m组回波时间,Geddy(t)表示涡流函数模型。
在实际应用中,可以通过最小二乘拟合,因此涡流成分参数的计算可通过求解以下公式实现:
其中,Ieddy[m]表示第m组涡流积分值,τobs[m]表示第m组观测时间,τecho[m]表示第m组回波时间,Geddy(t)表示涡流函数模型,A[n]表示指数函数的幅度,T[n]表示指数函数的衰减时间常数。最小二乘拟合目标为最小化涡流积分观测值与涡流积分函数期望值距离的平方和,最小化过程需要考虑针对涡流成分幅度A和衰减时间常数T的约束条件:一是所有A必须为正值,且其总和不超过测试梯度强度Gtest;二是所有T必须落在[Tmin,Tmax]范围内,该范围应根据磁共振设备梯度预加重模块提供的衰减时间常数的可行范围进行调整。数据拟合过程中,Ieddy[m]的具体数值由公式(1)给出,理论上,指数成分的数目N越大,拟合效果越好。然而在实际情况中,磁共振设备提供的预加重通道总数有限,N值若过大则无法实现所有指数成分的涡流补偿,故N的取值不应超过预加重通道总数。
本发明实施例采用经典的列文伯格-马夸尔特非线性最小二乘算法结合对数障碍函数法可高效、稳健地实现公式(4)的约束优化,求解出准确、可靠的涡流成分参数。需要说明的是,考虑约束条件的求解方法并不唯一,还可以是其它算法,如梯度下降法、牛顿法、信赖域法,无论采用何种算法求解公式(4)或其等效变型式,都可以进行计算,只要能达到计算求解的目的即可,本发明并不以此为限。
实际应用中,首先对采集到的复数信号进行取模值操作,得到的模值回波信号如图7所示,然后对每组回波信号主峰进行拟合和插值回归,得到高分辨的聚相时间δrph,图8是第一组回波信号,即τobs取0.1毫秒时的主峰拟合和插值回归示意图,之后通过聚相时间计算出回波时间τecho并计算回波时间内测试梯度涡流的积分值Ieddy。参考信号的聚相时间等于示例中的第八组信号的聚相时间(δrph[ref]=δrph[8])。此时,涡流积分函数模型中,除待求解的指数成分幅度A和衰减时间常数T外,其余变量全部确定。涡流积分的拟合效果如图9所示,可见本发明提供的方法可准确拟合实际梯度涡流情况。
本发明提供的用于磁共振成像的梯度涡流补偿方法,通过采集预设采样脉冲序列的回波信号,利用预设算法对回波信号进行数据处理,计算得到涡流成分参数,对磁共振成像梯度脉冲进行梯度涡流补偿,得到预期的目标梯度波形;测试与补偿的过程中无需额外硬件设备支持,噪声对确定回波位置产生的影响小,通过回波位置分析数据的方法抗噪能力较强,且不易受主磁场及射频场不稳定性等问题的影响,稳健性高;并且操作便捷,信号采集速度快,数据处理结果可靠,能够高效、准确地实现梯度涡流的测量与补偿;对测试样品的形状、材质没有严格要求,只需单组分、弛豫时间较长的液体样品;无需精确测量样品位置,样品一经放置,直至完成所有梯度方向的涡流测量前,无需再次移动其位置;整个样品所有质子均参与信号贡献,信号量大,信噪比高,有效解决了传统方法数据采集过程操作繁琐、耗时长等问题,克服了传统方法数据分析结果易受噪声等不确定性因素影响的问题。
本发明实施例还提供了一种用于磁共振成像的梯度涡流补偿系统,如图10所示,包括:
采集模块1,用于采集预设采样脉冲序列的回波信号。详细内容参见上述方法实施例中步骤S1的相关描述,在此不再赘述。
计算模块2,用于利用预设算法对回波信号进行数据处理,计算得到涡流成分参数。详细内容参见上述方法实施例中步骤S2的相关描述,在此不再赘述。
波形生成模块3,用于根据涡流成分参数生成预加重波形。详细内容参见上述方法实施例中步骤S3的相关描述,在此不再赘述。
校正模块4,用于通过预加重波形对磁共振成像梯度脉冲进行涡流补偿,得到预期的目标梯度波形。详细内容参见上述方法实施例中步骤S4的相关描述,在此不再赘述。
通过上述各个模块组成部分的协同合作,本发明提供的用于磁共振成像的梯度涡流补偿系统,通过采集预设采样脉冲序列的回波信号,利用预设算法对回波信号进行数据处理,计算得到涡流成分参数,对磁共振成像梯度脉冲进行梯度涡流补偿,得到预期的目标梯度波形;测试与补偿的过程中无需额外硬件设备支持,噪声对确定回波位置产生的影响小,通过回波位置分析数据的方法抗噪能力较强,且不易受主磁场及射频场不稳定性等问题的影响,稳健性高;并且操作便捷,信号采集速度快,数据处理结果可靠,能够高效、准确地实现梯度涡流的测量与补偿;对测试样品的形状、材质没有严格要求,只需单组分、弛豫时间较长的液体样品;无需精确测量样品位置,样品一经放置直至完成所有梯度方向的涡流测量前,无需再次移动其位置;整个样品所有质子均参与信号贡献,信号量大,信噪比高,有效解决了传统方法数据采集过程操作繁琐、耗时长等问题,克服了数据分析结果易受噪声等不确定性因素影响的问题。
本发明实施例提供一种计算机设备,如图11所示,包括:至少一个处理器401,例如CPU(Central Processing Unit,中央处理器),至少一个通信接口403,存储器404,至少一个通信总线402。其中,通信总线402用于实现这些组件之间的连接通信。其中,通信接口403可以包括显示屏(Display)、键盘(Keyboard),可选通信接口403还可以包括标准的有线接口、无线接口。存储器404可以是高速RAM存储器(Ramdom Access Memory,易挥发性随机存取存储器),也可以是非不稳定的存储器(non-volatile memory),例如至少一个磁盘存储器。存储器404可选的还可以是至少一个位于远离前述处理器401的存储装置。其中处理器401可以执行用于磁共振成像的梯度涡流补偿方法。存储器404中存储一组程序代码,且处理器401调用存储器404中存储的程序代码,以用于执行上述的用于磁共振成像的梯度涡流补偿方法。
其中,通信总线402可以是外设部件互连标准(peripheral componentinterconnect,简称PCI)总线或扩展工业标准结构(extended industry standardarchitecture,简称EISA)总线等。通信总线402可以分为地址总线、数据总线、控制总线等。为便于表示,图11中仅用一条线表示,但并不表示仅有一根总线或一种类型的总线。
其中,存储器404可以包括易失性存储器(英文:volatile memory),例如随机存取存储器(英文:random-access memory,缩写:RAM);存储器也可以包括非易失性存储器(英文:non-volatile memory),例如快闪存储器(英文:flash memory),硬盘(英文:hard diskdrive,缩写:HDD)或固降硬盘(英文:solid-state drive,缩写:SSD);存储器404还可以包括上述种类的存储器的组合。
其中,处理器401可以是中央处理器(英文:central processing unit,缩写:CPU),网络处理器(英文:network processor,缩写:NP)或者CPU和NP的组合。
其中,处理器401还可以进一步包括硬件芯片。上述硬件芯片可以是专用集成电路(英文:application-specific integrated circuit,缩写:ASIC),可编程逻辑器件(英文:programmable logic device,缩写:PLD)或其组合。上述PLD可以是复杂可编程逻辑器件(英文:complex programmable logic device,缩写:CPLD),现场可编程逻辑门阵列(英文:field-programmable gate array,缩写:FPGA),通用阵列逻辑(英文:generic arraylogic,缩写:GAL)或其任意组合。
可选地,存储器404还用于存储程序指令。处理器401可以调用程序指令,实现如本申请的用于磁共振成像的梯度涡流补偿方法。
本发明实施例还提供一种计算机可读存储介质,计算机可读存储介质上存储有计算机可执行指令,该计算机可执行指令可执行用于磁共振成像的梯度涡流补偿方法。其中,存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)、随机存储记忆体(Random Access Memory,RAM)、快闪存储器(Flash Memory)、硬盘(Hard Disk Drive,缩写:HDD)或固降硬盘(Solid-State Drive,SSD)等;存储介质还可以包括上述种类的存储器的组合。
显然,上述实施例仅仅是为清楚地说明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明创造的保护范围之中。
Claims (8)
1.一种用于磁共振成像的梯度涡流补偿方法,其特征在于,包括:
采集预设采样脉冲序列的回波信号,所述预设采样脉冲序列用于测量涡流数据,所述回波信号包括:基础回波信号和参考回波信号;
利用预设算法对所述回波信号进行数据处理,计算得到涡流成分参数;
根据所述涡流成分参数生成预加重波形;
通过所述预加重波形对磁共振成像梯度脉冲进行涡流补偿,得到预期的目标梯度波形;
其中,所述采集预设采样脉冲序列的回波信号的步骤,包括:
步骤S10:将测试样品放置于预设测试区域,完成预扫描,得到所述测试区域的校正参数;
步骤S20:预设置观测时间,利用所述观测时间、测试区域的校正参数及预设采样脉冲序列参数,运行预设采样脉冲序列,得到回波信号,所述预设采样脉冲序列参数包括:测试梯度数据、散相梯度数据及聚相梯度数据,其中所述测试梯度数据包括:测试梯度强度、测试梯度持续时间及测试梯度方向;
步骤S30:改变所述测试梯度方向,将预设采样脉冲序列参数替换为改变后的预设采样脉冲序列参数,重复步骤S20,直到完成所有梯度方向的回波信号的采集;
所述利用预设算法对所述回波信号进行数据处理,计算得到涡流成分参数的步骤,包括:
对所述回波信号进行取模值操作,得到模值回波信号;
获取所述模值回波信号的主峰,对所述模值回波信号的主峰进行拟合和插值回归,得到回波信号的聚相时间;
通过采样脉冲中心到采样开始时刻的时间加上对应的聚相时间得到回波时间,所述聚相时间包括:对应基础回波信号的聚相时间及对应参考回波信号的参考聚相时间;
利用所述回波时间、所述聚相时间和所述聚相梯度数据计算涡流成分参数。
2.根据权利要求1所述的用于磁共振成像的梯度涡流补偿方法,其特征在于,
所述预设置观测时间,利用所述观测时间、测试区域的校正参数及预设采样脉冲序列参数,运行预设采样脉冲序列,得到回波信号的步骤,包括:
步骤201:预设置第一观测时间;
步骤202:利用第一观测时间、测试区域的校正参数及第一预设采样脉冲序列参数,运行预设采样脉冲序列,得到基础回波信号;
步骤203:将所述第一预设采样脉冲序列参数中的测试梯度强度及测试梯度持续时间置为零,得到第二预设采样脉冲序列参数,根据第一观测时间运行预设采样脉冲序列,得到参考回波信号;
步骤204:将第二观测时间预设置为第一观测时间,并返回执行步骤S202,直到完成预设采样脉冲序列的所有回波信号的采集。
3.根据权利要求1所述的用于磁共振成像的梯度涡流补偿方法,其特征在于,所述利用所述回波时间、所述聚相时间和所述聚相梯度数据计算涡流成分参数的步骤,包括:
将回波时间内所述聚相时间与所述参考聚相时间之差确定为延迟时间;
利用所述聚相梯度数据中的聚相梯度强度及所述延迟时间,计算得到涡流积分值;
利用所述回波时间、涡流积分值、涡流积分函数模型,基于预设约束模型计算得到涡流成分参数,所述涡流成分参数包括:指数函数的幅度及指数函数的衰减时间常数。
4.根据权利要求3所述的用于磁共振成像的梯度涡流补偿方法,其特征在于,通过以下公式计算涡流积分值:
Ieddy[m]=Grph(δrph[m]-δrph[ref]),
其中,Ieddy[m]表示第m组涡流积分值,Grph表示聚相梯度强度,δrph[m]表示第m组聚相时间,δrph[ref]表示参考聚相时间。
5.根据权利要求3所述的用于磁共振成像的梯度涡流补偿方法,其特征在于,所述利用所述回波时间、涡流积分值及涡流积分函数模型,计算得到涡流成分参数的步骤,包括:
利用所述回波时间及涡流积分值,通过以下公式计算涡流函数模型:
其中,Geddy(t)表示涡流函数模型,n表示指数成分的序号,n=1,2,…,N,A[n]表示指数函数的幅度,T[n]表示指数函数的衰减时间常数;
利用所述观测时间、所述回波时间、所述涡流函数模型及所述涡流积分值,通过以下公式计算涡流成分参数:
其中,Ieddy[m]表示第m组涡流积分值,τobs[m]表示第m组观测时间,τecho[m]表示第m组回波时间,Geddy(t)表示涡流函数模型。
6.一种用于磁共振成像的梯度涡流补偿系统,其特征在于,包括:
采集模块,用于采集预设采样脉冲序列的回波信号,所述预设采样脉冲序列用于测量涡流数据,所述回波信号包括:基础回波信号和参考回波信号,其中采集预设采样脉冲序列的回波信号的步骤,包括:
步骤S10:将测试样品放置于预设测试区域,完成预扫描,得到所述测试区域的校正参数;
步骤S20:预设置观测时间,利用所述观测时间、测试区域的校正参数及预设采样脉冲序列参数,运行预设采样脉冲序列,得到回波信号,所述预设采样脉冲序列参数包括:测试梯度数据、散相梯度数据及聚相梯度数据,其中所述测试梯度数据包括:测试梯度强度、测试梯度持续时间及测试梯度方向;
步骤S30:改变所述测试梯度方向,将预设采样脉冲序列参数替换为改变后的预设采样脉冲序列参数,重复步骤S20,直到完成所有梯度方向的回波信号的采集;
计算模块,用于利用预设算法对所述回波信号进行数据处理,计算得到涡流成分参数,包括:
对所述回波信号进行取模值操作,得到模值回波信号;
获取所述模值回波信号的主峰,对所述模值回波信号的主峰进行拟合和插值回归,得到回波信号的聚相时间;
通过采样脉冲中心到采样开始时刻的时间加上对应的聚相时间得到回波时间,所述聚相时间包括:对应基础回波信号的聚相时间及对应参考回波信号的参考聚相时间;
利用所述回波时间、所述聚相时间和所述聚相梯度数据计算涡流成分参数;
波形生成模块,用于根据所述涡流成分参数生成预加重波形;
校正模块,用于通过所述预加重波形对磁共振成像梯度脉冲进行涡流补偿,得到预期的目标梯度波形。
7.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储计算机指令,所述计算机指令被处理器执行时实现如权利要求1-5中任一项所述的用于磁共振成像的梯度涡流补偿方法。
8.一种电子设备,其特征在于,包括:
存储器和处理器,所述存储器和所述处理器之间互相通信连接,所述存储器中存储有计算机指令,所述处理器通过执行所述计算机指令,从而执行如权利要求1-5中任一项所述的用于磁共振成像的梯度涡流补偿方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010524731.1A CN112881959B (zh) | 2020-06-10 | 2020-06-10 | 一种用于磁共振成像的梯度涡流补偿方法及系统 |
PCT/CN2020/119282 WO2021248752A1 (zh) | 2020-06-10 | 2020-09-30 | 一种用于磁共振成像的梯度涡流补偿方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010524731.1A CN112881959B (zh) | 2020-06-10 | 2020-06-10 | 一种用于磁共振成像的梯度涡流补偿方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112881959A CN112881959A (zh) | 2021-06-01 |
CN112881959B true CN112881959B (zh) | 2021-12-21 |
Family
ID=76043031
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010524731.1A Active CN112881959B (zh) | 2020-06-10 | 2020-06-10 | 一种用于磁共振成像的梯度涡流补偿方法及系统 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN112881959B (zh) |
WO (1) | WO2021248752A1 (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114252828A (zh) * | 2021-12-21 | 2022-03-29 | 安徽福晴医疗科技有限公司 | 一种磁共振成像系统及其预加重涡流校正方法和装置 |
CN114325531B (zh) * | 2021-12-31 | 2024-06-25 | 深圳市联影高端医疗装备创新研究院 | 磁共振系统延迟的校正方法、装置、计算机设备和介质 |
CN114509704A (zh) * | 2022-02-15 | 2022-05-17 | 湖南小快智造电子科技有限公司 | 一种安全用电智能监测仪 |
CN115656900B (zh) * | 2022-11-02 | 2023-06-23 | 佛山瑞加图医疗科技有限公司 | 一种用于降低系统误差对磁共振成像影响的方法及装置 |
CN115951281B (zh) * | 2022-12-08 | 2023-10-20 | 无锡鸣石峻致医疗科技有限公司 | 磁共振梯度延迟和梯度切换率同时检测方法、装置和介质 |
CN117233677B (zh) * | 2023-09-15 | 2024-06-04 | 安徽峻德医疗科技有限公司 | 非笛卡尔激发采样的磁共振梯度脉冲信号补偿方法及装置 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6127826A (en) * | 1999-01-27 | 2000-10-03 | Picker International, Inc. | EPI image based long term eddy current pre-emphasis calibration |
CN102103196A (zh) * | 2009-12-18 | 2011-06-22 | 东软飞利浦医疗设备系统有限责任公司 | 具有可调数字预加重的磁共振梯度波形产生系统及方法 |
CN103202694A (zh) * | 2012-01-13 | 2013-07-17 | 株式会社东芝 | 磁共振成像装置 |
CN104181479A (zh) * | 2013-05-23 | 2014-12-03 | 上海联影医疗科技有限公司 | 用于磁共振成像系统的涡流补偿方法 |
CN104614694A (zh) * | 2015-01-27 | 2015-05-13 | 华东师范大学 | 一种磁共振梯度涡流补偿方法 |
CN107677976A (zh) * | 2017-09-26 | 2018-02-09 | 中国科学院武汉物理与数学研究所 | 一种自适应核磁共振梯度预加重波形产生装置及方法 |
CN107907846A (zh) * | 2017-11-27 | 2018-04-13 | 深圳先进技术研究院 | 涡流校正方法、装置、移动终端及可读存储介质 |
CN108279393A (zh) * | 2017-01-05 | 2018-07-13 | 上海康达卡勒幅医疗科技有限公司 | 一种磁共振成像全自动预加重校正涡流的方法 |
CN108303665A (zh) * | 2018-02-27 | 2018-07-20 | 奥泰医疗系统有限责任公司 | 一种磁共振成像系统涡流补偿方法 |
CN109765512A (zh) * | 2019-01-18 | 2019-05-17 | 上海联影医疗科技有限公司 | 磁共振梯度系统及其涡流补偿方法及装置 |
CN109799472A (zh) * | 2019-03-08 | 2019-05-24 | 沈阳工业大学 | 一种基于深度学习的磁共振涡流补偿方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7112964B2 (en) * | 2004-08-02 | 2006-09-26 | General Electric Company | Eddy current measurement and correction in magnetic resonance imaging systems with a static phantom |
CN100339050C (zh) * | 2005-09-20 | 2007-09-26 | 华东师范大学 | 一种磁共振成像系统中减小梯度线圈涡流的方法 |
CN102096054B (zh) * | 2010-12-22 | 2013-03-13 | 华东师范大学 | 一种快速自旋回波脉冲序列的梯度系统延时校正方法 |
CN102298129B (zh) * | 2011-05-30 | 2014-07-09 | 苏州朗润医疗系统有限公司 | 一种用于核磁共振成像系统的涡流测量及补偿方法 |
CN105929350B (zh) * | 2016-05-05 | 2019-01-25 | 南京拓谱医疗科技有限公司 | 一种单次激发水脂分离成像误差校正系统及方法 |
US10324156B2 (en) * | 2017-07-19 | 2019-06-18 | Synaptive Medical (Barbados) Inc. | System and method to correct eddy current artifacts in magnetic resonance imaging |
-
2020
- 2020-06-10 CN CN202010524731.1A patent/CN112881959B/zh active Active
- 2020-09-30 WO PCT/CN2020/119282 patent/WO2021248752A1/zh active Application Filing
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6127826A (en) * | 1999-01-27 | 2000-10-03 | Picker International, Inc. | EPI image based long term eddy current pre-emphasis calibration |
CN102103196A (zh) * | 2009-12-18 | 2011-06-22 | 东软飞利浦医疗设备系统有限责任公司 | 具有可调数字预加重的磁共振梯度波形产生系统及方法 |
CN103202694A (zh) * | 2012-01-13 | 2013-07-17 | 株式会社东芝 | 磁共振成像装置 |
CN104181479A (zh) * | 2013-05-23 | 2014-12-03 | 上海联影医疗科技有限公司 | 用于磁共振成像系统的涡流补偿方法 |
CN104614694A (zh) * | 2015-01-27 | 2015-05-13 | 华东师范大学 | 一种磁共振梯度涡流补偿方法 |
CN108279393A (zh) * | 2017-01-05 | 2018-07-13 | 上海康达卡勒幅医疗科技有限公司 | 一种磁共振成像全自动预加重校正涡流的方法 |
CN107677976A (zh) * | 2017-09-26 | 2018-02-09 | 中国科学院武汉物理与数学研究所 | 一种自适应核磁共振梯度预加重波形产生装置及方法 |
CN107907846A (zh) * | 2017-11-27 | 2018-04-13 | 深圳先进技术研究院 | 涡流校正方法、装置、移动终端及可读存储介质 |
CN108303665A (zh) * | 2018-02-27 | 2018-07-20 | 奥泰医疗系统有限责任公司 | 一种磁共振成像系统涡流补偿方法 |
CN109765512A (zh) * | 2019-01-18 | 2019-05-17 | 上海联影医疗科技有限公司 | 磁共振梯度系统及其涡流补偿方法及装置 |
CN109799472A (zh) * | 2019-03-08 | 2019-05-24 | 沈阳工业大学 | 一种基于深度学习的磁共振涡流补偿方法 |
Non-Patent Citations (3)
Title |
---|
《Quantification and compensation of eddy-current-induced magnetic-field gradients》;William M. Spees等;《Journal of Magnetic Resonance》;20110930;全文 * |
《一种基于分段迭代的MRI成像涡流补偿方法》;郭亚平等;《生物医学工程研究》;20110630;全文 * |
《基于迭代TSVD的NMR二维谱反演算法》;周小龙等;《波谱学杂志》;20131231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112881959A (zh) | 2021-06-01 |
WO2021248752A1 (zh) | 2021-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112881959B (zh) | 一种用于磁共振成像的梯度涡流补偿方法及系统 | |
Buonincontri et al. | MR fingerprinting with simultaneous B1 estimation | |
Hargreaves et al. | Variable‐rate selective excitation for rapid MRI sequences | |
Addy et al. | Simple method for MR gradient system characterization and k‐space trajectory estimation | |
Brodsky et al. | Characterizing and correcting gradient errors in non‐cartesian imaging: are gradient errors linear time‐invariant (LTI)? | |
EP2667215B1 (en) | Method of performing an mri reference scan | |
GB2304901A (en) | MRI magnet shimming using 1-D projection FID signals | |
JPH0349453B2 (zh) | ||
CN105829906B (zh) | 具有水/脂肪分离的相位敏感的反转恢复mri | |
CN106646282B (zh) | 一种基于量化时延法提高fid信号测频精度的方法及电路 | |
US10379190B2 (en) | Method of determining a gradient impulse response function during execution of a magnetic resonance imaging or spectroscopy sequence | |
JP4564015B2 (ja) | 磁気共鳴撮影装置及び磁気共鳴撮影方法 | |
Robson et al. | Consequences of T2 relaxation during half‐pulse slice selection for ultrashort TE imaging | |
JP7023954B2 (ja) | プロペラmrイメージング | |
US11275143B2 (en) | Method for measuring eddy current fields in a magnetic resonance device, magnetic resonance device, computer program and electronically readable data carrier | |
Han et al. | Direct measurement of magnetic field gradient waveforms | |
Splitthoff et al. | SENSE shimming (SSH): a fast approach for determining B0 field inhomogeneities using sensitivity coding | |
Jordanova et al. | B1 estimation using adiabatic refocusing: BEAR | |
JP4456009B2 (ja) | Nmr測定方法および装置 | |
Fischer et al. | Monitoring and compensating phase imperfections in cine balanced steady‐state free precession | |
JP2000279390A (ja) | 磁気共鳴イメージング装置 | |
Renou et al. | Radio-frequency pulse calibration using the MISSTEC sequence | |
Andris et al. | Optimized measurement of magnetic field maps using nuclear magnetic resonance (NMR) | |
Elster et al. | Quantitative MRS: comparison of time domain and time domain frequency domain methods using a novel test procedure | |
CN114264990A (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 |