CN111542762A - 双分辨率Dixon磁共振成像 - Google Patents

双分辨率Dixon磁共振成像 Download PDF

Info

Publication number
CN111542762A
CN111542762A CN201980006440.9A CN201980006440A CN111542762A CN 111542762 A CN111542762 A CN 111542762A CN 201980006440 A CN201980006440 A CN 201980006440A CN 111542762 A CN111542762 A CN 111542762A
Authority
CN
China
Prior art keywords
magnetic resonance
resolution
point dixon
resonance imaging
pulse sequence
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.)
Granted
Application number
CN201980006440.9A
Other languages
English (en)
Other versions
CN111542762B (zh
Inventor
H·埃格斯
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of CN111542762A publication Critical patent/CN111542762A/zh
Application granted granted Critical
Publication of CN111542762B publication Critical patent/CN111542762B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4828—Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/243—Spatial mapping of the polarizing magnetic field
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/36—Electrical details, e.g. matching or coupling of the coil to the receiver
    • G01R33/3621—NMR receivers or demodulators, e.g. preamplifiers, means for frequency modulation of the MR signal using a digital down converter, means for analog to digital conversion [ADC] or for filtering or processing of the MR signal such as bandpass filtering, resampling, decimation or interpolation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image 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/5615—Echo 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]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image 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/5615—Echo 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/5616—Echo 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

