使用绝热射频脉冲测量射频B1场分布的磁共振成像方法与流程

本发明属于核磁共振成像技术领域，尤其是涉及一种使用绝热射频脉冲测量射频b1场分布的磁共振成像方法。该发明适用于测量磁共振成像中表面线圈、体线圈或相位阵列线圈发射时射频b1场在不同组织或样品中的分布，用于图像重建以及b1场校正等。

背景技术：

磁共振成像是诊断与评估疾病进展的重要临床方法，具有无电离辐射、非侵入性、高空间分辨率、任意层面成像、组织对比度高等优点，在肿瘤、脏器以及软组织病变等诊断方面具有不可替代的作用。

b1场表示磁共振线圈发射时射频脉冲的强度。磁共振成像中，因为组织不规则结构以及线圈的发射场不均匀，导致射频b1场在空间中呈现不均匀分布。随着目前磁共振谱仪的磁场强度越来越高，射频b1场不均匀问题越发严重。在目前磁共振中常用的并行成像以及化学交换饱和转移(chemicalexchangesaturationtransfer,cest)成像等方法中，也需要射频b1场分布进行图像重建和校正。发展测量射频b1场分布的磁共振方法，可以解决以上难题。

b0场表示磁共振成像时空间中某一点感应到的磁场大小，与核自旋频率紧密相关。用b0场偏移频率来表示射频激发时射频频率与核自旋频率的差值。b0场偏移会对b1场分布的测量产生影响，主要体现在b0场偏移较大时，偏共振效应下磁共振信号与b0和b1均有关。因此，测量b1场分布时不能忽略b0场偏移的影响。

现有的测量b1场分布的磁共振方法，主要是基于信号幅度的双翻转角方法，例如，inskoek等人[mappingoftheradiofrequencyfield.j.magn.reson.a.1993.45:p.82-85.]提出的基于不同翻转角的b1场测量方法，以及基于该方法的改进，例如，cunninghamch等人[sturateddouble-anglemethodforrapidb1+mapping.magn.reson.med.2006.55:p.1326-1333.]中提出的速度更快的基于不同翻转角b1场测量方法。然而，翻转角不同时脉冲的激发波形不是线性的，而且该方法对b0场偏移比较敏感。

yarnykhvl等人[actualflip-angleimaginginthepulsedsteadystate:amethodforrapidthree-dimensionalmappingofthetransmittedradiofrequencyfield.magn.reson.med.2007.57:p.192-200.]提出的afi(actualflip-angleimaging)方法，利用相同激发角下，两个不同恢复时间tr1和tr2的信号强度的比值来计算b1场分布，该方法速度快，但是，对b0场偏移以及运动比较敏感。

基于信号相位的方法，主要有sacolickli等人[b1mappingbybloch-siegertshift.magn.reson.med.2010.63:p.1315-1322.]提出的测量方法，bloch-siegertshift的相位变化与b1场平方成正比。相位方法中相位重建比较复杂，而且，可能伴随相位卷折的问题(相位超过360°)。

此外，其他方法包括patrickschuenke等人[simultaneousmappingofwatershiftandb1(wasabi)–applicationtofield-inhomogeneitycorrectionofcestmridata.magn.reson.med.2017.77:p.571-580.]提出的利用偏共振效应同时拟合b0和b1场分布的方法。该方法应用于化学交换饱和转移成像中，需要采集10幅或更多图像进行曲线拟合。

本发明提出一种使用绝热射频脉冲测量射频b1场分布的磁共振成像方法，是一种对b0场偏移不敏感、可广泛应用于不同脉冲序列中的b1测量方法，为射频b1场的测量以及图像重建校正等提供技术支持。

技术实现要素：

为解决现有技术存在的上述问题，本发明提供了一种使用绝热射频脉冲测量射频b1场分布的磁共振成像方法，本发明基于以下思路：

在旋转坐标系下，可以认为磁化矢量沿着有效b1场即beff场进动。其中，beff＝[(γb1)²+(△ω)²]^1/2，γ代表原子核的旋磁比，△ω表示b0场偏移的大小，beff表征beff场的强度，b1表征b1场的强度。进动频率ω＝γbeff。beff场与z轴(即主磁场方向)的夹角θ＝arctan(γb1/|△ω|)。

绝热脉冲是一种特殊条件下的射频脉冲。在绝热条件下，与初始beff场平行的磁化矢量将被绝热脉冲自旋锁定，并始终与beff场保持在同一直线上。

绝热脉冲需要满足绝热条件，定义则要求k>>1来保证绝热性，即绕beff场进动频率要远大于beff场与z轴夹角的变化。为了满足绝热条件，为射频脉冲设置一个大的频率偏置，定义为△ωrf,保证△ωrf>△ω。偏置值的加入既增加了beff值，又减小了θ随b1场的变化，保证了射频脉冲的绝热性。

选择一个初始强度为0，结束时强度达到最大的射频脉冲，设该射频脉冲的最大强度为b1max，在绝热条件下，即设置射频脉冲的频率偏置为△ωrf,保证△ωrf>△ω，射频脉冲开始时磁化矢量与z轴的夹角为0，射频脉冲结束时磁化矢量与z轴的夹角为θ1＝arctan[γb1max/|(△ω+△ωrf)|]。

设初始磁化矢量为m0，绝热脉冲结束后，磁化矢量与z轴方向的夹角为θ1＝arctan[γb1max/|(△ω+△ωrf)|]，施加损毁梯度场损毁掉横向磁化矢量，只保留初始磁化矢量m0在z轴方向上的投影mz1，则mz1＝m0cos(θ1)。同理，改变频率偏置为-△ωrf，则磁化矢量与z轴夹角θ2＝arctan[γb1max/|(△ω-△ωrf)|]，mz2＝m0cos(θ2)。

设b1,obs为计算的b1场的值，其计算公示如下：

公式2中，△ω²/(△ωrf²-△ω²)项为观测的b1,obs与实际b1max之间的偏差。在选定△ωrf为4000hz的情况下，△ω在-400到400hz范围内偏差小于等于1.01％。为保证脉冲的绝热性，以及测量结果对△ω不敏感，△ωrf的选值越大越好，2000hz或以上能满足实验条件。

使用绝热射频脉冲测量射频b1场分布的磁共振成像方法，包括以下步骤：

步骤1，设计绝热脉冲，绝热脉冲的起始强度为零，绝热脉冲逐渐增加到最大强度b1max。设置频率偏置△ωrf以及绝热脉冲时间；

步骤2，计算绝热脉冲的k值，保证k值满足绝热条件。k值的计算由公式给出，k>100情况下满足绝热条件。最大强度b1max越大，k值越小；△ωrf越大，k值越大；在满足k值的情况下脉冲时间尽量短，以减少驰豫带来的影响，其中，beff表征beff场的强度，γ为原子核的旋磁比，θ为beff场与z轴(即主磁场方向)的夹角，θ＝arctan(γb1/|△ω|)，b1表征b1场的强度，△ω为b0场偏移的大小；

步骤3，不施加绝热脉冲或将绝热脉冲的最大强度b1max设为零，采样获得初始磁化矢量大小m0；

步骤4，施加频率偏置为△ωrf的绝热脉冲以及损毁梯度，采样获得磁化矢量大小mz1；

步骤5，施加频率偏置为-△ωrf的绝热脉冲以及损毁梯度，采样获得磁化矢量大小mz2；

步骤6，利用以下公式计算b1,obs，

其中，m0为初始磁化矢量。

将b1,obs除以最大强度b1max进行归一化，得到归一化的射频b1场分布；

步骤7，依据所选择的频率偏置△ωrf，根据△ω²/(△ωrf²-△ω²)来计算由于b0场偏移给b1测量带来的偏差，△ωrf>△ω。

本发明相对于现有技术，具有以下有益效果：

1、本方法中绝热条件下绝热脉冲强度与信号强度呈现简单的三角函数关系，便于计算与分析。

2、本方法由于b0场偏移△ω带来的计算偏差为△ω²/(△ωrf²-△ω²)，通过施加大△ωrf值可以很好的避免b0场偏移对b1测量的影响。

3、本方法中绝热脉冲和损毁梯度施加于采样之前，可以集成为一个子模块，不受采样方法的限制。

附图说明

图1是本发明方法的流程图；

图2(a)展示了一个简单的线性增加的脉冲形状，b1max为40μt，时长为2ms，

图2(b)是图2(a)所示脉冲在偏置△ωrf为4000hz时k值随脉冲时间的变化，整个过程中k值远大于1，满足绝热条件。本发明所使用的脉冲不局限于该脉冲形状；

图3为使用绝热脉冲测量b1场分布磁共振成像方法的脉冲序列图；本图中仅介绍了使用自旋回波进行数据采集的脉冲序列图，本发明所使用的数据采集方式不受限制；

图4是绝热条件下磁化矢量随脉冲强度的演化图，初始磁化矢量m0被脉冲翻转至beff方向，此时在z轴方向上投影为mz1,2；

图5为3％琼脂糖溶液样品使用本发明的方法测量的归一化b1场分布图；

