CN107592912A - 利用rf线圈灵敏度映射的并行mr成像 - Google Patents

利用rf线圈灵敏度映射的并行mr成像 Download PDF

Info

Publication number
CN107592912A
CN107592912A CN201680027582.XA CN201680027582A CN107592912A CN 107592912 A CN107592912 A CN 107592912A CN 201680027582 A CN201680027582 A CN 201680027582A CN 107592912 A CN107592912 A CN 107592912A
Authority
CN
China
Prior art keywords
ste
stimulated echo
echo signal
pulses
signal
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
CN201680027582.XA
Other languages
English (en)
Other versions
CN107592912B (zh
Inventor
K·内尔克
P·博尔纳特
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 Electronics 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 Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN107592912A publication Critical patent/CN107592912A/zh
Application granted granted Critical
Publication of CN107592912B publication Critical patent/CN107592912B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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/5611—Parallel 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
    • 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
    • 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/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/246—Spatial mapping of the RF magnetic field B1
    • 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/5611—Parallel 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
    • G01R33/5612—Parallel RF transmission, i.e. RF pulse transmission using a plurality of independent transmission channels

Landscapes

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

Abstract

本发明涉及一种对对象(10)进行MR成像的方法。本发明的问题是提供一种改进的MR成像技术，其实现在并行成像中使用的RF接收天线(11、12、13)的空间灵敏度轮廓以及B1和/或B0映射的快速和鲁棒的确定。本发明的方法包括使对象(10)经受受激回波序列。采集两个或更多个受激回波信号(STE、STE*)，即直接受激回波信号(STE)和共轭受激回波信号(STE*)，其中，经由具有不同空间灵敏度轮廓的两个或更多个RF接收天线(11、12、13)的阵列并行地接收受激回波信号(STE、STE*)中的至少一个，并且其中，经由具有基本均匀的空间灵敏度轮廓的体RF线圈(9)接收受激回波信号(STE、STE*)中的至少另一个。通过将经由RF接收天线(11、12、13)的阵列接收的受激回波信号(STE、STE*)与经由体RF线圈(9)接收的受激回波信号(STE、STE*)进行比较，来导出指示阵列的个体RF接收天线(11、12、13)的空间灵敏度轮廓的灵敏度图。此外，本发明涉及MR设备(1)和用于MR设备(1)的计算机程序。

Description

利用RF线圈灵敏度映射的并行MR成像
技术领域
本发明涉及磁共振(MR)成像领域。其涉及对身体的至少部分进行MR成像方法。本发明还涉及MR设备和要在MR设备上运行的计算机程序。
背景技术
利用在磁场与核自旋之间的相互作用以便形成二维或三维图像的图像形成MR方法在当今得到广泛使用，尤其在医学诊断领域中，这是因为对于对软组织的成像，它们在许多方面优于其他成像方法，不需要电离辐射并且通常不是侵入性的。
根据通常的MR方法，要检查的患者的身体被布置在强的均匀磁场(B0场)中，所述强的均匀磁场的方向同时定义测量所基于的坐标系的轴(通常为z轴)。取决于磁场强度，磁场产生针对个体核自旋的不同的能级。能够通过应用定义的频率(所谓的Larmor频率或MR频率)的电磁交变场(RF场，也称为B1场)激励(自旋共振)这些能级。从宏观的角度来看，个体核自旋的分布产生总体磁化，所述总体磁化能够通过在磁场垂直于z轴延伸时应用适当频率的电磁脉冲(RF脉冲)而偏离平衡状态，使得磁化执行关于z轴的进动运动。进动运动描述了圆锥的表面，所述圆锥的孔径角被称为翻转角。翻转角的幅度取决于应用的电磁脉冲的强度和持续时间。在所谓的90°脉冲的情况下，自旋从z轴偏转到横向平面(翻转角90°)。
在RF脉冲终止之后，磁化弛豫回到原始的平衡状态，其中，z方向上的磁化以第一时间常数T1(自旋晶格或纵向弛豫时间)再次建立，并且在垂直于z方向的方向上的磁化以第二时间常数T2(自旋-自旋或横向弛豫时间)弛豫。可以借助于一个或多个接收RF线圈来检测磁化的变化，所述接收RF线圈以在垂直于z轴的方向上测量磁化的变化的方式被布置并取向在MR设备的检查体积内。在例如90°脉冲的应用之后，横向磁化的衰变伴随有从具有相同相位的有序状态到在其中所有相位角均匀地分布(失相)的状态的核自旋的过渡(由局部磁场不均匀性所诱导的)。失相可以借助于重聚焦脉冲(例如，180°脉冲)补偿。这在接收线圈中产生回波信号(自旋回波)。
为了实现身体中的空间分辨率，沿着三个主轴延伸的线性磁场梯度被叠加在均匀磁场上，从而引起自旋共振频率的线性空间依赖性。然后，在接收线圈中拾取的信号包含可以与身体中的不同位置相关联的不同频率的分量。经由RF线圈获得的MR信号数据对应于空间频域，并且被称为k空间数据。k空间数据通常包括利用不同的相位编码采集的多条线。每条线通过收集多个样本来数字化。k空间数据的集合借助于傅里叶变换被转换为MR图像。
已知的并行采集技术被广泛地用于MR成像中以加速MR信号采集。此类别中的一种方法是SENSE(灵敏度编码)。SENSE和其他并行采集技术使用从多个RF接收天线并行获得的欠采样k空间数据采集，其中，RF接收天线具有不同的空间灵敏度轮廓。在这些方法中，来自多个RF接收天线的(复)信号数据以这种方式与复加权组合：抑制最终重建的MR图像中的欠采样伪影(混叠)。该类型的复RF线圈阵列信号组合有时被称为空间滤波，并且包括在k空间域或图像域中(SENSE中)的组合以及作为混合的方法。
在SENSE成像中，通常根据由SENSE参考扫描获得的低分辨率参考数据估计线圈灵敏度轮廓。然后，在图像重建期间使用该线圈灵敏度信息，以使用直接反演算法在图像空间中“展开”混叠像素。
通常，用于给定诊断成像任务的MR设备根据选定的成像序列的类型和参数自动检测何时需要SENSE参考扫描。SENSE参考扫描被自动插入到要执行的序列的列表中，通常紧接在诊断成像序列之前。
SENSE参考扫描通常包括两个扫描，因为必须(i)经由要确定空间灵敏度轮廓的多个RF接收天线，并且(ii)经由具有基本均匀的空间灵敏度的体RF线圈来采集MR信号数据作为参考。为了将RF接收天线的阵列与体RF线圈解耦，必须分离地执行这两个扫描。
此外，对于许多MR成像应用(特别是在3特斯拉或更高的高主磁场强度下)，对发送的RF场的空间分布的准确测量通常是重要的，以支持适当的前瞻性(如果适用)和回顾性校正/补偿。这需要鲁棒而快速的B1映射技术。然而，大多数B1映射技术相对慢，使得整合到临床工作流程中变得困难。国际申请WO2013/05006提及了根据采集的FID和受激回波信号形成指示身体(10)的部分内的RF脉冲的RF场的空间分布的B1图。
有时，主磁场B0的分布也需要在实际诊断扫描之前的准备扫描中确定。这实现MR图像重建期间的B0匀场和/或B0不均匀性的补偿。
所有这些不同的准备扫描显著增加了总体扫描时间。
发明内容
从上述容易意识到，需要一种改进的MR成像技术，其使得能够快速和鲁棒地确定所使用的RF接收天线的灵敏度轮廓以及B1映射和/或B0映射。
根据本发明，公开了一种对被放置在MR设备的检查体积中的对象进行MR成像的方法。所述方法包括以下步骤：
-使所述对象经受RF脉冲和切变的磁场梯度的成像序列，所述成像序列是受激回波序列，所述受激回波序列包括：
i)在磁化准备时段期间朝向所述对象辐射的至少两个准备RF脉冲，以及
ii)在时间上在所述磁化准备时段之后的采集时段期间朝向所述对象辐射的一个或多个读取RF脉冲；
-在所述采集时段期间的每个读取RF脉冲之后采集两个或更多个受激回波信号，即直接受激回波信号和共轭受激回波信号，其中，经由具有不同空间灵敏度轮廓的两个或更多个RF接收天线的阵列并行地接收所述受激回波信号中的至少一个，并且其中，经由具有基本均匀的空间灵敏度轮廓的体RF线圈接收受激回波信号中的至少另一个；并且
-通过将经由RF接收天线的阵列接收的受激回波信号与经由体RF线圈接收的受激回波信号进行比较，来导出指示阵列的个体RF接收天线的空间灵敏度轮廓的灵敏度图。
通常，受激回波序列包括三个(例如60°或90°)RF脉冲，其中，前两个RF脉冲是磁化准备脉冲。第一准备RF脉冲激励磁共振并将纵向核磁化变换为横向核磁化。第二个准备RF脉冲沿着纵轴“存储”失相横向核磁化的分数。在90°RF脉冲的情况下，该分数几乎是失相横向磁化的一半。在采集时段期间应用第三RF脉冲，所述采集时段在时间上处于准备时段之后。第三RF脉冲(“读取RF脉冲”)又将存储的纵向核磁化变换为横向核磁化，从而生成所谓的受激回波。可以通过由一连串低翻转角读取RF脉冲代替第三RF脉冲，来加速基于受激回波的MR成像，其中，每个读取RF脉冲仅重聚焦在准备周期之后存储的纵向核磁化的小部分。
本发明的方法采用最近介绍的DREAM B1映射方法(Medical Resonance inMedicine，68,1517-1526，2012；Magnetic Resonance in Medicine，71,246-256，2014)，其允许在短时间间隔内的B1图的采集。DREAM方法包括受激回波型准备阶段，其中，准备阶段的RF脉冲的局部有效翻转角被编码为纵向磁化。在随后的采集阶段中获得该翻转角信息，其中，FID信号(自由感应衰变)和受激回波信号在单个读出梯度波瓣下准同时地被采集。然后根据在单个采集串中采集的FID和受激回波信号的幅度导出指示局部翻转角的空间分布的B1映射(即准备阶段的RF脉冲的局部幅度)。
根据本发明，在采集时段期间的每个读取RF脉冲之后采集两个受激回波信号，即直接受激回波信号和虚拟的“共轭”受激回波信号(其表示受激回波的复共轭对应物)。优选地，直接受激回波信号和共轭受激回波信号被采集为梯度重聚回波信号。这使得能够准同时采集根据其导出空间灵敏度图的受激回波信号。从而，减少不希望的相位包裹和T2效应。
根据本发明，如在已知的DREAM方法中，在受激回波序列的采集时段期间应用一个或多个读取RF脉冲，其中，在每个读取RF脉冲之后准同时地采集两个或更多个受激回波信号。DREAM序列的两个(直接和共轭)受激回波信号几乎相同，因为它们从相同的受激回波编码的纵向磁化中产生。化学位移效应可以通过序列的适当定时被中和。因此，可以通过经由RF接收天线的阵列和经由体RF线圈两者从受激回波信号采集MR信号数据，在单次扫描中获得SENSE参考数据。两个或更多个受激回波被隔开大约1-2ms的短时间间隔，其足够长以在体RF线圈和RF接收天线的阵列之间切换信号接收。
可以从经由该RF接收天线接收的受激回波信号针对给定RF接收天线重建一幅MR图像。可以从经由体RF线圈接收的受激回波信号重建另一MR图像。在MR图像重建之后，可以根据两幅MR图像的逐体素比率导出相应的RF接收天线的空间灵敏度图。这可以对阵列的每个个体RF接收天线执行。
因此，本发明的方法通过应用单个DREAM型准备序列实现用于并行成像的快速灵敏度映射。
在本发明的优选实施例中，在采集时段期间也采集一个或多个FID信号。如在常规DREAM方法中，然后可以根据FID信号和受激回波信号的逐体素强度比率导出B1图。FID信号和受激回波信号被准同时采集。可以根据FID信号重建MR图像，并且可以根据受激回波信号重建另外的MR图像。在MR图像重建之后，B1图可以根据分别从FID和受激回波信号重建的这两幅MR图像的逐体素强度比率导出。
因此，本发明的方法使得不仅能够采集RF接收天线的空间灵敏度图，而且能够在单次预扫描中采集B1图。
需要采集具有适当相位编码的多个FID信号和受激回波信号以生成完整的空间灵敏度图和/或完整的B1图。与本发明相结合，可有利地应用如EPI的高效采样方案用于此目的。
根据本发明的优选实施例，在每个读取RF脉冲之后采集FID和受激回波信号作为梯度重聚回波信号。如上所述，可以调节序列的定时(回波时间)，使得对于两个直接和共轭受激回波信号和/或FID而言磁化率和化学位移诱导效应基本相等。换言之，这意味着本发明的成像序列的参数(回波时间、梯度强度)被选择为使得来自水自旋和来自脂肪旋转的对根据其导出空间灵敏度轮廓的受激回波信号和/或对FID信号的贡献基本相同。例如，来自水自旋的信号贡献和来自脂肪自旋的信号贡献可以在所有受激回波信号中基本上同相。这使得该技术对化学位移效应鲁棒。
根据本发明的另一优选实施例，根据所采集的FID信号和受激回波信号导出指示检查体积内的主磁场的空间分布的B0图。可以通过使用成像序列的适当的参数从FID信号和受激回波信号的逐体素复合信号中不仅导出B1图，而且导出B0图。因此，可以在没有额外的测量步骤的单个预扫描中同时采集RF接收天线的空间灵敏度图、B1图和B0图。
本发明的方法是MR成像中的全局准备扫描的潜在候选者，因为其可用于获得用于并行成像的灵敏度信息和B1/B0图。本发明的基于受激回波的灵敏度映射技术允许在大约100ms的持续时间中以受激回波序列的单次拍摄采集灵敏度图，这使得该方法原则上是能够实时的。
迄今为止描述的本发明的方法可以借助于MR设备来执行，所述MR设备包括：至少一个主磁体线圈，其用于在检查体积内生成均匀、稳定的磁场；多个梯度线圈，其用于在检查体积内生成不同的空间方向上的切变的磁场梯度；至少一个RF线圈，其用于在检查体积内生成RF脉冲，用于接收来自被定位于检查体积中的患者的身体的MR信号；控制单元，其用于控制RF脉冲和切变的磁场梯度的时间演替；以及重建单元，其用于根据接收到的MR信号重建MR图像。本发明的方法优选地借助于MR设备的重建单元和/或控制单元的对应编程来实现。
本发明的方法可以有利地在目前临床使用的大多数MR设备中执行。为此，仅需要利用控制MR设备使其执行本发明的上述方法的步骤的计算机程序。计算机程序可以存在于数据载体上或存在于数据网络中，从而被下载以安装在MR设备的控制单元中。
附图说明
附图公开了本发明的优选实施例。然而，应当理解，附图仅出于说明的目的被设计，而不是作为本发明的限制的定义。在图中：
图1示意性地示出了用于执行本发明的方法的MR设备；
图2示出了图示根据本发明的成像序列的示意图。
具体实施方式
参考图1，示出了MR设备1。所述设备包括超导或电阻式主磁体线圈2，使得沿着z轴通过检查体积创建基本上均匀的、时间上恒定的主磁场B0。所述设备还包括(一阶、二阶以及-在适用的情况下-三阶)匀场线圈的集合2’，，其中，出于最小化检查体积内的B0偏差的目的，通过集合2’的个体匀场线圈的电流是可控制的。
磁共振生成和操纵系统应用一系列RF脉冲和切变的磁场梯度来反转或激励核磁自旋，诱导磁共振，重聚焦磁共振，操纵磁共振，空间地和以其他方式编码磁共振，使自旋饱和等，以执行MR成像。
更具体地，梯度脉冲放大器3沿着检查体积的x、y和z轴向全身梯度线圈4、5和6中的选定的线圈应用电流脉冲。数字RF频率发射器7经由发送/接收开关8将RF脉冲或脉冲包发送到体RF线圈9，以将RF脉冲发送到检查体积中。典型的MR成像序列包括彼此一起采取的短持续时间的RF脉冲分段的包，并且任何应用的磁场梯度实现核磁共振的选定的操纵。RF脉冲用于进行饱和，激励共振，反转磁化，重聚焦磁共振或操纵共振，并选择被定位在检查体积中的身体10的部分。MR信号还由体RF线圈9拾取，所述体RF线圈具有基本均匀的空间灵敏度。
为了借助于并行成像生成身体10的限制的区域的MR图像，具有不同空间灵敏度轮廓的局部阵列RF线圈11、12、13的集合被放置成与选择用于成像的区域相邻。阵列线圈11、12、13可以被用于接收由身体线圈RF发送引起的MR信号。在并行发送应用中，阵列RF线圈11、12、13也可以用于RF发送，例如用于RF匀场的目的。
得到的MR信号由体RF线圈9和/或由阵列RF线圈11、12、13拾取并且由优选地包括前置放大器(未示出)的接收器14解调。接收器14经由发送/接收开关8连接到RF线圈9、11、12和13。
主计算机15控制通过匀场线圈2’以及梯度脉冲放大器3和发射器7的电流，以生成多个MR成像序列中的任何，诸如回波平面成像(EPI)、回波体积成像、梯度和自旋回波成像、快速自旋回波成像等。对于选定的序列，接收器14在每个RF激励脉冲之后快速连续地接收单个或多个MR数据线。数据采集系统16执行对接收到的信号的模数转换，并将每条MR数据线转换成适于进一步处理的数字格式。在现代MR设备中，数据采集系统16是专用于采集原始图像数据的单独的计算机。
最终，数字原始图像数据由重建处理器17重建为图像表示，所述重建处理器应用傅里叶变换或其他适当的重建算法，诸如SENSE或SMASH。MR图像可以表示通过患者的平面切片、平行平面切片的阵列、三维体积等。然后将图像存储在图像存储器中，其中，其可以被访问，以将图像表示的切片、投影或其他部分转换为适当的格式以用于可视化，例如经由视频监测器18，视频监测器18提供所得到的MR图像的人可读显示。
图2示出了图示根据本发明的成像序列的图。
所描绘的成像序列是受激回波序列，其被细分为磁化准备时段21和采集时段22。在磁化准备时段21期间应用具有α的翻转角的两个准备RF脉冲。两个准备RF脉冲被由时间间隔Ts分离。在两个准备RF脉冲之间应用失相磁场梯度Gm2。在时间上位于准备时段21之后的采集时段22期间生成具有翻转角β的读取RF脉冲的序列。每个读取RF脉冲β创建FID信号，受激回波STE和共轭受激回波STE*，它们在单个读出梯度波瓣下准同时被采集为单独的梯度重聚回波。采集顺序由失相梯度Gm2的极性确定。三个梯度回波的定时由切变的梯度Gm1、Gm(由下相位图ΦGr图示的)控制。
紧接在受激回波准备时段21之后，纵向核磁化由下式给出：
Mz,FID＝cos2(α)·M0
其中，Mz,FID指代未准备(即同相)纵向磁化，并且Mz,STE和Mz，STE*指代两个镜像受激回波准备(即，失相)纵向磁化贡献。横向磁化分量(即来自第一RF脉冲α的FID，来自第二RF脉冲α的FID和自旋回波)由强破碎梯度破坏，并且将不再被考虑。
根据本发明，经由具有不同空间灵敏度轮廓的RF线圈11、12、13的阵列并行地接收受激回波STE。经由具有基本均匀的空间灵敏度的体RF线圈9接收共轭受激回波信号STE*。两个受激回波STE、STE*由约1-2ms的时间间隔分离，其足够长以借助于发送/接收开关8来在RF线圈11、12、13阵列与体RF线圈9之间切换信号接收。FID信号可以例如还经由阵列RF线圈11、12、13并行接收。
因此，成像序列的读取RF脉冲β生成三个横向信号贡献：
IFID＝SA·C(t)sin(β)Mz,FID
ISTE＝SA·C(t-TS)sin(β)Mz,STE
ISTE*＝SBC·C(t+TS)sin(β)Mz,STE*
其中,SA和SBC中的每个表示复合系统常数，其包括针对给定体素的相应阵列RF线圈(SA)和体RF线圈(SBC)的接收线圈灵敏度，并且β是读取RF脉冲的标称翻转角。此外，如上所述，TS是在准备阶段中分离两个RF脉冲的时间间隔，并且C描述了由于磁化率和化学位移效应的针对给定体素的静态信号失相：
其中，ρ(r)和ω(r)指代质子密度和偏共振频率偏移，并且积分描述给定体素上的总和。当STE信号ISTE重聚焦为受激回波时，STE*信号ISTE*进一步失相，并且因此通常在常规的受激回波实验中被丢弃。然而，图2所示的成像序列采用定制的切变的磁场梯度以在有意选择的回波时间处采集所有三个信号贡献作为单独的回波梯度回波。针对测量梯度Gm、重定相梯度Gm1和受激回波失相梯度Gm2的梯度面积A的关系为：
A(Gm1)＝-1.5A(Gm)
第一个方程确保FID的梯度回波在第二读出梯度Gm的中心处重聚焦。第二个方程确保直接受激回波STE和共轭受激回波STE*分别在第一和第三读出梯度Gm的中心处重聚焦。采集顺序(STE-FID-STE*或STE*-FID-STE)由受激回波失相梯度Gm2相对于读出梯度Gm的极性来确定。
在由例如采集带宽或RF和梯度功率限制产生的一般序列定时约束内，可以独立地选择梯度回波时间(即梯度回波顶部的时间)以获得用于不同回波的期望的频谱编码。例如，可以通过将三个读出梯度Gm级联成单个恒定梯度波瓣来应用针对三个梯度回波的等距定时方案。另外，如果STE/STE*与FID信号之间的时间间隔ΔT被选择为等于TS，则STE和STE*信号具有相同的频谱编码时间，即TEFID，并且仅在T2演进时间上不同，所述T2演进时间对于受激回波STE为TEFID+2TS，对于共轭回波STE*为TEFID。因此，从STE信号和STE*信号重建的MR图像是相同的，并且差别仅在于T2弛豫和它们各自的系统常数SA和SBC。如上所述，可以选择定时，使得两个受激回波STE、STE*通过仅1-2ms的时间间隔分离，使得T2弛豫可以忽略。因此阵列RF线11、12、13的空间灵敏度图可以直接从对于阵列的每个单独的RF接收线圈11、12、13计算的受激回波STE和STE*的逐体素比导出(ISTE/ISTE*，见上文)，其可以针对阵列的每个个体RF接收线圈11、12、13被计算。
另外，如在已知的DREAM方法中(见上文)，受激回波准备RF脉冲的翻转角α(以及因此B1图)可以根据下式从采集的受激回波和FID信号的比率导出：
如果TS设置为TS＝ΔT+TEFID，则可以根据两个信号的相位导出B0相位图：
这适用于在脂肪/水同相回波时间(例如，在3特斯拉处的2.3ms)处获得FID的规定。在这种情况下，FID和受激回波是在不同的水/脂肪同相频谱编码时间(例如，在3特斯拉处，对于STE*为4.6ms，对于FID为2.3ms，对于STE为0ms)采集的。因此，STE信号和STE*信号也会在T2*弛豫上不同，但T2 *和T2效应至少部分抵消，因为STE受激回波具有比STE*受激回波更强的T2加权和更弱的T2 *加权。

