一种基于压缩传感的快速高角分辨率扩散成像方法技术

本发明专利技术涉及一种基于压缩传感的快速高角分辨率扩散成像方法，包括：根据球面螺旋分布方式确定扩散敏感梯度磁场分布方向；为确定的每一个扩散敏感梯度磁场分布方向制定欠采样编码策略；采用全采样方式，采集第一个扩散敏感梯度磁场分布方向的一组k空间信号作为参考像；将其余扩散敏感梯度磁场分布方向及其欠采样编码策略输入磁共振扫描仪，依次采集其k空间信号；将得到的各组k空间信号分别与参考像k空间信号作差，并采用压缩传感方法重建得到差值像，再将差值像叠加到参考像上。本发明专利技术利用螺旋分布法、欠采样策略及HARDI数据特点，采集少量k空间信号就能很好地重建出未知图像，提高了采样及图像重建速度。

【技术实现步骤摘要】

本专利技术涉及磁共振成像
，特别是。
技术介绍
以下是本领域中的一些名词解释，其使用范围仅限于本专利技术k空间信号在磁共振成像过程，磁共振谱仪采集到的原始信号不同于图像像素所在的空间，而称为k空间，因此对应的原始信号称为k空间信号。相位编码梯度磁场为了对空间位置进行编码，就需要在不同的位置施加不同的磁场强度，这就要用到梯度磁场。当梯度的方向指向空间坐标系的y方向时，称这个方向为相位方向，对应的梯度磁场也就称为相位编码梯度磁场，简称相位梯度。扩散敏感梯度磁场在扩散张量成像中，除了通常施加的梯度磁场外，还要施加另外一种梯度磁场，这种梯度磁场能够提取出人体组织内的水分子扩散信息(大小，方向等)，因此称这种梯度磁场为扩散敏感梯度磁场。表观扩散系数是一种在宏观上度量水分子扩散速度的指标，这一指标是在观察范围内所有水分子扩散速度大小的平均，与扩散方向无关。扩散敏感梯度磁场的b值是衡量扩散敏感梯度磁场强弱的一个指标，由扩散敏感梯度磁场的强度和持续时间共同决定，近似相当于扩散敏感梯度磁场强度对时间的积分。现有技术中的磁共振成像技术主要有三个方面的发展一、磁共振成像技术和压缩传感技术磁共振成像(Magnetic Resonance Imaging, MRI)技术是目前广泛应用的医疗检查技术。其工作原理是，人体内的水分子会在核磁共振效应下发出电磁信号，经过特定的梯度磁场标记后，磁共振谱仪采集这些信号(即k空间信号)，然后经过傅里叶变换，重建出组织影像。这种重建方式需要采集的数据量与重建后的数据量是相等的，因此要想获得更高分辨率的图像，就要采集更多的数据。压缩传感(Compressed sensing, CS)理论使用了一种非线性重建方式,只需要采集少量k空间信号，就可得到与傅里叶重建几乎相同的图像质量。k空间信号通常是一行一行地采集的，每一行数据对应采集中的一个相位编码梯度磁场(简称相位梯度)，因此可以用相位梯度数目来衡量k空间信号的数目，相位梯度数目越少，k空间中采集的数据就越少。傅里叶重建需要对k空间全采样，CS重建则可以对k空间欠采样。通常情况下，CS重建只需采集1/10的k空间信号，就能很好的重建出图像。然而，过少的k空间信号，即便在CS重建下，也会使图像损失大量信息。但是如果两幅图像的差别比较小，并且已知其中一幅图像，那么即便只采集了很少的k空间信号，仍然可以很好的重建出未知图像。采集方式为采集非常少量的k空间信号，将这些信号与已知图像k空间中对应的信号作差，对得到的差值做CS重建，把重建得到的图像叠加到已知图像上，就会得到未知图像。使用这种方法可以使采集多组相似图像的过程大大缩短。而利用CS重建非常重要的步骤之一是制定欠采样策略，即随机选取相位梯度(一定要随机选取)。通常的做法是，在进行磁共振扫描前，把相位梯度选择好，待到扫描时，直接把选好的梯度输入到机器中，进行扫描。本专利技术即是基于这一思路。二、高角分辨率扩散成像及其与压缩传感的兼容性扩散张量成像(Diffusion Tensor Imaging, DTI)是磁共振扩散成像的进一步发展，在临床诊断及脑科学研究等方面得到广泛应用。它的原理是，在通常的磁共振扫描中增加扩散敏感梯度磁场，以使图像反映出组织内水分子的扩散信息。一次DTI成像需要在不同方向上施加扩散敏感梯度，每个方向下得到的图像称为一组扩散权重像(DiffusionWeight Image, DWI)。高角分辨率扩散成像(High Angular Resolution DiffusionImaging, HARDI)是在DTI的基础上施加更多扩散敏感梯度，得到的生理及病理信息远多于DTI，有重要临床意义及需求，但尚未应用于临床。主要原因是HARDI需要在数百甚至上千方向上采集DWI信号，导致其十分耗时，且极易受运动、畸变等因素影响。所以研究HARDI快速成像，解决其中的关键问题有重要的理论意义及应用前景。以平面自旋回波(Echo Planar Imaging, EPI)为代表的多种快速成像方法将磁共振成像的常规检查时间由数小时缩短到数分钟，但仍不能满足不断增长的临床需求。尤其是HARDI要求在尽量多的方向上测量水分子的表观扩散系数，在研究中甚至会多达上千方向，这必然导致扫描时间延长。让病人在漫长的扫描时间内保持静止是极困难的，而且还有呼吸、心跳及不自主运动等因素的影响。MRI数据的稀疏性使其适用基于CS理论的数据采集及图像重建方法，根据HARDI数据的特点改进当前单层欠采样方法可以大幅提高采集及图像重建速度。以脑部为例，HARDI数据中包括大量不同扩散敏感梯度磁场方向下得到的DWI图像，理论分析可知，其中只有白质部分的信号会有差异，水分子在脑中其他部位(如脑脊液，灰质)中的扩散是各向同性的，所以不同扩散敏感梯度磁场方向下这些部位的信号相同，由此可见不同扩散敏感梯度磁场方向的DWI图像之间差异很小。。因此HARDI符合使用CS的条件，可以设计一种方法，进一步把HARDI的扫描时间缩短。这也是本专利技术的重要思路之一。三、扩散敏感梯度磁场分布方式对于大量的梯度方向，如何合理的分布是一个重要的问题。目前使用的分布方式有正十二面体分布法和正二十面体分布法等，以正二十面体分布法为例，以正二十面体的中心为原点，由原点与正二十面体的顶点的连线确定方向，更多方向的分布可以通过进一步划分正二十面体的面来实现。这种方法能够均匀的分布大量方向，但是梯度对方向的数目有较严格的限制，不能把任意数量的梯度分布在各个方向上。因此，寻求一种通用的扩散敏感梯度磁场分布方法，也是磁共振成像需要解决的一个重要问题。
技术实现思路
本专利技术所要解决的技术问题是提供，用于解决现有技术中磁共振成像存在的图像采集质量低、采集速度较慢、扫描时间长、扩散敏感梯度磁场分布方法适应性低等问题。本专利技术解决上述技术问题的技术方案如下，包括步骤I，按球面螺旋分布方式在单位半球面上确定扩散敏感梯度磁场分布方向，并把半数的扩散敏感梯度磁场分布方向相对球心对称翻转到另一半球面，得到单位球面上所有的扩散敏感梯度磁场分布方向；步骤2，将所有扩散敏感梯度磁场分布方向按螺旋线前进方向编号，并为除第一个扩散敏感梯度磁场分布方向外的其余方向制定相应的欠采样编码策略；步骤3，采用全采样方式，采集第一个扩散敏感梯度磁场分布方向的一组k空间信号，对这组k空间信号作傅里叶重建，将重建后的结果作为参考像；步骤4，应用步骤2中的欠采样编码策略，依次采集除第一个扩散敏感梯度磁场分布方向外的其余方向的k空间信号；步骤5，将步骤4得到的一组欠采样的k空间信号与步骤3中的k空间信号作差，得到新的k空间欠采样信号；步骤6，采用压缩传感方法对步骤5得到的k空间欠采样信号进行重建，重建后得到该扩散敏感梯度磁场分布方向的差值像；步骤7，将步骤6得到的差值像叠加到步骤3得到的参考像上，得到该扩散敏感梯度磁场分布方向下的扩散权重像；步骤8，重复步骤5至步骤7，得到其余扩散敏感梯度磁场分布方向下的扩散权重像。在上述技术方案的基础上，本专利技术还可以做如下改进进一步，所述步骤I具体包括步骤11，选取单位球面的上顶点作为球面螺旋线初始点；步骤12，计算球面螺旋线在上半球的长度随环绕圈数变化的关系式；步骤13，计算在本文档来自技高网...

【技术保护点】
一种基于压缩传感的快速高角分辨率扩散成像方法，其特征在于，包括：步骤1，按球面螺旋分布方式在单位半球面上确定扩散敏感梯度磁场分布方向，并把半数的扩散敏感梯度磁场分布方向相对球心对称翻转到另一半球面，得到单位球面上所有的扩散敏感梯度磁场分布方向；步骤2，将所有扩散敏感梯度磁场分布方向按螺旋线前进方向编号，并为除第一个扩散敏感梯度磁场分布方向外的其余方向制定相应的欠采样编码策略；步骤3，采用全采样方式，采集第一个扩散敏感梯度磁场分布方向的一组k空间信号，对这组k空间信号作傅里叶重建，将重建后的结果作为参考像；步骤4，应用步骤2中的欠采样编码策略，依次采集除第一个扩散敏感梯度磁场分布方向外的其余方向的k空间信号；步骤5，将步骤4得到的一组欠采样的k空间信号与步骤3中的k空间信号作差，得到新的k空间欠采样信号；步骤6，采用压缩传感方法对步骤5得到的k空间欠采样信号进行重建，重建后得到该扩散敏感梯度磁场分布方向的差值像；步骤7，将步骤6得到的差值像叠加到步骤3得到的参考像上，得到该扩散敏感梯度磁场分布方向下的扩散权重像；步骤8，重复步骤5至步骤7，得到其余扩散敏感梯度磁场分布方向下的扩散权重像。

