发明内容
本发明要解决的技术问题是,提出一种用于确定开头提到的种类的磁共振控制序列的方法和相应的控制序列确定装置,利用其可以避免或至少减少在旋转对称的激励轮廓的二维选择性激励时通过梯度延迟时间引起的上述问题。
上述技术问题通过按照本发明的方法以及通过按照本发明的控制序列确定装置解决。
在按照本发明的方法中,如上所述,确定具有脉冲排列的磁共振控制序列,所述脉冲排列在至少两个空间的、优选互相正交的方向上选择性地起作用,以便激励在检查对象内部的受限的旋转对称的激励轮廓(也称为“激励体积”)。该脉冲排列具有二维选择性HF激励脉冲,其包括多个分开地、即时间上互相间隔开的HF子脉冲的序列。并行于HF子脉冲的序列,脉冲排列为了实现期望的选择作用而包括梯度脉冲,其在两个空间方向上协调地布置为,使得不同的HF子脉冲的HF能量输入(即HF能量的沉积)从发送k空间来看,分别在互相同心的、圆形的k空间发送轨迹上进行。即,各个HF子脉冲分别对应于一个自己的圆形发送轨迹(以下也称为“环形轨迹”),并且这些环形轨迹互相同心地位于一个平面中,该平面通过属于梯度的空间方向的k空间方向张开。在此(幅度调制的)HF子脉冲的包络线的幅度,以下也称为b1幅度,分别在所属的环形轨迹的一个遍历的持续时间期间(即在一个HF子脉冲的HF入射的持续时间T期间)是恒定的。
为了实现该轨迹形状,可以在静态的、至少近似均匀的B0场上分别在时间间隔T,即在发送轨迹的各自的圆形轨线上的环绕时间或者说HF子脉冲的发送时间期间,分别叠加第一梯度子脉冲(或者说第一梯度场)和同时第二梯度子脉冲(或者说第二梯度场),所述第一梯度子脉冲的幅度在时间间隔T期间正弦形地以周期T振荡并且其方向垂直于激励体积的对称轴,并且所述第二梯度子脉冲的幅度在时间间隔T期间余弦形地以周期T振荡并且其方向垂直于对称轴并且垂直于第一梯度场的方向。
基本上,这样的HF子脉冲类似于在US4812760中描述的空间选择性反转或重聚焦脉冲(“π脉冲”)来构造,其也具有矩形的包络线和沿着第一梯度方向的正弦形的梯度曲线以及沿着第二梯度方向的余弦形梯度曲线。然而,利用这样的“π脉冲”不能实现足够的空间定位。只有通过按照本发明将空间选择性HF脉冲由在互相同心的环形轨迹上的多个、即至少两个,但是优选更多个HF子脉冲组成才可以选择性地激励期望的圆柱形对称的激励轮廓。
通过如下,即,按照本发明对于选择性激励现在不是选择EPI或螺旋轨迹而是选择多个同心的环形轨迹,并且在此如后面还要示出的那样B1幅度在分别一个环形轨迹的遍历期间可以恒定,类似于一维选择性HF脉冲的情况(在随时间改变的B1幅度的情况下具有恒定的梯度幅度),尽管存在梯度延迟时间,还是可以实现选择轮廓的不寻常的鲁棒性。
在此已经指出,脉冲排列可以具有多个其他HF脉冲和梯度脉冲,其中也可以是选择性激励的脉冲。也就是不需要,所有子脉冲位于同心的环形轨迹上。例如,空间选择性作用的脉冲序列也可以附加地具有位于螺旋轨迹上的HF脉冲。特别地可以的是,例如k空间的外部区域在使用(快速)螺旋轨迹的情况下并且k空间的内部区域在使用(精确的)按照本发明同心的环形轨迹的情况下被激励。
按照本发明的用于确定这样的磁共振控制序列的控制序列确定装置需要至少一个输入接口装置,用于采集激励轮廓数据,其定义检查对象内部的受限的待激励的旋转对称的激励轮廓的空间伸展。这些激励轮廓数据例如可以是旋转对称的激励轮廓的位置数据,例如旋转轴的位置或方位,以及待激励的体积的直径。输入接口装置例如可以包括用户接口,以便例如手动地输入这些位置数据,特别是图形的用户接口,以便将期望的激励轮廓在检查对象的示出的图像数据中,例如定位片中标出或以其他方式标记。替换地或附加地,输入接口装置也可以具有用于(例如对于导航应用)自动产生对于检查对象的合适的激励轮廓数据或从相应的装置或存储器调用其的装置。
此外控制序列确定装置按照本发明还需要脉冲排列确定单元,用于确定至少一个脉冲排列,其在至少两个空间方向上选择性地作用,以便精确定义并限制地激励借助输入接口装置预先给出的激励轮廓。该脉冲排列确定单元这样构造,即,其这样确定脉冲排列的HF激励脉冲,使得其包括多个HF子脉冲的序列并且并行于HF子脉冲的序列这样协调地将梯度脉冲施加到两个空间方向上,使得在发送k空间中不同的HF子脉冲的HF能量输入分别在相互同心的圆形k空间发送轨迹(环形轨迹)上进行并且在此HF子脉冲的包络线的幅度分别在环形轨迹的一个遍历的持续时间期间是恒定的。
此外控制序列确定装置应当具有合适的控制序列输出接口,用于将序列传输到磁共振断层造影系统的其他控制单元。控制序列输出接口例如可以是如下的接口,其将序列传输到磁共振控制器,以便由此直接控制测量,但是也可以是如下的接口,其将数据经过网络发送和/或在存储器中为以后的使用而存储。
按照本发明的磁共振系统除了用于发送HF脉冲的高频发送装置之外还具有用于接通需要的梯度的梯度系统和控制装置,该控制装置构造为,用于为了基于预先给出的磁共振控制序列执行期望的测量而发送高频脉冲串(即HF脉冲)并且与之协调地经过梯度系统发送梯度脉冲串(即所属的梯度脉冲)。此外磁共振系统还具有上面描述的控制序列确定装置,以便以按照本发明的方式确定控制序列并且将其传输到控制装置。
相应地在按照本发明的用于运行磁共振系统的方法中按照前面描述的方法确定控制序列并且然后在使用控制序列的条件下运行磁共振系统。
如后面还要根据测试测量解释的,还可以在用于关于梯度延迟时间来调节和/或校准磁共振系统的方法的内部有利地采用按照本发明的控制方法。为此例如可以利用按照本发明的方法多次在人工改变梯度延迟时间的情况下,即在具有不同设置的梯度延迟时间的每个单个测量的情况下来控制、即运行磁共振系统。在此分别产生磁共振图像数据,例如确定的激励轮廓的简单的剖面图,其应当位于检查对象中的和/或在为了调节和/或校准测量而设置的模体等中的完全确定的位置。在这些重复的测量中优选反复地选择性地激励相同的激励轮廓。基于该磁共振图像数据的分析,然后可以关于梯度延迟时间进行磁共振系统的期望的调节或校准。该分析特别地可以例如视觉地通过用户进行,方式是,其在不同的拍摄之间反复地改变梯度延迟时间并且然后视觉地在图像数据中检查,激励轮廓是否位于期望的位置上。该调节和/或校准测量在改变梯度延迟时间的条件下一直重复,直到测量的激励轮廓精确地位于期望的位置上并且由此系统内在的(首先未知的)延迟时间被精确地补偿。由此确定的延迟时间也可以被记录并且例如在其他测量时被考虑。
控制序列确定装置的主要部分可以按照软件组件的形式构造。这特别地涉及脉冲排列确定单元。同样,提到的接口可以至少部分地以软件的形式构造并且可能援用现有的计算机的硬件接口。本发明由此也包括计算机程序或计算机程序产品,其可以直接加载到控制序列确定装置的存储器中,具有程序代码段,用于当程序在控制序列确定装置中运行时执行按照本发明的方法的所有步骤。这样的按照软件的实现具有优点,即,可以通过实现程序而以合适的方式修改为了确定控制序列而使用的迄今为止的装置,以便以按照本发明的方式确定最佳的控制序列。
从属权利要求以及以下的描述包含本发明的特别有利的扩展和构造,其中特别地一类权利要求也可以类似于另一类权利要求的从属权利要求来扩展并且不同的实施例的特征也可以组合以形成另外的实施例。
对于HF子脉冲在单个环形轨迹上的精确构造存在不同的可能性。可能的参数在此是为了遍历一个环形轨迹而所需的持续时间。
在第一变形中确保,圆形的环形轨迹的遍历的持续时间对于沿着圆形的环形轨迹作用的、HF激励脉冲的至少两个不同的HF子脉冲,在极端情况下对于所有HF子脉冲,是相同长的。
在此是最简单的可能性。工作方式类似于在二维螺旋脉冲情况下(即在沿着k空间中螺旋形轨迹发送HF脉冲的情况下)的恒定的角速度。
在一种替换的变形中,HF子脉冲和所属的梯度脉冲构造为,圆形的环形轨迹的遍历的持续时间对于至少两个不同的HF子脉冲分别单独地确定,即,不同的HF子脉冲可以特别是不同长的。通过传播时间的个性化,在以下意义上实现序列的优化,即,可以最小化HS激励脉冲的总持续时间并且尽管如此还可以可靠地遵循梯度脉冲的边界调节,诸如最大允许的转换速率(slew-rate)(梯度增长速率)和/或最大允许的梯度幅度。
为了实现这一点,优选地可以将每个HF子脉冲的圆形的k空间轨迹的遍历的持续时间选择为这样短,使得梯度脉冲(刚好)不超过预先给出的最大转换速率。这一点例如类似于具有恒定的梯度增长速率的螺旋形激励。
不管不同的HF子脉冲的脉冲持续时间是相同长的还是单独设置的,优选地单独地确定HF子脉冲的包络线的分别在圆形的发送轨迹的遍历的持续时间期间保持恒定的幅度(即b1幅度),即,HF子脉冲可以特别地具有不同的b1幅度。
为了实现特别快速的HF激励脉冲,即,所有的同心环形轨迹的遍历的总时间保持为尽可能小,可以将梯度脉冲优选地构造为,将与时间上相邻的HF子脉冲对应的两个同心的圆形k空间发送轨迹相反地在发送k空间中遍历。即,保证,例如将第一环形轨迹在顺时针方向遍历,将直接相邻的环形轨迹则逆时钟方向并且再下一个环形轨迹又顺时针遍历。如后面还要示出的,在该变形中通过在单个子脉冲之前和之后将预相位梯度和重聚相位梯度综合可以节省特别多的时间。
在HF激励脉冲的精确构造中与同心的环形轨迹有关的另一个参数是环形轨迹环形的间隔。
在第一变形中,这样构造梯度脉冲,使得圆形的k空间发送轨迹形成在发送k空间中的围绕共同的中心的等距的环。即在该变形中在两个环形轨迹之间的径向距离总是相同的。就此实现了能量特别均匀地输入到k空间中。
在替换的变形中,这样构造梯度脉冲,使得同心的圆形k空间发送轨迹形成这样的环,所述环的相互的径向距离是不同的。在特别优选的变形中同心的圆形k空间发送轨迹形成这样的环,所述环相互的径向距离朝着共同的中心变小。通过与k空间中心的更密的环形轨迹距离可以保证,特别地在具有关于激励轮廓的主信息的重要的频率范围中被沉积足够的能量。为此所需的时间于是可以在k空间的不太重要的边缘区域中被节省,方式是在那里在环形轨迹相互之间的距离更小。
另一个可设置的参数是HF激励脉冲的环形轨迹的数量。一方面有利的是,将环形轨迹相对密地设置,以便足够密地实现空间中的激励。另一方面随着环形轨迹数量增加,用于HF激励脉冲的总时间也必定延长。
为此特别有利的是,根据第一边激励的距离与激励轮廓的直径的比例来选择同心的、圆形的k空间发送轨迹的数量。因为在第一边激励的区域中不再期望任何激励并且由此边激励与激励轮廓的对称轴的距离和(圆柱形对称的)激励轮廓的直径是用于激励的关键参数,所以强调,足够的是,相应于第一边激励的该距离的双倍的值除以激励轮廓的直径来选择环形轨迹的数量。由此确保,环形轨迹相互的距离至少相应于在按照尼奎斯特理论的扫描情况下的最小距离。
另一个参数是在激励时在激励轮廓内部待实现的目标翻转角(目标激励角)或目标翻转角分布。
在第一变形中这样确定HF子脉冲的包络线的幅度,使得目标翻转角在圆柱形的激励轮廓的内部取决于与圆柱形激励轮廓的圆柱体轴的径向距离r。即,目标翻转角从内向外减小。特别优选地在此与径向距离的依赖关系是高斯形状的。
在替换的方法中这样确定HF子脉冲的包络线的幅度,使得目标翻转角在圆柱形激励轮廓内部基本上,即,在可能的容差内,是恒定的。在该变形中目标翻转角分布也就是关于激励轮廓是尽可能均匀的。
此外其他分布也是可能的,例如在内部区域中几乎恒定的或均匀的激励而在该均匀区域外部则特别是高斯形状的下降。
在此还要再次指出,上面提到的不同的变形对于各个参数,诸如在环形轨迹内部的传播时间、环形轨迹的距离、B1幅度或翻转角分布,可以任意组合。
具体实施方式
图1粗略示意性示出按照本发明的磁共振设备1。其一方面包括具有位于其中的检查空间8或患者通道8的本来的磁共振扫描仪2。卧榻7可以驶入该患者通道8中,从而卧于其上的患者O或受检者在检查期间可以在磁共振扫描仪2内部的特定位置上相对于布置于其中的磁体系统和高频系统放置或在测量期间可以在不同的位置之间移动。
磁共振扫描仪2的组件包括基本场磁体3、具有磁场梯度线圈以便产生在x、y和z方向上的磁场梯度的梯度系统4,以及全身高频线圈5。在x、y和z方向(空间坐标系)上的磁场梯度线圈互相独立地可控,从而通过预先给出的组合可以施加在任意逻辑空间方向上、例如在层选择方向上、在相位编码方向上或在读出方向上施加梯度,其不一定平行于空间坐标系的轴。在检查对象O中感应的磁共振信号的接收可以经由全身线圈5进行,利用所述全身线圈通常也发送高频信号以感应磁共振信号。但是通常这些信号利用例如带有置于检查对象O上或下的局部线圈(在此仅示出其中的一个)的局部线圈装置6来接收。所有这些组件对于专业人员是基本公知的并且由此在图1中仅粗略示意示出。
磁共振扫描仪2的组件可以由控制装置10控制。在此可以是控制计算机,其也可以由多个(可能空间上分离的和经过合适的电缆等互相连接的)单计算机组成。通过终端接口17,该控制装置10与终端20相连,通过所述终端,操作者可以控制整个设备1。在本情况中该终端20构造为具有键盘、一个或多个显示器以及其他输入设备例如鼠标等的计算机,从而可以为操作者提供图形的用户界面。
控制装置10尤其还具有梯度控制单元11,其又可以由多个部分组件组成。通过该梯度控制单元11利用控制信号按照梯度脉冲序列GS接通各个梯度线圈。在此如上所述是在测量期间在精确地预先给出的时间位置上并且以精确地预先给出的时间曲线设置的梯度线圈。
控制装置10还具有高频发送单元12,用于按照控制序列AS的预先给出的高频脉冲序列HFS将高频脉冲分别馈入到全身高频线圈5中。高频脉冲序列HFS包括上面提到的选择性激励脉冲。磁共振信号的接收于是借助局部线圈装置6进行,并且由此接收的原始数据RD由HF接收单元13读出并且处理。磁共振信号以数字的形式作为原始数据RD被传输到重建单元14,其由此重建图像数据BD并且将其存储在存储器16中和/或通过接口17传输到终端20,从而操作者可以对此进行观察。图像数据BD也可以通过网络NW存储到另外的位置和/或显示和分析。替换地,也可以经过局部线圈装置发送高频脉冲序列和/或磁共振信号可以由全身高频线圈接收(未示出)。
经过另外的接口18将控制命令传输到磁共振扫描仪2的其他组件,诸如卧榻7或基本场磁体3,或者接收测量值或其他信息。
梯度控制单元11、HF发送单元12和HF接收单元分别协调地通过测量控制单元15来控制。该测量控制单元通过相应的命令确保,发送期望的梯度脉冲序列GS和脉冲序列的高频脉冲序列HFS。此外必须确保,在合适的时刻将在局部线圈装置6的局部线圈上的磁共振信号通过HF接收单元13读出并进一步处理,即,必须例如通过将HF接收单元13的ADC接通到接收来设置读出窗。同样,测量控制单元15控制接口18。
这样的磁共振测量的基本流程和用于控制的提到的组件对于专业人员来说是公知的,从而在此不再详细讨论。此外这样的磁共振扫描仪2以及所属的控制装置10还具有多个另外的组件,所述组件在此同样不再详细解释。在此要指出,磁共振扫描仪2也可以另外地构造,例如具有侧面开口的患者空间,或者构造为更小的扫描仪,其中仅能放置一个身体部位。
为了开始测量,用户可以经过终端20从其中存储了用于不同的测量的多个控制协议P的存储器1中通常选择用于该测量设置的控制协议P。该控制协议P尤其包含用于各自的测量的不同的控制参数值SP。属于这些控制参数值SP的例如有序列类型、对于各个高频脉冲的目标磁化、回波时间、重复时间、不同的选择方向等。同样在此已经可以预先给出层厚、分辨率、层的数量或在3D激励的情况下预先给出厚片厚度或任意的激励轮廓的其他尺寸(即激励轮廓数据)。此外控制参数值SP还可以包括信息:相应的测量是否是具有导航的测量,如果是,可能已经是用于对于该导航应用所需的激励轮廓EP的激励轮廓数据。在图1示意性示出的应用中该激励轮廓EP可以是在z方向上在患者O的身体中延伸的圆柱形对称的所谓“铅笔束”,其通过患者O的横膈膜延伸。
所有这些控制参数值SP尤其经过控制序列确定装置22的输入接口24提供,由此其确定合适的控制序列AS。同样控制参数值SP也可以在调用该协议的情况下首先提供给用户用于接收,并且其可以将该值也借助用户界面任意改变并且匹配于当前的检查任务。特别地,用户在此也可以经过具有所属的图形用户界面的终端20的计算机21确定激励轮廓数据,例如直径d、旋转对称的激励轮廓EP的对称轴等的精确位置或者已经通过控制参数值SP改变在协议P中定义的激励轮廓数据。这一点在图1中通过接口装置23表示。
此外用户可以替代从存储器16也可以经过网络NW例如从磁共振设备的制造者处调用具有相应的控制参数值SP的控制协议,并且然后如下面描述的那样使用它们。
基于包括激励轮廓数据的控制参数值SP然后确定控制序列AS,按照所述控制序列最后通过测量控制单元15进行其余组件的控制。该控制序列于是尤其包含脉冲排列,以便选择性地激励选择的激励轮廓EP。控制序列AS如提到的那样在作为终端20的部分示出的控制序列确定装置22中被计算,并且经过控制序列输出接口25传输到磁共振扫描仪2的控制装置10。控制序列确定装置22尤其包含(其他组件在图1中为清楚起见没有示出)脉冲排列确定单元26,其确定用于激励轮廓EP的选择性激励的所述脉冲排列。整个控制序列确定装置22和其组件例如可以按照软件的形式在一个或多个合适的处理器上实现。控制序列确定装置22和其各个组件的精确工作方式在以下进一步解释,其中仅以上面提到的铅笔束激励为例,但是不将本发明限制于此。
为例设计二维HF脉冲,以激励圆柱形对称的激励轮廓EP的位置空间中的、与对称轴z垂直的x/y平面,必须除了该平面中的激励轮廓之外一般地首先选择k空间中所属的x/y平面内部的二维k空间发送轨迹(具有k空间中取决于时间的坐标(kx(t),ky(t));以下也简称为(kx,ky))。激励轮廓在此通过具有坐标值x和y的位置处的函数P(x,y)确定。函数P(x,y)说明了相对于均衡磁化的绝对值的、横向磁化的期望的绝对值,也就是没有单位。该轨迹的k空间覆盖(伸展)确定了空间分辨率,利用所述分辨率可以实现期望的选择轮廓。K空间中轨迹的扫描密度确定了位置空间中激励的第一边带的距离,以及FOV(Field of View=视野)。这样的边带由于离散的扫描而是不可避免的。
在高频入射期间二维梯度场(gx(t),gy(t))通过如下与所选的k空间轨迹相关联:
及
在此公式(1)中的负号是在RF脉冲的末端(t=T)开始积分的传统的结果。时间t说明直到RF脉冲末端的剩余时间。
HF脉冲(b1(t)在此是取决于时间的幅度)的所属的脉冲形状b1(t)通过期望的空间激励轮廓P(x,y)关于该期望的激励轮廓P(x,y)的、加权的二维傅里叶变换得到:
该公式从布洛赫等式(在小角度激励的限制中)导出(例如参见J.Pauly等的“A k-Space Analysis of Small-Tip-Angle Excitation”in Journal of Magn.Res.81,43至56页,1989和C.J.Hardy等的“Correcting for Nonuniform k-SpaceSampling in Two-Dimensional NMR Selective Excitation”in Journal of Magn.Res.87,639至645页,1990)。
公式(3)中第一加权系数
是k空间中的平移速度。
第二系数将k空间中的位置(kx,ky)处的k空间扫描的密度加权:
在公式(5)中ΔA(kx,ky)是k空间点(kx,ky)的周围中k空间的小的面积并且Δl(kx,ky)是由面积包围的轨迹的长度。按照公式(4)和(5)的系数直观地意味着,在k空间被缓慢地或以高的密度扫描的地方降低HF幅度,反之亦然。
根据按照本发明的方法为了激励轮廓的选择性激励应当选择具有至少两个同心的环形轨迹的k空间轨迹。在图2中示例性示出这样的k空间轨迹,其具有总共八个同心的环形轨迹TR1,TR2,…,TR8,其分别按照等距的径向距离Δkr关于k空间Sk的中点相互点对称地延伸。
通过在RF入射期间通过正弦或余弦函数调制两个选择梯度(即在选择性HF脉冲期间被发送的、用于选择的梯度脉冲),可以实现同心的环形轨迹:
在此n是取1和N之间的值的激励索引(或者说环形轨迹的索引),其中N是选择的轨迹的环的数量(也就是在按照图2的例子中适用N=8)。An是在第n个激励期间的梯度幅度的绝对值并且Tn是第n个激励(即在第n个环上的HF入射)的持续时间。公式(6)中的可能的共同的相位系数不失一般性省去以便将公式保持为尽可能简单。
公式(3)中的第一加权系数,即k空间平移速度,由此在环形轨迹的遍历期间是恒定的:
第二加权系数,即密度补偿系数
可以对于等距的环形轨迹通过在k空间中的相邻的环之间的径向距离Δkr的倒数值来估计。这一点在后面借助图5更详细解释。密度补偿系数在利用环形轨迹可变密度地占据k空间的情况下也就是环形轨迹非等距的距离的情况下才有意义。对于确定的环形轨迹n(其中1<n<N)则通过与两个最近的相邻环的平均的径向距离可以使用例如1作为密度补偿系数。
期望的激励轮廓P(x,y)按照本发明是旋转对称的并且,只要在极坐标中规定,仅取决于径向距离并且不取决于极角或方位角即,成立P(x,y)=P(r)。
通过公式(3)中的傅里叶变换从笛卡尔到极坐标的坐标变换首先获得:
在此是位置矢量的方位角并且是波矢量 的方位角并且Jm是第m阶贝塞耳函数。在假定期望的激励轮廓是旋转对称的情况下,可以将关于空间方位角的积分包括到和数中,并且仅剩下第0阶的项。由此可以将公式(3')简化写为:
因为在圆形轨线上波矢量k(t)=k的绝对值是恒定的,所以从该图示中直接得出,对于确定的圆形轨线(激励)在右侧上的积分也是恒定的并且由此b1幅度在单个环形轨迹期间也是恒定的。由此如开头解释的那样梯度线圈的接通延迟不再重要。
按照本发明的二维HF脉冲可以构造为具有多个这样的子脉冲的“综合的”HF脉冲(英语“composite RF-Pulse”),其中每个子脉冲由具有持续时间Tn的恒定幅度的“b1(n)”组成并且梯度在入射持续时间Tn期间正弦形或余弦形地以周期Tn这样振荡,使得在所属的激励k空间中遍历一个圆形轨线(即所属的第n个环形轨迹)。
具体地,为了激励定义的体积或激励轮廓,必须在单个子脉冲的设计之前首先更精确确定环形轨迹。以下在此作为例子又假定具有直径d的圆柱形棒体积,其例如可以在导航技术中被使用。
对于应用者有关的和由此由其规定的参数是棒的直径d和第一边激励SL(Sidelobe,旁瓣)的距离SSL。边激励(Seitenanregung)是由于k空间的离散扫描而不可避免的。应用者这样大地选择第一边带(Seitenband)SL的距离,使得其刚好处于检查对象的外部并且由此没有组织被激励。如果检查对象是成年人,则对于距离的直径的典型值是d~25 mm并且在上述导航应用中对于第一边带的距离是SSL~400mm。
从公式(7)中可以看出,复数值的b1场和期望的激励轮廓P(r)形成傅里叶变换对。
从特征P(r)=0对于r>SSL和尼奎斯特理论中根据以下得出k空间中扫描的最小距离
激励的规定的直径d通过扫描的k空间的直径2·krmax=2·N·Δkr确定环形轨迹的数量N:
因为期望的激励轮廓的傅里叶变换对于kr>krmax是零,所以其频谱被限制。由此可以应用尼奎斯特理论并且关系(9)是直接的结果。
从公式(8)和(9)中获得环形轨迹的最小数量N作为规定的参数d(激励轮廓的直径)和SSL(第一边激励的距离)的函数:
图3为了解释示出了通过具有直径d和具有按照距离SSL的边激励SL棒形激励轮廓EP的、一层的位置空间图像(即横截面)。位置空间图像是具有240mm直径的球形模体的图像。在FLASH序列中将传统的激励脉冲通过具有按照本发明的具有同心的环形轨迹的激励轨迹的二维选择性HF脉冲代替。环的数量设置为N=8并且第一边环的距离设置为等于64mm。这相应于在图2中示出的轨迹,即,图2和3也示出了两个傅里叶域中激励的关系(图3中的位置空间中的激励轮廓和图2中的所属的激励k空间Sk)。按照公式(10)由此预计为d=2·64/8mm=16mm的中央激励轮廓的直径。在按照图3的测量中确定的尺寸与该计算一致。第二边带具有128mm的半径或256mm直径并且由此已经位于球外部。在拍摄中使用具有在读出和相位编码方向上的各256个像素的正方形矩阵。其余的序列参数在图3中图像的建立中是TR=50ms,TE=6.3ms和翻转角30°,共振频率是123MHz。
当基于激励轮廓的直径d和旁瓣距离SSL确定了环形轨迹的数量N和环形轨迹n=1到N的半径krn=n·Δkr时,可以借助公式(1)和(2)确定对于每个环形轨迹的最大的梯度幅度An和脉冲持续时间Tn。通过公式(6)的、关于四分之一周期的积分获得在半径krn=n·Δkr和梯度幅度An的绝对值以及第n个环形轨迹的持续时间Tn之间的以下关系:
为了将激励脉冲的持续时间最小化,例如将公式(11)中Tn的选择为这样小,使得规定的最大梯度幅度Gmax和规定的最大允许梯度增长率Smax(英语“slewrate”)刚好不被超过。最大允许的梯度幅度Gmax通常通过梯度系统的性能确定,而在最大允许的梯度增长率Smax的规定时优选考虑梯度系统的性能和被检查的人的可能的神经刺激。
图4示出了单个环形轨迹的HF子脉冲RFn的设计的例子。图示在此以序列图或脉冲图的一个小的部分的形式进行。在该脉冲图中以通常的方式在上面的时间轴上关于时间t示出了高频脉冲以及在之下的不同时间轴上示出了与之协调地待接通的梯度脉冲。在此仅示出对于在该层内部的选择性激励有关的在x和y方向上的梯度脉冲GPx,n,GPy,n(梯度曲线的幅度An)。在梯度下面分别画出了通过涉及的梯度曲线对自旋起作用的、积累的0阶矩(在各自的梯度脉冲下面的面积的绝对值)F1,F2,根据作用方向具有不同的符号。
时间轴的零点与HF子脉冲RFn的中心重合,x梯度脉冲GPx,n相应于公式(6)具有正弦形曲线,y梯度脉冲GPy,n具有余弦形曲线,分别具有周期Tn(即环形轨迹的时间长度)。该选择可以在不失一般性的情况下任意进行。在时刻t=–Tn/2采用高频脉冲之前在两个梯度方向上分别示出一个预相位脉冲VPx,VPy。在时刻t=Tn/2高频脉冲RFn中断之后分别接通一个重聚相位脉冲RPx,RPy。
在y方向上的预相位脉冲VPy之所以是有意义的,是因为在时刻t=–Tn/2梯度场不能瞬间地从0上升到An。即其用于通过梯度线圈的电流的升高(英语“ramp up”)。相应地,在y方向上的重聚相位梯度RPy用于梯度电流的下降。两个梯度脉冲VPy,RPy的总矩对于每个HF子脉冲,即对于每个环形轨迹,是0(-F2+F2=0)。由此实现,由前面的子脉冲(预相位梯度)或者说由该子脉冲或前面的子脉冲(重聚相位梯度)在横向平面中激励的(静态的)自旋,由于这些梯度而没有获得任何相位。在具体在图4中示出的实施方式中选择对于预相位梯度或重聚相位梯度VPy,RPy的三角形(或梯形)曲线,因为由此梯度的持续时间在给出的最大幅度Gmax和增长率Smax的情况下可以被最小化。在此存在大的自由度。利用四分之三正弦波例如可以降低梯度噪音。此外可以从三个循环中选择梯度曲线,对于所述循环,一阶矩也消失,以便也可以不影响在y方向上流动的自旋。
在x方向上的预相位脉冲VPx用于将k空间轨迹对中心。其0阶矩F1等于在HF激励期间正弦形梯度脉冲GPx,n的四分之一周期之下的0阶矩F1。在x方向上的重聚相位梯度RPx具有相同的绝对面积F1和相反的符号。其将k空间轨迹带回到原点。在x方向上的梯度的总矩(具有和没有HF入射)是0。由此已经通过前面的子脉冲被激励的(静态)自旋保持不受影响。在图4中示出的实施方式中为了最小化所需的时间又选择三角形曲线。替换地,例如在此也可以是梯形曲线或半正弦波。如果二维选择性HF脉冲例如作为激励脉冲在扰相的梯度回波序列中被采用,其中所有自旋在激励之前处于纵向方向上,则例如可以将在x方向上的预相位梯度VPx对于激励的第一子脉冲省去,因为梯度对具有纵向取向的自旋没有影响。
如上所述,在每个单个子脉冲中b1幅度是恒定的。N个子脉冲的值b11,b12,…,b1N的计算在预先给出的、期望的旋转对称激励轮廓P(r)的情况下借助公式(7)进行。但是该公式(7)仅规定每个子脉冲的相对b1幅度或者说比例。但这是足够的。b1幅度的物理值(例如以μT为单位)可以由此如下确定:
共振频率位于HF(子)脉冲的带宽内部的自旋的磁化在脉冲的末端以如下翻转角从均衡位置翻转:
在此t0是子脉冲的接通时间并且回磁比(英语“gyromagnetic ratio”)γ是物理常数,其取决于被激励的核。对于质子,其值为γ=2π·42,57MHz/T。
因为b1幅度在子脉冲期间是恒定的,所以从公式(13)对于第n个子脉冲的翻转角值直接得到:
αn=γb′1nTn (13b)
在此b1n是借助公式(7)(具有等号替代比例符号)计算的用于第n个子脉冲的值,b1n′=c·b1n是第n个子脉冲的首先未知的物理的b1幅度。于是在假定HF脉冲的总持续时间相对于弛豫时间是短的(即自旋在HF入射期间没有在相关的周围范围内又弛豫)的情况下成立,综合的脉冲的翻转角α等于第n个子脉冲的翻转角值αn的和:
i在此仅是关于所有子脉冲的另一个动态变量,N仍是子脉冲的总数。
从公式(14)中按照如下获得第n个子脉冲的翻转角值αn:
最后的等号成立,因为在b1n′和b1n(b1n′=c·b1n)之间的比例常数c取决于子脉冲n。由此可以将在分母中的常数c提到和数之前并且约去。
公式(15)右边的所有项的计算在前面已经示出。从翻转角αn于是可以借助公式(13b)计算第n个子脉冲的幅度的所求的物理值:
如上讨论的,子脉冲的幅度是以两个系数,即,k空间横向速度(k-Raum-Transversalgeschwindigkeit)和密度补偿系数加权的、期望的径向轮廓的一维傅里叶变换。公式(7)中的积分在此可以根据期望的激励轮廓解析地或数值地计算。在实现的实施方式中选择高斯激励轮廓,
其中关于期望的直径d这样确定变量a,使得在高斯函数下面的总面积的90%位于直径d内部。
第一加权系数,即,k空间横向速度等于梯度幅度的绝对值An:
梯度幅度的绝对值An在整个环形轨迹期间是恒定的并且关于公式(8)至(11)作为预先给出的参数d,SSL的函数规定。
第二加权系数、即密度补偿系数的计算在以下以具有等距的环距离的同心环形轨迹为例解释。
根据公式(5),密度补偿系数可以通过与k空间扫描点(kx,ky)相关联的面积ΔA(kx,ky)和由该面积ΔA包围的轨迹的长度Δl(kx,ky)的商来近似。在每个实际的实现中HF幅度被数字化。Mn是第n个子脉冲的HF扫描点(英语“样本”)的数量。于是与每个扫描点相关联的第n个轨迹的长度Δln是第n个环形轨迹的直径除以扫描点的数量:
与单个样本相关联的面积ΔA如在图5(阴影面积)中可以理解的那样是该长度Δln乘以在两个相邻的环形轨迹之间的距离Δkr:
ΔAn=ΔlnΔkr(20)
对于第n个环形轨迹的密度补偿系数ρn由此是
在图(2)和(5)中示出的例子中(等距的环距离)密度补偿系数对于所有子脉冲是相同的。在不等距的环距离情况下(例如参见图10)可以对于一个环形轨迹近似地将公式(20)和(21)中的径向的环距离Δkr通过与两个相邻环形轨迹的平均距离代替。
旋转对称的激励轮廓的中心可以以离梯度系统的等中心的距离(δx,δy)被移开,方式是,利用与瞬时的梯度矩成线性比例的相位系数修改b1场:
在此b1n是对于在等中心中的激励的第n个子脉冲的取决于时间的b1幅度并且(kx,n(t),ky,n(t))是在遍历第n个环时在时刻t的波矢量。在此,HF脉冲与现有技术中已知的二维选择性HF脉冲并无不同。
对于按照本发明的脉冲排列的构造,其中二维选择性HF脉冲通过这样的上述子脉冲的综合来代替,存在不同的可能性。
在最简单的实施方式中简单地在时间上先后执行子脉冲。在单脉冲之间的时间距离在此在原理上可以自由选择。由于由前面的子脉冲影响的、自旋的不可避免的T1-和T2-弛豫,并且由于已经被激励的自旋的信号的去相位,例如由于在执行后面的子脉冲期间的B0非均匀性,在大多数情况下人们优选尽可能短的总脉冲持续时间。特别地优选试图将其中不进行HF入射的时间段最小化。
图6示出了为了采集图3中的图像(为此也参见图2中的环形轨迹)被使用的、按照本发明的二维选择性HF脉冲RFE的仿真的脉冲排列PA的脉冲图,以及所属的在x和y方向上的梯度脉冲GPx,GPy,(类似于图4中、但是现在是完整的二维选择性HF脉冲RFE)。在上面的轴上示出了以任意单位([a.u.])的HF子脉冲的b1幅度。在位于下面的轴上同样以任意单位示出了梯度幅度。时间轴的单位分别是μs。选择性HF脉冲可以被划分为八个子脉冲RF1,RF2,RF3,…,RF8。具有最大的k空间半径krmax的第一子脉冲的持续时间T这样短地选择,如在给定的最大允许的梯度幅度Gmax和给定的最大允许的梯度增长率Smax情况下所可能的那样。其他子脉冲RF2,RF3,…,RF8的持续时间T是非梯度限制的。而是用于遍历所有环形轨迹TR1,TR2,TR3,…,TR8的持续时间T(和由此各自的矩形的RF脉冲RF1,RF2,RF3,…,RF8的持续时间T)对于所有子脉冲相同地选择。借助在二维选择性螺旋脉冲中的术语,在此可以提及具有恒定的角速度的实施(英语“constant angular rate”)。
如上所述地计算各矩形的子脉冲RF1,RF2,RF3,…,RF8的(恒定的)幅度。
图7示出具有按照本发明的二维选择性HF脉冲RFE'的另一个实施变形的脉冲排列PA'的相应的(仿真的)脉冲图,以及在x和y方向上的所属的梯度脉冲GPx',GPy'。与按照图6的变形的主要区别在于,在此每个单脉冲RF1,RF2,RF3,…,RF8的持续时间被单独地计算,并且具体来说这样计算,使得在各自的第n个环形轨迹TR1,TR2,TR3,…,TR8的给定的k空间半径krn=(N-n+1)Δkr情况下给定的最大允许的梯度幅度Gmax和最大允许的梯度增长率Smax刚好不被超过。又借助二维选择性螺旋脉冲中的术语可以将这称作具有恒定的(梯度)增长率的实施方式(英语“constant slew rate”)。
在两个在图6和7中示出的实施方式中分别将第n个子脉冲的重聚相位梯度与第n+1个子脉冲的预相位梯度综合。两个梯度的这样的综合意味着,替代在时间上先后执行两个梯度,仅接通一个梯度(梯度脉冲),其0阶矩等于单梯度的0阶矩之和。综合的目的主要是节省时间。相邻的梯度脉冲的综合的另一个优点在于,由此降低梯度噪音、涡流和梯度线圈和梯度放大器的发热。
在按照图7的实施方式中,其中进行HF入射的时间间隔T1,T2,T3,…,T8是最佳的。每进一步的缩短导致超过最大允许的梯度增长率Smax或最大允许的梯度幅度Gmax。于是二维选择性RF脉冲RFE'或整个脉冲排列PA'的总持续时间的进一步缩短仅通过缩短没有HF入射的、其间进行重聚相位和预相位梯度的时间才可以。
图8为此示出了具有按照本发明的二维选择性HF脉冲RFE″的实施变形的脉冲排列PA″的相应(仿真的)脉冲图,以及在x和y方向上的所属的梯度脉冲GPx″,GPy″,其中在具有RF入射的间隔之间的时间再次降低。这通过如下实现,即,k空间中的相邻的环形轨迹以不同的旋转方向被遍历。例如具有奇数索引的子脉冲的旋转方向是逆时钟的而具有偶数索引的子脉冲的旋转方向是顺时针的。
该变形的优点可以根据图9来解释,其示出两个相邻的HF子脉冲RFn″,RFn+1″的放大部分和所属的在y方向上的梯度脉冲GPy″的部分(具体在图9中示出的子脉冲RFn″,RFn+1″例如可以相应于图8中的位于通过两条线表示的部分中的脉冲RF3,RF4,但是原该理对于所有其他时间上先后跟随的子脉冲RFn″,RFn+1″都是相同成立的)。
如在图9中明显可以看出的,在时间间隔Tn和时间间隔Tn+1之间的y梯度由两个三角形梯度组成。第一梯度具有幅度An和斜坡持续时间RTn,第二梯度具有幅度An+1和斜坡持续时间RTn+1。在此成立:
AnRTn+An-1RTn+1=0 (23)
两个三角形梯度的面积的绝对值(即0阶梯度矩的绝对值)由此相同。由此位于横向平面中在第n个激励脉冲之后的(静态)自旋由于这些梯度而没有获取任何相位。初始条件,即,第n个y重聚相位梯度的和第n+1个y预相位梯度的矩分别是零,可以在分别相反的延伸的相邻的环形轨迹的情况下由于幅度An和An+1的不同符号而通过更弱的条件代替,即,综合的梯度的总矩(参见等式(23))消失。这是在该方法变形中的附加的时间收益的原因。
在图10中示出了k空间Sk中按照本发明构造的发送轨迹的另一个实施例(类似于图2)。在此整个发送轨迹同样由围绕k空间中心的同心的环形轨迹组成。然而各个环形轨迹TR1,TR2,TR3,…不再具有与在按照图1的实施例的情况下相同的径向距离。而是环距离Δkr1,Δkr2,Δkr3,…朝着k空间中心逐渐减小。特别地k空间中心可以被这样过扫描,即,靠近中心将环距离选择为小于通过公式(8)右边给出的值。通过k空间中心的更密集扫描预计降低在期望的激励轮廓P(r)外部的不期望的激励。为了避免通过k空间中心的更密集扫描延长总脉冲,可以欠扫描k空间外围,即,将环距离Δkr1,Δkr2,Δkr3,…选择为大于通过公式(8)要求的。
在图11中示出了k空间Sk中按照本发明构造的发送轨迹的实施例(类似于图2中),其中发送轨迹仅在k空间Sk的内部区域中(例如在至最大半径krmax的一半的半径区域Δkri中)由围绕k空间中心的同心环形轨迹形成。在k空间Sk的外围半径区域Δkra中取而代之采用螺旋轨迹TRS。同心环形轨迹相对于常规的螺旋形轨迹的优点,如上所述和后面还要结合测试示出的,是更大的鲁棒性,特别是相对于梯度延迟。但是螺旋轨迹的一个优点是其快速性。为了利用这两个优点,由此在该变形中首先利用子脉冲经过具有(不完整的)螺旋轨迹的激励k空间的外围区域并且然后利用分别实现一个同心环形轨迹的子脉冲经过k空间的中央区域。
原理上当然也可以是相反的布置,即,外部使用同心的环形轨迹而内部使用螺旋轨迹。但是在图11中示出的变形就此而言是有利的,因为在外围k空间区域中的信息比来自于k空间中心的信息更不重要。这也是为什么在另外示出的变形中(参见图2和10)总是最后覆盖k空间中心的原因。
如上所示,在环形轨迹上被入射的每个子脉冲的b1幅度是恒定的,与这在公知的具有EPI或螺旋形轨迹的二维选择性HF脉冲的情况下不同。在入射的HF能量和实际上施加的梯度场之间的不期望的延迟由此仅导致仅在子脉冲的开始和仅在子脉冲的结束时的错误。这导致在存在梯度延迟时间的情况下选择轮廓的极好的鲁棒性。
图12至14示出该鲁棒性的试验验证。为此分别在二维选择性FLASH序列中将激励脉冲通过二维选择性HF脉冲代替。在所有三个图12至14的上面行中该二维选择性HF脉冲分别是具有EPI轨迹的传统的脉冲,在下面的行中是具有按照本发明的同心的环形轨迹的综合的HF脉冲。
中间的和右边栏中的图像在所有图中分别与左边栏中的所属的图像通过如下相区别,即将入射的RF场首先人工地以15μs或45μs相对于施加的梯度场延迟(延迟时间分别利用Del表示)。在测量中使用的模体是以模体液体填充的球。第一边带的距离在所有情况中是128mm并且由此在模体外部。
期望的激励轮廓在具有EPI轨迹的所有激励情况下分别是方形棒和在具有同心环形轨迹的所有激励情况下分别是圆柱体,以便分别将激励轮廓的形状匹配于轨迹的几何基本结构并且由此保证可比性。在按照本发明的激励情况下环形轨迹的数量N是8并且圆柱体的直径相应地是大约32mm。这样选择EPI脉冲的参数,使得方形棒的边长(至少理论上)也为32mm。
读出视野(英语“readout field of view”)在所有情况中分别是256mm并且垂直于棒或圆柱体轴定向。图像在西门子3T MAGNETOM Verio中被拍摄。
在图12中模体位于MR设备的等中心中。在图13和14中将模体以及激励的期望的位置以50mm在x方向或y方向上移动。棒或圆柱体轴在所有情况中分别在磁体的z方向上(也就是在B0场的方向上)被示出。
在所有三个实验中可以看出,EPI轨迹的轮廓在15μs延迟时间的情况下是显著的并且在45μs的情况下是完全无序的,而同心环形轨迹的轮廓很大程度地得到保持。只有在45情况下才看出轻微的拖尾效应。
在偏心激励(图13和14)情况下延迟时间引起激励相对于(位于读出视野中心中的)期望位置的轻微偏移。然而该误差在许多应用中、特别是在导航应用中是无害的,因为在此例如待激励的棒形体积例如不需要一定是定位在肝尖上,以正确探测呼吸运动。
此外该方法的特征也可以有利地用于调节和/或校准延迟时间。例如为此可以一直改变人工插入的附加延迟时间,直到测量的激励轮廓恰好位于期望的位置上并且由此系统内在的(首先未知的)延迟时间被精确补偿。该校准是必须一次性还是必须在活体内地进行,取决于延迟时间的原因。
为此的方法的简化流程图在图15中示意性示出。在步骤I中首先借助按照本发明的方法进行精确定义的激励轮廓的二维选择性激励。然后在步骤II中向用户显示基于该在步骤I中进行的测量重建的图像数据。用户可以判断,激励轮廓是否精确位于期望的位置上。如果不是,则他在步骤III中人工地改变梯度延迟时间。在此可以在不同的空间方向上改变梯度。然后在步骤I中进行具有相同的激励轮廓的新的测量并且在步骤II中又显示在此获得的图像。如果用户在步骤II中对结果满意,则磁共振系统是调整的(或按照值的相应协议进行校准)并且然后可以在步骤IV中进行本来的测量。替代通过用户进行的视觉监视可以在更合适的变形中也进行图像数据的自动分析。例如可以借助图像识别软件确定图像数据中激励轮廓的位置并且基于此然后自动地在步骤III中进行梯度延迟时间的调整等,直到最后按照图像数据补偿了延迟时间。
最后再次指出,前面描述的详细方法和装置是实施例并且原理也可以在宽的范围内由专业人员改变,而不脱离本发明的通过权利要求规定的范围。例如控制序列确定装置22替代在终端上也可以作为控制装置10本身的部分实现,特别是也可以作为测量控制单元15的组件。同样,控制序列确定装置也可以在例如经过网络NW与磁共振设备1相连的单独的计算系统上实现。为完整起见还要指出,不定冠词“一个”的使用不排除,涉及的特征也可以多重存在。同样概念“单元”或“模块”的使用也不排除,它们由多个必要时也可以在空间上分布的组件组成。