图6为使用绝热脉冲测量b1场分布的磁共振成像方法时，由于b0场偏移△ω带来的测量误差。

具体实施方式

为了便于本领域普通技术人员理解和实施本发明，下面结合实施例对本发明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不用于限定本发明。

实施例：

本实施例中所用样品为3％琼脂糖溶液，置于外径10mm的核磁样品管。使用bruker400m宽腔小动物成像仪，内径10mm的马鞍状体线圈，温度控制在300k。该条件下，样品的t1驰豫时间约为3400ms，t2驰豫时间约为40ms，脉冲序列tr时间设为5s。本方法测试对象可以是样品，也可以是动物或人。射频线圈可以是表面线圈，也可以是体线圈或相位阵列线圈等。

使用绝热射频脉冲测量射频b1场分布的磁共振成像方法，包括以下步骤：

步骤1、如图1流程所示，使用绝热脉冲测量b1场分布的磁共振成像方法首先要确定绝热脉冲的波形，使其满足初始强度为0，绝热脉冲结束前达到最大强度b1max。本实施例中，使用简单的线性变化的脉冲形状，如图4所示。设定最大强度b1max为40μt，设置频率偏置△ωrf为4000hz，绝热脉冲时间2ms。本方法使用绝热脉冲的脉冲形状不局限于线性脉冲，也可以是各种满足从0到b1max递增的形状脉冲。

步骤2、计算绝热度，也就是k值。首先，本实施例中绝热脉冲由1000个点构成。由公式beff＝[(γb1)²+(△ωrf)²]^1/2分别计算每个点对应的beff值，然后使用公式θ＝arctan[γb1/|(△ωrf)|]分别计算每个时间点对应的角度θ。两相邻点θ的差值与相邻点时间差的比值即为dθ/dt，最后可以使用公式计算k值。本实施例中所用脉冲的k值随时间分布如图2(b)所示，满足k>>1的绝热条件。

本实施例所用脉冲序列如图3所示，施加完步骤2中的绝热脉冲后，初始磁化矢量将被翻转到最终beff方向，如图4所示。此时立即施加1ms的x/y/z三个方向上的损毁梯度来损毁掉横向磁化矢量，然后开始数据采集。本实施例中使用的数据采集方式是快速自旋回波方式，本发明不局限采样方式，在施加完绝热脉冲和损毁梯度后，使用梯度回波或平面回波等采样方式同样可行。

步骤3，不施加绝热脉冲或将绝热脉冲的最大强度b1max设为零，采样获得初始磁化矢量大小m0；

步骤4，施加频率偏置为△ωrf的绝热脉冲以及损毁梯度，采样获得磁化矢量大小mz1；

步骤5，施加频率偏置为-△ωrf的绝热脉冲以及损毁梯度，采样获得磁化矢量大小mz2；

步骤6，利用以下公式计算b1,obs，

其中，m0为初始磁化矢量。

将b1,obs除以最大强度b1max进行归一化，得到归一化的射频b1场分布；

先采集不施加绝热脉冲(或脉冲最大强度为0)时的第一图像，然后分别采集施加射频偏置为△ωrf和-△ωrf的绝热脉冲的两幅图像，分别为第二图像和第三图像。

对第一图像、第二图像、第三图像分别进行背景分割，使第一图像、第二图像、第三图像的图像背景为零。然后对第一图像、第二图像、第三图像进行点对点计算。取第一图像上任意一点，其信号强度设为m0，取分别施加射频偏置为△ωrf和-△ωrf的第二图像和第三图像上与第一图像所取点相同位置的点，信号强度分别设为mz1、mz2。那么，测量到的b1,obs可以通过公式(1)即可计算得出。

通过步骤6将图像上每一点对应的b1,obs计算出来，得到b1分布图。然后将b1分布图每个点除以系统设定的b1max值，得到归一化的b1分布图。本实施例中计算的归一化后b1分布图如图5所示，样品管靠近线圈区域有局部的b1场不均匀。

步骤7，依据所选择的频率偏置△ωrf，根据△ω²/(△ωrf²-△ω²)来计算由于b0场偏移给b1测量带来的偏差，△ωrf>△ω。

b0场偏移△ω带来的计算误差为△ω²/(△ωrf²-△ω²)，图6显示了不同△ωrf下计算误差与△ω的关系。本实施例使用的△ωrf为4000hz，在△ω为400hz的情况下计算误差约1％，计算结果在一定b0场偏移下仍有较好的准确度。

需要指出的是，本发明中所描述的具体实施例仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例作各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或超越所附权利要求书所定义的范围。