CN108697369B - 磁共振成像设备及其方法 - Google Patents

磁共振成像设备及其方法 Download PDF

Info

Publication number
CN108697369B
CN108697369B CN201680082407.0A CN201680082407A CN108697369B CN 108697369 B CN108697369 B CN 108697369B CN 201680082407 A CN201680082407 A CN 201680082407A CN 108697369 B CN108697369 B CN 108697369B
Authority
CN
China
Prior art keywords
echo
data
space
echo data
mri apparatus
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
CN201680082407.0A
Other languages
English (en)
Other versions
CN108697369A (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.)
Samsung Electronics Co Ltd
Original Assignee
Samsung Electronics Co Ltd
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 Samsung Electronics Co Ltd filed Critical Samsung Electronics Co Ltd
Publication of CN108697369A publication Critical patent/CN108697369A/zh
Application granted granted Critical
Publication of CN108697369B publication Critical patent/CN108697369B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/56545Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by finite or discrete sampling, e.g. Gibbs ringing, truncation artefacts, phase aliasing artefacts
    • 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/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • 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/443Assessment of an electric or a magnetic field, e.g. spatial mapping, determination of a B0 drift or dosimetry
    • 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/561Image 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/5611Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • 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/561Image 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/5615Echo 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]
    • 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/561Image 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/5615Echo 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/5616Echo 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
    • 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/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/482MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
    • 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/56554Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in EPI
    • 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/56563Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)

Abstract

提供了一种用于基于多回波序列获取磁共振(MR)图像的磁共振成像(MRI)设备和所述MRI设备的方法。所述MRI设备包括:数据获取器,被配置为基于在第一回波时间产生的回波获取第一回波数据,并且基于在比第一回波时间晚的第二回波时间产生的回波获取第二回波数据,其中,第一回波数据包括k空间中与包括在第二回波数据中的部分重叠的部分。所述MRI设备还包括:图像处理器,被配置为基于第一回波数据和第二回波数据重建MR图像。

Description

磁共振成像设备及其方法
技术领域
与示例性实施例一致的设备和方法涉及磁共振成像(MRI)设备及其方法,更具体地,涉及用于通过使用多回波序列来获取磁共振(MR)图像的MRI设备以及所述MRI设备的方法。
背景技术
用于通过使用磁场对对象进行成像的磁共振成像(MRI)设备可以以期望的角度显示骨骼、腰椎间盘、关节、神经韧带、心脏等的立体图像。
MRI设备的优点在于MRI设备是非侵入性的,与计算机断层扫描(CT)设备相比表现出优异的组织对比度,并且没有由于骨组织而导致的伪像。而且,因为MRI设备可以在不移动对象的情况下以期望的方向捕捉各种横截面图像,所以MRI设备被广泛地与其它成像设备一起使用。
MRI设备可通过使用多回波序列来获取k空间数据。详细地,当使用多回波序列时,可通过激励一个射频(RF)脉冲然后通过使用多个产生的回波获取MR信号来减少捕捉图像所耗费的时间。
当MRI设备使用多回波序列时,MRI设备可使用通过在一个重复时间(TR)期间仅施加一次用于相位编码的梯度磁场来获取与多个回波时间分别相应的多条k空间数据的方法。
可选地,MRI设备可使用通过使用在一个TR期间产生的多个回波来获取一条k空间数据的方法。该方法被称为平面回波成像(EPI)方法。
多回波序列可基于梯度回波或自旋回波。如果多回波序列基于梯度回波,则由于读出的梯度磁场的符号必须连续反转,因此MRI设备可具有高性能来产生梯度磁场。随着近来已经改进了MRI设备的性能,不仅当使用梯度回波时可以减小TR,而且可获取具有期望对比度的图像。
发明内容
技术问题
示例性实施例可至少解决以上问题和/或缺点以及上面没有描述的其它缺点。此外,示例性实施例不必需克服上述缺点,并且可不克服上述任何问题。
示例性实施例可通过获取关于k空间中的在不同回波时间获取的第一回波数据和第二回波数据之间的重叠部分的附加数据来获取磁共振(MR)图像以提高重建的最终图像的质量。
示例性实施例可通过对在不同回波时间获取的第一回波数据和第二回波数据执行相位校正来减少在使用在不同回波时间获取的多条数据时在最终图像中可能发生的模糊或混叠的影响。
技术方案
根据示例性实施例的一方面,提供了一种用于基于多回波序列获取磁共振(MR)图像的磁共振成像(MRI)设备,MRI设备包括:数据获取器,被配置为基于在第一回波时间产生的回波获取第一回波数据,并且基于在比第一回波时间晚的第二回波时间产生的回波获取第二回波数据,其中,第一回波数据包括k空间中与包括在第二回波数据中的部分重叠的部分。MRI设备还包括:图像处理器,被配置为基于第一回波数据和第二回波数据来重建MR图像。
图像处理器还可被配置为:针对第一回波数据和第二回波数据中的一个或两个执行相位校正,并基于被执行了相位校正的第一回波数据和第二回波数据中的一个或两个来重建MR图像。
第一回波数据和第二回波数据中的一个或每一个可包括k空间的中心部分的数据。
图像处理器还可被配置为基于第一回波数据和第二回波数据中的k空间的中心部分的数据的相位来执行相位校正。
MRI设备还可包括:梯度磁场控制器,被配置为基于多回波序列控制施加的梯度磁场,并且在第一回波时间施加的读出梯度磁场的符号可与在第二回波时间施加的读出梯度磁场的符号相同。
数据获取器还可被配置为:响应于第一回波数据包括k空间的中心部分的数据并且第二回波数据不包括k空间的中心部分的数据,基于在第一回波时间产生的回波获取包括在第一回波数据中的所述重叠的部分的附加数据。
数据获取器还可被配置为:响应于第二回波数据包括k空间的中心部分的数据并且第一回波数据不包括k空间的中心部分的数据,基于在第二回波时间产生的回波获取包括在第二回波数据中的所述重叠的部分的附加数据。
数据获取器还可被配置为:确定将被获取的MR图像的特征,并且基于确定的MR图像的特征确定第一回波数据或第二回波数据将包括k空间的中心部分的数据。
图像处理器还可被配置为:基于第三回波数据和第四回波数据产生B0图,并基于产生的B0图对第一回波数据和第二回波数据中的一个或两个执行相位校正,其中,所述第三回波数据和第四回波数据是分别基于在第三回波时间产生的回波和在第四回波时间产生的回波而获取的。
MRI设备还可包括:梯度磁场控制器,被配置为基于多回波序列控制施加的梯度磁场,并且在第三回波时间施加的读出梯度磁场的符号可与在第四回波时间施加的读出梯度磁场的符号相同。
第一回波时间、第二回波时间、第三回波时间和第四回波时间可被包括在重复时间段中。
图像处理器还可被配置为:基于k空间中的除了包括在第一回波数据中的所述重叠的部分之外的第一回波数据的另一部分和第二回波数据来重建k空间数据,对包括在重建的k空间数据中的第二回波数据执行相位校正,并基于包括在第一回波数据中的所述重叠的部分和被执行了相位校正的第二回波数据来重建k空间数据。
根据另一示例性实施例的一方面,提供了一种用于基于多回波序列获取磁共振(MR)图像的磁共振成像(MRI)设备的方法,所述方法包括:基于在第一回波时间产生的回波获取第一回波数据,并且基于在比第一回波时间晚的第二回波时间产生的回波获取第二回波数据,其中,第一回波数据包括k空间中与包括在第二回波数据中的部分重叠的部分。所述方法还包括基于第一回波数据和第二回波数据重建MR图像。
所述方法还可包括:对第一回波数据和第二回波数据中的一个或两个执行相位校正,并且所述重建的步骤可包括基于第一回波数据和第二回波数据中的被执行了相位校正的一个或两个来重建MR图像。
第一回波数据和第二回波数据中的一个或每一个可包括k空间的中心部分的数据。
所述执行的步骤可包括:基于第一回波数据和第二回波数据中的k空间的中心部分的数据的相位来执行相位校正。
所述方法还可包括:基于多回波序列控制施加的梯度磁场,并且在第一回波时间施加的读出梯度磁场的符号可与在第二个回波时间施加的读出梯度磁场的符号相同。
所述方法还可包括:响应于第一回波数据包括k空间的中心部分的数据和第二回波数据不包括k空间的中心部分的数据,基于在第一个回波时间产生的回波,获取包括在第一回波数据中的所述重叠的部分的附加数据。
所述方法还可包括:响应于第二回波数据包括k空间的中心部分的数据和第一回波数据不包括k空间的中心部分的数据,基于在第二回波时间产生的回波获取包括在第二回波数据中的所述重叠的部分的附加数据。
所述方法还可包括:确定将被获取的MR图像的特征,并基于确定的MR图像的特征,确定第一回波数据或第二回波数据将包括k空间的中心部分的数据。
所述方法还可包括基于第三回波数据和第四回波数据产生B0图,并且基于所产生的B0图,针对第一回波数据和第二回波数据中的一个或两个执行相位校正,其中,第三回波数据和第四回波数据是分别基于在第三回波时间产生的回波和在第四回波时间产生的回波而获取的。
所述方法还可包括:基于多回波序列控制施加的梯度磁场,并且在第三回波时间施加的读出梯度磁场的符号可与在第四个回波时间施加的读出梯度磁场的符号相同。
第一回波时间、第二回波时间、第三回波时间和第四回波时间可被包括在重复时间段中。
所述重建可包括:基于k空间中除了包括在的第一回波数据中的所述重叠的部分之外的第一回波数据的另一部分和第二回波数据来重建k空间数据,对包括在重建的k空间数据中的第二回波数据执行相位校正,并基于包括在第一回波数据中的所述重叠的部分和被执行了相位校正的第二回波数据重新重建k空间数据。
一种非暂时性计算机可读存储介质可存储用于使计算机执行所述方法的程序。
根据另一示例性实施例的一方面,提供了一种用于基于多回波序列获取磁共振(MR)图像的磁共振成像(MRI)设备,MRI设备包括:数据获取器,被配置为:基于在第一回波时间产生的回波获取第一回波数据,基于在比第一回波时间晚的第二回波时间产生的回波获取第二回波数据,并基于在比第二回波时间晚的第三回波时间产生的回波获取第三回波数据,第二回波数据包括k空间中与包括在第三回波数据中的部分重叠的部分。MRI设备还包括:图像处理器,被配置为基于第一回波数据、第二回波数据和第三回波数据来重建MR图像。
图像处理器还可被配置为:针对第三回波数据执行相位校正;基于第二回波数据和被执行了相位校正的第三回波数据获取k空间数据,其中,第一回波数据包括k空间中与包括在k空间数据中的部分重叠的部分;针对k空间数据执行相位校正,并基于第一回波数据和被执行了相位校正的k空间数据重建k空间数据。
图像处理器还可被配置为:针对第一回波数据执行相位校正;基于第二回波数据和被执行了相位校正的第一回波数据获取k空间数据,其中,第三回波数据包括k空间中与包括在k空间数据中的部分重叠的部分;对k空间数据执行相位校正,并基于第三回波数据和被执行了相位校正的k空间数据重建k空间数据。
附图说明
从以下结合附图进行的示例性实施例的描述中,以上和/或其它方面将变得显而易见并且更容易理解,其中:
图1是示出根据示例性实施例的基于多回波序列由MRI设备执行以获取MR图像的处理的示图;
图2A是示出根据示例性实施例的基于多回波序列由MRI设备执行以获取k空间数据的过程的曲线图;
图2B是示出根据示例性实施例的基于梯度回波序列由MRI设备执行以获取MR图像的过程的曲线图;
图3A是根据示例性实施例的MRI设备的框图;
图3B是根据另一示例性实施例的MRI设备的框图;
图4A是根据示例性实施例的由图3A或图3B的MRI设备施加的回波链长度(ETL)为2的脉冲序列的脉冲序列图;
图4B是示出根据示例性实施例的根据图4A的脉冲序列图由图3A或图3B的MRI设备执行以获取MR图像的处理的示图;
图5A是根据示例性实施例的由图3A或图3B的MRI设备施加的ETL为2的脉冲序列的脉冲序列图;
图5B是示出根据示例性实施例的根据图5A的脉冲序列由图3A或图3B的MRI设备执行以获取MR图像的处理的示图;
图6A是根据示例性实施例的由图3A或图3B的MRI设备施加的ETL为3的脉冲序列的脉冲序列图;
图6B是示出根据示例性实施例的根据图6A的脉冲序列图由图3A或图3B的MRI设备执行以获取MR图像的处理的示图;
图7A是根据示例性实施例的由图3A或图3B的MRI设备施加的ETL为2的脉冲序列的脉冲序列图;
图7B是示出根据示例性实施例的根据图7A的脉冲序列图由图3A或图3B的MRI设备执行以获取MR图像的处理的示图;
图8A是根据示例性实施例的由图3A或图3B的MRI设备施加的ETL为3的脉冲序列的脉冲序列图;
图8B是示出根据示例性实施例的根据图8A的脉冲序列图由图3A或图3B的MRI设备执行以获取MR图像的处理的示图;
图9A是示出由图3A或图3B的MRI设备获取的第一回波数据和第二回波数据的示图;
图9B是示出根据示例性实施例的通过将广义自动校准部分并行采集(GRAPPA)方法应用于由图3A或图3B的MRI设备获取的第一回波数据和第二回波数据来重建k空间数据的处理的示图;
图10A是示出使用多频带方法由图3A或图3B的MRI设备获取的第一回波数据和第二回波数据的示图;
图10B是示出根据示例性实施例通过将GRAPPA方法应用于使用多频带方法由图3A或图3B的MRI设备获取的第一回波数据和第二回波数据来重建k空间数据的处理的示图;
图11是示出根据示例性实施例的MRI设备的方法的流程图;
图12是示出根据示例性实施例的MRI设备的另一方法的流程图;
图13是示出根据示例性实施例的MRI设备的另一方法的流程图;
图14是示出根据示例性实施例的MRI设备的另一方法的流程图;
图15是一般MRI系统的框图;以及
图16是根据示例性实施例的通信接口的框图。
具体实施方式
以下参照附图更详细地描述示例性实施例。
在以下描述中,即使在不同的附图中,相同的标号也用于相同的元件。说明书中定义的事项(诸如,详细的构造和元件)被提供以帮助全面理解示例性实施例。然而,显而易见的是,可以在没有那些具体定义的事项的情况下实践示例性实施例。而且,公知功能或构造可不被详细描述,因为它们会以不必要的细节模糊描述。
在下文中,将简要描述说明书中使用的术语,然后将详细描述本公开。
本说明书中使用的术语是考虑到关于本公开的功能而在本领域中当前广泛使用的那些通用术语,但是术语可以根据本领域普通技术人员的意图、先例或本领域中的新技术而变化。此外,一些术语可由申请人任意选择,并且在这种情况下,将在本说明书的详细描述中详细描述所选术语的含义。因此,这里使用的术语必须基于术语的含义以及整个说明书中的描述来定义。
当部件“包括”或“包含”元件时,除非存在与其相反的描述,否则该部件还可包括其它元件,而不排除其他元件。此外,示例性实施例中的术语“单元”表示诸如现场可编程门阵列(FPGA)或专用集成电路(ASIC)的软件组件或硬件组件,并且执行特定功能。然而,术语“单元”不限于软件或硬件。“单元”可以形成在可寻址存储介质中,或者可被形成为操作一个或更多个处理器。因此,例如,术语“单元”可表示诸如软件组件、面向对象的软件组件、类组件和任务组件的组件,并且可包括处理、功能、属性、过程、子例程、程序代码段、驱动器、固件、微代码、电路、数据、数据库、数据结构、表、数组或变量。由组件和“单元”提供的功能可以与更小数量的组件和“单元”相关联,或者可被划分为附加组件和“单元”。
现在将详细参照实施例,其中,实施例的示例在附图中示出。在以下描述中,可以不详细描述公知的功能或构造,以免以不必要的细节模糊示例性实施例。
在整个说明书中,“图像”可以表示由离散图像元素(例如,二维(2D)图像中的像素和三维(3D)图像中的体素)组成的多维数据。例如,图像可以是由X射线设备、计算机断层扫描(CT)设备、磁共振成像(MRI)设备、超声诊断设备或其它医学成像设备捕捉的对象的医学图像。
此外,在本说明书中,“对象”可以是人、动物或者人或动物的一部分。例如,对象可以是器官(例如,肝脏、心脏、子宫、大脑、乳房或腹部)、血管或它们的组合。所述对象可以是体模。体模表示具有与器官大致相同的密度、有效原子序数和体积的材料。例如,体模可以是具有与人体类似的特性的球形体模。
此外,在本说明书中,“用户”可以是但不限于例如医生的医疗专家、护士、医学实验室技术人员或修理医疗设备的技术人员。
此外,在本说明书中,“MR图像”表示通过使用核磁共振原理获取的对象的图像。
此外,在本说明书中,“脉冲序列”表示由MRI设备重复施加的信号的连续性。脉冲序列可包括射频(RF)脉冲的时间参数,例如重复时间(TR)或回波时间(TE)。
此外,这里使用的术语“脉冲序列图”可表示用于解释在MRI系统中施加的信号的序列的示图。例如,脉冲序列示意图可以是根据时间示出射频(RF)脉冲、梯度磁场、MR信号等的示图。
MRI系统是用于通过以对比度比较来表述MR信号的强度相对于在具有强度的磁场中产生的RF信号以获取对象的一部分的横截面图像的设备。例如,如果仅使原子核(例如,氢原子核)共振的RF信号瞬间朝向放置在强磁场中的对象发射,然后这种发射停止,则从原子核发射出MR信号,并因此MRI系统可接收MR信号并获取MR图像。MR信号表示从对象发射出的RF信号。可根据对象的预定原子(例如,氢)的密度、弛豫时间T1、弛豫时间T2和血液的流动等来确定MR信号的强度。
MRI系统包括与其它成像设备的特征不同的特征。与诸如根据检测硬件的方向获取图像的CT设备的成像设备不同,MRI系统可以获取朝向可选点的2D图像或3D体图像。与CT设备、X射线设备、位置发射断层扫描(PET)设备和单光子发射CT(SPECT)设备不同,MRI系统不使对象或检查者暴露于辐射,可以获取具有高软组织对比度的图像,并且可以获取用于精确捕获异常组织的神经图像、血管内图像、肌肉骨骼图像和肿瘤图像。
如这里所使用的,术语“和/或”包括相关联的列出的项中的一个或更多个的任意组合和所有组合。诸如“至少一个”的表述在元素列表之后时修饰整列元素而不修饰列中的个别元素。
图1是示出根据示例性实施例的基于多回波序列由MRI设备执行以获取MR图像的处理的示图。
详细地,参照图1,根据示例性实施例的MRI设备执行基于与第一回波时间TE1相应的第一回波数据E1和与第二回波时间TE2相应的第二回波数据E2重建最终图像的处理。
根据示例性实施例的MRI设备可获取基于多回波序列产生的MR信号。
多回波序列是指用于MRI设备通过使用在施加一次RF脉冲之后产生的多个回波来获取图像的脉冲序列。也就是说,多回波序列在一个TR期间可具有至少两个回波时间。
通过使用多回波序列获取的MR信号可被表示为k空间数据。
k空间数据是指通过将MR信号放置在k空间中而产生的信号,该MR信号是根据包括在高频多线圈中的通道从多个线圈中的每一个线圈接收的RF信号。
k空间数据可以是2D k空间数据或3D k空间数据。例如,2D k空间数据具有2D空间频率域,并且由与频率编码相应的kx轴和与相位编码相应的ky轴形成。此外,3D k空间数据由kx轴、ky轴和与空间中的前进方向相应的kz轴形成。kz轴与切片选择梯度相应。
可通过使用快速傅里叶变换(FFT)或傅里叶变换(FT)将k空间数据重建为MR图像。
当如上所述使用多回波序列时,MRI设备可使用以下方法,该方法通过在一个TR期间仅施加一次用于相位编码的梯度磁场来获取与多个回波时间相应的多条k空间数据。
使用多回波序列的MRI设备可使用通过使用在一个TR期间产生的多个回波来获取一条k空间数据的另一方法。该方法被称为平面回波成像(EPI)方法。
图1是示出根据示例性实施例的MRI设备使用获取分别与多个回波时间相应的多个k空间数据的方法的情况的示图。图1中所示的方块表示k空间。
参照图1,第一回波数据E1可以是通过使用在包括在一个TR中的多个回波时间中的第一回波时间TE1产生的回波而获取的数据。如图1所示,第一回波数据E1可包括k空间的中心部分的数据。此外,第一回波数据E1可以包括k空间中的ky坐标为0的数据。
第二回波数据E2可以是通过使用在包括在一个TR中的多个回波时间中的第二回波时间TE2产生的回波而获取的数据,第二回波时间TE2比第一回波时间TE1晚。如图1所示,第二回波数据E2可不包括k空间的中心部分的数据。此外,第二回波数据E2可不包括k空间中的ky坐标为0的数据。
第一回波数据E1和第二回波数据E2可以指原始数据或k空间数据。
原始数据可以是MR信号,该MR信号是通过MRI处理从包括在高频多线圈中的多个通道线圈中的每一个通道线圈接收的RF信号。
k空间数据可以是在k空间中通过对获取的MR信号进行采样而获取的数据。k空间数据可以是通过在k空间的所有点对获取的MR信号执行完全采样而获取的数据,或者可以是通过针对获取的MR信号执行在一些点执行采样而在其它点不执行采样的欠采样而获取的数据。可以通过使用广义自动校准部分并行采集(GRAPPA)方法或空间谐波同时采集(SMASH)方法来重建通过使用欠采样获取的多条不完整k空间数据中的非获取点的信号,其中,空间谐波同时采集(SMASH)方法使用具有附加线圈信息的图(例如,线圈灵敏度图)。
为方便起见,假设第一回波数据E1和第二回波数据E2是在k空间中采样的k空间数据来解释图1。k空间中的第一回波数据E1和第二回波数据E2的布置不限于图1中的布置,并且可以以各种方式来修改。
返回参照图1的获取MR图像的处理,根据示例性实施例的MRI设备可获取与第一回波时间TE1相应的第一回波数据E1,并且可获取与比第一回波时间TE1晚的第二回波时间TE2相应的第二回波数据E2。
在这种情况下,第一回波数据E1和第二回波数据E2可包括k空间中的重叠部分。例如,第一回波数据E1还可包括关于k空间中的重叠部分的附加数据15。附加数据15可占例如整个k空间的大约10%。详细地,当k空间包括256个行以重建具有256*256的分辨率的图像时,可以通过使用大约25个行或26个行来获取附加数据15。
如图1所示,第一回波数据E1可包括k空间中的k空间的中心部分的数据,并且可以包括附加数据15。在这种情况下,第一回波数据E1可以是“参考回波数据”。此外,为了MRI设备获取参考回波数据而产生的回波的回波时间可被称为“参考回波时间”。
可基于k空间的参考回波数据来确定将被获取的最终图像的特征。换句话说,参考回波时间可根据将被获取的最终图像的特征而变化。
例如,当T1加权图像将被获取时,根据示例性实施例的MRI设备可将第一回波时间TE1确定为参考回波时间。此外,当T2加权图像将被获取时,根据示例性实施例的MRI设备可将比第一回波时间TE1晚的第二回波时间TE2确定为参考回波时间。也就是说,根据将被获取的最终图像的特征,根据示例性实施例的MRI设备可将参考回波时间(即k空间的中心部分的数据的回波时间)确定为早期回波或晚期回波的时间。
根据示例性实施例的MRI设备可针对第二回波数据E2执行相位校正。相位校正是指将回波数据的相位校正为与参考回波时间相应的相位的处理。
当使用在不同回波时间获取的多条数据时,根据示例性实施例的MRI设备可受到回波时移、磁场不均匀性、自旋去相位和涡流的影响。
详细地,当根据示例性实施例的MRI设备使用在不同回波时间获取的多条数据时,多个回波时间可与预先设计的回波时间不同。此外,回波时间偏移的量可彼此不同。
根据示例性实施例,当MRI设备使用梯度回波脉冲序列时,与使用自旋回波脉冲序列的MRI设备相比,MRI设备可受到磁场不均匀性的更严重影响。
此外,当MRI设备使用梯度回波脉冲序列时,因为与使用自旋回波脉冲序列的MRI设备不同,MRI设备不使用180°脉冲,所以与使用自旋回波脉冲序列的MRI设备相比,MRI设备可受到自旋相移的更严重影响。
根据示例性实施例的MRI设备可执行相位校正以减小回波时移、磁场不均匀性和自旋去相位对信号强度的影响中的任何一个或任何组合。
也就是说,根据示例性实施例的MRI设备可通过执行相位校正来减少当使用在不同的回波时间获取的多条数据时可能出现的最终图像中的模糊或混叠的影响。
参照图1,当第一回波数据E1是参考回波数据时,根据示例性实施例的MRI设备可针对第二回波数据E2执行相位校正,并且可以获取基于第一回波时间TE1校正了相位的校正后的第二回波数据E2C。
根据示例性实施例的MRI设备可基于第一回波数据E1和校正后的第二回波数据E2C获取最终的k空间数据101。此外,MRI设备可基于最终的k空间数据101重建最终图像。
最终的k空间数据101可以是在k空间中重建的数据。例如,当第一回波数据E1和第二回波数据E2是欠采样的k空间数据时,最终的k空间数据101可以是通过使用GRAPPA方法或SMASH方法重建的数据。
由根据示例性实施例的MRI设备重建的最终k空间数据101的重叠数据19可以是通过将GRAPPA方法等应用于第一回波数据E1的附加数据15而重建的k空间数据。
根据示例性实施例的MRI设备可通过针对最终k空间数据101执行FT来获取最终图像。根据示例性实施例的MRI设备可通过将诸如灵敏度编码(SENSE)或局部灵敏度并行成像(PILS)的方法应用于最终k空间数据101来重建最终图像。
根据示例性实施例的MRI设备可通过获取关于第一回波数据E1和第二回波数据E2之间的在k空间中的重叠部分19的附加数据15来改善重建的最终图像的质量。
此外,当根据示例性实施例的MRI设备使用在不同回波时间获取的第一回波数据E1和第二回波数据E2时,MRI设备可通过组合多条回波数据的特征来增加最终图像的信噪比(SNR)和对比度。
此外,因为根据示例性实施例的MRI设备通过使用在不同的回波时间获取的第一回波数据E1和第二回波数据E2来获取最终的k空间数据101,所以MRI设备可减少获取最终图像耗费的时间。
图2A是示出根据示例性实施例的基于多回波序列由MRI设备执行以获取k空间数据的处理的示图。
参照图2A,根据示例性实施例的MRI设备使用的多回波序列可以是例如梯度回波序列。
根据示例性实施例的MRI设备可以激发一个RF脉冲,然后可以通过按照使得读出梯度磁场Gread具有相同的强度和交替的符号的方式应用读出梯度磁场Gread来产生多个梯度回波。
为了便于说明根据读出梯度磁场Gread产生的回波,用于相位编码的梯度磁场未在图2A中示出。
参照图2A,根据示例性实施例的MRI设备可在一个TR期间产生多个回波。也就是说,根据示例性实施例的MRI设备可以在一个TR期间在第一回波时间TE1、第二回波时间TE2、第三回波时间TE3和第四回波时间TE4产生回波。
参照图2A,根据示例性实施例的MRI设备可通过使用与第一回波时间TE1相应的读出梯度磁场,沿正方向211用数据填充k空间201的一个行。
此外,根据示例性实施例的MRI设备可以通过使用与第二回波时间TE2相应的读出梯度磁场,沿着负方向213用数据填充k空间203的一个行。
此外,根据示例性实施例的MRI设备可通过使用与第三回波时间TE3相应的读出梯度磁场,沿正方向215用数据填充k空间205的一个行。
此外,根据示例性实施例的MRI设备可通过使用与第四回波时间TE4相应的读出梯度磁场沿负方向217用数据填充k空间207的一个行。
根据示例性实施例的MRI设备可通过使用多个回波时间来获取具有不同特征的多条k空间数据。
例如,当通过使用在作为短回波时间的第一回波时间TE1获取的k空间数据来重建图像时,MRI设备可获取T2的影响或T2*的影响减小并且T1的影响增加的图像。
此外,当通过使用在作为长回波时间的第四回波时间TE4获取的k空间数据来重建图像时,MRI设备可获取T1的影响减小并且T2的影响或T2*的影响增加的图像。
根据示例性实施例的MRI设备可通过组合具有不同回波时间的多个k空间数据重建最终图像来调整组织对比度、包括在横截面图像中的静脉血液的清晰度、SNR等。
图2B是示出根据示例性实施例的基于梯度回波序列由MRI设备执行以获取MR图像的处理的曲线图。
根据示例性实施例的MRI设备可获取在回波时间具有最高强度的MR信号。可通过例如包括在MRI设备中的RF线圈来接收该MR信号。
根据示例性实施例的MRI设备可将与在回波时间获取的具有最高强度的MR信号相应的k空间数据放置在k空间中的kx坐标为0的位置。
例如,如参考图2A所述,当MRI设备通过使用与第一回波时间TE1相应的读出梯度磁场沿正方向211用数据填充空间201的一个行时,放置在k空间中的kx坐标为0的位置处的k空间数据可与回波信号的最高点相应。
参照图2B,根据示例性实施例的MRI设备可使用梯度回波序列来获取与在第一回波时间TE1、第二回波时间TE2、第三回波时间TE3和第四个回波时间TE4的回波信号的最高点相应的MR信号。
在这种情况下,如图2B所示,在一个TR期间,在根据示例性实施例的MRI设备获取的多个回波中可能发生回波时移。也就是说,可以在作为在一个TR期间获取的多个回波的多个回波时间的第一回波时间TE1、第二回波时间TE2、第三回波时间TE3和第四回波时间TE4获取与除回波信号的最高点以外的点相应的MR信号。
例如,如图2B所示,第一回波时间TE1可以向左偏移T1,第二回波时间TE2可以向右偏移T2。此外,第三回波时间TE3可以向左偏移T3,并且第四回波时间TE4可以向右偏移T4。
根据示例性实施例的MRI设备可以通过补偿回波时间偏移的量来以一维(1D)方式校正获取的k空间数据。此外,根据示例性实施例的MRI设备可以以2D或3D方式针对k空间中的多条数据执行相位校正。
例如,MRI设备可以使用B0图来执行相位校正。通过根据k空间中的回波时间差布置相位或频率差来获取B0图。详细地,MRI设备可通过使用具有回波时间差的两条k空间数据之间的相位差来获取k空间中的B0图。MRI设备可获取与映射到B0图的相位差相应的频率差。MRI设备可以产生B0图,其中,与根据回波时间差的频率相应的颜色被映射到该B0图。通过使用经由使用具有回波时间差的两条k空间数据之间的相位差而获取的B0图来执行相位校正的这种处理可被称为共轭相位重建。为了通过使用B0图执行相位校正,可以使用诸如共轭相位重建的各种公知方法中的任何一种。
根据示例性实施例的MRI设备可减小k空间中的数据的回波时移的影响,并且可以通过使用B0图执行相位校正来减小获取的最终图像中的模糊或混叠的影响。
图3A是根据示例性实施例的MRI设备300a的框图。
图3A的MRI设备300a是用于通过使用多回波序列获取MR图像的设备。详细地,MRI设备300a可以是用于通过使用多回波序列捕捉对象的MR图像的设备或用于处理通过使用多回波序列捕捉对象的MR图像而获取的数据的设备。
例如,MRI设备300a可以是用于通过包括在高频多线圈中的多个通道线圈将使用多回波序列的RF脉冲施加到对象并且通过使用通过所述多个通道线圈获取的MR信号重建MR图像的设备。
此外,MRI设备300a可以是用于提供将被施加于对象的多回波序列并通过使用获取的MR信号重建MR图像的服务器设备。服务器设备可以是执行MRI处理的医院等中的医疗服务器设备。
参照图3A,MRI设备300a可包括数据获取器310和图像处理器320。
数据获取器310可通过使用在第一回波时间TE1产生的回波来获取第一回波数据E1,并且可通过使用在比第一回波时间TE1晚的第二回波时间TE2产生的回波来获取第二回波数据E2。第一回波数据E1和第二回波数据E2可包括k空间中的重叠部分。
根据示例性实施例,第一回波数据和第二回波数据可以是原始数据,或者可以是与原始数据相应的k空间数据。
数据获取器310可获取分别与多个回波相应的多条不完整k空间数据。术语“不完整k空间数据”表示当MR信号在k空间中的至少一个点未被采样并且因此将在k空间中的所述至少一个点被重建时的k空间数据。可以通过使用获取的信号来重建来自“不完整k空间数据”的未获取的信号。
根据示例性实施例的数据获取器310可通过使用在参考回波数据的参考回波时间所产生的回波来获取关于重叠部分的附加数据。
参考回波数据可表示包括k空间的中心部分的数据的回波数据。此外,参考回波时间可表示为了获取参考回波数据而产生的回波的回波时间。
根据示例性实施例,第一回波数据和第二回波数据中的一个或每一个可包括k空间的中心部分的数据。
根据示例性实施例的数据获取器310可以确定将被获取的MR图像的特征。此外,数据获取器310可根据确定的MR图像的特征来确定第一回波数据或第二回波数据将包括k空间的中心部分的数据。确定MR图像的特征的步骤可包括确定例如MR图像是T1加权图像还是T2(或T2*)加权图像。
当将获取T1加权图像时,数据获取器310可以将最短回波时间确定为参考回波时间。在这种情况下,可通过使用在最短回波时间产生的回波来获取参考回波数据。例如,当将获取T1加权图像时,数据获取器310可确定第一回波数据将包括k空间的中心部分的数据。在这种情况下,参考回波时间可以被称为与早期回波相应。
此外,当将获取T2或T2*加权图像时,数据获取器310可将最长回波时间确定为参考回波时间。在这种情况下,可通过使用在最长回波时间产生的回波来获取参考回波数据。例如,当要获取T2或T2*加权图像时,数据获取器310可确定第二回波数据将包括k空间的中心部分的数据。在这种情况下,参考回波时间可以被称为与后期回波相应。
当第一回波数据包括k空间的中心部分的数据并且第二回波数据不包括k空间的中心部分的数据时,数据获取器310可通过使用在第一个回波时间产生的回波获取关于重叠部分的附加数据。
当第二回波数据包括k空间的中心部分的数据并且第一回波数据不包括k空间的中心部分的数据时,数据获取器310可通过使用在第二个回波时间产生的回波获取关于重叠部分的附加数据。
根据示例性实施例,已经解释了第一回波数据和第二回波数据是通过对在k空间中的从对象发出的磁MR信号进行采样而获取的k空间数据的情况。然而,数据获取器310可仅从RF线圈接收MR信号,并且图像处理器320可在k空间中执行采样处理。
图像处理器320可基于第一回波数据和第二回波数据来重建MR图像。
例如,第一回波数据和第二回波数据可以是k空间数据,即通过对MR信号进行欠采样而获取的不完整k空间数据。在这种情况下,可基于具有附加线圈信息(例如,线圈灵敏度图)的图(诸如,SMASH方法)或附加校准信号(诸如,GRAPPA方法)来重建不完整k空间数据。
此外,图像处理器320可针对第一回波数据和第二回波数据中的一个或两者执行相位校正,并且可以通过使用已经被执行了相位校正的第一回波数据和第二回波数据中的一个或两者来重建MR图像。
图像处理器320可基于第一回波数据和第二回波数据中的包括k空间的中心部分的数据的数据的相位来执行相位校正。
例如,当第一回波数据是参考回波数据时,图像处理器320可以基于第一回波数据的相位对第二回波数据执行相位校正。
图3B是根据另一示例性实施例的MRI设备300b的框图。
MRI设备300b可包括数据获取器315、图像处理器325、梯度磁场控制器330和输出接口340。
图3B的数据获取器315和图像处理器325可以执行与图3A的数据获取器310和图像处理器320的功能相同的功能,因此,将不给出数据获取器315和图像处理器325的重复解释。
数据获取器315可通过使用在第一回波时间TE1产生的回波获取第一回波数据E1,并且可通过使用在比第一回波时间TE1晚的第二回波时间TE2产生的回波来获取第二回波数据E2。第一回波数据E1和第二回波数据E2可包括k空间中的重叠部分。
根据示例性实施例的图像处理器325可通过使用第三回波数据和第四回波数据来产生B0图,其中,第三回波数据和第四回波数据是通过使用在第三回波时间和第四回波时间产生的回波来获取的。
根据示例性实施例,第一回波时间、第二回波时间、第三回波时间和第四回波时间可被包括在多回波序列的一个TR中。
此外,图像处理器325可通过使用产生的B0图对第一回波数据和第二回波数据中的一个或两者执行相位校正。
根据示例性实施例的梯度磁场控制器330可根据多回波序列控制梯度磁场被施加。
由梯度磁场控制器330在第一回波时间施加的第一读出梯度磁场和由梯度磁场控制器330在第二回波时间施加的第二读出梯度磁场可具有相同的符号。
此外,由梯度磁场控制器330在第三回波时间施加的第三读出梯度磁场和由梯度磁场控制器330在第四回波时间施加的第四读出梯度磁场可以具有相同的符号。
输出接口340可以将由图像处理器325产生的图像数据和重建的MR图像输出给用户。此外,输出接口340可以输出将被用于用户操作MRI设备300b的信息,诸如用户界面(UI)、用户信息或对象信息。
详细地,输出接口340可以在各种显示器(诸如,阴极射线管(CRT)显示器、液晶显示器(LCD)、等离子体显示板(PDP)显示器、有机发光二极管(OLED)显示器、场发射显示器(FED)、发光二极管(LED)显示器、真空荧光显示器(VFD)、数字光处理(DLP)显示器、平板显示器(FPD)、3D显示器或透明显示器)中的任何一个上显示预定的屏幕。
此外,输出接口340可显示用于显示关于MRI设备300b的MRI处理的信息的UI屏幕。
例如,输出接口340可以显示用于设计多回波序列的UI屏幕和用于显示MRI设备300b的成像环境和在MRI设备300b中使用的多回波序列的屏幕。
此外,输出接口340可以显示在重建MR图像的处理中产生的图像、关于当前正在执行的MRI协议的信息以及关于随后将被执行的MRI协议的信息。
此外,输出接口340可显示包括颜色被映射到的B0图的屏幕,以便用户可视地识别被产生用于执行相位校正的B0图。
MRI设备300b还可包括输入接口,并且可基于通过UI屏幕输入的信息确定将被用于重建MR图像的信息和将在MRI处理中使用的多回波序列。
详细地,MRI设备300b可基于通过UI屏幕输入的关于回波链长度(ETL)的信息来确定将在MRI处理中使用的多回波序列。此外,MRI设备300b可基于通过UI屏幕输入的关于参考回波数据的信息来确定参考回波数据。
此外,MRI设备300b可以基于通过UI屏幕输入的关于重建MR数据的方法的信息来确定重建方法。
此外,MR设备300b可基于通过UI屏幕输入的信息来确定被用于产生B0图的数据的类型。
图4A是根据示例性实施例的由MRI设备300a或300b施加的ETL为2的脉冲序列的脉冲序列图410。
参照如图4A,MRI设备300a或300b可根据脉冲序列图410中所示的多回波序列施加梯度磁场420。
详细地,脉冲序列图410可示出在一个TR期间由根据示例性实施例的MRI设备300a或300b施加的相位编码方向的梯度磁场PE1、梯度磁场PE2和梯度磁场PE2的翻转以及读出方向的梯度磁场420。读出方向可与k空间的kx方向相应,并且相位编码方向可与k空间的ky方向相应。
脉冲序列图410可示出多回波序列。用于在多回波序列的一个TR中获取数据的回波的数量可被称为ETL。根据图4A的脉冲序列图410的多回波序列,MRI设备300a和300b通过在第一回波时间TE1和第二回波时间TE2分别获取第一回波数据E1和第二回波数据E2来重建最终的k空间。也就是说,在脉冲序列图410中,ETL为2。
图4A的相位编码方向的梯度磁场PE1、梯度磁场PE2和梯度磁场PE2的翻转表示应用相位编码方向的梯度磁场的时间间隔。相位编码方向的梯度磁场PE1、梯度磁场PE2和梯度磁场PE2的翻转根据TR具有不同的强度。
参照图4A,如脉冲序列图410所示,MRI设备300a或300b可施加用于沿k空间的读出方向产生多个回波的梯度磁场420。MRI设备300a或300b可通过施加梯度磁场420在一个TR期间在第一回波时间TE1、第二回波时间TE2和第三回波时间TE1'产生回波。
在获取第一回波数据E1之前,MRI设备300a或300b可通过应用相位编码方向的梯度磁场PE,确定第一回波数据E1在k空间中的ky轴上的位置。
在应用相位编码方向的梯度磁场PE1之后,MRI设备300a或300b可通过使用与第一个回波时间TE1相应的读出梯度磁场沿正方向411填充k空间的第一回波数据E1中的一行数据。第一回波数据E1可以是用于重建最终的k空间数据的数据,并且可以是用于执行相位校正的数据。
MRI设备300a或300b可通过反转读出梯度磁场的符号通过使用与第三回波时间TE1'相应的读出梯度磁场通沿负方向413填充k空间的第三回波数据E1'中的一行数据。根据图4A的示例性实施例,第三回波数据E1'可以是用于执行相位校正的数据,而不是用于重建最终的k空间数据的数据。根据另一示例性实施例,第三回波数据E1'可包括用于重建最终的k空间数据的数据。
在获取第二回波数据E2之前,MRI设备300a或300b可通过应用相位编码方向的梯度磁场PE2来确定第二回波数据E2在k空间中的ky轴上的位置。
参照图4A,当MRI设备300a或300b施加相位编码方向的梯度磁场PE2时,k空间中的数据的位置可远离ky轴的原点。
MRI设备300a或300b可通过重新反转读出梯度磁场的符号通过使用与第二回波时间TE2相应的读出梯度磁场沿正方向415填充k空间的第二回波数据E2中的的一行数据。接下来,在进入下一个TR之前,MRI装置300a或300b可施加相位编码方向的梯度磁场PE2的翻转,其中,梯度磁场PE2的翻转具有与相位编码方向的梯度磁场PE2相同的强度和相反的方向。当MRI设备300a或300b施加相位编码方向的梯度磁场PE2的翻转时,将在所述下一个TR中获取的k空间的数据可被放置在ky轴的原点。
接下来,在所述下一个TR中,MRI设备300a或300b可通过应用相位编码方向的梯度磁场PE1来确定第一回波数据E1在k空间中的ky轴上的位置。
参照图4A所示,MRI设备300a或300b可通过在多个TR期间施加梯度磁场420来获取与第一回波时间TE1相应的第一回波数据E1、与第二回波时间TE2相应的第二回波数据E2以及与第三回波时间TE1'相应的第三回波数据E1'。
当MRI设备300a或300b获取第一回波数据E1和第二回波数据E2时在k空间中放置多条数据的方向411和415是相同的。也就是说,在第一回波时间TE1施加的读出梯度磁场和在第二回波时间TE2施加的读出梯度磁场具有相同的符号。
当读出梯度磁场具有相同的符号时,发生回波时移的方向可以是相同的,并且在这种情况下,MRI设备300a或300b可以减少重建的最终图像中的由回波时移引起的误差。
参照图4A,因为MRI设备300a或300b在获取第一回波数据E1之后不施加相位编码梯度磁场以获取第三回波数据E1',所以第一回波数据E1和第三回波数据E1'在k空间中的ky轴上的位置可以彼此相反。
当MRI设备300a或300b获取第一回波数据E1和第三回波数据E1'时在k空间中放置数据的方向411和413彼此不同。如果MRI设备300a或300b通过使用第一回波数据E1和第三回波数据E1'产生B0图,则MRI设备300a或300b被用于校正当回波时移发生的方向彼此不同时引起的误差。
如图4A所示,第一回波数据E1可包括k空间的中心部分的数据,第二回波数据E2可不包括k空间的中心部分的数据。
参照图4A,当将获取T1加权图像时,MRI设备300a或300b可将第一回波时间TE1确定为参考回波时间。在这种情况下,作为参考回波数据的第一回波数据E1可包括k空间中的与第二回波数据E2重叠的部分。
图4B是示出根据示例性实施例的根据图4的脉冲序列图410由MRI设备300a或300b执行以获取MR图像的处理的示图。
参照图4B,MRI设备300a或300b执行基于与第一回波时间TE1相应的第一回波数据E1和与第二回波时间TE2相应的第二回波数据E2重建最终图像的处理。
参照图4B,根据示例性实施例的MRI设备300a或300b可获取与第一回波时间TE1相应的第一回波数据E1,然后可获取与比第一回波时间TE1晚的第二回波时间TE2相应的第二回波数据E2。此外,MRI设备300a或300b可获取与第三回波时间TE1'相应的第三回波数据E1'。
参照图4B,参考回波数据可以是第一回波数据E1,参考回波时间可以是第一回波时间TE1。在这种情况下,第一回波数据E1可包括k空间中的与第二回波数据E2重叠的部分。例如,第一回波数据E1还可包括关于k空间中的与第二回波数据E2重叠的部分的附加数据425。附加数据425可以占例如整个k空间的大约10%。
根据示例性实施例的MRI设备300a或300b可针对第二回波数据E2执行相位校正。MRI设备300a或300b可通过使用第一回波数据E1和第三回波数据E1'来产生B0图,并且可以通过使用产生的B0图来执行相位校正。
如参照图4A所述,当在MRI设备300a或300b获取第一回波数据E1和第三回波数据E1'时在k空间中布置多条数据的方向411和413彼此不同时,MRI设备300a或300b被用于校正当回波时移发生的方向彼此不同时引起的误差以产生B0图。
MRI设备300a或300b可针对第二回波数据E2执行相位校正,并且可获取基于第一回波时间TE1校正了相位的第二回波数据E2C。由于相位校正,MRI设备300a或300b可减小当使用在不同回波时间获取的多条数据时可能发生的最终图像中的模糊或混叠的影响。
MRI设备300a或300b可以基于第一回波数据E1和校正后的第二回波数据E2C获取最终的k空间数据E1+E2C。如果第一回波数据E1和第二回波数据E2是欠采样的k空间数据,则最终的k空间数据E1+E2C可以是通过使用GRAPPA方法或SMASH方法重建的数据。
当重建最终的k空间数据E1+E2C时,MRI设备300a或300b可基于除了附加数据425之外的第一回波数据E1的数据和第二回波数据E2来重建k空间数据。MRI设备300a或300b可使用GRAPPA方法等来重建k空间数据。接下来,MRI设备300a或300b可基于包括附加数据425的第一回波数据E1和校正后的第二回波数据E2C重新重建k空间数据。在这种情况下,可以将GRAPPA方法等重新应用于附加数据425,并且MRI设备300a或300b可以获取最终的k空间数据E1+E2C。
根据图4A和图4B的示例性实施例,MRI设备300a或300b可通过使用在一个TR期间获取的回波数据来产生用于执行相位校正的B0图和最终的k空间数据。因此,MRI设备300a或300b可通过仅使用用于获取k空间数据的脉冲序列来获取将被用于相位校正的数据,而无需针对相位校正另外使用脉冲序列。
图5A是根据示例性实施例的由MRI设备300a或300b施加的ETL为2的另一脉冲序列的脉冲序列图510。
参照图5A,MRI设备300a或300b可根据脉冲序列图510中所示的多回波序列来施加梯度磁场520。
与图4A一样,脉冲序列图510可示出在一个TR期间由根据示例性实施例的MRI设备300a或300b施加的相位编码方向的梯度磁场PE1、梯度磁场PE2和梯度磁场PE2的翻转以及读出方向的梯度磁场520。
参照图5A,如脉冲序列图510所示,MRI设备300a或300b可施加用于沿k空间的读出方向产生多个回波的梯度磁场520。
MRI设备300a或300b可通过施加梯度磁场520在一个TR期间在第一回波时间TE1、第二回波时间TE2、第三回波时间TE1'和第四回波时间TE2'产生回波。也就是说,图5A的脉冲序列图510与图4A的脉冲序列图410的不同之处在于可另外产生与第四回波时间TE2'相应的回波。
MRI设备300a或300b可通过使用在一个TR中包括的第一回波时间TE1和第二回波时间TE2所产生的回波来获取第一回波数据E1和第二回波数据E2。MRI设备300a或300b通过获取第一回波数据E1和第二回波数据E2来重建最终的k空间。
根据图5A的示例性实施例,第三回波数据E1'和第四回波数据E2'可以是用于相位校正的数据,而不是用于重建最终的k空间数据的数据。根据另一示例性实施例,第三回波数据E1'和第四回波数据E2'可包括用于重建最终的k空间数据的数据。
首先,MRI设备300a或300b可通过在获取第一回波数据E1之前应用相位编码方向的梯度磁场PE1来确定第一回波数据E1在k空间中的ky轴上的位置。
在应用相位编码方向的梯度磁场PE1之后,MRI设备300a或300b可以通过使用与第一回波时间TE1相应的读出梯度磁场沿正方向511填充k空间的第一回波数据E1中的一行数据。
MRI设备300a或300b可通过反转读出梯度磁场的符号通过使用与第三回波时间TE1'相应的读出梯度磁场来沿负方向513填充k空间的第三回波数据E1'中的一行数据。
MRI设备300a或300b可以通过在获取第二回波数据E2之前应用相位编码方向的梯度磁场PE2来确定第二回波数据E2在k空间中的ky轴上的位置。
MRI设备300a或300b可通过重新反转读出梯度磁场的符号通过使用与第二回波时间TE2相应的读出梯度磁场来沿正方向515填充k空间的第二回波数据E2中的一行数据。
接下来,MRI设备300a或300b可施加相位编码方向的梯度磁场PE2的翻转,其中,梯度磁场PE2的翻转具有与相位编码方向的梯度磁场PE2相同的强度和相反的方向。接下来,MRI设备300a或300b可通过重新反转读出梯度磁场的符号通过使用与第四回波时间TE2'相应的读出梯度磁场来沿负方向517填充k空间的第四回波数据E2'中的一行数据。
当MRI设备300a或300b应用相位编码方向的梯度磁场PE2的翻转时,下一个TR的k空间中的数据可被放置在ky轴的原点。
接下来,在所述下一个TR中,MRI设备300a或300b可通过应用相位编码方向的梯度磁场PE1来确定k空间中的ky轴上的回波数据E1的位置。
参照图5A,MRI设备300a或300b可通过在多个TR期间施加梯度磁场520来获取与第一回波时间TE1相应的第一回波数据E1、与第二回波时间TE2相应的第二回波数据E2、与第三回波时间TE1'相应的第三回波数据E1'以及与第四回波时间TE2'相应的第四回波数据E2'。
参照图5A,第三回波数据E1'和第四回波数据E2'在k空间中的ky轴上的位置可以彼此相对。第一回波数据E1和第三回波数据E1'在ky轴上的位置可以彼此不同。
当MRI设备300a或300b获取第三回波数据E1'和第四回波数据E2'时在k空间中放置多条数据的方向513和517是相同的。如果MRI设备300a或300b通过使用第三回波数据E1'和第四回波数据E2'产生B0图,则MRI设备300a或300b可以最小化当回波时移发生的方向彼此不同时引起的误差。
如图5A所示,第一回波数据E1可包括k空间的中心部分的数据,第二回波数据E2可不包括k空间的中心部分的数据。
参照图5A,当将获取T1加权图像时,MRI设备300a或300b可将第一回波时间TE1确定为参考回波时间。在这种情况下,作为参考回波数据的第一回波数据E1可包括k空间中的与第二回波数据E2重叠的部分。
图5B是示出根据示例性实施例的根据图5A的脉冲序列图510由MRI设备300a或300b执行以获取MR图像的处理的示图。
参照图5B,MRI设备300a或300b执行基于与第一回波时间TE1相应的第一回波数据E1和与第二回波时间TE2相应的第二回波数据E2重建最终图像的处理。
根据示例性实施例的MRI设备300a或300b可获取与第一回波时间TE1相应的第一回波数据E1,并且可以获取与比第一回波时间晚的第二回波时间TE2相应的第二回波数据E2。此外,MRI设备300a或300b可获取与第三回波时间TE1'相应的第三回波数据E1'和与第四回波时间TE2'相应的第四回波数据E2'以执行相位校正。
参照图5B,参考回波数据可以是第一回波数据E1,参考回波时间可以是第一回波时间TE1。在这种情况下,第一回波数据E1可包括k空间中的与第二回波数据E2重叠的部分。例如,第一回波数据E1还可包括关于k空间中的与第二回波数据E2重叠的部分的附加数据525。附加数据525可占例如整个k空间的大约10%。
根据示例性实施例的MRI设备300a或300b可针对第二回波数据E2执行相位校正。如参考图5A所述,因为当MRI设备300a或300b获取第三回波数据E1'和第四回波数据E2'时在k空间中放置多条数据的方向513和517是相同的,所以MRI设备300a或300b可最小化由回波时移引起的误差。
MRI设备300a或300b可针对第二回波数据E2执行相位校正,并且可以获取基于第一回波时间TE1校正了相位的第二回波数据E2C。可以通过使用B0图来执行相位校正,其中,B0图是通过使用第三回波数据E1'和第四回波数据E2'产生的。
MRI设备300a或300b可基于第一回波数据E1和校正后的第二回波数据E2C获取最终的k空间数据E1+E2C。最终的k空间数据E1+E2C可以是在k空间中重建的k空间数据。当重建最终的k空间数据E1+E2C时,可以基于除了附加数据525之外的第一回波数据E1的数据和第二回波数据E2来首先重建k空间数据。MRI设备300a或300b可以使用GRAPPA方法等来重建k空间数据。接下来,MRI设备300a或300b可以基于包括附加数据525的第一回波数据E1和校正后的第二回波数据E2C重新重建k空间数据。在这种情况下,可以将GRAPPA方法等重新应用于附加数据525,并且MRI设备300a或300b可以获取最终的k空间数据E1+E2C。
根据图5A和图5B的示例性实施例,MRI设备300a或300b可通过使用通过施加具有相同符号的梯度磁场而获取的第一回波数据E1和第二回波数据E2来获取k空间数据。同样地,MRI设备300a或300b可通过使用通过施加具有相同符号的梯度磁场而获取的第三回波数据E1'和第四回波数据E2'来产生B0图。因此,MRI设备300a或300b可以最小化回波时间的误差。
此外,根据图5A和图5B的示例性实施例,MRI设备300a或300b可以通过使用在一个TR期间获取的回波数据来产生用于执行相位校正的B0图和最终的k空间数据。因此,MRI设备300a或300b可通过仅使用用于获取k空间数据的脉冲序列来获取将被用于相位校正的数据,而无需针对相位校正另外使用脉冲序列。
图6A是根据示例性实施例的由MRI设备300a或300b施加的ETL为3的脉冲序列的脉冲序列图610。
参照图6A,MRI设备300a或300b可根据脉冲序列图610中所示的多回波序列来施加梯度磁场620。
脉冲序列图610可以示出在一个TR期间由根据示例性实施例的MRI设备300a或300b施加的相位编码方向的梯度磁场PE1、梯度磁场PE2、梯度磁场PE2的翻转、梯度磁场PE2+PE3和梯度磁场PE的翻转以及读出方向的梯度磁场620。
参照图6A,MRI设备300a或300b可施加如脉冲序列图610所示的用于沿k空间的读出方向产生多个回波的梯度磁场620。
MRI设备300a或300b可通过施加梯度磁场在一个TR期间在第一回波时间TE1、第二回波时间TE2、第三回波时间TE1'、第四回波时间TE2'和第五回波时间TE3产生回波。也就是说,图6A的脉冲序列图610与图5A的脉冲序列图510的不同之处在于可另外产生与第五回波时间TE3相应的回波。
MRI设备300a或300b可以通过使用在一个TR中包括的第一回波时间TE1、第二回波时间TE2和第五回波时间TE3产生的回波来获取第一回波数据E1、第二回波数据E2和第五回波数据E3。MRI设备300a或300b通过获取第一回波数据E1、第二回波数据E2和第五回波数据E3来重建最终的k空间。
根据图6A的示例性实施例,第三回波数据E1'和第四回波数据E2'可以是用于相位校正的数据,而不是用于重建最终的k空间数据的数据。根据另一示例性实施例,第三回波数据E1'和第四回波数据E2'可包括用于重建最终的k空间数据的数据。
首先,MRI设备300a或300b可通过在获取第一回波数据E1之前施加相位编码方向的梯度磁场PE1来确定第一回波数据E1在k空间中的ky轴上的位置。
在施加相位编码方向的梯度磁场PE1之后,MRI设备300a或300b可通过使用与第一个回波时间TE1相应的读出梯度磁场沿正方向611填充k空间的第一回波数据E1中的一行数据。
MRI设备300a或300b可通过反转读出梯度磁场的符号通过使用与第三回波时间TE1'相应的读出梯度磁场来沿负方向613填充k空间的第三回波数据E1'中的一行数据。
MRI设备300a或300b可以通过在获取第二回波数据E2之前施加相位编码方向的梯度磁场PE2来确定第二回波数据E2在k空间中的ky轴上的位置。
MRI设备300a或300b可通过重新反转读出梯度磁场的符号通过使用与第二回波时间TE2相应的读出梯度磁场来沿正方向615填充k空间的第二回波数据E2中的一行数据。
接下来,MRI设备300a或300b可以施加相位编码方向的梯度磁场PE2的翻转,其中,梯度磁场PE2的翻转具有与相位编码方向的梯度磁场PE2相同的强度和相反的方向。接下来,MRI设备300a或30b可通过重新反转读出梯度磁场的符号通过使用与第四回波时间TE2'相应的读出梯度磁场来沿负方向617填充k空间的第四回波数据E2'中的一行数据。
MRI设备300a或300b可通过施加相位编码方向的梯度磁场PE2+PE3来确定第五回波数据E3在第k空间中的ky轴上的位置。
MRI设备300a或300b可以通过重新反转读出梯度磁场的符号通过使用与第五回波时间TE3相应的读出梯度磁场来沿着正方向619填充k空间的第五回波数据E3中的一行数据。
参照图6A,MRI设备300a或300b可通过在多个TR期间施加梯度磁场620来获取与第一回波时间TE1相应的第一回波数据E1、与第二回波时间TE2相应的第二回波数据E2、与第三回波时间TE1'相应的第三回波数据E1'、与第四回波时间TE2'相应的第四回波数据E2'以及与第五回波时间TE3相应的第五回波数据E3。
当MRI设备300a或300b获取第一回波数据E1、第二回波数据E2和第五回波数据E3时在k空间中放置多条数据的方向611、615和619是相同的。也就是说,在第一回波时间TE1施加的读出梯度磁场、在第二回波时间TE2施加的读出梯度磁场和在第五回波时间TE3施加的读出梯度磁场具有相同的符号。
参照图6A,第三回波数据E1'和第四回波数据E2'在k空间中的ky轴上的位置可以彼此相对。第一回波数据E1和第三回波数据E1'在ky轴上的位置可以彼此不同。
当MRI设备300a或300b获取第三回波数据E1'和第四回波数据E2'时在k空间中放置多条数据的方向613和617是相同的。如果MRI设备300a或300b通过使用第三回波数据E1'和第四回波数据E2'产生B0图,则可以最小化当回波时移的方向彼此不同时引起的误差。
如图6A所示,第一回波数据E1可包括k空间的中心部分的数据并且第二回波数据E2和第五回波数据E3可不包括k空间的中心部分的数据。
参照图6A,当将获取T1加权图像时,MRI设备300a或300b可将第一回波时间TE1确定为参考回波时间。在这种情况下,作为参考回波数据的第一回波数据E1可包括k空间中的与第二回波数据E2重叠的部分。此外,第二回波数据E2可包括k空间中的与第五回波数据E3重叠的部分。
图6B是示出根据示例性实施例的根据图6的脉冲序列图610由MRI设备300a或300b执行以获取MR图像的处理的示图。
参照图6B,MRI设备300a或300b执行基于与第一回波时间TE1相应的第一回波数据E1、与第二回波时间TE2相应的第二回波数据E2和与第五回波时间TE3相应的第五回波数据E3重建最终图像的处理。
根据示例性实施例的MRI设备300a或300b可获取与第一回波时间TE1相应的第一回波数据E1,可获取比第一回波时间TE1晚的第二回波时间TE2相应的第二回波数据E2,并且可获取与比第二回波时间TE2晚的第五回波时间TE3相应的第五回波数据E3。此外,MRI设备300a或300b可获取与第三回波时间TE1'相应的第三回波数据E1'和与第四回波时间TE2'相应的第四回波数据E2'以执行相位校正。
参照图6B,参考回波数据可以是第一回波数据E1,参考回波时间可以是第一回波时间TE1。在这种情况下,第一回波数据E1还可包括关于k空间中的与第二回波数据E2重叠的部分的附加数据625。此外,第二回波数据E2还可包括关于k空间中的与第三回波数据E3重叠的部分的附加数据635。
首先,根据示例性实施例的MRI设备300a或300b可针对第五回波数据E3执行相位校正。如参照图6A所述,因为当MRI设备300a或300b获取用于相位校正的第三回波数据E1'和第四回波数据E2'时在k空间中放置多条数据的方向613和617是相同的,所以可最小化由回波时移引起的误差。
MRI设备300a或300b可针对第五回波数据E3执行相位校正,并且可以获取基于第二回波时间TE2校正了相位的第五回波数据E3CE2。可通过使用B0图来执行相位校正,其中,B0图是通过使用第三回波数据E1'和第四回波数据E2'产生的。
MRI设备300a或300b针对第五回波数据E3执行第一相位校正,然后针对基于第二回波数据E2和校正的第五回波数据E3CE2获取的中间k空间数据E2+E3CE2执行第二相位校正。可以基于第二回波数据E2的相位来执行第一相位校正,并且可以基于第一回波数据E1的相位来执行第二相位校正。
MRI设备300a或300b可基于第二回波数据E2和校正后的第五回波数据E3CE2获取中间k空间数据E2+E3CE2。如果第二回波数据E2和第五回波数据E3是欠采样的k空间数据,则中间k空间数据E2+E3CE2可以是通过使用GRAPPA方法或SMASH方法重建的数据。
此外,当重建中间k空间数据E2+E3CE2时,MRI设备300a或300b可基于除了附加数据635之外的第二回波数据E2和第五回波数据E3来重建k空间数据。MRI设备300a或300b可使用GRAPPA方法等来重建k空间数据。接下来,MRI设备300a或300b可基于包括附加数据635的第二回波数据E2和校正后的第五回波数据E3CE2重新重建k空间数据。在这种情况下,可以将GRAPPA方法等重新应用于附加数据635,并且MRI设备300a或300b可获取中间k空间数据E2+E3CE2。
接下来,MRI设备300a或300b可针对基于第二回波数据E2和校正后的第五回波数据E3CE2获取的中间k空间数据E2+E3CE2执行第二相位校正。
MRI设备300a或300b可针对中间k空间数据E2+E3CE2执行相位校正,并且可获取基于第一回波时间TE1校正了相位的中间k空间数据(E2+E3CE2)CE1。可以通过使用B0图来执行相位校正,其中,该B0图是通过使用第三回波数据E1'和第四回波数据E2'产生的。
MRI设备300a或300b可基于第一回波数据E1和校正后的中间k空间数据(E2+E3CE2)CE1获取最终的k空间数据E1+(E2+E3CE2)CE1。最终的k空间数据E1+(E2+E3CE2)CE1可以是在k空间中重建的k空间数据。
由MRI设备300a或300b执行的通过使用附加数据625重建最终k空间数据E1+(E2+E3CE2)CE1的方法类似于参照图5B描述的方法,因此,将不给出其详细说明。
根据图6A和图6B的示例性实施例,MRI设备可通过使用通过施加具有相同符号的梯度磁场而获取的第一回波数据E1、第二回波数据E2和第五回波数据E3来获取k空间数据。同样地,MRI设备300a或300b可以通过使用通过施加具有相同符号的梯度磁场而获取的第三回波数据E1'和第四回波数据E2'来产生B0图。因此,MRI设备300a或300b可最小化回波时间的误差。
此外,根据图6A和图6B的示例性实施例,MRI设备300a或300b可通过使用在一个TR期间获取的回波数据来产生用于执行相位校正的B0图和最终的k空间数据。因此,MRI设备300a或300b可通过仅使用用于获取k空间数据的脉冲序列来获取将被用于相位校正的数据,而无需针对相位校正另外使用脉冲序列。
根据图4A至图6B的示例性实施例,当参考回波时间被确定为与早期回波相应时,MRI设备300a或300b可将第一回波数据确定为参考回波数据。
现在将参照图7A至图8B解释参考回波时间与后期回波相应的情况。
图7A是根据示例性实施例的由MRI设备300a或300b施加的ETL为2的脉冲序列的脉冲序列图710。
参照图7A,MRI设备300a或300b可根据脉冲序列图710中所示的多回波序列来施加梯度磁场720。
详细地,脉冲序列图710可示出在一个TR期间由根据示例性实施例的MRI设备300a或300b施加的相位编码方向的梯度磁场PE1、梯度磁场PE2和梯度磁场PE的翻转以及读出方向的梯度磁场720。
当图7A的脉冲序列图710和图5A的脉冲序列图510彼此进行比较时,读出方向的梯度磁场720可以与脉冲序列图510中的读出方向的梯度磁场520相同。在相位编码方向的梯度磁场PE1、梯度磁场PE2和梯度磁场PE的翻转方面,脉冲序列图710与脉冲序列图510不同,因此下面将关注于所述不同。
MRI设备300a或300b可通过在获取第一回波数据E1之前应用相位编码方向的梯度磁场PE1来确定第一回波数据E1在k空间中的ky轴上的位置。
在应用相位编码方向的梯度磁场PE1之后,MRI设备300a或300b可通过使用与第一回波时间TE1相应的读出梯度磁场来沿正方向711填充k空间的第一回波数据E1中的一行数据。第一回波数据E1可不包括k空间中的ky坐标为0的数据。
MRI设备300a或300b可通过在获取第三回波数据E1'之前应用相位编码方向的梯度磁场PE2来确定第三回波数据E1'在k空间中的ky轴上的位置。
MRI设备300a或300b可通过反转读出梯度磁场的符号通过使用与第三回波时间TE1'相应的读出梯度磁场来沿负方向713填充k空间的第三回波数据E1'中的一行数据。
MRI设备300a或300b可通过重新反转读出梯度磁场的符号通过使用与第二回波时间TE2相应的读出梯度磁场来沿正方向714填充k空间的第二回波数据E2中的一行数据。
接下来,MRI设备300a或300b可通过重新反转读出梯度磁场的符号通过使用与第四回波时间TE2'相应的读出梯度磁场来沿负方向715填充k空间的第四回波数据E2'中的一行数据。
接下来,MRI设备300a或300b可应用相位编码方向的梯度磁场PE的翻转。当MRI设备300a或300b应用相位编码方向的梯度磁场PE的翻转时,下一个TR中的k空间的数据可以被放置在ky坐标为0的位置。
如图7A所示,MRI设备300a或300b可在第一回波时间TE1获取不包括k空间的中心部分的数据的第一回波数据E1,并且可以在第二回波时间TE2获取包括k空间的中心部分的数据的第二回波数据E2。
参照图7A,当将获取T2或T2*加权图像时,MRI设备300a或300b可将第二回波时间TE2确定为参考回波时间。在这种情况下,作为参考回波数据的第二回波数据E2可包括k空间中的与第一回波数据E1重叠的部分。
图7B是示出根据示例性实施例的由MRI设备300a或300b执行的根据图7A的脉冲序列图710获取MR图像的处理的示图。
参照图7B,MRI设备300a或300b执行基于与第一回波时间TE1相应的第一回波数据E1和与第二回波时间TE2相应的第二回波数据E2重建最终图像的处理。
根据示例性实施例的MRI设备300a或300b可获取与第一回波时间TE1相应的第一回波数据E1,然后可获取与比第一回波时间TE1晚的第二回波时间TE2相应的第二回波数据E2。此外,MRI设备300a或300b可获取与第三回波时间TE1'相应的第三回波数据E1'和与第四回波时间TE2'相应的第四回波数据E2'。
参照图7B,参考回波数据可以是第二回波数据E2,参考回波时间可以是第二回波时间TE2。当图7B和图5B彼此进行比较时,参考回波数据和参考回波时间存在不同,而图7B和图5B中的其它方面可以是相同的。因此,以下将关注于所述不同。
参照图7B,参考回波数据可以是第二回波数据E2,并且在这种情况下,第二回波数据E2可包括k空间中的与第一回波数据E1重叠的部分。例如,第二回波数据E2还可包括关于k空间中的与第一回波数据E1重叠的部分的附加数据725。附加数据725可占例如整个k空间的大约10%。
根据示例性实施例的MRI设备300a或300b可针对第一回波数据E1执行相位校正。因为当MRI设备300a或300b获取第三回波数据E1'和第四回波数据E2'时在k空间中放置数据的方向713和715(参见图7A)是相同的,所以可最小化由回波时移引起的误差。
MRI设备300a或300b可针对第一回波数据E1执行相位校正,并且可以获取基于第二回波时间TE2校正了相位的第一回波数据E1C。可通过使用B0图来执行相位校正,其中,该B0图是通过使用第三回波数据E1'和第四回波数据E2'产生的。
MRI设备300a或300b可基于第一回波数据E1和校正后的第二回波数据E2C来获取最终的k空间数据E1C+2E。最终的k空间数据E1C+2E可以是在k空间中重建的k空间数据。当重建最终的k空间数据E1C+2E时,可以基于除了附加数据725之外的第二回波数据E2的数据和第一回波数据E1来重建k空间数据。MRI设备300a或300b可以使用GRAPPA方法等来重建k空间数据。接下来,MRI设备300a或300b可以基于包括附加数据725的第二回波数据E2和校正后的第一回波数据E1C重新重建k空间数据。在这种情况下,可以将GRAPPA方法等重新应用于附加数据725,并且MRI设备300a或300b可以获取最终的k空间数据E1C+2E。
图8A是根据示例性实施例的由MRI设备300a或300b施加的ETL为3的脉冲序列的脉冲序列图810。
参照图8A,MRI设备300a或300b可根据脉冲序列图810中所示的多回波序列来施加梯度磁场820。
脉冲序列图810示出了在一个TR期间由根据示例性实施例的MRI设备300a或300b施加的相位编码方向的梯度磁场PE1、梯度磁场PE2和梯度磁场PE的翻转以及读出方向的梯度磁场820。
参照图8A,如脉冲序列图810所示,MRI设备300a或300b可施加用于沿k空间的读出方向产生多个回波的梯度磁场820。图8A的脉冲序列图810与图7A的脉冲序列图710的不同之处在于可另外产生与第五回波时间TE3相应的回波。
当图8A的脉冲序列图810与图7A的脉冲序列图710相互比较时,读出方向的梯度磁场820的ETL为3,而脉冲序列图710中的读出方向的梯度磁场720的ETL为2,并且在读出方向的梯度磁场820和相位编码方向的梯度磁场PE1、梯度磁场PE2和梯度磁场PE的翻转方面存在不同,因此下面将关注于所述不同。
MRI设备300a或300b可以通过使用在一个TR中包括的第一回波时间TE1、第二回波时间TE2和第五回波时间TE3产生的回波来获取第一回波数据E1、第二回波数据E2和第五回波数据E3。MRI设备300a或300b通过获取第一回波数据E1、第二回波数据E2和第五回波数据E3来重建最终的k空间。
首先,MRI设备300a或300b可以通过在获取第一回波数据E1之前应用相位编码方向的梯度磁场PE1来确定第一回波数据E1在k空间中的ky轴上的位置。
在应用相位编码方向的梯度磁场PE1之后,MRI设备300a或300b可以通过使用与第一个回波时间TE1相应的读出梯度磁场沿正方向811填充k空间的第一回波数据E1中的一行数据。
MRI设备300a或300b可通过在获取第二回波数据E2之前应用相位编码方向的梯度磁场PE2来确定第二回波数据E2在k空间中的ky轴上的位置。
MRI设备300a或300b可通过使与回波时间TE2相应的读出梯度磁场的符号与用于获取第一回波数据E1的读出梯度磁场的符号相同来沿正方向813填充k空间的第二回波数据E2中的一行数据。
MRI设备300a或300b可通过在获取第三回波数据E1'之前应用相位编码方向的梯度磁场PE3来确定第三回波数据E1'在k空间中的ky轴上的位置。
MRI设备300a或300b可通过反转读出梯度磁场的符号通过使用与第三回波时间TE1'相应的读出梯度磁场沿负方向815填充k空间的第三回波数据E1'中的一行数据。
MRI设备300a或300b可通过重新反转读出梯度磁场的符号通过使用与第五回波时间TE3相应的读出梯度磁场沿正方向817填充k空间的第五回波数据E3中的一行数据。
接下来,MRI设备300a或300b可通过重新反转读出梯度磁场的符号通过使用与第四回波时间TE2'相应的读出梯度磁场来沿负方向819填充k空间的第四回波数据E2'中的一行数据。
当MRI设备300a或300b应用相位编码方向的梯度磁场PE的翻转时,下一个TR的k空间的数据可以被放置在ky坐标为0的位置。
参照图8A,MRI设备300a或300b可通过在多个TR期间施加梯度磁场820来获取与第一回波时间TE1相应的第一回波数据E1、与第二回波时间TE2相应的第二回波数据E2、与第三回波时间TE1'相应的第三回波数据E1'、与第四回波时间TE2'相应的第四回波数据E2'和与第五回波时间TE3相应的第五回波数据E3。
当MRI设备300a或300b获取第一回波数据E1、第二回波数据E2和第五回波数据E3时在k空间中放置多条数据的方向811、813和817是相同的。也就是说,在第一回波时间TE1施加的读出梯度磁场、在第二回波时间TE2施加的读出梯度磁场和在第五回波时间TE3施加的读出梯度磁场具有相同的符号。
当MRI设备300a或300b获取第三回波数据E1'和第四回波数据E2'时在k空间中放置多条数据的方向815和819是相同的。如果MRI设备300a或300b通过使用第三回波数据E1'和第四回波数据E2'产生B0图,则可以最小化当回波时移的方向彼此不同时引起的误差。
如图8A所示,第五回波数据E3可包括k空间的中心部分的数据,并且第一回波数据E1和第二回波数据E2可不包括k空间的中心部分的数据。
参照图8A,当将获取T2或T2*加权图像时,MRI设备300a或300b可以将第五回波时间TE3确定为参考回波时间。在这种情况下,作为参考回波数据的第五回波数据E3可包括k空间中的与第二回波数据E2重叠的部分。此外,第二回波数据E2可包括k空间中的与第一回波数据E1重叠的部分。
图8B是示出根据示例性实施例的根据图8的脉冲序列图810由MRI设备300a或300b执行以获取MR图像的处理的示图。
参照图8B,MRI设备300a或300b基于与第一回波时间TE1相应的第一回波数据E1、与第二回波时间TE2相应的第二回波数据E2和与第五回波时间TE3相应的第五回波数据E3执行重建最终图像的处理。
根据示例性实施例的MRI设备300a或300b可获取与第一回波时间TE1相应的第一回波数据E1,可获取与比第一回波时间TE1晚的第二回波时间TE2相应的第二回波数据E2,并且可以获取与比第二回波时间TE2晚的第五回波时间TE3相应的第五回波数据E3。此外,MRI设备300a或300b可以获取与第三回波时间TE1'相应的第三回波数据E1'和与第四回波时间TE2'相应的第四回波数据E2'以执行相位校正。
参照图8B,参考回波数据可以是第五回波数据E3,参考回波时间可以是第五回波时间TE3。在这种情况下,第五回波数据E3还可包括关于k空间中的与第二回波数据E2重叠的部分的附加数据825。此外,第二回波数据E2还可包括关于k空间中的与第一回波数据E1重叠的部分的附加数据835。
首先,根据示例性实施例的MRI设备300a或300b可针对第一回波数据E1执行相位校正。如参照图6A所述,因为当MRI设备300a或300b获取第三回波数据E1'和第四回波数据E2'时在k空间中放置多条数据的方向815和819是相同的,所以可最小化由回波时移引起的误差。
MRI设备300a或300b可针对第一回波数据E1执行相位校正,并且可以获取基于第二回波时间TE2校正了相位的第一回波数据E1CE2。可以通过使用B0图来执行相位校正,其中,该B0图是通过使用第三回波数据E1'和第四回波数据E2'产生的。
MRI设备300a或300b针对第一回波数据E1执行第一相位校正,然后针对基于第二回波数据E2和校正后的第一回波数据E1CE2获取的中间k空间数据E2+E1CE2执行第二相位校正。可基于第二回波数据E2的相位来执行第一相位校正,并且可以基于第五回波数据E3的相位来执行第二相位校正。
MRI设备300a或300b可基于第二回波数据E2和校正后的第一回波数据E1CE2获取中间k空间数据E2+E1CE2。如果第二回波数据E2和第五回波数据E3是欠采样的k空间数据,则中间k空间数据E2+E1CE2可以是通过使用GRAPPA方法或SMASH方法重建的数据。
此外,当重建中间k空间数据E2+E1CE2时,MRI设备300a或300b可以基于除了附加数据835之外的第二回波数据E2的数据和第五回波数据E3来重建k空间数据。MRI设备300a或300b可使用GRAPPA方法等来重建k空间数据。接下来,MRI设备可基于包括附加数据835的第二回波数据E2和校正后的第一回波数据E1CE2重新重建k空间数据。在这种情况下,可以将GRAPPA方法等重新应用于附加数据835,并且MRI设备300a或300b可获取中间k空间数据E2+E1CE2。
接下来,MRI设备300a或300b可针对基于第二回波数据E2和校正后的第一回波数据E1CE2获取的中间k空间数据E2+E1CE2执行第二相位校正。
MRI设备300a或300b可针对中间k空间数据E2+E1CE2执行相位校正,并且可获取基于第五回波时间TE3校正了相位的中间k空间数据(E2+E1CE2)CE3。可通过使用B0图来执行相位校正,其中,该B0图是通过使用第三回波数据E1'和第四回波数据E2'产生的。
MRI设备300a或300b可基于第一回波数据E1和校正后的中间k空间数据(E2+E1CE2)CE3获取最终的k空间数据E3+(E2+E1CE2)CE3。最终的k空间数据E3+(E2+E1CE2)CE3可以是在k空间中重建的k空间数据。
由MRI设备300a或300b执行的通过使用附加数据825重建最终k空间数据E3+(E2+E1CE2)CE3的方法与参照图7B描述的方法类似,因此,将不给出其解释。
图9A是示出由MRI设备300a或300b获取的第一回波数据E1和第二回波数据E2的示图。
详细地说,图9A示出了k空间中的欠采样的第一回波数据E1和欠采样的第二回波数据E2。
参照图9A,实际上获取的数据行由实线标记。此外,通过使用GRAPPA方法构建的数据行由虚线标记。此外,可在k空间中的ky坐标接近0的位置获取自动校准信号(ACS)。例如,由粗实线标记的行905可以是ACS行。
此外,图9A的第一核至第四核911、913、915和917中的每一个是具有大小为2*5的核。也就是说,MRI设备300a或300b可以通过使用第一核至第四核911、913、915和917中包括的2*5条数据之间的相关性来使用GRAPPA方法。
参照图9A,第一核911和第四核917可仅包括第一回波数据E1。此外,第二核913可仅包括第二回波数据E2。第三核915可包括第一回波数据E1和第二回波数据E2。
当MRI设备300a或300b基于包括通过使用不同回波获取的多条数据(即,第一回波数据E1和第二回波数据E2)的第三核915使用GRAPPA方法时的重建数据的准确度可比当MRI设备300a或300b基于包括通过使用相同回波获取的多条回波数据的第一核911、第二核913和第四核917使用GRAPPA方法时的重建数据的准确度低。
现在将参照图9B解释使用关于第一回波数据E1和第二回波数据E2之间的重叠部分的附加数据来提高基于第三核915的重建数据的准确度的情况。
图9B是示出根据示例性实施例的通过将GRAPPA方法应用于由MRI设备300a或300b获取的第一回波数据E1和第二回波数据E2来重建k空间数据的处理的示图。
在图9B中,示图910示出了欠采样的第一回波数据E1和欠采样的第二回波数据E2。图9B的示图910示出了除了关于k空间中的与第二回波数据E2重叠的部分的附加数据之外的欠采样的第一回波数据E1的一部分。此外,图9B的示图920示出了欠采样的第一回波数据E1的多条附加数据922和924。
首先,如图9B的示图930所示,MRI设备300a或300b可基于欠采样的第一回波数据E1和欠采样的第二回波数据E2获取重建的k空间数据。例如,MRI设备300a或300b可通过将GRAPPA方法应用于欠采样的第一回波数据E1和欠采样的第二回波数据E2来获取重建的k空间数据。
接下来,如图9B的示图940所示,MRI设备300a或300b可针对示图图930的第二回波数据E2执行相位校正,并且可以获取校正后的第二回波数据E2C。
接下来,如图9B的示图950所示,MRI设备300a或300b可将GRAPPA方法应用于多条附加数据952和954。
图9B的示图950中所示的多条附加数据952和954可与示图920的多条附加数据922和924相应。此外,图9B的示图950的多条附加数据952和954可与示图930的重建的k空间数据中具有低准确度的数据相应。如图9B的示图950所示,当基于多条附加数据952和954重新使用GRAPPA方法时,可提高重建的k空间数据的准确度。
图9B的示图950的附加数据952可以是欠采样的第一回波数据E1。也就是说,附加数据952可包括实际获取的数据行955和应用了GRAPPA方法的数据行953。尽管在图9B中多条附加数据952和954被示为两个行,但是多条附加数据952和954可以占整个k空间数据行的大约10%。与附加数据952一样,图9B的示图950的附加数据954也可以是欠采样的第一回波数据E1。
如示图950所示,当GRAPPA应用于多条附加数据952和954时,考虑到核的大小,可进一步使用位于多条附加数据952和954上边和下边的数据行。例如,图9B的示图940的数据行可以位于多条附加数据952和954的上边和下边。
例如,位于附加数据952下边的数据行951可与图9B的示图940的数据行941相应,位于附加数据952上边的数据行957可与图9B的示图940的数据行947相应。同样地,图9B的示图940的数据行也可以位于附加数据954的上边和下边。
参照图9B的示图960,MRI设备300a或300b可通过将GRAPPA方法应用于示图950的多条附加数据952和954来获取多条重建的附加数据962和964。在将GRAPPA方法应用于图9B的示图960的多条重建的附加数据962和964之后获取的数据行963和969可与示图950的多条附加数据952和954的数据行953和959相应。
参照图9B的示图970,MRI设备300a或300b可以通过用多条重建的附加数据962和964替换与图9B的示图940的数据行中的附加数据相应的第一回波数据942和944来获取最终的k空间数据。
图10A是示出使用多频带方法由MRI设备300a或300b获取的第一回波数据E1和第二回波数据E2的示图。
详细地,图10A示出了在ky轴和kz轴上的通过使用多频带方法获取的第一回波数据E1和第二回波数据E2。未示出kx轴的数据。通过使用多频带方法获取的回波数据是指关于具有不同z坐标的两个切片的k空间数据。例如,通过使用多频带方法获取的回波数据可以是关于两个切片的多条k空间数据的总和,或者是关于两个切片的多条k空间数据之间的差。当通过使用GRAPPA方法重建通过使用多频带方法获取的回波数据时,可以通过重建沿kz轴的两行数据来重建两个切片的最终图像。
参照图10A,部分O是获取数据的部分,部分X是不获取数据的部分。此外,部分A是ACS。
详细地,图10A的部分O是包括通过使用多频带方法获取的关于两个切片的两条数据的部分。
此外,图10A的部分X是必须通过使用GRAPPA方法获取数据的部分。可通过使用包括在例如2*1的核中的多条数据之间的相关性来获取图10A的部分X。
参照图10A,第一核1011可仅包括第一回波数据E1。此外,第二核1015可仅包括第二回波数据E2。第三核1019可包括第一回波数据E1和第二回波数据E2。
当MRI设备300a或300b基于包括通过使用不同回波获取的回波数据(即,第一回波数据E1和第二回波数据E2)的第三核1019使用GRAPPA方法时的重建数据的准确度可比当MRI设备300a或300b基于包括通过使用相同回波获取的回波数据的第一核1011和第二核1015使用GRAPPA方法时的重建数据的准确度低。
现在将参照图10B解释使用关于第一回波数据E1和第二回波数据E2之间的重叠部分的附加数据来提高基于第三核1019重建的数据的准确度的情况。
图10B是示出根据示例性实施例的通过将GRAPPA方法应用于由MRI设备300a或300b使用多频带方法获取的第一回波数据E1和第二回波数据E2来重建k空间数据的处理的示图。
图10B的示图1010示出了通过使用多频带方法获取的第一回波数据E1和第二回波数据E2。图10B的示图1010示出了除了关于k空间中的与第二回波数据E2重叠的部分的附加数据之外的通过使用多频带方法获取的第一回波数据E1的数据。此外,图10B的示图1020示出了通过使用多频带方法获取的第一回波数据E1的多条附加数据1022和1024。
首先,如图10B的示图1030所示,MRI设备300a或300b可基于通过使用多频带方法获取的第一回波数据E1和第二回波数据E2来获取重建的k空间数据。图10B中的部分O是包括通过使用多频带方法获取的关于两个切片的两条数据的部分。
例如,MRI设备300a或300b可通过将GRAPPA方法应用于通过使用多频带方法获取的第一回波数据E1和第二回波数据E2来获取重建的k空间数据。在这种情况下,MRI设备300a或300b可以通过使用尺寸为2*1的核1031来重建关于示图1010的部分X的k空间数据。参照图10B的示图1030,部分o是表示通过使用经由使用相同回波获取的回波数据而重建的k空间数据的部分。此外,部分c是表示通过使用经由使用不同回波获取的回波数据而重建的k空间数据的部分。
接下来,如图10B的示图1040所示,MRI设备300a或300b可针对示图图1030的第二回波数据E2执行相位校正,并且可以获取校正后的第二回波数据E2C。
接下来,如图10B的示图1050所示,MRI设备300a或300b可将GRAPPA方法应用于多条附加数据1052和1054。
图10B的示图1050中所示的多条附加数据1052和1054可与示图1020的多条附加数据1022和1024相应。此外,图10B的多条附加数据1052和1054可与示图1030的重建的k空间数据中具有低准确度的数据相应。如图10B的示图1050所示,当GRAPPA方法重新被应用于多条附加数据1052和1054时,可提高重建的k空间数据的准确度。
图10B的示图1050的多条附加数据1052和1054可以是通过使用多频带方法获取的第一回波数据E1。此外,多条附加数据1052和1054可包括将被应用GRAPPA方法的部分X和实际获取的部分O。尽管在图10B中多条附加数据1052和1054被示为一个行,但是多条附加数据1052和1054可与实际占整个k空间大约10%的数据行相应。
如示图1050所示,当GRAPPA被应用于多条附加数据1052和1054时,考虑到核的大小,可进一步使用位于多条附加数据1052和1054上边和下边的多条数据。例如,图10B的示图1040的多条数据可以位于多条附加数据1052和1054的上边和下边。
例如,位于附加数据1052下边的多条数据1051和1053可与图10B的示图1040的多条数据1041和1043相应,并且位于附加数据1052上边的数据1055可与图10B的示图1040的数据1045相应。同样地,图10B的示图1040的多条数据也可以位于附加数据1054的上边和下边。
参照图10B的示图1060,MRI设备300a或300b可获取通过将GRAPPA方法应用于示图1050的多条附加数据1052和1054而重建的多条附加数据1062和1064。如图10B的示图1060所示,可以通过使用包括在核1061中的多条数据之间的相关性来获取图10B的示图1060的多条重建的附加数据1062和1064。重建的附加数据1062可包括通过使用经由使用相同的回波获取的回波数据而重建的k空间数据,并且重建的附加数据1064可包括通过使用经由使用不同的回波获取的回波数据而重建的k空间数据。
参照图10B的示图1070,MRI设备300a或300b可以通过用多条重建的附加数据1062和1064至少替换图10B的示图1040的多条数据中的包括附加数据的多条数据1042和1044中的多个部分来获取最终的k空间数据。
图11是示出根据示例性实施例的MRI设备的方法的流程图。
根据示例性实施例的方法可由MRI设备300a或300b执行。此外,根据示例性实施例的方法可以是通过使用多回波序列获取MR图像的方法。
在操作S110中,MRI设备300a或300b可获取包括k空间中的重叠部分的第一回波数据E1和第二回波数据E2。
详细地,MRI设备300a或300b可通过使用在第一回波时间TE1产生的回波来获取第一回波数据E1,并且可通过使用在比第一回波时间TE1晚的第二回波时间TE2产生的回波来获取第二回波数据E2。
在操作S120中,MRI设备300a或300b可基于第一回波数据E1和第二回波数据E2来重建MR图像。
图12是示出根据示例性实施例的MRI设备的另一方法的流程图。
根据示例性实施例的方法可由MRI设备300a或300b执行。MRI设备300a或300b针对通过使用不同回波获取的回波数据执行相位校正。
在操作S210中,MRI设备300a或300b可以获取包括k空间中的重叠部分的第一回波数据E1和第二回波数据E2。
在操作S220中,MRI设备300a或300b可基于除了附加数据之外的第一回波数据E1的一部分和第二回波数据E2来重建k空间数据。
在操作S230中,MRI设备300a或300b可针对包括在重建的k空间数据中的第二回波数据E2执行相位校正。
在操作S240中,MRI设备300a或300b可基于第一回波数据E1的附加数据和被校正了相位的第二回波数据E2C重新重建k空间数据。
图13是示出根据示例性实施例的MRI设备的另一方法的流程图。
根据示例性实施例的方法可由MRI设备300a或300b执行。
在操作S310中,MRI设备300a或300b可确定MR图像的特征。
在操作S320中,MRI设备300a或300b可根据确定的MR图像的特征确定第一回波数据E1或第二回波数据E2将包括k空间的中心部分的数据。
在操作S330中,MRI设备300a或300b可获取包括k空间中的重叠部分的第一回波数据E1和第二回波数据E2。
详细地,MRI设备300a或300b可通过使用在第一回波时间TE1产生的回波来获取第一回波数据E1,并且可通过使用在比第一回波时间TE1晚的第二回波时间TE2产生的回波来获取第二回波数据E2。
在操作S340中,MRI设备300a或300b可基于第一回波数据E1和第二回波数据E2重建k空间数据。
图14是示出根据示例性实施例的MRI设备的另一方法的流程图。
根据示例性实施例的方法可由MRI设备300a或300b执行。
在操作S410中,MRI设备300a或300b可获取第一回波数据E1、第二回波数据E2和第五回波数据E3。
在操作S420中,MRI设备300a或300b可针对第五回波数据E3执行相位校正。
在操作S430中,MRI设备300a或300b可基于第二回波数据E2和相位被校正了的第五回波数据E3CE2获取中间k空间数据E2+E3CE2。
在操作S440中,MRI设备300a或300b可针对中间k空间数据E2+E3CE2执行相位校正以获取校正后的中间k空间数据(E2+E3CE2)CE1。
在操作S450中,MRI设备300a或300b可基于第一回波数据E1和校正后的中间k空间数据(E2+E3CE2)CE1来重建k空间数据。
图15是一般MRI系统的框图。参照图15,一般MRI系统可包括台架20、信号收发器30、监视器40、系统控制器50和操作部分60。
台架20防止由主磁体22、梯度线圈24和RF线圈26产生的电磁波的向外发射。在台架20的孔中形成静磁场和梯度磁场,并且RF信号朝向对象10发射。
主磁体22、梯度线圈24和RF线圈26可沿台架20的预定方向被布置。预定方向可以是同轴圆筒方向。对象10可被放置在能够沿着圆筒的水平轴插入圆筒中的工作台28上。
主磁体22产生静磁场或静态磁场,静磁场或静态磁场用于沿恒定方向排列对象10的原子核的磁偶极矩。由于主磁体22产生的磁场强且均匀,因此可以获取物体10的精确且准确的MR图像。
梯度线圈24包括X、Y和Z线圈,X、Y和Z线圈用于沿X、Y和Z轴方向产生以直角彼此相交的梯度磁场。梯度线圈24可通过根据对象10的区域不同地诱导共振频率来提供对象10的每个区域的位置信息。
RF线圈26可朝向患者发射RF信号并接收从患者发出的MR信号。详细地,RF线圈26可朝向包括在患者体内并具有进动运动的原子核发射具有与所述进动运动的频率相同频率的RF信号,停止发射RF信号,然后接收从包括在患者体内的原子核发出的MR信号。
例如,为了将原子核从低能量状态转换到高能量状态,RF线圈26可产生电磁波信号,并将所述电磁波信号施加到对象10,其中,所述电磁波信号是与原子核类型相应的RF信号。当由RF线圈26产生的电磁波信号施加到原子核时,原子核可以从低能量状态转换到高能量状态。然后,当由RF线圈26产生的电磁波消失时,被施加电磁波的原子核从高能量状态转换到低能量状态,从而发射具有拉莫频率的电磁波。换句话说,当停止向原子核施加电磁波信号时,原子核的能级从高能级变为低能级,因此原子核可以发射具有拉莫频率的电磁波。RF线圈26可接收来自包括在对象10中的原子核的电磁波信号。
RF线圈26可被实现为一个RF发送和接收线圈,该线圈具有产生电磁波的功能(其中,每个电磁波具有与原子核的类型相应的RF),并且具有接收从原子核发射的电磁波的功能。可选地,RF线圈26可被实现为具有产生电磁波的功能的发送RF线圈(其中,每个电磁波具有与原子核的类型相应的RF)和具有接收从原子核发射的电磁波的功能的接收RF线圈。
RF线圈26可被固定到台架20或可以是可拆卸的。当RF线圈26是可拆卸的时,RF线圈26可以是针对对象10的一部分的RF线圈,诸如头部RF线圈、胸部RF线圈、腿部RF线圈、颈部RF线圈、肩部RF线圈、手腕RF线圈或脚踝RF线圈。
RF线圈26可经由导线和/或无线地与外部装置进行通信,并且还可以根据通信频带执行双调谐通信。
根据结构,RF线圈26可以是鸟笼线圈、表面线圈或横向电磁(TEM)线圈。
根据发送和接收RF信号的方法,RF线圈26可以是发射专用线圈、接收专用线圈或发射和接收线圈。
RF线圈26可以是具有各种数量的通道(诸如16通道、32通道、72通道和144通道)的RF线圈。
台架20还可包括布置在台架20外部的显示器29和布置在台架20内部的显示器。台架20可通过分别布置在台架20外部和内部的显示器29和显示器向用户或对象10提供预定信息。
信号收发器30可根据预定的MR序列控制形成在台架20内部(即,在孔中)的梯度磁场,并控制RF信号和MR信号的发送和接收。
信号收发器30可包括梯度放大器32、发射和接收开关34、RF发射器36和RF接收器38。
梯度放大器32驱动包括在台架20中的梯度线圈24,并且可以在梯度磁场控制器54的控制下向梯度线圈24提供用于产生梯度磁场的脉冲信号。通过控制从梯度放大器32提供给梯度线圈24的脉冲信号,可将沿X、Y和Z轴方向的梯度磁场进行合成。
RF发射器36和RF接收器38可驱动RF线圈26。RF发射器36可以以拉莫频率向RF线圈26提供RF脉冲,并且RF接收器38可接收由RF线圈接收的MR信号。
发射和接收开关34可调整RF信号和MR信号的发射和接收方向。例如,发射和接收开关34可以在发射模式期间通过RF线圈26向对象10发射RF信号,并且在接收模式期间通过RF线圈26从对象10接收MR信号。可由RF控制器56输出的控制信号来控制发射和接收开关34。
图15的信号收发器30包括图3A的数据获取器310或图3B的数据获取器315。例如,图3A的数据获取器310或图3B的数据获取器315可被连接到包括在信号收发器30中的RF接收器38,并且可以从RF接收器38接收MR信号。
根据示例性实施例,信号收发器30可以通过使用通过RF接收器38接收的MR信号来获取第一回波数据E1和第二回波数据E2。
监视器40可以监视或控制台架20或安装在台架20上的装置。监视器40可包括系统监视器42、对象监视器44、工作台控制器46和显示器控制器48。
系统监视器42可监视和控制静磁场的状态、梯度磁场的状态、RF信号的状态、RF线圈26的状态、工作台28的状态、测量对象10的身体信息的装置的状态、供电状态、热交换器的状态和压缩机的状态。
对象监视器44监视对象10的状态。详细地,对象监视器44可以包括用于观察对象10的移动或位置的像机、用于测量对象10的呼吸的呼吸测量器、用于测量对象10的电活动的心电图(ECG)测量器、或用于测量对象10的温度的温度测量器。
工作台控制器46控制对象10所在的工作台28的移动。工作台控制器46可以根据系统控制器50的顺序控制来控制工作台28的移动。例如,在对象10的移动成像期间,工作台控制器46可以根据系统控制器50的顺序控制连续地或不连续地移动工作台28,因此可以以比台架20的视场(FOV)大的FOV对对象进行成像。
显示器控制器48控制设置在台架20外部的显示器29和设置在台架20内部的显示器。详细地,显示控制器48可以控制显示器29和显示器打开或关闭,并且可控制屏幕图像输出在显示器29和显示器上。此外,当扬声器位于台架20的内部或外部时,显示器控制器48可以控制扬声器打开或关闭,或者可以控制声音经由扬声器输出。
系统控制器50可包括:序列控制器52,用于控制在台架20中形成的信号序列;台架控制器58,用于控制台架20和安装在台架20上的装置。
序列控制器52可包括:梯度磁场控制器54,用于控制梯度放大器32;RF控制器56,用于控制RF发送器36、RF接收器38以及发射和接收开关34。序列控制器52可根据从操作部60接收的脉冲序列控制梯度放大器32、RF发送器36、RF接收器38和发射和接收开关34。这里,脉冲序列包括用于控制梯度放大器32、RF发送器36、RF接收器38以及发射和接收开关34的所有信息。例如,脉冲序列可包括关于施加到梯度线圈24的脉冲信号的强度、施加时间以及施加时序的信息。
在控制MRI系统的整体操作的同时,操作部60可请求系统控制器50发送脉冲序列信息。
操作部60可包括:图像处理器62,用于接收和处理由RF接收器38接收的MR信号;输出接口64;输入接口66。
图像处理器62可处理从RF接收器38接收的MR信号以产生对象10的MR图像数据。
图像处理器62接收由RF接收器38接收的MR信号,并对接收的MR信号执行各种信号处理(诸如,放大、频率变换、相位检测、低频放大和滤波)中的任何一个或任何组合。
图像处理器62可将k空间(例如,也称为傅立叶空间或频率空间)中的数字数据布置在存储器中,并且经由2D或3D傅里叶变换将数字数据重新排列为图像数据。
图像处理器62可针对图像数据执行合成处理或差分计算处理。合成处理可包括关于像素的加法处理或最大强度投影(MIP)处理。图像处理器62不仅可以在存储器或外部服务器中存储重新排列的图像数据,还可存储被执行了合成处理或差分计算处理的图像数据。
图像处理器62可并行地对MR信号执行任何信号处理。例如,图像处理器62可并行地针对由多通道RF线圈接收的多个MR信号执行信号处理,以将多个MR信号重新排列成图像数据。
图15的图像处理器62包括图3A的图像处理器320或图3B的图像处理器325。例如,图像处理器62可通过使用由信号收发器30获取的第一回波数据E1和第二回波数据E2来重建MR图像。
输出接口64可将由图像处理器62产生或重新排列的图像数据输出给用户。输出接口64还可输出被用于用户操纵MRI系统的信息,诸如UI、用户信息或对象信息。输出接口64可以是扬声器、打印机、CRT显示器、LCD、PDP显示器、OLED显示器、FED、LED显示器、VFD、DLP显示器、FPD、3D显示器、透明显示器或本领域普通技术人员公知的其他各种输出装置中的任何一个或任何组合。
图15的输出接口64可包括图3B的输出接口340。
用户可通过使用输入接口66输入对象信息、参数信息、扫描条件、脉冲序列或关于图像合成或差分计算的信息。输入接口66可以是键盘、鼠标、轨迹球、语音识别器、手势识别器、触摸屏或本领域普通技术人员公知的其它各种输入装置中的任何一个或任何组合。
信号收发器30、监视器40、系统控制器50和操作部60在图15中是分开的组件,但是,信号收发器30、监视器40、系统控制器50和操作部60的各自的功能可由另一组件执行。例如,图像处理器62将从RF接收器38接收的MR信号转换为数字信号,但是可选地,可以由RF接收器38或RF线圈26执行MR信号到数字信号的转换。
台架20、RF线圈26、信号收发器30、监视器40、系统控制器50和操作部60可以通过线缆或无线地彼此连接,并且当它们无线连接时,MRI系统可还包括用于在它们之间同步时钟信号的设备。台架20、RF线圈26、信号收发器30、监视器40、系统控制器50和操作部60之间的通信可通过使用高速数字接口(诸如,低压差分信令(LVDS))、异步串行通信(诸如,通用异步接收器发送器(UART))、低延迟网络协议(诸如,错误同步串行通信或控制器域网(CAN))或光通信来执行。
图16是根据示例性实施例的通信接口70的框图。操作图16,通信接口70可被连接到从图15的台架20、信号收发器30、监视器40、系统控制器50和操作部60中选择的至少一个。
通信接口70可以向医院服务器或医院中的其他医疗设备发送数据和从医院服务器或医院中的其他医疗设备接收数据(其中,所述医院服务器或医院中的其他医疗设备通过图片存档和通信系统(PACS)被连接),并根据医学数字影像和通讯(DICOM)标准执行数据通信。
如图16所示,通信接口70可通过线缆或无线地连接到网络80以与服务器92、医疗设备94或便携式装置96进行通信。
详细地,通信接口70可通过网络80发送和接收与对象的诊断有关的数据,并且还可以发送和接收由医疗设备94(诸如,CT设备、MRI设备或X射线设备)捕获的医学图像。另外,通信接口70可从服务器92接收对象的诊断历史或治疗方案,并使用该诊断历史或治疗方案来诊断对象。通信接口70不仅可与医院中的服务器92或医疗设备94进行数据通信,还可以与便携式装置96(诸如,移动电话、个人数字助理(PDA)或者医生或病人的笔记本计算机)进行数据通信。
此外,通信接口70可通过网络80向用户发送关于MRI系统的故障或关于医学图像质量的信息,并且从用户接收关于该信息的反馈。
通信接口70可包括能够与外部设备进行通信的至少一个组件。
例如,通信接口70可包括局域通信接口72、有线通信接口74和无线通信接口76。局域通信接口72表示用于与预定距离内的设备进行局域通信的接口。根据示例性实施例的局域通信技术的示例包括但不限于无线局域网(LAN)、Wi-Fi、蓝牙、ZigBee、Wi-Fi直连(WFD)、超宽带(UWB)、红外数据协会(IrDA)、蓝牙低能耗(BLE)和近场通信(NFC)。
有线通信接口74表示用于通过使用电信号或光信号执行通信的接口。根据示例性实施例的有线通信技术的示例包括使用双股电缆、同轴电缆和光纤电缆的有线通信技术以及其它公知的有线通信技术。
无线通信接口76向从移动通信网络中的基站、外部设备和服务器中选择的至少一个发送无线信号并且从移动通信网络中的基站、外部设备和服务器中选择的至少一个接收无线信号。这里,根据文本/多媒体消息的发送和接收,无线信号可以是语音呼叫信号、视频呼叫信号或以各种格式中的任何一种或任何组合的数据。
连接到MRI系统的服务器92、医疗设备94或便携式装置96可以是例如图3A或图3B的MRI设备300a或300b。也就是说,图16的通信接口70可被连接到MRI设备300a或300b。
根据一个或更多个实施例的MRI设备可通过获取k空间中的关于在不同回波时间获取的第一回波数据和第二回波数据之间的重叠部分的附加数据来改善重建的最终图像的质量。
此外,根据一个或多个实施例的MRI设备可通过对在不同回波时间获取的第一回波数据和第二回波数据执行相位校正来减少当使用在不同回波时间获取的多条数据时在最终图像中可能出现的模糊或混叠的影响。
本公开的上述实施例可以被编写为计算机程序,并且可使用计算机可读记录介质在执行所述程序的通用数字计算机中实现。
计算机可读记录介质的示例包括磁存储介质(例如,只读存储器(ROM)、软盘、硬盘等)、光学记录介质(例如,紧凑盘(CD-ROM)或数字通用盘(DVD)等)以及诸如互联网传输媒介的传输媒介。
虽然已经参照本公开的实施例示出并描述了本公开,但是本领域普通技术人员将理解,在不脱离本发明精神的情况下,可以在其中进行形式和细节上的各种改变,因此,上述实施例及其所有方面仅是示例而非限制的。

