一种磁共振快速成像方法及系统技术方案

本发明专利技术提出了一种磁共振快速成像方法及系统，其中，该方法包括：获取欠采样的磁共振数据；根据欠采样的磁共振数据，利用多个卷积核函数，计算多种待求的磁共振图像特征所对应的K空间数据；根据多种待求的磁共振图像特征所对应的K空间数据，利用压缩感知技术重建磁共振图像特征，获取重建的图像特征；将重建的图像特征进行合并，生成重建的磁共振图像。本发明专利技术提出的磁共振快速成像方法及系统首先利用压缩感知技术重建多个卷积核函数所对应的图像特征，这些图像特征或者更稀疏，或者具有更高的信噪比，因而更利于运用压缩感知技术进行重建；最后，图像特征的合并能够进一步地抑制重建的图像特征中的噪声和伪影，进而实现高质量的图像重建。

Magnetic resonance rapid imaging method and system

The invention provides a fast magnetic resonance imaging method and system, wherein, the method comprises: acquiring less magnetic resonance data sampling; according to the MRI data sampling, using multiple convolution kernel, K spatial data corresponding to magnetic resonance image characteristics calculation to solve the spatial data based on K; the corresponding image features a variety of magnetic resonance to solve the reconstruction using sensing magnetic resonance image feature compression, image feature reconstruction; image feature reconstruction combined with magnetic resonance image generation reconstruction. Fast magnetic resonance imaging method and system provided by the invention firstly uses compressed sensing reconstruction image features corresponding to multiple convolution kernel, these image features or more sparse, or has a higher signal-to-noise ratio, and thus more conducive to the use of reconstruction of compressed sensing technology; finally, combined with image features to further to suppress the noise in the reconstructed image features and artifacts, and realize the image reconstruction of high quality.

