CN103371820A - 借助磁共振技术采集呼吸的检查对象的测量数据组的方法 - Google Patents

借助磁共振技术采集呼吸的检查对象的测量数据组的方法 Download PDF

Info

Publication number
CN103371820A
CN103371820A CN2013101349354A CN201310134935A CN103371820A CN 103371820 A CN103371820 A CN 103371820A CN 2013101349354 A CN2013101349354 A CN 2013101349354A CN 201310134935 A CN201310134935 A CN 201310134935A CN 103371820 A CN103371820 A CN 103371820A
Authority
CN
China
Prior art keywords
excite
bunch
index
excites
space
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
CN2013101349354A
Other languages
English (en)
Other versions
CN103371820B (zh
Inventor
A.斯坦默
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.)
Siemens AG
Original Assignee
Siemens AG
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 Siemens AG filed Critical Siemens AG
Publication of CN103371820A publication Critical patent/CN103371820A/zh
Application granted granted Critical
Publication of CN103371820B publication Critical patent/CN103371820B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance 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/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/5676Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • G01R33/4826MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory in three dimensions

Landscapes

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

Abstract

不同的实施方式涉及用于借助三维磁共振技术确定呼吸的检查对象的磁共振图像的技术。检查对象的呼吸位置划分为多个预定义的簇。基于测量的检查对象的呼吸位置选择所述簇的特定的簇。待扫描的k空间区域划分为多个激发,所述激发分别借助特定的k空间轨迹采集多个k空间行。对于不同的簇一直采集不同的激发,直到至少两个紧挨的簇一起与测量数据组的所有激发对应。从与至少两个紧挨的簇对应的这些激发中确定MR图像。

Description

