[发明专利]k空间运动伪影矫正方法及装置

权利要求书

1.一种k空间运动伪影矫正方法，其特征在于，包括以下步骤：

S101，在成像区域连续采集多帧k空间数据，所述各帧内的k空间数据由至少一次射频激发采集得到，相邻帧数据的采集方向不同；

S102，将所述各帧中由一次射频激发采集得到的数据分为一组，所述各帧内至少有一个数据组；

S103，以前一帧中的一个数据组作为标准组，后一帧中的数据组作为待矫正组，利用所述标准组的数据对待矫正组的数据进行数据矫正，依次完成各帧数据组的数据矫正；

S104，将数据矫正后的各帧数据组填回相应帧的k空间中，得到各帧矫正完毕的k空间数据。

2.如权利要求1所述的k空间运动伪影矫正方法，其特征在于，所述步骤S101中采集多帧k空间数据时，奇数帧数据采集方向为kx方向，偶数帧数据采集方向为ky方向。

3.如权利要求1所述的k空间运动伪影矫正方法，其特征在于，所述步骤S103中，所述数据矫正包括：k空间中心矫正和/或旋转角度矫正和/或平移矫正。

4.如权利要求3所述的k空间运动伪影矫正方法，其特征在于，所述k空间中心矫正包括以下步骤：

将待矫正组k空间数据DataZ1变换到图像域得到数据ImageZ1；

将待矫正组k空间数据DataZ1加窗后变换到图像域得到数据ImageZ1Window；

取所述数据ImageZ1的模值为矫正后图像数据Image的模值；

取所述数据ImageZ1和所述数据ImageZ1Window的相位差为所述矫正后图像数据Image的相位；

将所述图像数据Image变换到k空间得到矫正完成数据DataZ1Corr。

5.如权利要求3所述的k空间运动伪影矫正方法，其特征在于，所述旋转角度矫正包括以下步骤：

将标准组数据和待矫正组数据分别取模值，所述标准组数据和待矫正组数据在ky和kx方向大小分别为K1和K2，所述标准组数据在ky方向步长为ΔK1，所述待矫正组数据在kx方向步长为ΔK2；

将取模值后的所述待矫正组数据在k空间内旋转α角度并插值到公共网格空间得到数据组MCi，所述公共网格空间在ky方向大小为K22，步长为ΔK22，在kx方向大小为K11，步长为ΔK11，所述MCi表示经过第i次旋转后插值到公共网格空间得到的数据组，i为正整数，α为0°到360°之间任一值；

将取模值后的所述标准组数据插值到所述公共网格空间得到数据组MB，计算MB与MCi的相关性；

找到最大相关性值对应的旋转角度α值，以所述旋转角度值为待矫正组数据作旋转角度矫正。

6.如权利要求3所述的k空间运动伪影矫正方法，其特征在于，所述平移矫正包括以下步骤：

将标准组数据和待矫正组数据插值到同一公共网格空间，所述标准组数据和待矫正组数据在ky和kx方向大小分别为K1和K2，所述标准组数据在ky方向步长为ΔK1，所述待矫正组数据在kx方向步长为ΔK2，所述公共网格空间在ky方向大小为K22，步长为ΔK22，在kx方向大小为K11，步长为ΔK11；

将公共网格空间内所述标准组数据与所述待矫正组数据的共轭的乘积变换到图像域得到图像，拟合找到图像最大值位置；

所述图像最大值位置与所述图像的中心位置的偏移量作为平移运动造成的偏移量，以所述偏移量为待矫正组作平移矫正。

7.如权利要求5或6所述的k空间运动伪影矫正方法，其特征在于，所述K11＝K1，K22＝K2，ΔK11＝ΔK1，ΔK22＝ΔK2。