Claims (12)

1.一种用于基于多回波序列获取磁共振MR图像的磁共振成像MRI设备,所述MRI设备包括:
数据获取器,被配置为:基于在第一回波时间产生的回波获取第一回波数据,并基于在比第一回波时间晚的第二回波时间产生的回波获取第二回波数据,其中,第一回波数据包括k空间中与包括在第二回波数据中的部分重叠的部分;以及
图像处理器,被配置为:
基于第三回波数据和第四回波数据产生B0图,其中,第三回波数据和第四回波数据是分别基于在第三回波时间产生的回波和在第四回波时间产生的回波而被获取的;
基于产生的B0图,对第一回波数据和第二回波数据中的至少一个回波数据执行相位校正;以及
基于第一回波数据和第二回波数据中的被执行了相位校正的至少一个回波数据重建MR图像。
2.如权利要求1所述的MRI设备,其中,第一回波数据和第二回波数据中的至少一个回波数据包括k空间的中心部分的数据。
3.如权利要求2所述的MRI设备,其中,图像处理器还被配置为基于第一回波数据和第二回波数据中的k空间的中心部分的数据的相位执行相位校正。
4.如权利要求1所述的MRI设备,还包括:梯度磁场控制器,被配置为基于多回波序列控制施加的梯度磁场,
其中,在第一回波时间施加的读出梯度磁场的符号与在第二回波时间施加的读出梯度磁场的符号相同。
5.如权利要求1所述的MRI设备,其中,数据获取器还被配置为:响应于第一回波数据包括k空间的中心部分的数据并且第二回波数据不包括k空间的中心部分的数据,基于在第一回波时间产生的回波,获取包括在第一回波数据中的所述重叠的部分的附加数据。
6.如权利要求1所述的MRI设备,其中,数据获取器还被配置为:响应于第二回波数据包括k空间的中心部分的数据并且第一回波数据不包括k空间的中心部分的数据,基于在第二回波时间产生的回波,获取包括在第二回波数据中的所述重叠的部分的附加数据。
7.如权利要求1所述的MRI设备,其中,数据获取器还被配置为:
确定将被获取的MR图像的特征;以及
基于确定的MR图像的特征,确定第一回波数据或第二回波数据将包括k空间的中心部分的数据。
8.如权利要求1所述的MRI设备,还包括:梯度磁场控制器,被配置为基于多回波序列控制施加的梯度磁场,
其中,在第三回波时间施加的读出梯度磁场的符号与在第四回波时间施加的读出梯度磁场的符号相同。
9.如权利要求1所述的MRI设备,其中,第一回波时间、第二回波时间、第三回波时间和第四回波时间被包括在重复时间段中。
10.如权利要求1所述的MRI设备,其中,图像处理器还被配置为:
基于k空间中的除包括在第一回波数据中的所述重叠的部分之外的第一回波数据的另一部分以及第二回波数据重建k空间数据;
对包括在重建的k空间数据中的第二回波数据执行相位校正;以及
基于包括在第一回波数据中的所述重叠的部分和执行了相位校正的第二回波数据,重建k空间数据。
11.一种用于基于多回波序列获取磁共振MR图像的磁共振成像MRI设备的方法,所述方法包括:
基于在第一回波时间产生的回波获取第一回波数据;
基于在比第一回波时间晚的第二回波时间产生的回波,获取第二回波数据,其中,第一回波数据包括k空间中与包括在第二回波数据中的部分重叠的部分;基于第三回波数据和第四回波数据产生B0图,其中,第三回波数据和第四回波数据是分别基于在第三回波时间产生的回波和在第四回波时间产生的回波而被获取的;
基于产生的B0图,对第一回波数据和第二回波数据中的至少一个回波数据执行相位校正;以及
基于第一回波数据和第二回波数据中的被执行了相位校正的至少一个回波数据来重建MR图像。
12.如权利要求11所述的方法,其中,第一回波数据和第二回波数据中的至少一个回波数据包括k空间的中心部分的数据。
CN201680082407.0A 2016-02-17 2016-12-19 磁共振成像设备及其方法 Active CN108697369B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
KR1020160018549A KR101819030B1 (ko) 2016-02-17 2016-02-17 자기 공명 영상 장치 및 자기 공명 영상 장치를 위한 방법
KR10-2016-0018549 2016-02-17
PCT/KR2016/014864 WO2017142180A1 (en) 2016-02-17 2016-12-19 Magnetic resonance imaging apparatus and method thereof

