专利名称：用于磁共振系统的激发脉冲序列生成器的制作方法
技术领域：
本发明涉及用于磁共振系统的激发脉冲序列生成器。
背景技术：
当人体组织等物质置于均勻静磁场(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)，是以厚平板内的单次均勻性激发为前提，利用三维傅里叶成像法进行磁共振成像，该方法能够在提高成像分辨率的同时保证分辨率的各项同性，大大提高了成像信噪比。
发明内容
本发明所要解决的问题是提供一种用于磁共振系统的激发脉冲序列生成器，它可生成实现磁共振系统的泛三维体选择性激励的序列。本发明解决其技术问题所采取的技术手段是该用于磁共振系统的激发脉冲序列生成器包括输入装置，所述输入装置用于接收来自所述磁共振系统的计算机的关于所述磁共振系统的成像空间内的期望三维激发目标数据、关于所述磁共振系统的激发脉冲序列配置方式的指令、关于所述磁共振系统的梯度驱动单元的梯度脉冲的最大幅值和最大切换率、 以及关于所述磁共振系统的射频驱动单元的射频脉冲的载波频率，并接收来自所述磁共振系统的B1+映射处理器的关于所述磁共振系统的射频线圈的各发射信道的三维空间敏感度分布数据，接收来自所述磁共振系统的扫描室接口单元的关于所述磁共振系统的被测试对象的定位信息数据，以及接收来自所述磁共振系统的脉冲序列时序控制器的关于所述磁共振系统的脉冲序列时序指令；第一处理器，所述第一处理器用于根据所述被测试对象的定位信息数据对所述期望三维激发目标数据进行修正，得到修正后的期望三维激发目标数据，并根据所述激发脉冲序列配置方式的指令，将所述修正后的期望三维激发目标数据发送给第二处理器或默认序列配置器；第二处理器，所述第二处理器用于根据所述三维空间敏感度分布数据和所述修正后的期望三维激发目标数据，得到相应的三维激发K空间轨迹容器的范围；第三处理器，所述第三处理器用于根据所述修正后的期望三维激发目标数据和所述三维激发K空间轨迹容器的范围的数据，确定相应的三维激发K空间轨迹；第四处理器，所述第四处理器用于根据所述三维激发K空间轨迹的数据和所述梯度脉冲的最大幅值和最大切换率，得到所述梯度驱动单元的第一梯度脉冲波形；第五处理器，所述第五处理器用于根据所述三维空间敏感度分布数据、所述修正后的期望三维激发目标数据和所述梯度脉冲波形的数据，得到所述射频驱动单元的关于所述各发射信道的射频脉冲的第一包络波形；默认序列配置器，所述默认序列配置器用于根据所述三维空间敏感度分布数据、 所述修正后的期望三维激发目标数据，通过已知的序列配置方式得到所述梯度驱动单元的第二梯度脉冲波形和所述射频驱动单元的关于各发射信道的射频脉冲的第二包络波形；输出装置，所述输出装置用于根据所述激发脉冲序列配置方式的指令和所述脉冲序列时序指令，将所述第一梯度脉冲波形或所述第二梯度脉冲波形发送给所述磁共振系统的梯度驱动单元、将所述射频脉冲的第一包络波形或所述射频脉冲的第二包络波形发送给所述磁共振系统的射频驱动单元，以及将所述射频脉冲的载波频率发送给所述磁共振系统的射频驱动单元。与现有技术相比，本发明的有益效果是(1)本发明通过对传统的激发脉冲序列生成器进行改进，实现了基于并行激发技术的泛三维体选择性激发的序列的生成。(2)第二处理器可用于确定K空间轨迹容器的范围，该范围能够直观地限定K空间轨迹的分布，有助于实现泛三维体选择性激发的序列的生成；不仅如此，对K空间轨迹的分布的限定可提高K 空间轨迹的采样效率，从而减少激发序列的时长、减小射频脉冲幅值。( 第三处理器根据第二处理器得到的K空间轨迹容器的范围，能够采用任意的K空间轨迹类型，从而得到K空间轨迹；第三处理器优先采用堆栈螺旋型轨迹来设计序列，能够大大减小激发所需要的射频脉冲包络波形的最大幅值。(4)第五处理器采用变速选择性激发方法或可变转换率螺旋型轨迹设计方法对梯度脉冲波形和计算出来的射频脉冲包络波形所组成的序列进行优化， 将进一步减小激发所需要的射频脉冲包络波形的最大幅值。减小激发所需要的射频脉冲包络波形的最大幅值，即减小射频能量，不仅能够降低SAR值，提高磁共振扫描的安全性，还有助于在射频能量强烈衰减的低信号区域内实现接近期望的激发目标。以上都是现有技术所不具备的特征。
图1是磁共振系统的结构示意图；图2是磁共振系统中的脉冲序列生成器的工作连接关系示意图；图3是激发脉冲序列生成器的工作连接关系示意图；图4是激发脉冲序列生成器的工作流程图；图5 (a)为传统的三维K空间轨迹中的回波立体型轨迹；图5(b)为本发明的K空间轨迹容器；图5 (c)为利用图5 (b)的K空间轨迹容器对图5 (a)的回波立体型轨迹进行优化后得到的轨迹；图6是图5 (c)的A-A平面映射图；图7是图5 (c)的B-B平面映射图；图8 (a)为传统的三维K空间轨迹中的堆栈螺旋型轨迹，图8 (b)为本发明的K空间轨迹容器，图8 (c)为利用图8 (b)的K空间轨迹容器对图8 (a)的堆栈螺旋型轨迹进行优化后得到的轨迹；图9是图8 (c)的A-A平面映射图；图10是图8(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 ；梯度场同时作用来完成对该磁共振信号的空间编码。用以生成该特定梯度场的梯度脉冲构成一个成像脉冲序列。一个激发过程和一个成像过程构成一个完整的单次扫描过程。而一个激发脉冲序列和一个成像脉冲序列构成一个完整的单次磁共振成像脉冲序列。参见图2，脉冲序列生成器116 —般主要包括激发脉冲序列生成器200、脉冲序列时序控制器202和成像脉冲序列生成器204。激发脉冲序列生成器200生成用以实现一个激发过程的脉冲序列，成像脉冲序列生成器204生成用以实现一个成像过程的脉冲序列， 脉冲序列时序控制器202通过发送指令给激发脉冲序列生成器200和成像脉冲序列生成器 204，来控制激发脉冲序列与成像脉冲序列的执行时间和顺序。激发脉冲序列生成器200与B1+映射处理器114相连。B1+映射处理器114用于计算射频线圈(体线圈108或阵列线圈110)各发射信道对应的多维空间敏感度分布情况。 激发脉冲序列生成器200将根据该多维空间敏感度分布情况的数据、以及来自计算机100 的关于磁共振系统的成像空间内的多维空间激发目标的数据，计算出能够实现该激发目标的梯度脉冲波形和射频线圈各发射信道对应的射频脉冲波形，并将波形数据分别发送给梯度驱动单元122和射频驱动单元124。射频驱动单元IM包括一个或多个射频脉冲发射器。在激发过程中，射频脉冲发射器依据来自于激发脉冲序列生成器200的数据和指令产生既定载波频率、既定包络波形的射频脉冲，并驱动射频线圈(如体线圈108或阵列线圈110)产生射频�。煌碧荻惹ピ谰堇醋杂诩し⒙龀逍蛄猩善�200的数据和指令产生既定的梯度波形，并驱动梯度线圈106产生梯度�。愿蒙淦党〗锌占浔嗦耄痪榷ǖ闹葱惺奔洌迪忠桓黾し⒐�。射频驱动单元124同时包含一个或多个射频信号接收通道，每个接收通道包含一个能够放大经由射频线圈108或110接收到的磁共振信号的射频信号放大器。在成像过程中，被测试对象(如图1中的人体)中的质子受到激励后衰减产生磁共振信号。同时梯度驱动单元依据来自于成像脉冲序列生成器204的数据和指令产生既定的梯度波形，并驱动梯度线圈106产生梯度�。愿么殴舱裥藕沤锌占浔嗦�。射频线圈(如体线圈108或阵列线圈110)接收该信号，并传递给射频驱动单元124，射频驱动单元IM依据来自于成像脉冲序列生成器204的数据和指令该对信号进行放大、解调、滤波及数字化处理。此外，在磁共振扫描过程中，脉冲序列时序控制器202同时通过生理信息采集控制器1 来实时的接收和分析被测试对象的相关生理信息。生理信息采集控制器1 通过一系列的传感器实时采集被测试对象心电信号、呼吸信号等生理信息。而脉冲序列时序控制器202通过分析这些信号，能够通过发送指令给激发脉冲序列生成器200和成像脉冲序列生成器204，实时的调整激发脉冲序列和成像脉冲序列，从而实现磁共振扫描与被测试对象的心跳和呼吸等生理运动的同步化，将扫描过程中伴随心跳和呼吸等生理运动的空间位移对磁共振成像结果的影响减至最小。激发脉冲序列生成器200还实时接收来自扫描室接口单元128的信息。扫描室接口单元1 与扫描室内用以探测被测试对象位置和磁体系统状况的传感器相连，并提供可操作平台，显示被测试对象位置和磁体系统运行状况。与扫描室接口单元1 相连的被测试对象定位系统130接收来自扫描室接口单元128的操作平台的数据，并依此控制扫描床 112的移动，将被测试对象平移到指定位置。扫描室接口单元1 将定位后的被测试对象的位置信息传送给激发脉冲序列生成器200。在成像后续过程中，经由射频驱动单元IM数字化后的磁共振信号被传递到数据采集单元118。数据采集单元118在计算机100的指令下配合成像脉冲序列生成器204传递的时序信息实时的接收磁共振数据。数据采集单元118提供足够的缓存来存储至少一个扫描周期内的数据。此后，数据采集单元118将磁共振数据传递给数据处理单元120用以进一步的处理和分析。数据处理单元120接收来自数据采集单元118的磁共振数据，并按照计算机100 的指令进行处理。该处理包括对原始K空间数据的傅里叶变换，对重建所得图像的滤波操作，用于功能磁共振成像的重计算，以及用于运动目标或流体成像的重计算等等。经由数据处理单元120得到的重建图像被传递给计算机100并存储其中。成像结果同时被传送到操作台显示器102或外置显示器132。计算机100也可以通过网络将图像传送到网络上的其他设备。更具体地，在激发过程中，激发脉冲序列生成器200依据计算机100的指令发送脉冲数据给梯度驱动单元122和射频驱动单元124，梯度驱动单元122据此生成梯度脉冲，并驱动梯度线圈106产生既定的梯度场。与此同时，射频驱动单元IM据此生成具有既定载波频率和包络的射频脉冲，并驱动体线圈108或阵列线圈110产生既定的射频场。通常该射频脉冲具有特定的载波频率和包络波形，从而配合梯度脉冲来实现对既定空间位置、既定空间范围内的质子的既定翻转角幅度的激发。在基于多发射信道射频线圈(如阵列线圈 110)的并行激发过程中，射频线圈各发射信道发送的射频脉冲具有同样的或相近的载波频率和时长，但是可以有不同的包络波形，结合各信道所对应的不同的射频场分布，最终产生期望中的激发结果。此结果包括被扫描体内受激发质子的空间范围、受激翻转角大小及其均勻性等等。因此，并行激发技术的实现，需首先获得射频线圈(如阵列线圈110)各发射信道所对应的射频场分布(或称各发射信道的空间敏感度)，并据此分布情况以及激发目标做脉冲序列设计。在本发明所采用的磁共振系统中，射频线圈(如阵列线圈110)各发射信道的空间敏感度的计算主要由B1+映射处理器114完成。进入B1+映射运行模式时，将首先通过多个特定的激发脉冲序列，得到多个与各发射信道空间敏感度相关的成像结果。更具体的，计算机100指令激发脉冲序列生成器200生成特定的激发脉冲序列波形，通过射频驱动单元IM和梯度驱动单元122驱动射频线圈(如阵列线圈110)的某一发射信道和梯度线圈106实现特定的激发；指令成像脉冲序列生成器204生成成像脉冲序列，并通过梯度驱动单元122驱动梯度线圈106进行空间编码、通过射频驱动单元IM控制具有均勻空间敏感度的接收通道的射频线圈(如体线圈108)接收经过空间编码的磁共振信号；指令数据采集单元118和数据处理单元120接收和处理数据从而得到成像结果，传递回计算机100并存储于计算机100 ；经过如上所述的多个特定的扫描过程，将该多个特定的激发脉冲序列激发得到的多个成像结果依次存储于计算机100中，并由计算机100传送给B1+映射处理器 114。B1+映射处理器114对这些成像结果做进一步的处理并计算出射频线圈各发射信道的空间敏感度情况。B1+映射处理器114计算完成后，将发送结束标记给计算机100，计算机 100据此退出B1+映射运行模式，并指令B1+映射处理器114将射频线圈所有发射信道的空间敏感度数据传送到激发脉冲序列生成器200。激发脉冲序列生成器200根据来自计算机 100的关于成像空间范围内的激发目标、关于梯度驱动单元122可生成的梯度脉冲的最大幅值和最大切换率、及关于射频驱动单元1 将生成的射频脉冲的载波频率，来自B1+映射处理器114的射频线圈各发射信道的空间敏感度数据，以及来自于扫描室接口单元128的关于被测试对象的定位信息数据，计算出梯度脉冲波形和射频线圈各发射信道对应的射频脉冲包络波形，并依据来自脉冲序列时序控制器202的脉冲序列时序指令，将计算结果及射频脉冲的既定载波频率等脉冲序列数据传送给梯度驱动单元122和射频驱动单元124 ； 由梯度驱动单元122和射频驱动单元IM分别生成梯度脉冲与射频脉冲，驱动梯度线圈106 和射频线圈(如阵列线圈110)来完成一个基于并行激发技术的激发过程。本发明改进了传统的激发脉冲序列生成器，能够实现基于并行激发技术的泛三维体选择性激发的序列的生成。在本发明中，如图3所示，激发脉冲序列生成器包括有输入装置300，第一处理器302、第二处理器304、第三处理器306、第四处理器308、第五处理器 310、默认序列配置器314和输出装置312。输入装置300用以接收来自于计算机100的关于成像空间内的期望三维激发目标数据、关于激发脉冲序列配置方式的指令、关于梯度驱动单元122可生成的梯度脉冲的最大幅值和最大切换率、关于射频驱动单元1 将生成的射频脉冲的载波频率，并接收来自于B1+映射处理器的关于射频线圈各发射信道的三维空间敏感度分布数据，接收来自于扫描室接口单元128的关于被测试对象的定位信息数据，以及接收来自于脉冲序列时序控制器202的脉冲序列时序指令。第一处理器302用以根据来自于输入装置300的被测试对象的定位信息数据对期望三维激发目标进行重定位，即根据被测试对象的定位点坐标，获得相对于期望三维激发目标的定位点坐标的三维偏移值，并将此偏移值计入期望三维激发目标，从而得到修正后的期望三维激发目标数据；然后根据来自于输入装置300的激发脉冲序列配置方式的指令，将修正后的期望三维激发目标数据发送给第二处理器304或默认序列配置器314 当激发脉冲序列配置方式的指令为非默认序列配置方式指令时，第一处理器302将修正后的期望三维激发目标数据发送给第二处理器304，当激发脉冲序列配置方式的指令为默认序列配置方式指令时，第一处理器302将修正后的期望三维激发目标数据发送给默认序列配置器 314。第二处理器304用以根据来自于输入装置300的关于射频线圈各发射信道三维空间敏感度分布数据和来自于第一处理器302的修正后的期望三维激发目标数据，得到出相应的三维激发K空间轨迹容器的范围。第三处理器306用以根据来自于第一处理器302的期望三维激发目标和来自于第二处理器304的三维激发K空间轨迹容器的范围的数据，确定相应的三维激发K空间轨迹。第四处理器308用以根据来自于第三处理器306的三维激发K空间轨迹的数据和来自输入装置300的梯度脉冲的最大幅值和最大切换率，得到梯度驱动单元122可生成的第一梯度脉冲波形。第五处理器310用以根据来自输入装置300的三维空间敏感度分布数据、来自于第一处理器302的修正后的期望三维激发目标数据和来自第四处理器308的梯度脉冲波形，得到对应于射频驱动单元124的关于射频线圈各发射信道的射频脉冲的第一包络波形。默认序列配置器314用于根据来自于输入装置300的关于射频线圈各发射信道的三维空间敏感度分布数据、来自于第一处理器302的修正后的期望三维激发目标数据，通过已知的序列配置方式确定出梯度驱动单元122的第二梯度脉冲波形和射频驱动单元IM 的关于射频线圈各发射信道的射频脉冲的第二包络波形。输出装置312用以根据来自于输入装置300的激发脉冲序列配置方式的指令和脉冲序列时序指令，将来自于第四处理器308的第一梯度脉冲波形数据或默认序列配置器 314的第二梯度脉冲波形数据发送给梯度驱动单元122、将来自于第五处理器310的射频脉冲的第一包络波形数据或默认序列配置器314的射频脉冲的第二包络波形数据和来自输入装置300的关于射频脉冲的载波频率数据发送给射频驱动单元124。当激发脉冲序列配置方式的指令为非默认序列配置方式指令时，输出装置312将来自于第四处理器308的第一梯度脉冲波形数据发送给梯度驱动单元122、将来自于第五处理器310的射频脉冲的第一包络波形数据和来自输入装置300的关于射频脉冲的载波频率数据发送给射频驱动单元124 ；当激发脉冲序列配置方式的指令为默认序列配置方式指令时，输出装置312将来自于默认序列配置器314的第二梯度脉冲波形数据发送给梯度驱动单元122、将来自于默认序列配置器314的射频脉冲的第二包络波形数据和来自输入装置300的关于射频脉冲的载波频率数据发送给射频驱动单元124。本发明激发脉冲序列生成器在非默认序列配置方式下的工作过程如下所述首先确定磁共振系统的成像空间内的期望三维激发目标，其参数包括激发目标在三维空间内的位置、形状、翻转角大小和分辨率要求，其中激发目标在三维空间内的位置将根据被测试对象的定位信息进行重定位；通过映射方法得到磁共振系统的射频线圈的各发射信道所对应的三维空间敏感度分布情况；根据期望三维激发目标、射频线圈各发射信道的三维空间敏感度分布情况，得到相应的三维激发K空间内的射频能量权重分布，并在此基础上，确定三维激发K空间轨迹容器的范围；将该K空间轨迹容器的范围作为限定K 空间轨迹分布的范围，并根据期望三维激发目标，选择K空间轨迹的类型、确定径向和纵向采样间隔，从而确定相应的三维激发K空间轨迹；根据三维激发K空间轨迹、以及梯度驱动单元122所能生成的梯度脉冲的最大幅值和最大切换率约束，得到出第一梯度脉冲波形； 根据期望三维激发目标、射频线圈各发射信道的三维空间敏感度分布情况、以及第一梯度脉冲波形，计算出射频线圈各发射信道所对应的射频脉冲的第一包络波形；由第一梯度脉冲波形、射频脉冲的第一包络波形和载波频率构成了实施三维空间选择性激励的非默认配置序列。梯度驱动单元122和射频驱动单元IM根据非默认配置序列数据生成相应的梯度脉冲和射频脉冲，并驱动梯度线圈和射频线圈将其同时施加于扫描空间内，从而实现期望的三维空间内选择性激发的目标。此外，本发明激发脉冲序列生成器在非默认序列配置方式下的工作过程还可包括对非默认配置序列采用变速选择性激发方法或可变转换率螺旋型轨迹设计方法进行优化，以进一步减小实现激发所需要的射频脉冲的第一包络波形的最大幅值。本发明激发脉冲序列生成器在非默认序列配置方式下的工作过程将在下文中详细描述。首先进行如下定义磁旋比Y，三维空间坐标向量r = [rx, ry, rJT,沿x、y、z各方向的梯度脉冲波形组合G(t) = [Gx (t), Gy (t), Gz(t) ]τ0三维K 空间轨迹 k (t) = [kx (t)，ky (t)，kz (t) ]τ。梯度脉冲波形与K空间轨迹的关系为众⑴=-γ fj G(t)出，这里T为脉冲时长。各发射信道的三维空间敏感度分布情况，1 = 1，2，. . .，L，L为信道总个数。各发射信道对应的射频脉冲包络波形B1,Jt) ,1 = 1,2,...,L0各发射信道对应的三维激发K空间的射频能量分布W1GOd = 1,2,...,L0横向磁化矢量Mxy = Mx+iMy,这里Mx为χ方向磁化矢量、My为y方向磁化矢量。纵向磁化矢量Mz，令质子受激发前或稳态时其值为M。。参见图4，在本发明激发脉冲序列生成器在非默认序列配置方式下的工作过程中，首先需要确定期望实现的三维激发目标、以及通过映射，得到关于射频线圈各发射信道的三维空间敏感度分布情况。这里的三维激发目标具体包括，该激发目标在三维空间内的位置、形状、分布范围，该激发目标内部的翻转角大小、均勻性，以及该激发目标的分辨率要求；其中激发目标在三维空间内的位置将根据被测试对象的定位信息进行重定位，即根据被测试对象的定位点坐标，获得相对于理想的期望三维激发目标的定位点坐标的三维偏移值，并将此偏移值计入理想的期望三维激发目标的坐标系统，从而得到修正后的期望三维激发目标的三维空间位置信息。这里的映射采用的是实翻转角成像法 (Actual Flip-Angle Imaging)来实现B1+的三维映射。实翻转角成像法在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 (r) = IYM0 Σ LSX (r) f W1 (k) eir ·kdk根据空间分辨率ΔΓ与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)分布的队元向量。&是表征第1个激励源产生的射频场空间敏感度的Nk元向量对角化后的NkXNk矩阵，A表达实际空间与K空间之间关系WNkXNk的矩阵且满足=。将上式重新组织如下
权利要求
1. 一种用于磁共振系统的激发脉冲序列生成器，其特征是，包括 输入装置，所述输入装置用于接收来自所述磁共振系统的计算机的关于所述磁共振系统的成像空间内的期望三维激发目标数据、关于所述磁共振系统的激发脉冲序列配置方式的指令、关于所述磁共振系统的梯度驱动单元的梯度脉冲的最大幅值和最大切换率、以及关于所述磁共振系统的射频驱动单元的射频脉冲的载波频率，并接收来自所述磁共振系统的B1+映射处理器的关于所述磁共振系统的射频线圈的各发射信道的三维空间敏感度分布数据，接收来自所述磁共振系统的扫描室接口单元的关于所述磁共振系统的被测试对象的定位信息数据，以及接收来自所述磁共振系统的脉冲序列时序控制器的关于所述磁共振系统的脉冲序列时序指令；第一处理器，所述第一处理器用于根据所述被测试对象的定位信息数据对所述期望三维激发目标数据进行修正，得到修正后的期望三维激发目标数据，并根据所述激发脉冲序列配置方式的指令，将所述修正后的期望三维激发目标数据发送给第二处理器或默认序列配置器；第二处理器，所述第二处理器用于根据所述三维空间敏感度分布数据和所述修正后的期望三维激发目标数据，得到相应的三维激发K空间轨迹容器的范围；第三处理器，所述第三处理器用于根据所述修正后的期望三维激发目标数据和所述三维激发K空间轨迹容器的范围的数据，确定相应的三维激发K空间轨迹；第四处理器，所述第四处理器用于根据所述三维激发K空间轨迹的数据和所述梯度脉冲的最大幅值和最大切换率，得到所述梯度驱动单元的第一梯度脉冲波形；第五处理器，所述第五处理器用于根据所述三维空间敏感度分布数据、所述修正后的期望三维激发目标数据和所述梯度脉冲波形的数据，得到所述射频驱动单元的关于所述各发射信道的射频脉冲的第一包络波形；默认序列配置器，所述默认序列配置器用于根据所述三维空间敏感度分布数据、所述修正后的期望三维激发目标数据，通过已知的序列配置方式得到所述梯度驱动单元的第二梯度脉冲波形和所述射频驱动单元的关于各发射信道的射频脉冲的第二包络波形；输出装置，所述输出装置用于根据所述激发脉冲序列配置方式的指令和所述脉冲序列时序指令，将所述第一梯度脉冲波形或所述第二梯度脉冲波形发送给所述磁共振系统的梯度驱动单元、将所述射频脉冲的第一包络波形或所述射频脉冲的第二包络波形发送给所述磁共振系统的射频驱动单元，以及将所述射频脉冲的载波频率发送给所述磁共振系统的射频驱动单元。
全文摘要
本发明提供一种用于磁共振系统的激发脉冲序列生成器，其生成的激发脉冲序列应用于磁共振系统能够实现优化的基于并行激发技术的三维空间选择性激励。本发明的激发脉冲序列生成器能够根据激发目标和射频线圈多个发射信道的空间敏感度情况优化的确定合适的激发K空间轨迹，从而确定梯度脉冲波形和各发射信道所对应的射频脉冲包络波形。磁共振系统的梯度驱动单元和射频驱动单元据此生成梯度脉冲和射频脉冲，并驱动梯度线圈和射频线圈将其同时施加于扫描空间内，从而实现期望的三维空间内选择性激发的目标。
文档编号G01R33/36GK102288929SQ201110121258
公开日2011年12月21日 申请日期2011年5月11日 优先权日2011年5月11日
发明者刘锋, 夏灵, 邵汀汀 申请人:浙江大学
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