【技术实现步骤摘要】
一种磁共振快速成像方法及系统
本专利技术涉及磁共振成像
，特别涉及一种磁共振快速成像方法及系统。
技术介绍
近年来，基于压缩感知技术的磁共振成像方法是快速磁共振成像领域的热点问题。压缩感知技术主要是利用信号的稀疏性以实现欠采样重建。磁共振成像的物理机制可以表示为：d＝Fρ+n(公式1)其中，d表示在磁共振仪上所采集的磁共振数据，F表示傅立叶编码矩阵，ρ为待求的磁共振图像，n通常假定为复高斯白噪声。压缩感知理论认为，如果信号在某一变换域内是稀疏的，那么通过一种与该变换域非相干的采样模式，就可以从少量的采样数据中精确的重建出真实的信号。并且，信号越稀疏，需要的采样数据就越少。基于压缩感知技术的磁共振成像方法可以表示为：其中，表示基于压缩感知技术重建的磁共振成像，ρ表示待求的磁共振图像，y为k空间的欠采样数据，Fu为欠采样的傅立叶编码矩阵，Ψ为稀疏变换矩阵(例如小波变换)，λ为正则化系数，||·||1表示求L1范数，||·||2表示求L2范数。然而，医学磁共振图像往往只是高度可压缩的，而并非严格稀疏的。这种稀疏性不足以及采样中的噪声会导致重建的图像中产生严重的伪影，降低了图像的成像质量，从而限制了压缩知感技术在加速磁共振成像中的实际应用。
技术实现思路
基于压缩感知的磁共振成像技术，在图像稀疏性不足以及有噪的情况下会产生比较严重的图像伪影。针对前述问题，本专利技术利用一系列卷积核函数提取图像的特征，这些图像特征需要比磁共振图像更稀疏或者具有更高的信噪比。具体的，该方法先利用压缩感知技术重建这些图像特征，再将这些重建得到的图像特征合并，以生成重建的图像。为达到上述目的，本专利技术提出了一种磁共振快速成像方法，包括：获取欠采样的磁共振数据；根据所述欠采样的磁共振数据，利用多个卷积核函数，计算多种待求的磁共振图像特征所对应的K空间数据；根据所述多种待求的磁共振图像特征所对应的K空间数据，利用压缩感知技术重建磁共振图像特征，获取重建的图像特征；将所述重建的图像特征进行合并，生成重建的磁共振图像。为达到上述目的，本专利技术还提出了一种磁共振快速成像系统，包括：磁共振数据获取模块，用于获取欠采样的磁共振数据；K空间数据计算模块，用于根据所述欠采样的磁共振数据，利用多个卷积核函数，计算多种待求的磁共振图像特征对应的K空间数据；图像特征重建模块，用于根据所述多种待求的磁共振图像特征所对应的K空间数据，利用压缩感知技术重建磁共振图像特征，获取重建的图像特征；重建图像生成模块，用于将所述重建的图像特征进行合并，生成重建的磁共振图像。本专利技术提出的磁共振快速成像方法及系统首先利用压缩感知技术重建多个卷积核函数所对应的图像特征，这些图像特征或者更稀疏，或者具有更高的信噪比，因而更利于运用压缩感知技术进行重建；最后，图像特征的合并能够进一步地抑制重建的图像特征中的噪声和伪影，进而实现高质量的图像重建。附图说明此处所说明的附图用来提供对本专利技术的进一步理解，构成本申请的一部分，并不构成对本专利技术的限定。在附图中：图1为本专利技术一实施例的磁共振快速成像方法流程图。图2为本专利技术一实施例的磁共振快速成像系统的结构示意图。具体实施方式以下配合图式及本专利技术的较佳实施例，进一步阐述本专利技术为达成预定专利技术目的所采取的技术手段。本专利技术提出了一种磁共振快速成像方法及系统，相比较现有的利用标准的压缩感知的磁共振成像技术的成像效果，本专利技术首先利用压缩感知技术重建多个卷积核函数所对应的图像特征；再将重建的图像特征进行合并，得到重建的磁共振图像，进而实现高质量的图像重建。图1为本专利技术一实施例的磁共振快速成像方法流程图。如图1所示，该方法包括：步骤1，获取欠采样的磁共振数据；步骤2，根据欠采样的磁共振数据，利用多个卷积核函数，计算多种待求的磁共振图像特征所对应的K空间数据；步骤3，根据多种待求的磁共振图像特征所对应的K空间数据，利用压缩感知技术重建磁共振图像特征，获取重建的图像特征；步骤4，将重建的图像特征进行合并，生成重建的磁共振图像。进一步的，在步骤2中，多种待求的磁共振图像特征所对应的K空间数据为：y1＝F{K1}⊙yy2＝F{K2}⊙yyp＝F{Kp}⊙y其中，y1、y2、…、yp表示多种待求的磁共振图像特征所对应的K空间数据，K1、K2、…、Kp表示多个卷积核函数，⊙为点乘运算符，F{·}表示傅立叶变换，y为欠采样的磁共振数据，p为卷积核函数的个数。其中，卷积核函数可以为边缘提取核函数或高斯核函数。利用边缘提取核函数重建的图像特征可以更稀疏，利用高斯核函数重建的图像特征具有更高的信噪比。因此，相比较
技术介绍
公式2的直接重建磁共振图像的方法而言，本专利技术利用压缩感知技术将能更精确的重建出这些图像特征。在步骤3中，基于压缩感知技术重建磁共振图像特征表示为：其中，表示重建的图像特征，ρ1、ρ2、…、ρp表示待求的磁共振图像特征，Fu为欠采样的傅立叶编码矩阵，λ1、λ2、…、λp为正则化系数，Ψ表示稀疏变换矩阵，||·||1表示求L1范数，||·||2表示求L2范数；其中，磁共振图像特征与磁共振图像的关系表示为：ρ1＝K1*ρρ2＝K2*ρρp＝Kp*ρ其中，*为卷积运算符，ρ为待求的磁共振图像。在步骤4中，主要是将步骤3中的重建图像特征进行合并，生成重建的磁共振图像。其中，本实施例中利用的合并方法是利用正则化约束的方法，即完整的重建图像由以下方式获得：其中，为重建的磁共振图像，函数R(·)用于约束磁共振图像的稀疏性，γ1为控制待求的磁共振图像稀疏性的正则化参数，γ2为控制待求的磁共振图像特征与重建的图像特征之间相似性的正则化参数，α1、α2、…、αp表示各个磁共振图像特征的权重。函数R(·)用于约束磁共振图像的稀疏性，其中，R(ρ)＝||Ψρ||1。在其他实施例中，步骤4中利用的合并方法并不仅限于此，上述实施例中给出的只是最优化的方法之一。通过上述步骤1-4便可以得到重建图像，该重建图像的质量要远优于直接利用压缩感知技术(公式2)获得的图像。基于同一专利技术构思，本专利技术实施例中还提供了一种磁共振快速成像系统，如下面的实施例。由于该系统解决问题的原理与上述方法相似，因此该系统的实施可以参见上述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图2为本专利技术一实施例的磁共振快速成像系统的结构示意图。如图2所示，该系统包括：磁共振数据获取模块1，用于获取欠采样的磁共振数据；K空间数据计算模块2，用于根据欠采样的磁共振数据，利用多个卷积核函数，计算多种待求的磁共振图像特征所对应的K空间数据；图像特征重建模块3，用于根据多种待求的磁共振图像特征所对应的K空间数据，利用压缩感知技术重建磁共振图像特征，获取重建的图像特征；重建图像生成模块4，用于将重建的图像特征进行合并，生成重建的磁共振图像。K空间数据计算模块2中，多种待求的磁共振图像特征所对应的K空间数据为：y1＝F{K1}⊙yy2＝F{K2}⊙yyp＝F{Kp}⊙y其中，y1、y2、…、yp表示多种待求的磁共振图像特征所对应的K空间数据，K1、K2、…、Kp表示多个卷积核函数，本文档来自技高网...