Landscapes

  • Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Signal Processing (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明提供了一种磁共振成像系统(100)。机器可执行指令(140)使控制所述磁共振成像系统的处理器：利用脉冲序列命令控制(200)所述磁共振成像系统采集两点Dixon磁共振数据和单点Dixon磁共振数据；使用所述两点Dixon磁共振数据来计算(202)第一分辨率磁场不均匀性图(148)；通过将所述第一分辨率磁性不均匀性图内插到所述第二分辨率来计算(204)第二分辨率磁场不均匀性图(154)；并且使用所述单点Dixon磁共振成像数据和所述第二分辨率磁场不均匀性图来计算(206)第二分辨率水图像(156)和第二分辨率脂肪图像(158)。所述第一分辨率低于所述第二分辨率。

Description

双分辨率Dixon磁共振成像
技术领域
本发明涉及磁共振成像，具体涉及Dixon磁共振成像。
背景技术
磁共振成像(MRI)扫描器使用大的静态磁场来对准原子的核自旋，作为用于产生患者体内的图像的流程的部分。该大的静态磁场被称为B0场。
在MRI扫描期间，由发射器线圈生成的射频(RF)脉冲引起对局部磁场的扰动，并且由核自旋发射的RF信号由接收器线圈进行检测。这些RF信号用于构建MRI图像。这些线圈也能够被称为天线。另外，发射器线圈和接收器线圈也能够被集成到执行两种功能的单个收发器线圈中。应理解，对术语收发器线圈的使用还指代使用单独的发射器线圈和接收器线圈的系统。所发射的RF场被称为B1场。
MRI扫描器能够构建切片或体积的图像。切片是仅一个体素厚的薄体积。体素是MRI信号在其上被平均的小体积，并且表示MRI图像的分辨率。体素在本文中也可以被称为像素。
磁共振成像的Dixon方法包括用于产生单独的水图像和脂质(脂肪)图像的一系列技术。各种Dixon技术(例如，但不限于，两点Dixon方法、三点Dixon方法和六点Dixon方法)在本文中被统称为Dixon技术或方法。对水图像和脂肪图像的重建依赖于对很大程度上由于B0不均匀性的相位误差的准确确定，以防止水信号和脂肪信号的交换。
Yu等人的期刊文章“SingleAcquisitionWater-FatSeparation:FeasibilityStudyforDynamicImaging”(MagneticResonanceinMedicine55:413–422(2006))公开了假定B0不均匀性的先验知识的单回波Dixon方法。
发明内容
本发明在独立权利要求中提供了一种磁共振成像系统、一种计算机程序产品以及一种方法。在从属权利要求中给出了实施例。
梯度回波Dixon成像通常与由于扫描时间、空间分辨率、回波时间和接收器带宽约束的高声压级相关联。
实施例可以通过将较低分辨率(第一分辨率)两点Dixon磁共振成像协议与较高分辨率(第二分辨率)单点Dixon磁共振成像协议合并来减少所生成的声学噪声量。两点Dixon磁共振成像协议用来确定第一分辨率磁场不均匀性图。第一分辨率磁场不均匀性图然后被内插到第二分辨率以创建第二分辨率磁场不均匀性图。第二分辨率磁场不均匀性图然后实现对单点Dixon磁共振成像协议的使用。
在一个方面中，本发明提供了一种磁共振成像系统，包括用于存储机器可执行指令和脉冲序列命令的存储器。如本文中所使用的脉冲序列命令包含命令或能够被转换成这种命令的数据，所述命令可以用来控制磁共振成像系统采集磁共振成像数据。例如脉冲序列命令可以采取图示各种动作何时由磁共振成像系统的各种部件执行的时序图的形式。所述脉冲序列命令被配置为根据两点Dixon磁共振成像协议从感兴趣区域采集两点Dixon磁共振数据。所述脉冲序列命令被配置为根据单点Dixon磁共振成像协议从所述感兴趣区域采集单点Dixon磁共振数据。所述两点Dixon磁共振成像协议被配置为从中心k空间区域对所述两点Dixon磁共振数据进行采样。
所述单点Dixon磁共振成像协议被配置为从扩展的k空间区域和所述中心k空间区域对所述单点Dixon磁共振数据进行采样。在一些情况下，特定数据点可以属于单点Dixon磁共振数据和两点Dixon磁共振数据两者。
所述扩展的k空间区域大于所述中心k空间区域。扩展的k空间区域与中心k空间区域的组合可以被视为被采样的完整k空间。所述扩展的k空间区域至少部分地围绕所述中心k空间区域。所述两点Dixon磁共振成像协议被配置用于以第一分辨率生成所述感兴趣区域的图像。所述单点Dixon磁共振成像协议被配置用于以第二分辨率生成所述感兴趣区域的图像。所述第二分辨率高于所述第一分辨率。
所述磁共振成像系统还包括用于控制所述磁共振成像系统的处理器。对所述机器可执行指令的执行使所述处理器利用所述脉冲序列命令控制所述磁共振成像系统采集所述两点Dixon磁共振数据和所述单点Dixon磁共振数据。对所述机器可执行指令的执行还使所述处理器使用所述两点Dixon磁共振数据来计算第一磁场不均匀性图。对所述机器可执行指令的执行还使所述处理器通过将所述第一分辨率磁性不均匀性图内插到所述第二分辨率来计算第二磁场不均匀性图。对所述机器可执行指令的执行还使所述处理器使用所述单点Dixon磁共振数据和所述第二分辨率磁场不均匀性图来计算第二分辨率水图像和第二分辨率脂肪图像。
该实施例可以是有益的，因为相比于以第二分辨率采集两点Dixon磁共振数据，以较低分辨率采集两点Dixon磁共振数据可以具有减少当两点Dixon磁共振数据被采集时出现的声学噪声量的效果。执行两点Dixon磁共振成像协议的优先是磁场不均匀性图能够被更可靠地确定。因此在该实施例中，两点Dixon磁共振成像协议以较低分辨率被执行，以降低噪声并且以第一分辨率获得磁场不均匀性图。这然后被内插，并且用作用于处理单点Dixon磁共振数据的输入。对第二磁场不均匀性图的使用降低体素被不正确地识别为主要包含水或脂肪的可能性。
在另一实施例中，所述脉冲序列命令被配置为使用双极性双回波脉冲序列从中心k空间区域对所述两点Dixon磁共振数据进行采样。所述第一回波是利用具有第一极性的第一读出梯度旁瓣来生成的，并且所述第二回波是利用具有第二极性的第二读出梯度旁瓣来生成的。
在另一实施例中，所述第一读出梯度旁瓣具有第一幅度和第一持续时间。所述第二读出梯度旁瓣具有第二幅度和第二持续时间。所述第一持续时间乘以所述第一幅度小于所述第二持续时间乘以所述第二幅度。该实施例具有使第一分辨率小于第二分辨率的效果。这例如可以是有益的，因为它可以减少当磁共振数据被采集时生成的噪声量。
如本文中所使用的术语“第一读出梯度旁瓣”和“第二读出梯度旁瓣”是用于特定读出梯度旁瓣的标签。“第一读出梯度旁瓣”可以依据实施方式而在“第二读出梯度旁瓣”之前或之后被执行。
在另一实施例中，第一幅度和第二幅度是相同的，这可以是有益的，因为它可以实现梯度幅度的最小化，这可以减少声学噪声的生成。
在另一实施例中，例如，第一幅度和第二幅度相等。
在另一实施例中，第一幅度小于第二幅度。这可以是有益的，因为它可以减少在两点Dixon磁共振数据的采集期间生成的声学噪声量。
在另一实施例中，第一持续时间乘以第一幅度乘以X小于或等于第二持续时间乘以第二幅度，其中，X是具有以下值中的任一个的数值：1.5、2、2.5、3、4和5。
在另一实施例中，所述第一读出梯度旁瓣具有第一幅度和第一持续时间。所述第二读出梯度旁瓣具有第二幅度和第二持续时间。第一持续时间乘以第一幅度乘以2等于第二持续时间乘以第二幅度。该实施例的优点之前已经被讨论过。
在另一实施例中，所述脉冲序列命令还被配置为对所述第二回波不对称地进行采样。这再次可以具有减少声学噪声量的益处。
在另一实施例中，所述第二分辨率在读出方向和/或相位编码方向上与所述第一分辨率的两倍一样高。该实施例可以是有益的，因为它可以减少在磁共振数据的采集期间生成的声学噪声量。
在另一实施例中，所述脉冲序列命令还被配置为使用所述双极性双回波脉冲序列从所述中心k空间区域和所述扩展的k空间区域对所述单点Dixon磁共振数据进行采样。在该实施例中，双极性双回波脉冲序列用来采集所有磁共振数据。再次，这可以具有减少所生成的声学噪声量的益处。
在另一实施例中，所述脉冲序列命令被配置为从所述双极性双回波脉冲序列的所述第二回波获得所述单点Dixon磁共振数据。例如，从第一回波和第二回波两者采集的数据可以用于两点Dixon磁共振协议，并且仅从第二回波采集的数据可以用作用于单点Dixon磁共振成像协议的输入。
在另一实施例中，所述脉冲序列命令被配置为使用单极性单回波脉冲序列从所述扩展的k空间区域对所述单点Dixon磁共振数据部分地进行采样。所述脉冲序列命令还被配置为使用所述双极性双回波脉冲序列从所述中心k空间区域对所述单点Dixon磁共振数据进行采样，并且对所述扩展的k空间区域部分地进行采样。在该范例中，单极性单回波脉冲序列与双回波脉冲序列的组合用来采集磁共振数据。这可以提供比单独使用双极性双回波脉冲序列更完整的对k空间的采样。这可以提供更好的图像质量。
在另一实施例中，所述第一读出梯度旁瓣具有第一幅度和第一持续时间。所述第二读出梯度旁瓣具有第二幅度和第二持续时间。所述第一持续时间乘以所述第一幅度小于所述第二持续时间乘以所述第二幅度。该实施例具有使第一分辨率小于第二分辨率的效果。这例如可以是有益的，因为它可以减少当磁共振数据被采集时生成的噪声量。
在另一实施例中，第一幅度和第二幅度是相同的，这可以是有益的，因为它可以实现梯度幅度的最小化，这可以减少声学噪声的生成。
在另一实施例中，例如，第一幅度和第二幅度可以是相同的。第一幅度的值然后能够被选择为减少噪声。
在另一实施例中，第一幅度小于第二幅度。这可以是有益的，因为它可以减少在对两点Dixon磁共振数据的采集期间生成的声学噪声量。
在另一实施例中，第一持续时间乘以第一幅度乘以X小于或等于第二持续时间乘以第二幅度。其中，X是具有以下值中的任一个的数值：1.5、2、2.5、3、4和5。
在另一实施例中，所述第一读出梯度旁瓣具有第一幅度和第一持续时间。所述第二读出梯度旁瓣具有第二幅度和第二持续时间。2乘以第一持续时间乘以第一幅度等于第二持续时间x第二幅度。该实施例的优点之前已经被讨论过。
在另一实施例中，所述单极性单回波脉冲序列具有第三读出梯度旁瓣，所述第三读出梯度旁瓣具有第三幅度和第三持续时间。第三持续时间乘以第三幅度大于第二持续时间乘以第二幅度。在这种情况下，单极性单回波脉冲序列用来对还未通过双极性双回波脉冲序列采样的k空间区域进行采样。当双极性双回波脉冲序列被使用时，扩展的k空间区域的邻近中心k空间区域的部分不被采样。该实施例可以用来对完全在中心k空间区域之外的k空间的线进行采样。对于单点Dixon协议，它提供了以相同的分辨率对k空间进行采样。
在另一实施例中，第三幅度等于或小于第二幅度。
在另一实施例中，所述单极性单回波脉冲序列具有第三读出梯度旁瓣，所述第三读出梯度旁瓣具有第三幅度和第三持续时间。第三持续时间乘以第三幅度等于第二持续时间乘以第二幅度乘以1.5倍。在这种情况下，单极性单回波脉冲序列用来对还未通过双极性双回波脉冲序列采样的k空间区域进行采样。当双极性双回波脉冲序列被使用时，扩展的k空间区域的邻近中心k空间区域的部分不被采样。该实施例可以用来对完全在中心k空间区域之外的k空间的线进行采样。对于单点Dixon协议，它提供了以相同的分辨率对k空间进行采样。
在另一实施例中，第三幅度和第二幅度是相同的。这可以是有益的，因为梯度脉冲的幅度可以被最小化，这可以减少声学噪声生成。
在另一实施例中，所述双极性双回波脉冲序列和所述单极性单回波脉冲序列具有完全相同的重复时间和翻转角。该实施例可以是有益的，因为对完全相同的重复时间和翻转角的使用意味着从两个脉冲序列采集的数据能够被组合，并且图像的对比度和其他质量将是相同的。
在另一实施例中，所述脉冲序列命令被配置为采集在时间上交错的所述两点Dixon磁共振成像数据和所述单点Dixon磁共振成像数据。例如，双极性双回波脉冲序列可以用来仅对中心k空间区域进行采样，并且单极性单回波脉冲序列可以用来对包括扩展的k空间区域和中心k空间区域两者的完整k空间进行采样。这两个序列然后可以被交错。在这种情况下，所述双极性双回波脉冲序列和所述单极性单回波脉冲序列具有完全相同的重复时间和翻转角，然而其不是严格必需的。
在另一实施例中，对所述第二分辨率水图像和所述第二分辨率脂肪图像的所述计算至少部分地使用所述两点Dixon磁共振成像数据来计算。以此方式，两点Dixon磁共振成像数据不仅用来计算磁场不均匀性图，而且用来改善例如得到的水图像和脂肪图像中的信噪比。
在另一方面中，本发明提供了一种操作磁共振成像系统的方法。所述方法包括利用脉冲序列命令控制所述磁共振成像系统采集两点Dixon磁共振数据和单点Dixon磁共振数据。所述脉冲序列命令被配置为根据两点Dixon磁共振成像协议从感兴趣区域采集两点Dixon磁共振数据。所述脉冲序列命令被配置为根据单点Dixon磁共振成像协议从所述感兴趣区域采集所述单点Dixon磁共振数据。所述两点Dixon磁共振成像协议被配置为从中心k空间区域采样所述两点Dixon磁共振数据。
所述单点Dixon磁共振成像协议被配置为从扩展的k空间区域对所述单点Dixon磁共振数据进行采样。所述扩展的k空间区域大于所述中心k空间区域。所述扩展的k空间区域至少部分地围绕所述中心k空间区域。所述两点Dixon磁共振成像协议被配置用于以第一分辨率生成所述感兴趣区域的图像。所述单点Dixon磁共振成像协议被配置用于以第二分辨率生成所述感兴趣区域的图像。所述第二分辨率高于所述第一分辨率。所述方法还包括使用所述两点Dixon磁共振数据来计算第一磁场不均匀性图。所述方法还包括通过将所述第一分辨率磁场不均匀性图内插到所述第二分辨率来计算第二磁场不均匀性图。所述方法还包括使用所述单点Dixon磁共振成像数据和所述第二分辨率磁场不均匀性图来计算第二分辨率水图像和第二分辨率脂肪图像。该实施例的优点之前已经被讨论过。
在另一方面中，本发明提供了一种计算机程序产品，包括用于由控制磁共振成像系统的处理器执行的机器可执行指令。对所述机器可执行指令的执行使所述处理器利用脉冲序列命令控制所述磁共振成像系统采集两点Dixon磁共振数据和单点Dixon磁共振数据。所述脉冲序列命令被配置为根据两点Dixon磁共振成像协议从感兴趣区域采集两点Dixon磁共振数据。所述脉冲序列命令被配置为根据单点Dixon磁共振成像协议从所述感兴趣区域采集所述单点Dixon磁共振数据。
所述两点Dixon磁共振成像协议被配置为从中心k空间区域对所述两点Dixon磁共振数据进行采样。所述单点Dixon磁共振成像协议被配置为从扩展的k空间区域对所述单点Dixon磁共振数据进行采样。所述扩展的k空间区域大于所述中心k空间区域。所述扩展的k空间区域至少部分地围绕所述中心k空间区域。所述两点Dixon磁共振成像协议被配置用于以第一分辨率生成所述感兴趣区域的图像。所述单点Dixon磁共振成像协议被配置用于以第二分辨率生成所述感兴趣区域的图像。所述第二分辨率高于所述第一分辨率。
对所述机器可执行指令的执行还使所述处理器使用所述两点Dixon磁共振数据来计算第一分辨率磁场不均匀性图。对所述机器可执行指令的执行还使所述处理器通过将所述第一分辨率磁场不均匀性图内插到所述第二分辨率来计算第二分辨率磁场不均匀性图。对所述机器可执行指令的执行还使所述处理器使用所述单点Dixon磁共振数据和所述第二分辨率磁场不均匀性图来计算第二分辨率水图像和第二分辨率脂肪图像。该实施例的优点之前已经被讨论过。
应当理解，只要组合的实施例不相互排斥，就可以组合本发明的前述实施例中的一个或多个。
如本领域的技术人员将认识到的，本发明的各方面可以被体现为一种装置、方法或计算机程序产品。因此，本发明的各方面可以采用以下形式：完全硬件实施例、完全软件实施例(包括固件、驻留软件、微代码等)或在本文中全部被通称为“电路”、“模块”或“系统”的组合软件方面和硬件方面的实施例。此外，本发明的各方面可以采用被体现在一个或多个计算机可读介质中的计算机程序产品的形式，所述一个或多个计算机可读介质具有被体现在其上的计算机可执行代码。
可以利用一个或多个计算机可读介质的任何组合。计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。如本文中所使用的“计算机可读存储介质”包括可以存储能由计算设备的处理器执行的指令的任何有形存储介质。计算机可读存储介质可以被称为计算机可读非瞬态存储介质。计算机可读存储介质还可以被称为有形计算机可读介质。在一些实施例中，计算机可读存储介质还能够存储能够由计算设备的处理器访问的数据。计算机可读存储介质的范例包括但不限于：软盘、磁硬盘驱动器、固态硬盘、闪速存储器、USB拇指驱动器、随机存取存储器(RAM)、只读存储器(ROM)、光盘、磁光盘以及处理器的寄存器文件。光盘的范例包括压缩盘(CD)和数字多用盘(DVD)，例如，CD-ROM、CD-RW、CD-R、DVD-ROM、DVD-RW或DVD-R盘。术语计算机可读存储介质还指能够由计算机设备经由网络或通信链路进行存取的各种类型的记录媒介。例如，可以在调制解调器上、在互联网上或在局域网上检索数据。可以使用任何适当的介质来发送在计算机可读介质上体现的计算机可执行代码，所述任何适当的介质包括但不限于：无线、有线、光纤线缆、RF等或前项的任何合适的组合。
计算机可读信号介质可以包括例如在基带中或作为载波的部分的、具有体现在其中的计算机可执行代码的传播数据信号。这种传播信号可以采用各种形式中的任何形式，包括但不限于：电磁、光学或其任何合适的组合。计算机可读信号介质可以是任何计算机可读介质：所述计算机可读介质不是计算机可读存储介质并且能够传递、传播或传输用于由指令执行系统、装置或设备使用的程序或与指令执行系统、装置或设备结合使用的程序。
“计算机存储器”或“存储器”是计算机可读存储介质的范例。计算机存储器是能由处理器直接访问的任何存储器。“计算机存储设备”或“存储设备”是计算机可读存储介质的另外的范例。计算机存储设备可以是任何非易失性计算机可读存储介质。在一些实施例中，计算机存储设备也可以是计算机存储器或反之亦然。
如本文中所使用的“处理器”包括能够执行程序或机器可执行指令或计算机可执行代码的电子部件。对包括“处理器”的计算设备的引用应被解读为可能包含多于一个处理器或处理核。处理器可以例如是多核处理器。处理器也可以指代在单个计算机系统内的或被分布在多个计算机系统之间的处理器的集合。术语计算设备也应被解读为可能指代计算设备的集合或网络，所述计算设备中的每个计算设备均包括一个处理器或多个处理器。计算机可执行代码可以由多个处理器来执行，这些处理器可以处于同一计算设备内，或者甚至可以跨多个计算设备而分布。
计算机可执行代码可以包括使处理器执行本发明的一方面的机器可执行指令或程序。用于执行针对本发明的各方面的操作的计算机可执行代码可以以一种或多种编程语言的任何组合来编写，所述一种或多种编程语言包括面向对象的编程语言(例如，Java、Smalltalk、C++等)和常规的程序编程语言(例如，“C”编程语言或类似的编程语言)，并且被编译成机器可执行指令。在一些实例中，计算机可执行代码可以是高级语言的形式或是预编译的形式，并且可以与解释器联合使用，所述解释器在运行中生成机器可执行指令。
计算机可执行代码可以完全在用户的计算机上、部分在用户的计算机上、作为独立软件包、部分在用户的计算机上并且部分在远程计算机上，或者完全在远程计算机或服务器上运行。在最后一种情况中，远程计算机可以通过任何类型的网络连接到用户的计算机，所述任何类型的网络包括局域网(LAN)或广域网(WAN)，或者可以(例如通过使用互联网服务提供商的互联网)对外部计算机进行连接。
参考根据本发明的实施例的方法、装置(系统)以及计算机程序产品的流程图图示和/或框图描述了本发明的各方面。应当理解，在适用时能够由计算机可执行代码形式的计算机程序指令来实施流程图、图示和/或框图的每个框或框的部分。还应当理解，当互不排斥时，可以对不同的流程图、图示和/或框图中的框的组合进行组合。可以将这些计算机程序指令提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器以产生机器，使得经由计算机或其他可编程数据处理装置的处理器运行的指令创建用于实施流程图和/或一个或多个框图框中指定的功能/动作的单元。
这些计算机程序指令也可以被存储在计算机可读介质中，所述计算机可读介质能够指导计算机、其他可编程数据处理装置或其他设备来以特定方式起作用，使得被存储在计算机可读介质中的指令产生一款制品，所述制品实施在流程图和/或一个或多个框图框中指定的功能/动作的指令。
计算机程序指令也可以被加载到计算机、其他可编程数据处理装置或其他设备上，以使一系列操作步骤被执行在计算机、其他可编程装置或其他设备上，从而产生计算机实施的过程，使得在计算机或其他可编程装置上运行的指令提供用于实施在流程图和/或一个或多个框图框中指定的功能/动作的过程。
如本文中所使用的“用户接口”是允许用户或操作者与计算机或计算机系统交互的接口。“用户接口”还可以被称为“人机接口设备”。用户接口可以向操作者提供信息或数据并且/或者从操作者接收信息或数据。用户接口可以使得来自操作者的输入能够被计算机接收，并且可以从计算机向用户提供输出。换言之，用户接口可以允许操作者控制或操纵计算机，并且接口可以允许计算机指示操作者的控制或操纵的效果。显示器或图形用户接口上对数据或信息的显示是向操作者提供信息的范例。通过键盘、鼠标、跟踪球、触控板、指点杆、图形输入板、操纵杆、游戏手柄、网络摄像头、头戴式设备、脚踏板、有线手套、遥控器以及加速度计来接收数据是使得能够从操作者接收信息或数据的用户接口部件的全部范例。
如本文中所使用的“硬件接口”涵盖使得计算机系统的处理器能够与外部计算设备和/或装置交互并且/或者控制外部计算设备和/或装置的接口。硬件接口可以允许处理器向外部计算设备和/或装置发送控制信号或指令。硬件接口还可以使得处理器能够与外部计算设备和/或装置交换数据。硬件接口的范例包括但不限于：通用串行总线、IEEE1394端口、并行端口、IEEE1284端口、串行端口、RS-232端口、IEEE-488端口、蓝牙连接、无线局域网连接、TCP/IP连接、以太网连接、控制电压接口、MIDI接口、模拟输入接口以及数字输入接口。
如本文中所使用的“显示器”或“显示设备”涵盖适于显示图像或数据的输出设备或用户接口。显示器可以输出视觉、听觉和或触觉数据。显示器的范例包括，但不限于：计算机监视器、电视屏幕、触摸屏、触觉电子显示器、盲文屏幕、阴极射线管(CRT)、存储管、双稳态显示器、电子纸、矢量显示器、平板显示器、真空荧光显示器(VF)、发光二极管(LED)显示器、电致发光显示器(ELD)、等离子显示面板(PDP)、液晶显示器(LCD)、有机发光二极管显示器(OLED)、投影仪以及头戴式显示器。
磁共振(MR)数据在本文中被定义为是在磁共振成像扫描期间由磁共振装置的天线所记录的对通过原子自旋发射的射频信号的测量结果。磁共振成像(MRI)图像或MR图像在本文中被定义为是对在磁共振成像数据内包含的解剖数据所重建的二维可视化或三维可视化。这种可视化能够使用计算机来执行。
附图说明
在下文中，将仅通过范例的方式并且参考附图来描述本发明的优选实施例，在附图中：
图1图示了磁共振成像系统的范例；
图2示出了图示操作权利要求1的磁共振成像系统的方法的流程图；
图3图示了脉冲序列的范例；
图4示出了k空间中的采样模式的范例；
图5图示了脉冲序列的另一范例；
图6示出了k空间中的采样模式的另一范例；
图7图示了脉冲序列的另一范例；
图8示出了k空间中的采样模式的另一范例；
图9图示了脉冲序列的另一范例；
图10图示了脉冲序列的另一范例；
图11示出了k空间中的采样模式的另一范例；以及
图12示出了k空间中的采样模式的另一范例。
附图标记列表
100 磁共振成像系统
104 磁体
106 磁体的膛
108 成像区
109 感兴趣区域
110 磁场梯度线圈
112 磁场梯度线圈电源
114 射频线圈
116 收发器
118 对象
120 对象支撑物
126 计算机系统
128 硬件接口
130 处理器
132 用户接口
134 计算机存储器
140 机器可执行指令
142 脉冲序列命令
144 两点Dixon磁共振数据
146 单点Dixon磁共振数据
148 第一分辨率磁场不均匀性图
150 第一分辨率水图像
152 第一分辨率脂肪图像
154 第二分辨率磁场不均匀性图
156 第二分辨率水图像
158 第二分辨率脂肪图像
200 利用脉冲序列命令控制磁共振成像系统采集两点Dixon磁共振数据和单点Dixon磁共振数据
202 使用两点Dixon磁共振数据来计算第一分辨率磁场不均匀性图
204 通过将第一分辨率磁性不均匀性图内插到第二分辨率来计算第二分辨率磁场不均匀性图
206 使用单点Dixon磁共振成像数据和第二分辨率磁场不均匀性图；奥计算第二分辨率水图像和第二分辨率脂肪图像
300 双极性双回波脉冲序列
302 读出梯度
304 相位编码梯度
306 切片选择梯度
308 射频发射和接收
310 第一读出梯度旁瓣
312 第一幅度
314 第一持续时间
316 第二读出梯度旁瓣
318 第二幅度
320 第二持续时间
322 RF脉冲
324 第一采样
326 第二采样
400 采样的k空间
402 中心k空间区域
404 扩展的k空间区域
406 在324期间采样的磁共振数据
408 在326期间采样的磁共振数据
500 单极性单回波脉冲序列
502 第三读出梯度旁瓣
504 第三幅度
506 第三持续时间
508 第三采样
600 在508期间采样的磁共振数据
700 双极性双回波脉冲序列
900 单极性单回波脉冲序列
1000 双极性双回波脉冲序列
1100 在324期间采样的磁共振数据
1102 在326期间采样的磁共振数据
1104 未采样的区域
1200 在508期间采样的磁共振数据
1202 在修改的脉冲序列900被应用的情况下的k空间区域
1204 在脉冲序列1000被应用的情况下的k空间区域
具体实施方式
在这些附图中，类似编号的元件是等效的元件或者执行相同的功能。如果功能等效，则在后面的附图中将不必讨论先前已经讨论过的元件。
图1示出了具有磁体104的磁共振成像系统100的范例。磁体104是超导圆柱型磁体，其具有穿过其的膛106。也可以使用不同类型的磁体；例如，也可以使用分裂式圆柱型磁体和所谓的开放式磁体。分裂式圆柱型磁体类似于标准的圆柱型磁体，除了低温恒温器已经被分成两个部分以允许进入磁体的等平面，这种磁体可以例如与带电粒子射束治疗结合使用。开放式磁体具有两个磁体部分，一个在另一个之上，其间具有足够大的空间以容纳对象：两个部分的区的布置类似于亥姆霍兹线圈的布置。开放式磁体很受欢迎，因为对象受限制较少。在圆柱型磁体的低温恒温器内部有一组超导线圈。在圆柱型磁体104的膛106内有成像区108，在该成像区中磁场足够强且均匀以执行磁共振成像。示出了成像区108内的感兴趣区域109。磁共振数据通常针对感兴趣区域来进行采集。对象118被示为由对象支撑物120支撑，使得对象118的至少部分处于成像区108和感兴趣区域109内。
在磁体的膛106内还有一组磁场梯度线圈110，该组磁场梯度线圈用于采集磁共振数据以对磁体104的成像区108内的磁自旋进行空间编码。磁场梯度线圈110被连接到磁场梯度线圈电源112。磁场梯度线圈110旨在是代表性的。通常，磁场梯度线圈110包含三组独立的线圈，这三组独立的线圈用于在三个正交空间方向上进行空间编码。磁场梯度电源向磁场梯度线圈提供电流。提供给磁场梯度线圈110的电流作为时间的函数而被控制并且可以是斜变的或脉冲的。
与成像区108相邻的是射频线圈114，射频线圈用于操纵成像区108内的磁自旋的取向并且用于接收也来自成像区108内的自旋的射频发射。射频天线可以包含多个线圈元件。射频天线也可以被称为通道或天线。射频线圈114被连接到射频收发器116。射频线圈114和射频收发器116可以由单独的发射线圈和接收线圈以及单独的发射器和接收器来替换。应当理解，射频线圈114和射频收发器116是代表性的。射频线圈114还旨在表示专用发射天线和专用接收天线。同样地，收发器116还可以表示单独的发射器和单独的接收器。射频线圈114还可以具有多个接收/发射元件，并且射频收发器116可以具有多个接收/发射通道。例如，如果执行诸如SENSE的并行成像技术，则射频114可以具有多个线圈元件。
收发器116和梯度控制器112被示为被连接到计算机系统126的硬件接口128。该计算机系统还包括与硬件系统128、存储器134和用户接口132通信的处理器130。存储器134可以是处理器130可访问的存储器的任何组合。这可以包括例如主存储器、高速缓存存储器以及诸如闪速RAM、硬盘驱动器的非易失性存储器或其他存储设备。在一些范例中，存储器134可以被认为是非瞬态计算机可读介质。
存储器134被示为包含机器可执行指令140。机器可执行指令140使得处理器130能够控制磁共振成像系统100的操作和功能。机器可执行指令140还可以使得处理器130能够执行各种数据分析和计算功能。计算机存储器134还被示为包含脉冲序列命令142。脉冲序列命令被配置用于控制磁共振成像系统100从对象采集磁共振数据。
计算机存储器134还示为包含当脉冲序列命令142被执行时从感兴趣区域109采集的两点Dixon磁共振数据144和单点Dixon磁共振数据146。存储器134还被示为包含从两点Dixon磁共振数据144计算的第一分辨率磁场不均匀性图148。存储器134还被示为可选地包含同样可以从第一分辨率磁场不均匀性图148计算的第一分辨率水图像150和第一分辨率脂肪图像152。
存储器134还被示为包含通过将第一分辨率磁场不均匀性图148内插到第二分辨率而创建或计算的第二分辨率磁场不均匀性图154。存储器134还被示为包含使用单点Dixon磁共振数据146和第二分辨率磁场不均匀性图154计算的第二分辨率水图像156和第二分辨率脂肪图像158。如果第一分辨率水图像150和第一分辨率脂肪图像152可获得，则它们也可以被使用。例如，第一分辨率水图像和第一分辨率脂肪图像也可以被内插到第二分辨率，并且用于一致性检查。
在一些情况下，在单点Dixon磁共振数据与两点Dixon磁共振数据144之间存在一定的交叠或冗余。例如，两点Dixon磁共振数据144中的一些可以用作单点Dixon磁共振数据146。
图2示出了图示操作图1的磁共振成像系统100的方法的流程图。首先，在步骤200中，利用脉冲序列命令142控制磁共振成像系统100采集两点Dixon磁共振数据144和单点Dixon磁共振数据146。两点Dixon磁共振数据144以第一分辨率被采集。单点Dixon磁共振数据146以第二分辨率被采集。第二分辨率高于第一分辨率。在步骤200之后，步骤202被执行。在步骤202中，第一分辨率磁场不均匀性图148使用两点Dixon磁共振数据144来计算。
接下来在步骤204中，第二分辨率磁场不均匀性图154通过将第一分辨率磁场不均匀性图148内插到第二分辨率来计算。最后在步骤206中，第二分辨率水图像156和第二分辨率脂肪图像158使用单点Dixon磁共振成像数据146和第二分辨率磁场不均匀性图154来计算。图2中的方法可以具有在对磁共振数据的采集期间减少声学噪声的效果。单点Dixon磁共振数据可以被采集，使得更少的声学噪声出现。然而，需要磁场不均匀性图。两点Dixon磁共振数据144以较低的第一分辨率被采集。这具有减少声学噪声量的效果。对两点Dixon磁共振数据144与单点Dixon磁共振数据146的采集的组合使得准确的水图像和脂肪图像能够在减少的声学噪声量的情况下被获得。
相比于标准梯度回波成像，梯度回波Dixon成像同时提供了优良的脂肪压缩和多个对比度。其通常利用双回波序列来实施，所述双回波序列在利用双极性读出梯度的每次激励之后生成两个回波。图3示出了覆盖三个重复时间的3D双梯度回波Dixon序列的示意性(简化)脉冲序列图。它描绘了双极性双回波脉冲序列300的表示，并未示出其所有细节。时间线图示了读出梯度302、相位编码梯度304、切片选择梯度306、以及射频发射和接收308。
读出梯度302具有第一读出梯度旁瓣310，第一读出梯度旁瓣具有第一幅度312和第一持续时间314。其还具有第二读出梯度旁瓣316，第二读出梯度旁瓣具有第二幅度318和第二持续时间320。第二幅度具有与第一幅度相反的极性。还存在RF脉冲322以及然后对应于第一读出梯度旁瓣310的第一采样324和对应于第二读出梯度旁瓣316的第二采样326。图3中描绘的脉冲序列300通常将用来针对两点Dixon磁共振成像协议采集磁共振数据。脉冲序列300可以例如被调整为以第二分辨率对两点Dixon磁共振数据进行采样。在下面的图4中图示了用于图3中的脉冲序列300的k空间采样模式。
图4示出了用来表示要通过图3的脉冲序列300采样的k空间的大正方形400。在采样的k空间400的中心内存在更小的正方形402。正方形402表示中心k空间区域。在正方形402与400之间是扩展的k空间区域404。扩展的k空间区域404至少部分地围绕中心k空间区域402。线406指示k空间中的第一轨迹，并且线408指示k空间中的第二轨迹。线406对应于在第一采样324期间采集的数据，并且线408对应于在图3的第二采样326期间采集的数据。在各种重复期间，轨迹406和408被垂直地移动以对整个正方形400进行采样。在该范例中，双极性双回波脉冲序列300用来对整个正方形进行采样。对这种采样模式的使用可能是不利的，因为大的声学噪声量可以被生成。
备选地，单回波序列可以被使用，所述单回波序列在每次激励之后仅采集一个回波，如下面的图5中示出的。它允许以脂肪压缩和信噪比的一致性为代价来达到更短的重复时间。下面的图5示出了覆盖三个重复时间的3D单梯度回波Dixon序列的示意性(简化)脉冲序列。它描绘了单极性单回波脉冲序列500的表示，该单极性单回波脉冲序列类似于脉冲序列300，除了在该实例中仅存在单个读出梯度旁瓣502。该单个读出梯度旁瓣在本文中被称为第三读出梯度旁瓣。术语第三读出梯度旁瓣用作将其与双极性双回波脉冲序列的第一读出梯度旁瓣和第二读出梯度旁瓣区别开的标签。第三读出梯度旁瓣502具有第三幅度504和第三持续时间506。再次存在RF脉冲322，并且存在对应于第三读出梯度旁瓣502的第三采样508。
图6再次示出了这次用来表示要通过图5的脉冲序列500采样的k空间的大正方形400。线600指示k空间中的轨迹，其对应于在采样508期间采集的数据。图5中的单极性单回波脉冲序列500可以用来对整个正方形进行采样。
梯度回波Dixon成像通常由于以下项与高声压级相关联：
例如在屏气的情况下在腹部成像中的扫描时间约束，
空间分辨率要求，
回波时间限制，以实现鲁棒的且在信噪比方面高效的水脂肪分离，以及
接收器带宽约束，以限制由于化学漂移和主场不均匀性的失真。
双回波序列需要读出梯度的额外反转，这促进了声学噪声。由于在估计主场可靠性的空间变化方面的根本困难，单回波序列使准确的水脂肪分离复杂化。
范例可以将低分辨率双回波与高分辨率单回波序列合并来减少梯度回波Dixon成像中的声学噪声，同时满足提到的约束。
在下面更详细地描述了两个范例。
下面的图7和图9中图示的第一实施例将具有降低的空间分辨率的双回波序列与具有全空间分辨率的单回波序列进行组合。空间分辨率在用于双回波序列的读出和相位编码方向上被减半，以降低梯度幅度并且因此降低声压级。由于相同的原因，并且为了保留对比度，对于单回波序列，重复时间和翻转角被保持相同。对于双回波序列，水信号和脂肪信号在回波时间处接近反相和同相，并且对于单回波序列，在回波时间处接近正交。以这种方式，水脂肪分离中的有利噪声传播被实现。
图7示出了修改的双极性双回波脉冲序列700。脉冲序列700类似于图3中描绘的脉冲序列300。然而，第一幅度312已经通过2的因子被降低，并且第二幅度318也已经通过2的因子被降低。用于相位编码304的另外的梯度旁瓣的幅度也已经通过2的因子被降低。这具有在读出和相位编码方向上通过2的因子降低分辨率的效果。这也具有减少当脉冲序列700被执行时生成的声学噪声量的效果。
图8再次示出了采样的k空间400。在该实例中，线406指示对应于在图7中的采样324期间采集的数据的k空间中的轨迹，并且线408指示对应于在图7中的采样326期间采集的数据的k空间中的轨迹。能够看出，仅在中心k空间区域402内对k空间进行采样。
减少声学噪声的一种方式将是首先使用脉冲序列700对中心k空间区域402进行采样。第一磁场不均匀性图然后可以被计算。单点Dixon磁共振数据然后可以使用脉冲序列500和与图6中的k空间中的轨迹类似的k空间中的轨迹来采集。当计算第二分辨率水图像和第二分辨率脂肪图像时，第一磁场不均匀性图可以如之前讨论的那样被使用。图5的脉冲序列500和图7的脉冲序列700也可以被交错。图5中的脉冲序列500和图7中的脉冲序列700可以被调整，使得它们具有完全相同的重复时间和翻转角。
图9示出了单极性单回波脉冲序列900的另一范例。脉冲序列900是图5的脉冲序列500的修改。它现在具有与图7中的脉冲序列700相同的重复时间。图7的脉冲序列700和图9的脉冲序列900然后能够被容易地交错，其中脉冲序列700沿着如图9中图示的k空间406和408中的轨迹进行采样，并且脉冲序列900沿着如图6中描绘的k空间600中的轨迹进行采样。
用于图7和图9的混合双/单回波序列的扫描时间比用于图3的双回波序列的扫描时间长大约50％。然而，它可以通过也在切片编码方向上或进一步在相位编码方向上降低空间分辨率来降低。此外，中心k空间的重复的采集可以被用来改善信噪比。备选地，任何重复的采集可以通过仅在相位编码或切片选择方向上降低空间分辨率来避免。
如果双回波序列和单回波序列在时间上交错，那么重复时间和翻转角也可以在这两个序列之间被改变以进一步以对比度为代价来降低声压级。
下面的图10中图示的第二实施例将低分辨率全回波采集与高分辨率部分回波采集进行组合。空间分辨率在用于第一回波的读出方向上被减半，以降低梯度幅度或回波时间。由于相同的原因，仅对第二回波的四分之三进行采样。
一般来说，低分辨率或部分回波采集促进梯度回波Dixon成像中的有利回波时间，尤其是当较低的读出梯度强度被用来减少声学噪声时。
图10图示了双极性双回波脉冲序列1000的另一范例。图10的脉冲序列1000类似于图7的脉冲序列700，除了第二持续时间320现在长于第一持续时间314。这具有对k空间不对称地进行采样的效果。相同的效果也可以通过相对于彼此调整幅度312和318或调整持续时间和幅度的组合来实现。然而，为了减少声学噪声量，可以有益的是降低两者的幅度并且仅仅使一个持续时间长于另一个。
图11图示了对k空间中的轨迹的影响。再次示出了完整的采样的k空间400。被标记1100的线对应于在图10的采样324期间采集的数据，并且被标记1102的线对应于在图10的采样326期间采集的数据。沿着轨迹1100采样的数据和沿着轨迹1102采样的数据的一部分能够是两点Dixon磁共振数据。例如，在中心k空间区域402外部沿着轨迹1102采样的数据能够被忽视。沿着轨迹1102采集到的所有数据能够是单点Dixon磁共振数据。对图10中的脉冲序列1000的使用然后将实现对两点Dixon磁共振数据和单点Dixon磁共振数据两者的同时采集。在该实例中，右侧的k空间区域1104保持不被采样。在图11中，中心k空间区域402已经被垂直地扩展，使得所采集的所有两点Dixon磁共振数据用于对第一分辨率磁场不均匀性图的计算。在其他范例中，中心k空间区域402可以是与图8中示出的尺寸相同的尺寸。
在图11中，因为“第一读出梯度旁瓣”可以依据实施方式而在“第二读出梯度旁瓣”之前或之后被执行，线1100和1102的读出的顺序可以被反转。箭头1100和1102的方向可以因此在一些范例中被反转。
图9的脉冲序列900能够被修改，使得它具有与图10中的脉冲序列1000相同的重复时间。图12图示了图10的脉冲序列1000和图9的修改的脉冲序列900如何能够被组合。被标记1202的k空间区域是使用修改的脉冲序列900进行采集的区域。被标记1204的区域是使用脉冲序列1000进行采集的k空间区域。能够看出，图12中的未采样的区域1104小于图11中的未采样的区域1104。为了如图12中图示的那样对k空间进行采样，图10的脉冲序列1000和图9的修改的脉冲序列900也可以被交错。
虽然已经在附图和前面的描述中详细说明和描述了本发明，但是这样的说明和描述应当被认为是说明性或示范性的而非限制性的；本发明不限于所公开的实施例。
通过研究附图、说明书以及权利要求书，本领域技术人员在实践请求保护的发明时能够理解并实现对所公开的实施例的其他变型。在权利要求中，“包括”一词不排除其他元件或步骤，并且词语“一”或“一个”不排除多个。单个处理器或其他单元可以实现在权利要求中记载的若干项的功能。虽然某些措施被记载在互不相同的从属权利要求中，但是这并不指示不能有利地使用这些措施的组合。计算机程序可以被存储/被分布在合适的介质上，例如，与其他硬件一起或作为其他硬件的部分供应的光学存储介质或固态介质，但是也可以以其他形式分布，例如，经由互联网或其他有线或无线的电信系统分布。权利要求中的任何附图标记都不应被解释为对范围的限制。

Claims (15)

1.一种磁共振成像系统(100)，包括：
存储器(134)，其用于存储机器可执行指令(140)和脉冲序列命令(142)，其中，所述脉冲序列命令被配置为根据两点Dixon磁共振成像协议(300、700、1000)从感兴趣区域采集两点Dixon磁共振数据，其中，所述脉冲序列命令被配置为根据单点Dixon磁共振成像协议(500、900、1000)从所述感兴趣区域采集单点Dixon磁共振数据(146)，其中，所述两点Dixon磁共振成像协议被配置为从中心k空间区域(402)对所述两点Dixon磁共振数据(144)进行采样(324、326)，其中，所述单点Dixon磁共振成像协议被配置为从扩展的k空间区域(404)和所述中心k空间区域对所述单点Dixon磁共振数据进行采样(326、508)，其中，所述扩展的k空间区域大于所述中心k空间区域，其中，所述扩展的k空间区域至少部分地围绕所述中心k空间区域，其中，所述两点Dixon磁共振成像协议被配置用于以第一分辨率采集所述感兴趣区域的图像，其中，所述单点Dixon磁共振成像协议被配置用于以第二分辨率采集所述感兴趣区域的图像，其中，所述第二分辨率高于所述第一分辨率；以及
处理器(130)，其用于控制所述磁共振成像系统，其中，对所述机器可执行指令的执行使所述处理器：
利用所述脉冲序列命令控制(200)所述磁共振成像系统采集所述两点Dixon磁共振数据和所述单点Dixon磁共振数据；
使用所述两点Dixon磁共振数据来计算(202)处于所述第一分辨率的第一分辨率磁场不均匀性图(148)；
通过将所述第一分辨率磁性不均匀性图内插到所述第二分辨率来计算(204)第二分辨率磁场不均匀性图(154)；并且
使用所述单点Dixon磁共振成像数据和所述第二分辨率磁场不均匀性图来计算(206)处于所述第二分辨率的第二分辨率水图像(156)和处于所述第二分辨率的第二分辨率脂肪图像(158)。
2.根据权利要求1所述的磁共振成像系统，其中，所述脉冲序列命令被配置为使用双极性双回波脉冲序列(300、700、1000)从中心k空间区域对所述两点Dixon磁共振数据进行采样，其中，所述第一回波是利用具有第一极性的第一读出梯度旁瓣(310)来生成的，并且所述第二回波是利用具有第二极性的第二读出梯度旁瓣(316)来生成的。
3.根据权利要求2所述的磁共振成像系统，其中，所述第一读出梯度旁瓣具有第一幅度(312)和第一持续时间(314)，其中，所述第二读出梯度旁瓣具有第二幅度(318)和第二持续时间(320)，并且其中，所述第一持续时间乘以所述第一幅度小于所述第二持续时间乘以所述第二幅度。
4.根据权利要求3所述的磁共振成像系统，其中，所述脉冲序列命令还被配置为对所述第二回波不对称地进行采样。
5.根据权利要求2、3或4所述的磁共振成像系统，其中，所述第二分辨率在读出方向和/或相位编码方向上与所述第一分辨率的两倍一样高。
6.根据权利要求2、3或4所述的磁共振成像系统，其中，所述脉冲序列命令还被配置为使用所述双极性双回波脉冲序列从所述中心k空间区域和所述扩展的k空间区域对所述单点Dixon磁共振数据进行采样。
7.根据权利要求6所述的磁共振成像系统，其中，所述脉冲序列命令被配置为从所述双极性双回波脉冲序列的所述第二回波(326)获得所述单点Dixon磁共振数据。
8.根据权利要求2所述的磁共振成像系统，其中，所述脉冲序列命令被配置为使用单极性单回波脉冲序列(500、900)从所述扩展的k空间区域对所述单点Dixon磁共振数据部分地进行采样，其中，所述脉冲序列命令还被配置为使用所述双极性双回波脉冲序列从所述中心k空间区域对所述单点Dixon磁共振数据进行采样，并且从所述扩展的k空间区域对所述单点Dixon磁共振数据部分地进行采样。
9.根据权利要求8所述的磁共振成像系统，其中，所述第一读出梯度旁瓣(310)具有第一幅度(312)和第一持续时间(314)，其中，所述第二读出梯度旁瓣(316)具有第二幅度(318)和第二持续时间(320)，其中，所述第一持续时间乘以所述第一幅度小于所述第二持续时间乘以所述第二幅度。
10.根据权利要求9所述的磁共振成像系统，其中，所述单极性单回波脉冲序列具有第三读出梯度旁瓣(502)，其中，所述第三读出梯度旁瓣具有第三幅度(504)和第三持续时间(506)，其中，所述第三持续时间乘以所述第三幅度大于所述第二持续时间乘以所述第二幅度。
11.根据权利要求8、9或10所述的磁共振成像系统，其中，所述双极性双回波脉冲序列和所述单极性单回波脉冲序列具有完全相同的重复时间和翻转角。
12.根据权利要求1、2或3所述的磁共振成像系统，其中，所述脉冲序列命令被配置为采集在时间上交错的所述两点Dixon磁共振成像数据和所述单点Dixon磁共振数据。
13.根据前述权利要求中的任一项所述的磁共振成像系统，其中，对所述第二分辨率水图像和所述第二分辨率脂肪图像的所述计算是至少部分地使用所述两点Dixon磁共振成像数据来计算的。
14.一种操作磁共振成像系统(100)的方法，其中，所述方法包括：
利用脉冲序列命令(142)控制(200)所述磁共振成像系统采集两点Dixon磁共振数据(144)和单点Dixon磁共振数据(146)，其中，所述脉冲序列命令被配置为根据两点Dixon磁共振成像协议(300、700、1000)从感兴趣区域(109)采集所述两点Dixon磁共振数据，其中，所述脉冲序列命令被配置为根据单点Dixon磁共振成像协议(500、900、1000)从所述感兴趣区域采集所述单点Dixon磁共振数据，其中，所述两点Dixon磁共振成像协议被配置为从中心k空间区域(402)对所述两点Dixon磁共振数据进行采样，其中，所述单点Dixon磁共振成像协议被配置为从扩展的k空间区域(404)和所述中心k空间区域对所述单点Dixon磁共振数据进行采样，其中，所述扩展的k空间区域大于所述中心k空间区域，其中，所述扩展的k空间区域至少部分地围绕所述中心k空间区域，其中，所述两点Dixon磁共振成像协议被配置用于以第一分辨率采集所述感兴趣区域的图像，其中，所述单点Dixon磁共振成像协议被配置用于以第二分辨率采集所述感兴趣区域的图像，其中，所述第二分辨率高于所述第一分辨率；
使用所述两点Dixon磁共振数据来计算(202)处于所述第一分辨率的第一分辨率磁场不均匀性图(148)；
通过将所述第一分辨率磁性不均匀性图内插到所述第二分辨率来计算(204)第二分辨率磁场不均匀性图(154)；以及
使用所述单点Dixon磁共振成像数据和所述第二分辨率磁场不均匀性图来计算(206)处于所述第二分辨率的第二分辨率水图像(156)和处于所述第二分辨率的第二分辨率脂肪图像(158)。
15.一种计算机程序产品，包括用于由控制磁共振成像系统(100)的处理器(130)执行的机器可执行指令(140)，其中，对所述机器可执行指令的执行使所述处理器：
利用脉冲序列命令(142)控制(200)所述磁共振成像系统采集两点Dixon磁共振数据(144)和单点Dixon磁共振数据(146)，其中，所述脉冲序列命令被配置为根据两点Dixon磁共振成像协议从感兴趣区域采集所述两点Dixon磁共振数据，其中，所述脉冲序列命令被配置为根据单点Dixon磁共振成像协议从所述感兴趣区域采集所述单点Dixon磁共振数据，其中，所述两点Dixon磁共振成像协议被配置为从中心k空间区域(402)对所述两点Dixon磁共振数据进行采样，其中，所述单点Dixon磁共振成像协议被配置为从扩展的k空间区域(404)和所述中心k空间区域对所述单点Dixon磁共振数据进行采样，其中，所述扩展的k空间区域大于所述中心k空间区域，其中，所述扩展的k空间区域至少部分地围绕所述中心k空间区域，其中，所述两点Dixon磁共振成像协议被配置用于以第一分辨率采集所述感兴趣区域的图像，其中，所述单点Dixon磁共振成像协议被配置用于以第二分辨率采集所述感兴趣区域的图像，其中，所述第二分辨率高于所述第一分辨率；
使用所述两点Dixon磁共振数据来计算(202)处于所述第一分辨率的第一分辨率磁场不均匀性图(148)；
通过将所述第一分辨率磁性不均匀性图内插到所述第二分辨率来计算(204)第二分辨率磁场不均匀性图(154)；并且
使用所述单点Dixon磁共振成像数据和所述第二分辨率磁场不均匀性图来计算(206)处于所述第二分辨率的第二分辨率水图像(156)和处于所述第二分辨率的第二分辨率脂肪图像(158)。
CN201980006440.9A 2018-01-11 2019-01-09 双分辨率Dixon磁共振成像 Active CN111542762B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP18151141.1 2018-01-11
EP18151141.1A EP3511725A1 (en) 2018-01-11 2018-01-11 Dual resolution dixon magnetic resonance imaging
PCT/EP2019/050385 WO2019137932A1 (en) 2018-01-11 2019-01-09 Dual resolution dixon magnetic resonance imaging

Publications (2)

Publication Number Publication Date
CN111542762A true CN111542762A (zh) 2020-08-14
CN111542762B CN111542762B (zh) 2024-03-15

Family

ID=60954928

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201980006440.9A Active CN111542762B (zh) 2018-01-11 2019-01-09 双分辨率Dixon磁共振成像

Country Status (5)

Country Link
US (1) US11280865B2 (zh)
EP (2) EP3511725A1 (zh)
JP (1) JP7337808B2 (zh)
CN (1) CN111542762B (zh)
WO (1) WO2019137932A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4043902A1 (en) * 2021-02-11 2022-08-17 Koninklijke Philips N.V. Dixon-type water/fat separation mr imaging
JP7221371B1 (ja) 2021-12-21 2023-02-13 ジーイー・プレシジョン・ヘルスケア・エルエルシー 磁気共鳴イメージング装置、血管画像生成方法、および記録媒体
US11927657B2 (en) * 2022-04-13 2024-03-12 Regents Of The University Of Michigan Multiresolution magnetic resonance fingerprinting systems and methods

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4871967A (en) * 1986-11-30 1989-10-03 Elscint Ltd. Spectral component separation in mr imaging
US5771893A (en) * 1994-09-06 1998-06-30 Kabushiki Kaisha Toshiba Method and apparatus for nuclear magnetic resonance imaging of physiological function information
US20080012565A1 (en) * 2006-04-25 2008-01-17 Mckenzie Charles A Self-calibration methods for parallel imaging and multipoint water-fat separation methods
US20080278162A1 (en) * 2007-05-02 2008-11-13 Bookwalter Candice A Cartesian continuous sampling with unequal gradients
WO2011080693A1 (en) * 2009-12-30 2011-07-07 Koninklijke Philips Electronics N.V. Water -fat separation in mri using partial echoes at different echo times
US20120268121A1 (en) * 2011-04-19 2012-10-25 Diego Hernando Method for r*2 quantification with magnetic resonance imaging
CN103300857A (zh) * 2012-03-12 2013-09-18 西门子公司 用于执行组合的磁共振－正电子发射断层造影的方法
US20160124064A1 (en) * 2013-06-06 2016-05-05 Koninklijke Philips N.V. Parallel mri with bo distortion correction and multi-echo dixon water-fat separation using regularised sense reconstruction
CN105659103A (zh) * 2013-08-30 2016-06-08 皇家飞利浦有限公司 狄克逊磁共振成像
CN105793721A (zh) * 2013-12-02 2016-07-20 皇家飞利浦有限公司 使用零回波时间脉冲序列的磁共振成像
US20160231405A1 (en) * 2013-09-16 2016-08-11 Koninklijke Philips N.V. Mr imaging with dixon-type water/fat separation
WO2016207035A1 (en) * 2015-06-26 2016-12-29 Koninklijke Philips N.V. Phase corrected dixon magnetic resonance imaging
US20170038446A1 (en) * 2014-04-24 2017-02-09 Junmin LIU Systems and methods for field mapping in magnetic resonance imaging
CN107072586A (zh) * 2015-08-06 2017-08-18 株式会社日立制作所 磁共振成像装置
CN107407714A (zh) * 2015-01-21 2017-11-28 皇家飞利浦有限公司 用于根据b0图和b1图来计算导出值的mri方法
CN107405502A (zh) * 2015-03-04 2017-11-28 皇家飞利浦有限公司 声辐射力成像

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002085376A (ja) * 2000-09-20 2002-03-26 Hitachi Medical Corp 核磁気共鳴イメージング装置および方法
WO2004097387A2 (en) * 2003-04-25 2004-11-11 Case Western Reserve University Dixon techniques in spiral trajectories with off-resonance correction
WO2005055136A2 (en) * 2003-11-26 2005-06-16 Board Of Regents, The University Of Texas System Method and apparatus for phase-sensitive magnetic resonance imaging

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4871967A (en) * 1986-11-30 1989-10-03 Elscint Ltd. Spectral component separation in mr imaging
US5771893A (en) * 1994-09-06 1998-06-30 Kabushiki Kaisha Toshiba Method and apparatus for nuclear magnetic resonance imaging of physiological function information
US20080012565A1 (en) * 2006-04-25 2008-01-17 Mckenzie Charles A Self-calibration methods for parallel imaging and multipoint water-fat separation methods
US20080278162A1 (en) * 2007-05-02 2008-11-13 Bookwalter Candice A Cartesian continuous sampling with unequal gradients
WO2011080693A1 (en) * 2009-12-30 2011-07-07 Koninklijke Philips Electronics N.V. Water -fat separation in mri using partial echoes at different echo times
US20120268121A1 (en) * 2011-04-19 2012-10-25 Diego Hernando Method for r*2 quantification with magnetic resonance imaging
CN103300857A (zh) * 2012-03-12 2013-09-18 西门子公司 用于执行组合的磁共振－正电子发射断层造影的方法
US20160124064A1 (en) * 2013-06-06 2016-05-05 Koninklijke Philips N.V. Parallel mri with bo distortion correction and multi-echo dixon water-fat separation using regularised sense reconstruction
CN105659103A (zh) * 2013-08-30 2016-06-08 皇家飞利浦有限公司 狄克逊磁共振成像
US20160231405A1 (en) * 2013-09-16 2016-08-11 Koninklijke Philips N.V. Mr imaging with dixon-type water/fat separation
CN105793721A (zh) * 2013-12-02 2016-07-20 皇家飞利浦有限公司 使用零回波时间脉冲序列的磁共振成像
US20170038446A1 (en) * 2014-04-24 2017-02-09 Junmin LIU Systems and methods for field mapping in magnetic resonance imaging
CN107407714A (zh) * 2015-01-21 2017-11-28 皇家飞利浦有限公司 用于根据b0图和b1图来计算导出值的mri方法
CN107405502A (zh) * 2015-03-04 2017-11-28 皇家飞利浦有限公司 声辐射力成像
WO2016207035A1 (en) * 2015-06-26 2016-12-29 Koninklijke Philips N.V. Phase corrected dixon magnetic resonance imaging
CN107072586A (zh) * 2015-08-06 2017-08-18 株式会社日立制作所 磁共振成像装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
H.W.PARK等: "Fast gradient-echo chemical-shift imaging", vol. 7, pages 340 - 345 *
HUANZHOU YU等: "413–422", vol. 55, pages 413 - 422 *
JINGFEI MA: "A single‐point dixon technique for fat‐suppressed fast 3D gradient‐echo imaging with a flexible echo time", vol. 27, pages 881 - 890, XP007919289, DOI: DOI 10.1002/jmri.21281 *
JINGFEI MA: "Dixon techniques for water and fat imaging", vol. 28, pages 543 - 558, XP002586398, DOI: 10.1002/JMRI.21492 *
WEITIAN CHEN等: "Fast Conjugate Phase Image Reconstruction Based on a Chebyshev Approximation to Correct for B0 Field Inhomogeneity and Concomitant Gradients", vol. 60, pages 1104 - 1111, XP055187444, DOI: 10.1002/mrm.21703 *

Also Published As

Publication number Publication date
US20200363486A1 (en) 2020-11-19
EP3692382A1 (en) 2020-08-12
WO2019137932A1 (en) 2019-07-18
US11280865B2 (en) 2022-03-22
EP3692382B1 (en) 2021-03-10
EP3511725A1 (en) 2019-07-17
JP7337808B2 (ja) 2023-09-04
JP2021510568A (ja) 2021-04-30
CN111542762B (zh) 2024-03-15

Similar Documents

Publication Publication Date Title
CN107923958B (zh) 相位校正的狄克逊磁共振成像
US11474182B2 (en) Motion detection in CEST magnetic resonance imaging based on Z-spectrum analysis
CN111542762B (zh) 双分辨率Dixon磁共振成像
WO2015197366A1 (en) Motion correction in magnetic resonance imaging
WO2018149737A1 (en) Phase error correction for bipolar read out gradients
JP2024056688A (ja) ３次元勾配インパルス応答関数の高次項の決定
CN111670370B (zh) 具有脂肪/水分离的mri
EP3019884A1 (en) Corrected magnetic resonance imaging using coil sensitivities
CN111263896B (zh) 在磁共振成像系统中对相位相关的伪影的数据驱动的校正
EP3699624A1 (en) Calculation of a b0 image using multiple diffusion weighted mr images
US11686800B2 (en) Correction of magnetic resonance images using simulated magnetic resonance images
EP4040179A1 (en) Adaptive water-fat shift in non-cartesian magnetic resonance imaging
EP4163657A1 (en) Motion correction using low resolution magnetic resonance images
JP7491670B2 (ja) ３次元勾配インパルス応答関数の高次項の決定
CN118103721A (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