[发明专利]基于磁共振图像的水脂分离方法、装置、设备及存储介质

权利要求书

1.一种基于磁共振图像的水脂分离方法，其特征在于，包括：

获取磁共振图像中的每个待处理像素点，分别计算每个所述待处理像素点的相位因子候选解，并将所述相位因子候选解中的各相位因子解分别划分至第一解集和第二解集中，其中，所述相位因子候选解包括全局最优解以及与所述全局最优解对应的分反解；

基于所述第一解集和所述第二解集中每个待处理像素点对应的相位因子解，确定出水脂转换区域，并计算出所述水脂转换区域内每个待处理像素点的目标相位因子解；

基于水脂转换区域以及所述水脂转换区域内每个待处理像素点的目标相位因子解确定出其余各所述待处理像素点的目标相位因子解；

根据所述磁共振图像中的每个待处理像素点的目标相位因子解提取所述磁共振图像中的水图和脂肪图。

2.根据权利要求1所述的方法，其特征在于，所述将所述相位因子候选解中的各相位因子解分别划分至第一解集和第二解集中，包括：

根据计算出的水脂分离结果将所述全局最优解以及所述分反解中的各相位因子解分别划分至第一解集和第二解集中。

3.根据权利要求2所述的方法，其特征在于，所述基于所述第一解集和所述第二解集中每个待处理像素点对应的相位因子解，确定出水脂转换区域，包括：

基于所述第一解集和所述第二解集中每个待处理像素点对应的相位因子解，分别计算每个待处理像素点与各邻域像素点在第一解集以及第二解集中的最大向量变化，根据第一解集以及第二解集对应的最大向量变化确定出水脂转换区域。

4.根据权利要求3所述的方法，其特征在于，所述基于所述第一解集和所述第二解集中每个待处理像素点对应的相位因子解，分别计算每个待处理像素点与各邻域像素点在第一解集中的最大向量变化，包括：

对于每一个待处理像素点，基于如下公式分别计算每个待处理像素点与各邻域像素点在第一解集中的最大向量变化：

其中，D_w(r)表示待处理像素点r与各邻域像素点在第一解集图中最大的向量变化；i表示各邻域像素点；abs(.)代表求取绝对值；angle(.)表示求取相位角；conj(.)表示求取复共轭；P_w(r_i)表示待处理像素点r的各邻域像素点i在第一解集中的相位因子解。

5.根据权利要求4所述的方法，其特征在于，所述根据第一解集以及第二解集对应的最大向量变化确定出水脂转换区域，包括：

如果当前待处理像素点在第一解集中的所述最大向量变化和在第二解集中的所述最大向量变化中的至少一个大于预设转换阈值时，则所述当前待处理像素点属于水脂转换区域，其中，所述预设转换阈值根据采样间隔时间内水脂信号之间的相位偏移确定。

6.根据权利要求3所述的方法，其特征在于，所述计算出所述水脂转换区域内每个待处理像素点的目标相位因子解，包括：

确定所述水脂转换区域中每个待处理像素点的各邻域像素点中所述最大向量变化最大的目标邻域像素点；

当采样间隔时间内水脂信号之间的相位偏移不等于180度的整数倍时，根据所述水脂转换区域中当前待处理像素点以及与当前待处理像素点对应的目标邻域像素点的相位因子解组合确定所述当前待处理像素点的目标相位因子解。

7.根据权利要求1所述的方法，其特征在于，所述基于水脂转换区域以及所述水脂转换区域内每个待处理像素点的目标相位因子解确定出其余各待处理像素点的目标相位因子解，包括：

根据所述第一解集和所述第二解集将除水脂转换区域之外的其余待处理像素点划分为第一数量的在空间上连续的待求解子区域，其中，每个待求解子区域内的所有待处理像素点都来自于同一个解集；

根据所述水脂转换区域内每个待处理像素点的目标相位因子解确定各所述待求解子区域中各待处理像素点的目标相位因子解。