借助磁共振技术采集呼吸的检查对象的测量数据组的方法
技术领域
本发明涉及一种用于借助磁共振技术采集呼吸的检查对象的测量数据组的方法、一种磁共振设备、一种计算机程序以及一种电子可读数据载体。
背景技术
磁共振技术(以下对于磁共振简称MR)是能够用来产生检查对象的内部的图像的公知技术。简言之,为此将检查对象在磁共振设备中定位于具有0.2特斯拉至7特斯拉和更大的场强的相对强的静态均匀基本磁场,也称为B0场中,使得其核自旋沿着基本磁场取向。为了触发核自旋共振,将高频激励脉冲(RF脉冲)入射到检查对象中,测量所触发的核自旋共振作为所谓的k空间数据并且在此基础上重建MR图像或确定光谱数据。为了测量数据的位置编码,与基本磁场重叠快速接通的磁梯度场。所记录的测量数据被数字化并且作为复数值被赋予k空间矩阵中。从被赋值的k空间矩阵中例如借助多维傅里叶变换可以重建所属的MR图像。
要借助MR检查的患者的呼吸运动,可以在磁共振成像中(英语“magneticresonance imaging”-MRI)主要在检查胸部和腹部的器官时,也就是通过患者呼吸运动影响的检查区域,导致所谓的鬼影(英语“ghosting”)、模糊(英语“blurring”)、以及在所产生的图像之间的配准错误。这些伪影可以使得例如通过医生基于这些图像来诊断变得困难,并且可以导致,例如忽视损伤。现有技术中存在用来降低呼吸运动引起的伪影的大量技术。一种这样的技术是呼吸门控。呼吸门控是在MR测量期间采集患者的呼吸并且与采集的测量数据对应的技术。在呼吸门控中,仅当所采集的呼吸运动满足特定的可预先给出的标准时才将测量数据用于重建。
在此,可以利用外部传感器,例如气体垫、或者利用MR信号,所谓的导航器来检测患者的呼吸。导航器通常是短的序列,其采集例如隔膜的或检查对象中的另一个信号源的MR信号,所述信号源的运动与患者的呼吸相关。通过隔膜或另一个信号源的位置可以检测呼吸运动。
在利用导航器的呼吸门控中,将导航器序列例如与成像的序列交错并且将利用导航器测量的隔膜位置然后与直接之后(或之前)采集的成像的数据对应。
区分回顾性和前瞻性的呼吸门控。
在回顾性的呼吸门控中在MR测量期间采集并且记录、但是不分析呼吸运动,而是多次测量待采集的k空间。为了重建,仅考虑测量的数据的一部分,优选是其中呼吸信号在特定的窗中在杰出的呼吸位置附近的那个部分。如果特定的、对于图像重建来说需要的k空间数据点在杰出的窗内部被多次测量,则可以将所述数据进行平均。相反如果一个数据点总是在窗外部被测量,则可以将与杰出的位置的偏差最小的那个数据点用于重建。
在前瞻性的呼吸门控中在测量期间分析借助呼吸传感器测量的、生理的呼吸信号(例如利用导航器序列测量的隔膜位置)并且基于所采集的生理信号控制MR测量。在最简单的实施方式中,即所谓的接受/拒绝算法(Acceptance/Rejection-Algorithmus,ARA),将成像的数据分组的测量(和必要时对应的导航器序列)一直重复直到生理的信号落在事先定义的接受窗中。
这样的接受/拒绝算法的例子在Todd S.Sachs,Craig H.Meyer,Bob S.Hu,Jim Kohli,Dwight G.Nishimura和Albert Macovski的文章“Real-Time MotionDetection in Spiral MRI Using Navigators”,MRM32:第639-645页(1994)中被描述。作者在螺旋序列中对于每个激励采集一个或多个导航器。这些导航器在此在采集图像数据之后被采集。不同的导航器通过其空间取向相区别。从每个导航器中借助互相关来计算沿着导航器的方向相对于参考导航器的空间位移。在第一成像扫描之后被采集的导航器扫描分别用作为参考。特定的成像扫描被一直重复直到利用导航器确定的关于参考的空间位移小于由使用者预先给出的阈值。也就是涉及基于一个或多个空间位移的接受/拒绝算法。
接受/拒绝算法的另一个例子由Wang et al.在“Navigator-Echo-basedReal-Time Respiratory Gating and Triggering for Reduction of Respiratory Effects inThree-dimensional Coronary MR Angiography”.Radiology198;第55-60页(1996)描述。在此生理的信号是利用导航器确定的、隔膜位置相对于参考状态的位移。与Sachs et al.的工作的区别是,在成像的扫描之前或之后分别采集一个导航器,并且在于,仅当利用两个导航器确定的位移小于阈值时才接受成像的扫描。
为了确定接受窗,通常对于每个患者进行一个所谓的预扫描,在预扫描中例如利用导航器序列采集呼吸运动,但是还不采集成像的数据。
前瞻性的呼吸门控通常比回顾性的呼吸门控更有效。对于前瞻性的呼吸门控的前提条件是目前通常给定的、MR设备的控制软件的实时能力。在此,实时能力被理解为,可以在序列运行期间分析利用序列(序列在此包括成像的和导航器序列)测量的数据,并且通过该分析的结果可以影响其它序列流程,其中在采集数据和影响其它流程之间的时间段比呼吸运动的典型的时间常数短(在此人的呼吸循环特别地为3和10秒之间)。
接受/拒绝算法的主要问题是,患者的呼吸通常在检查期间变化。呼吸运动的变化在此可以是这样的,即,呼吸位置在一次确定的接受窗内部很少或不再被检测。这导致延长的采集时间并且甚至可以导致,测量根本不是常规地结束。
关注该问题的、以距离最为重要的算法是“Phase Ordering With AutomaticWindow Selection”(PAWS),如在P.Jhooti,P.D.Gatehouse,J.Keegan,N.H.Bunce,A.M.Taylor,和D.N.Firmin的文章“Phase Ordering With Automatic WindowSelection(PAWS):A Novel Motion-Resistant Technique for3D Coronary Imaging”,Magnetic Resonance in Medicine43,第470–480页(2000)中和在美国专利US7039451B1中描述的那样。PAWS找到用于运行时间的最终的接受窗并且由此可以灵活地响应于变化的呼吸。PAWS的另一个目的是,确保一定程度的“相位编码排序”(或简称“相位排序”)。这意味着,k空间中相邻的行在类似的呼吸状态中被采集。特别地,在运动特别敏感的k空间中心附近采集时应当避免呼吸状态的变化。PAWS是对于三维笛卡尔采集技术被开发的。在此采用的ky-kz排序方案在每个导航器之后采集三维k空间的一个完整的kx-kz平面。沿着kz方向的k空间信号的调制由于在沿着kx-kz方向通过导航器(以及可能接通的准备脉冲、或对另一个生理信号的等待,诸如EKG触发器)中断静态的均衡之后的超越状态而由此是平滑的。在ky方向上由于保留的运动而可以形成不连续,所述不连续可以作为沿着第一相位编码方向ky的伪影和模糊在图像中显示。这一点仅当在k空间中心附近形成片段边界时才不成立。由呼吸传感器不能采集的蠕动运动也可以导致图像中的伪影。
PAWS存在于不同的所谓“bin”变形中。在PAWS中确定最终的接受窗的宽度。包括了该接受窗的呼吸位置与接受/拒绝算法不同自动地在运行时间被找到。k空间填充按照簇进行。一个簇(在原始工作中使用“容器(bin)”的概念而不是簇)的特征通过呼吸位置区域、接受区域表示,并且包括在测量了在与簇对应的呼吸位置区域中的呼吸位置之后已经被测量的所有的k空间行。在PAWS的n-bin变形中由n个相继跟随的簇覆盖一个呼吸位置区域,其宽度等于接受窗。
此外,为每个簇对应一个k空间中的开始位置,其中不同的开始位置的数量是n。将具有相邻的呼吸位置的簇对于n>1个不同的开始位置进行对应。如果利用导航器测量与一个簇对应的呼吸位置,则开始在该簇内部还没有测量的k空间行的测量。对于在此选择哪些还待测量的k空间行的决定一般地也包括由相邻的簇已经采集的k空间行。例如这样选择尚缺少的k空间行,使得n个相邻的簇的任意组尽可能完整,其中n个相邻的簇的任意组包含与当前测量的呼吸位置对应的簇;即,n个相邻的簇的组包括尽可能多的不同k空间行。只要n个相邻的簇的任意组包括所有待测量的k空间行,则结束测量,因为呼吸位置的总变化在该测量数据中由此限制到接受窗。
n个不同的开始点和PAWS的n-bin变形的簇通常导致k空间中的n个片段。在此每个片段由相邻的k空间行组成。利用导航器测量的呼吸位置的变化在一个片段内部,相应于与一个簇对应(在原始工作中提到“bin尺寸”)的位置区域,并且由此相应于采集窗的n分之一。关于整个k空间,呼吸位置的变化更大并且向上受到确定的接受窗限制。属于相同的片段的行,在类似的呼吸状态期间被测量。在片段边界处由此信号的调制随着呼吸改变。也就是产生片段边界处的位置跳变。PAWS的不同的bin变形的目的是,将片段边界从运动敏感的k空间中心移开。另一个目的是实现高的效率。
在已经提到的Jhooti et al.的文章中以及例如也在P.Jhooti,P.Gatehouse,J.Keegan,A.Stemmer,D.Firmin的随后文章“Phase ordering with Automatic WindowSelection(PAWS)with Half Fourier for Increased Scan Efficiency and ImageQuality”;Proc.Intl.Soc.Mag.Reson.Med.11(2004);第2146页中将1-bin、2-bin、3-bin和4-bin变形互相比较。该比较的结果是,PAWS的1-bin和2-bin变形是最有效的,即,对于接受窗的给定宽度,测量最快结束。1-bin变形被放弃,因为其不允许“相位排序”,4-bin变形(和更高)由于效率低而被放弃。3-bin变形比2-bin变形效率低。对此的原因是,具有在左和右k空间的边沿处的开始位置的簇的单向增长方向。只要在这两个***簇和中央的簇(具有在k空间中心的开始位置和双向的增长方向)之间的间隙被封闭,则其一直增长直到在另一个***的簇和中央的簇之间的缝隙也被封闭,只要测量呼吸位置并且与第一***簇对应。这通常导致在簇边界(片段边界)处多次采集的k空间行。该问题在2-bin变形中不存在。在该2-bin变形中每个第二个簇单向地从左k空间边沿经过k空间中心在右k空间边沿的方向上并且其余的簇单向地从右k空间边沿经过k空间中心在左k空间边沿的方向上增长。只要(具有相对的增长方向的)两个相邻的簇“相遇”,则测量结束。然而,在如利用2-bin变形那样对称地扫描k空间的情况下簇边界通常位于运动特别敏感的k空间中心的附近,这可以导致强的图像伪影。对于靠近k空间中心的簇边界的可能性在使用部分傅里叶(也就是k空间的不对称扫描)的情况下明显更小。
实际上重要的由此是PAWS的所谓的2-bin和3-bin变形,其中在对称扫描的情况下3-bin变形、在不对称扫描的情况下2-bin变形是优选的。该分析基于2-bin变形,在所述2-bin变形中相邻的簇的开始点在左和右k空间边沿之间交替。相应地,簇分别从与其对应的开始点出发首先在k空间中心的方向上增长。
还要提到,在一些工作中仅将唯一一个呼吸位置与一个簇对应。最终的接受窗的宽度于是为呼吸信号的分辨率的n倍。在该替换的表达中通过如下实现了接受窗的灵活选择,即,将利用传感器测量的呼吸位置首先这样粗化,使得n个相邻的产生的呼吸位置覆盖一个与接受窗的宽度相应的呼吸区域。
从Anthony Nuval,Thanh D.Nguyen,Richard Watts和Yi Wang:“Animproved real-time navigator gating algorithm for reducing motion effects incoronary magnetic resonance angiography”;Journal of X-Ray Science andTechnology11(2003),第115–123页和A.Nuval,T.Nguyen,Y.Wang:“RefinedPAWS Algorithms for3D Coronary MRA”.Proc.Intl.Soc.Mag.Reson.Med.11(2003),第1625页中公知3-bin PAWS算法的三种修改:
a)在原始的3-bin PAWS变形中具有在左k空间边沿、k空间中心和在右k空间边沿处的开始位置的簇周期地交替。在该修改的版本中,开始位置在左k空间边沿、k空间中心中、右k空间边沿和再次k空间中心之间周期地交替。相应地,为每个第二簇分配一个在k空间中心中的开始位置。利用该修改避免了在簇边界处的位置跳变,其是为一个簇所分配的接受区域的两倍大小。然而该修改也可以减少簇组合的数量,在所述簇组合中可以补充完整k空间。也就是降低了效率。
b)这样收紧中断标准,使得中央的簇必须已经对称地围绕k空间中心采集k空间的至少30%。该修改的目的是,避免在k空间中心附近的簇边界。该修改一般地也延长测量时间,并且由此降低效率。
c)借助预扫描建立出现的呼吸位置的直方图。在预扫描期间最通常出现的呼吸位置被分配给中央的簇。该修改也减小了靠近k空间中心的片段边界的可能性。但是通过现在需要的预扫描进一步降低了效率。此外借助预扫描获得的信息仅在规则的呼吸的情况下可以传输到实际的扫描。由此预扫描信息的建立与实际的PAWS思路背道而驰,所述思路是,相对于变化的呼吸模式要是稳健的。
PAWS最初是对于ky-kz排序方案开发的,其中分别在采集呼吸信号之后,采集具有第二相位编码梯度(在kz方向上)的特定值的所有的k空间行。相应地,“相位排序”也限制到笛卡尔k空间方向,这可以导致在该方向上的加强的保留运动伪影。
在近期发表的文章中将PAWS与所谓的径向相位编码(RPE)方案进行组合(Christoph Kolbitsch,Claudia Prieto,Jouke Smink,and Tobias Schaeffter:“HighlyEfficient Whole-Heart Imaging Using Radial Phase Encoding-Phase Ordering WithAutomatic Window Selection”;Magnetic Resonance in Medicine66(2011);第1008–1018页)。在此分别按照导航器采集的数据具有相同的运动敏感性。执行特殊的2-bin方案。在一个bin组中按照顺时针方向采集在k空间中的径向轮辐,而在另一个bin组中逆时针地采集。该方案的目的是,可以多次在不同的呼吸阶段中重建感兴趣区域(ROI;英语:“region of interest”)。
发明内容
本发明要解决的技术问题是,提出一种方法、一种磁共振设备、一种计算机程序以及一种电子可读数据载体,用来降低迄今为止公知的PAWS方法的保留的运动敏感性。
按照一个方面,本发明涉及一种用于借助三维MR技术确定呼吸的检查对象的MR图像的方法。检查对象的呼吸位置划分为多个预先定义的簇,所述簇分别覆盖检查对象的呼吸位置的一个预先定义的呼吸区域。测量数据组通过多个激发采集,所述激发沿着预先给出的三维k空间轨迹分别采集预先给定数量的k空间行。不同的激发采集不同的k空间行,从而通过多个激发扫描预定义的k空间区域。为激发分别对应一个激发索引,其说明对于各自的激发相对于检查对象的运动的敏感性。激发索引与激发的对应是这样的,即,紧挨的激发索引与相对于运动具有递增改变的敏感性的这样两个激发对应。该方法包括A)借助导航器技术测量检查对象的呼吸位置和B)选择特定的簇,测量的呼吸位置位于其覆盖的呼吸区域。该方法还包括C)选择一个待采集的激发作为具有对应的激发索引的那个激发,该激发索引与最后对于特定的簇采集的激发的激发索引紧挨。该方法还包括D)采集待采集的激发并且将激发与特定的簇对应。该方法还包括一直执行步骤A)-D),直到覆盖了检查对象的相邻的呼吸区域的至少两个紧挨的簇一起与测试数据组的所有激发对应。该方法还包括从与两个紧挨的簇对应的激发中确定MR图像。
K空间行也可以被称为View。
通过对于每个簇这样选择待采集的激发,使得分别选择紧挨的激发索引用于采集,可以实现,***地和关于相对于运动的递增改变的敏感性排序地填充簇。例如,不同的簇的填充可以这样进行,即,首先采集具有相对于运动的大的敏感性(相对于运动的小的敏感性)的这些簇,然后***地和排序地采集具有相对于运动增加的小敏感性(相对于运动增加的大的敏感性)的这些激发。
通过激发索引与激发的对应,从而紧挨的激发索引表示相对于运动的递增地改变的敏感性,可以实现:将激发索引按照一种排序的方式与激发对应。例如,可以首先将激发于其相对于运动的敏感性排序并且然后将这些排序的激发,例如按照连续的或上升的或下降的方式与激发索引对应。
例如可以这样划分预定义的簇,使得检查对象的呼吸位置的整个呼吸区域被无缝隙地覆盖。这样总是可以,在上面讨论的步骤B)中选择特定的簇,测量的呼吸位置位于其覆盖的呼吸区域中。
例如采集与k空间中心具有相对小的距离的一个或多个k空间行的这些激发可以具有相对于运动的相对大的敏感性。相应地,不采集或仅采集k空间中心附近的小数量的k空间行的这些激发可以具有相对于检查对象的运动的相对小的敏感性。通过这样对应激发索引,使得其可以表示相对于运动的敏感性,可以实现,激发的选择和采集在考虑相对于运动的敏感性的条件下进行并且例如可以就此优化。
例如,关于相对于运动的敏感性对激发的采集的优化,可以通过在相邻的簇之间相关地选择和采集激发来实现。可能不必要的是,对于相邻的簇的一个被采集的激发对于另一个相邻的簇来采集。这一点可以是这样,因为MR图像可以由至少两个紧挨的簇中确定,在所述簇中相关地选择激发索引。
一个簇可以具有多于一个的邻居,例如两个邻居。以下考察例如可以对于任何邻居成立。
例如,对于覆盖了检查对象的相邻的呼吸区域的两个紧挨的簇的第一个首先选择并采集具有与相对于运动的大的敏感性相应的对应的激发索引的这些激发。对于两个紧挨的簇的第二个可以首先选择和采集具有与相对于运动的小的敏感性相应的对应的激发索引的这些激发。
换言之,对于两个紧挨的簇的第一个来说待采集的激发的选择包括:选择至此对于两个紧挨的簇的第一个还没有被测量的并且具有相对于该至此还没有被测量的激发的运动的最大的敏感性的那个激发。对于两个紧挨的簇的第二个来说待采集的激发的选择可以包括:选择至此对于两个紧挨的簇的第二个来说还没有被测量的并且具有相对于该至此还没有测量的激发的运动的最小的敏感性的那个激发。对于两个紧挨的簇的一个还没有被测量的激发,也可以是对于两个相邻的簇的另一个还没有被测量的。
对于不同的激发索引的具体选择可以使用不同的技术。例如,连续增加的激发索引可以表示各自的激发相对于检查对象的运动的连续上升或下降的敏感性。但是激发索引的另外的变形也是可以的,例如在附加考虑各自的k空间轨迹的k空间组件的条件下。
在按照本发明的用于借助磁共振技术采集呼吸的检查对象的测量数据组的方法中,其中测量数据组通过多个激发采集,所述激发分别包括数量Nv个k空间轨迹,即所谓的View,并且其中在激发的测量之前利用合适的传感器测量当前的呼吸信号,换言之在不同的实施方式中测量数据组的采集可以包括以下步骤:
-A)为每个激发分配一个激发索引所述激发索引描述了激发在k空间中的相邻关系,其中,这样排序激发索引,使得激发相对于检查对象的运动的敏感性从激发索引ns=0到激发索引ns=ns0增加并且从激发索引ns=ns0至激发索引ns=Ns-1又降低,
-B)将在呼吸信号的测量之后待测量的激发与一个簇相应于测量的呼吸信号进行对应,其中,按照唯一的方式为一个簇对应一个呼吸信号区域,并且一个簇包括在呼吸信号的测量之后在簇的呼吸信号区域中已经被采集的激发的所有激发索引,
-C)根据由前面选择的簇和其直接的相邻的簇已经采集的激发索引选择对于待测量的激发的激发索引,和在呼吸信号的测量之后采集属于该索引的激发,其中簇的相邻关系相应于其所分配的呼吸区域进行定义,
-D)一直采集导航器测量和激发,直到相邻的簇共同包括具有所有Ns个激发索引的激发,其中,
对于待测量的激发的激发索引的选择可以这样进行,使得一个簇分别双向地从已经由该簇采集的激发索引增长,其中,该增长在索引区域的边缘周期地继续并且为相邻的簇交替地分配在具有相对于运动的高的敏感性的索引区域中的中央的开始位置,和在索引区域的边缘处的***的开始位置,所述开始位置然后作为对于待测量的激发的激发索引来选择,只要在与待测量的激发对应的该簇中至此还没有采集激发。
各个中央的开始位置,即例如具有相对于检查对象的运动的大的敏感性的,和***的开始位置,即例如具有相对于检查对象的运动的小的敏感性的,与相邻的簇的分配,以及两个簇的双向的增长方向,主动地避免了,在运动特别敏感的k空间中心附近形成簇边界并且由此降低了对于运动伪影的易发性。
通过上面提到的双向选择激发索引可以实现,交替地选择k空间轨迹,其位于与定义了k空间行的相位编码梯度方向垂直取向的平面的不同半边。
换言之,待采集的激发的选择还可以这样进行,即,尽可能均匀地沿着不同的方向和/或k空间区域的不同侧扫描k空间区域。
前面主要描述了如下的技术,所述技术涉及基于对于激发的预定义的集合的激发索引选择和采集特定的激发。然而,按照不同的实施方式也可以进行激发本身的集合的确定。特别地,可以这样确定激发或激发索引,使得它们具有相对于运动的特征性的敏感性。这一点例如可以通过将k空间区域划分为扇形来进行,所述扇形具有与k空间中心的明确定义的距离。该距离特别地可以表示相对于特定的采集的k空间行的运动的敏感性。
在一种实施方式中,在k空间中按照笛卡尔的行采集待采集的测量数据组,其中包括以下步骤:
-对每个激发选择数量Nv个View,
-确定数量Ns个激发,以便完整地填充待扫描的k空间,
-将View这样与k空间中的Nv个扇形对应,使得在每个扇形中布置近似相同数量的View,并且一个扇形中的所有View分别具有与k空间中心相似的距离,
-将每个扇形的每个View与Ns个激发中的一个相应于其在kz-ky平面中的取向进行对应,
-这样扫描每个激发的View,使得与相同的扇形和不同的激发对应的View分别占据激发内的相同的(时间)位置。
k空间行或View可以沿着相位编码梯度场方向布置,其中,多个激发的k空间轨迹或View分别主要径向地在垂直于相位编码梯度场方向取向的平面中布置。
主要径向地可以意味着,不同激发的k空间轨迹或View也可以具有切线的分量,但是所述分量相对于其径向分量来说相对小。
可以将k空间区域划分为扇形,所述扇形关于其与k空间中心的距离定义,其中对于多个激发的每一个来说,分别一个k空间行位于每个扇形中。
例如,每个激发的k空间轨迹可以按照顺序扫描不同的k空间行,所述顺序相应于k空间区域中的、不同的k空间行所位于其中的扇形的相邻关系。
k空间的这样的扫描相对于检查对象内部的运动、例如包括蠕动性运动来说是稳健的,因为运动沿着两个笛卡尔方向被“分散”,其方法是,一个激发的View相应于其在kz-ky平面中的取向而不是如迄今为止的那样沿着ky行并且由此沿着仅一个相位编码方向被获取。由此,该方法不太容易发生例如在现有技术中采用的分割的情况下沿着相位编码方向会形成的鬼影伪影。
以简单的方式可以通过kz-ky平面中的极坐标系中的View的方位角来确定取向。由此一个扇形的View与Ns个激发中的一个的对应,可以相应于其在kz-ky平面中的极坐标系中的方位角来进行。
与一个激发对应的View对于每个激发来说按照相同的顺序被扫描。也就是说,与一个共同的扇形对应的View在其激发中在各自的激发开始之后的相同的时间被获取。由此得到沿着激发的径向方向和由此与方位角方向正交的k空间的平滑调制,由此如上所述避免了伪影鬼影。
在每个激发中的View的扫描顺序例如可以相应于激发的View所对应的扇形来选择。因为在一个扇形中的View全部具有与k空间中心的相似距离,所以例如可以将该距离作为对于顺序的排序标准来考虑,由此用来扫描每个激发的View的顺序相应于扇形的相邻关系。
本发明基于以下思路:
门控技术在笛卡尔三维梯度回波序列的方面是特别重要的。这些序列通常在每个激励脉冲之后采集一个特定的k空间行,所述k空间行通过在第一相位编码方向的ky值和在第二相位编码方向上的kz值规定。该相位编码行一次或多次在不同的回波时间(例如在使用迪克松技术时)被读取。这样的单激励(包括信号编码和信号采集)的持续时间为仅数毫秒。由此在用于测量生理的信号诸如呼吸运动的唯一一个导航器序列之后和/或在执行用于抑制不期望的信号份额例如用于抑制脂肪信号的预接通模块之后通常采集多个相位编码行。
在特定的导航器序列之后被采集的所有相位编码行的集合,以下称为“激发(Shot)”。因为在通过导航器序列(和可能的其它预接通模块,例如用于脂肪饱和的)中断静态均衡之后的磁化处于超越状态中,所以在一个激发内ky-kz行的时间顺序确定k空间的调制和由此确定图像质量。
此外公知,中央的k空间行是运动最敏感的,并且一个特定的k空间行的运动敏感性随着其与k空间中心的距离而下降。运动敏感性在此理解为通过在检查的检查对象中的运动引起的伪影形成的容易程度。
PAWS最初是对于ky-kz排序方案开发的,其中在一个激发内采集具有第一相位编码方向ky的相位编码梯度的特定值的所有k空间行。在该排序方案中对每个激发的激励的数量等于在第二相位编码方向上的相位编码步骤Nz的数量。最终为了图像重建而被接受的激发的数量,等于在第一相位编码方向上的相位编码步骤Ny的数量。相应地在此得到,k空间的调制由于在中断静态均衡之后的超越状态而仅沿着第二相位编码方向延伸。相反,保留的运动伪影沿着第一相位编码方向出现。此外单个激发的运动敏感性通过第一相位编码方向ky的值来描述特征。
为了减小运动敏感性和此外为了达到PAWS的一般描述,以下假定:
a)每个激发的k空间行、或更一般地、“View”的数量是恒定的。View的概念也包括非笛卡尔k空间轨迹。例如一个View可以通过径向轮辐的方位角和在径向三维k空间轨迹中或螺旋k空间轨迹中(可选的还具有在与螺旋扫描的平面正交的方向上的笛卡尔扫描)与径向平面正交的相位编码方向上的kz坐标来描述。
b)为了待扫描的k空间的完整采集而需要的激发的数量是Ns。
c)可以为每个激发分配一个标尺,该标尺描述了k空间中的相邻关系。在[0,...,Ns-1]中的激发索引ns相应于该标尺排序。
d)存在具有最大运动敏感性的在[0,...,Ns-1]中的激发索引ns0的杰出的激发。相应地,运动敏感性在[0,...,ns0]的范围中增加并且在[ns0,...,Ns-1]范围中又下降。
在上面描述的(Jhooti et al.的)原始工作的ky-kz排序方案中激发的数量Ns等于Ny并且提供ky作为描述了相邻关系的标尺。在k空间的对称采集中ky由此取在-Ny/2和Ny/2-1之间的范围中的值,通过以下换算获得激发索引:
ns=ky+Ny/2。
具有最大的运动敏感性的激发索引ns0=Ny/2位于值域的中心附近。
一般性的描述允许PAWS应用于灵活的ky-kz排序方案。在此理解为笛卡尔k空间轨迹,其中各个k空间行不是沿着两个笛卡尔轴的一个被采集而是例如或多或少沿着径向线,由此如上面所描述的也减小了运动敏感性。
按照本发明的磁共振设备包括基本场磁体、梯度场***、高频天线和用于控制梯度场***和高频天线的控制装置,以及用于接收由高频天线记录的测量信号、用于分析测量信号和用于建立磁共振图像并且构造为用于执行所述方法的图像计算机。
按照本发明的计算机程序包括程序装置,当计算机程序在磁共振设备的控制装置中运行时,其执行所描述的方法的所有步骤。
按照本发明的电子可读数据载体包括在其上存储的电子可读的控制信息,其构造为,在磁共振设备的控制装置中使用所述数据载体的情况下执行所描述的方法。
关于方法给出的优点和实施方式也类似地地适用于磁共振设备、计算机程序产品和电子可读数据载体。
附图说明
本发明的其它优点和细节从以下描述的实施例以及结合附图得出。所举的例子不是本发明的限制。其中,
图1示意性示出了按照本发明的磁共振设备,
图2示出了View与扇形和激发按照本发明对应的例子,如在结合图5解释的一般性PAWS算法中可以采用的,
图3示意性示出了待扫描的k空间按照本发明按照扇形和激发排序的流程图,如在结合图5解释的一般性PAWS算法中可以采用的那样,
图4示意性示出了一般性的PAWS方法的流程图,
图5举例示出了根据ky坐标的原始工作的PAWS排序方案与对于3-bin-PAWS算法按照激发索引的一般化的比较,
图6举例示出了迄今为止的2-bin PAWS与使用在此介绍的一般化的PAWS方法的、按照本发明的优化的2-bin PAWS的比较,
图7示意性示出了优化的2-bin PAWS的流程图。
具体实施方式
图1示出了(磁共振成像或核自旋断层造影设备的)磁共振装置5的示意图。在此,基本场磁体1产生时间上恒定的强磁场以用于极化或对齐检查对象U的(例如位于检查台23上为了检查或测量而被移动到磁共振装置5中的人体的待检查的部位的)检查区域中的核自旋。对于核自旋共振测量所需的基本磁场的高的均匀性在典型地为球形的测量体积M中定义,人体的待检查的部位布置在该测量体积中。为了支持均匀性要求和特别是为了消除时间上不可变的影响,在合适的位置上安装由铁磁材料组成的所谓匀场片。时间上可变的影响通过匀场线圈2和用于匀场线圈2的控制器27消除。
在基本场磁体1中采用由三个子线圈组成的圆柱形的梯度线圈***3。由放大器24-26给每个子线圈提供电流以用于在笛卡尔坐标系的各个方向上产生线性梯度场。在此,梯度场***3的第一子线圈产生x方向上的梯度Gx,第二子线圈产生y方向上的梯度Gy,并且第三子线圈产生z方向上的梯度Gz。放大器24-26包括数模转换器(DAC),该数模转换器由用于时间正确地产生梯度脉冲的序列控制装置18控制。
高频天线4位于梯度线圈***3内,所述高频天线4将由高频功率放大器给出的高频脉冲转换为用于待检查的对象或者对象的待检查区域的核的激励以及核自旋的对齐的交变磁场。高频天线4由线圈的以例如环形的、线形或矩阵形布置的形式的一个或多个HF发送线圈和一个或多个HF接收线圈组成。高频天线4的HF接收线圈也将从进动的核自旋出发的交变场、即通常由一个或多个高频脉冲和一个或多个梯度脉冲组成的脉冲序列引起的核自旋回波信号,转换为电压(测量信号),该电压经由放大器7被传输到高频***22的高频接收通道8、8'。高频***22还包括发送信道9,在该发送信道9中产生用于激励核磁共振的高频脉冲。在此,将各个高频脉冲根据由设备计算机20预先给出的脉冲序列在序列控制装置18中数字地表示为复数的序列。该数列作为实部和虚部分别经由输入端12被传输到高频***22中的数模转换器(DAC)并且从该数模转换器被传输到发送信道9。在发送信道9中将脉冲序列加调制到高频载波信号上,其基频相应于测量空间中核自旋的共振频率。通过放大器28将调制的脉冲序列传输到高频天线4的HF发送线圈。
通过发送-接收转接器6进行发送运行到接收运行的切换。高频天线4的HF发送线圈将用于激励核自旋的高频脉冲入射到测量空间M,并且通过HF接收线圈扫描所得到的回波信号。相应获得的核共振信号在高频***22的接收信道的第一解调器8'中被相位敏感地解调到中间频率,并且在模拟-数字-转换器(ADC)中被数字化。该信号还被解调到频率0。到频率0的解调和到实部和虚部的分离在第二解调器8中在数字域中数字化之后进行,所述第二解调器8通过输出端11将解调的数据输出到图像计算机17。通过图像计算机17可以从这样获得的测量数据重建MR图像、特别是在使用按照本发明的方法的条件下,所述方法包括例如借助图像计算机17计算至少一个干扰矩阵和其逆。通过设备计算机20进行测量数据、图像数据和控制程序的管理。序列控制装置18根据利用控制程序的预定值来控制各个期望的脉冲序列的产生和k空间的相应扫描、特别地根据按照本发明的方法。在此,序列控制装置18特别地控制梯度的时间正确的接通、具有定义的相位振幅的高频脉冲的发送以及核共振信号的接收。由合成器19提供用于高频***22和序列控制装置18的时间基准。通过包括用于输入可能性的输入部件诸如键盘15和/或鼠标16和用于显示可能性的显示部件诸如显示屏14的终端13来选择用于产生MR图像的相应控制程序(所述MR图像例如存储在DVD21中)以及其它用户方面的输入(诸如应当一起覆盖期望的k空间的相邻的簇的期望的数量n),和显示所产生的MR图像。
图2示出了新的k空间排序方案的例子。示出了kz-ky平面,在该平面中按照笛卡尔布置沿着kx方向的k空间行,即所谓的(通过利用不同的图案填充的圆形表示的)View。
图3示意性示出了待扫描的k空间按照扇形和激发排序的流程图。
在此,首先选择对每个激发应当采集的数量Nv个View(块301)。这一点例如通过用户输入到磁共振设备5的终端13来进行。利用对每个激发对View的选择例如可以在使用用于确定呼吸信号作为生理信号的导航器的情况下,由此自由选择呼吸信号的时间分辨率,因为在导航器之后采集具有选择的数量的View的激发。从(通过又由用户选择的分辨率等确定的)待测量的View的总数和对每个激发的View的数量Nv中由此确定为了完整说明待扫描的k空间所需的激发的数量Ns(块303)。
在k空间中的View被划分到扇形S1,S2,S3,S4,S5,S6,S7,S8(块305)。属于相同的扇形的View利用相同的图案填充。此外扇形S1,S2,S3,S4,S5,S6,S7,S8的边界通过细虚线画出。
不同扇形的数量等于每个激发的View的数量并且例如是用户定义的参数。在示出的例子中扇形的和每个激发的View的数量等于8。
每个扇形的View的数量等于激发的数量Ns。在例子中在每个扇形中的View的数量和由此激发的数量等于49。与同一个扇形对应的View具有与k空间中心的相似的距离并且位于相同的半空间中(在图2的例子中第一半空间通过ky>0或(ky=0和kz≤0))定义)。从中在非对称采集k空间(部分傅里叶)的情况下得到优点。
每个激发现在对每个扇形采集一个View。在此特定扇形的View在导航器序列之后或在激发开始之后的相同的时刻被采集。
所有View与一个激发对应,其中在每个扇形中一个View与一个特定的激发对应(块307)。一个扇形的各个View与一个特定的激发的对应相应于其在kz-ky平面中的取向进行,例如相应于其在极坐标中的方位角进行。该方案产生沿着几乎径向的扫描方向的k空间的平滑调制(由于在中断静态均衡之后的超越状态)。
图2中举例通过先后顺序的实箭头表示三个激发。在采集测量数据时与相同的扇形和不同的激发对应的View分别占据激发内部的相同(时间)位置(块309)。如在图2的例子中可以看出的,对于每个激发采集的View的顺序分别相应于其扇形,在示出的情况中从S1至S2至S3至S4至S5至S6至S7至S8。激发由此在几乎径向的k空间轨迹上在此从待扫描的k空间的右边沿到左边沿引导。
该排序方案相对于迄今为止结合PAWS技术采用的排序方案来说具有优点,即,不太容易产生由于保留的运动导致的鬼影伪影,因为该运动按照方位角被分散。此外该方案允许对每个激发的View的自由选择和由此利用导航器序列采集的生理呼吸信号的时间分辨率。此外与所谓的椭圆形扫描兼容,在所述椭圆形扫描中k空间的特别是ky-kz平面的***的角中的View,以相对小的信息内容,为了更短的测量时间而不被采集,如在示出的例子中那样。此外k空间的所介绍到扇形的划分与k空间扫描的变化的密度兼容,如在利用自动校准的并行成像中出现的那样。
一般化的PAWS描述允许,与PAWS一起采用刚才描述的排序方案,方法是,如上所述将激发通过划分到扇形的View与满足上面给出的假定c)和d)的激发索引对应。
在此提供以下工作方式:
首先为每个激发分配一个例如在[–π,π]之间的方位角
Figure BDA00003066438400162
。为此,例如可以使用在第一半空间中的激发的View的平均方位角(arctan2(ky,kz))或者在扇形中的一个的,有利地位于k空间中心的附近中的扇形诸如图2中的扇形S4的View的方位角。该方位角是合适的标尺并且由此是激发的排序标准,其描述了激发的相邻关系。人们注意到,不要求激发→标尺的对应是可逆转的。有利地将激发索引ns0置为等于采集k空间中心的激发的激发索引。图2的例子中该激发具有方位角
Figure BDA00003066438400163
并且由此位于ns[0,Ns-1]的值域的中心中。
在真实的径向轨迹中所有激发具有相同的运动敏感性。尽管如此,可以应用一般化的PAWS描述。在ns0的选择中是自由的。
利用通过上面的假定a)至d)的确定,由此可以对任意的二维笛卡尔ky-kz排序方案和k空间轨迹应用PAWS算法,方法是将(Jhooti et al.的)原始工作的ky坐标通过激发索引ns来代替。这一点例如简单地通过为具有在左k空间边沿处(原始工作中的kymin)的开始位置的簇(原始工作中的“bin”)分配激发索引ns=0作为开始位置并且为具有在右k空间边沿处(原始工作中的kymax)的开始位置的簇分配激发索引ns=Ns-1作为开始位置,并且为具有在k空间中心中(在原始工作中)的开始位置的簇分配具有最大的运动敏感性的激励索引ns=ns0作为开始位置。在图5中示出了从根据ky坐标的原始工作的排序方案(上面)到按照对于3-bin PAWS算法的激发索引的排序方案(下面)的这样的转变的例子。在纵坐标中分别画出了呼吸位置AP。
图4中示意性示出了一般化的PAWS方法的流程图。
在此,首先如在PASW中通常的那样选择数量n,该数量说明了相邻的簇的数量,其一起应当完整填充k空间,以获得完整的测量数据组,所述测量数据组不超过在测量期间呼吸位置的预定的总变化(块401)。
Ns个激发的每个如上所述对应一个激发索引
Figure BDA00003066438400161
(块403),其中这样进行对应,使得激发索引ns这样排序,使得激发相对于检查对象的运动的敏感性从激发索引ns=0到激发索引ns=ns0增加并且从激发索引ns=ns0至激发索引ns=Ns-1又降低。
以用于确定呼吸信号和由此瞬时的呼吸位置的导航器测量开始测量(块405)。
在此,将在导航器测量之后要测量的激发以在PAWS方法中通常的方式相应于利用导航器测量所测量的呼吸信号与一个簇对应(块407)。
在通过导航器测量确定,以下要测量的激发与哪个簇对应之后,根据由前面选择的簇和其相邻的簇已经采集的激发索引来确定对于待测量的激发的激发索引(块409)。如果还没有采集与待测量的激发对应的簇中的激发,则选择具有相应于分配给所述簇的开始位置的激发索引的激发(块409)。
具有选择的激发索引的相应激发作为待测量的激发被采集(块411)。
如果在一个激发的最后采集之后在一个预先给出的数量n个相邻的簇中采集了具有所有Ns个激发索引的共同的激发(询问413),则测量结束(块415)。
相应地在一般化的PAWS方法中具有开始位置0的所谓的***的簇在至k空间中心的方向上增长,方法是其选择下一个还没有采集的较大的激发索引,并且具有开始位置Ns-1的***的簇增长,方法是其选择下一个更小的还没有采集的激发索引。
中央的簇(具有开始位置ns0)从n个可能的簇组合中选择最完整的那个,即,已经包括了具有不同的激发索引的大多数激发,并且然后分别根据从具有ns≤ns0的激发中形成的集合或从还没有由簇组合采集的具有ns≥ns0的激发中形成的集合是否具有多个元素而在至更小的或更大的激发索引的方向上增长。只要任意组的n个相邻的簇包括所有待测量的激发索引[0,...,Ns-1],则测量结束(块415),因为呼吸位置的总变化由此限制到接受窗。如果还没有在n个相邻的簇的簇组合中包括具有所有Ns激发索引的激发,则在块405中继续并且分别以新的随后的激发采集新的导航器测量。
在图2的例子中k空间中的相邻的View分别属于具有类似的方位角的激发并且由此属于具有相邻的激发索引的激发。由此一般化的PAWS算法负责,相邻的View在类似的呼吸状态期间被采集。在图2的例子中采集k空间中心的激发占据方位角
Figure BDA00003066438400171
并且由此位于激发索引的值域的中心。由此可以假定,在原始工作中获得的结果(关于效率和保留的运动倾向性)也可以直接传输到按照本发明的灵活的ky-kz排序方案。利用所描述的方法可以将PAWS原理应用于任意的笛卡尔ky-kz排序方案和任意的非笛卡尔k空间轨迹。
在采用用于采集运动的导航器的情况下通过执行导航器序列中断磁化的静态均衡。在按照Jhooti et al的原始工作的排序方案中该中断分别在采集Nz个TR间隔之后进行,其中Nz是在第二相位编码方向上的相位编码步骤的数量。呼吸信号的时间分辨率由此直接与沿着第二相位编码方向的成像的序列的空间分辨率相关联。然而利用导航器测量的呼吸信号仅一个有限的时间是有效的,该时间相对于呼吸间隔来说是短的。即,在Jhooti et al的原始工作中使用的排序方案固有地将最大的分辨率限制到两个笛卡尔方向中的一个。在采用介绍的具有排序方案的一般化的PAWS算法时,如上特别是参考图2所述,不存在这样的限制,因为每个激发的View的数量和由此呼吸信号的时间分辨率可以自由选择。该优点是特别重要的,因为呼吸门控的测量的目的通常是,避免由于有限的取决于患者的屏气时间而固有的、屏气测量的分辨率限制。
导航器信号的有限的时间有效性的问题部分地通过采用沿着Kz方向的一维中央的排序方案来避免。这样的排序方案在k空间中心中开始并且这样交替地采集具有kz的正的和负的值的View,使得相位编码梯度的绝对矩(也就是kz=0的距离)在分别两个相位编码步骤之后连续增长。但是该中央的排序方案具有缺陷,即,会导致在TR间隔之间通过大的相位编码跳变由于涡流引起的伪影。
替换通常的n-bin PAWS方法,以下介绍优化的2-bin PAWS方法。
图6举例示出了迄今为止的2-bin PAWS与按照本发明的优化的2-bin PAWS的比较,其中上面示出了迄今为止的2-bin PAWS并且下面示出了新的2-binPAWS。
如在(Jhooti et al.的)原始工作中描述的,在最初的2-bin PAWS变形中相邻的簇的开始位置在左和右k空间边沿之间交替。在图6中上面具有偶数索引的簇对应右开始位置并且具有奇数索引的簇对应左开始位置。这一点在此处刚刚描述的、一般化的显示中相应于在ns=0和ns=Ns-1之间的交替。相应地,具有开始位置ns=0的簇通过选择还没有被该簇采集的最小的激发索引而增长。具有开始位置ns=Ns-1的簇通过选择还没有由该簇采集的最大的激发索引来增长。以下始终提到激发索引ns,即使可以使用按照ky坐标的初始排序。与每个簇对应的呼吸位置区域的宽度通常相应于半个接受窗(Af)。只要两个任意的相邻的簇一起已经采集所有激发,测量结束。这一点在图6的例子中对于簇c4和c5是这样。图形地来说,“从左增长的”簇(开始位置ns=0)和两个相邻的“从右增长的”簇之一(开始位置ns=Ns-1)相遇,从而两个簇一起张开共同的值域[0,...,Ns-1]。只要这两个簇通常包括类似相同数量的激发,则在此通常在k空间中心附近的运动敏感的区域中得到一个簇边界(通过图6中的虚线椭圆标记)。
在这里介绍的、按照本发明的最佳的2-bin PAWS实施也是两个簇类型交替。一个簇类型具有带有最大的运动敏感性ns=0作为开始位置的激发并且在以下称为中央的簇。另一个簇类型没有明确的开始位置并且在以下称为***的簇。在图6的例子中具有奇数索引的簇是中央的簇并且具有偶数索引的簇是***的簇。***的簇的开始位置要么是ns=0要么是ns=Ns-1,或者说要么是ky坐标中的右k空间边沿要么是左k空间边沿,其中实际的开始位置在运行时间才决定。***的簇独立于其开始位置地增长,或者是从更大的还不属于该簇的激发索引在ns0或ky=0的方向上降序地、或者从最小的还不属于该簇的激发索引在ns0或ky=0的方向上升序地。关于哪个增长方向当前是优选的判断又在运行时间进行。这一点示意性在图7中的流程图中示出。
如果如上面已经描述的借助导航器测量测量呼吸位置(块701,相应于图4中的块405),其按照在PAWS中通常的对应(图4中的块407)位于***的簇cn的呼吸位置区域中,则首先询问,簇组合cn-c(n-1)或cn-c(n+1)是否接近完整(块705和707)。对此,首先对在簇c(n-1)和c(n+1)中已经采集的激发计数,并且存储这些各自的数量M-及M+(块705)。相邻的簇cx(***的簇cn与所述相邻的簇一起共同最接近完整)相应于该数量M-及M+来选择,其中当M-大于M+时选择簇c(n-1),并且当M+大于M-时选择簇c(n+1)。簇如通常的那样相应于其呼吸位置区域来标记(簇cn相应于第n个呼吸位置)。相应地c(n-1)和c(n+1)是中央的簇并且是***的簇cn的两个最相邻的簇。以这种方式当簇组合cn-c(n-1)更接近完整时选择簇c(n-1),否则选择簇c(n+1)。
接下来对还没有由两个簇(cn和选择的cx))采集的、具有在范围[0,...,ns]中的索引的激发Mlow的数量,和还没有由两个簇采集的、具有在范围[ns,...,Ns-1]中的索引的Mhigh的数量进行计数(块709)。如果与最后测量的呼吸位置对应的簇cn如在假定的情况中那样是***的簇(询问711中的“y”),则当Mlow大于Mhigh(询问713中的“y”)时,***的簇cn从其最小的还没有采集的激发索引在方向ns0或ky=0上增长,方法是采集最小的、还没有采集的激发索引(块715),否则从其最大的还没有采集的激发索引在方向ns0或ky=0上增长(询问713中的“n”),方法是采集最大的、还没有采集的激发索引(块717)。
也就是***的簇一般地张开两个关联的索引区域。一个在最小的激发索引ns=0(或左k空间边沿)的情况下开始并且在更大的激发索引的方向上增长。另一个在最大的激发索引ns=Ns-1(或右k空间边沿)的情况下开始并且在更小的激发索引的方向上增长。或者也可以说,对于***的簇的索引区域在区域边界处周期性或循环地继续。
当最后测量的呼吸位置位于中央的簇cn的呼吸位置区域中时,判断过程类似进行:
如刚才所述,询问簇组合cn-c(n-1)和cn-c(n+1)中的哪个更接近完整,并且选择该簇组合(块705和707)。接下来对还没有由两个簇采集的、具有在范围[0,...,ns]中的索引的激发Mlow的数量,以及还没有由两个簇采集的、具有在范围[ns,...,Ns-1]中的索引的Mhigh的数量进行计数(块709)。如果与最后测量的呼吸位置对应的簇cn如现在假定的那样是中央的簇(询问711中的“n”),则当Mlow大于Mhigh时中央的簇cn从其最小的已经采集的激发索引在方向ns=0或kmin上增长(询问719中的“n”),方法是采集最大的、小于最小的已经采集的激发索引的那个还没有被采集的激发索引(块723)。否则(Mhigh大于Mlow)(询问719中的“y”),簇cn从其最大的已经采集的激发索引在方向Ns-1或kmax上增长,方法是采集最小的、大于最大已经采集的激发索引的那个激发索引(块721)。
在任何情况下在块715、717、721和723中在一个激发的采集之后在询问725中检验,在选择的簇组合中所有Ns个期望的激发索引是否已经被采集。如果是(询问725中的“y”)则测量是完整的(块727)并且可以结束;如果否(询问725中的“n”)则以新的导航器测量继续。由此相对于PAWS的原始版本中断标准不变:只要两个相邻的簇(2-bin)的任意组已经采集了所有要测量的激发索引,则结束测量,因为呼吸位置的总变化限制到接受窗。
在图7中的总结了刚才描述的算法的流程图中,使用集合论的通常符号:
{.}...表示集合。
{xs|…}...表示“对于…成立”的所有激发索引xs的集合。
∈...表示“是…的元素”。
#{.}...表示集合的元素的数量。
∧...对于“与”的逻辑符号。
∨...对于“或”的逻辑符号。
PAWS的该最佳的2-bin版本统一了原始的2-bin PAWS版本的高效率与原始的3-bin变形的降低的伪影可能性。新的算法主动地将片段边界从运动特别敏感的k空间中心在k空间***的方向上移开。
在不同的PAWS变形的以下比较中假定,接受窗的总宽度是给定的。在n-bin变形的情况下该接受窗由n个相继跟随的簇的呼吸位置区域张开。例如为每个簇分配一个呼吸位置区域,其宽度相应于接受窗的1/n。这一点与上面已经引用的Jhooti et al.的MRM文章的附录A的比较稍有不同,其中一个簇的呼吸位置区域的宽度设置为等于导航器分辨率。在后一种工作方式中接受窗的总宽度是导航器分辨率的n倍并且随着bin的数量增加。这一点使得不同的bin变形的公正比较变得困难。
此处刚刚描述的新的2-bin变形的效率在如下意义上是最佳的,即,只要在由两个相邻的簇覆盖的区域中的呼吸位置被测量了Ns次则采集了所有Ns个激发并且由此可以结束测量。该特征由新的2-bin变形与原始的2-bin变形分享并且比原始的3-bin变形和上面引用的Nuval et al.的文章中的3-bin变形更突出。
与原始的2-bin变形不同,明显降低了在运动特别敏感的k空间中心附近形成簇边界的可能性。
这一点在图6中解释,在图6中如通常的那样每条线相应于一个簇cn。它们在垂直方向上相应于其呼吸位置区域布置。在水平方向上画出相位编码索引ky或在一般的图中画出激发索引ns。灰色阴影条示出由一个簇采集的ky行或激发。上面是已经在前面引用的Jhooti et al.的MRM文章的图11b,其示出了在测量的结束时原始的2-bin变形的相位编码线的选择。下面示出了新的2-bin变形的相应图示。注意到,多频繁地测量一个特定的呼吸位置的次数在两个图(Plot)中是相同的。在上面描述的原始PAWS方法中簇边界靠近k空间中心。在下面示出的新的变形中其明显在k空间***的方向上移动。在图6中分别利用虚线画出的椭圆来突出所述簇边界。
当分配给两个最终的簇的呼吸位置大约相同频繁地被测量并且中央的、运动特别敏感的具有激发索引ns0或ky=0的激发精确位于索引区域的中心时,该问题在原始版本中是最大的。在这种情况中簇边界精确与k空间中心(ky=0)一致。按照本发明的新的版本最佳地处理该特别重要的情况:簇边界处于值域的大约+25%和+75%并且由此与运动敏感的k空间中心距离最大。
在k空间的非对称扫描时在此介绍的新的2-bin变形也在如下意义上最佳地反应,即,在分摊到中央的簇的扫描的数量给定的情况下,片段边界与具有激发索引ns0或ky=0的中央的运动特别敏感的激发距离最大。由此新的版本实际地总是比原始的2-bin版本更好地反应。对此的原因是,激发的围绕中央的激发ns0的对称分布主动地嵌入到算法的判断过程中。
当且仅当分摊到最终的中央的簇的扫描的数量相对于分摊到最终的***的簇的扫描的数量来说是小的时,才可以在k空间中心附近形成簇边界。在该情况中该边界也如在原始版本中那样可以靠近k空间中心。但是该情况在接受窗的合适选择和呼吸位置在最可能的呼吸位置的附近的静态分布的情况下是极其不可能的并且在我们的利用新的2-bin变形的大量测量中没有观察到。通过类似于上面已经引用的Nuval et al.的文章中的修改b)的扩展,该情况甚至可以完全被避免:这样收紧中断标准,即,中央的最终的簇要么必须已经采集了所有激发Ns的最小百分比,要么***的最终的簇必须已经采集了所有Ns个激发。注意到,围绕k空间中心的对称分布对于新的2-bin算法来说是固有的并且不必要求(与现有技术中的3-bin变形不同)。
还要提到,在新的2-bin算法的(以及原始版本的)“没有呼吸”的边界情况中被最佳处理:所有激发由唯一一个簇采集,也就是不形成簇边界,而不管该簇是中央的还是***的。

