磁共振成像方法技术

本发明专利技术提供一种磁共振成像方法，该方法基于图像域和K空间域的迭代计算获得图像。根据本发明专利技术的方法包括以下步骤：S1:沿着指定编码方向对K空间进行有规律的欠采样；和S2：基于图像域和K空间域通过迭代计算获得未采集的K空间数据，以获得磁共振图像。本发明专利技术的该磁共振成像方法不仅能够用多通道来实现，也可以用单通道来实现，仅采集部分欠采样数据，通过在图像域和K空间域之间进行迭代计算实现快速的磁共振成像，使得图像的采集时间大大缩短，并可获得较高的成像质量。本发明专利技术的磁共振成像方法可用于缩短扫描时间，或在同样扫描时间的情况下改善信噪比。

【技术实现步骤摘要】
磁共振成像方法
本专利技术涉及磁共振成像技术，更具体地，涉及一种基于图像域和K空间域迭代的快速磁共振成像方法。
技术介绍
扫描速度是磁共振成像发展的瓶颈。如何有效地提高磁共振的成像速度一直是研发方面的重点项目。直接提高梯度线圈性能或主磁场强度，需要更高的硬件设备成本，并且受到物理、工程、生理、安全等诸多方面的限制。因此，改进扫描或重建算法等途径越来越受到国内外学者的重视。近年来，发展较活跃的快速成像方法主要包括压缩感知、并行成像以及SPEED(SkippedPhaseEncodingandEdgeDeghosting)方法等。压缩感知成像技术突破了香农-奈奎斯特的采样定理关于采样频率必须高于2倍信号带宽的限制，成为近年来研究的热点。虽然压缩感知在理论上能够实现快速成像，但是还有许多实际问题有待解决，例如缺乏稳定的抗噪声稀疏分解算法、图像临床应用的限制及重建时间的考虑等。并行成像技术利用线圈的空间敏感度信息，由多个接收线圈同时进行独立采样和编码，节约了空间定位所需梯度编码数目，加快了成像速度。目前主要的并行成像技术有SMASH、GRAPPA及SENSE等。并行成像很大程度上有赖于线圈敏感度的恰当分布和准确估测，面临难以避免的信噪比损失，并且无法实现单通道的扫描加速。SPEED算法利用信号的稀疏特性，通过采集若干组有规律跳跃的相位编码数据，利用线性方程逐点解析重建出图像。传统的SPEED方法可以单通道实现磁共振的快速扫描。然而，传统的SPEED方法需要沿相位编码(PE)方向分别采集三组欠采样数据，这三组欠采样数据具有相同的PE间隔和不同的PE偏移量d，通过求解最小二乘误差问题来实现混叠图像的分离。由于多组数据的采集以及逐点的解析计算，在加速效率和成像质量方面仍有待改进。因此，需要一种更加有效实用的快速磁共振扫描成像方法。
技术实现思路
为了解决上述问题，本专利技术提供一种基于图像域和K空间域迭代的磁共振成像方法，即全息高效成像方法HEIGHT(HighlyEfficientImagingwithGlobalHarmonicTransformations)。根据本专利技术的磁共振成像方法包括以下步骤：S1:沿着指定编码方向对K空间进行有规律的欠采样；和S2：基于图像域和K空间域通过迭代计算获得未采集的K空间数据，以获得磁共振图像，其中，欠采样步骤S1包括以下步骤：S10：沿着指定编码方向每N行跳跃采集一行的一组欠采样K空间数据S(1)，其余数据用0填充；S11：在K空间的中心区域完整采集低频信息，以获得中心数据S(C)；和S12：生成N-1组K空间数据S(2)，S(3)，…S(N)，分别以相同的跳跃间隔N，不同的编码偏移量1,2…N-1，拷贝采集K空间S(C)对应位置的低频数据，其它未采集部分全部用0填充；以及迭代计算步骤S2包括以下步骤：S20：将K空间数据S(1)沿指定编码方向高通滤波，并经过傅氏变换(FourierTransform)后获得边缘增强图像I(1)；S21：将N-1组K空间数据S(2)-S(N)中的第n组沿指定编码方向高通滤波，并经过傅氏变换后获得混叠图像I(n)；S22：将混叠图像I(n)的幅度替换为当前幅度和边缘增强图像I(1)的幅度值的均值；S23：将获得的混叠图像I(n)的幅度反傅氏变换为新的K空间数据S(n)'；S24：利用中心数据S(C)更新K空间S(n)'的中心数据，并在K空间S(n)对应的跳跃间隔处填0；S25：重复步骤S21-S24，直到所有的N-1组数据都达到收敛状态为止；和S26：合并所有的N组K空间数据，并对合并后的K空间数据进行高通滤波的逆运算，然后经过傅氏变换获得最终的磁共振图像。进一步地，指定编码方向为相位编码方向或者3D序列层编码方向。可选地，指定编码方向包括相位编码和3D序列层编码，从而实现相位编码方向和3D序列层编码方向上的共同加速。欠采样步骤S20的K空间S(1)低频信息直接跳跃拷贝S(C)对应的数据，无需重复采集。中心区域为在K空间中心的L行范围内的区域，L为整数。当采集256点的原始数据时，中心区域为在K空间中心的8-64行范围内的区域。加速因子通过以下公式获得：f＝1/N+L/M-L/M/N，其中，f为加速因子，N为跳跃采集的跳跃宽度，L为K空间中心区域完整采样行数，M为原始编码方向数据长度。根据本专利技术的该磁共振成像方法不仅能够用多通道来实现，也可以用单通道来实现，仅采集部分欠采样数据，通过在图像域和K空间域之间进行迭代计算实现快速的磁共振成像，从而缩短扫描时间以使图像的采集时间大大缩短，或者在同样扫描时间的情况下改善信噪比，并可获得较高的成像质量。附图说明本专利技术的上述及其它方面和特征将从以下结合附图对实施例的说明清楚呈现，其中：图1A和图1B显示沿编码方向以等间隔N＝4跳跃采集的两组欠采样数据重建的图像，其中这两组欠采样数据具有不同的编码偏移量；图1C和图1D分别显示对图1A和图1B所示的图像进行高通滤波后获得的边缘增强图像；图2E-2H是通过本专利技术的HEIGHT方法对小腿进行扫描和重建的结果图，其中图2E是显示伪影的通过傅氏变换重建得到的图像，图2F是完整K空间中间的图像，图2G是相位编码跳跃间隔N＝4时重建的图像，以及图2H是跳跃间隔加大到N＝8时重建的图像；以及图3是0.3T低场永磁系统通过本专利技术的HEIGHT方法采集单通道的图像数据的结果图，其中图3I和图3L是完整K空间重建的图像，图3J和图3M是通过根据本专利技术的方法重建的图像。具体实施方式下面参照附图详细描述本专利技术的说明性、非限制性实施例，对根据本专利技术的磁共振成像方法进行进一步说明。根据本专利技术的磁共振成像方法是基于图像域和K空间域迭代的快速磁共振成像技术，该方法沿编码方向对K空间有规律地采集欠采样数据，并将采集到的欠采样数据通过图像域和K空间域之间的迭代计算得到未采集的K空间数据以获得最终的图像。参照图1A-1D示例性地说明跳跃采集获得的图像。图1A和图1B分别由两组欠采样数据重建而来，这两组数据在编码方向均以等间隔N＝4跳跃采集，但具有不同的编码偏移量d。从图1A和图1B中可以看出，由于编码偏移量不同，图1A和图1B的对比度明显不同，这是因为图中重叠的鬼影具有不同的相位。此时，如果对图1A和图1B所示的图像施以高通滤波，则分别生成图1C和图1D中所示的两个边缘图像。由于边缘图像是稀疏的，因此有效地减少了图像内信号的混叠，从而可见这两个边缘图像具有很高的相似性。这种沿编码方向跳跃间隔为N而偏移量d不同的图像共有N个。所述N个图像的稀疏边缘图像具有几乎相等的幅度。这种图像固有的自然属性使得本专利技术的磁共振成像方法只需要采集一组欠采样数据，便可基于图像域和K空间域之间的迭代计算外推出其它N-1组数据，从而实现快速的磁共振扫描成像。接下来，将详细说明根据本专利技术的磁共振成像方法。本专利技术的磁共振成像方法首先沿指定编码方向每N行采集一组欠采样K空间数据，然后基于图像域和K空间域之间的迭代计算获得未采集的N-1组K空间数据，从而获得最终的图像。具体地，根据本专利技术的磁共振成像方法包括以下步骤：S1:沿着指定编码方向对K空间进行有规律的欠采样；和S2：基于图像域和K空间域通过迭代计算获得未采集的K空间数据本文档来自技高网...

