CN1846145A - MRI 中具有共享预脉冲的随机排序k-空间子集 - Google Patents

MRI 中具有共享预脉冲的随机排序k-空间子集 Download PDF

Info

Publication number
CN1846145A
CN1846145A CNA2004800256370A CN200480025637A CN1846145A CN 1846145 A CN1846145 A CN 1846145A CN A2004800256370 A CNA2004800256370 A CN A2004800256370A CN 200480025637 A CN200480025637 A CN 200480025637A CN 1846145 A CN1846145 A CN 1846145A
Authority
CN
China
Prior art keywords
profile
sub
profile sub
imaging
order
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.)
Pending
Application number
CNA2004800256370A
Other languages
English (en)
Inventor
G·M·贝克
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 CN1846145A publication Critical patent/CN1846145A/zh
Pending legal-status Critical Current

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/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
    • 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/563Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • 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
    • 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/56509Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • 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/56518Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field
    • 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/567Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
    • G01R33/5673Gating or triggering based on a physiological signal other than an MR signal, e.g. ECG gating or motion monitoring using optical systems for monitoring the motion of a fiducial marker

Landscapes

  • 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)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Vascular Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

一种磁共振扫描器(10)包括主磁体(12)、磁场梯度线圈(16)和一个或多个射频线圈(18)。序列发生器(60,62,64,66)通过将k-空间分成多个轮廓子集而产生抗伪影成像序列。该序列发生器(60,62,64,66)包括对每个轮廓子集内的轮廓进行排序以使每个子集内的轮廓顺序与其它轮廓子集顺序大体上不同的轮廓子集排序处理器(60,62,64),以及根据选定的长期轮廓排序对轮廓子集进行排序的轮廓子集采集排序处理器(60,66)。磁共振成像扫描器(10)执行抗伪影成像序列以生成成像数据。重建处理器(44)将生成的成像数据重建成基本上无伪影的重建图像。

Description

