基于两层紧框架稀疏模型的并行磁共振成像方法、装置及计算机可读介质制造方法及图纸

本发明专利技术涉及基于两层紧框架稀疏模型的并行磁共振成像方法、装置及计算机可读介质。提供了一种基于两层紧框架稀疏模型的并行磁共振成像方法，包括以下步骤：S1)基于欠采样所得的每个通道的图像数据，创建基于两层紧框架稀疏模型的压缩感知并行成像模型函数；S2)针对所创建的压缩感知并行成像模型函数进行优化迭代；以及S3)根据函数优化迭代过程所求得的目标图像数据重建图像。本发明专利技术的基于两层紧框架稀疏模型的并行磁共振成像方法解决了自适应力不高、算法复杂度大、稀疏表达没有充分利用等问题，从而提高成像的速度和精度。

Parallel magnetic resonance imaging method, device and computer readable medium based on two layer tight frame sparse model

The invention relates to a parallel magnetic resonance imaging method, device and computer readable medium based on a two layer tight frame sparse model. A parallel magnetic resonance imaging method based on the two layer tight frame sparse model is provided, including the following steps: S1) based on the image data of each channel under the undersampled data, the compression perception parallel imaging model function based on the two layer tight frame sparse model is created; the parallel imaging model function of the compressed sensing is created. Optimize iteration; and S3) reconstruct the image according to the target image data obtained by the function optimization iteration process. The parallel magnetic resonance imaging (MRI) method based on two layer tight frame sparse model solves the problems of low self adaptability, large algorithm complexity and lack of full use of sparse expression, so as to improve the speed and accuracy of imaging.

