CN115685032B - 一种幅条k空间下的快速自旋回波的校正方法及系统 - Google Patents
一种幅条k空间下的快速自旋回波的校正方法及系统 Download PDFInfo
- Publication number
- CN115685032B CN115685032B CN202211363046.0A CN202211363046A CN115685032B CN 115685032 B CN115685032 B CN 115685032B CN 202211363046 A CN202211363046 A CN 202211363046A CN 115685032 B CN115685032 B CN 115685032B
- Authority
- CN
- China
- Prior art keywords
- phase
- gradient
- echo
- correction
- readout
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 162
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000005284 excitation Effects 0.000 claims abstract description 46
- 230000009897 systematic effect Effects 0.000 claims abstract description 14
- 238000003384 imaging method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000010304 firing Methods 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种幅条k空间下的快速自旋回波的校正方法及系统,包括:在自旋回波的选层梯度方向的每个预设间隔上选取一个测试平面;根据预设角度旋转第一个测试平面,获取第一读出梯度和第一相位梯度的编码数值;根据自旋回波在第一聚相脉冲后对应形成的相位值,计算修正相位;根据回波的相位差值,获得第一读出梯度的修正值和第一相位梯度的修正值;使用第一聚相激发脉冲的修正相位、第一读出梯度的修正值和第一相位梯度的修正值,分别修正第一回波和第二回波之间的系统误差、所述第一读出梯度和第一相位梯度的编码数值。解决幅条k空间下由于梯度随着辐条角度不停在变化,一个激发脉冲不足以覆盖全部的梯度情况,导致成像出像伪影的问题。
Description
技术领域
本发明涉及核磁共振检测技术,具体涉及一种幅条k空间下的快速自旋回波的校正方法及系统。
背景技术
核磁共振成像中,快速自旋回波类序列,通常利用CPMG条件避免伪影的出现。由于理想的CPMG条件很难达到,所以会有一系列的补偿措施,来确保在回波链中各个级别回波的一致。各个补偿方法中,最灵活适应的通常采用预扫描的方式,在序列实施的选层梯度、读出梯度、相位梯度基本确立的情况下,确定回波之间的差异从而实施补偿。
梯度固定的情况下,误差通常是确定的,因此常规笛卡尔K空间,由于选层、读出梯度不变,针对这些梯度误差的预扫描流程相对简单,通常可以只采集一个激发下的误差来应用到全体k空间。如果采用辐条状k空间采集,由于梯度随着辐条角度不停在变化,一个激发不足以覆盖全部的梯度情况。
发明内容
针对上述问题,本发明提供一种幅条k空间下的快速自旋回波的校正方法,包括:
在自旋回波的选层梯度方向的每个预设间隔上选取一个测试平面;
根据预设角度旋转第一个测试平面,根据第一个测试平面的初始读出梯度、初始相位梯度的编码数值,获取第一读出梯度和第一相位梯度的编码数值;
根据自旋回波在第一聚相脉冲后对应形成的第一回波和第二回波在等中心处的相位值,计算所述第一聚相激发脉冲的修正相位;
根据第一回波和第二回波的相位差值,获得所述第一读出梯度的修正值和第一相位梯度的修正值;
使用所述第一聚相激发脉冲的修正相位、第一读出梯度的修正值和第一相位梯度的修正值,分别修正第一回波和第二回波之间的系统误差、所述第一读出梯度和第一相位梯度的编码数值。
进一步的,还包括:
根据预设角度继续旋转第一个测试平面,根据自旋回波在第二修正激发脉冲后对应形成的第一自旋回波和第二自旋回波在等中心处的相位值,计算所述第二修正激发脉冲的相位;
将第一个测试平面第一次旋转后获得的修正后的读出梯度、相位梯度作为初始读出梯度、初始相位梯度的编码数值,根据所述初始读出梯度、初始相位梯度的编码数值,获取第二读出梯度和第二相位梯度的编码数值;
根据第一回波和第二回波的相位差值,获得第二读出梯度的修正值和第一相位梯度的修正值;
用所述第二修正激发脉冲的相位、第二读出梯度的修正值和第二相位梯度的修正值,分别修正第一回波和第二回波之间的系统误差、修正所述第一个测试平面第二次旋转的读出梯度、相位梯度的编码数值和相位;根据上述步骤,依次修正第一个测试平面每次旋转的读出梯度、相位梯度的编码数值和相位,直至完成所述第一个测试平面所有梯度的修正。
进一步的,在完成所述第一个测试平面所有梯度的修正的步骤之后,还包括:
根据第一个测试平面的对相位、读出梯度和相位梯度的修正步骤,依次修正剩余的测试平面的相位、读出梯度和相位梯度。
进一步的,根据第一回波和第二回波的相位差值,获得第二读出梯度的修正值和第一相位梯度的修正值,包括:
根据第一回波和第二回波的相位差值,获得第二读出梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第二读出梯度的修正值;
根据初始相位梯度和第二相位梯度的相位斜率、磁旋比及初始相位到第二相位的持续时间,计算第二相位梯度的修正值。
进一步的,根据自旋回波在第一修正激发脉冲后对应形成的第一自旋回波和第二自旋回波在等中心处的相位值,计算所述第一修正激发脉冲的相位,包括:
在第一修正激发脉冲180。后,获取第一自旋回波和第二自旋回波在等中心处的相位值Φ1,Φ2,计算所述第一修正激发脉冲的相位Φ,具体为:
Φ=(Φ1-Φ2)/2
进一步的,根据第一回波和第二回波的相位差值,获得所述第一读出梯度的修正值和第一相位梯度的修正值,包括:
根据第一回波和第二回波的相位差值,获得第一读出梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第一读出梯度的修正值;
根据第一回波和第二回波的相位差值,获得第一相位梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第一相位梯度的修正值。
进一步的,根据初始读出梯度和第一读出梯度的相位斜率、磁旋比及初始梯度到第一梯度的持续时间,计算第一读出梯度的修正值,具体的为:
G=S/T
进一步的,根据初始相位梯度和第一相位梯度的相位斜率、磁旋比及初始相位到第一相位的持续时间,计算第一相位梯度的修正值,具体为:
G=S/T
进一步的,还包括:
若所述测试平面是xy平面,随着旋转角度α的改变,读出梯度和相位梯度将由x路和y路梯度合成:
Greadx=Gread*cos(α);Gready=Gread*sin(α);
Gphasex=Gphase*sin(α);Gphasey=Gphase*cos(α)。
本发明同时提供一种幅条k空间下的快速自旋回波的校正系统,包括:
测试平面选取模块,用于在自旋回波的选层梯度方向的每个预设间隔上选取一个测试平面;
第一编码数值获取模块,用于根据预设角度旋转第一个测试平面,根据第一个测试平面的初始读出梯度、初始相位梯度的编码数值,获取第一读出梯度和第一相位梯度的编码数值;
第一修正相位计算模块,用于根据自旋回波在第一聚相脉冲后对应形成的第一回波和第二回波在等中心处的相位值,计算所述第一聚相激发脉冲的修正相位;
第一修正值获取模块,用于根据第一回波和第二回波的相位差值,获得所述第一读出梯度的修正值和第一相位梯度的修正值;
第一修正模块,用于使用所述第一聚相激发脉冲的修正相位、第一读出梯度的修正值和第一相位梯度修正值,分别修正第一回波和第二回波之间的系统误差、修正所述第一读出梯度和第一相位梯度的编码数值。
本发明提供的一种幅条k空间下的快速自旋回波的校正方法及系统,构建幅条k空间下的快速自旋回波的校正模型,每个激发脉冲后,会重新计算梯度的补偿,解决幅条k空间下由于梯度随着辐条角度不停在变化,一个激发脉冲不足以覆盖全部的梯度情况,导致成像出像伪影的问题。
附图说明
图1是本发明提供的一种幅条k空间下的快速自旋回波的校正方法的流程示意图;
图2是本发明涉及的常规自旋回波的序列时序图;
图3是本发明涉及的修正梯度及修正相位的位置图;
图4是本发明涉及的幅条k空间成像平面的旋转角度的示意图;
图5是本发明涉及的幅条k空间下的快速自旋回波的校正模型的构建流程图;
图6是本发明涉及的一种幅条k空间下的快速自旋回波的校正系统的结构示意图。
具体实施方式
在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施的限制。
常规自旋回波的序列时序图如图2所示,在每一个180聚焦脉冲1后,形成一个自旋回波2,通过改变相位编码梯度3的幅度,采集K空间不同位置的编码信息。
现实中,CPMG条件很难满足,因此回波2以及后续回波存在着相位误差,这种误差会不可避免的在最终图像上出现伪影。
为了消除这种影响,本发明提供一种幅条k空间下的快速自旋回波的校正方法,如图1所示,包括如下步骤:
步骤S101,在自旋回波的选层梯度方向的每个预设间隔上选取一个测试平面。
步骤S102,根据预设角度旋转第一个测试平面,获取第一读出梯度和第一相位梯度的编码数值。
据预设角度旋转第一个测试平面,根据初始读出梯度、初始相位梯度的的编码数值,获得旋转后第一读出梯度和第一相位梯度的编码数值。
步骤S103,根据自旋回波在第一聚相脉冲后对应形成的第一回波和第二回波在等中心处的相位值,计算所述第一聚相激发脉冲的修正相位。
在图3所示的第一修正激发脉冲180。后,获取第一自旋回波和第二自旋回波在等中心处的相位值Φ1,Φ2,计算所述第一修正激发脉冲的相位Φ,具体为:
Φ=(Φ1-Φ2)/2
步骤S104,根据第一回波和第二回波的相位差值,获得所述第一读出梯度的修正值和第一相位梯度的修正值。
根据第一回波和第二回波的相位差值,获得第一读出梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第一读出梯度的修正值,具体的,根据初始读出梯度和第一读出梯度的相位斜率、磁旋比及初始梯度到第一梯度的持续时间,计算第一读出梯度的修正值,如图3中箭头4所示的梯度,具体的为:
G=S/T
根据第一回波和第二回波的相位差值,获得第一相位梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第一相位梯度的修正值,具体的,根据初始相位梯度和第一相位梯度的相位斜率、磁旋比及初始相位到第一相位的持续时间,计算第一相位梯度的修正值,如图3中箭头5所示的梯度,具体为:
G=S/T
第一读出梯度的修正值和第一相位梯度的修正值和第一修正激发脉冲的相位Φ的数值会随着选层位置,梯度大小的变化而变化。
步骤S105,使用所述第一聚相激发脉冲的修正相位、第一读出梯度的修正值和第一相位梯度的修正值,分别修正第一回波和第二回波之间的系统误差,修正所述第一读出梯度和第一相位梯度的编码数值。
在修正系统误差和第一读出梯度和第一相位梯度的编码数值后,根据预设角度继续旋转第一个测试平面,根据自旋回波在第二修正激发脉冲后对应形成的第一自旋回波和第二自旋回波在等中心处的相位值,计算所述第二修正激发脉冲的相位;
将第一个测试平面第一次旋转后获得的修正后的读出梯度、相位梯度作为初始读出梯度、初始相位梯度的编码数值,根据所述初始读出梯度、初始相位梯度的编码数值,获取第二读出梯度和第二相位梯度的编码数值;
根据第一回波和第二回波的相位差值,获得第二读出梯度的修正值和第一相位梯度的修正值;根据第一回波和第二回波的相位差值,获得第二读出梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第二读出梯度的修正值;根据初始相位梯度和第二相位梯度的相位斜率、磁旋比及初始相位到第二相位的持续时间,计算第二相位梯度的修正值。
用所述第二修正激发脉冲的相位、第二读出梯度的修正值和第二相位梯度的修正值,分别修正第一回波和第二回波之间的系统误差、修正所述第一个测试平面第二次旋转的读出梯度、相位梯度的编码数值和相位;根据上述步骤,依次修正第一个测试平面每次旋转的读出梯度、相位梯度的编码数值和相位,直至完成所述第一个测试平面所有梯度的修正。
在完成所述第一个测试平面所有梯度的修正后,根据第一个测试平面的对相位、读出梯度和相位梯度的修正步骤,依次修正剩余的测试平面的相位、读出梯度和相位梯度。从而构建幅条k空间下的快速自旋回波的校正模型。
进一步的,若所述测试平面是xy平面,随着旋转角度α的改变,读出梯度和相位梯度将由x路和y路梯度合成,如图4所示:
Greadx=Gread*cos(α);Gready=Gread*sin(α);
Gphasex=Gphase*sin(α);Gphasey=Gphase*cos(α)。
梯度输出的变化,导致系统状态的改变,修改正梯度4、5和修改正相位也会同时随着变化。为满足奈奎斯特采样定律,辐条k空间相位编码数为笛卡尔坐标系下的π/2倍,这意味了通过扫描时先预扫描计算修正值所花费的时间非常长。
幅条k空间下的快速自旋回波的校正模型的构建流程如图5所示,沿选层梯度方向每隔定长选取一个测试平面,确定读出梯度和相位梯度的编码数值,并按不大于10°的角度旋转,获取回波信号,根据回波信号计算对应的读出梯度的修正值G4,相位梯度的修正值G5,以及修正相位Φ,然后,使用上述三个修正值修正读出梯度和相位梯度的编码数值,以及回波的系统误差,接着,继续按不大于10°的角度旋转第一个测试平面,确定第二次旋转的读出梯度和相位梯度的编码数值,然后计算对应的读出梯度的修正值G4,相位梯度的修正值G5,以及修正相位Φ,使用上述三个修正值修正第二次旋转获得的读出梯度和相位梯度的编码数值,以及回波的系统误差,直至第一个测试平面的所有梯度修正完成。根据第一个测试平面的对相位、读出梯度和相位梯度的修正步骤,依次修正剩余的测试平面的相位、读出梯度和相位梯度。至此,校正模型构建完成。
基于同一发明构思,本发明同时提供一种条k空间下的快速自旋回波的校正系统,如图6所示,包括:
测试平面选取模块610,用于在自旋回波的选层梯度方向的每个预设间隔上选取一个测试平面;
第一编码数值获取模块630,用于根据预设角度旋转第一个测试平面,根据第一个测试平面的初始读出梯度、初始相位梯度的编码数值,获取第一读出梯度和第一相位梯度的编码数值;
修正相位计算模块640,用于根据自旋回波在第一聚相脉冲后对应形成的第一回波和第二回波在等中心处的相位值,计算所述第一聚相激发脉冲的修正相位;
第一修正值计算模块650,用于根据第一回波和第二回波的相位差值,获得所述第一读出梯度的修正值和第一相位梯度的修正值;
第一修正模块660,用于使用所述第一聚相激发脉冲的修正相位、第一读出梯度的修正值和第一相位梯度修正值,分别修正第一回波和第二回波之间的系统误差、修正所述第一读出梯度和第一相位梯度的编码数值。
进一步的,还包括:
第二修正相位计算模块,用于根据预设角度继续旋转第一个测试平面,根据自旋回波在第二修正激发脉冲后对应形成的第一自旋回波和第二自旋回波在等中心处的相位值,计算所述第二修正激发脉冲的相位;
第二修正值获取模块,用于将第一个测试平面第一次旋转后获得的修正后的读出梯度、相位梯度作为初始读出梯度、初始相位梯度的编码数值,根据所述初始读出梯度、初始相位梯度的编码数值,获取第二读出梯度和第二相位梯度的编码数值;
第二修正值获取模块,用于根据第一回波和第二回波的相位差值,获得第二读出梯度的修正值和第一相位梯度的修正值;
第二修正模块,用于使用所述第二修正激发脉冲的相位、第二读出梯度的修正值和第二相位梯度的修正值,分别修正第一回波和第二回波之间的系统误差、修正所述第一个测试平面第二次旋转的读出梯度、相位梯度的编码数值和相位;根据上述步骤,依次修正第一个测试平面每次旋转的读出梯度、相位梯度的编码数值和相位,直至完成所述第一个测试平面所有梯度的修正。
进一步的,还包括:
修正模块,用于根据第一个测试平面的对相位、读出梯度和相位梯度的修正步骤,依次修正剩余的测试平面的相位、读出梯度和相位梯度。
进一步的,第一修正值获取模块,包括:
第二读出梯度修正值获取子模块,用于根据第一回波和第二回波的相位差值,获得第二读出梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第二读出梯度的修正值;
第二相位梯度修正值获取子模块,用于根据初始相位梯度和第二相位梯度的相位斜率、磁旋比及初始相位到第二相位的持续时间,计算第二相位梯度的修正值。
进一步的,第一修正相位计算模块,包括:
在第一修正激发脉冲180。后,获取第一自旋回波和第二自旋回波在等中心处的相位值Φ1,Φ2,计算所述第一修正激发脉冲的相位Φ,具体为:
Φ=(Φ1-Φ2)/2
进一步的,第一修正值获取模块,包括:
第一读出梯度计算子模块,用于根据第一回波和第二回波的相位差值,获得第一读出梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第一读出梯度的修正值;
第一相位梯度计算子模块,用于根据第一回波和第二回波的相位差值,获得第一相位梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第一相位梯度的修正值。
进一步的,根据初始读出梯度和第一读出梯度的相位斜率、磁旋比及初始梯度到第一梯度的持续时间,计算第一读出梯度的修正值,具体的为:
G=S/T
进一步的,根据初始相位梯度和第一相位梯度的相位斜率、磁旋比及初始相位到第一相位的持续时间,计算第一相位梯度的修正值,具体为:
G=S/T
进一步的,还包括:
梯度合成模块,若所述测试平面是xy平面,随着旋转角度α的改变,读出梯度和相位梯度将由x路和y路梯度合成:
Greadx=Gread*cos(α);Gready=Gread*sin(α);
Gphasex=Gphase*sin(α);Gphasey=Gphase*cos(α)。
本发明提供的一种幅条k空间下的快速自旋回波的校正方法及系统,构建幅条k空间下的快速自旋回波的校正模型,每个激发脉冲后,会重新计算梯度的补偿,解决幅条k空间下由于梯度随着辐条角度不停在变化,一个激发脉冲不足以覆盖全部的梯度情况,导致成像出像伪影的问题。
最后应该说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。
Claims (9)
1.一种幅条k空间下的快速自旋回波的校正方法,其特征在于,包括:
在自旋回波的选层梯度方向的每个预设间隔上选取一个测试平面;
根据预设角度旋转第一个测试平面,根据第一个测试平面的初始读出梯度、初始相位梯度的编码数值,获取第一读出梯度和第一相位梯度的编码数值;
根据预设角度继续旋转第一个测试平面,根据自旋回波在第二修正激发脉冲后对应形成的第一自旋回波和第二自旋回波在等中心处的相位值,计算所述第二修正激发脉冲的相位;
将第一个测试平面第一次旋转后获得的修正后的读出梯度、相位梯度作为初始读出梯度、初始相位梯度的编码数值,根据所述初始读出梯度、初始相位梯度的编码数值,获取第二读出梯度和第二相位梯度的编码数值;
根据第一回波和第二回波的相位差值,获得第二读出梯度的修正值和第一相位梯度的修正值;
用所述第二修正激发脉冲的相位、第二读出梯度的修正值和第二相位梯度的修正值,分别修正第一回波和第二回波之间的系统误差、修正所述第一个测试平面第二次旋转的读出梯度、相位梯度的编码数值和相位;根据上述步骤,依次修正第一个测试平面每次旋转的读出梯度、相位梯度的编码数值和相位,直至完成所述第一个测试平面所有梯度的修正;
根据自旋回波在第一聚相激发脉冲后对应形成的第一回波和第二回波在等中心处的相位值,计算所述第一聚相激发脉冲的修正相位;
根据第一回波和第二回波的相位差值,获得所述第一读出梯度的修正值和第一相位梯度的修正值;
使用所述第一聚相激发脉冲的修正相位、第一读出梯度的修正值和第一相位梯度的修正值,分别修正第一回波和第二回波之间的系统误差、所述第一读出梯度和第一相位梯度的编码数值。
2.根据权利要求1所述的方法,其特征在于,在完成所述第一个测试平面所有梯度的修正的步骤之后,还包括:
根据第一个测试平面的对相位、读出梯度和相位梯度的修正步骤,依次修正剩余的测试平面的相位、读出梯度和相位梯度。
3.根据权利要求1所述的方法,其特征在于,根据第一回波和第二回波的相位差值,获得第二读出梯度的修正值和第一相位梯度的修正值,包括:
根据第一回波和第二回波的相位差值,获得第二读出梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第二读出梯度的修正值;
根据初始相位梯度和第二相位梯度的相位斜率、磁旋比及初始相位到第二相位的持续时间,计算第二相位梯度的修正值。
4.根据权利要求1所述的方法,其特征在于,根据自旋回波在第一修正激发脉冲后对应形成的第一自旋回波和第二自旋回波在等中心处的相位值,计算所述第一修正激发脉冲的相位,包括:
在第一修正激发脉冲180°后,获取第一自旋回波和第二自旋回波在等中心处的相位值Φ1,Φ2,计算所述第一修正激发脉冲的相位Φ,具体为:
Φ=(Φ1-Φ2)/2。
5.根据权利要求1所述的方法,其特征在于,根据第一回波和第二回波的相位差值,获得所述第一读出梯度的修正值和第一相位梯度的修正值,包括:
根据第一回波和第二回波的相位差值,获得第一读出梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第一读出梯度的修正值;
根据第一回波和第二回波的相位差值,获得第一相位梯度的修正梯度,根据所述修正梯度的持续时间,计算所述第一相位梯度的修正值。
8.根据权利要求1所述的方法,其特征在于,还包括:
若所述测试平面是xy平面,随着旋转角度α的改变,读出梯度和相位梯度将由x路和y路梯度合成:
Greadx=Gread*cos(α);Gready=Gread*sin(α);
Gphasex=Gphase*sin(α);Gphasey=Gphase*cos(α)。
9.一种幅条k空间下的快速自旋回波的校正系统,其特征在于,包括:
测试平面选取模块,用于在自旋回波的选层梯度方向的每个预设间隔上选取一个测试平面;
第一编码数值获取模块,用于根据预设角度旋转第一个测试平面,根据第一个测试平面的初始读出梯度、初始相位梯度的编码数值,获取第一读出梯度和第一相位梯度的编码数值;根据预设角度继续旋转第一个测试平面,根据自旋回波在第二修正激发脉冲后对应形成的第一自旋回波和第二自旋回波在等中心处的相位值,计算所述第二修正激发脉冲的相位;
将第一个测试平面第一次旋转后获得的修正后的读出梯度、相位梯度作为初始读出梯度、初始相位梯度的编码数值,根据所述初始读出梯度、初始相位梯度的编码数值,获取第二读出梯度和第二相位梯度的编码数值;
根据第一回波和第二回波的相位差值,获得第二读出梯度的修正值和第一相位梯度的修正值;
用所述第二修正激发脉冲的相位、第二读出梯度的修正值和第二相位梯度的修正值,分别修正第一回波和第二回波之间的系统误差、修正所述第一个测试平面第二次旋转的读出梯度、相位梯度的编码数值和相位;根据上述步骤,依次修正第一个测试平面每次旋转的读出梯度、相位梯度的编码数值和相位,直至完成所述第一个测试平面所有梯度的修正;
第一修正相位计算模块,用于根据自旋回波在第一聚相激发脉冲后对应形成的第一回波和第二回波在等中心处的相位值,计算所述第一聚相激发脉冲的修正相位;
第一修正值获取模块,用于根据第一回波和第二回波的相位差值,获得所述第一读出梯度的修正值和第一相位梯度的修正值;
第一修正模块,用于使用所述第一聚相激发脉冲的修正相位、第一读出梯度的修正值和第一相位梯度修正值,分别修正第一回波和第二回波之间的系统误差、修正所述第一读出梯度和第一相位梯度的编码数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211363046.0A CN115685032B (zh) | 2022-11-02 | 2022-11-02 | 一种幅条k空间下的快速自旋回波的校正方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211363046.0A CN115685032B (zh) | 2022-11-02 | 2022-11-02 | 一种幅条k空间下的快速自旋回波的校正方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115685032A CN115685032A (zh) | 2023-02-03 |
CN115685032B true CN115685032B (zh) | 2023-06-23 |
Family
ID=85047310
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211363046.0A Active CN115685032B (zh) | 2022-11-02 | 2022-11-02 | 一种幅条k空间下的快速自旋回波的校正方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115685032B (zh) |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07323021A (ja) * | 1994-05-31 | 1995-12-12 | Shimadzu Corp | Mrイメージング装置 |
JP2000157509A (ja) * | 1998-11-23 | 2000-06-13 | General Electric Co <Ge> | Mriシステムの較正方法、mriシステム及びリセット勾配発生方法 |
US6369568B1 (en) * | 1999-06-03 | 2002-04-09 | Ge Medical Systems Global Technology Company, Llc | Fast spin echo phase correction for MRI system |
JP3875479B2 (ja) * | 2000-10-20 | 2007-01-31 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴撮影装置 |
JP4122452B2 (ja) * | 2001-10-15 | 2008-07-23 | 株式会社日立メディコ | 磁気共鳴撮像装置 |
US6617850B2 (en) * | 2001-12-03 | 2003-09-09 | Mayo Foundation For Medical Education And Research | Motion correction of magnetic resonance images using phase difference of two orthogonal acquisitions |
WO2010074059A1 (ja) * | 2008-12-26 | 2010-07-01 | 株式会社 日立メディコ | 磁気共鳴イメージング装置及び読み出し傾斜磁場誤差補正方法 |
CN102096054B (zh) * | 2010-12-22 | 2013-03-13 | 华东师范大学 | 一种快速自旋回波脉冲序列的梯度系统延时校正方法 |
CN103278785B (zh) * | 2013-05-30 | 2015-08-19 | 华东师范大学 | 一种快速自旋回波脉冲序列中射频脉冲相位的优化方法 |
CN107024670B (zh) * | 2016-01-29 | 2020-08-04 | 上海联影医疗科技有限公司 | 磁共振系统的校正方法及装置 |
CN109946632B (zh) * | 2019-03-06 | 2023-10-20 | 佛山瑞加图医疗科技有限公司 | 剩磁相位差测量方法及剩磁测量方法 |
CN111352055B (zh) * | 2020-03-17 | 2022-05-06 | 无锡鸣石峻致医疗科技有限公司 | 一种前瞻性相位校正平面回波成像方法 |
CN114526690A (zh) * | 2022-01-19 | 2022-05-24 | 中科融合感知智能研究院(苏州工业园区)有限公司 | 一种用于结构光三维测量的相位校正方法、系统及介质 |
-
2022
- 2022-11-02 CN CN202211363046.0A patent/CN115685032B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN115685032A (zh) | 2023-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3330007B2 (ja) | 核共鳴信号の位相補正装置 | |
US10634753B2 (en) | MR imaging with motion detection | |
US7323873B2 (en) | Phase error measuring method, MR imaging method, and MRI system | |
US4833407A (en) | Scan time reduction using conjugate symmetry and recalled echo | |
US8154294B2 (en) | Magnetic resonance imaging apparatus and method for correcting error due to gradient magnetic field | |
US7847546B2 (en) | Magnetic resonance imaging apparatus | |
US8427153B2 (en) | Method for motion correction in magnetic resonance imaging using radio frequency coil arrays | |
US7619410B2 (en) | System and method for accelerated magnetic resonance parallel imaging | |
US8131047B2 (en) | Method and device for distortion correction in magnetic resonance imaging | |
KR101694610B1 (ko) | 자기 공명 장치를 이용하여 파라미터 맵을 획득하기 위한 방법 및 자기 공명 장치 | |
US20110093233A1 (en) | Through-time radial grappa calibration | |
US8497681B2 (en) | Magnetic resonance method and apparatus to determine phase correction parameters | |
US10393847B2 (en) | Method and apparatus for recording calibration data for a GRAPPA magnetic resonance imaging algorithm | |
US7154269B1 (en) | Iterative method for correction of geometric distortion resulting from phase evolution during segmented echo planar nuclear magnetic resonance imaging and apparatus therefor | |
CN109814058A (zh) | Caipirinha波space | |
US8542012B2 (en) | Through-time non-cartesian grappa calibration | |
CN113298902B (zh) | 卷褶视野磁共振图像的重建方法、计算机设备及存储介质 | |
JP2020523075A (ja) | 可変コントラストのスタック・オブ・スター収集を使用したmrイメージング | |
CN102841329B (zh) | 磁共振信号处理方法及装置 | |
US5680045A (en) | Grase-type MR pulse sequences | |
CN115685032B (zh) | 一种幅条k空间下的快速自旋回波的校正方法及系统 | |
CN113298901B (zh) | 卷褶视野磁共振图像的重建方法、计算机设备及存储介质 | |
McLeish et al. | Free‐breathing radial acquisitions of the heart | |
CN111142056B (zh) | 平面回波二维空间选择性脉冲的校正方法 | |
CN110412493B (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 |