Publications (2)

Publication Number Publication Date
CN108697369A CN108697369A (zh) 2018-10-23
CN108697369B true CN108697369B (zh) 2021-10-26

Family

ID=59562044

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680082407.0A Active CN108697369B (zh) 2016-02-17 2016-12-19 磁共振成像设备及其方法

Country Status (5)

Country Link
US (1) US10288713B2 (zh)
EP (1) EP3402399B1 (zh)
KR (1) KR101819030B1 (zh)
CN (1) CN108697369B (zh)
WO (1) WO2017142180A1 (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102017203936A1 (de) * 2017-03-09 2018-09-13 Siemens Healthcare Gmbh Leise Echo Planar Bildgebung
US10802096B2 (en) * 2017-12-26 2020-10-13 Uih America, Inc. Methods and systems for magnetic resonance imaging
US10884087B2 (en) * 2019-03-01 2021-01-05 Canon Medical Systems Corporation Apparatuses and methods to improve chemical species separation in magnetic resonance imaging
CN113391251B (zh) * 2020-03-12 2023-05-26 上海联影医疗科技股份有限公司 磁共振图像重建方法、装置和设备
JP2023046466A (ja) * 2021-09-24 2023-04-05 富士フイルムヘルスケア株式会社 磁気共鳴イメージング装置、および、画像処理方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5539311A (en) * 1992-12-03 1996-07-23 Hitachi Medical Corporation Method for reducing artifacts in magnetic resonance imaging
US5742163A (en) * 1996-04-26 1998-04-21 Picker International, Inc. Magnetic resonance scan calibration and reconstruction technique for multi-shot, multi-echo imaging
EP0845684A1 (en) * 1996-11-27 1998-06-03 Picker International, Inc. Phase Correction in Multiple-Echo Magnetic Resonance Imaging
CN101051074A (zh) * 2006-04-04 2007-10-10 株式会社东芝 磁共振成像装置和磁共振成像方法

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL121775A0 (en) * 1997-09-15 1998-02-22 Elscint Ltd Removing discontinuities in K-space data
DE10015068C2 (de) 2000-03-25 2002-06-27 Bruker Medical Gmbh Verfahren zur Erzeugung von Bildern der magnetischen Resonanz
JP4763142B2 (ja) 2001-03-02 2011-08-31 株式会社日立メディコ 磁気共鳴イメージング装置
JP2006061235A (ja) 2004-08-25 2006-03-09 Hitachi Medical Corp 磁気共鳴イメージング装置
JP5179182B2 (ja) 2005-07-27 2013-04-10 株式会社日立メディコ 磁気共鳴イメージング装置
WO2009081785A1 (ja) * 2007-12-25 2009-07-02 Hitachi Medical Corporation 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JP5591493B2 (ja) * 2008-07-17 2014-09-17 株式会社東芝 磁気共鳴イメージング装置
DE102009004896B4 (de) * 2009-01-16 2010-10-28 Siemens Aktiengesellschaft Magnetresonanzgerät und Verfahren zur Messung von Feldinhomogenitäten und Einstellung von Shim-Parametern
JP5828763B2 (ja) * 2009-07-16 2015-12-09 株式会社日立メディコ 磁気共鳴イメージング装置および磁気共鳴イメージング方法
JP5835989B2 (ja) 2010-08-11 2015-12-24 株式会社東芝 磁気共鳴イメージング装置および磁気共鳴イメージング方法
JP5925529B2 (ja) * 2011-03-31 2016-05-25 株式会社東芝 磁気共鳴イメージング装置
US9612307B2 (en) * 2011-04-21 2017-04-04 The Board Of Trustees Of The University Of Illinois Reducing artifacts in magnetic resonance images acquired using propeller echo planar imaging
CN103533887B (zh) * 2011-05-20 2016-01-20 株式会社日立医疗器械 磁共振成像装置、以及重构图像取得方法
DE102011077197B4 (de) * 2011-06-08 2013-05-16 Siemens Aktiengesellschaft Verzeichnungskorrektur bei einer Magnetresonanz-Bildgebung
DE102012208425B4 (de) 2012-05-21 2013-12-12 Siemens Aktiengesellschaft Kontinuierliches Korrigieren von Phasenfehlern einer multidimensionalen, ortsselektiven Magnetresonanz-Messsequenz
EP2762910A1 (en) * 2013-01-30 2014-08-06 Samsung Electronics Co., Ltd Susceptibility-weighted magnetic resonance imaging
DE102013205208B4 (de) 2013-03-25 2015-02-12 Siemens Aktiengesellschaft Multiecho-Magnetresonanz-Messsequenz mit erhöhter Auflösung
DE102013205785B4 (de) * 2013-04-02 2014-10-09 Siemens Aktiengesellschaft Ermittlung einer Magnetresonanzsystem-Ansteuersequenz auf Basis einer reduzierten Anzahl an Feldverteilungskarten
DE102013220012B4 (de) * 2013-10-02 2015-04-30 Siemens Aktiengesellschaft Bildrekonstruktion mit mehreren virtuellen Spulen
US9715726B2 (en) 2013-12-05 2017-07-25 Siemens Healthcare Gmbh Method and system for B0 drift and respiratory motion compensation in echo-planar based magnetic resonance imaging
US10234521B2 (en) * 2014-04-24 2019-03-19 Junmin LIU Systems and methods for field mapping in magnetic resonance imaging

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5539311A (en) * 1992-12-03 1996-07-23 Hitachi Medical Corporation Method for reducing artifacts in magnetic resonance imaging
US5742163A (en) * 1996-04-26 1998-04-21 Picker International, Inc. Magnetic resonance scan calibration and reconstruction technique for multi-shot, multi-echo imaging
EP0845684A1 (en) * 1996-11-27 1998-06-03 Picker International, Inc. Phase Correction in Multiple-Echo Magnetic Resonance Imaging
CN101051074A (zh) * 2006-04-04 2007-10-10 株式会社东芝 磁共振成像装置和磁共振成像方法

Also Published As

Publication number Publication date
US10288713B2 (en) 2019-05-14
EP3402399A1 (en) 2018-11-21
CN108697369A (zh) 2018-10-23
KR101819030B1 (ko) 2018-01-16
WO2017142180A1 (en) 2017-08-24
US20170234957A1 (en) 2017-08-17
EP3402399A4 (en) 2019-03-20
EP3402399B1 (en) 2023-11-22
KR20170096818A (ko) 2017-08-25

Similar Documents

Publication Publication Date Title
US10185013B2 (en) Magnetic resonance imaging (MRI) apparatus and method of generating MR image
US10254365B2 (en) Magnetic resonance imaging apparatus and image processing method thereof
US10151819B2 (en) Magnetic resonance imaging apparatus and method of scanning blood vessel
US10168405B2 (en) Method and apparatus for quantifying properties of an object through magnetic resonance imaging (MRI)
CN108697369B (zh) 磁共振成像设备及其方法
US10473742B2 (en) Magnetic resonance imaging apparatus and method of generating magnetic resonance image by using the same
US9478047B2 (en) Apparatus and method for reconstructing images by displaying user interface indicating image reconstruction modes
US10274563B2 (en) Magnetic resonance imaging apparatus and method
US20170234955A1 (en) Method and apparatus for reconstructing magnetic resonance image
US10466323B2 (en) Magnetic resonance imaging apparatus and method
US10213131B2 (en) Method of generating magnetic resonance image and medical imaging apparatus using the method
US10145916B2 (en) Magnetic resonance imaging apparatus and method of scanning magnetic resonance image using the same
US9927508B2 (en) Magnetic resonance imaging apparatus and method for operating the same
US9977109B2 (en) Magnetic resonance imaging apparatus and operating method for the same
US20150374247A1 (en) Method of measuring blood flow velocity performed by medical imaging apparatus, and the medical imaging apparatus
US10408900B2 (en) Method and apparatus for processing MRI images
US9523754B2 (en) Image processing method and medical imaging apparatus employing the method
EP3409192A1 (en) Magnetic resonance imaging apparatus and magnetic resonance image acquisition method thereof
US20180321348A1 (en) Magnetic resonance imaging apparatus and method therefor

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