专利名称：快速磁共振成像方法及系统的制作方法
快速磁共振成像方法及系统
技术领域：
本发明涉及磁共振技术，特别是涉及一种快速磁共振成像方法及系统。背景技术：
随着磁共振技术的发展，磁共振成像已经成为临床医学检查的重要手段，具有无辐射危害、多方位和多参数成像等优点，对软组织的检查非常敏感，不仅能够显示人体解剖结构的形态信息，而且还能反映人体组织的某些生理化信息。但是，磁共振成像的成像速度较慢，且在成像过程中被检者身体中的生理性运动将会使影像模糊，无法满足运动器官的动态成像要求。在传统的快速成像方法中，例如压缩感知成像方法，通过磁共振采集少量的数据重建出高分辨的图像，需要通过求解最小1。范数问题，其中，Itl为向量中非零元素的个数。求解最小Itl范数问题需要穷举信号中非零值的所有稀疏度中可能的排列，无法直接求解。而其它求得次最优解的方法，例如，最小Itl范数法、贪婪迭代匹配追踪算法等，在最小Itl范数法中的线性规划方法能够重建出磁共振图像，但是需要以较高的时间为代价，并不适用于动态成像，而贪婪迭代匹配追踪算法虽然能够较快得到最优解，但是是建立在稀疏度已知的基础上的，实际应用中信号的稀疏度往往是未知的，也不适用于高分辨率的动态成像。

发明内容基于此，有必要提供一种能提高动态成像中成像分辨率的快速磁共振成像方法。此外，还有必要提供一种能提高动态成像中成像分辨率的快速磁共振成像系统。一种快速磁共振成像方法，包括如下步骤进行变密度采样得到空间数据；对所述空间数据进行稀疏变换生成测量矩阵；通过正则化自适应匹配追踪算法对测量矩阵进行信号重构。优选地，所述进行变密度采样得到空间数据的步骤为根据能量分布确定变密度的随机欠采样轨迹；通过所述变密度的随机欠采样轨迹进行采样得到空间数据。优选地，所述对所述空间数据进行稀疏变换生成测量矩阵的步骤具体为运用小波变换对所述空间数据进行稀疏变换生成测量矩阵。优选地，所述通过正则化自适应匹配追踪算法对测量矩阵进行信号重构的步骤为通过余量筛选所述测量矩阵中的原子得到支撑集；通过最小二乘法进行信号逼近，并更新所述余量；判断所述更新后的余量与更新前的余量之间的绝对差值是否大于阶段转换阈值， 若否，则进一步判断所述更新余量的绝对值是否超出控制迭代次数，若是，则返回所述通过余量筛选所述测量矩阵中的原子得到支撑集的步骤，若是，则
停止迭代，并通过所述支撑集进行信号重建。优选地，所述通过余量筛选所述测量矩阵的步骤具体为通过余量与测量矩阵中各原子之间内积的绝对值计算得到相关系数；按照从大到小的顺序提取数量与转换阶段相匹配的相关系数，并将所述提取的相关系数对应的原子作为索引值存入索引集合中；在所述索引集合中通过正则化过程进行二次筛选得到支撑集。一种快速磁共振成像系统，包括采样模块，用于进行变密度采样得到空间数据；稀疏变换模块，用于对所述空间数据进行稀疏变换生成测量矩阵；重构模块，用于通过正则化自适应匹配追踪算法对测量矩阵进行信号重构。优选地，所述采样模块包括轨迹确定单元，用于根据能量分布确定变密度的随机欠采样轨迹；欠采样单元，用于通过所述变密度的随机欠采样轨迹进行采样得到空间数据。
优选地，所述稀疏变换模块还用于运用小波变换对所述空间数据进行稀疏变换生成测量矩阵。优选地，所述重构模块包括筛选单元，用于通过余量筛选所述测量矩阵中的原子得到支撑集；余量更新单元，用于通过最小二乘法进行信号逼近，并更新所述余量；迭代控制单元，用于判断所述更新后的余量与更新前的余量之间的绝对差值是否大于阶段转换阈值，若否，则进一步判断所述更新余量的绝对值是否超出控制迭代次数，若是，则通知所述筛选单元，若是，则通知所述信号重建单元；所述信号重建单元用于停止迭代，并通过所述支撑集进行信号重建。优选地，所述筛选单元包括计算单元，用于通过余量与测量矩阵中各原子之间内积的绝对值计算得到相关系数；提取单元，用于按照从大到小的顺序提取数量与转换阶段相匹配的相关系数，并将所述提取的相关系数对应的原子作为索引值存入索引集合中；正则化单元，用于在所述索引集合中通过正则化过程进行二次筛选得到支撑集。上述快速磁共振成像方法及系统，通过变密度采样得到少量的空间数据，并正则化自适应匹配追踪算法的作用下完成图像的高分辨率重建，变密度采样与正则化自适应匹配追踪算法的结合，使得成像速度大大提高，且成像质量也得到了提高，进而提高了动态成像中的成像分辨率。

