三维多次激发扩散加权成像的K空间采集方法及重建方法技术

本发明专利技术涉及三维多次激发扩散加权成像的K空间采集方法及重建方法，基于三维轮式放射状EPI采集方式，每次激发后通过EPI方式采集三维轮式放射状K空间中的一个面，然后旋转读出梯度，再下一个角度继续激发和采集，完成一个三维轮式放射状K空间的采集；重建时，首先进行相位校正，然后对每次激发采集的数据进行二维傅里叶重建并取其模值，丢弃相位，获得一个无相位的三维放射状投影空间，最后对该投影空间逐层进行重建，得到三维的扩散加权成像。本发明专利技术可大大缩短扫描时间，重建运算量小，在重建过程中完全不受多次激发引入的随机相位的影响，稳定性好且不限激发次数，可以获得三维各向同性高分辨率扩散加权图像以及高b值三维扩散加权成像等。

【技术实现步骤摘要】
三维多次激发扩散加权成像的K空间采集方法及重建方法
本专利技术涉及核磁共振成像
，尤其涉及三维多次激发扩散加权成像的K空间采集方法及重建方法。
技术介绍
扩散加权成像(DiffusionWeightedImaging,DWI)是一种在分子水平上无创地反映活体水分子扩散运动的成像方法，是目前活体测量水分子扩散运动的唯一影像手段。扩散加权成像主要依赖于水分子的运动而非组织的质子密度、T1或T2弛豫时间。扩散加权成像适用于活体细胞水平探测生物组织的微动态和微结构变化，在肿瘤的良恶性鉴别、疗效评估和预测起着举足轻重的作用。目前，临床上使用的扩散成像方法通常是二维单次激发平面回波成像(EchoPlanarImaging，简称EPI)。二维单次激发EPI成像的特点是扫描时间短，受被试者的运动影响较小。然而，二维单次激发成像技术也有它本身的不足。第一，由于沿着相位编码方向的采集带宽较小，在磁介质率相差较大的不同组织交界处会产生较严重的图像变形；第二，单次激发方式获取高分辨率图像时，需要长回波链，长回波链意味着大的T2衰减，会导致图像模糊以及极大降低性噪比；第三，二维激发方式在层面方向分辨率有限，一般在4mm以上，否则在获取高分辨率图像时性噪比极低。为了减小图像变形，提高图像分辨率和信噪比，近年来提出了多次激发三维扩散成像。多次激发通过减小每次激发采集的相位编码的数目，提高相位编码方向采集带宽，因而能有效减少图像变形，达到较大的采集矩阵，提高空间分辨率。不过，因为施加了扩散梯度，病人自主运动，生理运动以及机械振动，会导致质子无法完全聚相，致使每次激发过程中图像会产生一个随机的相位误差。目前可通过在图像域进行相位校正，以去除多次激发扩散成像中的相位误差。图像域的相位校正方法主要有两类，第一类是每次激发成像完毕前额外采集导航数据，获得每次激发的相位信息，从而在重建过程中把相位移除。该类方法存在的主要缺陷是，导航数据与成像数据采集带宽无法做到完全一样，从而引起导航数据和成像数据的畸变不一致，就会影响最终的图像质量。第二类称为多路复用灵敏度编码方法(Multiplexedsensitivityencoding，MUSE)，即将其中每一次激发采集的数据，利用并行成像方法恢复出全采集的图像，然后在图像域去除相位。该类方法会随着激发次数的增加，稳定性和性噪比急遽下降，因此难以获得三维高分辨率扩散成像，同时该类方法运算量极大，重建时间一般长达数十分钟至数小时。多次激发三维扩散成像的技术难点是如何消除扩散梯度引入的随机相位对重建的影响。采用放射状K空间采样和滤波反投影重建扩散加权成像，可以消除扩散梯度引入的随机相位，减小运动伪影。现有的该类技术中，一次激发只采集1到3个回波，得到放射状K空间中的1到3条放射状K空间线。若要增加回波数目，每个回波对应的读出梯度都需要梯度回绕并旋转一个角度，这样的读出梯度形式极其复杂，难以校准，从而难以保证图像质量，并且梯度效率低。现有的该类技术要填满整个K空间，二维成像需要上百个重复周期(TR)，三维成像需要上千个重复周期(TR)，扫描时间会长达数小时，因此没有得到实际应用。
技术实现思路
本专利技术旨在提供三维多次激发扩散加权成像的K空间采集方法及重建方法，在重建过程中不受多次激发引入的随机相位的影响，可稳定获得三维扩散加权成像。为达到上述目的，本专利技术采用的技术方案如下：三维多次激发扩散加权成像的K空间采集方法，基于三维轮式放射状EPI采集方式，每次激发后通过EPI方式采集三维轮式放射状K空间中的一个面，然后旋转读出梯度，在下一个角度继续激发和采集，直至完成一个三维轮式放射状K空间的采集。具体包括以下步骤，步骤1，射频通道依次发射激发脉冲、第一回聚脉冲和第二回聚脉冲，相应的在逻辑坐标系中依次施加第一层选梯度，第二层选梯度和第三层选梯度；在激发脉冲后施加正向第一扩散梯度；在第一回聚脉冲后，施加反向第二扩散梯度，紧接着施加正向的第三扩散梯度；在第二回聚脉冲后施加反向第四扩散梯度；步骤2，在施加完扩散梯度后，在Z轴方向施加层选编码梯度，在X，Y轴施加EPI读出梯度；在多次激发时，层选编码梯度保持不变，仅旋转读出梯度X和读出梯度Y，采集得到多个旋转角度的轮式放射状EPI数据。其中，扩散梯度根据扩散梯度编码方向，施加在三个轴中的任意一个轴或多个轴。进一步的，在逻辑坐标系中的Z轴依次施加第一层选梯度，第二层选梯度和第三层选梯度；或者在逻辑坐标系中的Z轴施加第一层选梯度，在X轴施加第二层选梯度，在Y轴施加第三层选梯度。三维多次激发扩散加权成像的重建方法，包括以下步骤；第一步,利用上述的K空间采集方法采集三维放射状K空间数据；第二步,根据EPI读出轨迹进行网格化插值；第三步，利用参考扫描数据对单次激发数据进行EPI相位矫正；第四步，对单次激发对应的数据进行二维傅里叶重建；第五步，将第四步的重建结果取模值，得到不含相位信息的投影数据集；第六步，逐层进行基于投影数据的二维重建，得到了三维的扩散加权图像。其中，所述第三步中，将参考扫描获得的数据进行拟合，获得0阶和1阶相位，并从实际扫描数据中将0阶和1阶相位减去。进一步的，第六步中对投影数据逐层进行滤波反投影重建，或迭代重建，或基于深度学习的反投影重建，得到三维的扩散加权成像。进一步的，所述第二步包括以下步骤，对于每一次激发，采集到的K空间数据，按照EPI读出轨迹，填充于该放射状角度所对应的一个平面内。其中，所述第三步中参考扫描数据的采集方法包括以下步骤，步骤1，射频通道依次发射激发脉冲、第一回聚脉冲和第二回聚脉冲，相应的在逻辑坐标系中的Z轴依次施加第一层选梯度，第二层选梯度和第三层选梯度；其中，激发脉冲，第一回聚脉冲，第二回聚脉冲，第一层选梯度，第二层选梯度和第三层选梯度，所有的参数均与第一步中保持一致；步骤2，在X，Y轴施加EPI读出梯度；选择在N个放射状角度下采集参考扫描，旋转读出梯度X和读出梯度Y，采集得到若干个旋转角度的轮式放射状EPI数据；其中，N＝2～n,n为填满K空间所需的放射状射线总数目,若参考扫描的旋转角度数目少于实际扫描的旋转角度，则通过余弦插值的方式插值得到实际扫描旋转角度下的参考数据。与现有技术相比，本专利技术具有以下有益效果：第一，相对于传统的放射状K空间采集方式，本专利技术可大大缩短扫描时间；第二、本专利技术在重建过程中完全不受多次激发引入的随机相位的影响，可稳定地获得三维扩散加权成像；其稳定性好，不限激发次数，从而可以获得三维各向同性高分辨率扩散加权图像以及高b值三维扩散加权成像等；第三，本专利技术只需要简单的傅里叶重建和基于投影数据的重建，相对于现有的相位消除技术，重建运算量小。附图说明图1是三维轮式放射状K空间采集方式的扫描梯度/射频序列示意图；图2是三维轮式放射状K空间参考扫描的扫描梯度/射频序列示意图；图3是小视野三维轮式放射状K空间采集方式的扫描梯度/射频序列示意图；图4是三维轮式放射状K空间数据填充方式示意图；图5是三维多次激发扩散加权成像的重建方法的流程图；图6是利用本专利技术获得的健康志愿者三维全脑扩散加权图像；其中，a为横断位视野中典型的层面，b为3平面视野；图7是利用本专利技术获得的健康志愿者前列腺小视野三维扩散加权图像及其ADC图。具体本文档来自技高网...

