一种改善并行磁共振重建图像质量的变密度数据采样方法与流程

1.本发明涉及并行磁共振成像技术领域，尤其涉及一种改善并行磁共振重建图像质量的变密度数据采样方法。


背景技术：

2.并行磁共振成像(pmri)是使用多个通道来接收线圈信息进行数据采集的技术，这些线圈可以“同时收集数据”。和传统的磁共振成像技术相比，pmri技术通过使用多个线圈采集数据，提高了数据采集的速度，因此在过去20多年里广泛应用于减少磁共振成像的时间。磁共振成像系统中的数据通常是以k空间形式进行存储的，磁共振成像系统先围绕患者采集数据，这些数据称之为k空间数据，然后将这些k空间数据通过重建算法就可以得到我们能看懂的磁共振图像。
3.一般情况下，医疗磁共振成像系统的所用图像的尺寸大小256
×
256，k空间数据的尺寸大小同样为256
×
256。磁共振系统采集数据的方式一般为：在每个相位编码方向每采集一次就填充一行数据，故在k空间尺寸为256
×
256的磁共振系统中，在256个相位编码方向分别采样，采样时每行填充256个k空间数据。pmri的成像时间取决于采集相位编码方向的数量，即采集的数据量，因此减少k空间采样的数据就能减少采集时间。但是减少采样数据会导致重建过程中出现的严重的叠影，这是必须要解决的问题。
4.在现实生活中，通过完整采样得到k空间采样数据，就可以通过重建算法还原得到完整的图像。但是由于完整采样耗费的时间太长，对于患者而言非常不方便，因此我们必须尽可能的减少采集数据所耗费的时间。为了达到加速成像的目的，我们会对k空间的数据进行欠采样。我们一般定义一个减小因子r，即在k空间进行数据采集时，在相位编码方向上每r行采集一次数据。当r＝2时，我们所消耗的时间只有完整采集的1/2。当r＝4时，我们所消耗的时间只有完整采集时的1/4。这种方式能大大缩短采集所消耗的时间。但是，这样会大大降低重建图像的质量，从而降低磁共振成像系统的实际应用价值。国内外研究者们围绕着这一方面做了大量的工作，并且取得了一定的成果。
5.目前将欠采样数据还原为图像并且尽可能保存完整的算法有很多，主流的方法大致分为两类：基于图像域的重建算法和基于k空间的成像算法。基于图像域的算法主要是先将欠采样数据重建为混叠图像，再对混叠图像进行处理以还原出清晰的图像，代表算法为sense；基于k空间的成像算法主要是先在k空间中估计出缺失的数据，然后再进行重建，代表算法有grappa,nl
‑
grappa等。这两类算法的采样模式都需要在k空间中心完整采样一定的数据，图像域重建算法是为了获得更多的图像信息，k空间重建算法是为了用这些数据估算缺失值。
6.有资料显示，目前运用的最广泛的就是均匀欠采样模式，均匀欠采样模式是欠采样时选择固定的减小因子orf以获取数据。grappa，nl
‑
grappa和其他方法的采样模式相似。在对欠采样数据进行采样时，需要在中央k空间中以全采样的形式获得一些acs数据。grappa，nl
‑
grappa使用这些acs线来计算线性拟合的系数。例如当orf＝4时，在数据采集
时，会在在采集空间的中心区域进行全采样，在外部区域每4行采集1次。
7.中国专利文献cn104063886b公开了一种“基于稀疏表示和非局部相似的核磁共振图像重建方法”。采用了具体步骤是：首先，对磁共振图像对应的傅里叶变换系数，采用变密度随机下采样的方式采样，对采样数据进行傅里叶反变换得到用于重建的初始参考图像；其次，对参考图像进行分块，得到的每类图像子块具有相似的结构特征，并得到每类图像子块对应的字典和图像子块的稀疏表示系数；最后，利用图像子块的非局部相似性估计原始图像，对图像子块的稀疏系数进行约束，结合图像在小波域的稀疏性，通过混合正则项求解模型进行迭代重建。上述技术方案采集部分数据以缩短采集时间导致图像质量下降。


技术实现要素：