Claims (9)

1.一种用于借助三维磁共振技术确定呼吸的检查对象的磁共振图像的方法,
其中,将检查对象的呼吸位置划分为多个预定义的簇,所述簇分别覆盖检查对象的呼吸位置的预定义的呼吸区域,
其中,所述测量数据组通过多个激发采集,所述激发分别采集预先给定数量的具有分别预定义的三维k空间轨迹的k空间行,
其中,不同的激发采集不同的k空间行,从而扫描预定义的k空间区域,
其中,所述激发分别对应一个激发索引(ns),该激发索引说明了各自的激发相对于检查对象的运动的敏感性,
其中,激发索引与激发的对应是这样的,即,紧挨的激发索引对应于这样两个相对于运动具有递增改变的敏感性的激发,
其中,该方法包括:
A)借助导航器技术测量检查对象的呼吸位置,
B)选择特定的簇,测量的呼吸位置位于其覆盖的呼吸区域,
C)选择待采集的激发作为具有对应的激发索引的那个激发,该激发索引与最后对于特定的簇采集的激发的激发索引紧挨,
D)采集待采集的激发并且将该激发与特定的簇对应,
-一直执行步骤A)-D),直到覆盖了检查对象的相邻的呼吸区域的至少两个紧挨的簇一起与测试数据组的所有激发对应,
-从与两个紧挨的簇对应的激发中确定MR图像。
2.根据权利要求1所述的方法,
其中,对于覆盖了检查对象的相邻的呼吸区域的两个紧挨的簇的第一个首先选择并采集具有与相对于运动的大的敏感性相应的对应的激发索引的这些激发,
其中,对于两个紧挨的簇的第二个首先选择和采集具有与相对于运动的小的敏感性相应的对应的激发索引的这些激发。
3.根据权利要求2所述的方法,
其中,对于两个紧挨的簇的第一个来说待采集的激发的选择包括:
-选择至此对于两个紧挨的簇的第一个还没有被测量的并且具有相对于该至此还没有被测量的激发的运动的最大的敏感性的那个激发,
其中,对于两个紧挨的簇的第二个来说待采集的激发的选择包括:
-选择至此对于两个紧挨的簇的第二个来说还没有被测量的并且具有相对于该至此还没有测量的激发的运动的最小的敏感性的那个激发。
4.根据上述权利要求中任一项所述的方法,
其中,待采集的激发的选择还这样进行,即,尽可能均匀地沿着不同的方向和/或在k空间区域的不同侧扫描k空间区域。
5.根据上述权利要求中任一项所述的方法,
其中,k空间行沿着相位编码梯度场方向布置,
其中,多个激发的k空间轨迹分别主要径向地在垂直于相位编码梯度场方向取向的平面中布置。
6.根据权利要求5所述的方法,
其中,k空间区域划分为扇形,所述扇形关于其与k空间中心的距离定义,
其中,对于多个激发中的每一个来说分别一个k空间行位于每个扇形中。
7.根据权利要求6所述的方法,
其中,每个激发的k空间轨迹按照顺序扫描不同的k空间行,所述顺序相应于k空间区域中的、不同的k空间行所位于其中的扇形的相邻关系。
8.一种磁共振设备(5),其构造为用于借助三维磁共振技术确定呼吸的检查对象的磁共振图像,
其中,检查对象的呼吸位置划分为多个预定义的簇,所述簇分别覆盖检查对象的呼吸位置的预定义的呼吸区域,
其中,所述测量数据组通过多个激发采集,所述激发分别采集预先给定数量的具有分别预定义的三维k空间轨迹的k空间行,
其中,不同的激发采集不同的k空间行,从而扫描预定义的k空间区域,
其中,所述激发分别对应一个激发索引,该激发索引说明了各自的激发相对于检查对象的运动的敏感性,
其中,激发索引与激发的对应是这样的,即,紧挨的激发索引对应于这样两个相对于运动具有递增改变的敏感性的激发,
其中,所述磁共振设备(5)包括序列控制器(18),该序列控制器构造为用于执行以下步骤:
A)借助导航器技术测量检查对象的呼吸位置,
B)选择特定的簇,测量的呼吸位置位于其覆盖的呼吸区域,
C)选择待采集的激发作为具有对应的激发索引的那个激发,该激发索引与最后对于特定的簇采集的激发的激发索引紧挨,
D)采集待采集的激发并且将该激发与特定的簇对应,
-一直执行步骤A)-D),直到覆盖了检查对象的相邻的呼吸区域的至少两个紧挨的簇一起与测试数据组的所有激发对应,
其中,所述磁共振设备(5)还包括图像计算机(17),其构造为用于执行以下步骤:
-从与两个紧挨的簇对应的激发中确定MR图像。
9.根据权利要求8所述的磁共振设备(5),其还构造为用于执行按照权利要求1-7中任一项所述的方法。
CN201310134935.4A 2012-04-20 2013-04-18 借助磁共振技术采集呼吸的检查对象的测量数据组的方法 Expired - Fee Related CN103371820B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
DE1020122065504 2012-04-20
DE102012206550.4 2012-04-20
DE102012206550A DE102012206550A1 (de) 2012-04-20 2012-04-20 Verfahren zur Akquisition eines Messdatensatzes eines atmenden Untersuchungsobjekts mittels Magnetresonanztechnik, Magnetresonanzanlage, Computerprogramm sowie elektronisch lesbarer Datenträger