Claims (9)

1.一种对被放置在MR设备(1)的检查体积中的对象(10)进行MR成像的方法，所述方法包括以下步骤：
-使所述对象(10)经受RF脉冲和切变的磁场梯度的成像序列，所述成像序列是受激回波序列，所述受激回波序列包括：
i)磁化准备时段(21)期间朝向所述对象(10)辐射的至少两个准备RF脉冲(α)，以及
ii)在时间上在所述磁化准备时段(21)之后的采集时段(22)期间朝向所述对象(10)辐射的一个或多个读取RF脉冲(β)；
-在所述采集时段(22)期间的每个读取RF脉冲(β)之后采集两个或更多个受激回波信号(STE、STE*)，即，直接受激回波信号(STE)和共轭受激回波信号(STE*)，其中，经由具有不同空间灵敏度轮廓的两个或更多个RF接收天线(11、12、13)的阵列并行地接收所述受激回波信号(STE、STE*)中的至少一个，并且其中，经由具有基本均匀的空间灵敏度轮廓的体RF线圈(9)接收所述受激回波信号(STE、STE*)中的至少另一个；并且
-通过将经由RF接收天线(11、12、13)的所述阵列接收的所述受激回波信号(STE、STE*)与经由所述体RF线圈(9)接收的所述受激回波信号(STE、STE*)进行比较，来导出指示所述阵列的个体RF接收天线(11、12、13)的所述空间灵敏度轮廓的灵敏度图。
2.根据权利要求1所述的方法，其中，在所述采集时段(22)期间采集一个或多个FID信号。
3.根据权利要求2所述的方法，其中，根据所述FID信号与所述受激回波信号的逐体素强度比率导出B1图。
4.根据权利要求1-3中的任一项所述的方法，其中，借助于多个连续的读取RF脉冲(β)来生成多个FID MR信号和受激回波(STE、STE*)MR信号。
5.根据权利要求1-4中的任一项所述的方法，其中，所述FID信号和/或所述两个或更多个受激回波信号(STE、STE*)被采集为梯度重聚回波信号。
6.根据权利要求1-5中的任一项所述的方法，其中，根据所采集的FID信号和受激回波信号(STE、STE*)导出指示所述检查体积内的主磁场的空间分布的B0图。
7.根据权利要求1-6中的任一项所述的方法，其中，所述成像序列的参数被选择为使得来自水自旋和来自脂肪自旋的对所述受激回波信号(STE、STE*)的贡献基本相同，所述空间灵敏度轮廓是根据所述受激回波信号导出的。
8.一种MR设备，包括：至少一个主磁体线圈(2)，其用于在检查体积内生成均匀、稳定的磁场；多个梯度线圈(4、5、6)，其用于在所述检查体积内生成不同空间方向上的切变的磁场梯度；至少一个RF线圈(9)，其用于在所述检查体积内生成RF脉冲和/或用于接收来自被定位于所述检查体积中的对象(10)的MR信号；控制单元(15)，其用于控制RF脉冲和切变的磁场梯度的时间演替；以及重建单元(17)，其用于根据接收到的MR信号重建MR图像；其中，所述MR设备(1)被布置为执行以下步骤：
-使所述对象(10)经受RF脉冲和切变的磁场梯度的成像序列，所述成像序列是受激回波序列，所述受激回波序列包括：
i)磁化准备时段(21)期间朝向所述对象(10)辐射的至少两个准备RF脉冲(α)，以及
ii)在时间上在所述磁化准备时段(21)之后的采集时段(22)期间朝向所述对象(10)辐射的一个或多个读取RF脉冲(β)；
-在所述采集时段(22)期间的每个读取RF脉冲(β)之后采集两个或更多个受激回波信号(STE、STE*)，即，直接受激回波信号(STE)和共轭受激回波信号(STE*)，其中，经由具有不同空间灵敏度轮廓的两个或更多个RF接收天线(11、12、13)的阵列并行地接收所述受激回波信号(STE、STE*)中的至少一个，并且其中，经由具有基本均匀的空间灵敏度轮廓的体RF线圈(9)接收所述受激回波信号(STE、STE*)中的至少另一个；并且
-通过将经由RF接收天线(11、12、13)的所述阵列接收的所述受激回波信号(STE、STE*)与经由所述体RF线圈(9)接收的所述受激回波信号(STE、STE*)进行比较，来导出指示所述阵列的个体RF接收天线(11、12、13)的所述空间灵敏度轮廓的灵敏度图。
9.一种要在MR设备(1)上运行的计算机程序，所述计算机程序包括用于以下操作的指令：
-生成RF脉冲和切变的磁场梯度的成像序列，所述成像序列是受激回波序列，所述受激回波序列包括：
i)磁化准备时段(21)期间辐射的至少两个准备RF脉冲(α)，以及
ii)在时间上在所述磁化准备时段(21)之后的采集时段(22)期间辐射的一个或多个读取RF脉冲(β)；
-在所述采集时段(22)期间的每个读取RF脉冲(β)之后采集两个或更多个受激回波信号(STE、STE*)，即，直接受激回波信号(STE)和共轭受激回波信号(STE*)，其中，经由具有不同空间灵敏度轮廓的两个或更多个RF接收天线(11、12、13)的阵列并行地接收所述受激回波信号(STE、STE*)中的至少一个，并且其中，经由具有基本均匀的空间灵敏度轮廓的体RF线圈(9)接收所述受激回波信号(STE、STE*)中的至少另一个；并且
-通过将经由RF接收天线(11、12、13)的所述阵列接收的所述受激回波信号(STE、STE*)与经由所述体RF线圈(9)接收的所述受激回波信号(STE、STE*)进行比较，来导出指示所述阵列的个体RF接收天线(11、12、13)的所述空间灵敏度轮廓的灵敏度图。
CN201680027582.XA 2015-05-13 2016-04-28 利用rf线圈灵敏度映射的并行mr成像 Active CN107592912B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP15167607.9 2015-05-13
EP15167607 2015-05-13
PCT/EP2016/059454 WO2016180640A1 (en) 2015-05-13 2016-04-28 Parallel mr imaging with rf coil sensitivity mapping

