一种磁共振成像全自动预加重校正涡流的方法技术

本发明专利技术涉及一种磁共振成像全自动预加重校正涡流的方法，在对涡流信号定量过程中使用了Levenberg‑Marquardt非线性正则化最小化拟合，准确而且全自动，可以对磁共振成像系统中涡流衰减函数进行全自动地准确地判定，当涡流衰减形式得到正确的定义后，在迭代预加重校正过程中又将E指数的时间常数固定，有效地将非线性问题转化为线性模型定量，然后通过SVD线性最小化方法对涡流幅值进行迭代校正，实现全过程完全自动化，不需要人工介入。

An automatic pre emphasis correction method for eddy current in magnetic resonance imaging

The present invention relates to a method of fully automatic preweighting correction of eddy current in magnetic resonance imaging, using Levenberg Marquardt nonlinear regularization minimization fitting in the quantitative process of eddy current signals, accurate and fully automatic, which can automatically and accurately determine eddy current attenuation function in magnetic resonance imaging system, when eddy current After the correct definition of the attenuation form, the time constant of the E exponent is fixed in the iterative pre accentuation process, and the nonlinear problem is converted to the linear model quantitatively. Then the eddy current amplitude is iteratively corrected by the SVD linear minimization method to realize the complete automation of the whole process without the need of artificial intervention.