【技术实现步骤摘要】
基于两层紧框架稀疏模型的并行磁共振成像方法、装置及计算机可读介质
本专利技术涉及磁共振成像
，尤其涉及基于两层紧框架稀疏模型的并行磁共振成像方法、装置及计算机可读介质。
技术介绍
目前，磁共振成像已经发展成为生物学工程领域重要技术之一。然而磁共振成像速度较慢，且在成像过程中人体器官的运动、被检查者的移动等都会影像成像的清晰度。因此，缩短成像时间对磁共振成像有着重要的意义。在传统的成像方法中，重建方法缺乏自适应能力和有效的稀疏表达，特别是在欠采倍数较大和加噪的情况下。压缩感应磁共振成像(Compressedsensingmagneticresonanceimaging，CSMRI)是用于加速MRI扫描的基于信号处理的技术。压缩感应CS利用了MR图像的稀疏性，并允许CSMRI从稀疏采样的K空间数据中恢复MR图像，CSMRI的经典公式可以写成：其中，u∈lQ×1和f∈lP×1分别表示MR图像及其相应的欠采样原始K空间数据，FP∈lQ×P表示欠采样的傅里叶编码矩阵，其中P＜＜Q，并且||Wu||1是一种分析模型，它在l1范数约束下利用变换W∈lQ×Q稀释图像。P和Q表示图像像素数量和测得数据的数量。然而，由于缺乏适应性或缺乏稀疏性，效率是以牺牲精度为代价的，特别是高度欠采样的噪声测量。为解决这个问题，目前已提出一些方法，主要包括以下三个方向。第一个方向主要是利用全局变换和冗余变换来分析稀疏变换的磁共振成像。经典的算法有基于压缩感知的全局变换和小波变换，这些方法重建模型较简单，其图像重建精度不高。第二个方向是利用自适应字典稀疏表达磁共振成像，其重建方法有DLMRI,BPFA,TBMDU等，它们利用字典来自适应学习图像的结构信息。这些方法的去噪能力较强，然而，它们的计算量非常复杂，稀疏性仍然直接限制在一个层上。第三个方向是在前面两个方向改善的，其目标是利用分析和合成稀疏模型的优点。如平衡的紧框架模型，它包含一个惩罚项来缩小分析和合成模型的距离。然而，尽管这种方法有很好的数学解释，但对于简单的一层，其稀疏的能力有待进一步提高。为了提高稀疏度，基于小波驱动的字典学习和基于全局变换的字典学习自适应从分析变换中代表稀疏系数。尽管重建效果提高了，但字典学习的计算量非常复杂。针对以上技术的缺点，本专利技术提出了一种基于两层紧框架稀疏模型的并行磁共振成像方法。它是在稀疏的图像中通过一个固定的紧框架和一个自适应的学习紧框架重建图像，它采用了三层的Bregman算法有效地解决重建模型的优化问题，以实现快速且高质量的磁共振成像。
技术实现思路
如上所述，目前磁共振重建的方法有：(1)基于压缩感知的全局变换和小波变换的磁共振成像，其重建模型较简单，其图像重建精度不高。(2)基于自适应字典稀疏表达磁共振成像，其计算量非常复杂，稀疏性仍然直接限制在一个层上。(3)基于分析和合成模型的磁共振成像，其稀疏性有待进一步提高。为了解决上述问题，针对目前基于压缩感知的并行磁共振成像方法的缺陷，本专利技术提出了一种基于两层紧框架稀疏模型的并行磁共振成像方法。该方法是在稀疏的图像中通过一个固定的紧框架和一个自适应的学习紧框架重建图像，采用了三层的Bregman算法有效地解决重建模型的优化问题，以实现快速且高质量的磁共振成像。本专利技术的方法解决了自适应力不高、算法复杂度大、稀疏表达没有充分利用等问题，从而提高成像的速度和精度。更具体而言，本专利技术提供了一种基于两层紧框架稀疏模型的并行磁共振成像方法，包括以下步骤：S1)基于欠采样所得的每个通道的图像数据，创建基于两层紧框架稀疏模型的压缩感知成像模型函数；S2)针对所创建的压缩感知成像模型函数进行优化迭代；以及S3)根据函数优化迭代过程所求得的目标图像数据来重建图像。在其中的一个实施例中，在步骤S1中创建的基于两层紧框架稀疏模型的压缩感知成像模型函数为：其中，J是采样通道个数，Wa,j表示第j通道的第一层固定的紧框架，Wb,j表示第j通道的第二层数据驱动的紧框架，Γ表示紧框架系统，vj∈lQ×1和fj∈lP×1分别表示第j通道的灵敏度编码图像和第j通道的原始k空间数据，FP∈lQ×P表示欠采样的傅里叶编码矩阵，其中P＜＜Q，并且其中，vj＝sju，u表示并行磁共振成像的目标图像，sj表示并行磁共振成像的第j通道的灵敏度。在其中的一个实施例中，所述压缩感知成像模型函数用下式的卷积模型近似：式中，J是采样通道个数，an,j是第j通道固定的核，bm,j表示第j通道学习的自适应核，vj∈lQ×1和fj∈lP×1分别表示第j通道的灵敏度编码图像和第j通道的原始k空间数据，FP∈lQ×P表示欠采样的傅里叶编码矩阵，其中P＜＜Q，并且其中，vj＝sju，u表示并行磁共振成像的目标图像，sj表示并行磁共振成像的第j通道的灵敏度。在其中的一个实施例中，采用三层Bregman迭代数值算法对所述压缩感知成像模型函数进行优化迭代，所述三层Bregman迭代数值算法的一个迭代循环包括如下步骤：S21)引入Bregman参数c，将第一层的Bregman迭代数值算法设为：其中，J是采样通道个数，k表示迭代次数，Wa,j表示第j通道的第一层固定的紧框架，Wb,j表示第j通道的第二层数据驱动的紧框架，Γ表示紧框架系统，vj∈lQ×1和fj∈lP×1分别表示第j通道的灵敏度编码图像和第j通道的原始k空间数据，FP∈lQ×P表示欠采样的傅里叶编码矩阵，其中P＜＜Q，并且其中，vj＝sju，u表示并行磁共振成像的目标图像，sj表示并行磁共振成像的第j通道的灵敏度；S22)定义第一辅助变量va,j＝Wa,jvj，得到第二层Bregman迭代：其中，α,μ表示权重参数，d为辅助变量，并且va,j的更新满足：所述紧框架Wa,j满足：归一化后全采样傅里叶编码矩阵F满足：FHF＝I由此得到Fvj(kx,ky)：式中Ωj表示第j通道的欠采样的K空间集；S23)引入第二辅助变量ηj＝Wb,jva,j来分解变量Wb,j和va,j之间的耦合，由此得到第三层Bregman迭代：其中，α,μ表示权重参数，e为引入的新的辅助变量，va,j的更新满足：临时固定数据驱动的紧框架变量Wb,j，按下式更新第二辅助变量ηj:式中shrink(x,a)＝sign(x)max(0,|x|-a)；固定第二辅助变量ηj，通过下式更新数据驱动紧框架变量Wb,j:在其中的一个实施例中，所述方法还包括利用在所述三层Bregman迭代数值算法的迭代过程中得到的结果，按照下式更新并行磁共振成像的目标图像：其中，J是采样通道个数，k表示迭代次数，u表示并行磁共振成像的目标图像，表示迭代优化更新的第j通道的图像。本专利技术的另一方面还提供了一种计算机可读介质，该计算机可读介质具有存储在其中的程序，该程序是计算机可执行的以使计算机执行上述基于两层紧框架稀疏模型的并行磁共振成像方法的各步骤。本专利技术的另一方面还提供了一种基于两层紧框架稀疏模型的并行磁共振成像装置，其包括：压缩感知成像模型函数创建模块，基于欠采样所得每个通道的图像数据来创建基于两层紧框架稀疏模型的压缩感知成像模型函数；优化迭代模块，针对所创建的压缩感知成像模型函数进行优化迭代；以及目标图像重建模块，根据函数优化迭代过程所求得的目标图像数据来重建图本文档来自技高网...

