[发明专利]径向重排方法

说明书

技术领域

本发明涉及磁共振成像(MRI)技术领域，特别是一种径向重排(radial reordering)方法，用于对MRI中回波的数据点进行重排。

背景技术

在重聚焦脉冲中使用非选择性的脉冲和调制的翻转角，可以提高3D快速自旋回波成像中的采样效率。有人提出了线性模式(linear mode)和中心径向模式(centric radial mode)这两种灵活重排方案，以避开分段3D成像中矩形k空间中采样标记矩阵的限制，提供灵活的k空间覆盖度。其中，线性模式主要适用于中等回波时间(Time to Echo，TE)成像，中心径向模式主要适用于短TE成像。

总体而言，中心径向模式可以提供较好的对比度纯度(contrast purity)，但是该模式只适用于短TE成像。线性模式适用于中等TE(约为40～50ms)成像，但是只能提供质量很差的对比度纯度。

在所提及的中心径向重排模式中，来自第一个回波或最后一个回波的数据填充到K空间中心，这种径向重排模式对于超短TE或超长TE成像来说效果不错。然而，中心径向重排模式不支持中等TE值，而中等TE值在一些躯干成像和整形外科成像的应用中非常有用。与线性重排模式相比，由于中等TE成像中的径向重排模式能够提供较好的对比度加权，降低k空间中部的亮度波动，以及降低部分傅立叶(Partial Fourier，PF)成像期间的图像模糊，所以有时需要优选采用中等TE成像中的径向重排模式。

因此，在MRI成像中急需一种能够支持中等TE成像的径向模式的重排技术。

发明内容

有鉴于此，本发明提供一种径向重排方法，其目的在于，本发明的径向重排方法能够支持中等TE成像。

因此，本发明提供了一种径向重排方法，该方法包括以下步骤：

计算k空间的Kz-Ky平面内所需采集的数据点的标记矩阵；

在所述平面中选择一参考轴，该参考轴经过k空间中心；

计算各数据点的半径和方位角，其中，如果任一数据点的方位角大于π，则将该数据点的半径改为负值，并将该数据点的方位角减去π；

按照半径升序模式，对所有数据点进行排序；

从所述标记矩阵中删除已排序序列中最先的数据点；

按照各数据点在已排序序列中的位置，为剩余各个数据点分配升序的回波号；

对于具有相同回波号的数据点，按照方位角升序模式，为这些数据点分配升序的回波链号。

所述参考轴沿Kz方向或者Ky方向。优选地，所述参考轴沿Kz方向和Ky方向中编码步较少的方向。

该方法进一步包括：利用公式f_PF＝(0.5L+M)/L计算部分傅立叶因子，其中f_PF为部分傅立叶因子，L为快速因子，M为k空间中心的回波号。

所删除的数据点的数量由部分傅立叶因子确定。

该方法进一步包括：将回波链的长度从L缩短为0.5L+M，其中L为快速因子，M为k空间中心的回波号。

本发明还提供了另一种径向重排方法，该方法包括以下步骤：

计算k空间的Ky-Kz平面内所需采集的数据点的标记矩阵；

在所述平面中选择一参考轴，该参考轴经过k空间中心；

删除所述标记矩阵中根据部分傅立叶采集所要省略的数据点；

计算剩余各数据点的半径和方位角，其中，如果任一数据点的方位角大于π，则将该数据点的半径改为负值，并将该数据点的方位角减去π；

按照半径升序模式，对所有数据点进行排序；

按照各数据点在已排序序列中的位置，为各个数据点分配升序的回波号；

对于具有相同回波号的数据点，按照方位角升序模式，为这些数据点分配升序的回波链号。

所述参考轴沿Kz方向或者Ky方向。优选地，所述参考轴沿Kz方向和Ky方向中编码步较少的方向。

该方法进一步包括：利用公式f_PF＝(0.5L+M)/L计算部分傅立叶因子，其中f_PF为部分傅立叶因子，L为快速因子，M为k空间中心的回波号。

该方法进一步包括：将回波链的长度从L缩短为0.5L+M，其中L为快速因子，M为k空间中心的回波号。

从上述方案中可以看出，由于本发明在k空间独特地对数据点重排，使得径向重排能够用于中等TE成像，并且提供了比中等TE成像中线性重排方案更好的图像质量。

附图说明