CN103679654A - 磁共振成像方法和系统 - Google Patents
磁共振成像方法和系统 Download PDFInfo
- Publication number
- CN103679654A CN103679654A CN201310633874.6A CN201310633874A CN103679654A CN 103679654 A CN103679654 A CN 103679654A CN 201310633874 A CN201310633874 A CN 201310633874A CN 103679654 A CN103679654 A CN 103679654A
- Authority
- CN
- China
- Prior art keywords
- spatial data
- adopt
- learning model
- degree
- deep learning
- 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空间数据;根据所述预设采样模式训练深度学习模型,使用训练获得的深度学习模型表征K空间中已采样点和未采样点间的映射关系,根据已获得的K空间数据和所述映射关系对所述K空间进行填充,获得完整K空间数据;对所述完整K空间数据进行逆傅立叶变换,获得磁共振图像。本发明提供的磁共振成像方法,由于利用的是欠采K空间数据,可以达到快速扫描、快速成像的目的;而且由于是对完整K空间数据进行逆傅立叶变换获得磁共振图像,因此避免了卷褶伪影,成像质量高。本发明还提供了一种磁共振成像系统。
Description
技术领域
本发明涉及图像技术领域,特别是涉及一种磁共振成像方法和系统。
背景技术
K空间为磁共振图像原始数据的填充存储空间格式,K空间数据经逆傅立叶变换可获得重建磁共振图像。数据采集需满足奈奎斯特采样定理,因此需要采集许多条的相位编码线,以满足重建磁共振图像的需要,这样就严重影响了扫描速度。
为了提高磁共振成像速度,快速成像技术广泛应用于磁共振成像中。该技术主要利用磁共振信号的稀疏性或接收线圈的敏感度信息,减少成像所必需的相位编码线数,获得更快的扫描速度。我们将减少了相位编码线数的K空间数据称为欠采K空间数据。然而若直接对欠采数据进行逆傅立叶变换,得到的磁共振图像有很强的卷褶伪影,严重影响成像质量。
发明内容
基于此,有必要针对使用欠采数据直接进行逆傅立叶变换,得到的磁共振图像有很强的卷褶伪影,影响成像质量的问题,提供一种磁共振成像方法和系统。
一种磁共振成像方法,所述方法包括:
使用预设采样模式获取对应K空间的欠采K空间数据;
根据所述预设采样模式训练深度学习模型,使用训练获得的深度学习模型表征K空间中已采样点和未采样点间的映射关系,根据已获得的K空间数据和所述映射关系对所述K空间进行填充,获得完整K空间数据;
对所述完整K空间数据进行逆傅立叶变换,获得磁共振图像。
一种磁共振成像系统,所述系统包括:
欠采K空间数据获取模块,用于使用预设采样模式获取对应K空间的欠采K空间数据;
完整K空间数据获取模块,用于根据所述预设采样模式训练深度学习模型,使用训练获得的深度学习模型表征K空间中已采样点和未采样点间的映射关系,根据已获得的K空间数据和所述映射关系对所述K空间进行填充,获得完整K空间数据;
逆傅立叶变换模块,用于对所述完整K空间数据进行逆傅立叶变换,获得磁共振图像。
上述磁共振成像方法,根据预设采样模式训练深度学习模型,从而使用训练获得的深度学习模型估计K空间中已采样点和未采样点之间的映射关系,从而根据给映射关系填充K空间,获得完整K空间数据。最后对完整K空间数据进行逆傅立叶变换获得磁共振图像。由于利用的是欠采K空间数据,可以达到快速扫描、快速成像的目的;而且由于是对完整K空间数据进行逆傅立叶变换获得磁共振图像,因此避免了卷褶伪影,成像质量高。
附图说明
图1为一个实施例中磁共振成像方法的流程示意图;
图2为一个实施例中K空间的示意图;
图3为一个实施例中根据所述预设采样模式训练深度学习模型,使用训练获得的深度学习模型表征K空间中已采样点和未采样点间的映射关系,根据已获得的K空间数据和映射关系对K空间进行填充,获得完整K空间数据的步骤的流程示意图;
图4为一个举例中用6个基向量表示一个一维向量的示意图;
图5为一个实施例中使用稀疏自编码器求解用以表示欠采K空间数据的基向量的示意图;
图6为一个实施例中磁共振成像系统的结构框图;
图7为一个实施例中完整K空间数据获取模块的结构框图;
图8为一个实施例中深度学习模型训练模块的结构框图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
除非上下文另有特定清楚的描述,本发明中的元件和组件,数量既可以单个的形式存在,也可以多个的形式存在,本发明并不对此进行限定。本发明中的步骤虽然用标号进行了排列,但并不用于限定步骤的先后次序,除非明确说明了步骤的次序或者某步骤的执行需要其他步骤作为基础,否则步骤的相对次序是可以调整的。可以理解,本文中所使用的术语“和/或”涉及且涵盖相关联的所列项目中的一者或一者以上的任何和所有可能的组合。
如图1所示,在一个实施例中,提供了一种磁共振成像方法,该方法包括:
步骤102,使用预设采样模式获取对应K空间的欠采K空间数据。
可通过磁共振扫描仪获取欠采的K空间数据。预设采样模式可为任意的欠采采样模式,比如变密度随机采样模式。K空间为磁共振图像原始数据的填充存储空间格式,K空间数据经逆傅立叶变换可获得重建磁共振图像。如图2所示,假设全采的K空间为M×N的矩阵,行方向为频率编码方向,列方向为相位编码方向,只在相位编码方向欠采,没有采集的部分用0来填充。假设加速因子为2,则采集的数据为M/2行。其中,K空间中心数据(如图2中202部分,可为预设面积或预设行数构成的区域)为全采,其余部分欠采,对采样模式没有限制。例如图2所采用的采样模式是变密度采样,图2中白色表示采集的部分,黑色表示未采集(欠采)的部分。
步骤104,根据预设采样模式训练深度学习模型,使用训练获得的深度学习模型表征K空间中已采样点和未采样点间的映射关系,根据已获得的K空间数据和映射关系对K空间进行填充,获得完整K空间数据。
K空间数据具有对称性,而且各条相位编码线具有相关性,因此可利用已采数据估计未采数据,从而获得完整K空间数据。具体地,可使用欠采K空间数据训练深度学习模型,并使用训练获得的深度学习模型估计已采样点和未采样点间的映射关系,从而根据该映射关系和已采样数据可获得K空间中的未采样数据,并填入K空间中,从而可获得完整的K空间数据。
在一个实施例中,根据预设采样模式训练深度学习模型的步骤包括:对K空间的全采中心数据按照预设采样模式进行欠采,获得欠采中心数据;用全采中心数据和欠采中心数据训练深度学习模型。
本实施例中,采用欠采K空间数据本身进行训练,因为K空间的中心数据是全采的,对全采中心数据按与采集欠采K空间数据时相同的采样模式进行欠采集,这样就获得了对应全采中心数据的欠采中心数据。然后用K空间的全采中心数据和获得的欠采中心数据训练深度学习模型,从而可利用全采中心数据和欠采中心数据的对应关系,估计出K空间中已采样点和未采样点之间的映射关系,从而便于后续根据该映射关系获得完整K空间数据。
在一个实施例中,根据预设采样模式训练深度学习模型的步骤包括:采集一组全采K空间数据;按照预设采样模式对全采K空间数据进行欠采,获得训练欠采数据;用全采K空间数据和训练欠采数据训练深度学习模型。
本实施例中,先采集一组全采K空间数据,该全采K空间数据经傅立叶逆变换可以获得重建后的图像。用于采集欠采K空间数据时相同的采样模式对全采K空间数据进行欠采集,这样就获得了对应全采K空间数据的训练欠采数据。然后用全采K空间数据和训练欠采数据训练深度学习模型,从而可利用全采K空间数据和训练欠采数据的对应关系,估计出欠采K空间数据对应的K空间中已采样点和未采样点之间的映射关系,从而便于后续根据该映射关系获得完整K空间数据。
深度学习(Deep Learning)是机器学习的一个分支,源于人工神经网络的研究,含有多隐层的多层感知器就是一种深度学习结构,深度学习通过组合低层特征形成更加抽象的高层表示属性类别或特征,以发现数据的分布式特征表示。深度学习模型需要多层来获得更抽象的特征表达。
具体地,深度学习模型包括若干层,每一层包括若干神经元单元,每个神经元单元为一个运算单元,在每一层从这些神经元单元的输入端输入数据,神经元单元的输出端输出数据。宏观上来讲,从深度学习模型的第一层输入数据,然后从第一层开始将每一层输出的特征作为下一层的输入,逐层抽象,最终从最后一层输出数据。通过逐层的特征变换,将样本在原空间的特征表示变换到一个新特征空间,从而使预测、估计更加容易。
如图3所示,在一个实施例中,步骤104包括:
步骤302,根据预设采样模式训练深度学习模型。
在一个实施例中,根据预设采样模式训练深度学习模型的步骤包括步骤1)~2):
1)使用超完备基向量表示训练样本,训练样本包括全采数据和按照预设采样模式对全采数据进行欠采获得的欠采数据。
训练样本包括全采数据和欠采数据,其中欠采数据是按照预设采样模式对全采数据进行欠采获得的。利用该全采数据和欠采数据可以估计出K空间中已采样点和未采样点之间的映射关系,从而利用该映射关系获得完整K空间数据。
对于任意的n维输入向量x,可以表示为一组基向量φi的线性组合,即其中φi为一组基向量,当基向量φi的总数k大于输入向量x维数n时,这组基向量就为超完备基向量。如图4为n=1,k=6的一个例子。用以表示欠采K空间数据的基向量φi与其系数di可以用以下方法求得。
定义代价函数(cost function)为其中m表示训练样本数,k表示基向量φi的总数。通过大量的训练样本xj,j=1,2,3,……m,使得代价函数的误差最小。由于di和φi都是未知,无法直接利用梯度下降算法去现实,但是可以通过一次固定一个参数如di,对φi进行梯度下降算法,来最小化误差。把得到的φi带回代价函数,此时固定φi,对di进行梯度下降法,使得代价函数最小,此为一次迭代,反复迭代直到误差小于设定的门限阈值,如0.0001,迭代终止,从而得到最终的di和φi。
进一步地,可定义s来表示φi的平均激活度,s为常数。例如s=0.1,表示对于100个x,每个φi只使用了10次。这样只需求φi,现在就可以利用稀稀疏自编码来求得φi了。
使用稀疏自编码器时,φi的表示如图5所示,对于输入x(r)(训练样本),有其中,r表示第r个训练样本,且样本总数为m。是激励函数,激励函数表示一个神经元单元输入与输出之间的函数关系,z表示自变量。hi表示用于计算φi的中间值,k表示基向量φi的总数。wij表示如图5所示的神经网络中连接第一层第j个神经元单元和第2层第i个神经元单元的参数;Bi表示如图5所示的神经网络中第2层第i个神经元单元的偏置项。由此可以看到如图5所示的神经网络的前半部分近似求得φi,后半部分网络输出由此我们重新定义目标函数而为了得到稀疏表示,我们先定义表示hi的平均激活度,求得与s比较,当的激活度小于或者大于s时给予惩罚,这样我们就能使得hi的平均激活度尽量接近s。
最小化目标函数来求得神经网络各层系数。如何刻画φi表示原始数据的能力,我们把φi前面的网络参数看成是对原始数据的一种编码方式,通过这组系数得到了φi,而右边的网络看作是对原始数据的译码。换句话说,通过左边的网络得到的一组φi必须有能力去恢复原始数据,所以我们只要把右边的网络输出为与输入相等。这样一个三层的神经网络可以提取一组超完备基向量,进而使用该一组超完备基向量表示训练样本。
2)将使用超完备基向量表示的训练样本输入深度学习模型的第一层,并从第一层开始,将深度学习模型每一层输出的特征作为下一层的输入从而对每一层进行训练,直至获得深度学习模型所有层的参数。
可利用梯度下降法来求得深度学习模型所有层的参数W和b,如何求下面会详细讲解。神经元单元的激励函数是当激励函数输出接近1时认为此神经元单元被激活,而激励函数输出为0是认为该神经元单元被抑止。
表示深度学习模型第2层第j个神经元单元的激活度,但是这一表示方法中并未明确指出哪一个输入x带来了这一激活度,所以可以用来表示神经元j在训练集上的平均活跃度。用自定常数ρ表示稀疏性,ρ可以为0.02。并且使得为了实现这一限制,即当与ρ不同时给予惩罚(减小权重)。
现在修改代价函数,其中,Sl表示深度学习模型中第l层的神经元单元个数(偏置神经元单元不计算在内)。β为权重系数,用于控制稀疏性惩罚因子的权重。项则也(间接地)取决于W,b,因为是隐藏神经元j的平均激活度,而隐藏层神经元的激活度取决于W和b.
目标是针对参数W和b来利用梯度下降法来最小化代价函数Jsparse(W,b)。
而整体代价函数Jsparse(W,b)的梯度可表示为
其中,m表示训练样本总数,x(i)表示第i个样本数据,y(i)是对应x(i)的深度学习模型的输出值。表示深度学习模型中连接第l层第j个神经元单元和第l+1层第i个神经元单元的参数;λ是权重衰减系数,用于控制的权重。表示第l+1层第i个神经元单元的偏置项,可以看作是连接第l层偏置神经元单元和第l+1层第i个神经元单元的参数。
计算梯度可以使用反向传播算法。
对于深度学习模型的第l层,令ΔW(l):=0,Δb(l):=0,表示将ΔW(l)和Δb(l)设置为全零矩阵或全零向量,其中符号:=表示赋值;ΔW(l)是一个与矩阵W(l)维度相同的矩阵,Δb(l)是一个与矩阵b(l)维度相同的矩阵。
计算
计算
更新权重参数:
其中,α表示学习速率。
重复梯度下降法的迭代步骤来减小代价函数的值,从而获得深度学习模型所有层的参数W(l)和b(l),完成训练。
步骤304,将欠采K空间数据输入训练获得的深度学习模型获得输出值,将输出值填入K空间中的对应位置,获得填充K空间数据。
具体地,由于训练获得的深度学习模型估计出了K空间中已采样点和未采样点之间的映射关系,将K空间中各个已采样点的数据输入训练获得的深度学习模型获得对应已采样点的输出值。该输出值即为估计的K空间中未采样点的值,将该输出值填入K空间中的对应位置,便可获得填充K空间数据。
步骤306,根据填充K空间数据和预设采样模式训练深度学习模型,继续执行将欠采K空间数据输入训练获得的深度学习模型获得输出值,将输出值填入K空间中的对应位置,获得填充K空间数据的步骤,直至满足迭代终止条件,将迭代终止后获得的填充K空间数据作为完整K空间数据。
由于一次填充估计出的未采样点的值可能未必精确,因此,需要再根据填充K空间数据和预设采样模式对深度学习模型继续进行训练。具体地,可将填充K空间数据作为全采数据,结合欠采K空间数据对深度学习模型继续进行训练。其中欠采K空间数据是按照预设采样模式采集的。然后将欠采K空间数据输入该继续训练所获得的深度学习模型获得输出值,将输出值填入K空间中的对应位置,获得填充K空间数据。这样进行迭代,直到满足迭代终止条件则终止迭代,并将迭代终止后获得的填充K空间数据作为完整K空间数据。
在一个实施例中,迭代终止条件包括:达到预设迭代次数或者本次迭代获得的填充K空间数据与上一次迭代获得的填充K空间数据间的误差小于误差阈值。
可以预先设置迭代次数,当达到预设迭代次数时,可认为获得的填充K空间数据已经很接近真实值,可停止迭代运算。
或者可通过比较每一次迭代获得的填充K空间数据与上一次迭代获得的填充K空间数据,判断两者的误差是否小于误差阈值,若是则说明获得的填充K空间数据已经很接近真实值,可停止迭代运算。
步骤106,对完整K空间数据进行逆傅立叶变换,获得磁共振图像。
K空间中每个采样点的数据都包含全层的信息,其中K空间中心数据的采样点数据主要决定磁共振图像的对比,采样时中心数据为全采,因此不会影响磁共振图像的对比信息。而K空间中中心数据之外的采样点数据主要决定磁共振图像的细节,由于使用深度学习模型估计出了欠采K空间中欠采的数据,获得了完整K空间数据,保留了磁共振图像的细节信息,避免出现卷褶伪影。
上述磁共振成像方法,根据预设采样模式训练深度学习模型,从而使用训练获得的深度学习模型估计K空间中已采样点和未采样点之间的映射关系,从而根据已获得的K空间数据和该映射关系填充K空间,获得完整K空间数据。最后对完整K空间数据进行逆傅立叶变换获得磁共振图像。由于利用的是欠采K空间数据,可以达到快速扫描、快速成像的目的;而且由于是对完整K空间数据进行逆傅立叶变换获得磁共振图像,因此避免了卷褶伪影,成像质量高。
除了上述基于深度学习的磁共振成像方法,也可以使用并行成像方法、压缩感知方法等方法,根据欠采K空间数据恢复出完整的原始信号,获得无伪影的图像。其中,并行成像方法是利用接收线圈的空间敏感度信息减少扫描的线数的方法,常用的并行成像方法(如GRAPPA,SENSE等)随加速因子增高,SNR随之下降,故加速因子一般小于3。压缩感知方法是利用随机测量矩阵把一个稀疏的高维信号投影到低维的空间上,这样的随机投影包含了待重建信号的足够信息,利用信号的稀疏性等先验知识,通过非线性重建算法能以很高的概率获得原始的稀疏信号。以压缩感知理论为基础,设计随机测量矩阵在K空间中稀疏采样,只须采集少量数据就可以经过非线性重建方法重建出高质量的磁共振图像。
但使用并行成像方法、压缩感知方法等方法在加速因子较高时(大于等于3)会出现较明显的卷褶伪影或SNR(信噪比)下降。上述磁共振成像方法利用深度学习的多层特征提取机制,非线性的对欠采的K空间数据进行填充,可获得更高的加速因子,进一步提高采集速度。而且并行成像方法和压缩感知方法对采样轨迹有一定要求,但上述磁共振成像方法对采样轨迹无具体要求,也可以方便的用于非笛卡尔采样模式,具有极大的应用范围。
如图6所示,在一个实施例中,提供了一种磁共振成像系统,该系统包括:欠采K空间数据获取模块602、完整K空间数据获取模块604和逆傅立叶变换模块606。
欠采K空间数据获取模块602用于使用预设采样模式获取对应K空间的欠采K空间数据。
欠采K空间数据获取模块602可用于通过磁共振扫描仪获取欠采的K空间数据。K空间为磁共振图像原始数据的填充存储空间格式,K空间数据经逆傅立叶变换可获得重建磁共振图像。如图2所示,假设全采的K空间为M×N的矩阵,行方向为频率编码方向,列方向为相位编码方向,只在相位编码方向欠采,没有采集的部分用0来填充。假设加速因子为2,则采集的数据为M/2行。其中,K空间中心数据(如图2中202部分)为全采,其余部分欠采,对采样模式没有限制。例如图2所采用的采样模式是变密度采样,图2中白色表示采集的部分,黑色表示未采集(欠采)的部分。
完整K空间数据获取模块604用于根据预设采样模式训练深度学习模型,使用训练获得的深度学习模型表征K空间中已采样点和未采样点间的映射关系,根据已获得的K空间数据和映射关系对K空间进行填充,获得完整K空间数据。
K空间数据具有对称性,而且各条相位编码线具有相关性,因此完整K空间数据获取模块604可用于利用已采数据估计未采数据,从而获得完整K空间数据。具体地,可使用欠采K空间数据训练深度学习模型,并使用训练获得的深度学习模型估计已采样点和未采样点间的映射关系,从而根据该映射关系和已采样数据可获得K空间中的未采样数据,并填入K空间中,从而可获得完整的K空间数据。
在一个实施例中,完整K空间数据获取模块604还用于对K空间的全采中心数据按照预设采样模式进行欠采,获得欠采中心数据;用于用全采中心数据和欠采中心数据训练深度学习模型。
本实施例中,采用欠采K空间数据本身进行训练,因为K空间的中心数据是全采的,对全采中心数据按与采集欠采K空间数据时相同的预设采样模式进行欠采集,这样就获得了对应全采中心数据的欠采中心数据。然后用K空间的全采中心数据和获得的欠采中心数据训练深度学习模型,从而可利用全采中心数据和欠采中心数据的对应关系,估计出K空间中已采样点和未采样点之间的映射关系,从而便于后续根据该映射关系获得完整K空间数据。
在一个实施例中,完整K空间数据获取模块604还用于采集一组全采K空间数据;用于按照预设采样模式对全采K空间数据进行欠采,获得训练欠采数据;用于用全采K空间数据和训练欠采数据训练深度学习模型。
本实施例中,先采集一组全采K空间数据,该全采K空间数据经傅立叶逆变换可以获得重建后的图像。用于采集欠采K空间数据时相同的预设采样模式对全采K空间数据进行欠采集,这样就获得了对应全采K空间数据的训练欠采数据。然后用全采K空间数据和训练欠采数据训练深度学习模型,从而可利用全采K空间数据和训练欠采数据的对应关系,估计出欠采K空间数据对应的K空间中已采样点和未采样点之间的映射关系,从而便于后续根据该映射关系获得完整K空间数据。
逆傅立叶变换模块606用于对完整K空间数据进行逆傅立叶变换,获得磁共振图像。
K空间中每个采样点的数据都包含全层的信息,其中K空间中心数据的采样点数据主要决定磁共振图像的对比,采样时中心数据为全采,因此不会影响磁共振图像的对比信息。而K空间中心数据之外的采样点数据主要决定磁共振图像的细节,由于使用深度学习模型估计出了欠采K空间中欠采的数据,获得了完整K空间数据,保留了磁共振图像的细节信息,避免出现卷褶伪影。
如图7所示,在一个实施例中,完整K空间数据获取模块604包括深度学习模型训练模块604a和K空间填充模块604b。
深度学习模型训练模块604a用于根据预设采样模式训练深度学习模型。
K空间填充模块604b用于将欠采K空间数据输入训练获得的深度学习模型获得输出值,将输出值填入K空间中的对应位置,获得填充K空间数据。
具体地,由于训练获得的深度学习模型估计出了K空间中已采样点和未采样点之间的映射关系,将K空间中各个已采样点的数据输入训练获得的深度学习模型获得对应已采样点的输出值。该输出值即为估计的K空间中未采样点的值,将该输出值填入K空间中的对应位置,便可获得填充K空间数据。
深度学习模型训练模块604a还用于根据填充K空间数据和预设采样模式训练深度学习模型,K空间填充模块604b用于执行将欠采K空间数据输入训练获得的深度学习模型获得输出值,将输出值填入K空间中的对应位置,获得填充K空间数据,直至满足迭代终止条件,将迭代终止后获得的填充K空间数据作为完整K空间数据。
由于一次填充估计出的未采样点的值可能未必精确,因此,深度学习模型训练模块604a用于再根据填充K空间数据和预设采样模式对深度学习模型进行训练,K空间填充模块604b用于将填充K空间数据输入训练获得的深度学习模型获得输出值,将输出值填入K空间中的对应位置,获得填充K空间数据。于是深度学习模型训练模块604a和K空间填充模块604b交替作用,直到满足迭代终止条件则终止迭代,并将迭代终止后获得的填充K空间数据作为完整K空间数据。
在一个实施例中,迭代终止条件包括:达到预设迭代次数或者本次迭代获得的填充K空间数据与上一次迭代获得的填充K空间数据间的误差小于误差阈值。
可以预先设置迭代次数,当达到预设迭代次数时,可认为获得的填充K空间数据已经很接近真实值,可停止迭代运算。
或者可通过比较每一次迭代获得的填充K空间数据与上一次迭代获得的填充K空间数据,判断两者的误差是否小于误差阈值,若是则说明获得的填充K空间数据已经很接近真实值,可停止迭代运算。
如图8所示,在一个实施例中,深度学习模型训练模块604a包括欠采K空间数据表示模块604a1和逐层训练模块604a2。
欠采K空间数据表示模块604a1,用于使用超完备基向量表示训练样本,训练样本包括全采数据和按照预设采样模式对全采数据进行欠采获得的欠采数据。
逐层训练模块604a2,用于将使用超完备基向量表示的训练样本输入深度学习模型的第一层,从第一层开始,将深度学习模型每一层输出的特征作为下一层的输入从而对每一层进行训练,直至获得深度学习模型所有层的参数。
上述磁共振成像系统,根据预设采样模式训练深度学习模型,从而使用训练获得的深度学习模型估计K空间中已采样点和未采样点之间的映射关系,从而根据给映射关系填充K空间,获得完整K空间数据。最后对完整K空间数据进行逆傅立叶变换获得磁共振图像。由于利用的是欠采K空间数据,可以达到快速扫描、快速成像的目的;而且由于是对完整K空间数据进行逆傅立叶变换获得磁共振图像,因此避免了卷褶伪影,成像质量高。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (12)
1.一种磁共振成像方法,所述方法包括:
使用预设采样模式获取对应K空间的欠采K空间数据;
根据所述预设采样模式训练深度学习模型,使用训练获得的深度学习模型表征K空间中已采样点和未采样点间的映射关系,根据已获得的K空间数据和所述映射关系对所述K空间进行填充,获得完整K空间数据;
对所述完整K空间数据进行逆傅立叶变换,获得磁共振图像。
2.根据权利要求1所述的方法,其特征在于,所述根据所述预设采样模式训练深度学习模型,包括:
对K空间的全采中心数据按照所述预设采样模式进行欠采,获得欠采中心数据;
用所述全采中心数据和所述欠采中心数据训练深度学习模型。
3.根据权利要求1所述的方法,其特征在于,所述根据所述预设采样模式训练深度学习模型,包括:
采集一组全采K空间数据;
按照所述预设采样模式对全采K空间数据进行欠采,获得训练欠采数据;
用所述全采K空间数据和所述训练欠采数据训练深度学习模型。
4.根据权利要求1所述的方法,其特征在于,所述根据所述预设采样模式训练深度学习模型,使用训练获得的深度学习模型表征K空间中已采样点和未采样点间的映射关系,根据已获得的K空间数据和所述映射关系对所述K空间进行填充,获得完整K空间数据,包括:
根据所述预设采样模式训练深度学习模型;
将所述欠采K空间数据输入训练获得的深度学习模型获得输出值,将所述输出值填入所述K空间中的对应位置,获得填充K空间数据;
根据所述填充K空间数据和预设采样模式训练所述深度学习模型,继续执行将所述欠采K空间数据输入训练获得的深度学习模型获得输出值,将所述输出值填入所述K空间中的对应位置,获得填充K空间数据的步骤,直至满足迭代终止条件,将迭代终止后获得的所述填充K空间数据作为完整K空间数据。
5.根据权利要求4所述的方法,其特征在于,所述迭代终止条件包括:
达到预设迭代次数或者本次迭代获得的填充K空间数据与上一次迭代获得的填充K空间数据间的误差小于误差阈值。
6.根据权利要求4所述的方法,其特征在于,所述根据所述预设采样模式训练深度学习模型,包括:
使用超完备基向量表示训练样本,所述训练样本包括全采数据和按照所述预设采样模式对所述全采数据进行欠采获得的欠采数据;
将所述使用超完备基向量表示的训练样本输入所述深度学习模型的第一层,并从所述第一层开始,将所述深度学习模型每一层输出的特征作为下一层的输入从而对每一层进行训练,直至获得所述深度学习模型所有层的参数。
7.一种磁共振成像系统,其特征在于,所述系统包括:
欠采K空间数据获取模块,用于使用预设采样模式获取对应K空间的欠采K空间数据;
完整K空间数据获取模块,用于根据所述预设采样模式训练深度学习模型,使用训练获得的深度学习模型表征K空间中已采样点和未采样点间的映射关系,根据已获得的K空间数据和所述映射关系对所述K空间进行填充,获得完整K空间数据;
逆傅立叶变换模块,用于对所述完整K空间数据进行逆傅立叶变换,获得磁共振图像。
8.根据权利要求7所述的系统,其特征在于,所述完整K空间数据获取模块还用于对K空间的全采中心数据按照所述预设采样模式进行欠采,获得欠采中心数据;用于用所述全采中心数据和所述欠采中心数据训练深度学习模型。
9.根据权利要求7所述的系统,其特征在于,所述完整K空间数据获取模块还用于采集一组全采K空间数据;用于按照所述预设采样模式对全采K空间数据进行欠采,获得训练欠采数据;用于用所述全采K空间数据和所述训练欠采数据训练深度学习模型。
10.根据权利要求7所述的系统,其特征在于,所述完整K空间数据获取模块包括深度学习模型训练模块和K空间填充模块:
所述深度学习模型训练模块用于根据所述预设采样模式训练深度学习模型;
所述K空间填充模块用于将所述欠采K空间数据输入训练获得的深度学习模型获得输出值,将所述输出值填入所述K空间中的对应位置,获得填充K空间数据;
所述深度学习模型训练模块还用于根据所述填充K空间数据和预设采样模式训练所述深度学习模型,所述K空间填充模块用于执行将所述欠采K空间数据输入训练获得的深度学习模型获得输出值,将所述输出值填入所述K空间中的对应位置,获得填充K空间数据,直至满足迭代终止条件,将迭代终止后获得的所述填充K空间数据作为完整K空间数据。
11.根据权利要求10所述的系统,其特征在于,所述迭代终止条件包括:
达到预设迭代次数或者本次迭代获得的填充K空间数据与上一次迭代获得的填充K空间数据间的误差小于误差阈值。
12.根据权利要求10所述的系统,其特征在于,所述深度学习模型训练模块包括:
欠采K空间数据表示模块,用于使用超完备基向量表示训练样本,所述训练样本包括全采数据和按照所述预设采样模式对所述全采数据进行欠采获得的欠采数据;
逐层训练模块,用于将所述使用超完备基向量表示的训练样本输入所述深度学习模型的第一层,从所述第一层开始,将所述深度学习模型每一层输出的特征作为下一层的输入从而对每一层进行训练,直至获得所述深度学习模型所有层的参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310633874.6A CN103679654B (zh) | 2013-11-29 | 2013-11-29 | 磁共振成像方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310633874.6A CN103679654B (zh) | 2013-11-29 | 2013-11-29 | 磁共振成像方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103679654A true CN103679654A (zh) | 2014-03-26 |
CN103679654B CN103679654B (zh) | 2017-02-15 |
Family
ID=50317110
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310633874.6A Active CN103679654B (zh) | 2013-11-29 | 2013-11-29 | 磁共振成像方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103679654B (zh) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106232003A (zh) * | 2014-04-21 | 2016-12-14 | 株式会社日立制作所 | 图像拍摄装置以及方法 |
WO2017113205A1 (zh) * | 2015-12-30 | 2017-07-06 | 中国科学院深圳先进技术研究院 | 一种基于深度卷积神经网络的快速磁共振成像方法及装置 |
CN106970343A (zh) * | 2017-04-11 | 2017-07-21 | 深圳先进技术研究院 | 一种磁共振成像方法及装置 |
CN107064845A (zh) * | 2017-06-06 | 2017-08-18 | 深圳先进技术研究院 | 基于深度卷积网的一维部分傅里叶并行磁共振成像方法 |
CN107993271A (zh) * | 2017-12-26 | 2018-05-04 | 上海交通大学 | 一种磁共振动态成像采样方法和图像重建方法 |
CN108335339A (zh) * | 2018-04-08 | 2018-07-27 | 朱高杰 | 一种基于深度学习和凸集投影的磁共振重建方法 |
CN108535675A (zh) * | 2018-04-08 | 2018-09-14 | 朱高杰 | 一种基于深度学习和数据自洽的磁共振多通道重建方法 |
CN109863512A (zh) * | 2016-09-01 | 2019-06-07 | 通用医疗公司 | 用于通过流形近似进行自动变换的系统和方法 |
CN110146834A (zh) * | 2019-05-16 | 2019-08-20 | 山东大学齐鲁医院 | 一种磁共振成像方法和装置 |
US10386435B2 (en) | 2013-11-06 | 2019-08-20 | The Research Foundation For The State University Of New York | Methods and systems for fast auto-calibrated reconstruction with random projection in parallel MRI |
CN110664378A (zh) * | 2019-10-28 | 2020-01-10 | 中国科学院深圳先进技术研究院 | 磁共振成像方法、装置、系统及存储介质 |
CN110988763A (zh) * | 2019-11-05 | 2020-04-10 | 深圳先进技术研究院 | 磁共振成像方法、装置、服务器及存储介质 |
CN110992433A (zh) * | 2018-10-02 | 2020-04-10 | 通用电气公司 | 用于利用空数据一致性来稀疏图像重建的系统和方法 |
WO2020135014A1 (zh) * | 2018-12-27 | 2020-07-02 | 深圳先进技术研究院 | 医学成像模型的建立方法、装置、设备及存储介质 |
CN111513716A (zh) * | 2019-02-05 | 2020-08-11 | 通用电气精准医疗有限责任公司 | 使用扩展灵敏度模型和深度神经网络进行磁共振图像重建的方法和系统 |
CN111624540A (zh) * | 2019-02-27 | 2020-09-04 | 上海东软医疗科技有限公司 | 磁共振成像方法和装置 |
CN111681297A (zh) * | 2020-05-28 | 2020-09-18 | 上海联影智能医疗科技有限公司 | 图像重建方法、计算机设备和存储介质 |
US10989779B2 (en) | 2017-09-29 | 2021-04-27 | Yonsei University, University - Industry Foundation (UIF) | Apparatus and method for reconstructing magnetic resonance image using learning, and under-sampling apparatus method and recording medium thereof |
CN113933773A (zh) * | 2021-09-30 | 2022-01-14 | 深圳先进技术研究院 | 基于深度学习的磁共振成像方法、系统、终端及存储介质 |
CN114325524A (zh) * | 2020-09-29 | 2022-04-12 | 上海联影医疗科技股份有限公司 | 磁共振图像重建方法、装置、系统及存储介质 |
WO2023050249A1 (zh) * | 2021-09-30 | 2023-04-06 | 深圳先进技术研究院 | 基于深度学习的磁共振成像方法、系统、终端及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102798829A (zh) * | 2012-08-14 | 2012-11-28 | 深圳先进技术研究院 | 基于机器学习的并行磁共振成像grappa方法 |
-
2013
- 2013-11-29 CN CN201310633874.6A patent/CN103679654B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102798829A (zh) * | 2012-08-14 | 2012-11-28 | 深圳先进技术研究院 | 基于机器学习的并行磁共振成像grappa方法 |
Non-Patent Citations (2)
Title |
---|
YUQIONG DING等: "Noise behavior of MR brain reconstructions using compressed sensing", 《ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY (EMBC), 2013 35TH ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE》 * |
翁卓等: "基于k空间加速采集的磁共振成像技术", 《中国生物医学工程学报》 * |
Cited By (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10386435B2 (en) | 2013-11-06 | 2019-08-20 | The Research Foundation For The State University Of New York | Methods and systems for fast auto-calibrated reconstruction with random projection in parallel MRI |
CN106232003B (zh) * | 2014-04-21 | 2019-06-07 | 株式会社日立制作所 | 图像拍摄装置以及方法 |
US11419555B2 (en) | 2014-04-21 | 2022-08-23 | Fujifilm Healthcare Corporation | Image capturing apparatus and method |
CN106232003A (zh) * | 2014-04-21 | 2016-12-14 | 株式会社日立制作所 | 图像拍摄装置以及方法 |
WO2017113205A1 (zh) * | 2015-12-30 | 2017-07-06 | 中国科学院深圳先进技术研究院 | 一种基于深度卷积神经网络的快速磁共振成像方法及装置 |
CN107182216A (zh) * | 2015-12-30 | 2017-09-19 | 中国科学院深圳先进技术研究院 | 一种基于深度卷积神经网络的快速磁共振成像方法及装置 |
CN107182216B (zh) * | 2015-12-30 | 2019-06-14 | 中国科学院深圳先进技术研究院 | 一种基于深度卷积神经网络的快速磁共振成像方法及装置 |
CN109863512B (zh) * | 2016-09-01 | 2023-10-20 | 通用医疗公司 | 用于通过流形近似进行自动变换的系统和方法 |
CN109863512A (zh) * | 2016-09-01 | 2019-06-07 | 通用医疗公司 | 用于通过流形近似进行自动变换的系统和方法 |
US11620772B2 (en) | 2016-09-01 | 2023-04-04 | The General Hospital Corporation | System and method for automated transform by manifold approximation |
CN106970343B (zh) * | 2017-04-11 | 2019-12-27 | 深圳先进技术研究院 | 一种磁共振成像方法及装置 |
CN106970343A (zh) * | 2017-04-11 | 2017-07-21 | 深圳先进技术研究院 | 一种磁共振成像方法及装置 |
CN107064845B (zh) * | 2017-06-06 | 2019-07-30 | 深圳先进技术研究院 | 基于深度卷积网的一维部分傅里叶并行磁共振成像方法 |
CN107064845A (zh) * | 2017-06-06 | 2017-08-18 | 深圳先进技术研究院 | 基于深度卷积网的一维部分傅里叶并行磁共振成像方法 |
US10989779B2 (en) | 2017-09-29 | 2021-04-27 | Yonsei University, University - Industry Foundation (UIF) | Apparatus and method for reconstructing magnetic resonance image using learning, and under-sampling apparatus method and recording medium thereof |
CN107993271A (zh) * | 2017-12-26 | 2018-05-04 | 上海交通大学 | 一种磁共振动态成像采样方法和图像重建方法 |
CN108535675B (zh) * | 2018-04-08 | 2020-12-04 | 朱高杰 | 一种基于深度学习和数据自洽的磁共振多通道重建方法 |
CN108535675A (zh) * | 2018-04-08 | 2018-09-14 | 朱高杰 | 一种基于深度学习和数据自洽的磁共振多通道重建方法 |
CN108335339A (zh) * | 2018-04-08 | 2018-07-27 | 朱高杰 | 一种基于深度学习和凸集投影的磁共振重建方法 |
CN110992433B (zh) * | 2018-10-02 | 2024-02-27 | 通用电气公司 | 用于利用空数据一致性来稀疏图像重建的系统和方法 |
CN110992433A (zh) * | 2018-10-02 | 2020-04-10 | 通用电气公司 | 用于利用空数据一致性来稀疏图像重建的系统和方法 |
US11175365B2 (en) | 2018-10-02 | 2021-11-16 | General Electric Company | System and method for sparse image reconstruction utilizing null data consistency |
WO2020135014A1 (zh) * | 2018-12-27 | 2020-07-02 | 深圳先进技术研究院 | 医学成像模型的建立方法、装置、设备及存储介质 |
CN111513716B (zh) * | 2019-02-05 | 2023-05-23 | 通用电气精准医疗有限责任公司 | 使用扩展灵敏度模型和深度神经网络进行磁共振图像重建的方法和系统 |
CN111513716A (zh) * | 2019-02-05 | 2020-08-11 | 通用电气精准医疗有限责任公司 | 使用扩展灵敏度模型和深度神经网络进行磁共振图像重建的方法和系统 |
CN111624540A (zh) * | 2019-02-27 | 2020-09-04 | 上海东软医疗科技有限公司 | 磁共振成像方法和装置 |
CN111624540B (zh) * | 2019-02-27 | 2022-11-22 | 上海东软医疗科技有限公司 | 磁共振成像方法和装置 |
CN110146834A (zh) * | 2019-05-16 | 2019-08-20 | 山东大学齐鲁医院 | 一种磁共振成像方法和装置 |
CN110664378A (zh) * | 2019-10-28 | 2020-01-10 | 中国科学院深圳先进技术研究院 | 磁共振成像方法、装置、系统及存储介质 |
CN110664378B (zh) * | 2019-10-28 | 2022-05-24 | 中国科学院深圳先进技术研究院 | 磁共振成像方法、装置、系统及存储介质 |
CN110988763B (zh) * | 2019-11-05 | 2021-12-17 | 深圳先进技术研究院 | 磁共振成像方法、装置、服务器及存储介质 |
CN110988763A (zh) * | 2019-11-05 | 2020-04-10 | 深圳先进技术研究院 | 磁共振成像方法、装置、服务器及存储介质 |
CN111681297A (zh) * | 2020-05-28 | 2020-09-18 | 上海联影智能医疗科技有限公司 | 图像重建方法、计算机设备和存储介质 |
CN111681297B (zh) * | 2020-05-28 | 2024-02-27 | 上海联影智能医疗科技有限公司 | 图像重建方法、计算机设备和存储介质 |
CN114325524B (zh) * | 2020-09-29 | 2023-09-01 | 上海联影医疗科技股份有限公司 | 磁共振图像重建方法、装置、系统及存储介质 |
CN114325524A (zh) * | 2020-09-29 | 2022-04-12 | 上海联影医疗科技股份有限公司 | 磁共振图像重建方法、装置、系统及存储介质 |
CN113933773A (zh) * | 2021-09-30 | 2022-01-14 | 深圳先进技术研究院 | 基于深度学习的磁共振成像方法、系统、终端及存储介质 |
WO2023050249A1 (zh) * | 2021-09-30 | 2023-04-06 | 深圳先进技术研究院 | 基于深度学习的磁共振成像方法、系统、终端及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103679654B (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103679654A (zh) | 磁共振成像方法和系统 | |
US10902651B2 (en) | Systems and methods for magnetic resonance image reconstruction | |
Tezcan et al. | MR image reconstruction using deep density priors | |
Akçakaya et al. | Compressed sensing with wavelet domain dependencies for coronary MRI: a retrospective study | |
CN105005012B (zh) | 基于压缩感知的腹部器官动态对比增强磁共振成像方法 | |
US10996306B2 (en) | MRI system and method using neural network for detection of patient motion | |
CN105518477B (zh) | 使用空间自适应正则化以用于图像重建的mri | |
CN110133556B (zh) | 一种磁共振图像处理方法、装置、设备及存储介质 | |
CN106997034A (zh) | 基于以高斯模型为实例整合重建的磁共振扩散成像方法 | |
CN109214989A (zh) | 基于多方向特征预测先验的单幅图像超分辨率重建方法 | |
EP2773985A1 (en) | Method for calibration-free locally low-rank encouraging reconstruction of magnetic resonance images | |
CN116402865B (zh) | 一种利用扩散模型的多模态影像配准方法、装置和介质 | |
US11776171B2 (en) | Systems and methods for magnetic resonance image reconstruction | |
CN104574456A (zh) | 一种基于图正则化稀疏编码的磁共振超欠采样k数据成像方法 | |
CN105678822A (zh) | 一种基于Split Bregman迭代的三正则磁共振图像重构方法 | |
CN103505208B (zh) | 一种磁共振成像方法 | |
CN106821381A (zh) | 一种血管成像方法 | |
Cheng et al. | Compressed sensing: From research to clinical practice with data-driven learning | |
Makhijani et al. | Accelerated 3D MERGE carotid imaging using compressed sensing with a hidden Markov tree model | |
Vedula et al. | 3D FLAT: Feasible learned acquisition trajectories for accelerated MRI | |
Awate et al. | Compressed sensing HARDI via rotation-invariant concise dictionaries, flexible K-space undersampling, and multiscale spatial regularity | |
FR2983999A1 (fr) | Procede de reconstruction d'un signal en imagerie medicale a partir de mesures experimentales perturbees, et dispositif d'imagerie medicale mettant en oeuvre ce procede | |
CN104504657A (zh) | 磁共振弥散张量去噪方法和装置 | |
CN104360295A (zh) | 基于字典学习的磁共振并行成像方法和装置 | |
Alush-Aben et al. | 3D FLAT: feasible learned acquisition trajectories for accelerated MRI |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |