CN106646303B - 一种欠采样磁共振波谱的快速重建方法 - Google Patents
一种欠采样磁共振波谱的快速重建方法 Download PDFInfo
- Publication number
- CN106646303B CN106646303B CN201611011513.8A CN201611011513A CN106646303B CN 106646303 B CN106646303 B CN 106646303B CN 201611011513 A CN201611011513 A CN 201611011513A CN 106646303 B CN106646303 B CN 106646303B
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- matrix
- resonance spectrum
- time
- lack sampling
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
一种欠采样磁共振波谱的快速重建方法,涉及磁共振波谱。首先对磁共振波谱的时间域信号进行欠采样,将采样数据填充在一个汉克尔矩阵当;然后利用矩阵的低秩特性重建出完整的时间域信号,再经傅里叶变换获得高质量波谱;在此过程中,利用矩阵因子分解方法,通过矩阵弗罗贝尼乌斯范数项的最小化来逼近矩阵核范数项的最小化,避免在求解核范数项时使用奇异值分解,从而达到加速重建的效果。效果优良,易于操作,可以应用在一维、二维以及更高维的欠采样的磁共振波谱重建中。
Description
技术领域
本发明涉及磁共振波谱,尤其是涉及基于时间域信号低秩汉克尔矩阵的一种欠采样磁共振波谱的快速重建方法。
背景技术
磁共振波谱被广泛地应用在临床医学和生物技术领域,其中高维磁共振波谱因能提供丰富信息而备受青睐。但对于高维核磁共振波谱来说,采样时间长是亟待解决的问题之一。用欠采样(相对于满足奈奎斯特准则下的全采样)来降低间接维的采样点数是一种缩短核磁共振采样时间的方法。但欠采样因为不满足奈奎斯特采样条件,容易造成谱峰重叠进而形成伪峰。为了获得高质量的谱,可以通过谱的自稀疏性从欠采样数据来重建这些谱(Xiaobo Qu,Xue Cao,Di Guo,Zhong Chen,"Compressed sensing for sparse magneticresonance spectroscopy,"International Society for Magnetic Resonance inMedicine 18th Scientific Meeting.Stockholm,Sweden,pp.3371,2010.)(Xiaobo Qu,DiGuo,Xue Cao,Shuhui Cai,Zhong Chen,"Reconstruction of self-sparse 2D NMRspectra from undersampled data in indirect dimension,"Sensors,vol.11,no.9,pp.8888-8909,2011.)。然而,谱峰宽度的增加将降低谱的自稀疏性,容易导致重建的宽谱峰发生畸变甚至完全丢失。屈小波等(Xiaobo Qu,Maxim Mayzel,Jian-Feng Cai,ZhongChen,Vladislav Orekhov."Accelerated NMR spectroscopy with low-rankreconstruction,"Angewandte Chemie International Edition,vol.54,no.3,pp.852-854,2015.)提出了一种基于低秩汉克尔矩阵的谱重建方法,该方法能同时重建不同宽度的谱峰,能够使用更少的欠采样数据点重建出高质量的谱。但该方法在迭代求解中最小核范数项时需要进行复杂度较高的奇异值分解,导致谱重建时间较长。
发明内容
本发明的目的在于提供效果优良、易于操作、重建时间短的欠采样磁共振波谱重建方法。
本发明包括以下步骤:
1)欠采样磁共振波谱的时间域信号,具体方法可为:
对磁共振信号进行欠采样,给定一组与预期的全采情况下信号的维度相同的且取值为0或1的数字,用于表示全采样磁共振波谱的特定时间点的信号是否采集,其中0表示不采样,1表示采样,根据0和1构成欠采样算子U;
2)构建时间域信号的汉克尔矩阵,具体方法可为:
假设x是待重建的磁共振波谱时间域信号,其维度与全采情况下数据的维度相同,然后通过线性算子R构建汉克尔矩阵,一维的磁共振波谱信号通过线性算子R构建成汉克尔矩阵,二维的磁共振波谱信号通过线性算子R构建成分块汉克尔矩阵,更高维的磁共振信号通过线性算子R构建成多维汉克尔矩阵;
3)设基于矩阵因子分解的低秩汉克尔矩阵的重建模型:
式(1)中,x为待重建磁共振波谱时间域信号,U为欠采样算子,y为欠采的磁共振波谱时间域信号,R是构建汉克尔矩阵的线性算子;||·||2为向量的二范数,λ是一个正则化参数,用于权衡和两项的重要性,符号“H”表示复共轭转置,||·||F表示矩阵的弗罗贝尼乌斯范数,P,Q是满足约束条件的任意矩阵;
4)采用基于矩阵因子分解的低秩汉克尔矩阵的重建算法:
式(1)的增广拉格朗日形式可以写成式(2):
其中D是拉格朗日乘子,其维度与Rx相同,β为大于零的参数,<·,·>表示内积;
然后通过式(3)进行迭代求解:
其中,下标k表示第k次的解,符号“-1”表示求矩阵的逆,τ是一个大于0的参数。当k=1时,Pk和Qk初始为随机矩阵,Dk是一个全为1的矩阵;当达到迭代停止准则时,迭代停止;
在步骤4)中,所述迭代停止的准则设定为达到最大迭代次数103或x在相邻两次迭代中的误差小于设定的阈值η,η设为10-6。
5)对补全的时间信号进行傅立叶变换即可得到谱图。
本发明利用弗罗贝尼乌斯范数项的最小化来逼近核范数项的最小化的性质(Srebro,Nathan."Learning with matrix factorizations,"Massachusetts Instituteof Technology,2004.),利用矩阵的因子分解来避免迭代过程中的奇异值分解,通过对磁共振波谱的时间信号进行汉克尔矩阵的低秩重建,并设计对应的迭代重建算法,实现包含多维磁共振波谱在内的欠采样信号的快速重建,进而得到高质量的磁共振谱。本发明首先对磁共振波谱的时间域信号进行欠采样,将采样数据填充在一个汉克尔矩阵当;然后利用矩阵的低秩特性重建出完整的时间域信号,再经傅里叶变换获得高质量波谱;在此过程中,利用矩阵因子分解方法,通过矩阵弗罗贝尼乌斯范数项的最小化来逼近矩阵核范数项的最小化,避免在求解核范数项时使用奇异值分解,从而达到加速重建的效果。效果优良,易于操作,可以应用在一维、二维以及更高维的欠采样的磁共振波谱重建中。
附图说明
图1是实施例中的采样模板。
图2是全采样的谱。
图3是本发明方法重建的谱。
具体实施方式
本实施例将重建二维磁共振波谱,直接维与间接维大小分别为N1=256和N2=116,包含87个谱峰。根据欠采样模板对磁共振波谱时间域信号进行欠采样,采集30%的数据,则本实施例中总的磁共振波谱数据点数为29696点,欠采样率为30%时得到的总采样数据点数为8909点。正则化参数λ=106,β=1。参见图1,本实施例的具体步骤如下:
1)欠采样磁共振波谱的时间域信号:对二维磁共振信号进行欠采样,给定一组与预期的全采情况下信号的维度相同的且取值为0或1的数字(如图2所示,黑色表示0,白色表示1),用于表示对全采样磁共振波谱的特定时间点的信号是否采集,其中0表示不采样,1表示采样。根据0和1构成欠采样算子U。
2)构建时间域信号的分块汉克尔矩阵:假设x是待重建的磁共振波谱时间域信号,其维度为256×116,然后通过线性算子R构建分块汉克尔矩阵,得到维度为7424×7611的分块汉克尔矩阵。
3)设基于矩阵因子分解的低秩汉克尔矩阵的重建模型:
式(1)中,x为待重建磁共振波谱时间域信号,U为欠采样算子,维度为256×116,y为欠采的磁共振波谱时间域信号,此实施例中为欠采样信号中丢失数据点进行填零后的欠采样数据,其维度为256×116,R是构建汉克尔矩阵的线性算子;||·||2为向量的二范数,正则化参数λ=106,符号“H”表示复共轭转置,||·||F表示矩阵的弗罗贝尼乌斯范数,P,Q是满足约束条件的矩阵,此实施例中P的维度为7424×400,Q的维度为7611×400。
4)基于矩阵因子分解的低秩汉克尔矩阵的重建算法:
式(1)的增广拉格朗日形式可以写成式(2):
其中D是拉格朗日乘子,其维度是7424×7611;β为大于零的参数,取值为1;<·,·>表示内积。
然后通过式(3)进行迭代求解如下:
其中,下标k表示第k次的解,符号“-1”表示求矩阵的逆,τ的取值为1。初值化算法(也就是k=1时)中,Pk和Qk初始为随机矩阵,Dk是一个全为1的矩阵。当达到迭代停止准则时,迭代停止。迭代停止准则设定为达到最大迭代次数103或x在相邻两次迭代中的误差小于设定的阈值η,η设为10-6。
5)数据后处理:对重建的时间信号进行傅立叶变换即可得到谱图(如图3所示)。作为参考,我们将原始的全采样时间信号做傅立叶变换得到谱图。可以看出,利用采集到的部分数据和本发明的数据补全方法,可以重建得到高质量的谱。
本发明的谱重建耗时是采用奇异值分解的谱重建(Xiaobo Qu,Maxim Mayzel,Jian-Feng Cai,Zhong Chen,Vladislav Orekhov."Accelerated NMR spectroscopy withlow-rank reconstruction,"Angewandte Chemie International Edition,vol.54,no.3,pp.852-854,2015.)耗时的7.6%,奇异值分解的谱重建时间与按照本发明方法的谱重建时间如表1所示。
表1
磁共振波谱重建方法 | 重建耗时(单位:s) |
基于奇异值分解的磁共振波谱重建 | 13758.4 |
本发明 | 1052.2 |
由表1可见,本发明可以达到快速重建波谱的目的。
Claims (4)
1.一种欠采样磁共振波谱的快速重建方法,其特征在于包括以下步骤:
1)欠采样磁共振波谱的时间域信号;
2)构建时间域信号的汉克尔矩阵;
3)设基于矩阵因子分解的低秩汉克尔矩阵的重建模型:
式(1)中,x为待重建磁共振波谱时间域信号,U为欠采样算子,y为欠采的磁共振波谱时间域信号,R是构建汉克尔矩阵的线性算子;||·||2为向量的二范数,λ是一个正则化参数,用于权衡和两项的重要性,符号“H”表示复共轭转置,||·||F表示矩阵的弗罗贝尼乌斯范数,P,Q是满足约束条件Rx=PQH的任意矩阵;
4)采用基于矩阵因子分解的低秩汉克尔矩阵的重建算法:
式(1)的增广拉格朗日形式写成式(2):
其中D是拉格朗日乘子,其维度与Rx相同,β为大于零的参数,<·,·>表示内积;
然后通过式(3)进行迭代求解:
其中,下标k表示第k次的解,符号“-1”表示求矩阵的逆,τ是一个大于0的参数;当k=1时,Pk和Qk初始为随机矩阵,Dk是一个全为1的矩阵;当达到迭代停止准则时,迭代停止;
5)对补全的时间信号进行傅立叶变换即得到谱图。
2.如权利要求1所述一种欠采样磁共振波谱的快速重建方法,其特征在于在步骤1)中,所述欠采样磁共振波谱的时间域信号的具体方法为:对磁共振信号进行欠采样,给定一组与预期的全采情况下信号的维度相同的且取值为0或1的数字,用于表示全采样磁共振波谱的特定时间点的信号是否采集,其中0表示不采样,1表示采样,根据0和1构成欠采样算子U。
3.如权利要求1所述一种欠采样磁共振波谱的快速重建方法,其特征在于在步骤2)中,所述构建时间域信号的汉克尔矩阵的具体方法为:假设x是待重建的磁共振波谱时间域信号,其维度与全采情况下数据的维度相同,然后通过线性算子R构建汉克尔矩阵,一维的磁共振波谱信号通过线性算子R构建成汉克尔矩阵,二维的磁共振波谱信号通过线性算子R构建成分块汉克尔矩阵,更高维的磁共振波谱信号通过线性算子R构建成多维汉克尔矩阵。
4.如权利要求1所述一种欠采样磁共振波谱的快速重建方法,其特征在于在步骤4)中,所述迭代停止的准则设定为达到最大迭代次数103或x在相邻两次迭代中的误差小于设定的阈值η,η设为10-6。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611011513.8A CN106646303B (zh) | 2016-11-17 | 2016-11-17 | 一种欠采样磁共振波谱的快速重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611011513.8A CN106646303B (zh) | 2016-11-17 | 2016-11-17 | 一种欠采样磁共振波谱的快速重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106646303A CN106646303A (zh) | 2017-05-10 |
CN106646303B true CN106646303B (zh) | 2018-12-18 |
Family
ID=58807928
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611011513.8A Active CN106646303B (zh) | 2016-11-17 | 2016-11-17 | 一种欠采样磁共振波谱的快速重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106646303B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108762072B (zh) * | 2018-05-21 | 2021-07-27 | 南京邮电大学 | 基于核范数子空间法和增广向量法的预测控制方法 |
CN108920423B (zh) * | 2018-05-24 | 2022-03-29 | 厦门理工学院 | 一种高保真谱重建方法 |
CN109191540B (zh) * | 2018-07-24 | 2022-06-24 | 厦门理工学院 | 一种基于截断核范数的磁共振波谱重建方法 |
CN108828482B (zh) * | 2018-08-03 | 2019-11-26 | 厦门大学 | 结合稀疏和低秩特性的欠采样磁共振扩散谱的重建方法 |
CN109165432B (zh) * | 2018-08-09 | 2022-12-13 | 厦门理工学院 | 一种基于部分奇异值和的磁共振波谱重建方法 |
CN109615675B (zh) * | 2018-12-04 | 2022-05-03 | 厦门大学 | 一种多通道磁共振成像的图像重建方法 |
CN110728624B (zh) * | 2019-09-29 | 2021-07-23 | 厦门大学 | 一种高分辨率扩散加权图像重建方法 |
CN110658484B (zh) * | 2019-10-17 | 2021-07-20 | 东北大学 | 一种磁共振波谱重建方法及系统 |
CN111324861B (zh) * | 2020-02-28 | 2022-05-03 | 厦门大学 | 一种基于矩阵分解的深度学习磁共振波谱重建方法 |
CN111538944B (zh) * | 2020-04-13 | 2022-06-24 | 厦门理工学院 | 一种基于子空间的磁共振波谱快速重建方法 |
CN111783631B (zh) * | 2020-06-29 | 2023-10-24 | 厦门大学 | 一种基于稀疏表示的深度学习磁共振波谱重建方法 |
CN111915007B (zh) * | 2020-07-29 | 2022-06-14 | 厦门大学 | 一种基于神经网络的磁共振谱降噪方法 |
CN112861074B (zh) * | 2021-03-09 | 2022-10-04 | 东北电力大学 | 基于Hankel-DMD的电力系统机电参数提取方法 |
CN114428222A (zh) * | 2022-01-28 | 2022-05-03 | 厦门大学 | 一种基于指数分解约束的磁共振谱去噪方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104739410B (zh) * | 2015-04-16 | 2017-03-15 | 厦门大学 | 一种磁共振图像的迭代重建方法 |
CN104793159B (zh) * | 2015-05-11 | 2017-06-23 | 厦门大学 | 一种高维核磁共振时域信号补全方法 |
CN105137373B (zh) * | 2015-07-23 | 2017-12-08 | 厦门大学 | 一种指数信号的去噪方法 |
CN105760682A (zh) * | 2016-02-29 | 2016-07-13 | 厦门理工学院 | 一种基于四元素汉克尔矩阵的四通道信号重建方法 |
CN105808869A (zh) * | 2016-03-16 | 2016-07-27 | 厦门理工学院 | 一种基于块Hankel矩阵的磁共振波谱重建方法 |
-
2016
- 2016-11-17 CN CN201611011513.8A patent/CN106646303B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106646303A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106646303B (zh) | 一种欠采样磁共振波谱的快速重建方法 | |
Ying et al. | Hankel matrix nuclear norm regularized tensor completion for $ n $-dimensional exponential signals | |
CN110378980A (zh) | 一种基于深度学习的多通道磁共振图像重建方法 | |
Makinen et al. | Lossless, scalable implicit likelihood inference for cosmological fields | |
CN105808869A (zh) | 一种基于块Hankel矩阵的磁共振波谱重建方法 | |
CN107423543B (zh) | 一种超复数磁共振波谱的快速重建方法 | |
CN111324861B (zh) | 一种基于矩阵分解的深度学习磁共振波谱重建方法 | |
CN104739410B (zh) | 一种磁共振图像的迭代重建方法 | |
CN111783631B (zh) | 一种基于稀疏表示的深度学习磁共振波谱重建方法 | |
Guo et al. | A fast low rank Hankel matrix factorization reconstruction method for non-uniformly sampled magnetic resonance spectroscopy | |
Wang et al. | A sparse model-inspired deep thresholding network for exponential signal reconstruction—Application in fast biological spectroscopy | |
CN103142228A (zh) | 压缩感知磁共振快速成像方法 | |
CN104013403B (zh) | 一种基于张量分解稀疏约束的三维心脏磁共振成像方法 | |
CN109191540B (zh) | 一种基于截断核范数的磁共振波谱重建方法 | |
CN109165432B (zh) | 一种基于部分奇异值和的磁共振波谱重建方法 | |
CN108828482B (zh) | 结合稀疏和低秩特性的欠采样磁共振扩散谱的重建方法 | |
CN113143243B (zh) | 一种基于谱分解的深度学习磁共振波谱重建方法 | |
Wu et al. | Fast NMR spectroscopy reconstruction with a sliding window based Hankel matrix | |
CN108537738A (zh) | 一种矩阵补全方法 | |
CN110658484B (zh) | 一种磁共振波谱重建方法及系统 | |
CN111538944B (zh) | 一种基于子空间的磁共振波谱快速重建方法 | |
CN106649201A (zh) | 一种基于指数信号的范德蒙分解的数据补全方法 | |
CN104932863A (zh) | 一种高维指数信号数据补全方法 | |
CN115170691A (zh) | 基于深度学习和交替方向乘子法的磁共振成像的重构方法 | |
CN108920423A (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 |