磁共振线圈合并系数计算方法、磁共振成像方法及其装置发明专利.docxVIP

磁共振线圈合并系数计算方法、磁共振成像方法及其装置 技术领域 本发明涉及磁共振成像领域，特别涉及一种磁共振线圈合并系数计算方法、磁共 振成像方法及其装置。 背景技术 在磁共振成像技术中，成像的速度是衡量成像方法的一个很重要标准。限制成像 速度的很重要因素是数据采集，以及k空间填充。一般的数据采集方式要采满k空间数据，然 后才能进行重建得到图像。磁共振并行采集重建技术，是利用线圈重组合并的方式，对欠采 样的数据进行填补，利用填补完整的k空间数据进行重建。利用这样的方式，可以根据需求， 只采集一部分k空间数据，不必采满整个k空间。因此这样的方法可以大大加快成像的速度。 比较常用的并行重建方法之一是GRAPPA。传统的GRAPPA算法如图1所示，黑色实点 代表为实际采集的k空间数据；白色空点为欠采样需要填补的数据；灰色实点表示为了计算 线圈合并参数，而适量全采的校准数据。GRAPPA算法认为，图中任意一个空心点可以表示为 周围黑色实点的线性叠加，相当于对多个线圈的数据进行了合并。而线圈合并系数nij(第i 个线圈，第j个位置，如图1)可以通过黑色的实点拟合灰色点来确定。线圈合并系数确定后， 其他白色空心点即可根据求得的合并系数以及黑色实点计算填补得到，从而重建得到完整 的k空间数据。 线圈合并系数又可称为卷积核，在传统方法中，卷积核的计算方向，只加入了相位 编码方向，以及通道方向。为了优化效果，近年来，很多方法引入了其他方向如频率编码方 向，以及k-t中时间方向等。这些维度的引入，使卷积核的信息增加，但是随着过多的引入， 信息会出现冗余，而且由于噪声的影响，使得计算的系数仍然受到影响；而卷积核过大，使 得要拟合的系数增多，也带来计算的不稳定性。而卷积核过小，则数据信息不足，计算过程 不够准确。 发明内容 本发明要解决的问题是提供一种磁共振并行成像线圈合并系数计算方法、磁共振 并行成像方法及其装置，解决在卷积核选取过程中的由于卷积核选取过大造成信息量增 加，信息冗余，并且噪声过多，拟合系数增多，计算不稳定的影响，以及在卷积核选取过小 时，计算数据不足，结果不准确的问题。 为了实现上述目的，本发明提供了一种磁共振并行成像线圈合并系数计算方法， 包括如下步骤： a)采集k空间数据，所述k空间数据包括采集数据和校准数据； b)对所述采集数据进行数据提炼得到提炼数据，所述采集数据在提炼方向上的维 数为n，所述提炼数据在所述提炼方向上的维数为t，所述t、n均为正整数且tn； c)由所述提炼数据和所述校准数据计算得到线圈合并系数。 优选的，所述步骤a)中k空间数据的数据方向包括以下任意一种或多种：频率编码 方向、相位编码方向、通道方向、时间方向、3d扫描中的第二相位编码方向或上述多个方向 合并形成的方向。 优选的，所述步骤b)中数据提炼方法为矩阵降维法。 优选的，所述步骤b)中矩阵降维法包括如下步骤： 将所述采集数据的第一方向向量作为基础向量，所述第一方向与提炼方向垂直； 计算所述采集数据的协方差矩阵，并选取所述协方差矩阵归一化正交特征矢量中 特征值较大的前t个值，组成提炼系数； 由所述提炼系数计算得到提炼数据。 优选的，所述步骤b)中矩阵降维法包括如下步骤： 将所述采集数据的第一方向向量作为基础向量，所述第一方向与提炼方向垂直； 按照相同数据方向将所述采集数据分类得到各类第一方向向量Obk_Sx1、Obk_ Sx2、…、Obk_Sxi； 由所述各类第一方向向量计算各类样本第一方向向量均值Obk_Sm1、Obk_Sm2、…、 Obk_Smi和总样本第一方向向量均值Obk_Sm，并进而计算得到样本类内离散度矩阵Obk_ Ss1、Obk_Ss2、…、Obk_Ssi以及样本间离散度矩阵Obk_Ssb； Obk_Ssi-1*Obk_Ssb的前t个最大特征值对应的特征向量分别为wi1、wi2、…、wit， 将wi1、wi2至wit组合成一个矩阵作为所述提炼数据。 进一步的，所述第一方向为所述采集数据的列向量方向或行向量方向。 本发明还提供了一种磁共振并行成像方法，包括在所述的磁共振并行成像线圈合 并系数计算方法计算得到线圈合并系数后，根据所述线圈合并系数和所述采集数据，填补 得到完整的k空间数据，将所述完整的k空间数据变换到图像域得到图像。 本发明还提供了一种磁共振并行成像装置，包括： 采集单元，适于采集k空间数据，所述k空间数据包括采集数据和校准数据，所述k 空间数据至少包含一个数据方向； 计算单元，适于对所述采集数据进行数据提炼得到提炼数据，所述采集数据在提 炼方向上的维数为n，所述提炼数据在所述提炼方向上的维数为t，所述t、n均为正整数且t n，并由所述提炼数据和所述校准数据计算得到线圈合并系数； 填补单元，适于根据所述线圈合并系数和所述采集数据，填