[发明专利]一种磁共振成像多回波磁场分布图拟合方法

权利要求书

1.一种磁共振成像多回波磁场分布图拟合方法，其特征在于，该方法包括以下步骤：

步骤1：利用多回波梯度回波序列采集得到磁共振原始数据；

步骤2：对所有回波的相位进行线性度检查，舍弃线性度差的数据点的部分或者全部回波信号；

步骤3：将经过线性度检查后的相位图利用非线性最小二乘法进行拟合得到初步的磁场分布图；

步骤4：利用周围可靠数据点即相位线性度好的点的局部磁场的平均值来替代回波信号全部舍去的点的局部磁场，最终得到完整的磁场分布图；其中：

步骤2中所述的相位进行线性度检查，具体步骤为：

第一步，对前三个回波进行线性拟合得到初始相位近似值和角频率近似值ω′，具体表述为如下线性方程：

其中，为初始相位理论值，ω为角频率理论值，TE_j为第j个回波的回波时间，为第j个回波的相位；

第二步，奇异点判断

采用拟合得到的初始相位近似值和角频率近似值计算图像上每一个像素点前三个回波的拟合相位与测量相位之间的平均差值如下式：

设置一个误差阈值th，如果则认为该像素点相位误差严重，在场图上为奇异点，对奇异点之后再做进一步处理；

第三步，非奇异点回波可靠性判断

利用下式计算第j(N≥j＞3)个回波的拟合相位与测量相位的差值δ_j：

如果δ_j＞th，则认为该回波相位可靠性低，舍弃该回波及其后面的所有回波信号。

2.根据权利要求1所述的磁共振成像多回波磁场分布图拟合方法，其特征在于，步骤3中所述利用非线性最小二乘法进行拟合得到初步的磁场分布图是利用最小二乘法进行多回波场图拟合，根据拟合值与测量值之间的误差不断进行迭代，求解出局部磁场即磁场分布图的最优解，如下式：

其中，s(r，TE_j)为信号测量值，TE_j为回波时间，N为回波个数，经过相位检查后线性度不同的数据点所使用的回波个数不同，a(r，TE_j)为每个回波信号的幅值，γ为旋磁比，γΔB(r)为所求场图，单位为频率(Hz)。