【技术实现步骤摘要】
本专利技术涉及磁共振成像，尤其是一种无导航多激发扩散磁共振成像的伪影校正方法及装置。
技术介绍
1、扩散磁共振成像(diffusion mri)是一种较新的医学成像技术，可通过检测活体组织中水分子的布朗运动受限制的方向和程度，获取相关的组织微观结构、神经纤维通路和其他重要的解剖结构或生理功能信息。扩散磁共振成像对运动非常敏感，非常容易受运动伪影影响，这在微观上是由扩散驱动的水分子位移引起的，在宏观上是由头部运动、血流和呼吸引起的相位改变。为了减少这种运动敏感性的影响，传统的dmri采用单次激发(single-shot)平面回波序列,它可以在一次激发后完成一幅图像的采集，成像速度很快，同时对受试的运动不敏感。但该序列存在明显的缺陷：由于相位编码方向带宽很窄，主磁场不均匀引起的图像畸变比较严重；由于信号衰减的影响，数据采集窗的长度受限于t2*值，无法进行高分辨率成像。采用多次激发(multi-shot)平面回波序列可以有效减小图像畸变，从而获得高空间分辨率扩散磁共振图像。但由于不同激发间的相位变化，又会导致更复杂的伪影，这会影响扩散加权磁共振成像的图像质量和后续的定量分析。
技术实现思路
1、有鉴于此，本专利技术实施例的主要目的在于提出一种无导航多激发扩散磁共振成像的伪影校正方法及装置，以提高扩散加权磁共振成像的图像质量，优化定量分析结果。
2、为实现上述目的，本专利技术实施例的一方面提供了一种无导航多激发扩散磁共振成像的伪影校正方法，所述方法包括：
3、获取扩
4、对所述扩散磁共振数据进行一维线性相位校正，得到初步校正后的扩散磁共振数据的b0数据以及初步校正后的扩散磁共振数据的b1数据；
5、将所述初步校正后的扩散磁共振数据的b1数据按激发数分割成若干个不同的k空间，得到分割后的b1数据的k空间数据，其中，所述分割后的b1数据的k空间数据包括每个激发对应的单次激发的b1数据的k空间数据；
6、根据所述初步校正后的扩散磁共振数据的b0数据的k空间数据以及所述分割后的b1数据的k空间数据，获取第一分块汉克尔矩阵；
7、根据所述分割后的b1数据的k空间数据，获取第二分块汉克尔矩阵；
8、对所述第一分块汉克尔矩阵以及所述第二分块汉克尔矩阵进行奇异值分解和奇异值截断，得到目标汉克尔矩阵；
9、根据所述目标汉克尔矩阵对未采样的所述分割后的b1数据的k空间数据进行估计，得到目标k空间数据；
10、将所述目标k空间数据进行逆傅里叶变换，并沿通道维度合并，得到重建图像。
11、在一些实施例中，所述根据所述初步校正后的扩散磁共振数据的b0数据的k空间数据以及所述分割后的b1数据的k空间数据，获取第一分块汉克尔矩阵，包括：
12、在所述初步校正后的扩散磁共振数据的b0数据的k空间数据中构建第三分块汉克尔矩阵；
13、在所述每个激发对应的单次激发的b1数据的k空间数据中构建第四分块汉克尔矩阵；
14、将所述第三分块汉克尔矩阵与多个所述第四分块汉克尔矩阵沿卷积核维度进行拼接，得到第一分块汉克尔矩阵。
15、在一些实施例中，所述根据所述分割后的b1数据的k空间数据，获取第二分块汉克尔矩阵，包括：
16、对所述每个激发对应的单次激发的b1数据的k空间数据与所述分割后的b1数据的k空间数据中的其他激发对应的单次激发的b1数据的k空间数据进行拟合，得到第二分块汉克尔矩阵。
17、在一些实施例中，所述对所述每个激发对应的单次激发的b1数据的k空间数据与所述分割后的b1数据的k空间数据中的其他激发对应的单次激发的b1数据的k空间数据进行拟合，包括：
18、获取在所述每个激发对应的单次激发的b1数据的k空间数据中构建的第四分块汉克尔矩阵；
19、将所述第四分块汉克尔矩阵沿通道维度进行拼接。
20、在一些实施例中，所述对所述第一分块汉克尔矩阵以及所述第二分块汉克尔矩阵进行奇异值分解和奇异值截断，得到目标汉克尔矩阵，包括：
21、设置目标秩；
22、对所述第一分块汉克尔矩阵进行奇异值分解，得到第一奇异值；
23、对所述第二分块汉克尔矩阵进行奇异值分解，得到第二奇异值；
24、根据所述目标秩，对所述第一奇异值进行奇异值截断，得到第一目标分块汉克尔矩阵；
25、根据所述目标秩，对所述第二奇异值进行奇异值截断，得到第二目标分块汉克尔矩阵；
26、其中，所述目标汉克尔矩阵包括所述第一目标分块汉克尔矩阵以及所述第二目标分块汉克尔矩阵。
27、在一些实施例中，所述根据所述目标汉克尔矩阵对未采样的所述分割后的b1数据的k空间数据进行估计，得到目标k空间数据，包括：
28、获取未采样的所述分割后的b1数据的k空间数据对应的若干个所述目标汉克尔矩阵的元素；
29、获取所述元素的平均值；
30、将所述元素的平均值作为中间k空间数据；
31、判断所述中间k空间数据是否收敛，若收敛，则将所述中间k空间数据作为目标k空间数据；
32、若不收敛，则返回根据所述初步校正后的扩散磁共振数据的b0数据的k空间数据以及所述分割后的b1数据的k空间数据，获取第一分块汉克尔矩阵这一步骤。
33、为实现上述目的，本专利技术实施例的另一方面提出了一种无导航多激发扩散磁共振成像的伪影校正装置，所述装置包括：
34、第一模块，用于获取扩散磁共振数据；
35、第二模块，用于对所述扩散磁共振数据进行一维线性相位校正，得到初步校正后的扩散磁共振数据的b0数据以及初步校正后的扩散磁共振数据的b1数据；
36、第三模块，用于将所述初步校正后的扩散磁共振数据的b1数据按激发数分割成若干个不同的k空间，得到分割后的b1数据的k空间数据，其中，所述分割后的b1数据的k空间数据包括每个激发对应的单次激发的b1数据的k空间数据；
37、第四模块，用于根据所述初步校正后的扩散磁共振数据的b0数据的k空间数据以及所述分割后的b1数据的k空间数据，获取第一分块汉克尔矩阵；
38、第五模块，用于根据所述分割后的b1数据的k空间数据，获取第二分块汉克尔矩阵；
39、第六模块，用于对所述第一分块汉克尔矩阵以及所述第二分块汉克尔矩阵进行奇异值分解和奇异值截断，得到目标汉克尔矩阵；
40、第七模块，用于根据所述目标汉克尔矩阵对未采样的所述分割后的b1数据的k空间数据进行估计，得到目标k空间数据；
41、第八模块，用于将所述目标k空间数据进行逆傅里叶变换，并沿通道维度合并，得到重建图像。
42、本专利技术实施例还提供了一种电子设备，该电子设备包括处理器以及存储器；存储器存储有程序；处理器执行程序以执行前述的一种无导航多激发扩散磁共振成像的伪影校正方法。
<本文档来自技高网...【技术保护点】
1.一种无导航多激发扩散磁共振成像的伪影校正方法，其特征在于，包括：
2.根据权利要求1所述的一种无导航多激发扩散磁共振成像的伪影校正方法，其特征在于，所述根据所述初步校正后的扩散磁共振数据的b0数据的k空间数据以及所述分割后的b1数据的k空间数据，获取第一分块汉克尔矩阵，包括：
3.根据权利要求1所述的一种无导航多激发扩散磁共振成像的伪影校正方法，其特征在于，
4.根据权利要求3所述的一种无导航多激发扩散磁共振成像的伪影校正方法，其特征在于，所述对所述每个激发对应的单次激发的b1数据的k空间数据与所述分割后的b1数据的k空间数据中的其他激发对应的单次激发的b1数据的k空间数据进行拟合，包括：
5.根据权利要求1所述的一种无导航多激发扩散磁共振成像的伪影校正方法，其特征在于，所述对所述第一分块汉克尔矩阵以及所述第二分块汉克尔矩阵进行奇异值分解和奇异值截断，得到目标汉克尔矩阵，包括：
6.根据权利要求1所述的一种无导航多激发扩散磁共振成像的伪影校正方法，其特征在于，所述根据所述目标汉克尔矩阵对未采样的所述分割后的b1数据的
7.一种无导航多激发扩散磁共振成像的伪影校正装置，其特征在于，包括：
8.一种电子设备，其特征在于，包括处理器以及存储器；
9.一种计算机可读存储介质，其特征在于，所述存储介质存储有程序，所述程序被处理器执行实现如权利要求1至6中任一项所述的方法。
10.一种计算机程序产品，包括计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6任一项所述的方法。
...【技术特征摘要】
1.一种无导航多激发扩散磁共振成像的伪影校正方法，其特征在于，包括：
2.根据权利要求1所述的一种无导航多激发扩散磁共振成像的伪影校正方法，其特征在于，所述根据所述初步校正后的扩散磁共振数据的b0数据的k空间数据以及所述分割后的b1数据的k空间数据，获取第一分块汉克尔矩阵，包括：
3.根据权利要求1所述的一种无导航多激发扩散磁共振成像的伪影校正方法，其特征在于，
4.根据权利要求3所述的一种无导航多激发扩散磁共振成像的伪影校正方法，其特征在于，所述对所述每个激发对应的单次激发的b1数据的k空间数据与所述分割后的b1数据的k空间数据中的其他激发对应的单次激发的b1数据的k空间数据进行拟合，包括：
5.根据权利要求1所述的一种无导航多激发扩散磁共振成像的伪影校正方法...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。