【技术特征摘要】
1.一种基于压缩传感的快速高角分辨率扩散成像方法，其特征在于，包括 步骤1，按球面螺旋分布方式在单位半球面上确定扩散敏感梯度磁场分布方向，并把半数的扩散敏感梯度磁场分布方向相对球心对称翻转到另一半球面，得到单位球面上所有的扩散敏感梯度磁场分布方向； 步骤2，将所有扩散敏感梯度磁场分布方向按螺旋线前进方向编号，并为除第一个扩散敏感梯度磁场分布方向外的其余方向制定相应的欠采样编码策略； 步骤3，采用全采样方式，采集第一个扩散敏感梯度磁场分布方向的一组k空间信号，对这组k空间信号作傅里叶重建，将重建后的结果作为参考像； 步骤4，应用步骤2中的欠采样编码策略，依次采集除第一个扩散敏感梯度磁场分布方向外的其余方向的k空间信号； 步骤5，将步骤4得到的一组欠采样的k空间信号与步骤3中的k空间信号作差，得到新的k空间欠采样信号； 步骤6，采用压缩传感方法对步骤5得到的k空间欠采样信号进行重建，重建后得到该扩散敏感梯度磁场分布方向的差值像； 步骤7，将步骤6得到的差值像叠加到步骤3得到的参考像上，得到该扩散敏感梯度磁场分布方向下的扩散权重像； 步骤8，重复步骤5至步骤7，得到其余扩散敏感梯度磁场分布方向下的扩散权重像。2.根据权利要求1所述的快速高角分辨率扩散成像方法，其特征在于，所述步骤I具体包括 步骤U，选取单位球面的上顶点作为球面螺旋线初始点； 步骤12，计算球面螺旋线在上半球的长度随环绕圈数变化的关系式； 步骤13，计算在指定扩散敏感梯度磁场分布方向数目下，以球面螺旋线间距为两相邻方向的球面距离，排列所有方向所需的螺旋线总长度随球面螺旋线环绕圈数变化的关系式； 步骤14，由步骤12及步骤13的两关系式确定球面螺旋线的环绕圈数，并构造上半球的球面螺旋线...

【专利技术属性】
技术研发人员：高嵩，王飞，王浩宇，
申请(专利权)人：北京大学，
类型：发明
国别省市：