【技术实现步骤摘要】
一种磁共振成像全自动预加重校正涡流的方法
本专利技术涉及磁共振成像领域，尤其涉及一种磁共振成像全自动预加重校正涡流的方法。
技术介绍
磁共振成像具有无损伤、软组织对比度高、任意方向断层等特点，目前已广泛用于医学临床诊断，而在磁共振成像过程中，系统会产生涡流，破坏主磁场和梯度场的均匀性，从而会在磁共振图像中产生伪影，尤其是快速自旋回波(FSE)，弥散成像(DWI)等成像序列会对磁共振成像系统的涡流极为敏感。因此准确有效地定量并且消除磁共振成像系统中的涡流是影响磁共振图像质量最为关键的因素之一。磁共振成像系统的涡流主要产生于施加成像梯度后，梯度电流会在梯度线圈和磁体中感应产生随时间呈自然指数形式衰减的感应电流，从而产生随时间衰减性的B0涡流场和梯度涡流场。涡流的物理形式为由单个或多个E指数组合而成的衰减函数，具体由几个E指数成分组成受系统特性、磁体和梯度线圈的设计生产、硬件设备安装等因素决定，不同的磁共振成像系统会有不同的涡流函数形式，即使在同一个磁共振成像系统中，在三个不同的磁体方向(X，Y和Z)的涡流衰减函数形式也会有不同。目前校正涡流最常用且有效的两种方法：分别是利用自屏蔽线圈和预加重补偿，前者属于被动的方法，是通过在梯度线圈设计和制造上实现对涡流场的自屏蔽作用来抑制梯度线圈中的涡流，但是此方法依赖于线圈制造商的线圈设计和制造工艺，而且并不能完全消除涡流的影响，并且往往会带来负极性的涡流；而预加重补偿则是属于主动的方法，是通过对所施加的梯度波形进行预加重修正，施加与涡流幅值和衰减时间常数大小相同，但极性相反的梯度场，与系统涡流相互抵消，从而消除涡流影响。虽然该方法可以有效校正系统涡流，但在实际的磁共振成像系统调试中，涡流的极性、时间参数以及衰减指数的个数受系统的设备，硬件条件状况变化差异很大，传统的方法利用非线性拟合数学工具对涡流衰减信号无法保证能准确定义涡流衰减模型。因此实际的定量涡流调试依赖于工程师的主观经验判断，往往并不准确，无法给出合适的涡流模型，或者通常的做法是将E指数的成分个数设成为一个固定值(3个或4个)，这样做往往会造成过拟合，因此传统的涡流预加重方法无法真正做到自动化地准确调试校正系统涡流。
技术实现思路
有鉴于此，本专利技术的目的是提供一种磁共振成像全自动预加重校正涡流的方法，以解决现有技术中的不足。为了达到上述目的，本专利技术的目的是通过下述技术方案实现的：提供一种磁共振成像全自动预加重校正涡流的方法，包括如下步骤：A)分别在X、Y或Z方向独立施加一固定幅值位G和持续时间为D的梯度，经过延迟时间τi后，发射射频，并采集梯度回波信号；B)利用对所有信号相位做定量计算得到各个延迟时间点τi的涡流值大小Ci，得到涡流随时间衰减的曲线；C)利用Levenberg-Marquardt算法，对定量得到的涡流衰减曲线Ci做多E指数衰减模型的非线性拟合，并且对所有E指数的幅度值做权重正则化控制，防止出现指数项过多的过拟合发生,自动剔除掉多余的指数项后，将所有E指数衰减常数固定Kj；D)对涡流E指数幅度做SVD线性拟合，根据拟合定量得到的对应于各个时间常数Kj的E指数涡流幅值Aj，并更新谱仪的涡流预加重参数，重新采集涡流信号；E)反复迭代进行以上D过程，当涡流值小于某一设定的阈值T之后，程序判定为最终迭代收敛，涡流校正完毕。上述磁共振成像全自动预加重校正涡流的方法，其中，步骤B)中利用对所有信号相位做定量计算得到各个延迟时间点τi的涡流值大小Ci的，表示为公式[1]:Ci＝mi/2πγTE[1]其中，mi表示第i次采集的傅里叶变换后回波信号的斜率(单位：弧度/毫米)，γ表示物理常量旋磁比(＝0.267513弧度/毫秒/微特斯拉)，TE表示回波时间(单位：微特斯拉)；其中，步骤C)中所用的多E指数衰减的函数表示为公式[2]:其中，n表示E指数成分的个数，τi表示梯度延迟时间，C(τi)表示第i个梯度延迟点采集计算得到的涡流大小，Aj和Kj分别表示第j个E指数衰减成分的幅值和时间常数；其中，对公式[2]进行拟合时所做的正则化控制可以表示为公式[3]:其中min{}表示最小化，λ表示正则化控制项的常数。与已有技术相比，本专利技术的有益效果在于：本专利技术在对涡流信号定量过程中使用了Levenberg-Marquardt非线性正则化最小化拟合，准确而且全自动，可以对磁共振成像系统中涡流衰减函数进行全自动地准确地判定，当涡流衰减形式得到正确的定义后，在迭代预加重校正过程中又将E指数的时间常数固定，有效地将非线性问题转化为线性模型定量，然后通过SVD线性最小化方法对涡流幅值进行迭代校正，实现全过程完全自动化，不需要人工介入。附图说明构成本专利技术的一部分的附图用来提供对本专利技术的进一步理解，本专利技术的示意性实施例及其说明用于解释本专利技术，并不构成对本专利技术的不当限定。在附图中：附图1是本专利技术方法中定量得到的X方向涡流随时间衰减的数值定量曲线；附图2是本专利技术方法中预加重校正前后系统X方向涡流曲线的效果对比图，实线所示为无预加重校正的系统X方向涡流定量曲线，虚线所示为有预加重校正的系统X方向涡流定量曲线；附图3是本专利技术方法中涡流预加重校正前后水模FSE图像验证的效果对比图：a为无涡流预加重校正的水模FSE图像，b为有涡流预加重校正的水模FSE图像。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。需要说明的是，在不冲突的情况下，本专利技术中的实施例及实施例中的特征可以相互组合。本专利技术在磁共振成像系统中的某一方向(X、Y或Z)施加一定量梯度，之后通过采集梯度回波信号来测量涡流值，通过测量不同延迟时间点的涡流值，得到涡流衰减曲线，然后利用Levenberg-Marquardt(列文伯格-马夸尔特法)的最小化算法对涡流衰减信号值做多E指数衰减的非线性拟合。并且在最小化过程中，对所有E指数的幅度值做权重正则化控制，防止出现指数项过多的过拟合发生。然后自动剔除掉多余的指数项后，将所有E指数衰减常数固定，将非线性模型转化为线性模型定量，对涡流E指数幅度做SVD线性拟合，根据拟合定量得到的涡流幅值和时间常数更新谱仪的涡流预加重参数，重新采集涡流信号，反复迭代进行以上过程，迭代过程中所有E指数衰减常数固定，只做涡流E指数幅度做SVD线性拟合，当涡流值小于某一设定的阈值之后，程序判定为最终迭代收敛，涡流校正完毕。全过程完全自动化，不需要人工介入。A)分别在X、Y或Z方向独立施加一固定幅值位G和持续时间为D的梯度，经过延迟时间τi后，发射射频，并采集梯度回波信号；B)利用对所有信号相位做定量计算得到各个延迟时间点τi的涡流值大小Ci，得到涡流随时间衰减的曲线；C)利用Levenberg-Marquardt算法，对定量得到的涡流衰减曲线Ci做多E指数衰减模型的非线性拟合，并且对所有E指数的幅度值做权重正则化控制，防止出现指数项过多的过拟合发生,自动剔除掉多余的指数项后，将所有E指数衰减常数固定Kj；D)对涡流E指数幅度做SVD线性拟合，根据本文档来自技高网...

【技术保护点】
1.一种磁共振成像全自动预加重校正涡流的方法，其特征在于，包括如下步骤：A)分别在X、Y或Z方向独立施加一固定幅值位G和持续时间为D的梯度，经过延迟时间τi后，发射射频，并采集梯度回波信号；B)利用对所有信号相位做定量计算得到各个延迟时间点τi的涡流值大小Ci，得到涡流随时间衰减的曲线；C)利用Levenberg‑Marquardt算法，对定量得到的涡流衰减曲线Ci做多E指数衰减模型的非线性拟合，并且对所有E指数的幅度值做权重正则化控制，防止出现指数项过多的过拟合发生,自动剔除掉多余的指数项后，将所有E指数衰减常数固定Kj；D)对涡流E指数幅度做SVD线性拟合，根据拟合定量得到的对应于各个时间常数Kj的E指数涡流幅值Aj，并更新谱仪的涡流预加重参数，重新采集涡流信号；E)反复迭代进行以上D过程，当涡流值小于某一设定的阈值T之后，程序判定为最终迭代收敛，涡流校正完毕。

【技术特征摘要】
1.一种磁共振成像全自动预加重校正涡流的方法，其特征在于，包括如下步骤：A)分别在X、Y或Z方向独立施加一固定幅值位G和持续时间为D的梯度，经过延迟时间τi后，发射射频，并采集梯度回波信号；B)利用对所有信号相位做定量计算得到各个延迟时间点τi的涡流值大小Ci，得到涡流随时间衰减的曲线；C)利用Levenberg-Marquardt算法，对定量得到的涡流衰减曲线Ci做多E指数衰减模型的非线性拟合，并且对所有E指数的幅度值做权重正则化控制，防止出现指数项过多的过拟合发生,自动剔除掉多余的指数项后，将所有E指数衰减常数固定Kj；D)对涡流E指数幅度做SVD线性拟合，根据拟合定量得到的对应于各个时间常数Kj的E指数涡流幅值Aj，并更新谱仪的涡流预加重参数，重新采集涡流信号；E)反复迭代进行以上D过程，当涡流值小于某一设定的阈值T之后，程序判...

【专利技术属性】
技术研发人员：裴孟超，姜小平，张成秀，李建奇，
申请(专利权)人：上海康达卡勒幅医疗科技有限公司，
类型：发明
国别省市：上海,31