[发明专利]一种磁共振化学位移编码成像方法、装置及设备

权利要求书

1.一种磁共振化学位移编码成像方法，其特征在于，包括：

在预定分辨率下的采样图像中，基于两点磁共振信号模型建立相位因子-拟合误差phasor-error谱；

将所述phasor-error谱中具有唯一的相位因子值且使所述phasor-error谱达到局部极小值的像素点确定为初始种子点；

根据所述初始种子点对待估像素点的相位因子值进行估计，获得在所述预定分辨率下的场图；

将所述预定分辨率下的场图在最高分辨率下分别进行映射，获得在所述最高分辨率下的多幅场图，并将所述多幅场图合并后获得重构场图；

从所述重构场图中确定重构种子点，并基于所述重构种子点进行估计以获得重构待估像素点的相位因子值；

根据所述重构种子点和所述重构待估像素点的相位因子值得到两种预定成分的分离图像。

2.如权利要求1所述的方法，其特征在于，将具有唯一相位因子值且使phasor-error谱达到局部极小值的像素点确定为初始种子点包括：

获得所述预定分辨率图像中的每一个像素点的phasor-error谱，并根据所述phasor-error谱确定使所述像素点达到局部极小值的相位因子值；

若所述像素点的信噪比大于预定值且具有唯一的相位因子值，则将所述像素点确定为初始种子点。

3.如权利要求2所述的方法，其特征在于，确定所述像素点的具有唯一的相位因子值包括：

若所述像素点只有一个局部极小值使得拟合误差小于预定值，则将所述局部极小值确定为所述像素点的唯一的相位因子值。

4.如权利要求3所述的方法，其特征在于，根据所述初始种子点对待估像素点的相位因子值进行估计包括：

若所述待估像素点包括有多个局部极小值，则将拟合误差小于预定值的局部极小值确定为所述待估像素点的多个相位因子值；

从所述待估像素点的多个相位因子值中选取与相邻种子点的相位因子值具有最大相似性的相位因子值作为所述待估像素点的相位因子值。

5.如权利要求4所述的方法，其特征在于，待估像素点的多个相位因子值与相邻种子点的相位因子值的最大相似性通过以下公式确定：

d＝argmax{D₁,D₂,...,D_S}

其中，D₁,D₂,…,D_s表示phasor相似性，并通过以下公式确定：

其中，K表示当前待估像素点的八邻域内的种子点个数；m_k表示第k个种子点的幅值，相应的相位因子值为d_k，d_s表示当前待估像素点的一个相位因子候选值；conj(.)表示取复共轭；angle(.)表示取角度。

6.如权利要求5所述的方法，其特征在于，根据所述初始种子点对待估像素点的相位因子值进行估计还包括：

选取幅值大于第一阈值的待估像素点建立第一待估像素点集合；

获取所述第一待估像素点集合中每个待估像素点的相邻种子点中两两之间的相位差，选取最大相位差小于第二阈值的待估像素点建立第二待估像素点集合；

将所述第二待估像素点集合中的待估像素点按所述相邻种子点的个数进行降序排列，选取预定数量的待估像素点建立第三待估像素点集合；

获取所述第三待估像素点集合中的待估像素点与所述相邻种子点之间最大的相位相似性，并将所述相似性大于第三阈值的待估像素点的相位因子值作为待估像素点的相位因子值；

将所述相似性大于第三阈值的待估像素点添加为种子点对所述第一待估像素点集合中的待估像素点进行估计，直至所述第一待估像素点集合中的待估像素点的相邻种子点数为零。