CN103048632B - 一种基于压缩传感的快速高角分辨率扩散成像方法 - Google Patents
一种基于压缩传感的快速高角分辨率扩散成像方法 Download PDFInfo
- Publication number
- CN103048632B CN103048632B CN201210556410.5A CN201210556410A CN103048632B CN 103048632 B CN103048632 B CN 103048632B CN 201210556410 A CN201210556410 A CN 201210556410A CN 103048632 B CN103048632 B CN 103048632B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- diffusion
- gradient distribution
- gradient
- distribution
- 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.)
- Expired - Fee Related
Links
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及一种基于压缩传感的快速高角分辨率扩散成像方法,包括:根据球面螺旋分布方式确定扩散敏感梯度磁场分布方向;为确定的每一个扩散敏感梯度磁场分布方向制定欠采样编码策略;采用全采样方式,采集第一个扩散敏感梯度磁场分布方向的一组k空间信号作为参考像;将其余扩散敏感梯度磁场分布方向及其欠采样编码策略输入磁共振扫描仪,依次采集其k空间信号;将得到的各组k空间信号分别与参考像k空间信号作差,并采用压缩传感方法重建得到差值像,再将差值像叠加到参考像上。本发明利用螺旋分布法、欠采样策略及HARDI数据特点,采集少量k空间信号就能很好地重建出未知图像,提高了采样及图像重建速度。
Description
技术领域
本发明涉及磁共振成像技术领域,特别是一种基于压缩传感的快速高角分辨率扩散成像方法。
背景技术
以下是本领域中的一些名词解释,其使用范围仅限于本发明:
k空间信号:在磁共振成像过程,磁共振谱仪采集到的原始信号不同于图像像素所在的空间,而称为k空间,因此对应的原始信号称为k空间信号。
相位编码梯度磁场:为了对空间位置进行编码,就需要在不同的位置施加不同的磁场强度,这就要用到梯度磁场。当梯度的方向指向空间坐标系的y方向时,称这个方向为相位方向,对应的梯度磁场也就称为相位编码梯度磁场,简称相位梯度。
扩散敏感梯度磁场:在扩散张量成像中,除了通常施加的梯度磁场外,还要施加另外一种梯度磁场,这种梯度磁场能够提取出人体组织内的水分子扩散信息(大小,方向等),因此称这种梯度磁场为扩散敏感梯度磁场。
表观扩散系数:是一种在宏观上度量水分子扩散速度的指标,这一指标是在观察范围内所有水分子扩散速度大小的平均,与扩散方向无关。
扩散敏感梯度磁场的b值:是衡量扩散敏感梯度磁场强弱的一个指标,由扩散敏感梯度磁场的强度和持续时间共同决定,近似相当于扩散敏感梯度磁场强度对时间的积分。
现有技术中的磁共振成像技术主要有三个方面的发展:
一、磁共振成像技术和压缩传感技术
磁共振成像(MagneticResonanceImaging,MRI)技术是目前广泛应用的医疗检查技术。其工作原理是,人体内的水分子会在核磁共振效应下发出电磁信号,经过特定的梯度磁场标记后,磁共振谱仪采集这些信号(即k空间信号),然后经过傅里叶变换,重建出组织影像。这种重建方式需要采集的数据量与重建后的数据量是相等的,因此要想获得更高分辨率的图像,就要采集更多的数据。
压缩传感(Compressedsensing,CS)理论使用了一种非线性重建方式,只需要采集少量k空间信号,就可得到与傅里叶重建几乎相同的图像质量。
k空间信号通常是一行一行地采集的,每一行数据对应采集中的一个相位编码梯度磁场(简称相位梯度),因此可以用相位梯度数目来衡量k空间信号的数目,相位梯度数目越少,k空间中采集的数据就越少。傅里叶重建需要对k空间全采样,CS重建则可以对k空间欠采样。通常情况下,CS重建只需采集1/10的k空间信号,就能很好的重建出图像。
然而,过少的k空间信号,即便在CS重建下,也会使图像损失大量信息。但是如果两幅图像的差别比较小,并且已知其中一幅图像,那么即便只采集了很少的k空间信号,仍然可以很好的重建出未知图像。采集方式为:采集非常少量的k空间信号,将这些信号与已知图像k空间中对应的信号作差,对得到的差值做CS重建,把重建得到的图像叠加到已知图像上,就会得到未知图像。使用这种方法可以使采集多组相似图像的过程大大缩短。而利用CS重建非常重要的步骤之一是制定欠采样策略,即随机选取相位梯度(一定要随机选取)。通常的做法是,在进行磁共振扫描前,把相位梯度选择好,待到扫描时,直接把选好的梯度输入到机器中,进行扫描。本发明即是基于这一思路。
二、高角分辨率扩散成像及其与压缩传感的兼容性
扩散张量成像(DiffusionTensorImaging,DTI)是磁共振扩散成像的进一步发展,在临床诊断及脑科学研究等方面得到广泛应用。它的原理是,在通常的磁共振扫描中增加扩散敏感梯度磁场,以使图像反映出组织内水分子的扩散信息。一次DTI成像需要在不同方向上施加扩散敏感梯度,每个方向下得到的图像称为一组扩散权重像(DiffusionWeightImage,DWI)。高角分辨率扩散成像(HighAngularResolutionDiffusionImaging,HARDI)是在DTI的基础上施加更多扩散敏感梯度,得到的生理及病理信息远多于DTI,有重要临床意义及需求,但尚未应用于临床。主要原因是HARDI需要在数百甚至上千方向上采集DWI信号,导致其十分耗时,且极易受运动、畸变等因素影响。所以研究HARDI快速成像,解决其中的关键问题有重要的理论意义及应用前景。
以平面自旋回波(EchoPlanarImaging,EPI)为代表的多种快速成像方法将磁共振成像的常规检查时间由数小时缩短到数分钟,但仍不能满足不断增长的临床需求。尤其是HARDI要求在尽量多的方向上测量水分子的表观扩散系数,在研究中甚至会多达上千方向,这必然导致扫描时间延长。让病人在漫长的扫描时间内保持静止是极困难的,而且还有呼吸、心跳及不自主运动等因素的影响。
MRI数据的稀疏性使其适用基于CS理论的数据采集及图像重建方法,根据HARDI数据的特点改进当前单层欠采样方法可以大幅提高采集及图像重建速度。以脑部为例,HARDI数据中包括大量不同扩散敏感梯度磁场方向下得到的DWI图像,理论分析可知,其中只有白质部分的信号会有差异,水分子在脑中其他部位(如脑脊液,灰质)中的扩散是各向同性的,所以不同扩散敏感梯度磁场方向下这些部位的信号相同,由此可见不同扩散敏感梯度磁场方向的DWI图像之间差异很小。
因此HARDI符合使用CS的条件,可以设计一种方法,进一步把HARDI的扫描时间缩短。这也是本发明的重要思路之一。
三、扩散敏感梯度磁场分布方式
对于大量的梯度方向,如何合理的分布是一个重要的问题。目前使用的分布方式有正十二面体分布法和正二十面体分布法等,以正二十面体分布法为例,以正二十面体的中心为原点,由原点与正二十面体的顶点的连线确定方向,更多方向的分布可以通过进一步划分正二十面体的面来实现。这种方法能够均匀的分布大量方向,但是梯度对方向的数目有较严格的限制,不能把任意数量的梯度分布在各个方向上。因此,寻求一种通用的扩散敏感梯度磁场分布方法,也是磁共振成像需要解决的一个重要问题。
发明内容
本发明所要解决的技术问题是提供一种基于压缩传感的快速高角分辨率扩散成像方法,用于解决现有技术中磁共振成像存在的图像采集质量低、采集速度较慢、扫描时间长、扩散敏感梯度磁场分布方法适应性低等问题。
本发明解决上述技术问题的技术方案如下:一种基于压缩传感的快速高角分辨率扩散成像方法,包括:
步骤1,按球面螺旋分布方式在单位半球面上确定扩散敏感梯度磁场分布方向,并把半数的扩散敏感梯度磁场分布方向相对球心对称翻转到另一半球面,得到单位球面上所有的扩散敏感梯度磁场分布方向;
步骤2,将所有扩散敏感梯度磁场分布方向按螺旋线前进方向编号,并为除第一个扩散敏感梯度磁场分布方向外的其余方向制定相应的欠采样编码策略;
步骤3,采用全采样方式,采集第一个扩散敏感梯度磁场分布方向的一组k空间信号,对这组k空间信号作傅里叶重建,将重建后的结果作为参考像;
步骤4,应用步骤2中的欠采样编码策略,依次采集除第一个扩散敏感梯度磁场分布方向外的其余方向的k空间信号;
步骤5,将步骤4得到的一组欠采样的k空间信号与步骤3中的k空间信号作差,得到新的k空间欠采样信号;
步骤6,采用压缩传感方法对步骤5得到的k空间欠采样信号进行重建,重建后得到该扩散敏感梯度磁场分布方向的差值像;
步骤7,将步骤6得到的差值像叠加到步骤3得到的参考像上,得到该扩散敏感梯度磁场分布方向下的扩散权重像;
步骤8,重复步骤5至步骤7,得到其余扩散敏感梯度磁场分布方向下的扩散权重像。
在上述技术方案的基础上,本发明还可以做如下改进:
进一步,所述步骤1具体包括:
步骤11,选取单位球面的上顶点作为球面螺旋线初始点;
步骤12,计算球面螺旋线在上半球的长度随环绕圈数变化的关系式;
步骤13,计算在指定扩散敏感梯度磁场分布方向数目下,以球面螺旋线间距为两相邻方向的球面距离,排列所有方向所需的螺旋线总长度随球面螺旋线环绕圈数变化的关系式;
步骤14,由步骤12及步骤13的两关系式确定球面螺旋线的环绕圈数,并构造上半球的球面螺旋线;
步骤15,由球面螺旋线的初始点出发,沿着球面螺旋线前进的方向,依次计算出与初始点间的球面螺旋线长度为球面螺旋线间距整数倍的点的直角坐标,并按照先后顺序对计算出的点进行编号;
步骤16,将所有点按照编号分为奇数组和偶数组,把所有偶数组的点相对球心作对称翻转,翻转到下半球上,并删除原来的点;
步骤17,将球心与球面上的所有点的连线方向作为扩散敏感梯度磁场分布方向。
进一步,所述步骤2具体包括:
步骤21,向磁共振扫描仪随机输入一幅全采样图像的相位编码梯度数目及层数;
步骤22,将磁共振扫描仪的最大相位编码梯度强度平均划分为若干个等级,并将每个等级作为一个梯度的对应强度;
步骤23,按照高斯随机分布的方式,选择所有划分出的等级中的若干个等级的梯度作为梯度组;
步骤24,按照同样的方式分别为全采样图像的其他层面选择梯度组;
步骤25,将所有梯度组作为一个扩散敏感梯度磁场分布方向下的欠采样策略;
步骤26,按上述步骤为除螺旋初始点之外的扩散敏感梯度磁场分布方向制定欠采样编码策略。
进一步,所述步骤4中依次采集除第一个扩散敏感梯度磁场分布方向外的其余方向的k空间信号具体包括:
步骤41,设定扩散敏感梯度磁场的b值,并将除第一个扩散敏感梯度磁场分布方向外的其余扩散敏感梯度磁场分布方向及其对应的欠采样编码策略制作成文本文件,读入磁共振扫描仪的控制端;
步骤42,选择适合多层扫描的磁共振采集序列,按步骤41的文本文件修改磁共振采集序列的相位编码梯度,使磁共振扫描仪按照欠采样策略中的相位编码梯度采集;
步骤43,按照扩散敏感梯度磁场分布方向的编号,依次在所有扩散敏感梯度磁场分布方向下采集k空间信号。
本发明的有益效果是:本发明的一种基于压缩传感的快速高角分辨率扩散成像方法,主要有以下几个方面的优点:
一、利用欠采样策略,只需采集少量k空间信号就能很好地重建出未知图像,缩短了采集过程;
二、利用HARDI数据的特点改进当前单层欠采样方法,大幅度提高了采样及图像重建速度,同时减短了扫描时间;
三、采用螺旋分布法来选择扩散敏感梯度磁场分布方向,该方向更为通用,能把任意数量的梯度分布在各个方向上。
四、快速的HARDI技术使得被观察者在磁共振扫描仪中扫描的时间缩短,减少了运动伪影;
五、压缩传感的重建方式能够部分抑制图像中的高斯噪声,有助于图像信噪比的提高。
附图说明
图1为本发明的基于压缩传感的快速高角分辨率扩散成像方法的流程示意图;
图2为本发明中采用螺旋分布方式确定扩散敏感梯度磁场分布方向的流程示意图;
图3为本发明中制定欠采样编码策略的流程示意图;
图4为本发明采用欠采样编码策略采集扩散敏感梯度磁场分布方向下k空间信号的流程示意图;
图5为本发明采用压缩传感方法重建差值像的流程示意图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
如图1所示,实施例一是一种基于压缩传感的快速高角分辨率扩散成像方法,其实现方法如下所述。
一、步骤1,在球坐标系中构造一个单位球面(半径r=1,方位角phi∈[0,2π),仰角theta∈[0,π)),按球面螺旋分布方式在单位半球面上确定扩散敏感梯度磁场分布方向,并把半数的扩散敏感梯度磁场分布方向相对球心对称翻转到另一半球面,得到单位球面上所有的扩散敏感梯度磁场分布方向。
如图2所示,步骤1的实施过程具体包括:
步骤11,选取单位球面的上顶点作为螺旋初始点,坐标为(r=1,phi=0,theta=0)。
步骤12,计算球面螺旋线在上半球的长度s1随环绕圈数n变化的关系式,关系式为完全椭圆积分函数s1=EllipticE(-16n^2);
步骤13,计算在指定扩散敏感梯度磁场分布方向数目m下,以球面螺旋线间距d为两相邻方向的球面距离,排列所有方向所需的总长度s2随球面螺旋线环绕圈数n变化的关系式,关系式为s2=(pi*m)/2n;
步骤14,根据步骤12及步骤13的两关系式的交点即可确定球面螺旋线的环绕圈数n,并构造上半球的球面螺旋线;
步骤15,由球面螺旋线的初始点出发,沿着球面螺旋线前进的方向,依次计算出与初始点间的球面螺旋线长度为球面螺旋线间距整数倍的点的直角坐标,按照先后顺序对计算出的点进行编号;
步骤16,将所有点按照编号分为奇数组和偶数组,把所有偶数组的点相对球心作对称翻转,翻转到下半球上。只把偶数组的方向翻转到下半球,奇数组的方向留在上半球,是为了把所有方向在球面上均匀分布,可防止磁共振采集中由于方向相隔太近而出现涡流效应。
步骤17,将球面上的所有点与球心的连线方向作为扩散敏感梯度磁场分布方向。
二、步骤2,将所有扩散敏感梯度磁场分布方向按螺旋线前进方向编号,并为除第一个扩散敏感梯度磁场分布方向外的其余方向制定相应的欠采样编码策略。
如图3所示,步骤2的实施过程具体包括:
步骤21,向磁共振扫描仪随机输入一幅全采样图像的相位编码梯度数目及层数,以相位编码梯度数目为256、层数为128为例。
步骤22,将磁共振扫描仪的最大相位编码梯度强度平均划分为256个等级(包含正负),并将256个等级作为256个梯度的对应强度。
步骤23,按照高斯随机分布的方式,从256个梯度中选择少量梯度作为梯度组,这里以选择5个梯度为例。
步骤24,按照同样的方式分别为全采样图像的其他层面选择127个梯度组。
步骤25,将所有梯度组作为一个扩散敏感梯度磁场分布方向下的欠采样策略,可知欠采样策略是通过设定相位编码梯度的数目和大小来决定采集哪些数据的。
步骤26,按上述步骤为除螺旋初始点(即第一个扩散敏感梯度磁场分布方向)之外的扩散敏感梯度磁场分布方向制定欠采样编码策略。
三、步骤3,采用全采样方式,并通过磁共振扫描仪采集第一个扩散敏感梯度磁场分布方向的一组k空间信号,将这组k空间信号作为参考像。在这之前,在不加任何扩散敏感梯度的情况下(即扩散敏感梯度的b值为零),使用全采样的方式采集一组k空间数据。这里,扫描不加扩散敏感梯度磁场下的k空间数据是HARDI中必做的一部分,是为了比较有扩散敏感梯度磁场和没有扩散敏感梯度磁场的差别。
四、步骤4,将除第一个扩散敏感梯度磁场分布方向(即螺旋线初始点的方向)外的其余扩散敏感梯度磁场分布方向信息及其对应的欠采样编码策略输入磁共振扫描仪,依次采集对应扩散敏感梯度磁场分布方向下k空间信号。
如图4所示,步骤4的实施过程具体包括:
步骤41,设定扩散敏感梯度磁场的b值,并将除第一个扩散敏感梯度磁场分布方向外的其余扩散敏感梯度磁场分布方向及其对应的欠采样编码策略制作成文本文件,读入磁共振扫描仪的控制端;
步骤42,选择适合多层扫描的磁共振采集序列,按步骤41的文本文件修改磁共振采集序列的相位编码梯度,使磁共振扫描仪按照欠采样策略中的相位编码梯度采集;
步骤43,按照扩散敏感梯度磁场分布方向的编号,依次在所有扩散敏感梯度磁场分布方向下采集k空间信号。
五、步骤5,如图5所示,具体包括:
步骤51,将步骤4得到的一组欠采样的k空间信号与步骤3中的k空间信号作差,得到新的k空间欠采样信号;
步骤52,采用压缩传感方法对步骤51得到的k空间欠采样信号进行重建,重建后得到图像近似为该扩散敏感梯度磁场分布方向下的扩散权重像与参考像的差值像;
步骤53,将步骤52得到的差值像叠加到步骤3得到的参考像上,得到该扩散敏感梯度磁场分布方向下的扩散权重像;
步骤54,重复步骤51至步骤53,得到其余扩散敏感梯度磁场分布方向下的扩散权重像。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于压缩传感的快速高角分辨率扩散成像方法,其特征在于,包括:
步骤1,按球面螺旋分布方式在单位半球面上确定扩散敏感梯度磁场分布方向,并把半数的扩散敏感梯度磁场分布方向相对球心对称翻转到另一半球面,得到单位球面上所有的扩散敏感梯度磁场分布方向;
所述步骤1具体包括:
步骤11,选取单位球面的上顶点作为球面螺旋线初始点;
步骤12,计算球面螺旋线在上半球的长度随环绕圈数变化的关系式;
步骤13,计算在指定扩散敏感梯度磁场分布方向数目下,以球面螺旋线间距为两相邻方向的球面距离,排列所有方向所需的螺旋线总长度随球面螺旋线环绕圈数变化的关系式;
步骤14,由步骤12及步骤13的两关系式确定球面螺旋线的环绕圈数,并构造上半球的球面螺旋线;
步骤15,由球面螺旋线的初始点出发,沿着球面螺旋线前进的方向,依次计算出与初始点间的球面螺旋线长度为球面螺旋线间距整数倍的点的直角坐标,并按照先后顺序对计算出的点进行编号;
步骤16,将所有点按照编号分为奇数组和偶数组,把所有偶数组的点相对球心作对称翻转,翻转到下半球上,并删除原来的点;
步骤17,将球心与球面上的所有点的连线方向作为扩散敏感梯度磁场分布方向;
步骤2,将所有扩散敏感梯度磁场分布方向按螺旋线前进方向编号,并为除第一个扩散敏感梯度磁场分布方向外的其余方向制定相应的欠采样编码策略;
步骤3,采用全采样方式,采集第一个扩散敏感梯度磁场分布方向的一组k空间信号,对这组k空间信号作傅里叶重建,将重建后的结果作为参考像;
步骤4,应用步骤2中的欠采样编码策略,依次采集除第一个扩散敏感梯度磁场分布方向外的其余方向的k空间信号;
步骤5,将步骤4得到的一组欠采样的k空间信号与步骤3中的k空间信号作差,得到新的k空间欠采样信号;
步骤6,采用压缩传感方法对步骤5得到的k空间欠采样信号进行重建,重建后得到该扩散敏感梯度磁场分布方向的差值像;
步骤7,将步骤6得到的差值像叠加到步骤3得到的参考像上,得到该扩散敏感梯度磁场分布方向下的扩散权重像;
步骤8,重复步骤5至步骤7,得到其余扩散敏感梯度磁场分布方向下的扩散权重像。
2.根据权利要求1所述的快速高角分辨率扩散成像方法,其特征在于,所述步骤2具体包括:
步骤21,向磁共振扫描仪随机输入一幅全采样图像的相位编码梯度数目及层数;
步骤22,将磁共振扫描仪的最大相位编码梯度强度平均划分为若干个等级,并将每个等级作为一个梯度的对应强度;
步骤23,按照高斯随机分布的方式,选择所有划分出的等级中的若干个等级的梯度作为梯度组;
步骤24,按照同样的方式分别为全采样图像的其他层面选择梯度组;
步骤25,将所有梯度组作为一个扩散敏感梯度磁场分布方向下的欠采样策略;
步骤26,按上述步骤为除螺旋初始点之外的扩散敏感梯度磁场分布方向制定欠采样编码策略。
3.根据权利要求1所述的快速高角分辨率扩散成像方法,其特征在于,所述步骤4中依次采集除第一个扩散敏感梯度磁场分布方向外的其余方向的k空间信号具体包括:
步骤41,设定扩散敏感梯度磁场的b值,并将除第一个扩散敏感梯度磁场分布方向外的其余扩散敏感梯度磁场分布方向及其对应的欠采样编码策略制作成文本文件,读入磁共振扫描仪的控制端;
步骤42,选择适合多层扫描的磁共振采集序列,按步骤41的文本文件修改磁共振采集序列的相位编码梯度,使磁共振扫描仪按照欠采样策略中的相位编码梯度采集;
步骤43,按照扩散敏感梯度磁场分布方向的编号,依次在所有扩散敏感梯度磁场分布方向下采集k空间信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210556410.5A CN103048632B (zh) | 2012-12-19 | 2012-12-19 | 一种基于压缩传感的快速高角分辨率扩散成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210556410.5A CN103048632B (zh) | 2012-12-19 | 2012-12-19 | 一种基于压缩传感的快速高角分辨率扩散成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103048632A CN103048632A (zh) | 2013-04-17 |
CN103048632B true CN103048632B (zh) | 2015-11-25 |
Family
ID=48061340
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210556410.5A Expired - Fee Related CN103048632B (zh) | 2012-12-19 | 2012-12-19 | 一种基于压缩传感的快速高角分辨率扩散成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103048632B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104280705B (zh) * | 2014-09-30 | 2017-01-11 | 深圳先进技术研究院 | 基于压缩感知的磁共振图像重建方法和装置 |
WO2018082027A1 (zh) * | 2016-11-04 | 2018-05-11 | 深圳先进技术研究院 | 用于三维动态磁共振成像的圆柱形k空间采集方法和系统 |
CN106997034B (zh) * | 2017-04-25 | 2019-07-26 | 清华大学 | 基于以高斯模型为实例整合重建的磁共振扩散成像方法 |
WO2019000348A1 (en) * | 2017-06-29 | 2019-01-03 | Shanghai United Imaging Healthcare Co., Ltd. | SYSTEM AND METHOD FOR ACCELERATING MAGNETIC RESONANCE IMAGING |
CN107703467A (zh) * | 2017-08-21 | 2018-02-16 | 深圳先进技术研究院 | 三维动态磁共振成像的采集方法、装置、设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101283288A (zh) * | 2005-09-22 | 2008-10-08 | 威斯康星校友研究基金会 | 扩散加权mri中的高度约束背投重构过程 |
CN102298128A (zh) * | 2011-07-08 | 2011-12-28 | 哈尔滨工业大学 | 基于自适应扩散基函数分解的组织纤维束结构信息提取方法 |
CN102309328A (zh) * | 2011-10-19 | 2012-01-11 | 中国科学院深圳先进技术研究院 | 弥散张量成像方法及系统 |
CN102389309A (zh) * | 2011-07-08 | 2012-03-28 | 首都医科大学 | 基于压缩感知理论的磁共振图像重建的方法 |
CN102663701A (zh) * | 2011-12-12 | 2012-09-12 | 中国科学院深圳先进技术研究院 | 一种磁共振参数重建方法及系统 |
CN102736046A (zh) * | 2011-04-06 | 2012-10-17 | 西门子公司 | 产生扩散信息的方法和相应的磁共振设备 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8170311B2 (en) * | 2008-05-02 | 2012-05-01 | The Uwm Research Foundation, Inc. | Toeplitz random encoding for reduced acquisition in compressed sensing magnetic resonance imaging |
US7659718B1 (en) * | 2008-07-31 | 2010-02-09 | The Board Of Trustees Of The Leland Stanford Junior University | Blip design for random sampling compressed sensing of flyback 3D-MRSI |
US9121915B2 (en) * | 2010-12-09 | 2015-09-01 | The Board Of Trustees Of The Leland Stanford Junior University | Multi-dimensional cardiac and respiratory imaging with MRI |
DE102011077197B4 (de) * | 2011-06-08 | 2013-05-16 | Siemens Aktiengesellschaft | Verzeichnungskorrektur bei einer Magnetresonanz-Bildgebung |
-
2012
- 2012-12-19 CN CN201210556410.5A patent/CN103048632B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101283288A (zh) * | 2005-09-22 | 2008-10-08 | 威斯康星校友研究基金会 | 扩散加权mri中的高度约束背投重构过程 |
CN102736046A (zh) * | 2011-04-06 | 2012-10-17 | 西门子公司 | 产生扩散信息的方法和相应的磁共振设备 |
CN102298128A (zh) * | 2011-07-08 | 2011-12-28 | 哈尔滨工业大学 | 基于自适应扩散基函数分解的组织纤维束结构信息提取方法 |
CN102389309A (zh) * | 2011-07-08 | 2012-03-28 | 首都医科大学 | 基于压缩感知理论的磁共振图像重建的方法 |
CN102309328A (zh) * | 2011-10-19 | 2012-01-11 | 中国科学院深圳先进技术研究院 | 弥散张量成像方法及系统 |
CN102663701A (zh) * | 2011-12-12 | 2012-09-12 | 中国科学院深圳先进技术研究院 | 一种磁共振参数重建方法及系统 |
Non-Patent Citations (1)
Title |
---|
压缩传感理论在磁共振成像技术中的应用;王飞 等;《中国医学物理学杂志》;20121130;第29卷(第6期);第3755-3758、3833页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103048632A (zh) | 2013-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103048632B (zh) | 一种基于压缩传感的快速高角分辨率扩散成像方法 | |
CN106264529B (zh) | 用于三维动态磁共振成像的圆柱形k空间采集方法和系统 | |
CN104101849B (zh) | 磁共振系统控制序列的确定 | |
CN103238082B (zh) | 使用多点Dixon技术和低分辨率校准的MR成像 | |
CN111505553B (zh) | 磁共振成像系统和方法 | |
CN102772210B (zh) | 弥散加权磁共振成像 | |
CN103477238A (zh) | 采用来自先验采集的约束的压缩感测mr图像重建 | |
CN106997034B (zh) | 基于以高斯模型为实例整合重建的磁共振扩散成像方法 | |
CN103597370B (zh) | 空间编码的相位对比磁共振成像 | |
CN102914754A (zh) | 用于产生磁共振照片的方法 | |
DE102008007048B4 (de) | Dynamische Verzeichnungskorrektur bei EPI-Messungen in der medizinischen Magnet-Resonanz-Bildgebung | |
CN104337516A (zh) | 磁共振控制序列的确定 | |
CN103323805B (zh) | 基于小波域稀疏表示的speed快速磁共振成像方法 | |
CN110573896B (zh) | 用于磁共振成像加速的系统和方法 | |
CN102772207B (zh) | 磁共振成像方法和装置 | |
CN101975936A (zh) | 一种基于cs压缩感知技术的快速磁共振成像方法 | |
CN102654569B (zh) | 用于处理复数图像数据的方法和装置 | |
CN1683939A (zh) | 在磁共振断层造影中融合加权成像的运动校正的多点方法 | |
CN106793971A (zh) | 磁共振成像装置 | |
US20130278257A1 (en) | Mri pulse sequence based on q-space trajectory technique | |
CN106574954A (zh) | 针对epi的具有奈奎斯特伪影校正的并行mr成像 | |
CN106264530B (zh) | 用于三维动态磁共振成像的笛卡尔k空间采集方法和系统 | |
CN108459289A (zh) | 一种基于数据一致性的多次激发弥散加权磁共振成像方法 | |
CN107037385A (zh) | 数字mri图谱的构建方法及设备 | |
CN105249967B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151125 Termination date: 20161219 |
|
CF01 | Termination of patent right due to non-payment of annual fee |