[发明专利]一种磁共振成像全自动预加重校正涡流的方法

权利要求书

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之后，程序判定为最终迭代收敛，涡流校正完毕。

2.如权利要求1所述磁共振成像全自动预加重校正涡流的方法，其特征在于，

其中，步骤B)中利用对所有信号相位做定量计算得到各个延迟时间点τi的涡流值大小Ci的，表示为公式[1]:

C_i＝m_i/2πγTE [1]

其中，mi表示第i次采集的傅里叶变换后回波信号的斜率(单位：弧度/毫米)，γ表示物理常量旋磁比(＝0.267513弧度/毫秒/微特斯拉)，TE表示回波时间(单位：微特斯拉)；

其中，步骤C)中所用的多E指数衰减的函数表示为公式[2]:

其中，n表示E指数成分的个数，τ_i表示梯度延迟时间，C(τi)表示第i个梯度延迟点采集计算得到的涡流大小，Aj和Kj分别表示第j个E指数衰减成分的幅值和时间常数；

其中，对公式[2]进行拟合时所做的正则化控制可以表示为公式[3]:

其中min{}表示最小化，λ表示正则化控制项的常数。