CN108303665B - 一种磁共振成像系统涡流补偿方法 - Google Patents

一种磁共振成像系统涡流补偿方法 Download PDF

Info

Publication number
CN108303665B
CN108303665B CN201810162124.8A CN201810162124A CN108303665B CN 108303665 B CN108303665 B CN 108303665B CN 201810162124 A CN201810162124 A CN 201810162124A CN 108303665 B CN108303665 B CN 108303665B
Authority
CN
China
Prior art keywords
magnetic resonance
eddy current
gradient
waveform
axis
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
Application number
CN201810162124.8A
Other languages
English (en)
Other versions
CN108303665A (zh
Inventor
罗海
王文周
朱高杰
周翔
王超
陈梅泞
刘霞
吴子岳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ALLTECH MEDICAL SYSTEMS LLC
Original Assignee
ALLTECH MEDICAL SYSTEMS LLC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by ALLTECH MEDICAL SYSTEMS LLC filed Critical ALLTECH MEDICAL SYSTEMS LLC
Priority to CN201810162124.8A priority Critical patent/CN108303665B/zh
Priority to PCT/CN2018/080112 priority patent/WO2019165656A1/zh
Publication of CN108303665A publication Critical patent/CN108303665A/zh
Application granted granted Critical
Publication of CN108303665B publication Critical patent/CN108303665B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56518Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field

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

本发明公开一种磁共振成像系统涡流补偿方法,属于磁共振成像技术领域,包括以下步骤:步骤a、磁共振谱仪根据扫描序列参数产生梯度磁场的波形并输出到梯度放大器,驱动磁共振扫描仪工作,并由磁共振谱仪采集磁共振数据;步骤b、根据扫描序列参数和涡流补偿参数,进行梯度轨迹校正,即计算出受涡流影响后的梯度磁场随时间变化的波形,进而得到校正后的K空间轨迹;步骤c、将步骤a采集到的磁共振数据按照步骤b得到的校正后的K空间轨迹进行网格化重建或者利用非均匀傅里叶变换重建出磁共振图像。本发明基于软件进行涡流补偿,避免了硬件投入,提高了成像质量。

Description

一种磁共振成像系统涡流补偿方法
技术领域
本发明涉及磁共振成像技术领域,尤其涉及一种磁共振成像系统涡流补偿方法。
背景技术
在磁共振成像过程中,梯度磁场与主磁场及射频磁场共同组成磁共振成像的三大要素。梯度磁场由梯度线圈产生,梯度线圈中的电流随时间快速切换,会在周围导体结构中产生涡流,涡流会产生一个空间和时间上都不断变化的磁场,使成像区域内的梯度场发生畸变,进而影响成像质量,如图1所示,图中实线为理想梯度磁场波形,虚线为受涡流影响的波形(横轴为时间,纵轴为梯度磁场强度)。因此需要对脉冲序列发出的梯度波形进行校正,使得实际产生的梯度场接近于理想的波形,该过程称为“涡流补偿”。
传统的涡流补偿方法通过梯度预加重(Gradient Pre-emphasis),即依据涡流的数学物理模型,预先加大梯度磁场驱动输入,使得实际产生的梯度磁场逼近期望值,如图2所示,图中虚线为输入梯度波形即梯度预加重后的波形,实线为实际产生的梯度磁场波形(横轴为时间,纵轴为梯度磁场强度)。
基于梯度预加重涡流补偿技术的磁共振成像系统如图3所示,在系统安装校正阶段,根据涡流的数学物理模型测出涡流补偿参数并设置到磁共振扫描工作站;在进行磁共振扫描时,磁共振谱仪根据扫描序列参数产生梯度磁场的波形,该波形通过磁共振谱仪中的梯度预加重模块进行梯度预加重,并输出到梯度放大器,驱动磁共振扫描仪工作,采集到的磁共振信号可直接重建出图像。
但是,这种传统的涡流补偿方案,主要存在以下几方面问题。第一、磁共振谱仪中需要增加梯度预加重模块,增加了磁共振谱仪设计复杂度和硬件成本;第二、梯度预加重使得实际驱动电流增加,需要牺牲5%~10%的性能来保证预加重后的梯度电流不超过梯度放大器的工作区间;第三、梯度预加重模块有数字补偿和模拟补偿两种方案,受系统实时性以及硬件规模的限制,梯度预加重的指数项数目有限,一般商业磁共振系统仅支持4~8个指数项用于涡流补偿;第四,磁共振系统中由梯度放大器,滤波器,梯度线圈等组成的梯度系统,其频率响应带宽有限,一般为5~10KHz,对于超短时涡流,例如时间常数小于0.1ms的涡流,其频率响应会大于10KHz,梯度系统不能很好的响应这种超短时涡流的梯度预加重,使得涡流补偿达不到预期效果。
发明内容
本发明旨在提供一种磁共振成像系统涡流补偿方法,用于解决现有的涡流补偿方式硬件成本高、降低系统性能、使用限制大、响应频段窄的问题。
为达到上述目的,本发明采用的技术方案如下:一种磁共振成像系统涡流补偿方法,包括以下步骤:
步骤a、磁共振谱仪根据扫描序列参数产生梯度磁场的波形并输出到梯度放大器,驱动磁共振扫描仪工作,并由磁共振谱仪采集磁共振数据;
步骤b、根据理想梯度波形和涡流补偿参数进行梯度轨迹校正,即计算出受涡流影响后的梯度磁场随时间变化的波形,进而得到校正后的K空间轨迹;
步骤c、将步骤a采集到的磁共振数据按照步骤b得到的校正后的K空间轨迹进行网格化重建或者利用非均匀傅里叶变换重建出磁共振图像。
进一步的,在步骤b中,包括以下步骤:
步骤b1、根据理想梯度波形和涡流的物理数学模型计算出校正后的物理梯度波形,即X、Y、Z轴受涡流影响后的实际梯度波形;
步骤b2、将校正后的物理梯度波形从物理轴映射到逻辑轴,得到校正后的逻辑梯度波形;
步骤b3、根据校正后的逻辑梯度波形,得出校正后的K空间轨迹。
进一步的,在步骤b1中,X轴受涡流影响后的实际波形计算公式如下:
其中GX(t)、GY(t)、GZ(t)分别为X轴、Y轴、Z轴的理想波形,G′X(t)为X轴受涡流影响的实际波形,为第i个时间常数对应的自相关涡流补偿卷积核函数,为第i个交叉项涡流补偿卷积核函数,即Y轴和Z轴对X轴的影响,卷积核函数表示为如下形式:
其中,分别为系统校正阶段测量到的第i项涡流补偿参数,即根据涡流物理数学模型得到的补偿增益和时间常数,
上述公式计算出了X轴受涡流影响后的实际梯度波形,按照相同原理,分别计算出Y轴和Z轴受涡流影响后的实际梯度波形G′Y(t)和G′Z(t),即得到了校正后的物理轴梯度波形。
进一步的,在步骤b2中,映射关系如下:
其中,R为旋转矩阵,即由切片层面所在坐标系转换到物理坐标系的变换矩阵,R-1为旋转矩阵的逆矩阵,GR(t)、GP(t)、GS(t)分别对应逻辑坐标系中的频率编码方向校正后的逻辑梯度波形、相位编码方向校正后的逻辑梯度波形、层面选择方向校正后的逻辑梯度波形。
进一步的,在步骤b3中,对于第n步相位编码,其在相位编码方向的K空间坐标按下式计算:
其中,γ为氢核的磁旋比,GP(t,n)为第n步相位编码校正后的梯度波形,te xcite表示射频激发时刻,techc表示采集到回波信号中心时刻;
频率编码方向的K空间坐标按下式计算:
其中,tm表示采集回波信号中第m个数据点时对应的时刻。
优选的,所述涡流补偿参数是根据涡流的数学物理模型测出。
优选的,所述理想梯度波形由磁共振扫描序列参数得到。
本发明相对于现有技术,有效降低了硬件成本,提高了系统性能,操作简便,成像效果好。
附图说明
图1为梯度磁场波形受涡流影响示意图;
图2为梯度预加重后梯度磁场波形受涡流影响示意图;
图3为基于梯度预加重涡流补偿技术的磁共振成像系统示意图;
图4为本发明涡流补偿方法磁共振成像系统示意图;
图5为本发明K空间轨迹校正及成像流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。
如图4、图5所示,一种磁共振成像系统涡流补偿方法,包括以下步骤:
步骤a、磁共振谱仪根据扫描序列参数产生梯度磁场的波形并输出到梯度放大器,驱动磁共振扫描仪工作,并由磁共振谱仪采集磁共振数据;
步骤b、根据理想梯度波形和涡流补偿参数进行梯度轨迹校正,即计算出受涡流影响后的梯度磁场随时间变化的波形,进而得到校正后的K空间轨迹;
步骤c、将步骤a采集到的磁共振数据按照步骤b得到的校正后的K空间轨迹利于傅里叶变换重新网格化到笛卡尔坐标系中,或者利用非均匀傅里叶变换重建出磁共振图像。
所述涡流补偿参数是根据涡流的数学物理模型测出。
所述理想梯度波形由磁共振扫描序列参数得到。
在步骤b中,还包括以下步骤:
步骤b1、根据X轴、Y轴、Z轴(分别对应三个梯度线圈)的理想梯度波形和涡流补偿参数分别计算出X、Y、Z轴受涡流影响后的实际梯度波形,即可得到校正后的物理梯度波形;
步骤b2、将校正后的物理梯度波形从物理轴映射到逻辑轴,得到校正后的逻辑梯度波形;
步骤b3、根据校正后的逻辑梯度波形,得出校正后的K空间轨迹。
在步骤b1中,X轴受涡流影响后的实际波形计算公式如下:
其中GX(t)、GY(t)、GZ(t)分别为X轴、Y轴、Z轴的理想波形,G′X(t)为X轴受涡流影响的实际波形,为第i个时间常数对应的自相关涡流补偿卷积核函数,为第i个交叉项涡流补偿卷积核函数,即Y轴和Z轴对X轴的影响,卷积核函数表示为如下形式:
其中,分别为系统校正阶段测量到的第i项涡流补偿参数,即根据涡流物理数学模型得到的补偿增益和时间常数,
上述公式计算出了X轴受涡流影响后的实际梯度波形,按照相同原理,分别计算出Y轴和Z轴受涡流影响后的实际梯度波形G′Y(t)和G′z(t),即得到了校正后的物理轴梯度波形。
在步骤b2中,映射关系如下:
其中,R为旋转矩阵,即由切片层面所在坐标系转换到物理坐标系的变换矩阵,R-1为旋转矩阵的逆矩阵,GR(t)、GP(t)、GS(t)分别对应逻辑坐标系中的频率编码方向校正后的逻辑梯度波形、相位编码方向校正后的逻辑梯度波形、层面选择方向校正后的逻辑梯度波形。
在步骤b3中,对于第n步相位编码,其在相位编码方向的K空间坐标按下式计算:
其中,γ为氢核的磁旋比,Gp(t,n)为第n步相位编码校正后的梯度波形,te xcite表示射频激发时刻,techc表示采集到回波信号中心时刻;
频率编码方向的K空间坐标按下式计算:
其中,tm表示采集回波信号中第m个数据点时对应的时刻。
当然,本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (6)

1.一种磁共振成像系统涡流补偿方法,其特征在于:包括以下步骤:
步骤a、磁共振谱仪根据扫描序列参数产生梯度磁场的波形并输出到梯度放大器,驱动磁共振扫描仪工作,并由磁共振谱仪采集磁共振数据;
步骤b、根据理想梯度波形和涡流补偿参数进行梯度轨迹校正,即计算出受涡流影响后的梯度磁场随时间变化的波形,进而得到校正后的K空间轨迹,包括以下步骤:
步骤b1、根据理想梯度波形和涡流补偿参数计算出校正后的物理梯度波形,
步骤b2、将校正后的物理梯度波形从物理轴映射到逻辑轴,得到校正后的逻辑梯度波形,
步骤b3、根据校正后的逻辑梯度波形,得出校正后的K空间轨迹;
步骤c、将步骤a采集到的磁共振数据按照步骤b得到的校正后的K空间轨迹进行网格化重建或者利用非均匀傅里叶变换重建出磁共振图像。
2.根据权利要求1所述的磁共振成像系统涡流补偿方法,其特征在于:在步骤b1中,X轴受涡流影响后的实际波形计算公式如下:
其中GX(t)、GY(t)、GZ(t)分别为X轴、Y轴、Z轴的理想波形,G′X(t)为X轴受涡流影响的实际波形,为第i个时间常数对应的自相关涡流补偿卷积核函数,为第i个交叉项涡流补偿卷积核函数,即Y轴和Z轴对X轴的影响,卷积核函数表示为如下形式:
其中,分别为系统校正阶段测量到的第i项涡流补偿参数,即根据涡流物理数学模型得到的补偿增益和时间常数,
上述公式计算出了X轴受涡流影响后的实际梯度波形,按照相同原理,分别计算出Y轴和Z轴受涡流影响后的实际梯度波形G′Y(t)和G′Z(t),即得到了校正后的物理轴梯度波形。
3.根据权利要求2所述的磁共振成像系统涡流补偿方法,其特征在于:在步骤b2中,映射关系如下:
其中,R为旋转矩阵,即由切片层面所在坐标系转换到物理坐标系的变换矩阵,R-1为旋转矩阵的逆矩阵,GR(t)、GP(t)、GS(t)分别对应逻辑坐标系中的频率编码方向校正后的逻辑梯度波形、相位编码方向校正后的逻辑梯度波形、层面选择方向校正后的逻辑梯度波形。
4.根据权利要求3所述的磁共振成像系统涡流补偿方法,其特征在于:在步骤b3中,对于第n步相位编码,其在相位编码方向的K空间坐标按下式计算:
其中,γ为氢核的磁旋比,GP(t,n)为第n步相位编码校正后的梯度波形,texcite表示射频激发时刻,techo表示采集到回波信号中心时刻;
频率编码方向的K空间坐标按下式计算:
其中,tm表示采集回波信号中第m个数据点时对应的时刻。
5.根据权利要求1所述的磁共振成像系统涡流补偿方法,其特征在于:所述涡流补偿参数是根据涡流的数学物理模型测出。
6.根据权利要求1所述的磁共振成像系统涡流补偿方法,其特征在于:所述理想梯度波形由磁共振扫描序列参数得到。
CN201810162124.8A 2018-02-27 2018-02-27 一种磁共振成像系统涡流补偿方法 Active CN108303665B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201810162124.8A CN108303665B (zh) 2018-02-27 2018-02-27 一种磁共振成像系统涡流补偿方法
PCT/CN2018/080112 WO2019165656A1 (zh) 2018-02-27 2018-03-23 一种磁共振成像系统涡流补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810162124.8A CN108303665B (zh) 2018-02-27 2018-02-27 一种磁共振成像系统涡流补偿方法

Publications (2)

Publication Number Publication Date
CN108303665A CN108303665A (zh) 2018-07-20
CN108303665B true CN108303665B (zh) 2019-12-13

Family

ID=62848754

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810162124.8A Active CN108303665B (zh) 2018-02-27 2018-02-27 一种磁共振成像系统涡流补偿方法

Country Status (2)

Country Link
CN (1) CN108303665B (zh)
WO (1) WO2019165656A1 (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108303665B (zh) * 2018-02-27 2019-12-13 奥泰医疗系统有限责任公司 一种磁共振成像系统涡流补偿方法
CN111329479B (zh) * 2018-12-19 2023-09-29 西门子保健有限责任公司 用于提供mr图像的方法、磁共振成像设备和介质
CN109752681B (zh) * 2019-01-22 2021-10-29 奥泰医疗系统有限责任公司 倾斜平面回波成像方法及核磁共振成像系统
CN109799472B (zh) * 2019-03-08 2022-06-28 沈阳工业大学 一种基于深度学习的磁共振涡流补偿方法
CN112881959B (zh) * 2020-06-10 2021-12-21 苏州纽迈分析仪器股份有限公司 一种用于磁共振成像的梯度涡流补偿方法及系统
CN112630712B (zh) * 2020-11-25 2024-01-02 哈尔滨医科大学 一种简便的核磁共振梯度波形失真预矫正方法
CN113075603B (zh) * 2021-03-25 2022-07-29 上海交通大学 磁共振成像脉冲序列设计方法
CN113296037A (zh) * 2021-05-21 2021-08-24 电子科技大学 一种高场磁共振梯度控制器
CN114124723A (zh) * 2021-11-19 2022-03-01 上海电气(集团)总公司智惠医疗装备分公司 一种磁共振谱仪的组网方法和磁共振谱仪
CN115629347B (zh) * 2022-10-20 2023-09-19 无锡鸣石峻致医疗科技有限公司 一种磁共振成像系统中获得梯度轨迹的方法、装置和介质
CN117233677B (zh) * 2023-09-15 2024-06-04 安徽峻德医疗科技有限公司 非笛卡尔激发采样的磁共振梯度脉冲信号补偿方法及装置

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1427267A (zh) * 2001-12-21 2003-07-02 中国科学院电工研究所 磁共振成像的方法和装置
JP3930439B2 (ja) * 2003-02-06 2007-06-13 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 渦電流補正方法および磁気共鳴撮影装置
JP5611661B2 (ja) * 2009-06-04 2014-10-22 株式会社東芝 磁気共鳴イメージング装置
CN102298129B (zh) * 2011-05-30 2014-07-09 苏州朗润医疗系统有限公司 一种用于核磁共振成像系统的涡流测量及补偿方法
KR101473872B1 (ko) * 2013-02-05 2014-12-18 삼성전자 주식회사 자기공명영상장치 및 그 제어방법
CN104181479B (zh) * 2013-05-23 2015-07-01 上海联影医疗科技有限公司 用于磁共振成像系统的涡流补偿方法
DE102014222495A1 (de) * 2014-11-04 2016-05-04 Siemens Aktiengesellschaft Korrektur von Magnetresonanz-Bildaufnahmen mit Hilfe konvertierter Feldkarten
CN107677976B (zh) * 2017-09-26 2019-07-16 中国科学院武汉物理与数学研究所 一种自适应核磁共振梯度预加重波形产生装置及方法
CN108303665B (zh) * 2018-02-27 2019-12-13 奥泰医疗系统有限责任公司 一种磁共振成像系统涡流补偿方法

Also Published As

Publication number Publication date
CN108303665A (zh) 2018-07-20
WO2019165656A1 (zh) 2019-09-06

Similar Documents

Publication Publication Date Title
CN108303665B (zh) 一种磁共振成像系统涡流补偿方法
JP6037652B2 (ja) 拡散強調磁気共鳴データの生成方法、磁気共鳴システムおよびコンピュータ読み取り可能な記憶媒体
JP5611661B2 (ja) 磁気共鳴イメージング装置
US8890530B2 (en) Magnetic resonance imaging method, magnetic resonance imaging apparatus, and control device of magnetic resonance imaging apparatus
CN107367703B (zh) 磁共振扫描方法、系统、装置及存储介质
DE19821780A1 (de) Korrektur von durch Maxwell-Terme bei einer Schnitt-Verschiebungs-Echo-Planar-Abbildung verursachten Artefakten
DE19801808B4 (de) Vorrichtung und Verfahren zur Verringerung von Maxwell-Term-Artefakten bei schnellen Spinecho-Magnetresonanzbildern
CN112763958B (zh) 一种基于神经网络的多次激发平面回波磁共振成像方法
CN109696646B (zh) 一种快速自旋回波脉冲序列中射频脉冲的优化方法
JPH05192308A (ja) 非対称なnmrエコー取得とともに短いパルスシーケンスを使用する高分解能イメージング装置
DE19814677A1 (de) Korrektur einer durch Maxwell-Terme verursachten Verschlechterung eines Axial-Bild-Signals
CN110320485A (zh) 测量梯度延时和一阶场不均匀度的装置、方法及存储介质
CN104665831A (zh) 图像伪影去除方法及图像伪影去除装置
JP5675044B2 (ja) 磁気共鳴イメージング装置
JP2000023943A (ja) 磁気共鳴方法及び装置
US10809341B1 (en) Readout-segmented echo planar imaging with k-space averaging
CN108577841B (zh) 一种propeller技术中抑制非刚性运动的权重计算方法
US20200326398A1 (en) Phase Correction for Echo-Planar Imaging
US5412321A (en) Magnetic resonance imaging method and apparatus
JPH0622928A (ja) 超高速mrイメージング方法
JPH0838459A (ja) Mr画像再生方法
JP3507586B2 (ja) 磁気共鳴イメージング装置
US4939462A (en) Method and apparatus for MRI
US6753683B2 (en) System and method for generating spectra
US6268729B1 (en) Correction of magnetic resonance signals in a magnetic resonance imaging method, and a magnetic resonance device

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
PP01 Preservation of patent right
PP01 Preservation of patent right

Effective date of registration: 20230905

Granted publication date: 20191213

PD01 Discharge of preservation of patent
PD01 Discharge of preservation of patent

Date of cancellation: 20231117

Granted publication date: 20191213