CN108154484B - 一种基于自适应联合稀疏编码的并行磁共振成像方法、装置及计算机可读介质 - Google Patents

一种基于自适应联合稀疏编码的并行磁共振成像方法、装置及计算机可读介质 Download PDF

Info

Publication number
CN108154484B
CN108154484B CN201711246873.0A CN201711246873A CN108154484B CN 108154484 B CN108154484 B CN 108154484B CN 201711246873 A CN201711246873 A CN 201711246873A CN 108154484 B CN108154484 B CN 108154484B
Authority
CN
China
Prior art keywords
representing
image
channel
formula
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711246873.0A
Other languages
English (en)
Other versions
CN108154484A (zh
Inventor
王珊珊
梁栋
谭莎
刘新
郑海荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201711246873.0A priority Critical patent/CN108154484B/zh
Publication of CN108154484A publication Critical patent/CN108154484A/zh
Application granted granted Critical
Publication of CN108154484B publication Critical patent/CN108154484B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/73—Deblurring; Sharpening
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing
    • G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明提供了一种基于自适应联合稀疏编码的并行磁共振成像方法、装置以及计算机可读介质。本发明方法通过求解一个l2‑lF‑l2，1最小目标函数，l2范数是数据拟合项，lF范数代表稀疏表达误差，l2，1混合范数代表每个通道之间的联合稀疏约束；再采用分而治之的方法和对应的算法更新稀疏矩阵、字典和K空间数据，最后，通过所有通道的均方根和求得重建图像。本方发明方法利用了l2，1范数开发通道的联合稀疏，开发信息稀疏性的同时，可以免校准，且该方法有较强的鲁棒性。

Description

一种基于自适应联合稀疏编码的并行磁共振成像方法、装置 及计算机可读介质
技术领域
本发明涉及磁共振成像领域，尤其涉及一种基于自适应联合稀疏编码的并行磁共振成像方法、装置以及计算机可读介质。
背景技术
目前，为了提高磁共振成像速度，基于压缩感知的并行磁共振成像比较流行。基于如何使用灵敏度信息，可把重建方法大概分为三类：精确校准，自动校准和免校准。免校准方法结合先验信息，在灵敏度信息上具有较强的鲁棒性。自动校准不依赖于准确的预估计灵敏度信息。K空间自校准灵敏度方法有GRAPPA，L1-SPIRiT等，这些方法的图像重建质量较好。但随着加速倍数提高，这些算法的性能随之降低。针对该缺陷，提出了 JSENSE，SparseBLIP等算法，在估计灵敏度信息的同时还能估计欠采图像。然而同时估计是一个非凸问题，对灵敏度的初始值非常敏感，且计算量非常复杂。免校准的重建方法不需要特殊的校准信息来估计灵敏度。它们在 K空间或图像域通过通道间的关系来重建图像。然而，现在提出的免校准重建方法只是直接从转换域利用联合稀疏性。
发明内容
针对现有技术的上述问题，本发明提出了一种基于自适应联合稀疏编码的并行磁共振成像方法及装置，其基于块的稀疏性，能够提高稀疏性，且能在进一步提高加速倍数的同时，保证重建质量。
本发明提供了一种基于自适应联合稀疏编码的并行磁共振成像方法，包括以下步骤：
步骤a：构建一个基于免校准的联合稀疏编码的并行磁共振成像模型即重建模型，所述模型定义为：
Figure BDA0001490895220000021
式中，V代表重建的图像，
Figure BDA0001490895220000022
代表过完备(P＞＞M)的字典，θ表示目标函数，
Figure BDA0001490895220000023
代表稀疏矩阵，
Figure BDA0001490895220000024
表示从第j-th通道的图像提取第l-th图像块，FM代表傅里叶变换，vj表示每个通道的图像，所有通道的图像vj合并在一起组成的矩阵即为V，fj代表K空间数据，
Figure BDA0001490895220000025
代表块提取矩阵，M表示原子的尺寸，N表示图像的向量化尺寸，P表示原子的个数，下标F表示傅里叶变换，λ表示数据拟合项权重，β表示联合稀疏正则化权重， | | · | |2，1表示联合稀疏项，
Figure BDA0001490895220000026
式中，
Figure BDA0001490895220000027
表示第j-th通道的第l-th提取块中第 p-th像素；
步骤b：临时固定X,采用梯度下降的方法求解字典D，
Figure BDA0001490895220000028
式中γ表示学习率，k表示迭代次数，
由于
Figure BDA0001490895220000029
该公式是重建目标函数对字典D的倒数，H 表示共轭转置，
因此获得下式：
Dk+1＝Dk-γ(DX-RlV)XH (3)
步骤c：更新联合稀疏系数X；
步骤d：通过下式更新vj，
Figure BDA0001490895220000031
步骤e：通过下述公式更新K空间数据：
Figure BDA0001490895220000032
其中，
Figure BDA0001490895220000033
表示第j个通道的第k次迭代更新值；
步骤f：对K空间数据进行傅里叶反变换再次更新vj，获得
Figure BDA0001490895220000034
以及
步骤g：根据所有通道的图像vj，获得更新后的图像。
优选地，步骤a中，构建重建模型即求解一个
Figure BDA0001490895220000038
最小目标函数，其中，
Figure BDA0001490895220000039
范数是数据拟合项，
Figure BDA00014908952200000310
范数代表稀疏表达误差，
Figure BDA00014908952200000311
混合范数代表每个通道之间的联合稀疏约束。
优选地，步骤a中，包括初始化
Figure BDA0001490895220000035
的步骤。
优选地，步骤b中，对原子进行归一化处理，以避免字典原子尺度模糊。
优选地，步骤d中，更新vj进一步包括如下步骤：
步骤d1：根据最小二乘法获得下式：
Figure BDA0001490895220000036
式中，
Figure BDA0001490895220000037
是一个对角矩阵，对角线上元素值等于对应位置图像像素重叠的块的个数；当假设这些块映射在图像的边缘，则对角线上元素值全部相等，
Figure BDA0001490895220000041
特殊情况下，当填充的距离为单位距离1时，ω＝M；
步骤d2：图像域通过归一化全采傅里叶编码矩阵转换到傅里叶域，获得：
Figure BDA0001490895220000042
式中，
Figure BDA0001490895220000043
是0和1的对角矩阵，1代表K空间采到的点的位置，ω表示像素重复的个数，IN表示N*N的单位矩阵，
Figure BDA0001490895220000044
表示第j-th通道零填充的傅里叶测量值；
步骤d3：将步骤d2中的公式变换成：
Figure BDA0001490895220000045
式中，Fvj(kx,ky)表示K空间位置(kx,ky)的更新值，
Figure BDA0001490895220000046
表示第 j-th通道的零填充K空间测量值；η表示一个小标量，Ω表示K空间已经采到的点的集合，
Figure BDA0001490895220000047
步骤d4：通过步骤c3的公式(9)中频域中的数据通过逆傅里叶变换求得vj。
优选地，步骤g中，对所有通道的图像vj的平方进行求和，然后开方，获得更新后的图像。
本发明的另一个方面提供了一种基于自适应联合稀疏编码的并行磁共振成像装置，包括：模型构建模块，构建一个基于免校准的联合稀疏编码的并行磁共振成像模型即重建模型，所述模型定义为：
Figure BDA0001490895220000051
式中，V代表重建的图像，
Figure BDA0001490895220000052
代表过完备(P＞＞M)的字典，θ表示目标函数，
Figure BDA0001490895220000053
代表稀疏矩阵，
Figure BDA0001490895220000054
表示从第j-th通道的图像提取第l-th图像块，FM代表傅里叶变换，vj表示每个通道的图像，所有通道的图像vj合并在一起组成的矩阵即为V，fj代表K空间数据，
Figure BDA0001490895220000055
代表块提取矩阵，M表示原子的尺寸，N表示图像的向量化尺寸，P表示原子的个数，下标F表示傅里叶变换，λ表示数据拟合项权重，β表示联合稀疏正则化权重，||·||2，1表示联合稀疏项，
Figure BDA0001490895220000056
式中，
Figure BDA0001490895220000057
表示第j-th通道的第l-th提取块中第 p-th像素；
更新模块，包括：更新联合稀疏编码X的模块、更新字典D的模块、更新每一个通道的图像vj的模块以及更新K空间数据fj的模块；其中，
更新字典D的模块临时固定X,采用梯度下降的方法求解字典D，
Figure BDA0001490895220000058
式中γ表示学习率，k表示迭代次数，
由于
Figure BDA0001490895220000059
该公式是重建目标函数对字典D的倒数，H 表示共轭转置，
因此获得下式：
Dk+1＝Dk-γ(DX-RlV)XH (3)
更新每一个通道的图像vj的模块通过下式更新vj，
Figure BDA0001490895220000061
更新K空间数据fj的模块通过下述公式更新K空间数据：
Figure BDA0001490895220000062
其中，
Figure BDA0001490895220000063
表示第j个通道的第k次迭代更新值；
并且对K空间数据进行傅里叶反变换再次更新vj，获得
Figure BDA0001490895220000064
成像模块，根据所有通道的图像vj，获得更新后的图像。
优选地，该成像装置还包括：初始化模块，用于初始化
Figure BDA0001490895220000065
优选地，构建重建模型即求解一个
Figure BDA0001490895220000067
最小目标函数，其中，
Figure BDA0001490895220000068
范数是数据拟合项，
Figure BDA0001490895220000069
范数代表稀疏表达误差，
Figure BDA00014908952200000610
混合范数代表每个通道之间的联合稀疏约束。
优选地，所述更新模块在更新字典D时，对原子进行归一化处理，以避免字典原子尺度模糊。
优选地，更新每一个通道的图像vj的模块进一步执行如下步骤：
步骤d1：根据最小二乘法获得下式：
Figure BDA0001490895220000066
式中，
Figure BDA0001490895220000071
是一个对角矩阵，对角线上元素值等于对应位置图像像素重叠的块的个数；当假设这些块映射在图像的边缘，则对角线上元素值全部相等，
Figure BDA0001490895220000072
特殊情况下，当填充的距离为单位距离1时，ω＝M；
步骤d2：图像域通过归一化全采傅里叶编码矩阵转换到傅里叶域，获得：
Figure BDA0001490895220000073
式中，
Figure BDA0001490895220000074
是0和1的对角矩阵，1代表K空间采到的点的位置，ω表示像素重复的个数，IN表示N*N的单位矩阵，
Figure BDA0001490895220000075
表示第j-th通道零填充的傅里叶测量值；
步骤d3：将步骤d2中的公式变换成：
Figure BDA0001490895220000076
式中，Fvj(kx,ky)表示K空间位置(kx,ky)的更新值，
Figure BDA0001490895220000077
表示第 j-th通道的零填充K空间测量值；η表示一个小标量，Ω表示K空间已经采到的点的集合，
Figure BDA0001490895220000078
步骤d4：通过步骤d3的公式(9)中频域中的数据通过逆傅里叶变换求得vj。
优选地，所述成像模块对所有通道的图像vj的平方进行求和，然后开方，获得更新后的图像。
本发明还提供了一种计算机可读介质，该计算机可读介质具有存储在其中的程序，该程序是计算机可执行的以使计算机执行上述的基于自适应联合稀疏编码的并行磁共振成像方法的各步骤。
有益效果
本方发明方法利用了l2,1范数开发通道的联合稀疏，开发信息稀疏性的同时，可以免校准，且该方法有较强的鲁棒性。
附图说明
图1是本发明的基于自适应联合稀疏编码的并行磁共振成像方法的流程图。
图2是本发明的基于自适应联合稀疏编码的并行磁共振成像装置的框图。
具体实施方式
在下列说明中，为了提供对本发明的彻底了解而提出许多具体细节。本发明可在不具有部分或所有这些具体细节的情况下实施。在其他情况下，为了不对本发明造成不必要的混淆，不详述众所周知的过程操作。虽然本发明将结合具体实施例来进行说明，但应当理解的是，这并非旨在将本发明限制于这些实施例。
下面根据图1对本发明的一个实施例的基于自适应联合稀疏编码的并行磁共振成像方法进行详细说明。
首先，步骤a，构建一个基于免校准的联合稀疏编码的并行磁共振成像模型即重建模型，所述模型定义为：
Figure BDA0001490895220000091
式中，V代表重建的图像，
Figure BDA0001490895220000092
代表过完备(P＞＞M)的字典，θ表示目标函数，
Figure BDA0001490895220000093
代表稀疏矩阵，
Figure BDA0001490895220000094
表示从第j-th通道的图像提取第l-th图像块，FM代表傅里叶变换，vj表示每个通道的图像，所有通道的图像vj合并在一起组成的矩阵即为V，fj代表K空间数据，
Figure BDA0001490895220000095
代表块提取矩阵，M表示原子的尺寸，N表示图像的向量化尺寸，P表示原子的个数，下标F表示傅里叶变换，λ表示数据拟合项权重，β表示联合稀疏正则化权重，||·||2，1表示联合稀疏项，
Figure BDA0001490895220000096
式中，
Figure BDA0001490895220000097
表示第j-th通道的第l-th提取块中第 p-th像素。
这里，步骤a中，构建重建模型即求解一个
Figure BDA0001490895220000099
最小目标函数，其中，
Figure BDA00014908952200000910
范数是数据拟合项，
Figure BDA00014908952200000911
范数代表稀疏表达误差，
Figure BDA00014908952200000912
混合范数代表每个通道之间的联合稀疏约束。步骤a中，还包括初始化X0,V0,fj 0的步骤。
接着进行步骤b：临时固定X,采用梯度下降的方法求解字典D，
Figure BDA0001490895220000098
式中γ表示学习率，k表示迭代次数，
由于
Figure BDA0001490895220000101
该公式是重建目标函数对字典D的倒数，H表示共轭转置，因此获得下式：
Dk+1＝Dk-γ(DX-RlV)XH (3)
这里，步骤b中，对原子进行归一化处理，以避免字典原子尺度模糊。
接着进行步骤c：更新联合稀疏系数X，直到内层迭代达到规定次数，
Figure BDA0001490895220000102
其中，q是内循环又嵌套的一层循环的次数，B表示辅助变量，α表示辅助参数。
接下来进行步骤d：通过下式更新每一个通道的图像vj，
Figure BDA0001490895220000103
这里，步骤d中，更新vj进一步包括如下步骤：
步骤d1：根据最小二乘法获得下式：
Figure BDA0001490895220000104
式中，
Figure BDA0001490895220000105
是一个对角矩阵，对角线上元素值等于对应位置图像像素重叠的块的个数；当假设这些块映射在图像的边缘，则对角线上元素值全部相等，
Figure BDA0001490895220000111
特殊情况下，当填充的距离为单位距离1时，ω＝M；
步骤d2：图像域通过归一化全采傅里叶编码矩阵转换到傅里叶域，获得：
Figure BDA0001490895220000112
式中，
Figure BDA0001490895220000113
是0和1的对角矩阵，1代表K空间采到的点的位置，ω表示像素重复的个数，IN表示N*N的单位矩阵，
Figure BDA0001490895220000114
表示第j-th通道零填充的傅里叶测量值；
步骤d3：将步骤d2中的公式变换成：
Figure BDA0001490895220000115
式中，Fvj(kx,ky)表示K空间位置(kx,ky)的更新值，
Figure BDA0001490895220000116
表示第 j-th通道的零填充K空间测量值；η表示一个小标量，Ω表示K空间已经采到的点的集合，
Figure BDA0001490895220000117
步骤d4：通过步骤d3的公式(9)中频域中的数据通过逆傅里叶变换求得vj。
接下来进行步骤e：通过下述公式更新K空间数据：
Figure BDA0001490895220000118
其中，
Figure BDA0001490895220000121
表示第j个通道的第k次迭代更新值。
接着进行步骤f：对K空间数据进行傅里叶反变换再次更新vj，获得
Figure BDA0001490895220000122
最后进行步骤g：对所有通道的图像vj的平方进行求和，然后开方，获得更新后的图像。
现有的并行磁共振成像方法中，一般利用二维随机欠采样，该采样模式会使得直接反傅里叶变换的图像有类噪声的伪影且细节丢失与模糊。而本发明提出的方法可以通过整体模型与算法实现4倍加速的高质量磁共振图像重建。具体地，联合稀疏范数可以较好地抑制噪声，字典更新可以自适应地获取图像的结构信息，而整体算法的精确迭代循环求解可以实现最小误差的重建质量。
此外，本发明还提供了一种计算机可读介质，该计算机可读介质具有存储在其中的程序，该程序是计算机可执行的，以使计算机执行包括上述各步骤的处理。
本发明还提供了一种用于上述方法的基于自适应联合稀疏编码的并行磁共振成像装置。该装置包括了不同模块用于实现上述方法中提及的各个步骤。
下面根据图2说明本发明的一个实施例的成像装置。例如，该装置可以包括模型构建模块、更新模块和成像模块。其中，模型构建模块，构建一个基于免校准的联合稀疏编码的并行磁共振成像模型即重建模型。更新模块，包括更新联合稀疏编码X的模块、更新字典D的模块、更新每一个通道的图像vj的模块以及更新K空间数据fj的模块。成像模块，根据所有通道的图像vj，获得更新后的图像。优选地，该成像装置还进一步包括：初始化模块，用于初始化
Figure BDA0001490895220000123
上述模型构建模块中，构建重建模型即求解一个
Figure BDA0001490895220000131
最小目标函数，其中，
Figure BDA0001490895220000132
范数是数据拟合项，
Figure BDA0001490895220000133
范数代表稀疏表达误差，
Figure BDA0001490895220000134
混合范数代表每个通道之间的联合稀疏约束。
优选地，上述更新模块在更新字典D时，需要对原子进行归一化处理，以避免字典原子尺度模糊。
优选地，上述成像模块对所有通道的图像vj的平方进行求和，然后开方，从而获得更新后的图像。
尽管已经根据优选的实施方案对本发明进行了说明，但是存在落入本发明范围之内的改动、置换以及各种替代等同方案。还应当注意的是，存在多种实现本发明的方法和系统的可选方式。因此，意在将随附的权利要求书解释为包含落在本发明的主旨和范围之内的所有这些改动、置换以及各种替代等同方案。

Claims (13)

1.一种基于自适应联合稀疏编码的并行磁共振成像方法，其特征在于，包括以下步骤：
步骤a：构建一个基于免校准的联合稀疏编码的并行磁共振成像模型即重建模型，所述重建模型定义为：
Figure FDA0003305101510000011
式中，V代表重建的图像，
Figure FDA0003305101510000012
代表过完备的字典，且P＞＞M，θ表示目标函数，
Figure FDA0003305101510000013
代表稀疏矩阵，
Figure FDA0003305101510000014
表示从第j-th通道的图像提取第l-th图像块，FM代表傅里叶变换，vj表示每个通道的图像，所有通道的图像vj合并在一起组成的矩阵即为V，fj代表K空间数据，
Figure FDA0003305101510000015
代表块提取矩阵，M表示原子的尺寸，N表示图像的向量化尺寸，P表示原子的个数，下标F表示傅里叶变换，λ表示数据拟合项权重，β表示联合稀疏正则化权重，||·||2,1表示联合稀疏项，
Figure FDA0003305101510000016
式中，
Figure FDA0003305101510000017
表示第j-th通道的第l-th提取块中第p-th像素；
步骤b：临时固定X,采用梯度下降的方法求解字典D，
Figure FDA0003305101510000018
式中γ表示学习率，k表示迭代次数，
由于
Figure FDA0003305101510000019
该公式是重建目标函数对字典D的导数，H表示共轭转置，
因此获得下式：
Dk+1＝Dk-γ(DX-RlV)XH (3)
步骤c：更新稀疏矩阵X；
步骤d：通过下式更新vj，
Figure FDA0003305101510000021
步骤e：通过下述公式更新K空间数据：
Figure FDA0003305101510000022
其中，
Figure FDA0003305101510000023
表示第j个通道的第k次迭代更新值；
步骤f：对K空间数据进行傅里叶反变换再次更新vj，获得
Figure FDA0003305101510000024
以及
步骤g：根据所有通道的图像vj，获得更新后的图像。
2.根据权利要求1所述的方法，其特征在于，步骤a中，构建重建模型即求解一个
Figure FDA0003305101510000025
最小目标函数，其中，
Figure FDA0003305101510000026
范数是数据拟合项，
Figure FDA0003305101510000027
范数代表稀疏表达误差，
Figure FDA0003305101510000028
混合范数代表每个通道之间的联合稀疏约束。
3.根据权利要求1所述的方法，其特征在于，步骤a中，包括初始化
Figure FDA00033051015100000211
的步骤。
4.根据权利要求1所述的方法，其特征在于，步骤b中，对原子进行归一化处理，以避免字典原子尺度模糊。
5.根据权利要求1所述的方法，其特征在于，步骤d中，更新vj进一步包括如下步骤：
步骤d1：根据最小二乘法获得下式：
Figure FDA00033051015100000210
式中，
Figure FDA0003305101510000031
是一个对角矩阵，对角线上元素值等于对应位置图像像素重叠的块的个数；当假设这些块映射在图像的边缘，则对角线上元素值全部相等，
Figure FDA0003305101510000032
特殊情况下，当填充的距离为单位距离1时，ω＝M；
步骤d2：图像域通过归一化全采傅里叶编码矩阵转换到傅里叶域，获得：
Figure FDA0003305101510000033
式中，
Figure FDA0003305101510000034
是0和1的对角矩阵，1代表K空间采到的点的位置，ω表示像素重复的个数，IN表示N*N的单位矩阵，
Figure FDA0003305101510000035
表示第j-th通道零填充的傅里叶测量值；
步骤d3：将步骤d2中的公式变换成：
Figure FDA0003305101510000036
式中，Fvj(kx,ky)表示K空间位置(kx,ky)的更新值，
Figure FDA0003305101510000037
表示第j-th通道的零填充K空间测量值；η表示一个小标量，Ω表示K空间已经采到的点的集合，
Figure FDA0003305101510000038
步骤d4：通过步骤d3的公式(9)中频域中的数据通过逆傅里叶变换求得vj。
6.根据权利要求1所述的方法，其特征在于，步骤g中，对所有通道的图像vj的平方进行求和，然后开方，获得更新后的图像。
7.一种基于自适应联合稀疏编码的并行磁共振成像装置，其特征在于，包括：
模型构建模块，构建一个基于免校准的联合稀疏编码的并行磁共振成像模型即重建模型，所述重建模型定义为：
Figure FDA0003305101510000041
式中，V代表重建的图像，
Figure FDA0003305101510000042
代表过完备的字典，且P＞＞M，θ表示目标函数，
Figure FDA0003305101510000043
代表稀疏矩阵，
Figure FDA0003305101510000044
表示从第j-th通道的图像提取第l-th图像块，FM代表傅里叶变换，vj表示每个通道的图像，所有通道的图像vj合并在一起组成的矩阵即为V，fj代表K空间数据，
Figure FDA0003305101510000045
代表块提取矩阵，M表示原子的尺寸，N表示图像的向量化尺寸，P表示原子的个数，下标F表示傅里叶变换，λ表示数据拟合项权重，β表示联合稀疏正则化权重，||·||2,1表示联合稀疏项，
Figure FDA0003305101510000046
式中，
Figure FDA0003305101510000047
表示第j-th通道的第l-th提取块中第p-th像素；
更新模块，包括：更新稀疏矩阵X的模块、更新字典D的模块、更新每一个通道的图像vj的模块以及更新K空间数据fj的模块；其中，
更新字典D的模块临时固定X,采用梯度下降的方法求解字典D，
Figure FDA0003305101510000051
式中γ表示学习率，k表示迭代次数，
由于
Figure FDA0003305101510000052
该公式是重建目标函数对字典D的导数，H表示共轭转置，
因此获得下式：
Dk+1＝Dk-γ(DX-RlV)XH (3)
更新每一个通道的图像vj的模块通过下式更新vj，
Figure FDA0003305101510000053
更新K空间数据fj的模块通过下述公式更新K空间数据：
Figure FDA0003305101510000054
其中，
Figure FDA0003305101510000055
表示第j个通道的第k次迭代更新值；
并且对K空间数据进行傅里叶反变换再次更新vj，获得
Figure FDA0003305101510000056
成像模块，根据所有通道的图像vj，获得更新后的图像。
8.根据权利要求7所述的成像装置，其特征在于，还包括：初始化模块，用于初始化
Figure FDA00033051015100000512
9.根据权利要求7所述的成像装置，其特征在于，构建重建模型即求解一个
Figure FDA0003305101510000058
最小目标函数，其中，
Figure FDA0003305101510000059
范数是数据拟合项，
Figure FDA00033051015100000510
范数代表稀疏表达误差，
Figure FDA00033051015100000511
混合范数代表每个通道之间的联合稀疏约束。
10.根据权利要求7所述的成像装置，其特征在于，所述更新模块在更新字典D时，对原子进行归一化处理，以避免字典原子尺度模糊。
11.根据权利要求7所述的成像装置，其特征在于，更新每一个通道的图像vj的模块进一步执行如下步骤：
步骤d1：根据最小二乘法获得下式：
Figure FDA0003305101510000061
式中，
Figure FDA0003305101510000062
是一个对角矩阵，对角线上元素值等于对应位置图像像素重叠的块的个数；当假设这些块映射在图像的边缘，则对角线上元素值全部相等，
Figure FDA0003305101510000063
特殊情况下，当填充的距离为单位距离1时，ω＝M；
步骤d2：图像域通过归一化全采傅里叶编码矩阵转换到傅里叶域，获得：
Figure FDA0003305101510000064
式中，
Figure FDA0003305101510000065
是0和1的对角矩阵，1代表K空间采到的点的位置，ω表示像素重复的个数，IN表示N*N的单位矩阵，
Figure FDA0003305101510000066
表示第j-th通道零填充的傅里叶测量值；
步骤d3：将步骤d2中的公式变换成：
Figure FDA0003305101510000067
式中，Fvj(kx,ky)表示K空间位置(kx,ky)的更新值，
Figure FDA0003305101510000071
表示第j-th通道的零填充K空间测量值；η表示一个小标量，Ω表示K空间已经采到的点的集合，
Figure FDA0003305101510000072
步骤d4：通过步骤d3的公式(9)中频域中的数据通过逆傅里叶变换求得vj。
12.根据权利要求7所述的成像装置，其特征在于，所述成像模块对所有通道的图像vj的平方进行求和，然后开方，获得更新后的图像。
13.一种计算机可读介质，该计算机可读介质具有存储在其中的程序，该程序是计算机可执行的以使计算机执行权利要求1-6中任一项所述的基于自适应联合稀疏编码的并行磁共振成像方法的各步骤。
CN201711246873.0A 2017-12-01 2017-12-01 一种基于自适应联合稀疏编码的并行磁共振成像方法、装置及计算机可读介质 Active CN108154484B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711246873.0A CN108154484B (zh) 2017-12-01 2017-12-01 一种基于自适应联合稀疏编码的并行磁共振成像方法、装置及计算机可读介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711246873.0A CN108154484B (zh) 2017-12-01 2017-12-01 一种基于自适应联合稀疏编码的并行磁共振成像方法、装置及计算机可读介质

Publications (2)

Publication Number Publication Date
CN108154484A CN108154484A (zh) 2018-06-12
CN108154484B true CN108154484B (zh) 2021-12-17

Family

ID=62465935

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711246873.0A Active CN108154484B (zh) 2017-12-01 2017-12-01 一种基于自适应联合稀疏编码的并行磁共振成像方法、装置及计算机可读介质

Country Status (1)

Country Link
CN (1) CN108154484B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11668777B2 (en) 2017-12-01 2023-06-06 Shenzhen Institutes Of Advanced Technology Adaptive joint sparse coding-based parallel magnetic resonance imaging method and apparatus and computer readable medium
CN110414530A (zh) * 2019-07-22 2019-11-05 上海交通大学 基于黎曼流形优化的图像去混合冲激噪声方法及系统
CN111091604B (zh) * 2019-11-18 2022-02-01 中国科学院深圳先进技术研究院 快速成像模型的训练方法、装置及服务器

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8014616B2 (en) * 2007-11-02 2011-09-06 Siemens Aktiengesellschaft System and method for fixed point continuation for total variation based compressed sensing imaging
CN102959426A (zh) * 2010-07-02 2013-03-06 皇家飞利浦电子股份有限公司 使用欠采样线圈数据来估计线圈灵敏度的并行磁共振成像
CN104360295A (zh) * 2014-11-04 2015-02-18 中国科学院深圳先进技术研究院 基于字典学习的磁共振并行成像方法和装置
CN105676156A (zh) * 2015-12-31 2016-06-15 深圳先进技术研究院 一种基于多通道协作编码的磁共振成像重建方法和装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8014616B2 (en) * 2007-11-02 2011-09-06 Siemens Aktiengesellschaft System and method for fixed point continuation for total variation based compressed sensing imaging
CN102959426A (zh) * 2010-07-02 2013-03-06 皇家飞利浦电子股份有限公司 使用欠采样线圈数据来估计线圈灵敏度的并行磁共振成像
CN104360295A (zh) * 2014-11-04 2015-02-18 中国科学院深圳先进技术研究院 基于字典学习的磁共振并行成像方法和装置
CN105676156A (zh) * 2015-12-31 2016-06-15 深圳先进技术研究院 一种基于多通道协作编码的磁共振成像重建方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MR Image Reconstruction From Highly;Saiprasad Ravishankar;《IEEE TRANSACTIONS ON MEDICAL IMAGING》;20110531;1-14 *
基于自适应稀疏表达的并行磁共振重建方法;谭莎;《集成技术》;20160531;1-6 *

Also Published As

Publication number Publication date
CN108154484A (zh) 2018-06-12

Similar Documents

Publication Publication Date Title
Tezcan et al. MR image reconstruction using deep density priors
US10902651B2 (en) Systems and methods for magnetic resonance image reconstruction
Uecker et al. ESPIRiT—an eigenvalue approach to autocalibrating parallel MRI: where SENSE meets GRAPPA
Pierre et al. Multiscale reconstruction for MR fingerprinting
CN108154484B (zh) 一种基于自适应联合稀疏编码的并行磁共振成像方法、装置及计算机可读介质
US8335955B2 (en) System and method for signal reconstruction from incomplete data
US9430854B2 (en) System and method for model consistency constrained medical image reconstruction
US20120081114A1 (en) System for Accelerated MR Image Reconstruction
US9939509B2 (en) Variable density incoherent spatiotemporal acquisition (VISTA) for highly accelerated magnetic resonance imaging
JP2007503903A (ja) 局所適応的な非線形ノイズ削減
CN110133556B (zh) 一种磁共振图像处理方法、装置、设备及存储介质
US20140340083A1 (en) Parallel acquisition image reconstruction method and device for magnetic resonance imaging
US11035919B2 (en) Image reconstruction using a colored noise model with magnetic resonance compressed sensing
US11867786B2 (en) Parameter map determination for time domain magnetic resonance
CN109633502A (zh) 磁共振快速参数成像方法及装置
CN112991483B (zh) 一种非局部低秩约束的自校准并行磁共振成像重构方法
US11776171B2 (en) Systems and methods for magnetic resonance image reconstruction
CN109343060A (zh) 基于深度学习时频分析的isar成像方法及系统
WO2019104702A1 (zh) 一种基于自适应联合稀疏编码的并行磁共振成像方法、装置及计算机可读介质
CN109471053B (zh) 一种基于双约束的介电特性迭代成像方法
CN104360295B (zh) 基于字典学习的磁共振并行成像方法和装置
Huang et al. Sparsity-promoting orthogonal dictionary updating for image reconstruction from highly undersampled magnetic resonance data
CN112581385B (zh) 基于多先验约束扩散峰度成像张量估计方法、介质和设备
CN109543717B (zh) 基于自适应邻域及字典的联合协作表达高光谱分类方法
CN115356672B (zh) 多维磁共振成像方法、系统及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant