专利名称：用于磁共振成像的三维空间选择性激励的序列设计方法
技术领域：
本发明涉及用于磁共振成像的三维空间选择性激励的序列设计方法。
背景技术：
当人体组织等物质置于均勻静磁场(Btl场，即主磁场)下时，组织内的质子磁矩方向将倾向于与该Btl场一致，且沿着该方向以拉莫频率进动，令此时的磁化矢量方向为ζ轴正向，幅值为Mtl，称组织内的质子处于平衡态。如果此时施加方向与Btl场垂直(即xy平面内)、且接近拉莫频率的射频场(B1场)，组织内部的磁化矢量将向xy平面旋转，其轨迹处于与B1场垂直的平面内，其旋转的角度称为受激翻转角，此时的B1场称为激发场。撤去激发场后，组织内部的磁化矢量的xy分量将逐渐减弱、ζ分量将逐渐增强，并逐渐恢复为原平衡态磁化矢量Mtl，该过程称为弛豫。在弛豫过程中，之前受激发的质子将释放出以拉莫频率为载波频率的信号，从而被磁共振系统采集并处理成像。由此可见，人们赖以分析组织结构、功能状态的信息来自于这个受激区域内质子的弛豫信号。只有当受激区域内的所有质子从同样的状态(即同样的翻转角)开始弛豫， 该区域内的组织结构和功能状态才能在磁共振成像中得以准确的体现。此为受激翻转角的均勻性要求。影响翻转角激发均勻性的主要因素为激发射频场的均勻性。随着人们对磁共振成像高分辨率的追求，主磁场场强由ο. 35T逐渐提高到3T甚至11T，相应的射频场频率随之由15MHz提高至128MHz乃至470MHz，射频场的波长逐渐与人体尺寸相近，场的干涉和衰减效应显著，射频线圈的设计与改进已经不足以改善加载负载(即人体)后的射频场均勻性。 因此激发脉冲序列的设计成为必需。裁剪脉冲法是利用小翻转角前提下、目标场分布与射频场之间的傅里叶变换关系，来计算对应的激发脉冲序列的方法。它可以补偿射频场不均勻的影响，并且均勻激发任意形状的目标区域。但是该方法用在多维、尤其是三维上的脉冲序列的设计时，所得到的脉冲序列常常占用过长的激发时间，使得该激发过程极易受到偏共振效应的影响，而无法被实际应用。近期发展的并行激发技术(parallel transmission或parallel excitation)，采用多个射频线圈或同一个射频线圈上的多个发射信道发射射频脉冲，该技术利用各射频线圈或各发射信道的射频场分布的变化性所带来的K空间内的点扩散作用，来实现对K空间的欠采样，从而加速对目标场的激发，减小脉冲序列时长。根据现有技术，基于并行激发技术的三维裁剪脉冲法主要应用于选层激发。它使用立体回波(echo-volumar)轨迹(或称spoke或fast-kz轨迹)。在这里，设定对应于层选方向的K空间方向为“kz”，垂直于kz方向的的K空间平面为kx_ky平面；立体回波轨迹在 kx-ky平面上对应有多个相位编码点(或称spoke点)，每个spoke点确定一条等长的沿kz 方向穿越该点的spoke线。K空间走线方式为从原点出发，以最短距离经过所有的spoke线， 最终返回原点。spoke线的长度通常由激发目标的分辨率决定。传统的设计方法的主要目的在于寻找1^_1^平面上的spoke点的最小集合，在确保激发准确度的同时缩短激发时间。 传统方法有四种，分别是傅里叶方法(Fourier-based method)、逆方法(inversion-based method)、强制稀疏方法(sparsity-enforced method),禾口联合设计方法(joint design method) 0傅里叶方法直接对层内(xy平面上)的目标场分布进行傅里叶变换，并把spoke 点安置于傅里叶变换结果的最大幅值区域；该方法操作简便直接，但是无法计入射频场不均勻性的影响(或称空间敏感度，spatial sensitivity)。逆方法基于激发K空间内的射频场能量分布与质子受激情况之间存在的傅里叶关系，导出以寻找令激发结果最接近目标场的spoke点分布为目标的优化问题，在该优化问题的结果集中，寻找其中权重相对较高的点作为最终的spoke点集合；该方法弥补了傅里叶方法的不足，能够在优化问题中加入空间敏感度因素的考量，但是其结果通常是一个密集分布的高权值集合，而遍历所有这些高权值的spoke点，将导致过长的激发时间。强制稀疏法与逆方法计算原理类似，但是在计算过程中加以改进、促使计算结果的稀疏化，从而找到一个数量适中的高权值点集合。联合设计法不同于前述的三种方法之处在于，其在计算过程中同时考虑spoke点集选取的优化和射频脉冲设计的优化。强制稀疏法和联合设计法能够得到最优的spoke点集，从而最大程度的优化激发K空间轨迹。但是，上述所有方法都无法直接扩展到除了选层激发外的其他三维体选择性激励的应用上。目前尚没有针对泛三维体选择性激励的基于并行激发技术的K空间轨迹设计策略及其序列设计方法，以及可实现泛三维体选择性激励的激发脉冲序列生成装置。此外，上述所有方法都局限于立体回波类型的K空间轨迹。但是，高场下针对泛三维体的选择性激励将是非常有意义的应用。首先在于，小范围视场(或称小范围FOV，FOV即field of view)成像中，由于成像K空间的欠采样，虽然缩短了成像时间，但是也带来了包含其他部位信息的伪影；通过三维体选择性激励，选择性的激发该小范围FOV内的质子，则成像结果中将不再存在来自其他部位受激质子干扰的伪影了。同样的，在磁共振血管造影术中，通过选择性的激发或者反转来实现那些流入待考察的血管之前的那部分血液内的氢质子的自旋标记，则可以在不加造影剂的前提下得到关于这个待考察的血管的很好对比度的成像结果，从而避免因使用造影剂而造成对人体的伤害。此外，在脑功能成像中应用的单平板成像(single slab imaging)，是以厚平板内的单次均勻性激发为前提，利用三维傅里叶成像法进行磁共振成像，该方法能够在提高成像分辨率的同时保证分辨率的各项同性，大大提高了成像信噪比。
发明内容
本发明所要解决的问题是提供一种用于磁共振成像的三维空间选择性激励的序列设计方法，通过它可得到实现磁共振系统的泛三维体选择性激励的序列。本发明解决其技术问题所采取的技术手段是本发明一种用于磁共振成像的三维空间选择性激励的序列设计方法，其特征是，包括如下步骤 (1)确定磁共振系统的成像空间内的期望三维激发目标，并且通过&+映射方法计算出所述磁共振系统的射频线圈的各发射信道所对应的三维空间敏感度分布情况；(2)根据所述期望三维激发目标和所述三维空间敏感度分布情况，计算出相应的三维激发K空间轨迹容器的范围；(3)根据所述期望三维激发目标和所述三维激发K空间轨迹容器的范围，确定相应的三维激发K空间轨迹；(4)根据所述三维激发K空间轨迹和所述磁共振系统中的梯度驱动单元的梯度脉冲的最大幅值和最大切换率，计算出所述梯度驱动单元的梯度脉冲波形；(5)根据所述期望三维激发目标、所述三维空间敏感度分布情况和所述梯度脉冲波形，计算出所述射频线圈的各发射信道的射频脉冲的包络波形；由所述梯度脉冲波形、所述射频脉冲的包络波形和所述射频脉冲的载波频率构成所述三维空间选择性激励的序列。进一步地，本发明对步骤( 所得到的序列采用变速选择性激发方法或可变转换率螺旋型轨迹设计方法进行优化。进一步地，本发明所述步骤( 中，按以下步骤计算出相应的三维激发K空间轨迹容器的范围1)根据所述期望三维激发目标，生成用于表征所述期望三维激发目标的第一向量；根据所述射频线圈的各发射信道的三维空间敏感度分布情况，生成满足以下关系式 ⑴的第一矩阵Afull = [S1A,…，SlA] (I)式(I)中，Afull即所述第一矩阵，S1是表示所述射频线圈的第1个发射信道的三维空间敏感度分布情况的对角化矩阵，A是表达实际空间与K空间之间关系的矩阵且满足以下关系式(II) A(m, η) = iye'r-'k"(II)式(II)中，A(m，n)表达矩阵A的第m行、第η列的元素的数值，i是虚数符号，γ 是磁旋比，^是对实际空间离散化得到的第m个点的空间坐标，kn是对K空间离散化得到的第η个点的K空间坐标；2)根据所述第一向量和所述第一矩阵，计算得到用于表征所述射频线圈的所有发射信道所对应的激发K空间的能量分布的第二向量；3)将所述第二向量根据以下关系式(III)，分离为用于表达所述射频线圈的各发射信道所对应的激发K空间的能量分布的向量组Wfull = [W1, . . . , ffL]T (III)式(III)中，Wfull即所述第二向量，W1是表达所述射频线圈的第1个发射信道所对应的激发K空间的能量分布，T是矩阵转置符号；4)将所述向量组的幅值加和，并将该加和结果归一化为结合了所述射频线圈的各发射信道所对应的激发K空间的能量分布情况的统一的能量权重分布；5)根据所述能量权重分布，确定出所述K空间轨迹所需遍历的高权重区域；6)根据所述高权重区域内的主体部分，确定一个规则封闭的形状并将该形状作为 K空间轨迹容器的参考范围；7)根据所述三维空间敏感度分布情况，计算出所述射频线圈的各发射信道所对应的三维激发K空间内的点扩散作用范围；根据所述点扩散作用范围削减所述K空间轨迹容器的参考范围，得到K空间轨迹容器的范围。进一步地，本发明所述步骤(3)中，按以下步骤确定相应的三维激发K空间轨迹a)根据所述期望三维激发目标，确定三维K空间轨迹的类型；b)根据所述期望三维激发目标、所述三维K空间轨迹的类型，确定三维K空间轨迹的径向采样间隔、纵向采样间隔；c)根据所述三维K空间轨迹的类型、径向采样间隔和纵向采样间隔，以及所述K空间轨迹容器的范围，确定相应的三维激发K空间轨迹。进一步地，本发明所述步骤a)中，所述三维K空间轨迹的类型为立体回波型或堆栈螺旋型。与现有技术相比，本发明的有益效果是(1)由于引入了用于限定三维K空间轨迹分布范围的K空间轨迹容器，本发明的序列设计方法能够实现包括选层激发的泛三维体选择性激发的序列设计；与此同时，该分布范围的限定提高了 K空间轨迹的采样效率，能够有效的减少激发序列的时长、以及减小射频脉冲幅值。(2)在引入的K空间轨迹容器的范围内设计K空间轨迹，允许任意的K空间轨迹类型的使用，尤其是堆栈螺旋型轨迹的使用，能够大大减小激发所需要的射频脉冲包络波形的最大幅值(3)在引入的K空间轨迹容器的范围内设计K空间轨迹，能够充分利用并行激发技术的加速因素，能够在减少激发序列的时长的同时，维持激发效果；利用。( 利用变速选择性激发方法或可变转换率螺旋型轨迹设计方法对序列进行优化，能够进一步减小激发所需要的射频脉冲包络波形的最大幅值。减小激发所需要的射频脉冲包络波形的最大幅值，即减小射频能量，不仅能够降低SAR值，提高磁共振扫描的安全性，还有助于在射频能量强烈衰减的低信号区域内实现接近期望的激发目标。以上都是现有技术所不具备的特征。
图1是磁共振系统的结构示意图；图2是本发明的序列设计方法的流程图；图3 (a)为传统的三维K空间轨迹中的回波立体型轨迹，图3 (b)为本发明的K空间轨迹容器，图3 (c)为利用图3 (b)的K空间轨迹容器对图3 (a)的回波立体型轨迹优化后的轨迹；图4是图3 (c)的A-A平面映射图；图5是图3 (c)的B-B平面映射图；图6 (a)为传统的三维K空间轨迹中的堆栈螺旋型轨迹，图6 (b)为本发明的K空间轨迹容器，图6 (c)为利用图6 (b)的K空间轨迹容器对图6 (a)的堆栈螺旋型轨迹优化后的轨迹；图7是图6 (c)的A-A平面映射图；图8是图6 (c)的B-B平面映射图。
具体实施例方式本发明的序列设计方法适用于包含有多发射信道射频线圈(如阵列线圈)的磁共振系统。参见图1，本发明所适用的磁共振系统一般包括一个计算机100、磁共振扫描系统 102、可调整被测试对象位置的扫描床112、Bl+映射处理器114、脉冲序列生成器116、数据采集单元118、数据处理单元120、梯度驱动单元122、射频驱动单元124、生理信息采集控制器126、扫描室接口单元128、被测试对象定位系统130以及外置显示器132。磁共振扫描系统102主要包括用以产生主磁场的主磁体104、用以产生梯度磁场的梯度线圈106和用以产生射频场的射频线圈(例如体线圈108或阵列线圈110)、以及可调整被测试对象位置的扫描床112。计算机100则包含正当商业途径可购买到的处理器和存储器，以及一个正当商业途径可购买到的操作系统。该计算机100拥有与磁共振系统中其他各模块或设备相连的接口，使得计算机100的指令得以在该磁共振系统中传送并得以执行。计算机经由该接口与磁共振系统中的四个模块单元进行实时的数据交换，该四个模块单元分别是=B1+映射处理器114、脉冲序列生成器116、数据采集单元118和数据处理单元 120。脉冲序列生成器116获取计算机100的指令并依据指令来控制梯度驱动单元122 和射频驱动单元124。脉冲序列生成器116依据指令产生用于实现预期扫描的梯度脉冲波形和射频脉冲波形的具体的数字指令，并将它传送给梯度驱动单元122和射频驱动单元 124。梯度驱动单元122依据指令产生梯度脉冲并施加在梯度线圈106上，以生成实现空间编码的梯度磁场(^、 和&。在一个激发过程中，射频驱动单元IM依据指令产生射频激发脉冲并施加在射频线圈(例如体线圈108或阵列线圈110)上、以生成期望中的射频磁场； 射频场与梯度磁场共同作用来完成对既定空间位置、既定空间范围内的既定翻转角幅度的激发。用以生成该特定射频场与梯度场的射频脉冲和梯度脉冲构成一个激发脉冲序列。在一个成像过程中，射频驱动单元1 依据指令通过射频线圈(体线圈108或阵列线圈110) 接收磁共振信号，并将数字化后的信号数据传递给数据采集单元118 ；梯度场同时作用来完成对该磁共振信号的空间编码。用以生成该特定梯度场的梯度脉冲构成一个成像脉冲序列。一个激发过程和一个成像过程构成一个完整的单次扫描过程。而一个激发脉冲序列和一个成像脉冲序列构成一个完整的单次磁共振成像脉冲序列。射频驱动单元IM包括一个或多个射频脉冲发射器。在激发过程中，射频脉冲发射器依据来自于脉冲序列生成器116的数据和指令来产生既定载波频率、既定包络波形的射频脉冲，并驱动射频线圈(如体线圈108或阵列线圈110)产生射频场；同时梯度驱动单元依据来自于脉冲序列生成器116的数据和指令产生既定的梯度波形，并驱动梯度线圈 106产生梯度场，对该射频场进行空间编码；经过既定的执行时间，实现一个激发过程。射频驱动单元124同时包含一个或多个射频信号接收通道，每个接收通道包含一个能够放大经由射频线圈108或110接收到的磁共振信号的射频信号放大器。在成像过程中，被测试对象(如图1中的人体)中的质子受到激励后衰减产生磁共振信号。同时梯度驱动单元依据来自于脉冲序列生成器116的数据和指令产生既定的梯度波形，并驱动梯度线圈106产生梯度场，对该磁共振信号进行空间编码。射频线圈(如体线圈108或阵列线圈110)接收该信号，并传递给射频驱动单元124，射频驱动单元IM依据来自于脉冲序列生成器116的数据和指令来对信号进行放大、解调、滤波及数字化处理。此外，在磁共振扫描过程中，脉冲序列生成器116同时通过生理信息采集控制器 126来实时的接收和分析被测试对象的相关生理信息。生理信息采集控制器1 通过一系列的传感器实时采集被测试对象的心电信号、呼吸信号等生理信息。而脉冲序列生成器116 通过分析这些信号，能够实时的调整脉冲序列，来实现磁共振扫描与被测试对象的心跳和呼吸等生理运动的同步化，将扫描过程中伴随心跳和呼吸等生理运动的空间位移对磁共振成像结果的影响减至最小。脉冲序列生成器116还实时接收来自扫描室接口单元1 的信息。扫描室接口单元1 与扫描室内用以探测被测试对象位置和磁体系统状况的传感器相连，并提供可操作平台，显示被测试对象位置和磁体系统运行状况。与扫描室接口单元1 相连的被测试对象定位系统130接收来自扫描室接口单元128的操作平台的数据，并依此控制扫描床112 的移动，将被测试对象平移到指定位置。扫描室接口单元1 将定位后的被测试对象的位置信息传送给脉冲序列生成器116。在成像后续过程中，经由射频驱动单元IM数字化后的磁共振信号被传递到数据采集单元118。数据采集单元118在计算机100的指令下配合脉冲序列生成器116传递的时序信息实时的接收磁共振数据。数据采集单元118提供足够的缓存来存储至少一个扫描周期内的数据。此后，数据采集单元118将磁共振数据传递给数据处理单元120用以进一步的处理和分析。数据处理单元120接收来自数据采集单元118的磁共振数据，并按照计算机100 的指令进行处理。该处理包括对原始K空间数据的傅里叶变换，对重建所得图像的滤波操作，用于功能磁共振成像的重计算，以及用于运动目标或流体成像的重计算等等。经由数据处理单元120得到的重建图像被传递给计算机100并存储其中。成像结果同时被传送到操作台显示器102或外置显示器132。计算机100也可以通过网络将图像传送到网络上的其他设备。更具体地，在激发过程中，脉冲序列生成器116依据计算机100的指令发送脉冲数据给梯度驱动单元122和射频驱动单元124，梯度驱动单元122据此生成梯度脉冲，并驱动梯度线圈106产生既定的梯度场。与此同时，射频驱动单元IM据此生成具有既定载波频率和包络的射频脉冲，并驱动体线圈108或阵列线圈110产生既定的射频场。通常该射频脉冲具有特定的载波频率和包络波形，从而配合梯度脉冲来实现对既定空间位置、既定空间范围内的质子的既定翻转角幅度的激发。在基于多发射信道射频线圈(如阵列线圈110) 的并行激发过程中，射频线圈各发射信道发送的射频脉冲具有同样的或相近的载波频率和时长，但是可以有不同的包络波形，结合各信道所对应的不同的射频场分布，最终产生期望中的激发结果。此结果包括被扫描体内受激发质子的空间范围、受激翻转角大小及其均勻性等等。因此，并行激发技术的实现，需首先获得射频线圈(如阵列线圈110)各发射信道所对应的射频场分布(或称各发射信道的空间敏感度)，并据此分布情况以及激发目标做脉冲序列设计。在本发明所采用的磁共振系统中，射频线圈(如阵列线圈110)各发射信道的空间敏感度的计算主要由B1+映射处理器114完成。进入B1+映射运行模式时，将首先通过多个特定的激发脉冲序列，得到多个与各发射信道空间敏感度相关的成像结果。更具体的，计算机100指令脉冲序列生成器116生成特定的激发脉冲序列波形，通过射频驱动单元 IM和梯度驱动单元122驱动射频线圈(如阵列线圈110)的某一发射信道和梯度线圈106 实现特定的激发；指令脉冲序列生成器116生成成像脉冲序列波形，通过梯度驱动单元122 驱动梯度线圈106、射频驱动单元IM控制具有均勻空间敏感度的接收通道的射频线圈(如体线圈108)接收磁共振信号；指令数据采集单元118和数据处理单元120接收和处理数据从而得到成像结果，传递回计算机100并存储于计算机100 ；经过如上所述的多个特定的扫描过程，将该多个特定的激发脉冲序列激发得到的多个成像结果依次存储于计算机100 中。并由计算机100传送给B1+映射处理器114。B1+映射处理器114对这些成像结果做进一步的处理并计算出射频线圈各发射信道的空间敏感度情况。B1+映射处理器114计算完成后，将发送结束标记给计算机100，计算机100据此退出B1+映射运行模式，并指令B1+映射处理器114将射频线圈所有发射信道的空间敏感度数据传送到脉冲序列生成器116。脉冲序列生成器116根据来自计算机100的关于成像空间范围内的激发目标、关于梯度驱动单元122可生成的梯度脉冲的最大幅值和最大切换率、关于射频驱动单元IM将生成的射频脉冲的载波频率、以及来自B1+映射处理器114的射频线圈各发射信道的空间敏感度数据， 以及来自于扫描室接口单元128的关于被测试对象的定位信息数据，计算出梯度脉冲波形和射频线圈各发射信道对应的射频脉冲包络波形，并依据来自脉冲序列生成器116的脉冲序列时序指令，将计算结果及射频脉冲的既定载波频率等脉冲序列数据传送给梯度驱动单元122和射频驱动单元124 ；由梯度驱动单元122和射频驱动单元IM分别生成梯度脉冲与射频脉冲、并驱动梯度线圈106和射频线圈(如阵列线圈110)来完成一个基于并行激发技术的激发过程。本发明的三维空间选择性激励的序列设计方法由脉冲序列生成器116采用，能够生成实现基于并行激发技术的泛三维体选择性激发的脉冲序列。该脉冲序列应用于所述的基于并行激发技术的激发过程中，将能够实现磁共振系统的成像空间内的期望三维激发目标，即实现泛三维体的选择性激发。本发明的三维空间选择性激励的冲序列设计方法包括以下步骤确定磁共振系统的成像空间内的期望三维激发目标，其参数包括激发目标在三维空间内的位置、形状、翻转角大小、及分辨率要求，其中激发目标在三维空间内的位置将根据被测试对象的定位信息进行重定位；通过B1+映射方法得到磁共振系统的射频线圈的各发射信道所对应的三维空间敏感度分布情况；根据所述的期望三维激发目标、射频线圈各发射信道的三维空间敏感度分布情况，得到相应的三维激发K空间内的射频能量权重分布，并在此基础上，确定三维激发K空间轨迹容器的范围；将该K空间轨迹容器的范围作为限定K空间轨迹分布的范围，并根据所述的期望三维激发目标，选择K空间轨迹的类型、确定径向和纵向采样间隔，从而确定相应的三维激发K空间轨迹；根据所述的三维激发K空间轨迹、以及梯度驱动单元122所能生成的梯度脉冲的最大幅值和最大切换率约束，得到梯度脉冲波形；根据所述的期望三维激发目标、射频线圈各发射信道的三维空间敏感度分布情况、以及所述的梯度脉冲波形，计算出射频线圈各发射信道所对应的射频脉冲包络波形； 由所述的梯度脉冲波形、射频脉冲的包络波形和载波频率构成了实施三维空间选择性激励的序列。梯度驱动单元122和射频驱动单元IM根据该序列数据生成相应的梯度脉冲和射频脉冲，并驱动梯度线圈和射频线圈将其同时施加于扫描空间内，从而实现期望的三维空间内选择性激发的目标。此外，本发明的三维空间选择性激励的冲序列设计方法还包括，对所述序列采用变速选择性激发方法或可变转换率螺旋型轨迹设计方法进行优化，以进一步减小实现激发所需要的射频脉冲包络波形的最大幅值。本发明的三维空间选择性激励的冲序列设计方法将在下文中详细描述。首先完成如下定义磁旋比Y。
三维空间坐标向量r = [rx, ry, rz]T。沿x、y、z各方向的梯度脉冲波形组合G(t) = [Gx (t), Gy (t), Gz(t) ]τ0三维K 空间轨迹 k(t) = [kx(t)，ky(t)，kz(t)]T。梯度脉冲波形与K空间轨迹的关系为众⑴=-γ fj G(t)出，这里T为脉冲时长。各发射信道的三维空间敏感度分布情况，1 = 1，2，. . .，L，L为信道总个数。各发射信道对应的射频脉冲包络波形B1, 1 (t)，1 = 1，2，. . .，L。各发射信道对应的三维激发K空间的射频能量分布W1GOd = 1,2,...,L0横向磁化矢量Mxy = Mx+iMy,这里Mx为χ方向磁化矢量、My为y方向磁化矢量。纵向磁化矢量Mz，令质子受激发前或稳态时其值为M。。参见图2，本发明的三维空间选择性激励的冲序列设计方法中，首先需要确定期望实现的三维激发目标、以及通过B1+映射得到关于射频线圈各发射信道的三维空间敏感度分布情况。这里的三维激发目标具体包括，该激发目标在三维空间内的位置、形状、 分布范围，该激发目标内部的翻转角大小、均勻性，以及该激发目标的分辨率要求；其中激发目标在三维空间内的位置将根据被测试对象的定位信息进行重定位，即根据被测试对象的定位点坐标，获得相对于理想的期望三维激发目标的定位点坐标的三维偏移值，并将此偏移值计入理想的期望三维激发目标的坐标系统，从而得到修正后的期望三维激发目标的三维空间位置信息。这里的B1+映射采用的是实翻转角成像法(Actual Flip-Angle Imaging)来实现的三维映射。实翻转角成像法在V. Yarnykh发表的论文中有具体描述(“Actual flip-angle imaging in the pulsed steady state :A method for rapid three-dimensional mapping of the transmitted radiofrequency field， “ Magnetic Resonance in Medicine，vol. 57，pp.192-200，2007)。然后，将根据所述的期望实现的三维激发目标和射频线圈各发射信道的三维空间敏感度分布情况，得到相应的三维激发K空间内的射频能量权重分布，将其中的高权重分布范围的主体部分作为用以“装载”和“约束”之后设计的K空间轨迹的容器范围、即三维激发K空间轨迹容器的范围。具体的方法是，首先在多发射信道并行激发条件下，根据小翻转角前提下的选择性射频脉冲设计理论，可有Mxy(T) = I7M,YjlS1(T)I Wl^dk根据空间分辨率ΔΓ与K空间最大值kmax之间的对应关系kmax= l/2/Δι·，以及视场FOV与K空间分辨率Δ k之间的对应关系Ak= 1/2/F0V,将空间和K空间离散化为rm， m = 1，· · ·，Nk 和 kn，η = 1，· · ·，Νκ，，可有M = Σ LS1 · A · W1其中M是表征空间范围FOX内的mxy(r)分布的Nk元向量，即表征激发目标的第一向量。&是表征第1个激励源产生的射频场空间敏感度的K元向量对角化后的KXNk矩阵，A表达实际空间与K空间之间关系的NkXNk的矩阵且满足= iYM0eir"-k"。将上式重新组织如下
权利要求
1.一种用于磁共振成像的三维空间选择性激励的序列设计方法，其特征是，包括如下步骤(1)确定磁共振系统的成像空间内的期望三维激发目标，并且通过B1+映射方法计算出所述磁共振系统的射频线圈的各发射信道所对应的三维空间敏感度分布情况；(2)根据所述期望三维激发目标和所述三维空间敏感度分布情况，计算出相应的三维激发K空间轨迹容器的范围；(3)根据所述期望三维激发目标和所述三维激发K空间轨迹容器的范围，确定相应的三维激发K空间轨迹；(4)根据所述三维激发K空间轨迹和所述磁共振系统中的梯度驱动单元的梯度脉冲的最大幅值和最大切换率，计算出所述梯度驱动单元的梯度脉冲波形；(5)根据所述期望三维激发目标、所述三维空间敏感度分布情况和所述梯度脉冲波形， 计算出所述射频线圈的各发射信道的射频脉冲的包络波形；由所述梯度脉冲波形、所述射频脉冲的包络波形和所述射频脉冲的载波频率构成所述三维空间选择性激励的序列。
2.根据权利要求1所述的用于磁共振成像的三维空间选择性激励的序列设计方法，其特征是对步骤(5)所得到的序列采用变速选择性激发方法或可变转换率螺旋型轨迹设计方法进行优化。
3.根据权利要求1所述的用于磁共振成像的三维空间选择性激励的序列设计方法，其特征是，所述步骤(2)中，按以下步骤计算出相应的三维激发K空间轨迹容器的范围1)根据所述期望三维激发目标，生成用于表征所述期望三维激发目标的第一向量；根据所述射频线圈的各发射信道的三维空间敏感度分布情况，生成满足以下关系式(I)的第一矩阵Afull = [ S1A,…，S1A ](I)式(I)中，^^即所述第一矩阵，&是表示所述射频线圈的第1个发射信道的三维空间敏感度分布情况的对角化矩阵，A是表达实际空间与K空间之间关系的矩阵且满足以下关系式(II)A(m,M) = iyeir"'k(II)式(II)中dfe/ )表达矩阵J的第《行、第列的元素的数值，i是虚数符号，r是磁旋比，4是对实际空间离散化得到的第 个点的空间坐标，先是对K空间离散化得到的第/7 个点的K空间坐标；2)根据所述第一向量和所述第一矩阵，计算得到用于表征所述射频线圈的所有发射信道所对应的激发K空间的能量分布的第二向量；3)将所述第二向量根据以下关系式(III)，分离为用于表达所述射频线圈的各发射信道所对应的激发K空间的能量分布的向量组Wfull = [ Wp …，WlY(III)式(ΙII)中，即所述第二向量，W1是表达所述射频线圈的第/个发射信道所对应的激发κ空间的能量分布，上标τ是矩阵转置符号；4)将所述向量组的幅值加和，并将该加和结果归一化为表征了包含所述射频线圈所有发射信道所对应的激发K空间的能量分布情况的能量权重分布；5)根据所述能量权重分布，确定出所述K空间轨迹所需遍历的高权重区域；6)根据所述高权重区域内的主体部分，确定一个规则封闭的形状并将该形状作为K空间轨迹容器的参考范围；7)根据所述三维空间敏感度分布情况，计算出所述射频线圈的各发射信道所对应的三维激发K空间内的点扩散作用范围；根据所述点扩散作用范围削减所述K空间轨迹容器的参考范围，得到K空间轨迹容器的范围。
4.根据权利要求1所述的用于磁共振成像的三维空间选择性激励的序列设计方法，其特征是，所述步骤(3)中，按以下步骤确定相应的三维激发K空间轨迹a)根据所述期望三维激发目标，确定三维K空间轨迹的类型；b)根据所述期望三维激发目标、所述三维K空间轨迹的类型，确定三维K空间轨迹的径向采样间隔、纵向采样间隔；c)根据所述三维K空间轨迹的类型、径向采样间隔和纵向采样间隔，以及所述K空间轨迹容器的范围，确定相应的三维激发K空间轨迹。
5.根据权利要求4所述的用于磁共振成像的三维空间选择性激励的序列设计方法，其特征是，所述步骤a)中，所述三维K空间轨迹的类型为立体回波型或堆栈螺旋型。
全文摘要
本发明公开一种用于磁共振成像的三维空间选择性激励的序列设计方法。本发明的序列设计方法，能够根据激发目标和射频线圈多个发射信道的空间敏感度情况优化的确定合适的激发K空间轨迹，从而确定梯度脉冲波形和各发射信道所对应的射频脉冲包络波形。磁共振系统中的梯度驱动单元和射频驱动单元据此生成梯度脉冲和射频脉冲，并驱动梯度线圈和射频线圈将其同时施加于扫描空间内，从而实现期望的三维空间内选择性激发的目标。
文档编号G01R33/32GK102283649SQ201110121218
公开日2011年12月21日 申请日期2011年5月11日 优先权日2011年5月11日
发明者刘锋, 夏灵, 邵汀汀 申请人:浙江大学