MRI中具有共享预脉冲的随机排序k-空间子集
以下涉及每个采集周期或激发(shot)获得多个轮廓(profile)的诊断成像技术。其尤其可应用于采用共享预脉冲的磁共振成像中,且将具体参照上述应用进行描述。然而,其还可应用于导航装置选通和跟踪序列或如心脏选通(cardiac gating)成像的其它同步成像序列中,以及可应用于有利地采用共享脉冲的其它磁共振应用中,如磁共振谱等。
在磁共振成像中,采集共享一个预脉冲的多个轮廓是减少比吸收率(SAR)和减少成像时间的已知技术。这些优点在基本上任何类型的磁共振成像中都是有价值的,且在诸如3特斯拉扫描器的高场强扫描器中具有特别的优点。该预脉冲可以是检测身体运动部分的当前位置的导航装置预脉冲或定时脉冲,如指示心动周期中出现选定心脏状态的心脏选通定时脉冲,而不是操纵磁共振的施加脉冲,如脂肪抑制预脉冲或磁化传递对比(MTC)预脉冲等。
通常用于这类应用的轮廓顺序是一种分段k-空间方法[例如,参见“Strategies to Improve Contrast in TurboFlash Imaging:Reordered Phase Encoding and K-Space Segmentation”,D.Chien,D.A.Atkinson和R.A.Edelman J.Magn.Res.Im., 1p.63(1991)]。通常用于该类成像(共享预脉冲)的分段采集方法产生良好的图像质量且不会导致较大的信号恶化或图像伪影。然而,采用分段采集方法限制了对采用为进行对比加权或运动抑制等而选择的轮廓排序的选择权。
例如,在对比增强MRA中,优选采用“中央”、“CENTRA”或其它长期轮廓排序在造影剂团通过感兴趣动脉时快速采集造影剂团的成像数据。中央、CENTRA或类似轮廓排序有利地首先采集k-空间中心周围的成像数据,接着采集远离k-空间中心的k-空间数据。类似地,选定的轮廓采集排序可在重建图像中提供有利的运动抑制或脂肪抑制。由于分段采集方法实行特殊的轮廓排序以减少伪影或其它假象,该排序不能同时可选成实现中央、CENTRA或其它有利的排序。
本发明提出一种克服上述及其它缺陷的改进的装置和方法。
根据一个方面,提供一种磁共振成像方法。将k-空间分成多个轮廓子集。为每个轮廓子集选择轮廓排序。采集多个轮廓子集的轮廓。按为每个轮廓子集选择的轮廓排序采集该子集的轮廓。按选定的轮廓子集采集顺序采集轮廓子集。将所采集的多个轮廓子集的轮廓重建成重建图像。
根据另一方面,公开一种磁共振成像装置。提供将k-空间分成多个轮廓子集的装置。提供为每个轮廓子集选择轮廓排序的装置。提供采集该多个轮廓子集的轮廓的成像数据采集装置,每个子集的轮廓按为该轮廓子集选择的轮廓排序采集,且按所选定的轮廓子集采集顺序采集轮廓子集。提供将所采集的多个轮廓子集的轮廓重建成重建图像的装置。
本发明的一个优点在于能够对轮廓子集进行长期排序选择,以实施对比加权、运动抑制或其它技术,同时通过在每个轮廓子集内选择短期排序大大减少伪影或其它图像假像。每个轮廓子集的大体不同的排序,如随机排序,减少了重建图像中的伪影,这些重建图像采用共享预脉冲和伪影减少与图像切片或三维图像体积的采集相结合。
另一个优点在于采用大体中央长期轮廓子集排序或另一种长期子集排序提供对比增强或呼吸保持MR成像,所述排序有利地在首先采集k-空间中心区域的同时,在每个轮廓子集内提供随机的或其它抑制伪影的排序。
在阅读下面优选实施例的详细描述后,本领域普通技术人员将更清楚本发明的许多其它优点和有益效果。
本发明可采用各种元件和元件布置,以及各种处理操作和处理操作安排的形式。附图仅用于举例说明优选实施例而不应解释为限制本发明。
图1示意性示出采用轮廓子集内具有随机或以其它方式大体不对称以及不相关轮廓排序的共享预脉冲进行磁共振成像的磁共振成像系统。
图2示意性示出采用三次激发的磁共振成像脉冲序列,每次激发包括共享预脉冲和每次激发中共享轮廓的随机排序,以及轮廓子集的交替低高长期采集排序。
图3示出由图2的脉冲序列进行的数据采集k-空间图,其包括预脉冲和每个轮廓共享顺序的指示。在右侧表示这些激发,示出了由这些激发采集的轮廓子集的交替低高长期采集排序。
图4示出采用多个轮廓子集的三维体积数据采集,其中每个轮廓子集的轮廓共享预脉冲,且按照大体中央采集顺序采集这些轮廓子集。
图5示出图4的轮廓子集的示例性随机共享排序。
参考图1,磁共振成像扫描器10包括主磁体线圈12,其优选为超导线圈,当然也可采用电阻主磁体线圈或永久性磁体。主磁体线圈12受激励,在检查区域14内产生基本上均匀的主磁场。磁场梯度线圈16在选定的空间方向上产生梯度以对由激励射频线圈18产生的磁共振进行空间编码。在图1中,示出了全身射频线圈18;然而,可采用局部线圈,如头部线圈、相控射频线圈阵列和SENSE线圈等代替或与全身射频线圈18协作以激发磁共振和/或检测磁共振回波。
磁共振序列控制器30调整和控制与全身射频线圈18或另一个射频线圈耦合以激发磁共振回波的射频发射器34,并且控制与梯度线圈16耦合以对被激发的磁共振回波进行空间编码的磁场梯度控制器32。与全身射频线圈18或另一个射频线圈耦合的一个或多个射频接收器36检测、解调磁共振回波且使其数字化,并将数字磁共振采样存储在k-空间存储器40内。重建处理器44进行基于傅里叶变换的重建或另一种类型的图像重建以由存储的k-空间磁共振采样生成一个或多个重建图像。
重建图像存储在图像存储器46内,由视频处理器50处理并显示在用户界面52上,通过局部计算机网络或因特网发送,或者进行其它处理。优选地,用户界面52包括显示器、打印机或其它使放射师或其它操作者能够观察、回顾或对重建图像进行其它操作的输出设备。此外,用户界面52优选使放射师或其它操作者能够与磁共振序列控制器30通信以建立磁共振成像序列、修改成像序列、执行成像序列或者控制磁共振成像扫描器10。
在采用共享预脉冲进行成像时,序列发生器60优选生成合适的成像序列,其中包括优选地将轮廓分组成轮廓子集。每个轮廓子集共享一个或多个预脉冲。预脉冲和下列轮廓子集在此也被称为激发。下面的方程(1)提供轮廓总数(N轮廓)、激发数(N激发)和共享因数(r)之间的关系:
Figure A20048002563700081
其中共享因数r指每个激发内的轮廓数。
为显著减少伪影,每个激发内的轮廓优选相对于其它激发的轮廓排序,使得该轮廓顺序与其它激发中的轮廓顺序大体不同且基本上不相关,且相对于其它激发中的轮廓顺序不对称。在一种方法中,序列发生器60访问随机数发生器62以在每个激发内产生随机轮廓顺序。应当理解,随机数发生器62可以是在计算机系统中通常采用的随机数发生器或伪随机数发生器。在另一种方法中,访问共享排序数据库64以检索每个轮廓子集的轮廓顺序,使得每个检索到的轮廓子集顺序与其它检索到的轮廓顺序大体不同或基本上不相关。在又一种方法中,序列发生器60从用户界面52接收合适的不同轮廓子集顺序的用户输入。
虽然每个轮廓子集或激发内的轮廓顺序优选是随机化的或以其它方式被选择成使伪影或其它伪像减少,优选从激发排序数据库66选择激发或轮廓子集的长期排序。激发或轮廓子集排序限定图像采集的长期轮廓排序,且能够被选作交替低高长期排序、中央排序或其它合适的排序。例如,在对比增强MRA中有利地采用中央、CENTRA或在采集远端k-空间数据之前首先优选采集大体位于中心的k-空间数据的其它排序。选自激发排序数据库66的中央、CENTRA或类似轮廓子集排序首先有利地采集大体位于k-空间中心周围的k-空间数据,而后采集离k-空间中心更远处的k-空间数据。类似地,选自激发排序数据库66的轮廓子集采集排序可在重建的图像中提供有利的运动抑制或脂肪抑制。
一旦由序列发生器60建立包括每个激发或子集内选定的短期轮廓顺序和选定的长期激发顺序在内的合适成像序列,磁共振序列控制器30就执行建成的序列以采集成像数据。成像数据由重建处理器44重建,而且得到的图像显示在用户界面52上或以其它方式使用。
继续参考图1并进而参考图2,其描述了由序列发生器60建立的示例性成像序列。该成像序列具有共享因数r=5,也就是说,每个预脉冲有5个共享轮廓。在示例性序列中,k-空间以单位步幅(unitarysteps)在ky=-7至ky=7(包含ky=-7和ky=7)的整数步中成像。轮廓数是15,且采用三个激发采集15个轮廓。这三个激发在图2中标记为“激发1”,“激发2”,“激发3”。每个激发包括预脉冲70和轮廓子集72。
预脉冲可以是基本上任何产生选定磁共振特性的磁共振激励或操作。例如,预脉冲70可以是脂肪饱和预脉冲、磁化传递对比(MTC)预脉冲等。子集72的每个轮廓具有其自己的激励脉冲,该激励脉冲不被轮廓子集的其它轮廓共享。或者,激励脉冲可由相应轮廓子集的少量轮廓共享。代替激励或操作磁共振或除此之外,预脉冲70可对应于能够检测身体运动部分当前位置的导航装置预脉冲,或者定时信号,例如使轮廓子集72的采集与心动周期同步的心脏选通信号。
预脉冲70引起或指示大体随时间变化的磁共振响应。在图2中,每个预脉冲70生成或指示随时间增加的与时间相关的响应。然而,应当理解,也可出现其它的磁共振响应对时间的其它依赖性,如随时间衰减的与时间有关的响应,这取决于选定的预脉冲。对于生理选通预脉冲,与时间相关的响应74是由心脏运动和血流变化等产生的空间变化磁共振。
与时间相关的响应74对于每个预脉冲70是基本上相同的。结果,每个轮廓子集74内的每个轮廓位置的与时间相关的响应74的值是基本上相同的。例如,在图2中,示例性增加的与时间相关的响应74对于每个轮廓子集74首先采集的轮廓具有最低响应值,且对于每个轮廓子集74的最后采集的轮廓具有最高响应值。优选地,预脉冲70和轮廓子集72在时间上相对设置以为轮廓子集74的中心轮廓提供最佳预脉冲性能(例如,在MTC预脉冲情形中的最佳加权,或者在心脏选通预脉冲情形中的最低程度运动的心肌)。对于包含5个轮廓的示例性轮廓子集74,优选对所采集的第三轮廓进行最佳化。
对于中心轮廓以外子集的轮廓,预脉冲为在中心轮廓之前或之后采集的轮廓提供次于最佳的性能。如果每个轮廓子集72内的轮廓排序基本对称或相关,则该非均匀响应会在重建图像中共同导致伪影或其它图像恶化。例如,如果每个轮廓子集内的轮廓从最低的ky值运行至最大的ky值,则具有对应于每个激发内预脉冲共享轮廓数r的重复周期的周期性通常被引入所采集的成像数据。当由重建处理器44处理时,k-空间内的这种周期性可导致图像伪影。
继续参考图1和图2且进而参考图3,为抑制这种图像恶化,每个轮廓子集内的轮廓具有大体不同的排序,该排序相对于其它轮廓子集的排序是不对称和不相关的。如图3中所最清楚看出的,示例性磁共振序列在相位编码(ky)和读出(kx)的方向为矩形平面切片采集k-空间数据,且每个轮廓对应于固定相位编码ky值处的读出kx值线。在该示范性实施例中,作为第1次激发“激发1”采集中心轮廓子集ky∈{-2,-1,0,1,2},其中括号{}表示无序集。“激发2”采集包含更高正性轮廓的轮廓子集ky∈{3,4,5,6,7}。“激发3”采集包含更高负性轮廓的轮廓子集ky∈{-7,-6,-5,-4,-3}。在图2和3的实施例中,轮廓子集的分布为交替型低高激发排序,从而首先采集k-空间中心部分,其在如造影剂辅助成像之类的时间敏感成像中是有利的,在造影剂辅助成像中在应用造影剂后在有限时间内可获得高对比成像数据。
在每次激发内,该轮廓子集的轮廓具有与其它轮廓子集的排序大体不同和不相关的排序。具体对于本实例,“激发1”具有轮廓排序ky=[0,2,-2,1,-1],其中括号[]表示有序集。“激发2”具有轮廓排序ky=[6,3,7,5,4],且“激发3”具有轮廓排序ky=[-5,-7,-4,-3,-6]。为使这些排序之间基本上不相关或不对称,对k-空间中的每个轮廓标示出共享顺序是有用的,如图3所示。在“激发1”中,ky=0被首先采集且其共享顺序=1,ky=2被第二个采集且其共享顺序=2,ky=-2被第三个采集且其共享顺序=3,ky=1被第四个采集且其共享顺序=4,ky=-1被最后采集且其共享顺序=5。按照ky值增加的顺序将共享排序列表,给出“激发1”的共享排序[3,5,1,4,2]。通过类似分析,“激发2”的共享排序为[2,5,4,1,3]且“激发3”的共享排序为[2,5,1,3,4]。这些共享排序表示在图3中“共享顺序”栏下。通常,每个轮廓子集的共享顺序值的范围为从1至r,其中r是共享因数且表示每个轮廓子集中的轮廓数。
在优选实施例中,序列发生器60分割k-空间,并且通过访问随机数据发生器62为每个轮廓子集计算合适的排序以计算应用于每个轮廓子集的大体不同的随机共享排序。通常,随机共享排序提供具有基本上相对于彼此不相关或不对称的排序的轮廓子集。当k-空间具有更多数据线时,这在大多数现代成像过程的普遍的,k-空间的每一段可填入在两次或更多次激发中采集的数据。
可选地,在为轮廓子集计算随机共享排序后,序列发生器进行每对轮廓子集的随机共享排序比较以检查不想要的相关性。例如,这种检查可指示“激发2”和“激发3”的共享排序之间的有疑问的相关性,因为这两个共享排序都以序列[2,5,...]开始。如果认为这种相关性是不想要的,该比较应会对这两个激发中的一个触发新的随机共享顺序计算。
代替当创建磁共振成像序列时使用序列发生器60计算随机排序,序列发生器60也可适当地访问共享排序数据库64以获得应用于每个轮廓子集的不同的存储共享排序。所存储的共享排序可采用随机数发生器62生成。在另一种方法中,手动选择基本上不相关且不对称的共享排序并将它们输入到共享排序数据库64内。在又一种方法中,对各种共享排序选择采用图像模拟,而且选择使模拟图像中的伪影极大减少的共享排序并将其存储到共享排序数据库64内。在一种经验性方法中,采用各种具有不同共享排序的轮廓子集的成像序列对基本上均匀的模型成像,而后选择产生基本无伪影的重建图像的序列的共享排序并将其存储到共享排序数据库64中。
每个轮廓子集的共享排序优选与其它子集的共享排序大体不同。优选地,各个轮廓子集大体不同的排序是唯一排序;也就是说,每个轮廓子集优选具有与该脉冲序列中全部其它轮廓子集不同的唯一共享顺序。
在某些情况下,大体不同的排序可包括一个或几个类似或同样的共享排序。考虑图2中的实例,每个轮廓子集72中的5个轮廓可以以5!=120种不同方式进行排序。因此,如果激发数超过120,某些排序将会重复。当在共享排序之间存在偶然随机重复性相似或存在共享排序的偶然重复时,仍可实现显著的伪影减少。因此,即使在轮廓子集之间存在偶然的共享顺序的相似或重复,仍可认为共享排序是大体上彼此不同的。
在优选实施例中,轮廓排序是与随机数发生器62协作随机生成的。然而,还可预见基于脉冲序列的基本原理分析以及图像重建过程计算轮廓子集的基本上非相关和不对称轮廓顺序。这种计算出的轮廓顺序不是随机的,但大体上彼此不同以避免显著的伪影。
磁共振序列可以基本上是任何类型的序列。大体平衡类型的序列获益于避免涡流伪影的子集划分。在一个优选实施例中,该序列是平衡场回波序列。
参考图4,其公开了对轮廓子集采用大体不同的轮廓排序以减少图像恶化,如伪影,该方法可用于非平面或非基于切片的三维成像。图4示出k-空间78的ky-kz维。在该示出的实例中,ky-kz维分量被分割成大体方形或矩形轮廓部分或子集80,每一个包含九个布置成在ky和kz方向上排列成3×3结构的轮廓。每个轮廓沿与ky和kz方向垂直的kx方向的读出方向采集,且每个轮廓子集80的九个轮廓共享同一个预脉冲(共享因数r=9)。优选沿轮廓子集80的大体中央排序82采集轮廓子集80,其中首先采集k-空间78的中心处的轮廓子集,而后轮廓子集采集顺序从k-空间78的中心向外盘旋。然而,也可采用其它的轮廓子集80的采集排序。
在优选实施例中,由y-和x-方向的视场比确定不对称因数。每个轮廓子集80内轮廓限定具有由不对称因数确定的y-z平面内各向异性的k-空间邻接段。
在一个可替换实施例中,共享因数可减少至r=3。三个预脉冲每个之后三个轮廓,在九个轮廓之间随机排序。在共享因数为r=5的另一个可替换实施例中,第二激发的第五个轮廓成为k-空间的下一个3×3部分的第一个轮廓。遗留的轮廓,如每个部分内的轮廓,可以是九个轮廓中随机选择的一个。
与参照图2和3所描述的平面切片采集相比,图4的中央采集由于在轮廓子集80的轮廓排序的对称性和相关性而可能遭受伪影或其它图像恶化的影响。通过采用序列发生器60为每个轮廓子集80产生不同的随机的或者以其它方式不相关和不对称的轮廓顺序而基本上抑制这种图像恶化。在这种情形下,每个共享排序由3×3矩阵限定,该矩阵的元素具有在1至r=9范围内的共享排序值。
参考图5,示出了图4的ky-kzk-空间78,其中为每个轮廓子集80显示随机计算的3×3共享顺序。每个轮廓子集80的共享顺序值为从1至r=9,且表示采集该轮廓子集80的激发内轮廓采集顺序。图5的共享顺序采用伪随机数据发生器计算,无需进行任何计算后比较来检测相似或相同的共享顺序。对于共享顺序r=9,应当理解,可能的共享排序数为9!=362,880。因此,两个轮廓子集80之间的共享排序恰好重复的可能性很低。
参照图4和5描述的排序采用每个3×3轮廓子集80内的随机排序,并且采用轮廓子集80的中央排序82以首先采集k-空间78的中心区域。例如当采集造影剂辅助成像数据时,该中央排序82是有利的,这是由于引入造影剂的对比度随时间衰减。然而应当理解,可以采用除示例性中央排序82以外其它的轮廓子集排序。对于基本上任何这类轮廓子集排序,每个轮廓子集80内的随机或以其它方式不对称和不相关排序有利地减少伪影或其它图像恶化。
已参照优选实施例描述了本发明。显然,在阅读和理解前述具体描述后,其它人员可进行修改和变更。本发明解释为包括全部这些修改和变更,只要它们落入附加的权利要求书及其等同描述的范围内。

Claims (18)

1、一种磁共振成像方法包括:
将k-空间分成多个轮廓子集;
为每个轮廓子集选择轮廓排序;
采集该多个轮廓子集的轮廓,每个子集的轮廓按为该轮廓子集选择的轮廓排序被采集,按选定的轮廓子集采集顺序采集这些轮廓子集;和
将所采集的多个轮廓子集的轮廓重建成重建图像。
2、根据权利要求1所述的磁共振成像方法,其中将k-空间分成多个轮廓子集满足:
其中N轮廓为被分割的k-空间的轮廓数,N激发为包括多个轮廓子集的轮廓子集总数,以及r是包含在每个轮廓子集内的轮廓总数。
3、根据权利要求1所述的磁共振成像方法,其中为每个轮廓子集选择轮廓排序包括:
选择每个轮廓子集内轮廓的顺序,其中每个轮廓子集内所选定的轮廓顺序大体与其它轮廓子集的所选定的顺序不同。
4、根据权利要求3所述的磁共振成像方法,其中轮廓采集包括为每个轮廓子集:
生成共享预脉冲;和
在该共享预脉冲的生成后采集该轮廓子集序列的轮廓。
5、根据权利要求4所述的磁共振成像方法,其中生成共享预脉冲包括以下之一:
生成操纵磁共振的预脉冲,
生成检测运动状态的导航装置预脉冲,和
生成指示心脏状态的心脏选通信号。
6、根据权利要求3所述的磁共振成像方法,其中:
每个轮廓对应于沿读出方向的一维线;和
每个轮廓子集内的轮廓限定了k-空间的邻接段。
7、根据权利要求6所述的磁共振成像方法,其中选定的轮廓子集采集顺序包括:
首先采集靠近k-空间中心的轮廓子集;和
在采集靠近k-空间中心的轮廓子集后,采集远离k-空间中心的轮廓子集。
8、根据权利要求1所述的磁共振成像方法,其中:
k-空间是三维k-空间且每个轮廓对应于三维k-空间的一维线;
通过y-和z-方向的视场比确定不对称因数;和
每个轮廓子集内的轮廓限定了具有由不对称因数确定的y-z平面内各向异性的k-空间邻接段。
9、根据权利要求1所述的磁共振成像方法,其中为每个轮廓子集选择轮廓排序包括:
选择每个子集内的轮廓顺序以使重建图像中的伪影最小化。
10、根据权利要求1所述磁共振成像方法,其中为每个轮廓子集选择轮廓排序包括以下之一:
随机排序每个轮廓子集内的轮廓,
基于为每个轮廓子集存储的共享顺序选择顺序,其中不同的存储共享排序用于选择每个轮廓子集内的轮廓的顺序,和
选择每个子集内的轮廓顺序,该顺序相对于其它轮廓子集的选定顺序大体不对称和不相关。
11、根据权利要求1所述的磁共振成像方法,其中采集轮廓包括:
对每个轮廓,激发共振并产生单个回波,在该回波期间采集该轮廓。
12、根据权利要求1所述磁共振成像方法,其中:
轮廓采集采用平衡场回波序列
13、一种磁共振成像装置包括:
用于将k-空间分成多个轮廓子集的装置(60);
用于为每个轮廓子集选择轮廓排序的装置(60,62,64);
用于采集该多个轮廓子集的轮廓的成像数据采集装置(10),其中每个子集的轮廓按为该轮廓子集选择的轮廓排序被采集,且按所选定的轮廓子集采集顺序采集轮廓子集;和
用于将所采集的多个轮廓子集的轮廓重建成重建图像的装置(44)。
14、根据权利要求13所述的磁共振成像装置,其中图像采集装置(10)生成单个共享预脉冲并以为相应轮廓子集选定的轮廓顺序采集该子集的多个轮廓。
15、根据权利要求13所述的的磁共振成像装置,其中用于为每个轮廓子集选择轮廓排序的装置(60,62,64)包括:
随机数发生器(62);和
用于访问随机数发生器(62)以为每个轮廓子集生成随机轮廓顺序的序列发生装置(60)。
16、一种磁共振成像装置包括:
磁共振成像扫描器(10),其包括:
主磁体(12),在成像区域(14)内生成主磁场,
磁场梯度线圈(16),横越成像区域(14)产生选定磁场梯度,和
一个或多个射频线圈(18),生成射频激励脉冲并检测磁共振信号;
通过将k-空间分成多个轮廓子集而生成抗伪影成像序列的序列发生器(60,62,64,66),该序列发生器(60,62,64,66)包括:
轮廓子集排序处理器(60,62,64),其对每个轮廓子集内的轮廓排序,使得每个子集内的轮廓顺序大体与其它轮廓子集的顺序不同,和
轮廓子集采集排序处理器(60,66),其根据选定的长期轮廓排序对轮廓子集进行排序;
扫描器控制器(30),其控制磁共振成像扫描器(10)执行抗伪影成像序列以生成成像数据,和
重建处理器(44),其将生成的成像数据重建成基本上无伪影的重建图像。
17、根据权利要求16所述的磁共振成像装置,其中主磁场为至少3特斯拉。
18、根据权利要求16所述的磁共振成像装置,其中扫描器控制器(30)通过使射频线圈(18)激发共振和控制射频线圈(18)和梯度线圈(16)中的至少一个以形成共振回波,且射频线圈(18)读出所形成的回波期间的相应轮廓,来为每个轮廓产生成像数据。
CNA2004800256370A 2003-09-08 2004-08-30 MRI 中具有共享预脉冲的随机排序k-空间子集 Pending CN1846145A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US50112203P 2003-09-08 2003-09-08
US60/501,122 2003-09-08

Publications (1)

Publication Number Publication Date
CN1846145A true CN1846145A (zh) 2006-10-11

Family

ID=34273019

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2004800256370A Pending CN1846145A (zh) 2003-09-08 2004-08-30 MRI 中具有共享预脉冲的随机排序k-空间子集

Country Status (5)

Country Link
US (1) US7362099B2 (zh)
EP (1) EP1664822A1 (zh)
JP (1) JP2007504877A (zh)
CN (1) CN1846145A (zh)
WO (1) WO2005024449A1 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806874A (zh) * 2009-02-12 2010-08-18 美国西门子医疗解决公司 用于排序代表mr图像数据的频域分量的系统
CN108369264A (zh) * 2015-12-16 2018-08-03 皇家飞利浦有限公司 用于对磁共振成像(mri)系统进行无线通信同步的系统和方法
CN108711448A (zh) * 2018-04-03 2018-10-26 复旦大学附属中山医院 一种扫描序列的处理方法、系统、存储设备及移动终端
CN113100741A (zh) * 2020-01-13 2021-07-13 上海联影医疗科技股份有限公司 一种磁共振扫描方法、设备及存储介质

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7265547B2 (en) * 2005-09-16 2007-09-04 General Electric Company Method and apparatus for acquiring MR data with a segmented multi-shot radial fan beam encoding order
WO2008041208A2 (en) * 2006-10-03 2008-04-10 Koninklijke Philips Electronics N.V. Segmented multi-shot mri involving magnetization preparation
WO2009047690A2 (en) * 2007-10-09 2009-04-16 Koninklijke Philips Electronics N.V. Segmented multi-shot mri involving magnetization preparation
EP3080633A1 (en) 2013-12-12 2016-10-19 Koninklijke Philips N.V. Mr imaging using multi-echo segmented k-space acquisition
DE102014202649B4 (de) * 2014-02-13 2015-12-10 Siemens Aktiengesellschaft Leise MR-Bildgebung durch eine variable Anzahl von Pulssequenzabschnitten zwischen zwei Vorpulsen
DE102014202606B4 (de) 2014-02-13 2015-10-01 Siemens Aktiengesellschaft Erstellung von Magnetresonanzaufnahmen unter Verwendung von Aufnahmeblöcken aus unterschiedlichen Schichten in unterschiedlichen Reihenfolgen
JP6596020B2 (ja) * 2014-04-22 2019-10-23 コーニンクレッカ フィリップス エヌ ヴェ 前パルス及びナビゲータを使用する運動補正を伴う磁気共鳴撮像
CN113099207B (zh) * 2021-03-31 2022-09-02 吉林工程技术师范学院 基于相位调制的微透镜阵列式深度学习三维鬼成像方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5122747A (en) 1990-11-21 1992-06-16 Mayo Foundation For Medical Education And Research Spiral three-dimensional fourier transform NMR scan
JPH05176912A (ja) * 1991-05-08 1993-07-20 Toshiba Corp 磁気共鳴イメージング装置および方法
US5499629A (en) 1994-02-01 1996-03-19 Board Of Trustees Of The Leland Stanford Junior University Slice profile stabilization for segmented k-space magnetic resonance imaging
JP3559597B2 (ja) * 1994-12-21 2004-09-02 株式会社東芝 Mri装置
JPH08196522A (ja) * 1995-01-25 1996-08-06 Toshiba Corp 磁気共鳴イメージング装置
US5912557A (en) 1997-02-19 1999-06-15 Mayo Foundation For Medical Education And Research Centric phase encoding order for 3D NMR data acquisition
US6025714A (en) 1997-04-25 2000-02-15 Toshiba America Mri, Inc. Magnetic resonance imaging (MRI) using fast spin echo (FSE) imaging process
US6144201A (en) 1997-12-26 2000-11-07 Kabushiki Kaisha Toshiba MR imaging utilizing ECG gating technique
EP1071971A4 (en) 1998-04-13 2004-04-21 Univ Pennsylvania SINGLE-PULSE RARE SEQUENCE PHASE NON-SENSITIVE PREPARATION FOR BROADCAST IMAGE FORMATION
JP3699304B2 (ja) * 1999-08-13 2005-09-28 ジーイー横河メディカルシステム株式会社 磁気共鳴撮像装置
DE19962476B4 (de) 1999-12-24 2004-04-08 Forschungszentrum Jülich GmbH Verfahren zur bildgebenden Untersuchung einer Probe mittels einer Aufnahmesequenz und Umordnung von Echosignalen
US6486668B1 (en) 2000-11-20 2002-11-26 Ge Medical Systems Global Technology Company, Llc Preparatory pulse sequence for suppression of artifacts in MR images
US6639211B1 (en) * 2000-11-22 2003-10-28 Koninklijke Philips Electronics, N.V. Contrast-enhanced MRA including an effective zero-latency method of bolus detection
WO2004079386A1 (en) 2003-03-03 2004-09-16 Koninklijke Philips Electronics N.V. Magnetic resonance imaging method

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806874A (zh) * 2009-02-12 2010-08-18 美国西门子医疗解决公司 用于排序代表mr图像数据的频域分量的系统
CN101806874B (zh) * 2009-02-12 2014-07-30 美国西门子医疗解决公司 用于排序代表mr图像数据的频域分量的系统
US9442961B2 (en) 2009-02-12 2016-09-13 Siemens Medical Solutions Usa, Inc. System for ordering frequency domain components representing MR image data
CN108369264A (zh) * 2015-12-16 2018-08-03 皇家飞利浦有限公司 用于对磁共振成像(mri)系统进行无线通信同步的系统和方法
CN108369264B (zh) * 2015-12-16 2021-08-10 皇家飞利浦有限公司 用于对磁共振成像(mri)系统进行无线通信同步的系统和方法
CN113419201A (zh) * 2015-12-16 2021-09-21 皇家飞利浦有限公司 用于对磁共振成像(mri)系统进行无线通信同步的系统和方法
CN108711448A (zh) * 2018-04-03 2018-10-26 复旦大学附属中山医院 一种扫描序列的处理方法、系统、存储设备及移动终端
CN113100741A (zh) * 2020-01-13 2021-07-13 上海联影医疗科技股份有限公司 一种磁共振扫描方法、设备及存储介质
CN113100741B (zh) * 2020-01-13 2024-02-27 上海联影医疗科技股份有限公司 一种磁共振扫描方法、设备及存储介质

