[发明专利]一种DWI-EPI图像奈奎斯特重影校正方法

权利要求书

1.一种DWI-EPI图像奈奎斯特重影校正方法，其特征在于，包括以下步骤：

获取磁共振原始数据，对所述磁共振原始数据进行分割，提取用于参考扫描校正的参考扫描数据，将剩余数据填充相应的k空间位置，形成处于(Kx,Ky)域的k空间数据；

将处于(Kx,Ky)域的k空间数据和参考扫描数据均变换为(x,Ky)混合域数据；

对所述(x,Ky)混合域数据执行参考扫描校正，得到第一校正数据；

将处于(x,Ky)混合域的第一校正数据变换为(Kx,Ky)域数据；

对所述(Kx,Ky)域数据执行无参考扫描校正法校正，得到第二校正数据；

对所述第二校正数据对零填充处理，获得第三校正数据，对所述第三校正数据进行二维逆傅里叶变换得到无伪影图像。

2.根据权利要求1所述的DWI-EPI图像奈奎斯特重影校正方法，其特征在于，对所述磁共振原始数据进行分割具体为：

根据序列梯度信息、采集频率及系统延迟信息，舍弃噪声点和读梯度爬坡处采集的信号数据，只取循环采集中处于读梯度平稳段数据点，并剔除无相位编码数据，使数据矩阵尺寸与扫描时设置的矩阵尺寸一致。

3.根据权利要求1所述的DWI-EPI图像奈奎斯特重影校正方法，其特征在于，所述参考扫描数据为扫描采集的数据矩阵中Ky＝0的数据。

4.根据权利要求1所述的DWI-EPI图像奈奎斯特重影校正方法，其特征在于，填充相应的k空间位置时，将k空间数据中的负度梯度数据左右颠倒使得k空间数据度越方向都一致。

5.根据权利要求1所述的DWI-EPI图像奈奎斯特重影校正方法，其特征在于，分别将所述k空间数据和参考扫描数据沿Kx方向作反傅里叶变换将每个不同相位编码的信号变换到(x,Ky)混合域。

6.根据权利要求1所述的DWI-EPI图像奈奎斯特重影校正方法，其特征在于，所述参考扫描校正具体为：

(x,Ky)混合域的所述参考扫描数据包含三行数据，将第1行和第3行对应位置取平均得到插值行，将所述插值行乘以(x,Ky)混合域的所述参考扫描数据第2行的共轭，得到复数行数据；

使用一阶直线拟合方法基于所述复数行数据拟合获得相位差信息行数据；

将处于(x,Ky)混合域的图像数据中的正度越行剔除所述相位差信息行数据，获得所述第一校正数据。

7.根据权利要求6所述的DWI-EPI图像奈奎斯特重影校正方法，其特征在于，所述一阶直线拟合方法为最小二乘法。

8.根据权利要求1所述的DWI-EPI图像奈奎斯特重影校正方法，其特征在于，对处于(x,Ky)混合域的所述第一校正数据的x方向使用傅里叶变换，获得(Kx,Ky)域数据。

9.根据权利要求1所述的DWI-EPI图像奈奎斯特重影校正方法，其特征在于，所述无参考扫描校正法为改进信息熵校正法，校正过程包括：

生成信息熵校正法的常参数b和一阶线性参数k的参数范围，在所述参数范围内随机生成多个参数组合，每个参数组合对一条一阶线性相位误差；

对所述(Kx,Ky)域数据变换至(x,Ky)混合域，分别利用每条所述一阶线性相位误差对混合域数据进行相位误差剔除操作，得到多份校正后数据；

将多份所述校正后数据变换到(x,y)图像域，计算对应图像的熵值，获得熵值最小的一个参数组合，作为初始参数；

基于所述初始参数搜寻获得最优参数对，以所述最优参数对进行信息熵校正，获得所述第二校正数据。

10.根据权利要求9所述的DWI-EPI图像奈奎斯特重影校正方法，其特征在于，采用无导数法搜寻获得所述最优参数对。