[发明专利]一种利用时空信息的磁共振动脉自旋标记序列部分容积校正方法

说明书

技术领域

本发明属于医学影像的部分容积效应校正方法技术领域，特别涉及一种利用时空信息的磁共振动脉自旋标记序列的部分容积校正方法。

背景技术

磁共振(Magnatic Resonance Imaging,MRI)灌注成像是用来反映组织的微血管分布及血流灌注情况的MRI检查技术。动脉自旋标记技术(Arterial Spin Labeling,ASL)，通过磁化血液中水分子作为内源性对比剂，来评价脑血流灌注情况，可完全无创检测，可重复性高，且可获得绝对定量的脑血流(Cerebral Blood Flow,CBF)，完全不受血脑屏障的影响。

ASL技术以磁标记的动脉血作为内源性对比剂，在成像平面上游，利用反脉冲标记动脉血的质子；延迟一段时间后，待标记的血液进入组织，血液与组织进行物质交换后成像(即标记图像，label)。label图像包括静态组织和标记动脉血的信号。为了消除静态组织的信号，进行另外一次未标记血成像(即控制图像，control)，control图像只包括静态组织信号。将label图像与control图像减影，所得的差值图像只与流入成像平面的标记动脉血有关。ASL图像的缺点在于它对噪声敏感。一般说来，label图像与control图像的差值大约是control像灰度值的1～2％。因此，一对label/control图像的差值不能够有效的反映灌注情况，需要大约60对数据来获得合理的信噪比(Signal to Noise Ratio,SNR)。

同时，作为一组配对图像，label和control图像需要快速、连续获取。所以，快速成像方法(如EPI方法)普遍应用于ASL成像，但这类方法的使用却降低了图像的空间分辨率，使图像受部分容积(Partial Volume,PV)效应影响严重。目前，针对ASL图像的PV校正算法很少，Asllani等(Asllani I,et al.Magn Reson Med,2008；60(6):1362-1371.)首先使用局部区域的线性回归(Linear Regression,LR)方法对图像进行校正，但该方法却使局部区域平滑，易引起CBF计算的误差。之后，Liang等(Liang X,et al.Magn Reson Med,2013；69(2):531-537.)采用最小截平方和方法对平滑效应进行后处理。同时，Chappell等(Chappell MA,et al.Magn Reson Med,2011；65(4):1173-1183.)对多个反转时间的ASL序列的PV进行校正。

虽然上述方法已可以在空间域对ASL数据进行一定程度的PV校正，但ASL序列中的时间信息并未被充分使用。

发明内容

针对目前ASL序列的PV校正方法中，时间信息利用率低的问题，本发明提供了一种利用时空信息的ASL序列PV校正方法，该方法充分利用ASL序列中的时间和空间信息，对序列进行有效PV校正，从而大幅提高ASL图像质量。

本发明是通过以下技术方案来实现：

一种利用时空信息的磁共振动脉自旋标记序列的部分容积校正方法，包括以下步骤：

(1)采集被试者的MRI数据，包括结构像和动脉自旋标记序列；

(2)将结构像和动脉自旋标记序列进行配准；

(3)利用SPM软件对结构像进行分割，获取灰质、白质和脑脊液的概率分布图像；

(4)利用线性回归方法对动脉自旋标记序列进行空间部分容积校正，并利用校正结果作为最大期望算法的初始值，在时间方向上进行部分容积校正，充分利用动脉自旋标记序列的空间和时间信息，获得准确的校正结果。

步骤(1)中，采集的结构像为T1或T2序列，以及动脉自旋标记序列。

步骤(2)所述将结构像和动脉自旋标记序列进行配准，是利用SPM提供的MNI坐标系作为中间值来进行配准。

步骤(4)利用动脉自旋标记的四维数据，沿着时间轴，将具有相同三维空间坐标i的体素构成一个时间向量；假定该向量中所有元素相互独立，且其中的灰质和白质符合高斯分布，利用最大期望算法估计混合组织中灰质和白质灌注信号，具体求解方法如下：

对于某一空间体素i，磁矩表示为：

ΔM_i＝P_iGMΔM_iGM+P_iWMΔM_iWM(1)