Publications (2)

Publication Number Publication Date
CN103371820A true CN103371820A (zh) 2013-10-30
CN103371820B CN103371820B (zh) 2016-02-03

Family

ID=49290179

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310134935.4A Expired - Fee Related CN103371820B (zh) 2012-04-20 2013-04-18 借助磁共振技术采集呼吸的检查对象的测量数据组的方法

Country Status (3)

Country Link
US (1) US9261575B2 (zh)
CN (1) CN103371820B (zh)
DE (1) DE102012206550A1 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105182264A (zh) * 2014-06-05 2015-12-23 西门子公司 磁共振技术中参数图的产生
CN108567422A (zh) * 2018-03-02 2018-09-25 上海东软医疗科技有限公司 一种血管成像方法和装置
CN110706201A (zh) * 2019-09-04 2020-01-17 上海东软医疗科技有限公司 磁共振成像方法、装置、设备及***
CN112315449A (zh) * 2020-11-24 2021-02-05 上海联影医疗科技股份有限公司 磁共振成像方法、***、电子装置和存储介质

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8970217B1 (en) 2010-04-14 2015-03-03 Hypres, Inc. System and method for noise reduction in magnetic resonance imaging
JP6440980B2 (ja) * 2013-07-10 2018-12-19 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
DE102015107347A1 (de) * 2015-05-11 2016-11-17 Universitätsspital Basel Ein magnetresonanztomographie-verfahren mit asymmetrischer radialer akquisition von k-raum-daten
US10267883B2 (en) * 2015-12-11 2019-04-23 Siemens Healthcare Gmbh System and method for motion resolved MRI

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7039451B1 (en) * 1999-08-20 2006-05-02 Imperial College Innovations Limited Phase ordering with automatic window selection (PAWS) for motion resistant MRI
CN101144853A (zh) * 2006-09-13 2008-03-19 西门子公司 用于三维磁共振成像的测量序列以及磁共振设备
CN102078192A (zh) * 2009-11-27 2011-06-01 西门子公司 借助磁共振技术采集呼吸着的检查对象的测量数据的方法
CN102078196A (zh) * 2009-11-27 2011-06-01 西门子公司 借助磁共振采集相应于对象的呼吸运动的信号数据的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102011007850B4 (de) * 2011-04-21 2015-05-13 Siemens Aktiengesellschaft Verfahren zur getriggerten Akquisition von Messdaten mittels eines Magnetresonanzgeräts unter Unterdrückung von unerwünschten Signalen, Magnetresonanzgerät, Computerprogramm und elektronisch lesbarer Datenträger

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7039451B1 (en) * 1999-08-20 2006-05-02 Imperial College Innovations Limited Phase ordering with automatic window selection (PAWS) for motion resistant MRI
CN101144853A (zh) * 2006-09-13 2008-03-19 西门子公司 用于三维磁共振成像的测量序列以及磁共振设备
CN102078192A (zh) * 2009-11-27 2011-06-01 西门子公司 借助磁共振技术采集呼吸着的检查对象的测量数据的方法
CN102078196A (zh) * 2009-11-27 2011-06-01 西门子公司 借助磁共振采集相应于对象的呼吸运动的信号数据的方法

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105182264A (zh) * 2014-06-05 2015-12-23 西门子公司 磁共振技术中参数图的产生
CN105182264B (zh) * 2014-06-05 2018-11-16 西门子公司 磁共振技术中参数图的产生
CN108567422A (zh) * 2018-03-02 2018-09-25 上海东软医疗科技有限公司 一种血管成像方法和装置
CN108567422B (zh) * 2018-03-02 2021-03-05 上海东软医疗科技有限公司 一种血管成像方法和装置
CN110706201A (zh) * 2019-09-04 2020-01-17 上海东软医疗科技有限公司 磁共振成像方法、装置、设备及***
CN110706201B (zh) * 2019-09-04 2022-08-09 上海东软医疗科技有限公司 磁共振成像方法、装置、设备及***
CN112315449A (zh) * 2020-11-24 2021-02-05 上海联影医疗科技股份有限公司 磁共振成像方法、***、电子装置和存储介质
CN112315449B (zh) * 2020-11-24 2023-08-08 上海联影医疗科技股份有限公司 磁共振成像方法、***、电子装置和存储介质