Publications (2)

Publication Number Publication Date
CN107592912A true CN107592912A (zh) 2018-01-16
CN107592912B CN107592912B (zh) 2020-12-01

Family

ID=53180566

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680027582.XA Active CN107592912B (zh) 2015-05-13 2016-04-28 利用rf线圈灵敏度映射的并行mr成像

Country Status (7)

Country Link
US (1) US10591565B2 (zh)
EP (1) EP3295194B1 (zh)
JP (1) JP6458170B2 (zh)
CN (1) CN107592912B (zh)
BR (1) BR112017024032A2 (zh)
RU (1) RU2683605C1 (zh)
WO (1) WO2016180640A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114325525A (zh) * 2020-09-29 2022-04-12 西门子医疗有限公司 省时的基于具有激励的回波的双回波序列的b0图的产生

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6445470B2 (ja) * 2013-03-15 2018-12-26 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． パラレルのマルチスライスｍｒ画像法
DE102017203082B4 (de) * 2017-02-24 2019-05-09 Siemens Healthcare Gmbh Parallele Bildgebung mit verbessertem Referenzdatensatz zur Bestimmung der Gewichtungsmatrix
US11354586B2 (en) 2019-02-15 2022-06-07 Q Bio, Inc. Model parameter determination using a predictive model
DE102020212272A1 (de) * 2020-09-29 2022-03-31 Siemens Healthcare Gmbh Echo-spezifische k-Raum-Abtastung bei Multi-Echo-Sequenzen
EP4040180A1 (en) * 2021-02-05 2022-08-10 Siemens Healthcare GmbH System and method for mri coil sensitivity estimation and reconstruction
WO2022271286A1 (en) * 2021-06-21 2022-12-29 Q Bio, Inc. Maxwell parallel imaging

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103649765A (zh) * 2011-04-11 2014-03-19 皇家飞利浦有限公司 具有b1绘制的mr成像
WO2014053927A1 (en) * 2012-10-02 2014-04-10 Koninklijke Philips N.V. Metal resistant mr imaging reference scan
CN104081219A (zh) * 2012-01-12 2014-10-01 皇家飞利浦有限公司 具有b1映射的mr成像
WO2015036340A1 (en) * 2013-09-10 2015-03-19 Koninklijke Philips N.V. Metal resistant mr imaging
WO2015059000A1 (en) * 2013-10-21 2015-04-30 Koninklijke Philips N.V. Mr imaging with temperature mapping

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4983921A (en) * 1989-08-18 1991-01-08 The Regents Of The University Of California Rapid calibration of nutation angles in MRI
US8334692B2 (en) * 2005-06-24 2012-12-18 Koninklijke Philips Electronics N.V. Simultaneous multinuclear magnetic resonance imaging
RU2451946C2 (ru) * 2006-12-22 2012-05-27 Конинклейке Филипс Электроникс Н.В. Рч катушка для использования в мр системе формирования изображения
US8198891B2 (en) * 2009-06-15 2012-06-12 General Electric Company System, method, and apparatus for magnetic resonance RF-field measurement
GB2492549B (en) 2011-07-04 2015-08-26 Pelsis Ltd Deterrent Device
DE102013201671B3 (de) * 2013-02-01 2014-07-31 Siemens Aktiengesellschaft Verfahren zum Erfassen von MR-Daten und zur Bestimmung eines B1-Magnetfelds sowie entsprechend ausgestaltete Magnetresonanzanlage
US10156625B2 (en) * 2013-08-12 2018-12-18 Koninklijke Philips N.V. MR imaging with B1 mapping

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103649765A (zh) * 2011-04-11 2014-03-19 皇家飞利浦有限公司 具有b1绘制的mr成像
CN104081219A (zh) * 2012-01-12 2014-10-01 皇家飞利浦有限公司 具有b1映射的mr成像
WO2014053927A1 (en) * 2012-10-02 2014-04-10 Koninklijke Philips N.V. Metal resistant mr imaging reference scan
WO2015036340A1 (en) * 2013-09-10 2015-03-19 Koninklijke Philips N.V. Metal resistant mr imaging
WO2015059000A1 (en) * 2013-10-21 2015-04-30 Koninklijke Philips N.V. Mr imaging with temperature mapping

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KAY NEHRKE等: "Volumetric B1+ Mapping of the Brain at 7T Using DREAM", 《MAGNETIC RESONANCE IN MEDICINE》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114325525A (zh) * 2020-09-29 2022-04-12 西门子医疗有限公司 省时的基于具有激励的回波的双回波序列的b0图的产生
CN114325525B (zh) * 2020-09-29 2024-04-19 西门子医疗有限公司 用于产生用于磁共振检查的b0图的方法和磁共振设备

Also Published As

Publication number Publication date
JP2018512977A (ja) 2018-05-24
BR112017024032A2 (pt) 2018-07-17
US20180120402A1 (en) 2018-05-03
CN107592912B (zh) 2020-12-01
US10591565B2 (en) 2020-03-17
EP3295194A1 (en) 2018-03-21
RU2683605C1 (ru) 2019-03-29
WO2016180640A1 (en) 2016-11-17
JP6458170B2 (ja) 2019-01-23
EP3295194B1 (en) 2018-09-19

Similar Documents

Publication Publication Date Title
CN104685368B (zh) 用于耐金属mr成像的方法和设备
US8405395B2 (en) Method for simultaneous multi-slice magnetic resonance imaging
US6603989B1 (en) T2 contrast in magnetic resonance imaging with gradient echoes
CN104204839B (zh) 使用apt对比增强和多回波时间采样的mr成像
CN107592912A (zh) 利用rf线圈灵敏度映射的并行mr成像
CN105143906B (zh) 金属抗mr成像
CN104204838B (zh) 具有Dixon类型水/脂肪分离和关于主磁场的不均匀性的先验知识的MRI
CN105531597B (zh) 耐金属的mr成像
US20140070805A1 (en) Mr imaging with b1 mapping
EP2902794A1 (en) MR imaging using multi-channel rf excitation
US8076939B2 (en) Method for fast magnetic resonance radiofrequency coil transmission profile mapping
CN106574954A (zh) 针对epi的具有奈奎斯特伪影校正的并行mr成像
CN102713658B (zh) 磁化率梯度绘图
CN102365559A (zh) 使用并行信号采集的mr成像
US20190041481A1 (en) Massively parallel magnetic resonance imaging wherein numerous off-surface coils are used to acquire partially under-sampled magnetic resonance signal data
CN105814449A (zh) 具有水/脂肪分离的零回波时间mr成像
CN103675735A (zh) 用于控制磁共振系统的方法和控制装置
CN102378910A (zh) 利用液体衰减反转恢复（flair）的双对比度mr成像
CN107548464A (zh) 自旋回波mr成像
CN107076818A (zh) 零回波时间mr成像
CN106133545A (zh) 利用对k空间中心的采样的零回波时间MR成像
CN107810425B (zh) 消除非t2加权信号贡献的t2加权mr成像
US6046588A (en) Magnetic resonance apparatus
CN109983358A (zh) Propeller mr成像
CN109716155A (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