CN101975935A - 一种基于部分回波压缩感知的快速磁共振成像方法 - Google Patents
一种基于部分回波压缩感知的快速磁共振成像方法 Download PDFInfo
- Publication number
- CN101975935A CN101975935A CN2010102720899A CN201010272089A CN101975935A CN 101975935 A CN101975935 A CN 101975935A CN 2010102720899 A CN2010102720899 A CN 2010102720899A CN 201010272089 A CN201010272089 A CN 201010272089A CN 101975935 A CN101975935 A CN 101975935A
- Authority
- CN
- China
- Prior art keywords
- data
- signal
- image space
- random mask
- space
- 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
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于部分回波压缩感知的快速磁共振成像方法。现有的成像方法速度比较慢,硬件成本高。本发明方法首先采集随机变密度部分回波数据;具体方法为在k-空间中心区域高密度采集数据,四周随机稀疏采集,生成二维随机掩模;将频率编码轴上需采集的数据点处都添加二维随机掩模,形成三维随机掩模;按照生成的三维随机掩模采集k-空间数据;然后基于小波域软域值去噪的凸集投影方法重建;最后基于有限差分变换的最小化L1范数的非线性重建,具体方法对图像空间信号x进行稀疏变换,确立优化目标,求解优化目标。本发明方法将部分回波和压缩感知技术结合起来应用在MRI的数据采集中,在减小回波时间的同时,缩减了数据采集时间。
Description
技术领域
本发明属于磁共振的图像成像领域,涉及一种基于部分回波压缩感知的快速磁共振成像方法。
背景技术
磁共振成像MRI(Magnetic Resonance Imaging)是一种利用体外测量的核磁共振信号产生体内器官物理和化学特性的断层图像成像技术。由于MRI对机体没有不良影响,具有较高的软组织分辨能力等优点,目前已在临床疾病检测中得到广泛应用。但是,MRI在应用中还经常受到数据采集时间过长的限制。
MRI的数据采集是在k-空间(频率空间)进行的。k空间编码由kz方向的层编码、ky方向的相位编码和kx方向的频率编码组成。在相位和层编码平面(ky_kz)上减少数据采集点数,可以有效地缩短数据采集时间。在磁共振血管成像MRA(Magnetic Resonance Angiography)中常用部分回波(Partial Echo)采集来缩短回波时间,减小高阶流动伪影。部分回波采集,就是只采集频率编码kx方向的一部分数据。
近年来,压缩感知(Compressed Sensing)技术突破了香浓(Shannon)采样定理关于采样速率必须高于2倍信号带宽的极限,从而在信号处理等多个领域成为新的研究热点(Candès et al.,Robust Uncertainty Principles:ExactSignal Reconstruction from Highly Incomplete Frequency Information,2004)。压缩感知技术利用信号自身或其在变换域中的稀疏性,只需随机采集少量的数据点,就可通过非线性重建等算法恢复出原始的信号来。然而,MRI受硬件扫描系统有规律往返轨迹的限制,不能做到完全的随机采集。
2007年,M.Lustig等人提出了适合MRI硬件扫描系统的伪随机数据采集方式,首次将压缩感知技术成功应用于MRI领域(M.Lustig,et.al.,SparseMRI:The Application of Compressed Sensing for Rapid MR Imaging,MagneticResonance in Medicine,58:1182-1195,2007)。他们利用变密度随机采集的方式采集数据,采用小波变换等方式作为图像的稀疏变换,用基于L1范数最小化的非线性优化方法重建图像。此后,有人将压缩感知技术应用到MRI的多个应用领域。
目前已申请的有关压缩感知技术方面的发明专利有:宽带移动通信中利用压缩感知减少导频数的信道估计方法(200910079441.4),提出一种利用压缩感知技术的原理,降低系统估计信道时所需的导频符号个数和保证信道估计性能的方法。一种稀疏信号的重构方法(200910023785.3)提出一种稀疏信号的重构方法,主要解决从观测向量重构原始稀疏信号速率低的问题。基于压缩感知技术的分布式信源编码的方法(200910242622.4),提出利用压缩感知技术的优点和视频图像的稀疏特性,形成一种新的分布式信源编码方法。基于正交匹配追踪的具有压缩感知过程的信号检测方法(201010032485.4),提出一种基于正交匹配追踪的具有压缩感知过程的信号检测方法,解决了利用现有的压缩感知实现信号检测时需要重构信号所带来的资源浪费的问题。传感器网络物理信号协同压缩感知系统及感知方法(200910198056.1),提出一种传感器网络物理信号协同压缩感知系统及感知方法。以上专利主要应用于宽带移动、无线电等信号检测和处理领域。
以上申请的发明专利表明,还没有一种专门针对MRI领域的压缩感知技术方面的专利。
发明内容
部分回波技术是有规律的采集k-空间,而压缩感知技术则是随机的采集,本发明针对现有技术的不足,将部分回波和压缩感知技术结合起来应用在MRA的数据采集中,在缩短回波时间,抑制高阶流动伪影的同时,缩减数据采集的时间,提供一种适用于硬件的随机变密度部分回波数据采集方法,并提供一种基于小波域软域值去噪的凸集投影方法POCS(Projection onto Convex Sets)重建和基于有限差分变换的最小化L1范数的非线性重建相结合的方法来重建图像。
本发明主要包括三大步骤:随机变密度部分回波数据采集、基于小波软域值去噪的POCS重建、基于有限差分变换的最小化L1范数的非线性重建。
随机变密度部分回波数据采集包括产生三维随机掩模、根据三维随机掩模采集数据两个步骤:
1-1、产生三维随机掩模
在相位编码和层编码平面产生二维的随机掩模,具体是:在信息量集中的k-空间中心区域高密度采集数据,在中心区域四周进行随机稀疏采集,生成二维随机掩模;
频率编码轴上需采集的数据点处都添加二维随机掩模,由此形成基于多层二维随机掩模的三维随机掩模。需采集的数据点在三维随机掩模中以“1”表示,不采集的数据点以“0”表示。
1-2、根据三维随机掩模采集数据
按照步骤1-1生成的三维随机掩模采集k-空间数据。
基于小波域软域值去噪的POCS重建包括初始化、利用低频信息产生相位约束、傅立叶反变换至图像空间、离散小波变换、小波域的软域值法去噪、离散小波反变换、图像空间数据合成、经过傅立叶变换转换到k-空间和k-空间数据合成:
2-1、初始化
初始化序号i和k-空间数据yi,使得i=0,yi=y;其中y代表采集到的欠采样数据;
2-2、利用低频信息产生相位约束
取步骤1-1的k-空间数据yi的中心低频部分,做傅立叶反变换得到图像空间信号,取出该图像空间信号中的相位信号φ作为三维POCS迭代的约束条件。
2-3、傅立叶反变换至图像空间
对步骤2-1的yi做傅立叶反变换到图像空间,形成图像空间信号,以xi表示。
2-4、离散小波变换
用正交小波对步骤2-3的xi进行离散小波变换。
2-5、小波域的软域值法去噪
用软域值(Soft Thresholding)的方法对步骤2-4生成的小波系数进行去噪处理。
其中,ψ表示步骤2-4中的小波变换,λ为域值,median表示中值滤波,α是常数,xres×yres为图像的像素点数。
2-6、离散小波反变换
将步骤2-5中去噪后的小波系数经过离散小波反变换至图像空间。
2-7、图像空间数据合成
取出步骤2-6产生的图像空间信号xi的幅值信号|xi|,与步骤2-2产生的相位信号φ组合成新的图像空间数据xi′。
xi ′=|xi|eiφ (公式3)
2-8、经过傅立叶变换转换到k-空间
将步骤2-7产生的xi′数据,经傅立叶变换到k-空间形成数据wi。
2-9、k-空间数据合成
用实际采集到的y信号替换步骤(8)产生的wi中相应位置的信号,得到信号yi′。
2-10、重复迭代
将步骤2-9中生成的yi′信号作为yi送至步骤2-3,重复步骤2-3至步骤2-10n次,其中n为自然数且n≥4。
2-11、将yi作傅立叶反变换至图像空间形成图像空间信号x。
基于有限差分变换的最小化L1范数的非线性重建包括有限差分变换(Finite Difference)、确立最小化L1范数优化目标和优化目标的求解三个步骤。
3-1、有限差分变换
采用有限差分变换对图像空间信号x进行稀疏变换。
3-2、确立最小化L1范数优化目标
假设Fu代表欠采样数据的傅立叶变换,非线性凸优化重建需解决以下问题:
其中,需重建的图像由优化变量x表示,ψ为有限差分变换,优化目标为x在有限差分变换域的L1范数最小化,域值参数ε是期望的噪声水平,用以控制重建的保真度。
3-3、优化目标的求解
采用复共轭梯度法进行优化目标的求解,得到磁共振图像。
采用本发明方法可以将部分回波和压缩感知技术结合起来应用在MRI的数据采集中,在减小回波时间的同时,缩减数据采集时间,同时本发明具有以下特点:
(1)将压缩感知技术中的随机数据采集与有规律的部分回波采集相结合,解决完全随机采样不适用于实际硬件的矛盾。通过频率编码方向的欠采样缩短回波时间,减小高阶的流动伪影,通过相位编码与层编码平面的随机欠采样,缩短数据采集时间。
(2)用已采数据和已采数据的低频相位约束,通过POCS算法在图像域和k-空间域反复迭代,可以恢复大部分未采集的k-空间数据。
(3)利用小波域软域值去噪的方法抑制随机采样带来的变换域中的噪声状伪影。
(4)针对MRA图像在有限差分域中具有较高的稀疏表现,采用基于有限差分变换的最小化L1范数的非线性重建进一步优化MRA图像重建质量。
附图说明
图1是相位编码和层编码平面的二维随机掩模示意图;
图2是三维随机掩模在相位编码和频率编码截面上的示意图;
图3是POCS算法的二维图解;
图4(a)是常规全采样参考数据经过傅立叶重建的图像;
图4(b)是欠采样MRA数据的傅立叶重建结果图;
图4(c)是欠采样MRA数据用本发明方法的重建结果图。
具体实施方式
以下结合附图对本发明作进一步说明。
一种基于部分回波压缩感知的快速磁共振成像方法包括三大步骤:随机变密度部分回波数据采集、基于小波软域值去噪的POCS重建、基于有限差分变换的最小化L1范数的非线性重建。
随机变密度部分回波数据采集包括产生三维随机掩模、根据三维随机掩模采集数据两个步骤:
1-1、产生三维随机掩模
在相位编码和层编码平面产生二维的随机掩模,具体是:在信息量集中的k-空间中心区域高密度采集数据,在中心区域四周进行随机稀疏采集,生成二维随机掩模;
频率编码轴上需采集的数据点处都添加二维随机掩模,由此形成基于多层二维随机掩模的三维随机掩模。需采集的数据点在三维随机掩模中以“1”表示,不采集的数据点以“0”表示。
1-2、根据三维随机掩模采集数据
按照步骤1-1生成的三维随机掩模采集k-空间数据。
基于小波域软域值去噪的POCS重建包括初始化、利用低频信息产生相位约束、傅立叶反变换至图像空间、离散小波变换、小波域的软域值法去噪、离散小波反变换、图像空间数据合成、经过傅立叶变换转换到k-空间和k-空间数据合成:
2-1、初始化
初始化序号i和k-空间数据yi,使得i=0,yi=y;其中y代表采集到的欠采样数据;
2-2、利用低频信息产生相位约束
取步骤1-1的k-空间数据yi的中心低频部分,做傅立叶反变换得到图像空间信号,取出该图像空间信号中的相位信号φ作为三维POCS迭代的约束条件。
2-3、傅立叶反变换至图像空间
对步骤2-1的yi做傅立叶反变换到图像空间,形成图像空间信号,以xi表示。
2-4、离散小波变换
用正交小波对步骤2-3的xi进行离散小波变换。
2-5、小波域的软域值法去噪
用软域值(Soft Thresholding)的方法对步骤2-4生成的小波系数进行去噪处理。
其中,ψ表示步骤2-4中的小波变换,λ为域值,median表示中值滤波,α是常数,xres×yres为图像的像素点数。
2-6、离散小波反变换
将步骤2-5中去噪后的小波系数经过离散小波反变换至图像空间。
2-7、图像空间数据合成
取出步骤2-6产生的图像空间信号xi的幅值信号|xi|,与步骤2-2产生的相位信号φ组合成新的图像空间数据xi′。
xi′=|xi|eiφ (公式3)
2-8、经过傅立叶变换转换到k-空间
将步骤2-7产生的xi′数据,经傅立叶变换到k-空间形成数据wi。
2-9、k-空间数据合成
用实际采集到的y信号替换步骤(8)产生的wi中相应位置的信号,得到信号yi′。
2-10、重复迭代
将步骤2-9中生成的yi′信号作为yi送至步骤2-3,重复步骤2-3至步骤2-10n次,其中n为自然数且n≥4。
2-11、将步骤2-10产生的yi′(此处的yi′=yi)信号作傅立叶反变换至图像空间形成图像空间信号x。
基于有限差分变换的最小化L1范数的非线性重建包括有限差分变换(Finite Difference)、确立最小化L1范数优化目标和优化目标的求解三个步骤。
3-1、有限差分变换
采用有限差分变换对图像空间信号x进行稀疏变换。
3-2、确立最小化L1范数优化目标
假设Fu代表欠采样数据的傅立叶变换,非线性凸优化重建需解决以下问题:
其中,需重建的图像由优化变量x表示,ψ为有限差分变换,优化目标为x在有限差分变换域的L1范数最小化,域值参数ε是期望的噪声水平,用以控制重建的保真度。
3-3、优化目标的求解
采用复共轭梯度法进行优化目标的求解。
以下结合MRA图像的部分回波压缩感知欠采集和重建进行实例说明。假设要采集的MRA图像的三维矩阵大小为kx×ky×kz=512×384×64。首先,产生部分回波随机变密度三维掩模,如图1所示,三维掩模在相位编码与层编码平面二维的随机变密度欠采样率为70%,即该平面只采集384×64×30%=7372个点。如图2所示,在频率编码方向部分回波的欠采样率为30%,只采集362个点。欠采样总量高达79%。
用三维梯度回波序列按三维掩模采集MRA数据。采集的数据先进行三维POCS重建,图3为三维POCS重建的二维示意图。如图3所示:初始化时,i=0,令yi=y;取yi的低频信号,在其周围填零后,做傅立叶反变换,取出变换后的相位信号φ作为约束条件;同时对yi做傅立叶反变换,变换到图像空间,以xi表示;再对xi进行基于正交小波基的小波变换;在小波域,用公式(1)所示的软域值算法对小波系数去噪;去噪后的小波系数经过小波反变换到图像空间;取出图像空间信号的幅值信号|xi|,与相位约束信号φ组合成新的xi′信号;xi′信号经过傅立叶变换到k-空间,形成wi信号;然后用实际采集到的y信号替换wi信号中相应位置的信号,形成yi′,将yi′作为下一次迭代的yi信号;重复迭代5次,对yi′作傅立叶反变换形成图像空间信号x。
POCS重建后恢复出的k-空间数据yi′作为参考数据,按照公式(4)所示,进行基于有限差分变换的最小化L1范数的非线性重建。需重建的图像为非线性重建中需优化的变量x,以x在有限差分域中的L1范数最小为优化目标,以x在k-空间中相对参考数据yi′的连续性为约束条件,用复共轭梯度的方法求解得到最优的重建图像
图4(a)为常规全采样参考数据采集经过傅立叶重建的图像,图4(b)为欠采样的MRA数据经过直接傅立叶重建的图像,图4(c)为本发明所用方法的重建图像。从放大的图中可以看到,欠采样的MRA数据经过直接傅立叶重建的图像质量较差,如箭头所示处,很多中小血管信息丢失了,而它们却在图4(c)中却清晰可见。
相比直接傅立叶重建方法,用本发明方法重建的图像质量更好。在高达79%的欠采样下,相比参考图像,采用本发明方法重建图像中的小动脉可见性仅略有下降。高欠采样带来数据采集时间的缩短,因此,部分回波压缩感知方法为临床MRA数据采集缩减时间提供了可能。
Claims (1)
1.一种基于部分回波压缩感知的快速磁共振成像方法,其特征在于该方法包括如下步骤:
步骤(1)采集随机变密度部分回波数据;具体方法为:
1-1、在信息量集中的k-空间中心区域高密度采集数据,在k-空间中心区域四周随机稀疏采集,生成二维随机掩模;将频率编码轴上需采集的数据点处都添加二维随机掩模,由此形成基于多层二维随机掩模的三维随机掩模;
1-2、按照步骤1-1生成的三维随机掩模采集k-空间数据;
步骤(2)基于小波域软域值去噪的凸集投影方法重建;具体方法为:
2-1、初始化序号i和k-空间数据yi,使得i=0,yi=y;其中y为采集到的欠采样数据;
2-2、取步骤2-1的k-空间数据yi的中心低频部分,做傅立叶反变换得到图像空间信号,取出该图像空间信号中的相位信号φ作为三维凸集投影方法迭代的约束条件;
2-3、对yi做傅立叶反变换到图像空间,形成图像空间信号xi;
2-4、用正交小波对xi进行离散小波变换;
2-5、用软域值法对步骤2-4中所生成的小波系数进行去噪处理;
2-6、将去噪后的小波系数经过离散小波反变换至图像空间;
2-7、将步骤2-6产生的图像空间信号xi的幅值信号|xi|与相位信号φ组合成新的图像空间数据xi′,xi′=|xi|eiφ;
2-8、将xi′经傅立叶变换到k-空间形成数据wi;
2-9、用采集到的y信号替换步骤2-8产生的wi中相应位置的信号,得到信号yi′;
2-10、令yi=yi′,重复步骤2-3至步骤2-10n次,其中n为自然数且n≥4;
2-11、将yi作傅立叶反变换至图像空间形成图像空间信号x;
步骤(3)基于有限差分变换的最小化L1范数的非线性重建;具体方法为:
3-1、采用有限差分变换对图像空间信号x进行稀疏变换;
3-3、采用复共轭梯度法求解优化目标,得到磁共振图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010272089 CN101975935B (zh) | 2010-09-03 | 2010-09-03 | 一种基于部分回波压缩感知的快速磁共振成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010272089 CN101975935B (zh) | 2010-09-03 | 2010-09-03 | 一种基于部分回波压缩感知的快速磁共振成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101975935A true CN101975935A (zh) | 2011-02-16 |
CN101975935B CN101975935B (zh) | 2012-12-05 |
Family
ID=43575837
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010272089 Expired - Fee Related CN101975935B (zh) | 2010-09-03 | 2010-09-03 | 一种基于部分回波压缩感知的快速磁共振成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101975935B (zh) |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102309328A (zh) * | 2011-10-19 | 2012-01-11 | 中国科学院深圳先进技术研究院 | 弥散张量成像方法及系统 |
CN103027682A (zh) * | 2011-12-12 | 2013-04-10 | 深圳先进技术研究院 | 动态对比度增强磁共振成像方法和系统 |
CN103033784A (zh) * | 2012-12-12 | 2013-04-10 | 厦门大学 | 一种射频脉冲控制的压缩感知磁共振成像方法 |
CN103150707A (zh) * | 2013-01-10 | 2013-06-12 | 华东师范大学 | 一种消除磁共振图像中尖峰噪声的方法 |
CN103142228A (zh) * | 2012-12-14 | 2013-06-12 | 中国科学院深圳先进技术研究院 | 压缩感知磁共振快速成像方法 |
CN103300859A (zh) * | 2013-05-31 | 2013-09-18 | 王勇 | 一种混合范数的高质量快速cs-mri成像方法 |
CN103323805A (zh) * | 2013-05-29 | 2013-09-25 | 杭州电子科技大学 | 基于小波域稀疏表示的speed快速磁共振成像方法 |
CN103356193A (zh) * | 2013-07-19 | 2013-10-23 | 哈尔滨工业大学深圳研究生院 | 压缩感知框架下的快速弥散张量成像方法及系统 |
CN104027113A (zh) * | 2014-07-01 | 2014-09-10 | 中国科学院武汉物理与数学研究所 | 一种基于先验知识和稀疏采样的肺部快速磁共振成像方法 |
CN104107044A (zh) * | 2014-06-27 | 2014-10-22 | 山东大学(威海) | 一种基于tv范数和l1范数的压缩感知磁共振图像重构方法 |
CN104665831A (zh) * | 2013-11-29 | 2015-06-03 | 株式会社日立制作所 | 图像伪影去除方法及图像伪影去除装置 |
CN104825161A (zh) * | 2015-06-04 | 2015-08-12 | 中国科学院武汉物理与数学研究所 | 基于过完备字典与先验知识的高质量肺部磁共振成像方法 |
CN104939828A (zh) * | 2014-03-28 | 2015-09-30 | 阿尔贝特-路德维希斯弗赖堡大学 | 磁共振成像方法 |
CN105578183A (zh) * | 2015-12-16 | 2016-05-11 | 西安交通大学 | 一种基于高斯混合模型的压缩传感视频编解码方法 |
CN107072575A (zh) * | 2014-10-03 | 2017-08-18 | 波士顿科学医学有限公司 | 用于映射心脏组织的医疗装置 |
CN108825205A (zh) * | 2018-04-09 | 2018-11-16 | 中国石油大学(北京) | 井下核磁共振波谱信号压缩感知采集方法与装置 |
CN108872292A (zh) * | 2018-05-25 | 2018-11-23 | 中国石油大学(北京) | 基于稀疏采样重建的多维拉普拉斯磁共振方法与装置 |
WO2019000348A1 (en) * | 2017-06-29 | 2019-01-03 | Shanghai United Imaging Healthcare Co., Ltd. | SYSTEM AND METHOD FOR ACCELERATING MAGNETIC RESONANCE IMAGING |
CN109143135A (zh) * | 2017-06-28 | 2019-01-04 | 圣纳普医疗(巴巴多斯)公司 | 在梯度运动传感成像用途中校正梯度不均匀性的方法 |
CN109596659A (zh) * | 2018-12-06 | 2019-04-09 | 上海交通大学 | 二维以及三维快速自旋回波成像方法 |
CN110461270A (zh) * | 2017-02-14 | 2019-11-15 | 阿特雷塞斯有限责任公司 | 具有压缩和/或cmos开窗的高速光学跟踪 |
WO2020098047A1 (zh) * | 2018-11-14 | 2020-05-22 | 深圳先进技术研究院 | 磁共振电影成像方法、装置、设备和存储介质 |
CN111325808A (zh) * | 2018-12-15 | 2020-06-23 | 深圳先进技术研究院 | 一种非线性梯度成像的图像重建方法及其相关设备 |
CN111640080A (zh) * | 2020-06-05 | 2020-09-08 | 郑州轻工业大学 | 一种基于高光谱全变差的cs图像去噪重建方法 |
CN112763958A (zh) * | 2020-12-10 | 2021-05-07 | 复旦大学 | 一种基于神经网络的多次激发平面回波磁共振成像方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5821751A (en) * | 1995-11-24 | 1998-10-13 | Siemens Aktiengesellschaft | Method for magnetic resonance image acquisition and reconstruction of the basic of wavelet encoding |
US20070083114A1 (en) * | 2005-08-26 | 2007-04-12 | The University Of Connecticut | Systems and methods for image resolution enhancement |
CN101051075A (zh) * | 2007-04-24 | 2007-10-10 | 骆建华 | 基于复奇异谱分析的磁共振部分k数据图像重建方法 |
CN101810480A (zh) * | 2010-04-16 | 2010-08-25 | 上海交通大学 | 基于缺失数据重构的磁共振图像截断伪影消除方法 |
-
2010
- 2010-09-03 CN CN 201010272089 patent/CN101975935B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5821751A (en) * | 1995-11-24 | 1998-10-13 | Siemens Aktiengesellschaft | Method for magnetic resonance image acquisition and reconstruction of the basic of wavelet encoding |
US20070083114A1 (en) * | 2005-08-26 | 2007-04-12 | The University Of Connecticut | Systems and methods for image resolution enhancement |
CN101051075A (zh) * | 2007-04-24 | 2007-10-10 | 骆建华 | 基于复奇异谱分析的磁共振部分k数据图像重建方法 |
CN101810480A (zh) * | 2010-04-16 | 2010-08-25 | 上海交通大学 | 基于缺失数据重构的磁共振图像截断伪影消除方法 |
Non-Patent Citations (1)
Title |
---|
MICHAEL LUSTIG ET AL: "Sparse MRI : The Application of Compressed Sensing for Rapid MR Imaging", 《MAGNETIC RESONANCE IN MEDICINE》 * |
Cited By (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102309328B (zh) * | 2011-10-19 | 2012-11-14 | 中国科学院深圳先进技术研究院 | 弥散张量成像方法及系统 |
CN102309328A (zh) * | 2011-10-19 | 2012-01-11 | 中国科学院深圳先进技术研究院 | 弥散张量成像方法及系统 |
CN103027682A (zh) * | 2011-12-12 | 2013-04-10 | 深圳先进技术研究院 | 动态对比度增强磁共振成像方法和系统 |
CN103033784B (zh) * | 2012-12-12 | 2015-01-21 | 厦门大学 | 一种射频脉冲控制的压缩感知磁共振成像方法 |
CN103033784A (zh) * | 2012-12-12 | 2013-04-10 | 厦门大学 | 一种射频脉冲控制的压缩感知磁共振成像方法 |
CN103142228A (zh) * | 2012-12-14 | 2013-06-12 | 中国科学院深圳先进技术研究院 | 压缩感知磁共振快速成像方法 |
CN103150707A (zh) * | 2013-01-10 | 2013-06-12 | 华东师范大学 | 一种消除磁共振图像中尖峰噪声的方法 |
CN103323805B (zh) * | 2013-05-29 | 2015-04-29 | 杭州电子科技大学 | 基于小波域稀疏表示的speed快速磁共振成像方法 |
CN103323805A (zh) * | 2013-05-29 | 2013-09-25 | 杭州电子科技大学 | 基于小波域稀疏表示的speed快速磁共振成像方法 |
CN103300859A (zh) * | 2013-05-31 | 2013-09-18 | 王勇 | 一种混合范数的高质量快速cs-mri成像方法 |
CN103356193A (zh) * | 2013-07-19 | 2013-10-23 | 哈尔滨工业大学深圳研究生院 | 压缩感知框架下的快速弥散张量成像方法及系统 |
CN103356193B (zh) * | 2013-07-19 | 2016-12-28 | 哈尔滨工业大学深圳研究生院 | 压缩感知框架下的快速弥散张量成像方法及系统 |
CN104665831A (zh) * | 2013-11-29 | 2015-06-03 | 株式会社日立制作所 | 图像伪影去除方法及图像伪影去除装置 |
CN104939828A (zh) * | 2014-03-28 | 2015-09-30 | 阿尔贝特-路德维希斯弗赖堡大学 | 磁共振成像方法 |
CN104107044A (zh) * | 2014-06-27 | 2014-10-22 | 山东大学(威海) | 一种基于tv范数和l1范数的压缩感知磁共振图像重构方法 |
CN104027113A (zh) * | 2014-07-01 | 2014-09-10 | 中国科学院武汉物理与数学研究所 | 一种基于先验知识和稀疏采样的肺部快速磁共振成像方法 |
CN107072575A (zh) * | 2014-10-03 | 2017-08-18 | 波士顿科学医学有限公司 | 用于映射心脏组织的医疗装置 |
CN104825161A (zh) * | 2015-06-04 | 2015-08-12 | 中国科学院武汉物理与数学研究所 | 基于过完备字典与先验知识的高质量肺部磁共振成像方法 |
CN104825161B (zh) * | 2015-06-04 | 2017-06-09 | 中国科学院武汉物理与数学研究所 | 基于过完备字典与先验知识的高质量肺部磁共振成像方法 |
CN105578183A (zh) * | 2015-12-16 | 2016-05-11 | 西安交通大学 | 一种基于高斯混合模型的压缩传感视频编解码方法 |
CN105578183B (zh) * | 2015-12-16 | 2019-10-11 | 西安交通大学 | 一种基于高斯混合模型的压缩传感视频编解码方法 |
CN110461270A (zh) * | 2017-02-14 | 2019-11-15 | 阿特雷塞斯有限责任公司 | 具有压缩和/或cmos开窗的高速光学跟踪 |
CN109143135A (zh) * | 2017-06-28 | 2019-01-04 | 圣纳普医疗(巴巴多斯)公司 | 在梯度运动传感成像用途中校正梯度不均匀性的方法 |
CN109143135B (zh) * | 2017-06-28 | 2022-04-05 | 圣纳普医疗公司 | 在梯度运动传感成像用途中校正梯度不均匀性的方法 |
CN110573896A (zh) * | 2017-06-29 | 2019-12-13 | 上海联影医疗科技有限公司 | 用于磁共振成像加速的系统和方法 |
CN110573896B (zh) * | 2017-06-29 | 2022-05-24 | 上海联影医疗科技股份有限公司 | 用于磁共振成像加速的系统和方法 |
WO2019000348A1 (en) * | 2017-06-29 | 2019-01-03 | Shanghai United Imaging Healthcare Co., Ltd. | SYSTEM AND METHOD FOR ACCELERATING MAGNETIC RESONANCE IMAGING |
US11327132B2 (en) | 2017-06-29 | 2022-05-10 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for magnetic resonance imaging acceleration |
CN108825205A (zh) * | 2018-04-09 | 2018-11-16 | 中国石油大学(北京) | 井下核磁共振波谱信号压缩感知采集方法与装置 |
CN108825205B (zh) * | 2018-04-09 | 2020-09-22 | 中国石油大学(北京) | 井下核磁共振波谱信号压缩感知采集方法与装置 |
CN108872292B (zh) * | 2018-05-25 | 2019-11-08 | 中国石油大学(北京) | 基于稀疏采样重建的多维拉普拉斯磁共振方法与装置 |
CN108872292A (zh) * | 2018-05-25 | 2018-11-23 | 中国石油大学(北京) | 基于稀疏采样重建的多维拉普拉斯磁共振方法与装置 |
WO2020098047A1 (zh) * | 2018-11-14 | 2020-05-22 | 深圳先进技术研究院 | 磁共振电影成像方法、装置、设备和存储介质 |
CN109596659B (zh) * | 2018-12-06 | 2021-07-30 | 上海交通大学 | 二维以及三维快速自旋回波成像方法 |
CN109596659A (zh) * | 2018-12-06 | 2019-04-09 | 上海交通大学 | 二维以及三维快速自旋回波成像方法 |
CN111325808A (zh) * | 2018-12-15 | 2020-06-23 | 深圳先进技术研究院 | 一种非线性梯度成像的图像重建方法及其相关设备 |
CN111325808B (zh) * | 2018-12-15 | 2021-02-26 | 深圳先进技术研究院 | 一种非线性梯度成像的图像重建方法及其相关设备 |
CN111640080A (zh) * | 2020-06-05 | 2020-09-08 | 郑州轻工业大学 | 一种基于高光谱全变差的cs图像去噪重建方法 |
CN111640080B (zh) * | 2020-06-05 | 2021-03-16 | 郑州轻工业大学 | 一种基于高光谱全变差的cs图像去噪重建方法 |
CN112763958A (zh) * | 2020-12-10 | 2021-05-07 | 复旦大学 | 一种基于神经网络的多次激发平面回波磁共振成像方法 |
CN112763958B (zh) * | 2020-12-10 | 2022-06-21 | 复旦大学 | 一种基于神经网络的多次激发平面回波磁共振成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101975935B (zh) | 2012-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101975935B (zh) | 一种基于部分回波压缩感知的快速磁共振成像方法 | |
CN102389309B (zh) | 基于压缩感知理论的磁共振图像重建的方法 | |
Quinsac et al. | Compressed sensing of ultrasound images: Sampling of spatial and frequency domains | |
CN101975936A (zh) | 一种基于cs压缩感知技术的快速磁共振成像方法 | |
Cheng et al. | Highly scalable image reconstruction using deep neural networks with bandpass filtering | |
KR101000786B1 (ko) | 움직임 추정 및 움직임 보상을 이용한 압축 센싱 기반 동적자기공명 영상기법 | |
CN106485764B (zh) | Mri图像的快速精确重建方法 | |
US10739432B2 (en) | Dynamic magnetic resonance imaging | |
CN109115820B (zh) | 基于平面回波成像的磁共振水脂分离和定量方法及装置 | |
Dar et al. | Synergistic reconstruction and synthesis via generative adversarial networks for accelerated multi-contrast MRI | |
CN109345473B (zh) | 一种基于自适应快速迭代收缩阈值算法的图像处理方法 | |
CN103027681A (zh) | 用于重构并行获取的mri图像的系统 | |
Baker et al. | Translational-invariant dictionaries for compressed sensing in magnetic resonance imaging | |
CN105678822A (zh) | 一种基于Split Bregman迭代的三正则磁共振图像重构方法 | |
Elahi et al. | Compressively sampled MR image reconstruction using generalized thresholding iterative algorithm | |
Yiasemis et al. | On retrospective k-space subsampling schemes for deep MRI reconstruction | |
Xiao et al. | SR-Net: a sequence offset fusion net and refine net for undersampled multislice MR image reconstruction | |
US20230380714A1 (en) | Method and system for low-field mri denoising with a deep complex-valued convolutional neural network | |
Yu et al. | Compressed sensing MRI using Singular Value Decomposition based sparsity basis | |
Shah et al. | Compressively sampled MR image reconstruction using hyperbolic tangent-based soft-thresholding | |
CN111415315B (zh) | 一种放射状采集扩散加权成像运动伪影校正方法 | |
Zhao et al. | JoJoNet: Joint-contrast and Joint-sampling-and-reconstruction Network for Multi-contrast MRI | |
CN102579045B (zh) | 磁共振超欠采样k数据的稀疏化成像方法 | |
CN112669400B (zh) | 基于深度学习预测和残差架构的动态mr重建方法 | |
He et al. | Compressed sensing-based simultaneous recovery of magnitude and phase MR images via dual trigonometric sparsity |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121205 Termination date: 20150903 |
|
EXPY | Termination of patent right or utility model |