Also Published As

Publication number Publication date
US7362099B2 (en) 2008-04-22
JP2007504877A (ja) 2007-03-08
EP1664822A1 (en) 2006-06-07
WO2005024449A1 (en) 2005-03-17
US20060273791A1 (en) 2006-12-07

Similar Documents

Publication Publication Date Title
Zhang et al. Accelerating parameter mapping with a locally low rank constraint
US10996306B2 (en) MRI system and method using neural network for detection of patient motion
CN103124517B (zh) Mri图像的多维迭代相位循环重构
US9897675B2 (en) Magnetic resonance fingerprinting (MRF) with simultaneous multivolume acquisition
CN103969611B (zh) 用于运行磁共振系统的方法和控制装置
US6794869B2 (en) Moving table MRI with frequency-encoding in the z-direction
CN105005011B (zh) 在三维快速自旋回波中的混合采集的磁共振成像方法
US10162037B2 (en) Navigator-based data correction for simultaneous multislice MR imaging
US20030011369A1 (en) Moving table MRI with frequency-encoding in the z-direction
US8155419B2 (en) MRI acquisition using sense and highly undersampled fourier space sampling
CN106981051B (zh) 用于运动解析的mri的系统和方法
EP2656093A1 (en) Parallel mri method using calibration scan, coil sensitivity maps and navigators for rigid motion compensation
CN104412118B (zh) 改进的磁共振采集的方法和系统
JPH06217960A (ja) 映画磁気共鳴撮像方法及び装置
EP1344081A2 (en) Signal analysis for navigated magnetic resonance imaging
CN1846145A (zh) MRI 中具有共享预脉冲的随机排序k-空间子集
US11696700B2 (en) System and method for correcting for patient motion during MR scanning
CN1102320A (zh) 梯度和自旋回波磁共振成像中的重影抑制
CN109814058A (zh) Caipirinha波space
EP1293794A2 (en) MRI acquisition method using the SENSE technique
CN100582760C (zh) 磁共振成像设备和磁共振图像产生方法
US9329251B2 (en) System and method for magnetic resonance imaging using multiple spatial encoding magnetic fields
KR20160011012A (ko) 자기 공명 영상 장치 및 그 제어 방법
WO2002042787A2 (en) Multishot echo planar imaging fluoroscopic technique
US20180292494A1 (en) Multi-contrast images from a magnetic resonance imaging scan

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Open date: 20061011