CN103049923A - 磁共振快速成像的方法 - Google Patents
磁共振快速成像的方法 Download PDFInfo
- Publication number
- CN103049923A CN103049923A CN201210528522XA CN201210528522A CN103049923A CN 103049923 A CN103049923 A CN 103049923A CN 201210528522X A CN201210528522X A CN 201210528522XA CN 201210528522 A CN201210528522 A CN 201210528522A CN 103049923 A CN103049923 A CN 103049923A
- Authority
- CN
- China
- Prior art keywords
- image
- gradient
- dictionary
- alpha
- updating
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000003384 imaging method Methods 0.000 title claims abstract description 26
- 238000005070 sampling Methods 0.000 claims description 10
- 102000002274 Matrix Metalloproteinases Human genes 0.000 claims description 5
- 108010000684 Matrix Metalloproteinases Proteins 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 4
- 238000006243 chemical reaction Methods 0.000 abstract description 3
- 238000012545 processing Methods 0.000 abstract description 2
- 238000007634 remodeling Methods 0.000 abstract 2
- 239000011159 matrix material Substances 0.000 description 7
- 230000009466 transformation Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明属于图像处理技术领域,尤其涉及一种磁共振快速成像的方法。本发明的共振快速成像的方法包括:步骤A:在图像梯度域的水平梯度图像和垂直梯度图像上进行字典学习,建立图像模型;步骤B:利用重建算法交替更新图像块的稀疏表示,恢复水平梯度和垂直梯度,然后在这两个方向梯度重建图像。本发明实施例磁共振快速成像的方法通过自适应学习字典的引入可以克服固定的有限差分变换导致目标图像的块状效应,可以处理结构更复杂的图像,从而产生更精确的重建;另外本发明实施例对梯度图像进行处理,其比原图像更稀疏,因而字典学习可以变得更精确和鲁棒,从而更稀疏地表示图像,且具有更好的保真度,且可以恢复出更多的细节。
Description
技术领域
本发明属于图像处理技术领域,尤其涉及一种磁共振快速成像的方法。
背景技术
为了缩短磁共振图像采集时间,压缩感知理论被成功应用到磁共振成像中。压缩感知理论利用信号在某个基的稀疏性,实现了在非相干采样矩阵下,只需少量采样(远少于奈奎斯特采样理论所需的采样)即可高质量重建原始信号。而且,若信号在某个基越稀疏,那么所需要的采样量则越少。因此,在压缩感知理论中,一个重要的问题就是稀疏基的选取。在过去几年中,从欠采样K空间信号来重建磁共振图像一般使用如下全变分模型(TV模型):
其中, 称为各向异性离散总变分正则化方程, 定义为水平方向和垂直方向的差分算子;μ>0为罚参数;第二项作为保真项,用l2-norm控制误差;我们定义为需要重建的图像,表示傅里叶欠采样信号,两者的关系为Fpu=f,其中表示傅里叶欠采样矩阵。
TV模型有着很好的保存图像边缘的能力。但是在较大的欠采样下可能会产生块状效应,因此把该模型应用于磁共振快速成像上时,图像重建质量会受到一定的影响。为了提高图像重建质量,可以在TV模型基础上添加其它稀疏约束项(例如小波变换),例如近年来由Yang等人提出的RecPF快速重建方法模型如下:
其中,ψ表示小波变换,μ1,μ2>0,用于权衡前两个正则化项和保真项。
除了引入基于已知固定的稀疏变换的正则化来提高稀疏性,近年来,基于图像块字典学习的图像稀疏表示方法得到了越来越多的研究。其典型的方法描述如下:
对一幅的图像u,用一个图像块提取算子R对图像进行分块,得到包含L个图像块的集合R(u)=[R1u,R2u,…,RLu],定义为大小为的图像块的向量形式。Elad等人提出,对于所有的图像块Rlu,都可以在字典D*上稀疏的表示,其模型如下:
其中,αl表示第l个图像块在字典D*上的表示系数,T0控制表示系数的稀疏度,||·||0定义为向量中非零元素的个数。
Ravishanker等人把字典学习模型应用到K空间欠采样的磁共振图像重建上,提出了DLMRI模型:
其中,Γ=[α1,α2,…,αL]定义为所有图像块对应的稀疏系数矩阵。前一项保证图像块在自适应的学习字典上稀疏表示,后一项为在K空间上的信号保真项,而正则化参数v平衡这两项间的权重。求解该模型一般用两步迭代交替更新的方法:第一步,应用K-SVD算法训练稀疏表示字典来消除混淆和噪声;第二步,通过学习字典和稀疏系数来重建图像。
现有技术的缺陷在于:在TV模型和RecPF模型中,TV正则化及小波稀疏约束等作为非自适应的固定变换,并不能理想地稀疏表示所有的图像;而Elad等人和Ravishanker等人提出的字典学习方法都是在图像域上训练字典,不能使图像更稀疏的表示,导致需要采样的信号多,另外上述的方法在图像域上学习反应细节较差。
发明内容
本发明提供了一种磁共振快速成像的方法,旨在解决现有的磁共振成像方法不能理想地稀疏表示所有的图像,导致需要采样信号多且成像反应细节差的技术问题。
本发明提供的技术方案为:一种磁共振快速成像的方法,包括:
步骤A:在图像梯度域的水平梯度图像和垂直梯度图像上进行字典学习,建立图像模型;
步骤B:利用重建算法交替更新图像块的稀疏表示,恢复水平梯度和垂直梯度,然后在这两个方向梯度重建图像。
本发明的技术方案还包括:所述步骤A还包括:在图像域上进行字典学习,即在图像域和梯度域的水平梯度图像和垂直梯度图像上同时进行字典学习,建立图像模型。
本发明的技术方案还包括:在所述步骤B中,恢复图像及其水平、垂直梯度图像,然后根据图像以及水平、垂直方向梯度重建图像。
本发明的技术方案还包括:在所述步骤A中,建立的图像模型为:
本发明的技术方案还包括:所述步骤B的重建算法包括:引入辅助变量w(i),i=1,2,应用布雷格曼技术,定义
本发明的技术方案还包括:所述步骤B的重建算法还包括:更新梯度图像变量w(i),i=1,2;更新梯度图像块的稀疏表示;交替更新字典D(i)和系数矩阵在稀疏编码阶段,固定字典D(i),通过贪婪算法正交匹配追踪来更新在字典更新阶段,固定系数通过奇异值分解来逐列更新字典的每一列,最小化近似误差。
本发明的技术方案还包括:所述建立的图像模型为:
本发明的技术方案还包括:所述重建算法包括:引入辅助变量v,应用布雷格曼技术,将图像模型变更为:
本发明的技术方案还包括:所述重建算法还包括:更新u,固定其余的变量;更新变量v;更新字典及系数:B和xl,D和αl,l=1,2,…,L,交替更新字典D(i)和系数矩阵在稀疏编码阶段,固定字典D(i),通过贪婪算法正交匹配追踪来更新在字典更新阶段,则固定系数通过奇异值分解来逐列更新字典的每一列,最小化近似误差。
本发明的技术方案具有如下优点或有益效果:本发明实施例磁共振快速成像的方法通过自适应学习字典的引入可以克服固定的有限差分变换导致目标图像的块状效应,可以处理结构更复杂的图像,从而产生更精确的重建;另外本发明实施例对梯度图像进行处理,其比原图像更稀疏,因而字典学习可以变得更精确和鲁棒,从而更稀疏地表示图像,且具有更好的保真度,且可以恢复出更多的细节。
附图说明
附图1是本发明第一实施例的磁共振快速成像的方法的流程图;
附图2是本发明第二实施例的磁共振快速成像的方法的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
请参阅图1,为本发明第一实施例的磁共振快速成像的方法的流程图。本发明第一实施例的磁共振快速成像的方法包括:
步骤100:在图像梯度域的水平梯度图像和垂直梯度图像上进行字典学习,建立图像模型;
其中,自适应学习字典的引入可以克服有限差分变换导致重建图像的块状效应,同时也可以更好地表示结构复杂的图像,从而产生更精确的重构,在稀疏的变换域上学习字典,其稀疏表示效果越好,重建的效果也因此更好。对于图像的一阶梯度,其水平梯度图像和垂直梯度图像比原图像更稀疏。因此,为了同时利用字典学习的优点和磁共振图像在梯度域上的稀疏性质,本发明第一实施例的磁共振快速成像在图像梯度域(包含两个方向,水平方向和垂直方向)上进行字典学习的方法,建立如下新模型(称为TVDL模型):
步骤110:利用重建算法交替更新图像块的稀疏表示,恢复水平梯度和垂直梯度,然后在这两个方向梯度重建图像。
对于本发明第一实施例提出的模型(5),可以应用Bregman技术来求解,重建算法称为GradDLRec,该算法交替更新图像块的稀疏表示,恢复水平梯度和垂直梯度,然后由这两个方向梯度重建图像。算法的具体描述如下:
引入辅助变量w(i),i=1,2,方程(5)可以改写如下:
通过应用Bregman方法,定义 可以得到如下子问题:
其中v2定义为正罚参数。应用交替方向法(ADM,Alternating DirectionMethod)来求解问题(7)。问题(7)求解过程如下:
由于(9)是一个简单的最小二乘问题,可以直接求出变量u的解析解:
步骤72:更新梯度图像变量w(i),i=1,2(分别对应水平、垂直方向);
由于w(1)和w(2)是分离的,因此可以分别求解:
该最小二乘问题的解析解如下:
由于问题(7)中关于梯度图像在水平和垂直方向上的字典和系数是分离的,因此可以分别求解两个方向上对应的字典及系数:
求解(13)的方法与模型(3)使用的K-SVD方法一样,即交替更新字典D(i)和系数矩阵在稀疏编码阶段,固定字典D(i),通过贪婪算法正交匹配追踪(OMP)来更新在字典更新阶段,则固定系数通过奇异值分解(SVD)来逐列更新字典的每一列,从而最小化近似误差。
以上即为求解问题(7)的过程。综上所述,本发明实施例提出的完整的GradDLRec算法可归纳如下:
算法1:GradDLRec算法,包括:
1:初始化:
2:For k=1,2,…直到满足终止准则:
3:
4:更新
5:
6:
7:End
8:输出uk+1
请参阅图2,为本发明第二实施例的磁共振快速成像的方法的流程图。本发明第二实施例的磁共振快速成像的方法包括:
步骤200:在图像域和梯度域的水平梯度图像和垂直梯度图像上同时进行字典学习,建立图像模型;
在步骤200中,在图像及其水平和垂直梯度图像上同时进行字典学习,建立如下图像模型:
由于 因此提出如下新模型(称为HTVDL模型):
其中,D,B=[B1,B2]分别为原图像块和梯度图像块对应的字典, 为相应的系数;τ,μ1为权重参数。
步骤210:利用重建算法交替更新图像块的稀疏表示,恢复图像及其水平、垂直梯度图像,然后根据图像以及水平、垂直方向梯度重建图像。
引入辅助变量v,问题(15)可以改写成如下形式:
应用Bregman技术,问题(16)可以通过如下方法求解:
其中,μ2为正罚参数。
问题(17)的求解过程如下:
步骤171:更新u,固定其余的变量,问题(17)关于变量u的方程如下:
该最小二乘问题的解满足如下等式:
令定义为正规化的全采样傅里叶矩阵,即FTF=1N。对(20)左右两边做傅里叶变换:
应用傅里叶变换的卷积定理,可得:
步骤172:更新变量v;
问题(17)关于变量v的方程如下:
解为:
步骤173:更新字典及系数:B和xl,D和αl,l=1,2,…,L
求解(25)(26)的方法与K-SVD、DLMRI一致。
以上即为求解问题(17)的过程。
综上所述,本发明实施例提出的完整的HDLRec算法可归纳如下:
算法2:HDLRec算法,包括:
1:初始化:
2:For k=1,2,…直到满足终止准则:
3:更新
5:更新
6:
7:
8:End
9:输出uk+1
本发明实施例磁共振快速成像的方法提出两种新的自适应字典学习模型,第一种为稀疏梯度域自适应字典学习模型(简称TVDL模型),该模型在图像的水平梯度和垂直梯度上分别进行字典学习,可以看做TV模型的自适应策略的推广。第二种模型为图像域及其一阶梯度域自适应字典学习模型(简称HTVDL模型),该模型在图像域及其一阶梯度域的两个方向上同时进行字典学习,可以看做是图像域字典学习的高阶模型,该模型内含了TVDL模型,因而保存了TVDL模型的优点。
本发明实施例磁共振快速成像的方法通过自适应学习字典的引入可以克服固定的有限差分变换导致目标图像的块状效应,可以处理结构更复杂的图像,从而产生更精确的重建;另外本发明实施例对梯度图像进行处理,其比原图像更稀疏,因而字典学习可以变得更精确和鲁棒,从而更稀疏地表示图像,且具有更好的保真度,且可以恢复出更多的细节。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种磁共振快速成像的方法,包括:
步骤A:在图像梯度域的水平梯度图像和垂直梯度图像上进行字典学习,建立图像模型;
步骤B:利用重建算法交替更新图像块的稀疏表示,恢复水平梯度和垂直梯度,然后在这两个方向梯度重建图像。
2.根据权利要求1所述的磁共振快速成像的方法,其特征在于,所述步骤A还包括:在图像域上进行字典学习,即在图像域和梯度域的水平梯度图像和垂直梯度图像上同时进行字典学习,建立图像模型。
3.根据权利要求2所述的磁共振快速成像的方法,其特征在于,在所述步骤B中,恢复图像及其水平、垂直梯度图像,然后根据图像以及水平、垂直方向梯度重建图像。
4.根据权利要求1所述的磁共振快速成像的方法,其特征在于,在所述步骤A中,建立的图像模型为:
其中,第一项为梯度图像在字典上稀疏表示,第二项保证重建结果与K空间采样信号相匹配;权重v1=(λ/σ),σ为测量噪声的标准差,λ为正常数。
5.根据权利要求1所述的磁共振快速成像的方法,其特征在于,所述步骤B的重建算法包括:引入辅助变量w(i),i=1,2,应用布雷格曼技术,定义
7.根据权利要求2所述的磁共振快速成像的方法,其特征在于,所述建立的图像模型为: 其中,
D,B=[B1,B2]分别为原图像块和梯度图像块对应的字典, 为相应的系数;τ,μ1为权重参数。
8.根据权利要求7所述的磁共振快速成像的方法,其特征在于,所述重建算法包括:引入辅助变量v,应用布雷格曼技术,将图像模型变更为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210528522.XA CN103049923B (zh) | 2012-12-10 | 2012-12-10 | 磁共振快速成像的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210528522.XA CN103049923B (zh) | 2012-12-10 | 2012-12-10 | 磁共振快速成像的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103049923A true CN103049923A (zh) | 2013-04-17 |
CN103049923B CN103049923B (zh) | 2015-12-09 |
Family
ID=48062551
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210528522.XA Active CN103049923B (zh) | 2012-12-10 | 2012-12-10 | 磁共振快速成像的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103049923B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103472419A (zh) * | 2013-08-30 | 2013-12-25 | 深圳先进技术研究院 | 磁共振快速成像方法及其系统 |
CN104376198A (zh) * | 2014-11-04 | 2015-02-25 | 中国科学院深圳先进技术研究院 | 自适应磁共振并行成像方法和装置 |
CN104899906A (zh) * | 2015-06-12 | 2015-09-09 | 南方医科大学 | 基于自适应正交基的磁共振图像重建方法 |
CN105654527A (zh) * | 2015-12-31 | 2016-06-08 | 中国科学院深圳先进技术研究院 | 一种基于结构化字典学习的磁共振成像重建方法和装置 |
CN106056554A (zh) * | 2016-06-01 | 2016-10-26 | 南昌大学 | 一种梯度域卷积稀疏编码的磁共振快速成像方法 |
CN106537168A (zh) * | 2014-02-11 | 2017-03-22 | 通用医疗公司 | 用于磁共振指纹识别中的自适应字典匹配的系统和方法 |
WO2019104702A1 (zh) * | 2017-12-01 | 2019-06-06 | 深圳先进技术研究院 | 一种基于自适应联合稀疏编码的并行磁共振成像方法、装置及计算机可读介质 |
CN110678899A (zh) * | 2017-06-02 | 2020-01-10 | 三菱电机株式会社 | 成像系统和图像处理方法 |
WO2020019684A1 (zh) * | 2018-07-27 | 2020-01-30 | 深圳先进技术研究院 | 低剂量pet图像重建方法、装置、设备及存储介质 |
CN115032453A (zh) * | 2022-08-08 | 2022-09-09 | 四川大学 | 一种多频动态相量测量方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106056644B (zh) * | 2016-05-24 | 2020-03-17 | 深圳先进技术研究院 | Ct扫描的数据处理方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090141995A1 (en) * | 2007-11-02 | 2009-06-04 | Siemens Corporate Research, Inc. | System and Method for Fixed Point Continuation for Total Variation Based Compressed Sensing Imaging |
CN102540116A (zh) * | 2011-12-08 | 2012-07-04 | 中国科学院深圳先进技术研究院 | 磁共振成像方法和系统 |
CN102708576A (zh) * | 2012-05-18 | 2012-10-03 | 西安电子科技大学 | 基于结构字典的分块图像压缩感知重建方法 |
-
2012
- 2012-12-10 CN CN201210528522.XA patent/CN103049923B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090141995A1 (en) * | 2007-11-02 | 2009-06-04 | Siemens Corporate Research, Inc. | System and Method for Fixed Point Continuation for Total Variation Based Compressed Sensing Imaging |
CN102540116A (zh) * | 2011-12-08 | 2012-07-04 | 中国科学院深圳先进技术研究院 | 磁共振成像方法和系统 |
CN102708576A (zh) * | 2012-05-18 | 2012-10-03 | 西安电子科技大学 | 基于结构字典的分块图像压缩感知重建方法 |
Non-Patent Citations (2)
Title |
---|
SAIPRASAD RAVISHANKAR ET AL: "MR Image Reconstruction From Highly Undersampled k-Space Data by Dictionary Learning", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》 * |
VISHAL M. PATEL ET AL: "Gradient-Based Image Recovery Methods From Incomplete Fourier Measurements", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103472419A (zh) * | 2013-08-30 | 2013-12-25 | 深圳先进技术研究院 | 磁共振快速成像方法及其系统 |
CN106537168A (zh) * | 2014-02-11 | 2017-03-22 | 通用医疗公司 | 用于磁共振指纹识别中的自适应字典匹配的系统和方法 |
CN104376198B (zh) * | 2014-11-04 | 2017-06-09 | 中国科学院深圳先进技术研究院 | 自适应磁共振并行成像方法和装置 |
CN104376198A (zh) * | 2014-11-04 | 2015-02-25 | 中国科学院深圳先进技术研究院 | 自适应磁共振并行成像方法和装置 |
CN104899906A (zh) * | 2015-06-12 | 2015-09-09 | 南方医科大学 | 基于自适应正交基的磁共振图像重建方法 |
CN104899906B (zh) * | 2015-06-12 | 2019-02-12 | 南方医科大学 | 基于自适应正交基的磁共振图像重建方法 |
CN105654527A (zh) * | 2015-12-31 | 2016-06-08 | 中国科学院深圳先进技术研究院 | 一种基于结构化字典学习的磁共振成像重建方法和装置 |
CN106056554A (zh) * | 2016-06-01 | 2016-10-26 | 南昌大学 | 一种梯度域卷积稀疏编码的磁共振快速成像方法 |
CN106056554B (zh) * | 2016-06-01 | 2019-04-12 | 南昌大学 | 一种梯度域卷积稀疏编码的磁共振快速成像方法 |
CN110678899A (zh) * | 2017-06-02 | 2020-01-10 | 三菱电机株式会社 | 成像系统和图像处理方法 |
WO2019104702A1 (zh) * | 2017-12-01 | 2019-06-06 | 深圳先进技术研究院 | 一种基于自适应联合稀疏编码的并行磁共振成像方法、装置及计算机可读介质 |
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 |
WO2020019684A1 (zh) * | 2018-07-27 | 2020-01-30 | 深圳先进技术研究院 | 低剂量pet图像重建方法、装置、设备及存储介质 |
CN115032453A (zh) * | 2022-08-08 | 2022-09-09 | 四川大学 | 一种多频动态相量测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103049923B (zh) | 2015-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103049923B (zh) | 磁共振快速成像的方法 | |
CN103400402B (zh) | 基于低秩结构稀疏的压缩感知mri图像重建方法 | |
CN103472419B (zh) | 磁共振快速成像方法及其系统 | |
CN107064845A (zh) | 基于深度卷积网的一维部分傅里叶并行磁共振成像方法 | |
CN106373167A (zh) | 一种基于深度神经网络的压缩传感核磁共振成像方法 | |
CN104899906B (zh) | 基于自适应正交基的磁共振图像重建方法 | |
CN104751420B (zh) | 一种基于稀疏表示和多目标优化的盲复原方法 | |
CN105513026A (zh) | 一种基于图像非局部相似的压缩感知重构方法 | |
CN104933683A (zh) | 一种用于磁共振快速成像的非凸低秩重建方法 | |
CN103077544B (zh) | 一种磁共振参数匹配方法、装置及医疗图像处理设备 | |
CN109658468B (zh) | 磁共振参数成像方法、装置、设备及存储介质 | |
CN102217934A (zh) | 磁共振成像方法及系统 | |
CN112991483B (zh) | 一种非局部低秩约束的自校准并行磁共振成像重构方法 | |
CN110675318B (zh) | 一种基于主结构分离的稀疏表示图像超分辨率重建方法 | |
US20200072933A1 (en) | Method for removing ghost artifact of echo planar imaging by using neural network and apparatus therefor | |
CN104739410A (zh) | 一种磁共振图像的迭代重建方法 | |
CN104574456A (zh) | 一种基于图正则化稀疏编码的磁共振超欠采样k数据成像方法 | |
CN105118078A (zh) | 欠采样的ct图像重建方法 | |
CN107767337A (zh) | 一种图像超分辨率重建方法、装置、移动终端及存储介质 | |
CN104915935B (zh) | 基于非线性压缩感知与字典学习的压缩光谱成像方法 | |
CN102163338B (zh) | 一种压缩感知系统中的高效重建方法 | |
CN108122262B (zh) | 基于主结构分离的稀疏表示单帧图像超分辨率重建算法 | |
Huang et al. | Sparsity-promoting orthogonal dictionary updating for image reconstruction from highly undersampled magnetic resonance data | |
CN104360295B (zh) | 基于字典学习的磁共振并行成像方法和装置 | |
CN109934884A (zh) | 一种基于变换学习和联合稀疏性的迭代自一致性并行成像重构方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |