CN111133327B - Dixon型水/脂肪分离MR成像 - Google Patents

Dixon型水/脂肪分离MR成像 Download PDF

Info

Publication number
CN111133327B
CN111133327B CN201880061999.7A CN201880061999A CN111133327B CN 111133327 B CN111133327 B CN 111133327B CN 201880061999 A CN201880061999 A CN 201880061999A CN 111133327 B CN111133327 B CN 111133327B
Authority
CN
China
Prior art keywords
space
water
echo
fat
echo signals
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
CN201880061999.7A
Other languages
English (en)
Other versions
CN111133327A (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 CN111133327A publication Critical patent/CN111133327A/zh
Application granted granted Critical
Publication of CN111133327B publication Critical patent/CN111133327B/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/4828Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • 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/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
    • 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

Abstract

本发明涉及一种Dixon型MR成像的方法。本发明的目的是提供一种使能高效和可靠的水/脂肪分离的方法。本发明的方法包括以下步骤:使对象(10)经受成像序列,其包括至少一个激励RF脉冲和切换的磁场梯度,其中,两个回波信号(第一回波信号和第二回波信号)在不同的回波时间(TE1、TE2)处生成;从所述对象(10)采集所述回波信号;根据所述回波信号重建水图像和/或脂肪图像,其中,来自水和脂肪的对回波信号的贡献使用k空间的第一区域中的两点Dixon技术和k空间的第二区域中的单点Dixon技术来分离,其中,所述第一区域与所述第二区域不同。换句话说,本发明提出应用于第一和第二回波信号两者的用于水/分离的两点Dixon技术与应用于两个回波信号之一(即,第一回波信号数据或第二回波信号数据)的单点Dixon技术之间的自适应切换,这取决于k空间中的位置。此外,本发明涉及MR设备(1)和要在MR设备(1)上运行的计算机程序。

Description

Dixon型水/脂肪分离MR成像
技术领域
本发明涉及磁共振(MR)成像的领域。其涉及一种置于MR设备的检查体积中的对象的MR成像的方法。本发明还涉及一种MR设备和一种要在MR设备上运行的计算机程序。
背景技术
利用磁场和核自旋之间的相互作用以便形成二维或三维图像的图像形成MR方法现在被广泛使用,特别是在医学诊断领域中,因为对于对软组织的成像,它们在许多方面优于其他成像方法,不需要电离辐射并且通常是无创的。
根据一般的MR方法,要检查患者的身体被布置在强、均匀的磁场B0中,磁场B0的方向同时限定测量基于的坐标系的轴(通常地,z轴)。磁场B0根据可以通过施加限定频率(所谓的拉莫尔频率,或者MR频率)的电磁交变场(RF场)所激励(自旋共振)的磁场强度产生针对个体核自旋的不同的能级。从宏观的角度来看,个体核自旋的分布产生整体磁化,所述整体磁化可以通过施加垂直于z轴的适当频率(RF脉冲)的电磁脉冲而偏离平衡状态,使得所述磁化关于z轴执行进动运动。进动运动描述了圆锥的表面,所述圆锥的孔径角被称为翻转角。翻转角的幅度取决于所施加的电磁脉冲的强度和持续时间。在所谓的90°脉冲的情况下,自旋从z轴偏转到横向平面(翻转角度90°)。
在RF脉冲的终止之后,磁化弛豫回到原始平衡状态,其中,z方向上的磁化以第一时间常数T1(自旋-晶格或纵向弛豫时间)再次建立,并且所述磁化强度在垂直于z方向的方向上以第二时间常数T2(自旋-自旋或横向弛豫时间)弛豫。可以借助于接收RF线圈来检测磁化的变化,所述接收RF线圈以这样的方式在MR设备的检查体积内被布置和取向:在垂直于z轴的方向上测量磁化的变化。在施加例如90°脉冲之后,横向磁化的衰变伴随着核自旋(由局部磁场不均匀性引起)从具有相同相位的有序状态到所有相位角均匀分布(失相)的状态的转变。失相可以借助于重聚焦脉冲(例如180°脉冲)来补偿。这在接收线圈中产生回波信号。
为了实现身体中的空间分辨率,沿着三个主轴延伸的恒定磁场梯度叠加在均匀磁场B0上,从而导致自旋共振频率的线性空间相关性。然后,在接收线圈中拾取的信号包含可以与身体中的不同位置相关联的不同频率的分量。经由接收线圈获得的信号数据对应于空间频率域并且被称为k空间数据。k空间数据通常包括利用不同相位编码采集的多条线。通过收集多个样本来对每条k空间线进行数字化。例如借助于傅里叶变换将一组k空间数据转换为MR图像。
在MR成像中,常常期望获得关于水和脂肪对总体信号的相对贡献的信息,以抑制它们中的一个的贡献或者分离地或共同地分析它们两者的贡献。如果组合来自在不同回波时间处采集的两个或更多个对应回波的信息,则可以计算这些贡献。这可以被认为是化学位移编码,其中,通过在略微不同的回波时间处采集两幅或更多幅MR图像来定义和编码额外的维度,化学位移维度。对于水-脂肪分离,这些类型的实验常常被称为Dixon型测量。借助于Dixon MR成像或Dixon水/脂肪MR成像,通过计算来自在不同回波时间处采集的两个或更多个对应的回波的水和脂肪的贡献来实现水-脂肪分离。通常这样的分离是可能的,因为在脂肪和水中存在氢的已知的进动频率差异。在其最简单的形式中,水和脂肪图像通过“同相”和“异相”数据集的相加或相减来生成。
近年来已经提出了若干Dixon型MR成像方法。除了用于水/脂肪分离的不同策略之外,已知技术的主要特征在于它们采集的回波(或“点”)的特定数量以及它们施加于使用的回波时间上的约束。常规的所谓两点和三点方法需要同相和反相回波时间,其中,水和脂肪信号分别在复平面中平行和反平行。三点方法已经逐渐被推广以允许灵活的回波时间。因此,它们不再将回波时间处的水和脂肪信号之间的角度或相位限制于某个值。以这种方式,它们在成像序列设计中提供了更大的自由度,并且具体地实现来自采集的信噪比(SNR)增益与分离中的SNR损失之间的权衡。另一方面,期望仅采样两个而不是三个回波以减少扫描时间。Eggers等人(Magn.Reson.Med.,65:96-107,2011)已经提出双回波灵活Dixon型MR成像方法。使用具有灵活回波时间的这样的Dixon型MR成像方法,同相和反相图像不再必要地采集,而是任选地由水和脂肪图像合成。
各种策略可以使用在Dixon成像中以用于采集两个不同回波时间处的回波信号,包括:(i)双通策略,其中,每个回波信号使用正幅度读出磁场梯度在激励RF脉冲之后分离地采集;(ii)飞回策略,其中,这两个回波信号使用正幅度读出磁场梯度与负幅度重绕机磁场梯度组合在相同激励RF脉冲之后采集;以及(iii)双极策略,其中,这两个回波信号在相同激励RF脉冲之后被采集,一个回波使用正幅度读出磁场梯度来采集,并且另一回波使用负幅度读出磁场梯度来采集。
双极策略提供许多有吸引力的优点,诸如较短的扫描时间、较高的SNR(信噪比)效率、更鲁棒的场图估计、降低的运动引起的伪影、以及对短T2*的较低的灵敏度。然而,交变读出磁场梯度引起若干问题,包括延迟效应和成像误配准,其限制现有的Dixon水/脂肪分离技术的直接应用。Lu等人(“Water-fat separation with bipolar multi-echosequences”,Magn.Reson.Med.2008,60,198-209)提出在水/脂肪分离之前应用对所采集的回波信号的后处理。该后处理包括移位k空间数据以校正k空间回波未对准并且基于低分辨率场图使图像翘曲以移除场不均匀性引起的误配准。实际的水/脂肪分离然后在k空间中执行,由此分离的水和脂肪图像中的化学位移引起的伪影被消除。
对应于选定的回波时间,适用于水/脂肪分离的Dixon技术的逆问题通常在k空间的中心部分中是良态的。然而,由于有限采样率,逆问题可能在k空间的更外围的部分是病态的。为此目的,Lu等人还提出通过对角加载的正则化以防止这些k空间区域中的过度噪声放大。然而,该方法引入水和脂肪信号贡献的估计中的偏置。
因此,本发明的目标是提供一种使用多回波采集技术(诸如双回波TSE或双回波GRASE采集或双极采集策略)使能高效和可靠的Dixon水/脂肪分离的方法。
发明内容
根据本发明,公开了一种对被放置于MR设备的检查体积中的对象(例如,患者的身体)进行MR成像的方法。所述方法包括以下步骤:
使所述对象经受成像序列,所述成像序列包括至少一个激励RF脉冲和切换的磁场梯度,其中,两个回波信号(第一回波信号和第二回波信号)在不同的回波时间处生成,
从所述对象采集所述回波信号,并且
根据所述回波信号重建水图像和/或脂肪图像,其中,来自水和脂肪的对所述回波信号的贡献使用k空间的第一区域中的两点Dixon技术和k空间的第二区域中的单点Dixon技术分离,其中,所述第一区域与所述第二区域不同。
本发明提出应用于第一和第二回波信号两者的水/分离的两点Dixon技术与应用于所述两个回波信号之一(即,第一回波信号数据或第二回波信号数据)的单点Dixon技术之间的自适应切换,这取决于k空间中的位置。备选地,单点Dixon技术可以个体地应用于这两个回波,并且结果可以然后组合以降低噪声。
优选地,两个回波信号根据本发明使用双极读出磁场梯度来采集。所述两个回波信号使用具有相反极性的一对时间相邻读出磁场梯度来采集。第一回波使用正幅度读出磁场梯度在第一回波时间处采集,并且第二回波使用负幅度读出磁场梯度在第二回波时间处采集(反之亦然)。
由于随着双极采集发生的延迟效应和图像误配准,所述两个回波信号优选地在重建水和脂肪图像之前针对k空间未对准来校正。所述两个回波信号还应当在重建所水和脂肪图像之前针对由于主磁场不均匀性的相位误差和图像空间未对准来校正。关于这些校正步骤的细节并且关于执行k空间参考中的水/脂肪分离由上文记载的Lu等人的文章实现。
如上所述,成像参数将通常被选择,使得应用于水/脂肪分离的两点Dixon技术的逆问题在k空间的中心部分中是良态的。然而,利用相同成像参数,两点Dixon分离的逆问题在k空间的更外围的部分中能够是病态的。因此,应用两点Dixon技术的第一k空间区域应当被确定以覆盖k空间的中心部分(其中,对应的逆问题是良态的)。本发明的洞察力在于,尽管两点Dixon技术变为病态,但是在该范例中朝向k空间的周围,应用于第一和/或第二回波信号数据的单点Dixon技术的逆问题变为(或者保持)足够良态。因此,应用单点Dixon技术的第二k空间区域优选地被确定为覆盖根据本发明的k空间的外围部分。
通常,本发明提出根据相应的k空间区域中的两点和单点Dixon技术的逆问题的预测调节确定第一和第二k空间区域。还可以考虑相应的Dixon分离技术的预测的水/脂肪偏置。最后,水图像和/或脂肪图像在水和脂肪贡献的分离之后使用k空间数据来重建。
根据本发明的优选的实施例,来自水和脂肪的贡献被假定为那些k空间区域中的k空间中的埃尔米特,其中,相应的两点或单点Dixon技术的逆问题是病态的。如之前所述,Dixon技术的逆问题的调节可以通过正则化改进。然而,该正则化引入水和脂肪信号贡献的估计的偏置。本发明的洞察力在于,逆问题的调节可以通过每当(无约束的)逆问题是病态的时将图像空间中的水和脂肪信号贡献的初始相位限于相等来改进(在不使由正则化造成的水和脂肪贡献的估计偏置的情况下)。水和脂肪信号贡献然后被假定为在图像空间中是实的,并且因此在k空间中是埃尔米特。逆问题的调节可以通过相应的水/脂肪分离算法的噪声放大来确定。逆问题可以被认为是病态的,例如,如果噪声放大的水平高于给定阈值。
Dixon算法假定来自水和脂肪的信号贡献是k空间中的埃尔米特的本发明的方法提供在两个回波信号上操作的k空间中的单点和两点水/脂肪分离。其可以局部用于选定的k空间区域(根据相应的逆问题的调节)或者甚至全局用于所有k空间区域。与公共复合水/脂肪分离不同,其常常在没有正则化的情况下保持良态。
因此到目前为止所描述的本发明的方法可以借助于MR设备来执行,所述MR设备包括:至少一个主磁体线圈,其用于在检查体积内生成实质上均匀的静磁场B0;多个梯度线圈,其用于在检查体积内的不同空间方向上生成切换的磁场梯度;至少一个身体RF线圈,其用于在检查体积内生成RF脉冲和/或用于接收来自被定位于检查体积中的患者的身体的MR信号,控制单元,其用于控制RF脉冲和切换的磁场梯度的时间演替;以及重建单元,其用于根据接收到的MR信号重建MR图像。可以通过MR设备的重建单元和/或控制单元的对应的编程来实现本发明的方法。
本发明的方法可以有利地在目前临床使用的大多数MR设备上执行。为此,仅需要利用计算机程序,通过所述计算机程序控制MR设备,使得其执行本发明的上述方法步骤。计算机程序可以存在于数据载体上或存在于数据网络中,从而下载以安装在MR设备的控制单元中。
附图说明
附图公开了本发明的优选实施例。然而,应该理解,附图仅出于说明的目的而被设计,而不是作为对本发明的限制的定义。在附图中:
图1示出了用于执行本发明的方法的MR设备;
图2示出了根据本发明的使用双极读出磁场梯度的Dixon成像序列的示意性(简化)脉冲序列图;并且
图3示出了将由k空间中的Dixon水/脂肪分离造成的噪声放大图示为k空间位置的函数的示图。
具体实施方式
参考图1,MR设备1被示出为框图。所述设备包括超导或常导主磁体线圈2,使得沿z轴通过检查体积创建基本上均匀的、时间上恒定的主磁场B0。设备还包括一组(1阶、2阶,以及-在适用的情况下-3阶)匀场线圈2’,其中,流过该组2’的个体匀场线圈的电流是可控的,以用于使检查体积内的B0偏差最小化的目的。
磁共振生成和操纵系统应用一系列RF脉冲和切换的磁场梯度以反转或激励核磁自旋,诱导磁共振,重聚焦磁共振,操纵磁共振,空间地和以其他方式编码磁共振,使自旋饱和等,以执行MR成像。
更具体地,梯度脉冲放大器3沿检查体积的x、y和z轴将电流脉冲施加到全身梯度线圈4、5和6中的选定梯度线圈。数字RF频率发射器7经由发送/接收开关8将RF脉冲或脉冲包发射到体RF线圈9,以将RF脉冲发射到检查体积中。典型的MR成像序列包括短持续时间的RF脉冲分段的包,其与任何施加的磁场梯度一起实现核磁共振的选定的操纵。RF脉冲用于使共振饱和,激励共振,反转磁化,重聚焦共振或操纵共振,并且选择被定位于检查体积中的身体10的部分。MR信号也通过身体RF线圈9来拾取。
为了生成身体10的有限区域的MR图像,一组本地阵列RF线圈11、12、13邻近针对成像所选择的区域被放置。阵列线圈11、12、13可以被用于接收通过经由身体RF线圈的RF发射引起的MR信号。
得到的MR信号由身体RF线圈9和/或阵列RF线圈11、12、13拾取,并由接收器14解调,接收器14优选地包括前置放大器(未示出)。接收器14经由发送/接收开关8连接到RF线圈9、11、12和13。
主计算机15控制匀场线圈2’以及梯度脉冲放大器3和发射器7,以生成本发明的成像序列。对于选定的序列,接收器14在每个RF激励脉冲之后快速连续地接收单条或多条MR数据线。数据采集系统16对接收到的信号执行模数转换,并将每条MR数据线转换成数字格式以用于进一步处理。在现代MR设备中,数据采集系统16是专用于采集原始图像数据的单独计算机。
最终,数字原始图像数据由重建处理器17重建成图像表示,重建处理器17应用傅立叶变换或其他适当的重建算法,诸如SENSE。MR图像可以表示通过患者的平面切片、平行平面切片的阵列、三维体积等。图像然后被存储在图像存储器中,在图像存储器中可以访问图像以将图像表示的切片、投影或其他部分转换成适当的格式以用于可视化,例如经由视频监测器18,视频监测器18提供得到的MR图像的人可读显示。
主计算机15和重建处理器17被编程为执行如上文并且在以下中所描述的本发明的方法。
在图2中,描绘了根据本发明的使用双极读出磁场梯度的Dixon成像序列的示意性脉冲序列图。示图示出了频率编码方向(M)、相位编码方向(P)和切片选择方向(S)上的切换的磁场梯度。此外,示图示出了RF激励脉冲以及由ACQ1和ACQ2指定的在其期间采集回波信号的时间间隔。示图覆盖一对回波信号的采集。多个这样的回波信号对通过使用不同相位编码(P)的描绘的序列的多次重复来采集以完全地覆盖k空间的所要求的区域。每对回波信号使用具有相反极性的对应的一对读出磁场梯度(M)采集。双极读出梯度的定时和幅度被选择以移位回波信号的采集窗口ACQ1、ACQ2,使得提供了不同回波时间TE1和TE2和相应地来自水质子和脂肪质子的信号贡献的不同相位偏移。这些信号贡献的Dixon型分离基于水和/或脂肪图像的重建的最后步的这些相位偏移。
在重建之前,两个回波信号针对k空间未对准和相位误差以及由于主磁场不均匀性的图像空间未对准被校正(参见例如Lu等人的“Water-fat separation with bipolarmulti-echo sequences”,Magn.Reson.Med.2008,60,198-209)。
以下重建涉及来自水和脂肪的对k空间中的回波信号的贡献的分离。为此目的,两点Dixon技术(参见例如Eggers等人“Dual-echo Dixon imaging with flexible choiceof echo times”,Magn.Reson.Med.2011,65,96-107)应用于k空间的第一区域中,并且单点Dixon技术(参见例如“A single-point Dixon technique for fat-suppressed fast 3Dgradient-echo imaging with a flexible echo time”,J.Magn.Reson.Imaging 2008,27,881-890)应用于k空间的第二区域中。单点Dixon技术可以应用于两个回波信号之一,即,第一回波信号数据或者第二回波信号数据。备选地,单点Dixon技术可以个体应用于这两个回波,并且结果可以然后组合以降低噪声。
与两点Dixon技术不同,单点Dixon技术将图像空间中的水和脂肪信号建模为实变量。在k空间中,这对应于共轭对称水和脂肪信号。为了利用该性质,水/脂肪分离联合在k空间中有利地被执行用于对应的负和正频率。备选地,单点Dixon技术要应用于的k空间中的回波信号可以被变换到水/脂肪分离的图像空间,因为脂肪位移在个体回波信号数据中是单极的并且可以简单地在水/脂肪分离之后校正。
本发明提出根据相应的k空间区域中的两点和单点Dixon技术的逆问题的预测调节确定第一和第二k空间区域。这被图示在图3中。示图示出了作为k空间位置(样本号)的函数的由k空间中的相应的水/脂肪分离算法造成的噪声放大,一次在没有(曲线31的虚线部分)和有(曲线31的实线部分)自适应正则化的情况下使用两点Dixon技术,如在Lu等人的上文记载的文章中,并且一次将单点Dixon技术应用到第一回波(实曲线32)、第二回波(虚曲线32),并且个体地应用到两者和结果的平均(曲线33)。在该范例中,假定具有1.5T处的-0.5ms和1.5ms的回波移位和0.75的部分回波因子的Dixon快速自旋回波(TSE)序列。接近k空间的中心(样本号256),采用的两点Dixon技术的噪声放大接近于0.5的最佳值,其对应于两倍平均。然而,朝向k空间的外围,噪声放大急剧增加(除非应用正则化)。如由实线31所指示的,在该范例中,噪声放大限于2.0的最大值。然而,该正则化引入水和脂肪信号贡献的估计的偏置。尽管噪声放大针对两点Dixon技术朝向k空间的外围变得不利(指示变为病态的逆问题),但是其对于应用于第一或第二回波的单点Dixon技术变为更有利。这从曲线32、33是明显的。不管哪里两个曲线32、33之一保持在曲线31下面,依赖于单点Dixon技术而不是两点Dixon技术是有利的。这表征本发明的意义内的第二k空间区域。噪声放大是沿着曲线31的最小值的中心k空间区域形成本发明的意义内的第一k空间区域。在某些k空间位置处(也在第二区域内),个体应用到第一和第二回波信号数据的单点Dixon技术的结果的平均提供噪声放大中的另外的降低。
最后,水图像和/或脂肪图像在水和脂肪贡献的前述分离之后使用k空间数据来重建。
两个复合回波信号S1和S2(在主磁场不均匀性引起的相位的估计和解调之后)可以通过以下来建模:
Figure BDA0002423417900000091
其中,
Figure BDA0002423417900000092
其中,W和F指代k空间中的水和脂肪信号,t指代采集时间,并且w和Δf指代脂肪的采用的谱模型的峰的相对幅度和频率。针对W和F的解由以下给出:
Figure BDA0002423417900000093
其中,
Figure BDA0002423417900000094
其中,λ指代正则化参数。
如上文参考图1-3所解释的,λ的适合的选择允许将噪声放大因子限于期望值。然而,矩阵CHC的对角加载将偏置引入水/脂肪分离中。为了避免这一点,本发明提出如果矩阵CHC的逆是病态的,则将图像空间中的水和脂肪信号的初始相位限于相等。水和脂肪信号贡献然后被假定为在图像空间中是实的,并且因此在k空间中是埃尔米特。仅考虑一个回波,S通过以下内容描述:
Figure BDA0002423417900000101
其中,下标RI指代实部和虚部,并且k2=-k1。包括第二回波也导致:
Figure BDA0002423417900000102
其构成八个方程的过定线性系统。在这两种情况下,W和F如上文导出,而无正则化。该受限的Dixon技术提供k空间中的一点和两点水/脂肪分离,其在k空间样本对上操作。其可以局部用于选定的k空间区域或者全局用于所有k空间区域。作为逆问题的调节的量度,图4的示图示出了作为k空间位置(样本号)的函数的由k空间中的相应的水/脂肪分离算法造成的噪声放大,一次在没有(曲线41的虚线部分)和有(曲线41的实线部分)正则化的情况下使用常规复合两点Dixon技术,并且在没有正则化的情况下使用提出的受约束的水/脂肪分离,仅考虑第一回波(曲线42)和两者回波(曲线43)。如从图4可以看到的,与常规的复合分离(曲线41)不同,受约束于实值水和脂肪信号贡献的分离(曲线42、43)在没有正则化的情况下保持良态。其噪声放大仅呈现跨k空间的微小变化并且保持在合理的界限内。
最后,参考图5示范性地说明了提出的受约束的Dixon水/脂肪分离。志愿者的骨盆利用3特斯拉MR设备上的3D T1加权的扰相双梯度回波序列成像。选择350(AP)x 350(RL)x200(FH)mm3的视场、0.7 x 0.7 x 2.5mm3的分辨率、2.0/3.7/5.5ms的TE1/TE2/TR和10°的翻转角。像素带宽总计为对应于单回波图像中的0.65像素的脂肪谱的主峰的信号的移位的650Hz。水和脂肪信号的主磁场不均匀性引起的相位和初始相位首先根据单回波图像确定并且消除,这忽略了任何误配准。在没有和有正则化的情况下根据常规Dixon算法的复合分离以及在没有正则化的情况下的提出的“混合复实”受约束分离然后应用于k空间中。图5示出了针对单个切片的结果。图5示出了在没有(左)和有(右)正则化的情况下使用常规的复合Dixon水/脂肪分离(顶行)和在没有正则化的情况下使用提出的混合复实分离(底部)重建的水图像。尽管正则化在常规复分离中是实质的以防止高频噪声击穿,但是其对于由本发明提出的方法是明显不需要的。常规观察为正则化的不利效应的图像锐度的损耗可以通过本发明的方法避免,而不遭受过度的噪声放大。

Claims (18)

1.一种对被放置在MR设备(1)的检查体积中的对象(10)进行MR成像的方法,所述方法包括以下步骤:
使所述对象(10)经受成像序列,所述成像序列包括至少一个激励RF脉冲和切换的磁场梯度,其中,在不同的回波时间(TE1、TE2)处生成两个回波信号,所述两个回波信号是第一回波信号和第二回波信号,
从所述对象(10)采集所述两个回波信号,
根据所述两个回波信号来重建水图像和/或脂肪图像,其中,使用k空间的第一区域中的两点Dixon技术和k空间的第二区域中的单点Dixon技术来分离来自水和脂肪的对所述回波信号的贡献,其中,所述第一区域与所述第二区域不同。
2.根据权利要求1所述的方法,其中,所述两个回波信号是使用双极读出磁场梯度采集的。
3.根据权利要求1或2所述的方法,其中,所述单点Dixon技术被应用于所述第一回波信号数据或者所述第二回波信号数据。
4.根据权利 要求1或2所述的方法,其中,所述单点Dixon技术在两个个体计算步骤中分别被应用于所述第一回波信号数据和所述第二回波信号数据,其中,来自水和脂肪的所述贡献是通过组合所述两个个体计算步骤的结果来计算的。
5.根据权利要求1-2中的任一项所述的方法,其中,在重建所述水图像和所述脂肪图像之前针对k空间未对准对所述两个回波信号进行校正。
6.根据权利要求1-2中的任一项所述的方法,其中,在重建所述水图像和所述脂肪图像之前针对由于主磁场不均匀性导致的相位误差和图像空间未对准对所述两个回波信号进行校正。
7.根据权利要求1-2中的任一项所述的方法,其中,所述k空间的第一区域覆盖k空间的中心部分,而所述k空间的第二区域覆盖k空间的外围部分。
8.根据权利要求1-2中的任一项所述的方法,其中,所述k空间的第一区域和所述k空间的第二区域是根据对相应的k空间区域中的所述两点Dixon技术和所述单点Dixon技术的逆问题的调节来确定的。
9.根据权利要求1-2中的任一项所述的方法,其中,来自水和脂肪的所述贡献被假定为其中相应的两点Dixon技术或单点Dixon技术的逆问题为病态的k空间区域中的k空间中的埃尔米特。
10.根据权利要求1-2中的任一项所述的方法,其中,所述两点Dixon技术和/或所述单点Dixon技术在没有正则化的情况下被应用。
11.一种对被放置在MR设备(1)的检查体积中的对象(10)进行MR成像的方法,所述方法包括以下步骤:
使所述对象(10)经受成像序列,所述成像序列包括至少一个激励RF脉冲和切换的磁场梯度,其中,在不同的回波时间(TE1、TE2)处生成两个回波信号,所述两个回波信号是第一回波信号和第二回波信号,
从所述对象(10)采集所述两个回波信号,
根据所述两个回波信号来重建水图像和/或脂肪图像,其中,使用两点Dixon技术和/或单点Dixon技术来分离来自水和脂肪的对所述回波信号的贡献,其中,来自水和脂肪的所述贡献被假定为其中相应的两点Dixon技术或单点Dixon技术的逆问题为病态的那些k空间区域中的k空间中的埃尔米特。
12.根据权利要求11所述的方法,其中,所述两个回波信号是使用双极读出磁场梯度采集的。
13.根据权利 要求11或12所述的方法,其中,所述单点Dixon技术在两个个体计算步骤中分别被应用于所述第一回波信号数据和所述第二回波信号数据,其中,来自水和脂肪的所述贡献是通过组合所述两个个体计算步骤的结果来计算的。
14.根据权利要求11-12中的任一项所述的方法,其中,在重建所述水图像和所述脂肪图像之前针对k空间未对准对所述两个回波信号进行校正。
15.根据权利要求11-12中的任一项所述的方法,其中,在重建所述水图像和所述脂肪图像之前针对由于主磁场不均匀性导致的相位误差和图像空间未对准对所述两个回波信号进行校正。
16.根据权利要求11-12中的任一项所述的方法,其中,所述两点Dixon技术和/或所述单点Dixon技术在没有正则化的情况下被应用。
17.一种MR设备,包括:至少一个主磁体线圈(2),其用于在检查体积内生成均匀的静磁场B0;多个梯度线圈(4、5、6),其用于在所述检查体积内的不同的空间方向上生成切换的磁场梯度;至少一个RF线圈(9),其用于在所述检查体积内生成RF脉冲和/或用于接收来自被定位于所述检查体积中的对象(10)的MR信号;控制单元(15),其用于控制RF脉冲和切换的磁场梯度的时间演替;以及重建单元(17),其用于根据接收到的MR信号来重建MR图像,其中,所述MR设备(1)被配置为执行根据权利要求1-16中的任一项所述的方法的步骤。
18.一种计算机可读介质,其上存储有要在MR设备上运行的计算机程序,所述计算机程序包括用于执行根据权利要求1-16中的任一项所述的方法的指令。
CN201880061999.7A 2017-08-24 2018-08-17 Dixon型水/脂肪分离MR成像 Active CN111133327B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
EP17187711.1 2017-08-24
EP17187711.1A EP3447517A1 (en) 2017-08-24 2017-08-24 Dixon-type water/fat separation mr imaging
US201862644714P 2018-03-19 2018-03-19
US62/644,714 2018-03-19
PCT/EP2018/072295 WO2019038192A1 (en) 2017-08-24 2018-08-17 MAGNETIC RESONANCE IMAGING (MRI) WITH WATER SEPARATION / DIXON TYPE GREASE

Publications (2)

Publication Number Publication Date
CN111133327A CN111133327A (zh) 2020-05-08
CN111133327B true CN111133327B (zh) 2022-09-02

Family

ID=59699591

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201880061999.7A Active CN111133327B (zh) 2017-08-24 2018-08-17 Dixon型水/脂肪分离MR成像

Country Status (5)

Country Link
US (1) US11327133B2 (zh)
EP (2) EP3447517A1 (zh)
JP (1) JP6814325B2 (zh)
CN (1) CN111133327B (zh)
WO (1) WO2019038192A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10823806B2 (en) * 2017-11-22 2020-11-03 Siemens Healthcare Gmbh Magnetic resonance imaging using dataset undersampling
CN112763955B (zh) * 2020-12-31 2022-12-20 苏州朗润医疗系统有限公司 一种基于图割算法的磁共振图像水脂分离方法
CN114137462B (zh) * 2021-11-19 2023-10-20 中国科学院深圳先进技术研究院 一种用于低场磁共振的多对比度成像方法和设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102692612A (zh) * 2011-03-22 2012-09-26 株式会社东芝 磁共振成像系统以及方法
CN102779327A (zh) * 2011-05-12 2012-11-14 上海联影医疗科技有限公司 一种基于欠采样k-空间数据的水脂分离方法
CN103257333A (zh) * 2012-02-17 2013-08-21 西门子(深圳)磁共振有限公司 一种磁共振成像中的水脂分离成像方法及装置
CN103513202A (zh) * 2012-06-16 2014-01-15 上海联影医疗科技有限公司 一种磁共振成像中的dixon水脂分离方法
JP2014511750A (ja) * 2011-04-21 2014-05-19 コーニンクレッカ フィリップス エヌ ヴェ 脂肪抑制のために化学シフト・エンコードを用いるコントラスト増強磁気共鳴血管造影法
WO2015091956A1 (en) * 2013-12-19 2015-06-25 Koninklijke Philips N.V. Mri with dixon-type water/fat separation with estimation of the main magnetic field variations
CN105556326A (zh) * 2013-09-16 2016-05-04 皇家飞利浦有限公司 具有dixon类型的水/脂肪分离的mr成像

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7042215B2 (en) * 2003-04-25 2006-05-09 Case Western Reserve University Three point dixon techniques in MRI spiral trajectories with off-resonance correction where each TE is a multiple of 2.2 milliseconds
WO2009012115A1 (en) * 2007-07-13 2009-01-22 Board Of Regents, The University Of Texas System Methods for efficient and improved phase-sensitive mri
US7777486B2 (en) * 2007-09-13 2010-08-17 The Board Of Trustees Of The Leland Stanford Junior University Magnetic resonance imaging with bipolar multi-echo sequences
JP2010019588A (ja) * 2008-07-08 2010-01-28 Nissan Motor Co Ltd 車両用ナビゲーションシステム及び車両用ナビゲーションシステムにおける位置情報の補正方法、並びに、情報配信サーバ及び車載ナビゲーション装置
JP5676559B2 (ja) * 2009-03-30 2015-02-25 コーニンクレッカ フィリップス エヌ ヴェ エコー時間の柔軟な選択を有する2ポイントDixon技術
RU2557334C2 (ru) * 2010-06-24 2015-07-20 Конинклейке Филипс Электроникс Н.В. Динамическая контрастная улучшенная мр визуализация с реконструкцией сжатого измерения
CN103620437B (zh) * 2011-06-27 2017-10-27 皇家飞利浦有限公司 使用具有fid和多梯度回波采集的超短回波时间脉冲序列以及水‑脂肪分离处理的骨mri
US8957681B2 (en) * 2012-02-20 2015-02-17 Wisconsin Alumni Research Foundation System and method for magnetic resonance imaging water-fat separation with full dynamic range using in-phase images
US9449386B2 (en) * 2012-02-28 2016-09-20 Board Of Regents, The University Of Texas System Method and apparatus for extended phase correction in phase sensitive magnetic resonance imaging
JP6467344B2 (ja) * 2012-09-04 2019-02-13 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. ディクソン水脂肪分離を用いるプロペラ
CN105103001B (zh) * 2013-04-03 2018-10-19 皇家飞利浦有限公司 使用高snr同相图像和较低snr至少部分地异相图像的dixon型水/脂肪分离mri
DE102013217651B4 (de) 2013-09-04 2016-03-24 Siemens Aktiengesellschaft Mehrpunkt Dixon-Technik
US9519040B2 (en) * 2014-04-04 2016-12-13 Board Of Regents, The University Of Texas System Methods and apparatus for phase correction in magnetic resonance imaging
WO2015190508A1 (ja) 2014-06-12 2015-12-17 株式会社 日立メディコ 磁気共鳴イメージング装置及び水脂肪分離画像作成方法
DE102015202062A1 (de) 2015-02-05 2016-08-11 Siemens Healthcare Gmbh Rekonstruktion von Magnetresonanzbilddaten für mehrere chemische Substanzarten bei Multi-Echo-Bildgebungsverfahren
EP3286574B1 (en) 2015-04-24 2022-11-09 Institut National de la Santé et de la Recherche Médicale (INSERM) Method for fat characterization using mri images acquired using a multiple gradient-echo sequence with bipolar gradients

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102692612A (zh) * 2011-03-22 2012-09-26 株式会社东芝 磁共振成像系统以及方法
JP2014511750A (ja) * 2011-04-21 2014-05-19 コーニンクレッカ フィリップス エヌ ヴェ 脂肪抑制のために化学シフト・エンコードを用いるコントラスト増強磁気共鳴血管造影法
CN102779327A (zh) * 2011-05-12 2012-11-14 上海联影医疗科技有限公司 一种基于欠采样k-空间数据的水脂分离方法
CN103257333A (zh) * 2012-02-17 2013-08-21 西门子(深圳)磁共振有限公司 一种磁共振成像中的水脂分离成像方法及装置
CN103513202A (zh) * 2012-06-16 2014-01-15 上海联影医疗科技有限公司 一种磁共振成像中的dixon水脂分离方法
CN105556326A (zh) * 2013-09-16 2016-05-04 皇家飞利浦有限公司 具有dixon类型的水/脂肪分离的mr成像
WO2015091956A1 (en) * 2013-12-19 2015-06-25 Koninklijke Philips N.V. Mri with dixon-type water/fat separation with estimation of the main magnetic field variations

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Dual-echo Dixon imaging with flexible choice of echo times;Holger Eggers 等;《Magnetic Resonance in Medicine 65》;20100921;第96-107页 *
Robust Spectral Denoising for Water-Fat Separation in Magnetic Resonance Imaging;Felix Lugauer 等;《Medical Image Computing and Computer-Assisted Intervention -- MICCAI 2015 》;20151231;第667-674页 *
water-fat separation with bipolar multiecho sequences;Wenmiao Lu 等;《Magnetic Resonance in Medicine》;20080625;第198-209页 *
基于多分辨率和区域增长的两点Dixon水脂分离技术;刘镖水 等;《南方医科大学学报》;20171231;第37卷(第2期);第245-250页 *

Also Published As

Publication number Publication date
RU2020111830A (ru) 2021-09-24
EP3447517A1 (en) 2019-02-27
CN111133327A (zh) 2020-05-08
EP3673281B1 (en) 2024-04-24
WO2019038192A1 (en) 2019-02-28
EP3673281A1 (en) 2020-07-01
JP6814325B2 (ja) 2021-01-13
US20200355772A1 (en) 2020-11-12
RU2020111830A3 (zh) 2021-10-25
US11327133B2 (en) 2022-05-10
JP2020527089A (ja) 2020-09-03

Similar Documents

Publication Publication Date Title
RU2592039C2 (ru) Формирование магнитно-резонансного изображения с использованием многоточечного способа диксона
US8866476B2 (en) MR imaging using parallel signal acquisition
US11041926B2 (en) Dixon-type water/fat separation MR imaging
CN111758041A (zh) 使用多梯度回波序列进行的Dixon MR成像
CN110720047B (zh) 双回波Dixon型水/脂肪分离MR成像
CN111133327B (zh) Dixon型水/脂肪分离MR成像
WO2011080693A1 (en) Water -fat separation in mri using partial echoes at different echo times
CN109716155B (zh) 具有迪克逊型水/脂肪分离的mr成像
EP2503348A1 (en) MRI using a multi-point Dixon technique and a low resolution calibration scan
WO2018114554A1 (en) Dixon-type water/fat separation mr imaging
CN109983358B (zh) Propeller mr成像
CN111164444B (zh) 具有经改进的脂肪位移校正的Dixon型水/脂肪分离MR成像
US20220057467A1 (en) Epi mr imaging with distortion correction
RU2773486C2 (ru) Магнитно-резонансная визуализация с разделением вода/жир методом диксона
US20230366962A1 (en) Dixon-type water/fat separation mr imaging
CN114402214A (zh) 狄克逊型水/脂肪分离mr成像

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