【技术保护点】
一种磁共振快速成像方法，其特征在于，包括：获取欠采样的磁共振数据；根据所述欠采样的磁共振数据，利用多个卷积核函数，计算多种待求的磁共振图像特征所对应的K空间数据；根据所述多种待求的磁共振图像特征所对应的K空间数据，利用压缩感知技术重建磁共振图像特征，获取重建的图像特征；将所述重建的图像特征进行合并，生成重建的磁共振图像。

【技术特征摘要】
1.一种磁共振快速成像方法，其特征在于，包括：获取欠采样的磁共振数据；根据所述欠采样的磁共振数据，利用多个卷积核函数，计算多种待求的磁共振图像特征所对应的K空间数据；根据所述多种待求的磁共振图像特征所对应的K空间数据，利用压缩感知技术重建磁共振图像特征，获取重建的图像特征；将所述重建的图像特征进行合并，生成重建的磁共振图像。2.如权利要求1所述的磁共振快速成像方法，其特征在于，根据所述欠采样的磁共振数据，利用多个卷积核函数，计算得到多种待求的磁共振图像特征所对应的K空间数据，包括：多种待求的磁共振图像特征所对应的K空间数据为：y1＝F{K1}⊙yy2＝F{K2}⊙yyp＝F{Kp}⊙y其中，y1、y2、…、yp表示多种待求的磁共振图像特征所对应的K空间数据，K1、K2、…、Kp表示多个卷积核函数，⊙为点乘运算符，F{·}表示傅立叶变换，y为欠采样的磁共振数据，p为卷积核函数的个数。3.如权利要求2所述的磁共振快速成像方法，其特征在于，所述卷积核函数为边缘提取核函数或高斯核函数。4.如权利要求2所述的磁共振快速成像方法，其特征在于，根据所述多种待求的磁共振图像特征所对应的K空间数据，利用压缩感知技术重建磁共振图像特征，获取重建的图像特征，包括：基于压缩感知技术重建磁共振图像特征表示为：其中，表示重建的图像特征，ρ1、ρ2、…、ρp表示待求的磁共振图像特征，Fu为欠采样的傅立叶编码矩阵，λ1、λ2、…、λp为正则化系数，Ψ表示稀疏变换矩阵，||·||1表示求L1范数，||·||2表示求L2范数；其中，磁共振图像特征与磁共振图像的关系表示为：ρ1＝K1*ρρ2＝K2*ρρp＝Kp*ρ其中，*为卷积运算符，ρ为待求的磁共振图像。5.如权利要求4所述的磁共振快速成像方法，其特征在于，将所述重建的图像特征进行合并，生成重建的磁共振图像，包括：重建的磁共振图像由以下方式获得：其中，为重建的磁共振图像，函数R(·)用于约束磁共振图像的稀疏性，γ1为控制待求的磁共振图像稀疏性的正则化参数，γ2为控制待求的磁共振图像特征与重建的图像特征之间相似性的正则化参数，α1、α2、…、αp分别表示各个所述磁共振图像特征的权重。6.如权利要求5所述的磁共振快速成像方法，其特征在于，所述函数R(·)用于...

【专利技术属性】
技术研发人员：彭玺，梁栋，王珊珊，安一硕，刘新，郑海荣，
申请(专利权)人：中国科学院深圳先进技术研究院，
类型：发明
国别省市：广东,44