【技术保护点】
一种磁共振成像方法，包括以下步骤：S1:沿着指定编码方向对K空间进行有规律的欠采样；和S2：基于图像域和K空间域通过迭代计算获得未采集的K空间数据，以获得磁共振图像。

【技术特征摘要】
1.一种磁共振成像方法，包括以下步骤：S1:沿着指定编码方向对K空间进行有规律的欠采样；和S2：基于图像域和K空间域通过迭代计算获得未采集的K空间数据，以获得磁共振图像，其中，所述欠采样步骤S1包括以下步骤：S10：沿着所述指定编码方向每N行跳跃采集一行的一组欠采样K空间数据S(1)，其余数据用0填充；S11：在K空间的中心区域完整采集低频信息，以获得中心数据S(C)；和S12：生成N-1组K空间数据S(2)，S(3)，…S(N)，分别以相同的跳跃间隔N，不同的编码偏移量1,2…N-1，拷贝采集K空间S(C)对应位置的低频数据，其它未采集部分全部用0填充；以及所述迭代计算步骤S2包括以下步骤：S20：将所述K空间数据S(1)沿所述指定编码方向高通滤波，并经过傅氏变换后获得边缘增强图像I(1)；S21：将所述N-1组K空间数据S(2)-S(N)中的第n组沿所述指定编码方向高通滤波，并经过傅氏变换后获得混叠图像I(n)；S22：将所述混叠图像I(n)的幅度替换为当前幅度和所述边缘增强图像I(1)的幅度值的均值；S23：将获得的所述混叠图像I(n)的幅度反傅氏变换为新的K空间数据S(n)'；S24：利用所述中心数据S(C)更新K空间S(n)'的中心数据，并...

【专利技术属性】
技术研发人员：向清三，
申请(专利权)人：包头市稀宝博为医疗系统有限公司，
类型：发明
国别省市：内蒙古;15