CN101750596B - 磁共振测量中用于确定预定信号振幅的方法 - Google Patents

磁共振测量中用于确定预定信号振幅的方法 Download PDF

Info

Publication number
CN101750596B
CN101750596B CN200910254208.5A CN200910254208A CN101750596B CN 101750596 B CN101750596 B CN 101750596B CN 200910254208 A CN200910254208 A CN 200910254208A CN 101750596 B CN101750596 B CN 101750596B
Authority
CN
China
Prior art keywords
target
frequency
pulse
magnetization
phase
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.)
Expired - Fee Related
Application number
CN200910254208.5A
Other languages
English (en)
Other versions
CN101750596A (zh
Inventor
汉斯-彼得·福茨
弗朗兹·施米特
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 CN101750596A publication Critical patent/CN101750596A/zh
Application granted granted Critical
Publication of CN101750596B publication Critical patent/CN101750596B/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/58Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
    • G01R33/583Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency
    • 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/24Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/246Spatial mapping of the RF magnetic field B1
    • 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/543Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription

Abstract

本发明涉及一种用于在磁共振测量中确定检查对象(12)的预定信号振幅的方法,其中在脉冲序列中将具有小于检查对象的T1时间的重复时间TR的多个高频脉冲以脉冲序列的形式入射到检查对象中,具有以下步骤:在入射各高频脉冲之后的预定时刻确定对于脉冲序列的基本上所有高频脉冲的目标磁化;根据相应的目标磁化和现有的磁化针对检查对象的不同区域对于基本上所有高频脉冲确定目标翻转角和目标相位;对于基本上所有高频脉冲确定振幅和相位的时间变化,以在入射高频脉冲之后产生相应的目标磁化;以及,入射具有相应确定的振幅和相位关系的高频脉冲。

Description

磁共振测量中用于确定预定信号振幅的方法
技术领域
本发明涉及一种用于在磁共振测量中确定检查对象的预定信号振幅的方法以及为此的磁共振设备,在该磁共振测量中在脉冲序列中将多个高频脉冲按照脉冲顺序入射到检查对象中。本发明特别应用于降低MR拍摄中的空间图像非均匀性,该图像非均匀性的原因是检查对象中高频场分布的空间变化。
背景技术
特别在全身成像时在具有高场强例如3特斯拉的MR设备的情况下在图像中出现伪影,迄今为止这种伪影妨碍了这种检查的广泛使用。随着使用的场强B0的增加,图像伪影也增强。在更高的场强下图像伪影也更多更强地出现,并且在这种场强下也越来越影响头部的成像。
由于不均匀的B1场,即入射的高频场,在MR成像和光谱学中的伪影和不一致长久以来在磁共振技术中是公知的。利用常规的方法不可能直接影响HF场的B1均匀性,从而常规的方法很大程度上基于对B1非均匀性尽可能的不敏感性。例如使用组合(所谓的Composite)脉冲和绝热脉冲,然而该脉冲关于可达到的翻转角、在用于层选择的使用中的相位关系、脉冲时间和在将HF功率入射到身体中时的吸收率,具有受限制的可应用性。由于这个原因,这样的脉冲通常被用于磁化的准备,但是在成像序列中在磁化的激励和再聚焦中不能采用。
此外公知其本身对于翻转角变化不太敏感的成像序列以及磁化准备,以减少随后的成像序列的敏感性(参见参考文献(1))。此外公知的是,通过HF脉冲和梯度脉冲同时对自旋系统的影响实现所产生的横向磁化的空间调制。二维及三维脉冲的可达到的均匀性原则上是不受限制的,然而,该调制导致非常长的脉冲时间。利用多个HF信道并行发送的可能性,可以缩短该脉冲时间(参见参考文献(2)和(3))。但是,可达到的脉冲时间比其在通常的成像序列中代替迄今为止通常的层选择性的或者非选择性的脉冲,还是太长。同样公知用于补偿单个激励的B1场非均匀性的方法,这些方法使用少数几个子轨迹,这些子轨迹在成像中仅扫描少数几个k空间点(参见参考文献(4)和(5))。
通过利用多个空间上分离的HF发送线圈或者说HF信道在时间上同时入射HF脉冲可以直接影响HF场。通过在多个并行运行的HF发送器中匹配相位值和幅度值可以在空间上调制所产生的HF场,即,也以RF补偿(Shimming)公知的方法。在此,可到达的均匀性主要受到可用的并行发送信道数量的限制。并行发送的方法具有如下优点:该方法可以直接应用于所有通常的成像方法,而无需在时间上改变成像序列。
在MR成像中几乎仅使用周期性的成像序列,其中以特定的时间间隔入射具有预定翻转角和相位角的HF脉冲。在现代的快速成像序列中,HF脉冲如此快速地相继跟随,以至于通过HF脉冲产生的横向和纵向磁化直到随后的HF脉冲还没有重新弛豫。在其它的成像序列、例如多自旋回波序列或特定的梯度回波序列中,HF脉冲以短间隔相继跟随,以产生多个不同相位编码的MR信号。在这种情况下自旋系统的演变(Evolution)非常复杂并且利用布洛赫(Bloch)公式可能在少数几个脉冲之后数值上只能非常麻烦地计算。扩展的相位图形算法(Extended Phase Graph Algorithm EPG)是对于置于硬脉冲的串中的自旋演变的布洛赫(Bloch)公式的k空间类似描述(例如参见参考文献(7)和(8))。此处借助不同的去相位状态来描述自旋系统,并且可能的状态的数量以HF脉冲的三倍数量增加。在一个回波、即实际的MR信号中,根据序列仅读出一个状态。该状态的总体(Population)、即回波的信号强度,从多个可能的回波路径馈入,这些回波路径在HF序列的期间根据相应的翻转角和应用的脉冲被占满。可以唯一地从脉冲的翻转角和相位来确定形成回波的磁化。在考虑弛豫时还必须已知弛豫时间。
反问题、即翻转角和相位的计算(该计算导致具有预先定义的振幅的回波串),是不唯一的。同样也未知对于整个回波串确定一个非唯一解的一般方法。然而,存在所谓的预测方法(look-Ahead-Verfahren),该方法从磁化状态出发计算所需的翻转角,以利用一个或少数几个脉冲得到预定义的信号振幅。已表明,能够准备(praeparieren)取决于翻转角的平衡状态。利用具有开始值α(m)和结束值α(n+m)的m个翻转角的连续序列可以在属于开始翻转角以及结束翻转角的平衡状态之间交替,而不产生强信号波动(参见参考文献(8)-(12))。
迄今为止使用利用可变的翻转角产生预先定义的磁化的可能性,以例如稳定回波串中的信号振幅和避免信号波动(参见参考文献(14))。此外使用该可能性,以降低身体中的能量沉积、即信号吸收或者SAR(信号吸收率)。在其中记录了外部k空间区域的信号的回波串的片段中,使用低翻转角,而逐渐增大翻转角,以产生对于平均k空间区域的信号(参见参考文献(10)、(12)、(13))。同样,使用该可能性,以延缓沿着回波串的信号衰减,并使用于三维快速自旋回波记录的较长的回波串是可能的(参见参考文献(8)、(9)、(11)、(12))。
此外还可以降低快速自旋回波序列对B1场非均匀性的灵敏性。在此,准备磁化状态,从该磁化状态出发随后产生的回波振幅尽可能少地依赖于应用的脉冲使用的翻转角(参见参考文献(1))。
但在MR拍摄中特别是在较高场强的情况下,降低B1非均匀性引起的伪影的问题依然存在。
发明内容
因此本发明要解决的技术问题是,进一步降低B1引起的磁场非均匀性和由此产生的伪影。
按照本发明的一个方面,提供一种用于在MR测量中确定检查对象的预定信号振幅的方法,其中在脉冲序列的情况下将多个HF脉冲入射到检查对象中。在此,在入射各个HF脉冲之后的预定时刻对于脉冲序列的基本上所有HF脉冲确定目标磁化。在下一步骤中,根据各个目标磁化对于基本上所有HF脉冲对于检查对象的不同空间区域确定目标翻转角和目标相位,该目标磁化是产生各个HF脉冲之后应该形成的并且取决于在不同区域中存在的磁化。此外为了实现目标翻转角和目标相位,对于基本上所有HF脉冲确定振幅关系和相位关系在时间上的变化,并且将具有各个确定的振幅关系和相位关系的HF脉冲入射到检查对象中。
本发明基于以下思想,通过超过脉冲序列相应于目标磁化在空间上改变HF脉冲的空间上非均匀的翻转角特性和相位特性,产生磁化的空间调制。这意味着,根据位置每个磁化经历翻转角和相位的不同的序列,然而考虑,对于给定的时刻,在脉冲序列中能够同时读出预先定义的信号振幅,即,目标磁化。
然后,这些目标翻转角和目标相位导致预先定义的磁化,即目标磁化。对于空间目标翻转角和目标相位的确定,可以对于每个图像点以及检查对象中的每个区域使用从现有技术中公开的方法(例如上面提到的预测方法),以计算产生预先定义的磁化的翻转角和相位。优选利用多个HF发送信道同时入射HF脉冲,在此对于所有发送信道确定振幅和相位关系,以实现检查对象的不同区域中空间上预先给定的目标翻转角。优选对于脉冲序列的每个HF脉冲确定对于所有发送信道的振幅关系和相位关系。这意味着,对于每个HF信号和对于脉冲序列的每个HF脉冲,在入射前确定振幅和相位关系。
按照本发明的优选实施方式,在入射具有n个脉冲的脉冲序列中的第i个HF脉冲之后(i=1至n),计算磁化,并且在第i+1个HF脉冲时对于检查对象的不同区域在确定目标翻转角和目标相位时将该磁化作为开始磁化使用。这意味着,计算所有i个脉冲对磁化的作用(i小于i+1),并且作为现有的磁化被使用,以在步骤i+1中为了达到目标磁化计算目标翻转角和目标相位的空间图。优选对于开始磁化的确定,考虑在第i个和第i+1个HF脉冲之间磁化的演变。
这例如可以通过数据记录和同时接通梯度而形成。考虑存在的磁化,因为检查的组织成分的T2时间大于或等于重复时间TR,从而必须考虑从前面的脉冲还存在的磁化。检查对象可以有具有不同T2时间的组织。然而T2时间这样位于一个范围中,使得在重复时间TR之后下一个HF脉冲跟随之前,至少几个组织成分的横向磁化还没有衰减。然后在计算目标翻转角和目标相位时考虑该存在的磁化。
值得期待的是,对于脉冲序列的所有HF脉冲来说目标磁化基本上相同。当达到,在脉冲序列期间在不同的相位编码步骤中的信号探测时分别探测到相同的磁化,总的来说这意味着良好的图像质量。例如可以通过每个HF脉冲之后的相同的目标磁化来实现保持相同的磁化。
另一方面值得期待的是,特别是对一个或多个时刻,例如在读出k空间中心时,达到预定的磁化。在这种情况下具有优势的是,在不同的相位编码步骤中各个目标磁化与该目标一致。即,对于各个HF脉冲的目标磁化可以完全不同,以实现对各个时刻与目标磁化的精确一致。
同样还可以,在确定目标翻转角和目标相位时除了目标磁化之外还考虑其它设计目标。特别具有优势的是,这样确定目标翻转角和目标相位,使得超过图像拍摄最小特殊吸收率(SAR)。一般地,在对于各个HF脉冲确定目标翻转角和目标相位时可以考虑所有从非位置选择性的翻转角设计已知的并且在第3页第2和3段提到的目标(信号稳定、SAR最小化、信号衰减和对比度的修改、磁化准备等等)。
此外还可以,在确定目标翻转角和目标相位时必须考虑由脉冲序列引起的边界条件。根据成像序列的不同,不同的边界条件可以适用于目标翻转角和相位的计算。在CMPG(Car Purcell Meiboom Gill)序列的情况下例如相位不是沿着回波串可自由选择的。
此外,HF脉冲的局部相位可以经过脉冲序列改变,并且具体来说与位置有关。由此改变的发送相位不再全局地通过匹配的接收相位补偿。然而发送相位关于线圈的灵敏性图和对于每个位置的脉冲设计是已知的。由此在图像再现中可以回顾地考虑数据的相位校正。
为了实现HF序列的从HF脉冲到HF脉冲的不同的激励特性和再聚焦特性,优选通过从不同的HF发送信道同时入射HF脉冲来实现在检查对象的不同区域中的目标翻转角和目标相位。
为了实现检查对象的不同区域中的目标翻转角和目标相位,还单独地匹配在HF脉冲期间接通的磁场梯度。
这意味着,对于每个HF脉冲匹配具有振幅和相位的HF脉冲变化以及梯度脉冲变化,以实现在不同身体区域中期望的磁化。为了从目标翻转角和目标相位计算HF脉冲和梯度脉冲的单独匹配的时间变化,可以使用从文献中公知的方法(参见参考文献(2)-(5))。同样,在该步骤中具有优势的是,考虑其它设计目标特别是各个HF脉冲SAR的最小化。
本发明还涉及一种用于在MR测量中确定检查对象的预定信号振幅的磁共振设备,其中设置HF发送单元和用于确定目标磁化的单元,该单元在入射各个HF脉冲之后的预定时刻对于脉冲序列的基本上所有HF脉冲确定目标磁化。计算单元根据各个目标磁化对于检查对象的不同区域、例如对于后面的MR图像的每个图像点,确定对于基本上所有HF脉冲的目标翻转角和目标相位。HF控制单元对于基本上所有HF脉冲分别确定振幅和相位关系,用于在入射各个HF脉冲之后产生相应的目标磁化。然后HF发送单元入射具有相应确定的振幅和相位关系的各个HF脉冲。优选HF发送单元具有多个发送信道,并且HF控制单元在所有HF脉冲的情况下对所有HF发送信道确定振幅和相位关系,以实现各个HF脉冲的目标翻转角和目标相位。
附图说明
以下参考附图详细解释本发明。在此
图1示意性示出了按照本发明的MR设备,
图2示出了用于产生具有预定的信号振幅的MR图像的流程图,
图3示出了按照现有技术的具有发送信道的成像,
图4示出了根据现有技术和根据本发明的脉冲序列的比较,
图5示出了根据现有技术和根据本发明的另一个脉冲序列,并且
图6示出了根据现有技术和根据本发明的具有多个发送信道的脉冲序列。
具体实施方式
图1示出了按照本发明的MR设备。该MR设备具有用于产生极化场B0的磁铁10。位于卧榻11上的检查对象12被驶入磁铁的同心,以便在那里通过入射HF脉冲和梯度场产生MR图像。如何能够通过HF脉冲的这样的脉冲序列来激励给定的磁化并且通过同时接通磁场梯度来建立MR图像,对专业人员来说是公知的并且在此不作详细解释。MR设备与中央的控制单元13相连,利用该控制单元能够控制MR检查的过程。控制单元13可以具有脉冲序列控制装置14,其控制何时必须产生HF脉冲以及何时必须产生梯度场的过程。同样,设置HF发送单元15,其优选具有多个HF发送信道,这些发送信道能够被单独控制并且被提供以HF功率。为了控制HF发送单元15,设置HF控制单元16,其对于不同的发送信道确定振幅和相位关系,以便在HF脉冲的情况下达到预定的磁化,例如目标磁化。此外设置用于在各个HF脉冲之后确定目标磁化的单元17。目标磁化例如可以经过输入单元18在对于系统中的操作人员的显示屏19上被输入,或者其可以由系统根据成像序列预先给出。计算单元20根据对于检查区域确定的目标磁化计算检查对象不同区域的目标翻转角和目标相位的图。然后图像计算器21可以由探测的磁化,如已知的,计算MR图像,然后该MR图像例如可以在显示单元19上被显示。在中央控制单元13中设置的单元作为单独的单元示出。当然这些单元不必构成为分离的单元。例如不同的单元可以组合在一个综合的单元中,或者其可以具有处理器,这些处理器借助命令执行属于单元的功能。
图2示意性示出了一个流程图,利用该流程图能够实现具有预定信号振幅的MR信号。在步骤100中对于第一HF脉冲开始该方法之后,在步骤110中定义对于脉冲序列的n个HF脉冲的目标磁化,其中,在最简单的情况下目标磁化都是相同的。此外,在下一步骤120中确定存在的磁化。在第一迭代步骤中,这例如可以是在B0场方向上的输出磁化,其中(M(x),M(y),M(z))=(0,0,M(z)0)。在下一步骤130中对于不同的空间区域,在使用各个存在的目标磁化和存在的磁化的条件下计算目标翻转角和目标相位的空间图。然后在步骤140中从目标翻转角和目标相位的空间图中计算对于每个发送信道的HF脉冲。这例如可以利用现有技术的方法进行,该方法进行所谓的RF补偿或者基于k空间的二维或三维脉冲设计。当HF发送单元具有多个发送信道时,对每个发送信道进行该计算。在步骤130和140中还可以考虑在图像序列设计中起作用的其它预定值。例如可以在考虑特殊吸收率(SAR)的条件下确定目标翻转角和目标相位,使SAR最小化。在步骤150中在每个HF脉冲之后计算磁化,即,计算脉冲对磁化的作用,并且在步骤160中计算脉冲之间磁化的演变并且作为存在的磁化反馈到步骤120,以便利用下一个磁化和在步骤160中计算的磁化又对下一个HF脉冲确定目标翻转角和目标相位的空间图。因为在入射下一个HF脉冲之前还存在现有的横向磁化,所以在计算目标翻转角和目标相位的空间图时应该考虑该现有的横向磁化。当对于脉冲序列的所有HF脉冲确定了振幅和相位关系时,该方法结束。在图2示出的实施方式的变形中,可以使用在步骤160中计算的磁化,以便在每个迭代中根据在步骤160中计算的实际上存在的磁化附加地更新目标磁化。
在另一个实施方式中,整个成像序列可以视为一个优化问题,这原则上意味着,可以使用对于脉冲的优化的现有求解方法,以便根据相同的标准优化整个脉冲序列。整个脉冲序列被作为一个HF脉冲处理,在此必须满足特定的边界条件:例如必须在数据采集期间的特定时间记录回波,或者在该时刻现有的磁化必须满足特定的边界条件。此外用于数据采集所必须的梯度是不可自由选择的。
在大多数成像序列中在不同的相位编码步骤中多次读出磁化,其中在每个相位编码步骤中实现一个或两个HF脉冲用于产生信号/回波。通过计算HF脉冲沿着脉冲序列、即典型地沿着相位编码的空间翻转角和相位分布,可以通过预先给定在检查对象的不同区域的预定信号振幅来减小由B1引起的信号非均匀性。
为了理解本发明,以下将根据现有技术和根据本发明的脉冲序列进行互相比较。
在图3中示出了按照现有技术的脉冲序列,在此利用发送信道产生HF脉冲。成像序列是HF脉冲31-34的序列,其中在每个HF脉冲之后进行信号读出35,直到信号记录以期望的分辨率结束。为了稳定信号振幅例如可以使用不同的翻转角和相位。例如以翻转角α1和相位
Figure G2009102542085D00081
来入射HF脉冲31,以翻转角α2和相位
Figure G2009102542085D00082
入射第二HF脉冲等等。
在图4的上部同样示出了按照现有技术的方法,在该方法中优化HF脉冲。为了横向磁化的空间调制,在发送信道中以时间上可变的梯度入射具有空间上变化的相位和空间上变化的翻转角的HF脉冲40,如通过振幅曲线46和梯度曲线47示意性示出的。在现有技术中该振幅、相位和梯度曲线对于所有的HF脉冲都是相同的。在根据本发明的方法中例如第一HF脉冲40具有空间可变的翻转角α1(x)和空间可变的相位
Figure G2009102542085D00083
第二HF脉冲41具有空间可变的翻转角α2(x)和相应的相位
Figure G2009102542085D00084
等等。同样示出了通过预先给定不同的翻转角计算出的HF脉冲的和梯度演变的时间上的曲线。第一HF脉冲40的振幅演变48与HF脉冲41的振幅演变49以及与HF脉冲42的振幅演变50不同。同样,HF脉冲40、41、42的梯度演变51、52、53是不同的。如从图4的下部可以看出的,对于不同的拍摄步骤,HF脉冲改变,因为对于每个步骤重新计算振幅和相位关系。在第一相位编码步骤中例如使用具有振幅48和预定相位的HF脉冲40,在第二相位编码步骤中使用具有振幅49和预定相位的HF脉冲,而在第三相位编码步骤中使用具有另一相位的振幅50。分别入射具有梯度曲线51-53的HF脉冲48-50。
图5示出,在多个并行运行的HF发送器中相位和振幅的校准如何能够导致空间调制的HF场。该方法作为RF补偿是公知的。图5的上部是根据现有技术公知的RF补偿,其中在每个HF脉冲54-56的情况下在每个信道中的相位和振幅保持相同。因此根据现有技术在第一信道中振幅为α1(x)并且相位为在此该振幅和相位在HF脉冲的曲线中不变。然而根据本发明的方法,对于每个脉冲重新计算振幅和相位权重,即,为了获得翻转角和相位,在四个信道中应用振幅a11-a41和相位而在第二HF脉冲55的情况下对于翻转角α2和相位
Figure G2009102542085D00087
使用振幅a12-a42和相位与图4类似的例子在具有多个发送信道的图6中示出。在图6的上部示出了根据现有技术的具有多个发送信道的空间选择的激励,其中为了达到翻转角和相位,HF脉冲57-59在不同的信道中具有不同的振幅曲线,如分别通过HF脉冲曲线60-63以及示例性的梯度曲线64示出的。这些示例性示出的振幅曲线和梯度曲线对于所有的HF脉冲57-59都是相同的。
然而按照本发明,利用第一HF脉冲57设置具有空间可变的相位的空间可变的翻转角α1(x),其在示出的例子中通过在不同的信道60-63中HF脉冲曲线和梯度曲线64达到,如上面的例子中那样。
然而,为了达到具有相位
Figure G2009102542085D00092
的翻转角α2,根据本发明,在第二HF脉冲时改变在不同的信道中的振幅关系和相位关系,如通过振幅曲线65-68和梯度曲线69示例性示出的。这使得可以考虑自旋系统经过整个脉冲序列的演变并且作为对于优化的其它自由度允许HF脉冲和梯度脉冲对于不同的记录步骤的变化。
参考文献
(1)Madhuranthakam AJ,Busse RF,Brittain JH,Rofsky NM,Alsop DC.B1-insensitive fast spin echo using adiabatic square wave enabling of the echo train(SWEET)excitation.Magn Reson Med 59(6)1386-1393,2008.
(2)Katscher U,
Figure G2009102542085D00101
P.Leussler C,van den Brink JS.Transmit SENSE.Magn Reson MEd 49(1)144-150,2003.
(3)Zhu Y.Parallel excitation with an array of transmit coils.Magn Reson Med51(4)775-784,2004.
(4)Saekho S,Yip C,Noll DC,Boada F E,Stenger,VA.Fast-kzthree-dimensional tailored radiofrequency pulse for reduced B1 inhomogeneity.Magn Reson Med 55(4),719-724,2006.
(5)Setsompop K,Wald LL,Alagappan V,Gagoski B,Hebrank F,Fontius U,Schmitt F,Adalsteinsson E.Parallel RF transmission with eight channels at 3 Tesla.Magn Reson Med 56(5),1163-1171,2006
(6)Scheffler K.A pictorial description of steady-states in rapid magneticresonance imaging.Magn Reson Imag 1999;11:291-304.
(7)Hennig J.Echoes-how to generate,recognize,use or avoid them inMR-imaging sequences.Conc Magn Reson 1991;3:125-143.
(8)Alsop DC.The sensitivity of low flip angle RARE imaging.Magn ResonMed 1997;37:176-184.
(9)Mugler JP,Kiefer B,Brookeman JR.Three-dimensional T2-weightedimaging of the brain using very long spin-echo trains.In:Proceedings of the 8thAnnual Meeting of ISMRM,Denver,CO,USA,2000.p.687.
(10)Busse RF,Riederer SJ.Reducing SAR in real-time SSFSE imaging withvariable-flip hard refocusing RF pulses.In:proceedings of the 9th Annual Meeting ofISMRM,Glasgow,Scotland,2001.p.1790
(11)Mugler  JP.Practical Implementation of optimized tissue-specificprescribed signal evolutions for improved turbo-spin-echo imaging.In:Proceedingsof the 11th Annual Meeting of ISMRM,Toronto,Canada,2003.p.203.
(12)Hennig J,Weigel M,Scheffler K.Multiecho sequences with variablerefocusing flip angles;optimization of signal behavior using smooth transitionsbetween pseudo steady states(TRAPS).Magn Reson Med 2003;49:527-535.
(13)Hennig J.Calculation of flip angles for echo trains with predefinedamplitudes with the extended phase graph(EPG)-algorithm:principles andapplications to hyperecho and TRAPS sequences.Magn Reson Med 2004;51:68-80.
(14)LeRoux P,Hinks RS.Stabilization of echo amplitudes in FSE sequences.Magn Reson Med 1993;30:183-191.

Claims (10)

1.一种用于在磁共振测量中确定检查对象(12)的预定信号振幅的方法,其中,在脉冲序列中将具有小于检查对象的T2时间的重复时间TR的多个高频脉冲以脉冲序列的形式入射到检查对象中,该方法具有以下步骤:
对在入射各高频脉冲之后的预定时刻确定对于该脉冲序列的所有高频脉冲的目标磁化;
根据相应的目标磁化和现有的磁化对检查对象的不同区域对所有高频脉冲确定目标翻转角和目标相位;
对于所有高频脉冲分别确定振幅和相位的时间变化,以在入射各个高频脉冲之后产生各个目标磁化;以及
入射具有相应确定的振幅和相位关系的高频脉冲。
2.根据权利要求1所述的方法,其特征在于,利用多个高频发送信道同时入射所述高频脉冲,其中,对于所有高频发送信道确定振幅和相位关系。
3.根据权利要求1所述的方法,其特征在于,在入射脉冲序列的第i个高频脉冲之后计算磁化,并在第i+1个高频脉冲时在确定对于检查对象不同区域的目标翻转角和目标相位时将该磁化用作开始磁化。
4.根据权利要求3所述的方法,其特征在于,对于开始磁化的确定考虑第i个和第i+1个高频脉冲之间磁化的演变。
5.根据权利要求1至4中任一项所述的方法,其特征在于,所述目标磁化对于脉冲序列的所有高频脉冲基本相同。
6.根据权利要求1至4中任一项所述的方法,其特征在于,在确定目标翻转角和目标相位时考虑由脉冲序列引起的边界条件。
7.根据权利要求1至4中任一项所述的方法,其特征在于,对于脉冲序列的每个高频脉冲同时接通磁场梯度,该磁场梯度是根据每个高频脉冲的相应目标翻转角和相应目标相位针对不同的区域选择的。
8.根据权利要求1至4中任一项所述的方法,其特征在于,在考虑脉冲序列的最小特殊吸收率的条件下确定高频脉冲的目标翻转角和目标相位。
9.一种用于在磁共振测量中确定检查对象的预定信号振幅的磁共振设备,其中在脉冲序列中将多个高频脉冲以脉冲序列的形式入射到检查对象中,该磁共振设备具有:
高频发送单元(15);
单元(17),用于在入射脉冲序列的各个高频脉冲之后的预定时刻对脉冲序列的所有高频脉冲确定目标磁化;
计算单元(20),用于根据各个目标磁化对于检查对象的不同区域对所有高频脉冲自动确定目标翻转角和目标相位;
高频控制单元(16),对所有高频脉冲分别确定振幅和相位关系,以在入射各个高频脉冲之后产生相应的目标磁化,其中,高频发送单元将具有相应确定的振幅和相位关系的高频脉冲入射到检查对象中。
10.根据权利要求9所述的磁共振设备,其特征在于,所述高频发送单元具有多个发送信道,
其中,所述高频控制单元对所有高频发送信道和所有高频脉冲确定振幅和相位关系,以实现各个高频脉冲的目标翻转角和目标相位。
CN200910254208.5A 2008-12-10 2009-12-10 磁共振测量中用于确定预定信号振幅的方法 Expired - Fee Related CN101750596B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102008061455A DE102008061455B4 (de) 2008-12-10 2008-12-10 Verfahren und Vorrichtung zur Bestimmung einer vorbestimmten Signalamplitude bei MR-Messungen
DE102008061455.6 2008-12-10

Publications (2)

Publication Number Publication Date
CN101750596A CN101750596A (zh) 2010-06-23
CN101750596B true CN101750596B (zh) 2014-03-12

Family

ID=42193984

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910254208.5A Expired - Fee Related CN101750596B (zh) 2008-12-10 2009-12-10 磁共振测量中用于确定预定信号振幅的方法

Country Status (3)

Country Link
US (1) US8427157B2 (zh)
CN (1) CN101750596B (zh)
DE (1) DE102008061455B4 (zh)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0619269D0 (en) * 2006-09-29 2006-11-08 King S College London Method of mapping a magnetic field for use in magnetic resonance imaging
DE102008061455B4 (de) * 2008-12-10 2011-03-17 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Bestimmung einer vorbestimmten Signalamplitude bei MR-Messungen
DE102010027312B4 (de) * 2010-07-16 2017-10-12 Siemens Healthcare Gmbh Verfahren zur Erstellung von MR-Bildern sowie entsprechend ausgestaltete Magnetresonanzanlage
DE102010063565B4 (de) * 2010-12-20 2014-09-04 Siemens Aktiengesellschaft Verfahren und Einrichtung zur Ermittlung einer Magnetresonanzsystem-Ansteuersequenz
CN102654568A (zh) 2011-03-01 2012-09-05 西门子公司 用来确定对于磁共振成像的激励参数的方法和装置
DE102011078849B4 (de) * 2011-07-08 2013-02-21 Siemens Aktiengesellschaft Festlegen von Anregungsparametern für eine parallele MR-Bildgebung unter Verwendung von Elementen einer k-Raum-Kovarianzmatrix für Signalrauschen in k-Raum-Daten
DE102011005174A1 (de) * 2011-03-07 2012-09-13 Friedrich-Alexander-Universität Erlangen-Nürnberg Verfahren und Einrichtung zur Ermittlung einer Magnetresonanzsystem-Ansteuersequenz
DE102011006149B4 (de) * 2011-03-25 2012-10-31 Siemens Aktiengesellschaft Ermittlung einer Magnetresonanzsystem-Ansteuersequenz unter Optimierung der Sendebandbreite
DE102011006151B4 (de) 2011-03-25 2014-04-30 Siemens Aktiengesellschaft Ermittlung einer hinsichtlich eines HF-Energieparameters optimierten HF-Pulslänge eines Magnetresonanzsystems
DE102011007825B4 (de) 2011-04-20 2012-11-22 Bruker Biospin Mri Gmbh Verfahren zur Bestimmung der räumlichen Verteilung von Magnetresonanzsignalen in Subvolumen eines Untersuchungsobjektes
DE102011007823B4 (de) * 2011-04-20 2012-11-15 Bruker Biospin Mri Gmbh Verfahren zur Erzeugung eines gewünschten zeitlichen Verlaufs des Magnetisierungszustandes in einem Untersuchungsobjekt im Rahmen eines Experimentes der magnetischen Resonanz
DE102011081509B4 (de) * 2011-08-24 2013-03-28 Siemens Aktiengesellschaft Ermittlung einer Ansteuerungssequenz für eine Magnetresonanzeinrichtung aus einer Mehrzahl von Rechtecksubpulsen
DE102013205785B4 (de) * 2013-04-02 2014-10-09 Siemens Aktiengesellschaft Ermittlung einer Magnetresonanzsystem-Ansteuersequenz auf Basis einer reduzierten Anzahl an Feldverteilungskarten
JP6621978B2 (ja) * 2013-12-16 2019-12-18 キヤノンメディカルシステムズ株式会社 Mri装置
US10359486B2 (en) * 2016-04-03 2019-07-23 Q Bio, Inc. Rapid determination of a relaxation time
US10222441B2 (en) * 2016-04-03 2019-03-05 Q Bio, Inc. Tensor field mapping
JP7152956B2 (ja) * 2016-05-31 2022-10-13 キュー バイオ インコーポレーション テンソル場マッピング
US11354586B2 (en) 2019-02-15 2022-06-07 Q Bio, Inc. Model parameter determination using a predictive model
CN111157932B (zh) * 2020-01-02 2022-08-30 华东师范大学 一种快速自旋回波脉冲序列中射频脉冲的优化方法
US11614508B1 (en) 2021-10-25 2023-03-28 Q Bio, Inc. Sparse representation of measurements

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2831147A (en) * 1948-04-06 1958-04-15 Weber Joseph Electronic frequency analyzer device
US4684891A (en) * 1985-07-31 1987-08-04 The Regents Of The University Of California Rapid magnetic resonance imaging using multiple phase encoded spin echoes in each of plural measurement cycles
US5270653A (en) * 1992-04-10 1993-12-14 The Board Of Trustees Of The Leland Stanford Junior University Selective three dimensional excitation for magnetic resonance applications
DE19511794B4 (de) * 1995-03-30 2005-12-29 Siemens Ag Verfahren zur Gewinnung von Bilddaten in einem Kernspintomographiegerät und Kernspintomographiegerät zur Durchführung des Verfahrens
DE19543079A1 (de) * 1995-11-18 1997-05-22 Philips Patentverwaltung Verfahren zum Bestimmen der räumlichen und/oder spektralen Verteilung der Kernmagnetisierung
JP3753668B2 (ja) * 2002-03-12 2006-03-08 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Rfパルスチューニング装置
US7053618B2 (en) * 2003-11-26 2006-05-30 General Electric Company Method and apparatus to generate an RF excitation consistent with a desired excitation profile using a transmit coil array
US6989673B2 (en) * 2003-11-26 2006-01-24 General Electric Company Method and apparatus to reduce RF power deposition during MR data acquisition
WO2007098011A2 (en) * 2006-02-17 2007-08-30 Regents Of The University Of Minnesota High field magnetic resonance
US7385396B2 (en) * 2006-04-20 2008-06-10 General Electric Company SAR reduction in MR imaging with parallel RF transmission
DE102006058162B3 (de) * 2006-12-09 2008-06-12 Bruker Biospin Mri Gmbh Verfahren zur Gewinnung von Amplituden- und Phasenverläufen von HF-Pulsen für die räumlich-selektive Anregung
DE102007004620B4 (de) * 2007-01-30 2010-02-04 Siemens Ag Verbessertes dreidimensionales schichtselektives Mehrschicht-Anregungsverfahren in der MRT-Bildgebung
DE102007044463B4 (de) * 2007-09-18 2009-05-14 Bruker Biospin Mri Gmbh Verfahren zur Bestimmung der räumlichen Verteilung von Magnetresonanzsignalen durch mehrdimensionale HF-Anregungspulse
DE102008029175B4 (de) * 2008-06-19 2010-09-30 Siemens Aktiengesellschaft Verfahren zur Ermittlung einer Pulssequenz zur Ansteuerung einer Hochfrequenz-Sendespule
DE102008061455B4 (de) * 2008-12-10 2011-03-17 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Bestimmung einer vorbestimmten Signalamplitude bei MR-Messungen
DE102009002112B3 (de) * 2009-04-01 2010-08-19 Bruker Biospin Mri Gmbh Verfahren zur Gewinnung von Amplituden- und Phasenverläufen von HF-Pulsen für die räumlich-selektive Anregung
DE102009016341B4 (de) * 2009-04-06 2011-02-10 Siemens Aktiengesellschaft Verfahren zur Bestimmung von k-Raumpositionen für eine Modellierung von HF-Pulsen für Magnetresonanzanregungen sowie Magnetresonanzgerät und Computerprogramm zur Durchführung des Verfahrens

Also Published As

Publication number Publication date
US20100141252A1 (en) 2010-06-10
CN101750596A (zh) 2010-06-23
US8427157B2 (en) 2013-04-23
DE102008061455A1 (de) 2010-06-24
DE102008061455B4 (de) 2011-03-17

Similar Documents

Publication Publication Date Title
CN101750596B (zh) 磁共振测量中用于确定预定信号振幅的方法
US8698495B2 (en) Flip angle imaging with improved B1 mapping for multi-RF transmit systems
Katscher et al. Transmit sense
Vannesjo et al. Gradient system characterization by impulse response measurements with a dynamic field camera
US8228063B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
US6369568B1 (en) Fast spin echo phase correction for MRI system
CN102293650B (zh) 磁共振成像装置以及方法
JP6018401B2 (ja) 拡散強調エコープラナー撮像法において高次渦電流に誘発された歪みを予測補正するためのシステムおよび方法
US9041394B2 (en) Magnetic resonance imaging apparatus executing phase-corrected imaging pulse sequence based on data obtained from multiple pulse sub-sequences executed without readout direction gradients but instead using phase or slice encoding direction gradients during readout
US10401456B2 (en) Parallel MR imaging with Nyquist ghost correction for EPI
US10274566B2 (en) Dixon-type water/fat separation MRI using high-SNR in-phase image and lower-SNR at least partially out-of-phase image
US9494668B2 (en) Systems and methods for measuring nuclear magnetic resonance spin-lattice relaxation time T1 and spin-spin relaxation time T2
CN103278785B (zh) 一种快速自旋回波脉冲序列中射频脉冲相位的优化方法
Liu et al. DWI using navigated interleaved multishot EPI with realigned GRAPPA reconstruction
Malik et al. Spatially resolved extended phase graphs: modeling and design of multipulse sequences with parallel transmission
US10101428B2 (en) Method and magnetic resonance apparatus for generating at least one combination image dataset
RU2015116879A (ru) Опорное сканирование при устойчивой к металлам мр визуализации
US7239140B1 (en) Method and apparatus for vibration-related artifact reduction
CN105473067A (zh) 用于降低mr成像中声音噪声水平的系统和方法
Malik et al. Direct signal control of the steady‐state response of 3D‐FSE sequences
CN104515960B (zh) 磁共振设备中三维梯度回波序列的噪声发展的优化
US7239137B2 (en) Method and apparatus for fast spin echo (FSE) prescan phase correction
Corbin et al. Robust 3D Bloch‐Siegert based mapping using multi‐echo general linear modeling
KR20180007327A (ko) 자기 공명 기술에서의 아티팩트 감소
Nussbaum et al. Thermal variation in gradient response: measurement and modeling

Legal Events

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

Granted publication date: 20140312

Termination date: 20201210