【技术保护点】
1.三维多次激发扩散加权成像的K空间采集方法，其特征在于:基于三维轮式放射状EPI采集方式，每次激发后通过EPI方式采集三维轮式放射状K空间中的一个面，然后旋转读出梯度，在下一个角度继续激发和采集，直至完成一个三维轮式放射状K空间的采集。

【技术特征摘要】
1.三维多次激发扩散加权成像的K空间采集方法，其特征在于:基于三维轮式放射状EPI采集方式，每次激发后通过EPI方式采集三维轮式放射状K空间中的一个面，然后旋转读出梯度，在下一个角度继续激发和采集，直至完成一个三维轮式放射状K空间的采集。2.根据权利要求1所述的方法，其特征在于：包括以下步骤，步骤1，射频通道依次发射激发脉冲、第一回聚脉冲和第二回聚脉冲，相应的在逻辑坐标系中依次施加第一层选梯度，第二层选梯度和第三层选梯度；在激发脉冲后施加正向第一扩散梯度；在第一回聚脉冲后，施加反向第二扩散梯度，紧接着施加正向的第三扩散梯度；在第二回聚脉冲后施加反向第四扩散梯度；步骤2，在施加完扩散梯度后，在Z轴方向施加层选编码梯度，在X，Y轴施加EPI读出梯度；在多次激发时，层选编码梯度保持不变，仅旋转读出梯度X和读出梯度Y，采集得到多个旋转角度的轮式放射状EPI数据。3.根据权利要求2所述的方法，其特征在于：扩散梯度根据扩散梯度编码方向，施加在三个轴中的任意一个轴或多个轴。4.根据权利要求2所述的方法，其特征在于：在逻辑坐标系中的Z轴依次施加第一层选梯度，第二层选梯度和第三层选梯度。5.根据权利要求2所述的方法，其特征在于：在逻辑坐标系中的Z轴施加第一层选梯度，在X轴施加第二层选梯度，在Y轴施加第三层选梯度。6.三维多次激发扩散加权成像的重建方法，其特征在于：包括以下步骤；第一步,利用权利要求1-5中任一项所述的方法采集三维放射状K空间数据；第二步,根据EPI读出轨迹进行网格化插值；第三步，利用参考扫描数据对单次激发数据...

【专利技术属性】
技术研发人员：罗海，周翔，朱高杰，王超，陈梅泞，吕蓓，蒋辉，吴子岳，
申请(专利权)人：奥泰医疗系统有限责任公司，
类型：发明
国别省市：四川,51