图1为一个实施例中快速磁共振成像方法的流程图；图2为一个实施例中进行变密度采样得到空间数据的方法流程图；图3为一个实施例中通过正则化自适应匹配追踪算法对测量矩阵进行信号重构的方法流程图；图4为图3中通过余量筛选测量矩阵中的原子得到支撑集的方法流程图5为一个实施例中快速磁共振成像系统的结构示意图；图6为一个实施例中采样模块的结构示意图；图7为一个实施例中重构模块的结构示意图；图8为图7中筛选单元的结构示意图。
具体实施方式在一个实施例中，如图1所示，一种快速磁共振成像方法，包括如下步骤步骤S10，进行变密度采样得到空间数据。本实施例中，在K空间中进行变密度采样，以得到较为少量的空间数据，降低采样所耗费的时间。在一个实施例中，如图2所示，上述步骤SlO的具体过程包括步骤S110，根据能量分布确定变密度的随机欠采样轨迹。本实施例中，K空间中能量分布是各不相同的，即可在K空间的中间部分进行稠密采样，在外围区域进行稀疏采样，进而生成变密度的随机欠采样轨迹。步骤S130，通过变密度的随机欠采样轨迹进行采样得到空间数据。本实施例中，在基于笛卡尔坐标的变密度随机欠采样轨迹可对频率编码方向上的数据进行全采样，对相位编码方向和层编码方向上进行变密度的随机欠采样，也可以是在频率编码方向上的数据进行全采样，对相位编码方向或层编码方向进行变密度的随机欠采样，然后将得到的数据坐标通过插值的方法网格化到笛卡尔坐标系中，即，采集到的信号编码后转化成数据，并将数据放入笛卡尔坐标系这一二维数据空间中。步骤S30，对空间数据进行稀疏变换生成测量矩阵。本实施例中，压缩感知理论的前提条件是具有稀疏性或可压缩性，为简化模型，只考虑长度为N的空间数据X，记为χ (η)。信号数据χ用一组基Ψτ= [Ψι，Ψ2，…，Ψη]的线性组合表示，其中，Ψ的转置，则
权利要求
1.一种快速磁共振成像方法，包括如下步骤 进行变密度采样得到空间数据；对所述空间数据进行稀疏变换生成测量矩阵； 通过正则化自适应匹配追踪算法对测量矩阵进行信号重构。
2.根据权利要求1所述的快速磁共振成像方法，其特征在于，所述进行变密度采样得到空间数据的步骤为根据能量分布确定变密度的随机欠采样轨迹； 通过所述变密度的随机欠采样轨迹进行采样得到空间数据。
3.根据权利要求1所述的快速磁共振成像方法，其特征在于，所述对所述空间数据进行稀疏变换生成测量矩阵的步骤具体为运用小波变换对所述空间数据进行稀疏变换生成测量矩阵。
4.根据权利要求1所述的快速磁共振成像方法，其特征在于，所述通过正则化自适应匹配追踪算法对测量矩阵进行信号重构的步骤为通过余量筛选所述测量矩阵中的原子得到支撑集； 通过最小二乘法进行信号逼近，并更新所述余量；判断所述更新后的余量与更新前的余量之间的绝对差值是否大于阶段转换阈值，若否，则进一步判断所述更新余量的绝对值是否超出控制迭代次数，若是，则返回所述通过余量筛选所述测量矩阵中的原子得到支撑集的步骤，若是，则停止迭代，并通过所述支撑集进行信号重建。
5.根据权利要求4所述的快速磁共振成像方法，其特征在于，所述通过余量筛选所述测量矩阵的步骤具体为通过余量与测量矩阵中各原子之间内积的绝对值计算得到相关系数； 按照从大到小的顺序提取数量与转换阶段相匹配的相关系数，并将所述提取的相关系数对应的原子作为索引值存入索引集合中；在所述索引集合中通过正则化过程进行二次筛选得到支撑集。
6.一种快速磁共振成像系统，其特征在于，包括 采样模块，用于进行变密度采样得到空间数据；稀疏变换模块，用于对所述空间数据进行稀疏变换生成测量矩阵； 重构模块，用于通过正则化自适应匹配追踪算法对测量矩阵进行信号重构。
7.根据权利要求6所述的快速磁共振成像系统，其特征在于，所述采样模块包括 轨迹确定单元，用于根据能量分布确定变密度的随机欠采样轨迹；欠采样单元，用于通过所述变密度的随机欠采样轨迹进行采样得到空间数据。
8.根据权利要求6所述的快速磁共振成像系统，其特征在于，所述稀疏变换模块还用于运用小波变换对所述空间数据进行稀疏变换生成测量矩阵。
9.根据权利要求6所述的快速磁共振成像系统，其特征在于，所述重构模块包括 筛选单元，用于通过余量筛选所述测量矩阵中的原子得到支撑集；余量更新单元，用于通过最小二乘法进行信号逼近，并更新所述余量； 迭代控制单元，用于判断所述更新后的余量与更新前的余量之间的绝对差值是否大于阶段转换阈值，若否，则进一步判断所述更新余量的绝对值是否超出控制迭代次数，若是，则通知所述筛选单元，若是，则通知所述信号重建单元；所述信号重建单元用于停止迭代，并通过所述支撑集进行信号重建。
10.根据权利要求9所述的快速磁共振成像系统，其特征在于，所述筛选单元包括 计算单元，用于通过余量与测量矩阵中各原子之间内积的绝对值计算得到相关系数； 提取单元，用于按照从大到小的顺序提取数量与转换阶段相匹配的相关系数，并将所述提取的相关系数对应的原子作为索引值存入索引集合中；正则化单元，用于在所述索引集合中通过正则化过程进行二次筛选得到支撑集。
全文摘要
一种快速磁共振成像方法，包括如下步骤进行变密度采样得到空间数据；对所述空间数据进行稀疏变换生成测量矩阵；通过正则化自适应匹配追踪算法对测量矩阵进行信号重构。上述快速磁共振成像方法及系统，通过变密度采样得到少量的空间数据，并正则化自适应匹配追踪算法的作用下完成图像的高分辨率重建，变密度采样与正则化自适应匹配追踪算法的结合，使得成像速度大大提高，且成像质量也得到了提高，进而提高了动态成像中的成像分辨率。
文档编号G01R33/56GK102565737SQ201210013858
公开日2012年7月11日 申请日期2012年1月17日 优先权日2011年12月12日
发明者冯翔, 刘新, 史彩云, 谢国喜, 邱本胜 申请人:中国科学院深圳先进技术研究院