Also Published As

Publication number Publication date
US20130281828A1 (en) 2013-10-24
DE102012206550A1 (de) 2013-10-24
US9261575B2 (en) 2016-02-16
CN103371820B (zh) 2016-02-03

Similar Documents

Publication Publication Date Title
CN103371820A (zh) 借助磁共振技术采集呼吸的检查对象的测量数据组的方法
CN103371821A (zh) 借助磁共振技术采集呼吸的检查对象的测量数据组的方法
CN103376434A (zh) 借助磁共振技术采集呼吸的检查对象的测量数据组的方法
CN104101853B (zh) 用于确定共振频率偏差的方法和磁共振设备
US6516210B1 (en) Signal analysis for navigated magnetic resonance imaging
US10634753B2 (en) MR imaging with motion detection
CN102736046B (zh) 产生扩散信息的方法和相应的磁共振设备
CN1653348B (zh) 产生时间连续的磁共振图像的磁共振成像方法
CN107750338A (zh) 使用星形叠层采集的mr成像
CN102870000A (zh) 使用单个和多个通道接收器线圈的用于同时多切片磁共振成象的方法
CN103371823A (zh) 借助磁共振技术采集呼吸的检查对象的测量数据组的方法
CN101843487A (zh) 借助并行采集技术建立图像的方法、磁共振设备和程序
CN102743172A (zh) 借助磁共振设备抑制不期望的信号触发地采集测量数据
CN108693492A (zh) 用于相位循环的磁共振指纹(phc-mrf)的***和方法
CN103091655B (zh) 用于建立磁共振图像的方法以及相应的磁共振设备
JP2020523075A (ja) 可変コントラストのスタック・オブ・スター収集を使用したmrイメージング
CN106772157A (zh) 一种磁共振成像的方法和装置
CN106796274A (zh) 具有伪迹抑制的propeller‑mr成像
CN107548464A (zh) 自旋回波mr成像
US11852705B2 (en) MR imaging using a 3D radial or spiral acquisition with soft motion gating
JP2021512671A (ja) 固有運動補正を有するスタック・オブ・スター収集を用いたmr撮像
US5766128A (en) Respiratory motion compensation using segmented k-space magnetic resonance imaging
CN102565736A (zh) 确定层关于相对其运动的部位的位置的方法和磁共振设备
CN116724246A (zh) Mr***摄影

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160203

Termination date: 20200418

CF01 Termination of patent right due to non-payment of annual fee