8.本发明主要解决原有的技术方案采集部分数据以缩短采集时间导致图像质量下降的技术问题，提供一种改善并行磁共振重建图像质量的变密度数据采样方法，充分利用“图像的大部分信息都包含在k空间中心”的这一理论，在k空间中心采集更多的数据，或者还原出更准确的数据，提升图像的质量，具体采用保持采集acs数据的数量与传统方法一致，增大最外围k空间的orf值，以减少acs区域两侧orf的值，由于在acs区域使用较低orf进行采集，还原出来的数据会比传统方法准确，并将改进的采样方法使用nl
‑
grappa进行还原，在保持了总采集时间不变，采集acs数据不变的情况下，重建后能获得更好的效果，重建后的伪影也明显更轻。
9.本发明的上述技术问题主要是通过下述技术方案得以解决的：本发明包括以下步骤：
10.s1进行磁共振成像的数据采集；k空间中心区域和以往采集模式一样，采集同样数量的acs数据。中心区域外，越靠近减小因子r越小，越远离减小因子r越大。
11.s2使用nl
‑
grappa利用全采样的acs线计算出所使用不同r的拟合系数；再使用这些拟合系数去重建对应r区域的数据。
12.s2进行逆傅立叶变换重建图像。得到所需核磁共振图像。
13.作为优选，所述的步骤s1在k空间中心区域进行完整采样，采集acs数据；k空间中心区域外，越靠近中心区域，使用越小的减小因子r进行采集；越远离中心区域，使用越大的减小因子r进行采集。
14.作为优选，所述的步骤s1具体包括：
15.s1.1进行采集数据参数定义；
16.s1.2使用变密度采集模式进行采集数据；
17.s1.3在分割完的区域内使用不同的r进行欠采样。
18.作为优选，所述的步骤s1.1包括：
19.定义标准减小因子r,r∈(1,2,3
…
)，在不使用r时，磁共振系统每行都采集数据，即全采样；在使用r时，磁共振系统每r行采集一次数据，即欠采样；
20.定义acs数据量n
acs
，表示k空间中心区域进行全采样，即r＝1，获得acs数据的行数；
21.定义磁共振系统采集数据的总量为sum，sum＝n
acs
n
r
，其中，n
r
表示使用以往采集模式使用r欠采样时，在acs区域外采集到的数据行数，
[0022][0023]
其中，n
full
为整个k空间相位编码的行数。
[0024]
作为优选，所述的步骤s1.2具体包括在变密度采集模式中，在acs外的空间使用多种不同的r进行采样，使用n
r_more
代表使用n种r欠采样时采集到的数据行数，则总采集行数变为sum＝n
acs
n
r_more
；
[0025]
首先将k空间的总数据量n
full
分为除acs区域外的n个区域，在n个不同区域中数据行使用不同的r进行采样，公式表示为：
[0026][0027]
其中表示使用第i种r欠采样时该区域k空间相位编码的总行数，这些区域是除了n
acs
外的所有k空间相位编码的行数，然后能够得到每个区域的欠采样数据量：
[0028][0029]
其中表示使用第i种r欠采样时，在该区域采集到的数据行数，最终能够得到acs外欠采样的数据量：
[0030][0031]
作为优选，所述的步骤s1.3将以往采集模式使用的单一减小因子r定义为r
nom
，由以往采集模式sum＝n
acs
n
r
和本方案sum＝n
acs
n
r_more
得到n
r_more
＝n
r
，即本方案在acs区域外欠采样的数据量与以往采集模式相同，令靠近k空间中心区域的r减小时，将其他区域的r增大，保持欠采样总数据量n
r_more
不变。根据图像的大部分信息都存储在k空间中心的理论，越靠近k空间中心区域(即acs区域)的r应尽可能小，即小于r
nom
，目的是为了获得更多的图像信息。因此将靠近k空间中心区域的r减小时，需要将其他区域的r增大，以保持欠采样总数据量n
r_more
不变。距k空间中心区域越远，图像包含的信息越少，所以将k空间最外围区域的r增大。即减小靠近k空间中心的r，增大远离k空间中心的r。
[0032]
本方案是在采集总数据量sum和acs线数量n
acs
保持不变的基础上提高图像质量。
[0033]
作为优选，所述的步骤s2具体包括：
[0034]
s2.1使用nl
‑
grappa算法的估算权重方法，通过采集到的acs数据估算出各个欠采样区域r中的权重；
[0035]
s2.2使用估算出来的权重估算出缺失的数据，并存入相应的位置，将k空间中缺失的数据补全完整后得到相应的完整k空间数据。
[0036]
作为优选，所述的步骤s3具体包括：对于直角坐标系下的二维k空间数据，利用逆傅立叶变换方法，得到磁共振图像，计算方法如下所示:
[0037][0038]
其中，s(k
x
,k
y
)为k空间采样数据，ρ(x,y)为经过重建后得到的图像。
[0039]
本发明的有益效果是：充分利用“图像的大部分信息都包含在k空间中心”的这一理论，在k空间中心采集更多的数据，或者还原出更准确的数据，提升图像的质量，具体采用保持采集acs数据的数量与传统方法一致，增大最外围k空间的orf值，以减少acs区域两侧orf的值，由于在acs区域使用较低orf进行采集，还原出来的数据会比传统方法准确，并将改进的采样方法使用nl
‑
grappa进行还原，在保持了总采集时间不变，采集acs数据不变的情况下，重建后能获得更好的效果，重建后的伪影也明显更轻。
附图说明
[0040]
图1是本发明的一种流程图。
[0041]
图2是本发明的一种以往r＝4的数据采集模式图。
[0042]
图3是本发明的一种数据采集模式图
[0043]
图4是本发明的一种k空间欠采样数据图。
[0044]
图5是本发明的一种采集模式示意图。
[0045]
图6是本发明的一种k空间数据图。
[0046]
图7是本发明的一种还原完整后的k空间数据图。
[0047]
图8是本发明的一种使用逆傅里叶变换二维k空间数据图。
具体实施方式
[0048]
下面通过实施例，并结合附图，对本发明的技术方案作进一步具体的说明。
[0049]
实施例：本实施例的一种改善并行磁共振重建图像质量的变密度数据采样方法，如图1所示，包括
[0050]
步骤1、对于磁共振成像的数据采集模式，我们在k空间中心区域和以往采集模式一样，采集同样数量的acs数据。但在k空间中心区域外，越靠近中心区域，使用较小减小因子r进行采集；远离中心区域，使用较大减小因子r进行采集，以保持采集时间不变。
[0051]
步骤2、使用nl
‑
grappa利用全采样的acs线计算出所使用不同r的拟合系数，再使用这些拟合系数去重建对应r区域的数据
[0052]
步骤3、使用逆傅立叶变换将k空间数据变为图像，得到所需核磁共振图像。
[0053]
步骤1所述的k空间中心区域和以往采集模式一样，采集同样数量的acs数据。中心区域外，越靠近减小因子r越小，越远离减小因子r越大，具体操作如下：
[0054]1‑
1.定义标准减小因子r,r∈(1,2,3
…
)。在不使用r时，磁共振系统每行都采集数据，即全采样；在使用r时，磁共振系统每r行采集一次数据，即欠采样。定义acs数据量n
acs
，表示k空间中心区域进行全采样(即r＝1)获得acs数据的行数。如图2所示能形象说明以往的采集模式，虚线为未采集数据，实现为采集数据。在图中r＝4，在acs区域外每4行采集一次数据。
[0055]
定义磁共振系统采集数据的总量为sum，sum＝n
acs
n
r
，其中，n
r
表示使用以往采集模式使用r欠采样时，在acs区域外采集到的数据行数。
[0056][0057]
其中，n
full
为整个k空间相位编码的行数。
[0058]1‑
2.使用本技术方案中的变密度采集模式进行采集数据，在变密度采集模式中，我们在acs外的空间使用多种不同的r进行采样。使用n
r_more
代表我们使用n种r欠采样时采集到的数据行数。则总采集行数变为sum＝n
acs
n
r_more
。
[0059]
我们首先需要将k空间的总数据量n
full
分为除acs区域外的n个区域，在n个不同区域中数据行使用不同的r进行采样，公式可以表示为：
[0060][0061]
其中表示使用第i种r欠采样时该区域k空间相位编码的总行数。这些区域是除了n
acs
外的所有k空间相位编码的行数。然后我们能够得到每个区域的欠采样数据量：
[0062][0063]
其中表示使用第i种r欠采样时，在该区域采集到的数据行数。最终能够得到acs外欠采样的数据量：
[0064][0065]
例：以k空间尺寸为256
×
256，r
nom
＝4，n
acs
＝16为例。其中n
full
＝256。其中外围区域欠采样总行数n
r
计算方式为：
[0066][0067]
本方案采集模式以使用3种减小因子为例，在靠近k空间中心的区域我们将r减少为r
nom
＝2，在远离k空间中心的区域我们将r
nom
增大为r3＝6，在中间区域r
nom
保持不变，r2＝4。保持外围采集总数据量与n
r
一致，则此时需要
[0068]
n
r_more
＝n
r1
n
r2
n
r3
＝n
r
＝60
[0069]
同时这3个区域的总行数为除去acs外的所有k空间区域，即
[0070][0071]
这可以得到许多种分配的方案，我们拟定即r1、r2、r3分别采集相同的数据量。如图5所示表示了以往的采集模式(a)和本方案的采集模式(b),其中白色行为采集到的数据行，黑色为未采集到的数据。
[0072]
采集到的k空间数据如图6所示。(a)为完整的k空间数据，(b)为以往的采集方式，(c)为本方案采集模式后的k空间数据。
[0073]1‑
3.在分割完的区域内使用不同的r进行欠采样。我们将以往采集模式使用的单一减小因子r定义为r
nom
。根据图像的大部分信息都存储在k空间中心的理论，越靠近k空间中心区域(即acs区域)的r应尽可能小，即小于r
nom
，目的是为了获得更多的图像信息。
[0074]
本方案是在采集总数据量sum和acs线数量n
acs
保持不变的基础上提高图像质量。由以往采集模式sum＝n
acs
n
r
和本方案sum＝n
acs
n
r_more
得到n
r_more
＝n
r
。即本方案在acs区域
外欠采样的数据量应与以往采集模式相同。
[0075]
因此将靠近k空间中心区域的r减小时，我们需要将其他区域的r增大，以保持欠采样总数据量n
r_more
不变。距k空间中心区域越远，图像包含的信息越少，所以我们将k空间最外围区域的r增大。即减小靠近k空间中心的r，增大远离k空间中心的r。
[0076]
如图3所示能形象说明本方案的采集模式。和图1的以往以单一的r
nom
＝4采集模式相比，在靠近acs的区域我们将r减少为r
nom
＝2，在远离acs的区域我们将r
nom
增大为r3＝6，在中间区域r
nom
保持不变，r2＝4。图中表示使用各自的r采集的区域。实线为最终采集到的数据，虚线为未采集到的数据。
[0077]
步骤2所述使用nl
‑
grappa利用全采样的acs线计算出所使用不同r的拟合系数具体操作如下：
[0078]2‑
1.使用nl
‑
grappa算法的估算权重方法，通过采集到的acs数据估算出各个欠采样区域r中的权重。
[0079]2‑
2.使用估算出来的权重估算出缺失的数据，并存入相应的位置，将k空间中缺失的数据补全完整后得到相应的完整k空间数据。
[0080]
以步骤1中的例子为例，通过nl
‑
grappa使用采集到的16行acs数据分别计算r1＝2、r2＝4、r3＝6的权重，然后使用权重去分别还原对应区域的欠采样k空间数据。还原完整后的k空间数据如图7所示。
[0081]
根据步骤2重建图2本方案采集到的k空间欠采样数据如图3所示。
[0082]
步骤3所述的逆傅立叶变换重建图像具体如下：
[0083]3‑
1.对于直角坐标系下的二维k空间数据，可以利用逆傅立叶变换方法，可得到磁共振图像。计算方法如下所示:
[0084][0085]
其中，s(k
x
,k
y
)为k空间采样数据，ρ(x,y)为经过重建后得到的图像。
[0086]
将步骤2中得到的二维k空间数据使用逆傅里叶变换后得到的图像如图8(c)所示。其中图(a)为原始图像，图(b)为使用以往采集模式r
nom
＝4，n
acs
＝16情况下通过nl
‑
grappa重建后得到的图像，图(c)为使用本方案的采集方法，与以往采集模式采集相同数据量下通过nl
‑
grappa重建后得到的图像。
[0087]
grappa:一种基于k空间进行数据还原的算法。
[0088]
nl
‑
grappa：非线性grappa。grappa算法的改进版，效果更佳准确。
[0089]
k空间：k空间是寻常空间在傅利叶转换下的对偶空间，主要应用在磁振造影的成像分析。磁共振仪器采集到的的原始数据，就是k空间数据。一个256x256的k空间数据矩阵，经过逆傅立叶变换得到一个256x256尺寸的图像。
[0090]
acs：在k空间中心，全采样用以计算还原权重的数据。
[0091]
全采样：即磁共振系统采集数据并填入k空间时，每行数据都进行采集采集。
[0092]
欠采样：采集数据时，多行才采集一次。
[0093]
orf：减小因子。用于确定欠采样时所间隔的行数。当orf＝4时，每4行采集一次。
[0094]
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替
代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
[0095]
尽管本文较多地使用了k空间、减小因子等术语，但并不排除使用其它术语的可能性。使用这些术语仅仅是为了更方便地描述和解释本发明的本质；把它们解释成任何一种附加的限制都是与本发明精神相违背的。