【技术保护点】
1.一种基于两层紧框架稀疏模型的并行磁共振成像方法，其特征在于，包括以下步骤：S1)基于欠采样所得的每个通道的图像数据，创建基于两层紧框架稀疏模型的压缩感知成像模型函数；S2)针对所创建的压缩感知成像模型函数进行优化迭代；以及S3)根据函数优化迭代过程所求得的目标图像数据来重建图像。

【技术特征摘要】
1.一种基于两层紧框架稀疏模型的并行磁共振成像方法，其特征在于，包括以下步骤：S1)基于欠采样所得的每个通道的图像数据，创建基于两层紧框架稀疏模型的压缩感知成像模型函数；S2)针对所创建的压缩感知成像模型函数进行优化迭代；以及S3)根据函数优化迭代过程所求得的目标图像数据来重建图像。2.根据权利要求1所述的基于两层紧框架稀疏模型的并行磁共振成像方法，其特征在于，在步骤S1中创建的基于两层紧框架稀疏模型的压缩感知成像模型函数为：其中，J是采样通道个数，Wa,j表示第j通道的第一层固定的紧框架，Wb,j表示第j通道的第二层数据驱动的紧框架，Γ表示紧框架系统，vj∈lQ×1和fj∈lP×1分别表示第j通道的灵敏度编码图像和第j通道的原始k空间数据，FP∈lQ×P表示欠采样的傅里叶编码矩阵，其中P<<Q，并且其中，vj＝sju，u表示并行磁共振成像的目标图像，sj表示并行磁共振成像的第j通道的灵敏度。3.根据权利要求2所述的基于两层紧框架稀疏模型的并行磁共振成像方法，其特征在于，所述压缩感知成像模型函数用下式的卷积模型近似：其中，J是采样通道个数，an,j是第j通道固定的核，bm,j表示第j通道学习的自适应核，vj∈lQ×1和fj∈lP×1分别表示第j通道的灵敏度编码图像和第j通道的原始k空间数据，FP∈lQ×P表示欠采样的傅里叶编码矩阵，其中P<<Q，并且其中，vj＝sju，u表示并行磁共振成像的目标图像，sj表示并行磁共振成像的第j通道的灵敏度。4.根据权利要求2所述的基于两层紧框架稀疏模型的并行磁共振成像方法，其特征在于，采用三层Bregman迭代数值算法对所述压缩感知成像模型函数进行优化迭代，所述三层Bregman迭代数值算法的一个迭代循环包括如下步骤：S21)引入Bregman参数c，将第一层的Bregman迭代数值算法设为：其中，J是采样通道个数，k表示迭代次数，Wa,j表示第j通道的第一层固定的紧框架，Wb,j表示第j通道的第二层数据驱动的紧框架，Γ表示紧框架系统，vj∈lQ×1和fj∈lP×1分别表示第j通道的灵敏度编码图像和第j通道的原始k空间数据，FP∈lQ×P表示欠采样的傅里叶编码矩阵，其中P<<Q，并且其中，vj＝sju，u表示并行磁共振成像的目标图像，sj表示并行磁共振成像的第j通道的灵敏度；S22)定义第一辅助变量va,j＝Wa,jvj，得到第二层Bregman迭代：其中，α,μ表示权重参数，d为辅助变量，并且va,j的更新满足：所述紧框架Wa,j满足：归一化后全采样傅里叶编码矩阵F满足：FHF＝I由此得到Fvj(kx,ky)：式中Ωj表示第j通道的欠采样的K空间集；S23)引入第二辅助变量ηj＝Wb,jva,j来分解变量Wb,j和va,j之间的耦合，由此得到第三层Bregman迭代：其中，α,μ表示权重参数，e为引入的新的辅助变量，va,j的更新满足：临时固定数据驱动的紧框架变量Wb,j，按下式更新第二辅助变量ηj:式中shrink(x,a)＝sign(x)max(0,|x|-a)；固定第二辅助变量ηj，通过下式更新数据驱动紧框架变量Wb,j:5.根据权利要求4中任一项所述的基于两层紧框架稀疏模型的并行磁共振成像方法，其特征在于，利用在所述三层Bregman迭代数值算法的迭代过程中得到的结果，按照下式更新并行磁共振成像的目标图像：其中，J是采样通道个数，k表示迭代次数，u表示并行磁共振成像的目标图像，表示迭代优化更新的第j通道的图像。6.一种基于两层紧框架稀疏模型的并行磁共振成像装置，其特征在于，包括：压缩感知成像模型函数创建模块，基于欠采...

【专利技术属性】
技术研发人员：王珊珊，梁栋，谭莎，刘新，郑海荣，
申请(专利权)人：中国科学院深圳先进技术研究院，
类型：发明
国别省市：广东,44