一种磁共振化学位移编码成像方法、装置及设备，属于磁共振成像技术领域。所述方法包括：在基于两点磁共振信号模型建立的phasor‑error谱中将具有唯一相位因子值且使该谱达到局部极小值的像素点确定为初始种子点；根据初始种子点对待估像素点的相位因子值进行估计获得场图；将场图在最高分辨率下进行映射及合并后获得重构场图；从重构场图中确定重构种子点，并通过对重构种子点进行估计以获得重构待估像素点的相位因子值；根据重构种子点和重构待估像素点的相位因子值得到两种预定成分的分离图像。所述方法将同时含有两种成分的区域识别为种子点，排除了最高分辨率下由于相位因子值跳变造成的偏差，保证最终种子点选择的正确性。
1.一种磁共振化学位移编码成像方法，其特征在于，包括：
在预定分辨率下的采样图像中，基于两点磁共振信号模型建立相位因子-拟合误差phasor-error谱；
将所述phasor-error谱中具有唯一的相位因子值且使所述phasor-error谱达到局部极小值的像素点确定为初始种子点；
根据所述初始种子点对待估像素点的相位因子值进行估计，获得在所述预定分辨率下的场图；
将所述预定分辨率下的场图在最高分辨率下分别进行映射，获得在所述最高分辨率下的多幅场图，并将所述多幅场图合并后获得重构场图；
从所述重构场图中确定重构种子点，并基于所述重构种子点进行估计以获得重构待估像素点的相位因子值；
根据所述重构种子点和所述重构待估像素点的相位因子值得到两种预定成分的分离图像。
2.如权利要求1所述的方法，其特征在于，将具有唯一相位因子值且使phasor-error谱达到局部极小值的像素点确定为初始种子点包括：获得所述预定分辨率图像中的每一个像素点的phasor-error谱，并根据所述phasor-error谱确定使所述像素点达到局部极小值的相位因子值；
若所述像素点的信噪比大于预定值且具有唯一的相位因子值，则将所述像素点确定为初始种子点。
3.如权利要求2所述的方法，其特征在于，确定所述像素点的具有唯一的相位因子值包括：若所述像素点只有一个局部极小值使得拟合误差小于预定值，则将所述局部极小值确定为所述像素点的唯一的相位因子值。
4.如权利要求3所述的方法，其特征在于，根据所述初始种子点对待估像素点的相位因子值进行估计包括：若所述待估像素点包括有多个局部极小值，则将拟合误差小于预定值的局部极小值确定为所述待估像素点的多个相位因子值；
从所述待估像素点的多个相位因子值中选取与相邻种子点的相位因子值具有最大相似性的相位因子值作为所述待估像素点的相位因子值。
5.如权利要求4所述的方法，其特征在于，待估像素点的多个相位因子值与相邻种子点的相位因子值的最大相似性通过以下公式确定：d＝argmax{D1,D2,...,DS}
其中，D1,D2,…,Ds表示phasor相似性，并通过以下公式确定：其中，K表示当前待估像素点的八邻域内的种子点个数；mk表示第k个种子点的幅值，相应的相位因子值为dk，ds表示当前待估像素点的一个相位因子候选值；conj(.)表示取复共轭；angle(.)表示取角度。
6.如权利要求5所述的方法，其特征在于，根据所述初始种子点对待估像素点的相位因子值进行估计还包括：选取幅值大于第一阈值的待估像素点建立第一待估像素点集合；
获取所述第一待估像素点集合中每个待估像素点的相邻种子点中两两之间的相位差，选取最大相位差小于第二阈值的待估像素点建立第二待估像素点集合；
将所述第二待估像素点集合中的待估像素点按所述相邻种子点的个数进行降序排列，选取预定数量的待估像素点建立第三待估像素点集合；
获取所述第三待估像素点集合中的待估像素点与所述相邻种子点之间最大的相位相似性，并将所述相似性大于第三阈值的待估像素点的相位因子值作为待估像素点的相位因子值；
将所述相似性大于第三阈值的待估像素点添加为种子点对所述第一待估像素点集合中的待估像素点进行估计，直至所述第一待估像素点集合中的待估像素点的相邻种子点数为零。
7.如权利要求6所述的方法，其特征在于，根据所述初始种子点对待估像素点的相位因子值进行估计还包括：通过将所述第一阈值减小预定量以确定第四待估像素点集合，并确定所述第四待估像素点集合中的待估像素点的相位因子值，直至所述第一阈值减为零；
将所述第一待估像素点集合中未确定相位因子值的待估像素点的相位因子值确定为相邻种子点的相位因子值的平均值。
8.如权利要求2所述的方法，其特征在于，所述phasor-error谱通过以下公式确定：其中，d表示phasor-error谱；I表示4×4矩阵；上标“+”表示取矩阵违逆，A+＝(ATA)-1AT；
连接符 表示点乘；err(d)表示phasor d与拟合误差之间的关系； 和 均表示预定行向量；S为磁共振信号模型矩阵，A为S的系数矩阵。
9.如权利要求1所述的方法，其特征在于，将多个预定分辨率下的场图在最高分辨率下分别进行映射包括：获取最高分辨率下每一个重构待估像素点的所有相位因子值与对应的预定分辨率下的像素点的相位因子值之间的相位相似性，将最高相似性的相位因子值确定为所述重构待估像素点的相位因子值。
10.如权利要求9所述的方法，其特征在于，确定所述重构待估像素点的相位因子值包括：若所述重构待估像素点在所有的场图中均具有相同的相位因子值，则将所述重构待估像素点确定为新的种子点以及将所述相同的相位因子值确定为所述重构待估像素点的相位因子值；
若所述重构待估像素点在所有的场图中具有不同的相位因子值，则将所述重构待估像素点重新设定为待估像素点以及将所述不同的相位因子值确定为相位因子候选值。
11.如权利要求10所述的方法，其特征在于，确定所述重构待估像素点的相位因子值还包括：以所述新的种子点为起点在最高分辨率下对所述相位因子候选值进行估计，以确定所述重构待估像素点的相位因子值。
12.如权利要求1所述的方法，其特征在于，根据所述重构种子点和所述重构待估像素点的相位因子值得到两种预定成分的分离图像包括：将所述重构场图进行平滑处理后分别获得两种成分的图像。
13.一种磁共振化学位移编码成像装置，其特征在于，所述装置包括：相位误差谱建立模块，用于在预定分辨率下的采样图像中，基于两点磁共振信号模型建立相位因子-拟合误差phasor-error谱；
种子点选取模块，用于将所述phasor-error谱中具有唯一的相位因子值且使所述phasor-error谱达到局部极小值的像素点确定为初始种子点；
估计模块，用于根据所述初始种子点对待估像素点的相位因子值进行估计，获得在所述预定分辨率下的场图；
场图重构模块，用于将所述预定分辨率下的场图在最高分辨率下分别进行映射，获得在所述最高分辨率下的多幅场图，并将所述多幅场图合并后获得重构场图；
相位因子确定模块，用于从所述重构场图中确定重构种子点，并基于所述重构种子点进行估计以获得重构待估像素点的相位因子值；
成像模块，用于根据所述重构种子点和所述重构待估像素点的相位因子值得到两种预定成分的分离图像。
14.如权利要求13所述的装置，其特征在于，在所述种子点选取模块中包括：极小值确定子模块，获得所述预定分辨率图像中的每一个像素点的phasor-error谱，并根据所述phasor-error谱确定使所述像素点达到局部极小值的相位因子值；
种子点确定子模块，用于若所述像素点的信噪比大于预定值且具有唯一的相位因子值，则将所述像素点确定为初始种子点。
15.如权利要求14所述的装置，其特征在于，所述种子点确定子模块用于若所述待估像素点包括有多个局部极小值，则将拟合误差小于预定值的局部极小值确定为所述待估像素点的多个相位因子值，并从所述待估像素点的多个相位因子值中选取与相邻种子点的相位因子值具有最大相似性的相位因子值作为所述待估像素点的相位因子值。
16.如权利要求15所述的装置，其特征在于，在所述估计模块中包括：第一集合建立子模块，用于选取幅值大于第一阈值的待估像素点建立第一待估像素点集合；
第二集合建立子模块，用于获取所述第一待估像素点集合中每个待估像素点的相邻种子点中两两之间的相位差，选取最大相位差小于第二阈值的待估像素点建立第二待估像素点集合；
第三集合建立子模块，用于将所述第二待估像素点集合中的待估像素点按所述相邻种子点的个数进行降序排列，选取预定数量的待估像素点建立第三待估像素点集合；
第一相位因子确定子模块，用于获取所述第三待估像素点集合中的待估像素点与所述相邻种子点之间最大的相位相似性，并将所述相似性大于第三阈值的待估像素点的相位因子值作为待估像素点的相位因子值；
第一估计子模块，用于将所述相似性大于第三阈值的待估像素点作为新的种子点对所述第一待估像素点集合中的待估像素点进行估计，直至所述第一待估像素点集合中的待估像素点的相邻种子点数为零。
17.如权利要求16所述的装置，其特征在于，在所述估计模块中还包括：第四集合建立模块，用于通过将所述第一阈值减小预定量以确定第四待估像素点集合；
第二相位因子确定子模块，用于确定所述第四待估像素点集合中的待估像素点的相位因子值，直至所述第一阈值减为零；
第三相位因子确定子模块，用于将所述第一待估像素点集合中未确定相位因子值的待估像素点的相位因子值确定为相邻种子点的相位因子值的平均值。
18.如权利要求13所述的装置，其特征在于，在所述场图重构模块中包括：第四相位因子确定子模块，用于获取最高分辨率下每一个重构待估像素点的所有相位因子值与对应的预定分辨率下的待估像素点的相位因子值之间的相位相似性，将最高相似性的相位因子值确定为所述重构待估像素点的相位因子值。
19.如权利要求18所述的装置，其特征在于，在所述相位因子确定模块中包括：第五相位因子确定子模块，用于若所述重构待估像素点在所有的场图中均具有相同的相位因子值，则将所述重构待估像素点确定为新的种子点以及将所述相同的相位因子值确定为所述重构待估像素点的相位因子值；
第六相位因子确定子模块，用于若所述重构待估像素点在所有的场图中具有不同的相位因子值，则将所述重构待估像素点重新设定为待估像素点以及将所述不同的相位因子值确定为相位因子候选值；
第七相位因子确定子模块，用于以所述新的种子点为起点在最高分辨率下对所述相位因子候选值进行估计，以确定所述重构待估像素点的相位因子值。
20.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质包括计算机指令，当所述计算机指令被执行时，如权利要求1至12任意一项所述的方法被执行。
包括存储有计算机可读指令的存储器，所述处理器被配置来通过执行所述计算机可读指令以执行如权利要求1至12任意一项所述的方法。
一种磁共振化学位移编码成像方法、装置及设备
技术领域
[0001] 本发明涉及一种磁共振化学位移编码成像方法、装置及设备，属于磁共振成像技术领域。
背景技术
[0002] 磁共振化学位移编码成像是一种基于组织中各成分之间化学位移差异的成像方法，通过在不同的回波时间采集信号，并通过拟合混合信号模型得到。临床上最常用的化学位移编码成像是水脂分离成像，它主要用于脂肪抑制和脂肪定量等应用中。在分离过程中，通过多回波化学位移编码成像可对B0局部场进行校正，并同时得到纯水和纯脂肪的图像。但是，在B0场偏较大或成像组织空间分离的情况下，传统的多回波化学位移编码成像在估计B0场偏时容易收敛到一个错误的局部极小值，进而造成水脂分反的结果。
[0003] 多回波化学位移编码成像的原理是假设组织信号由多种不同成分被同时激发得到，并已知各成分相对于水的化学位移，然后拟合信号与各成分、B0场图和回波时间(Echo Time，TE)之间的数学模型。实际应用最多的化学位移编码成像是水脂分离成像，即通过多回波信号分离水和脂肪两种成分。
[0004] 常用多回波水脂分离成像根据其TE个数的不同可分为两点分离技术和多点(大于2)分离技术。两点分离技术通常只能分离具有两种化学位移的成分，但是在一些要求快速成像的应用中，两点分离技术具有更快的采集及成像速度。
[0005] 在现有技术中有一种两点水脂分离技术，该技术假设场图分布均匀，选取两个特定的回波时间分别采集水和脂肪质子相位同向(In-phase，IP)和反向(Opposed-phase，OP)的图像，通过简单的相加和/或相减操作得到水图和脂肪图。为了解决B0场不均匀造成的影响，在现有技术中还有一种基于局部增长的两点水脂分离技术，该技术利用待估像素点和种子点之间的角度差确定局部增长的路径，并根据待估像素点候选场图值和其邻域范围内所有种子点的场图加和的差值来确定该待估像素点的场图，但该方法只能得到分离的纯水和纯脂肪图像，需要加入额外的信息才能识别出二者中的纯水图和纯脂肪图。
[0006] 另外，在现有技术中还有一种两点水脂分离技术，该技术不再限制第二幅采集图像中水和脂肪相位反向，而是采集水脂相位部分反向(partially-opposed-phase，POP)图像，利用余弦定理得到两组水和脂肪的幅值图像，在此基础上解得两个场图候选值然后将原始图像分割为多个不交叠的子块，根据子块内场图分布的光滑特性，选取出每一个子块内具有最大场图相似性的场图集合作为初始场图，利用区域迭代场图提取算法从场图候选解中将整个原始图像空间内的场图估计出来，最后利用最小二乘法得到分离的水和脂肪图像。但该方法的缺陷是在场图提取过程中容易出现选取矛盾，即从不同方向进行场图选取时得到不同的选取结果，造成最后的场图中存在跳变。
发明内容
[0007] 本发明为解决现有的化学位移编码成像方法存在的只能分离两种成分信号却不能单独识别其中一种成分以及场图中存在跳变的问题，进而提出了一种磁共振化学位移编码成像方法、装置及设备，具体包括如下的技术方案：
[0008] 一种磁共振化学位移编码成像方法，包括：
[0009] 在预定分辨率下的采样图像中，基于两点磁共振信号模型建立phasor-error谱；
[0010] 将所述phasor-error谱中具有唯一的相位因子值且使所述phasor-error谱达到局部极小值的像素点确定为初始种子点；
[0011] 根据所述初始种子点对待估像素点的相位因子值进行估计，获得在所述预定分辨率下的场图；
[0012] 将所述预定分辨率下的场图在最高分辨率下分别进行映射，获得在所述最高分辨率下的多幅场图，并将所述多幅场图合并后获得重构场图；
[0013] 从所述重构场图中确定重构种子点，并基于所述重构种子点进行估计以获得重构待估像素点的相位因子值；
[0014] 根据所述重构种子点和所述重构待估像素点的相位因子值得到两种预定成分的分离图像。
[0015] 一种磁共振化学位移编码成像装置，所述装置包括：
[0016] 相位误差谱建立模块，用于在预定分辨率下的采样图像中，基于两点磁共振信号模型建立phasor-error谱；
[0017] 种子点选取模块，用于将所述phasor-error谱中具有唯一的相位因子值且使所述phasor-error谱达到局部极小值的像素点确定为初始种子点；
[0018] 估计模块，用于根据所述初始种子点对待估像素点的相位因子值进行估计，获得在所述预定分辨率下的场图；
[0019] 场图重构模块，用于将所述预定分辨率下的场图在最高分辨率下分别进行映射，获得在所述最高分辨率下的多幅场图，并将所述多幅场图合并后获得重构场图；
[0020] 相位因子确定模块，用于从所述重构场图中确定重构种子点，并基于所述重构种子点进行估计以获得重构待估像素点的相位因子值；
[0021] 成像模块，用于根据所述重构种子点和所述重构待估像素点的相位因子值得到两种预定成分的分离图像。
[0022] 本发明的有益效果是：通过在低分辨率下进行场图估计，可在低分辨率下将同时含有两种成分的区域识别为种子点，增加了种子点的数量和分布范围，利用自检验机制将多幅高分辨率的场图进行合并，排除了最高分辨率下由于相位因子值跳变造成的偏差，保证了最终种子点选择的正确性，从而能够对两者成分的信号进行准确的识别。
附图说明
[0023] 图1以示例的方式示出了磁共振化学位移编码成像方法的流程图。
[0024] 图2以示例的方式示出了基于多分辨率最安全路径局部增长的自检验场图估计方法的流程图。
[0025] 图3以示例的方式示出了一个种子点的phasor-error谱图，其中的横坐标表示相位因子(phasor)，纵坐标表示对应的拟合误差(err(d)，无量纲)，图中的直线表示阈值、曲线表示拟合误差。
[0026] 图4以示例的方式示出了最安全路径局部增长算法的流程图。
[0027] 图5以示例的方式示出了磁共振化学位移编码成像装置的结构图。
[0028] 图6以示例的方式示出了计算机设备的结构图。
[0029] 图7为实施例一提出的验证所述磁共振化学位移编码成像方法的颈部水脂分离实例图，图中从左向右依次为原始图像、分离脂肪图、分离水图和估计得到的场图。
[0030] 图8为实施例一提出的验证所述磁共振化学位移编码成像方法的腹部水脂分离实例图，图中从左向右依次为原始图像、分离脂肪图、分离水图和估计得到的场图。
具体实施方式
[0031] 在现有技术中，研究者在现有的两点分离方法的基础上提出了一种两点分离技术，该技术不再限制第一幅采集图像中水和脂肪相位同向，只是限制第一幅采集图像中的水脂相位差θ1和第二幅图像中的水脂相位差θ2之间的差或和不能等于2π或其整数倍。另外，还有一种不对TE限制的两点分离技术，该技术在求得两组场图候选解的基础上将场图估计转化为全局能量最小化的问题，并采用树消息传播算法对其进行估计。
[0032] 但上述现有算法均存在其缺点：上述第一种两点分离技术在不同的子块得到各自的场图集合，将其直接拼接后作为初始场图进行区域迭代场图提取，这样做的缺陷是在场图提取过程中容易出现选取矛盾，即从不同方向进行场图选取时得到不同的选取结果，造成最后的场图中存在跳变；上述第二种两点分离技术在处理IP/OP图像时，因其两组场图候选解的相位相差π，使得正确的场图和反转的场图具有相同的能量，造成树消息传播算法容易陷入全局场图反转的情况。
[0033] 针对上述缺点，本发明提出一种两点水脂分离方法，该方法是一种基于多分辨率最安全路径局部增长策略和自校验场图估计的方法，能够实现正确的场图估计和水脂分离。
[0034] 结合图1所示，本发明提出的磁共振化学位移编码成像方法包括：
[0035] 步骤11，在预定分辨率下的采样图像中，基于两点磁共振信号模型建立phasor-error谱。
[0036] 其中，本步骤可将二维的处于最高分辨率的原始图像进行低通滤波并根据预先定义的欠采样系数进行欠采样，得到不同的低分辨率下的二维原始图像，低分辨率下的二维原始图像中的一个像素对应最高分辨图像中的多个像素，将这些像素作为一个整体描绘其phasor-error(相位因子-拟合误差)谱，即为低分辨下一个像素的phasor-error谱。
[0037] 步骤12，将phasor-error谱中具有唯一的相位因子值使该phasor-error谱达到局部极小值的像素点确定为初始种子点。
[0038] 其中，本步骤根据phasor-error谱确定使每个像素点达到局部极小值的相位因子值。若该像素点的信噪比大于预定值且具有唯一的相位因子值，则可将该像素点确定为种子点。
[0039] 可选的，将相位因子值达到局部极小值的像素点确定为初始种子点包括：
[0040] 获得预定分辨率图像中的每一个像素点的phasor-error谱，并根据该phasor-error谱确定使像素点达到局部极小值的相位因子值；
[0041] 若该像素点的信噪比大于预定值且具有唯一的相位因子值，则将该像素点确定为初始种子点。
[0042] 其中，确定该像素点的具有唯一的相位因子值的过程可以包括：
[0043] 若该像素点只有一个局部极小值使得拟合误差小于预定值，则将该局部极小值确定为该像素点的唯一的相位因子值。
[0044] 步骤13，根据初始种子点对待估像素点的相位因子值进行估计，获得在预定分辨率下的场图。
[0045] 其中，对待估像素点的相位因子值进行估计的过程包括：
[0046] 若待估像素点包括有多个局部极小值，则将拟合误差小于预定值的局部极小值确定为该待估像素点的多个相位因子值；
[0047] 从该待估像素点的多个相位因子值中选取与相邻种子点的相位因子值具有最大相似性的相位因子值作为该待估像素点的相位因子值。
[0048] 可选的，可通过最安全路径局部增长算法确定该待估像素点的相位因子值，其过程包括：
[0049] 选取幅值大于第一阈值的待估像素点建立第一待估像素点集合；
[0050] 获取第一待估像素点集合中每个待估像素点的相邻种子点中两两之间的相位差，选取最大相位差小于第二阈值的待估像素点建立第二待估像素点集合；
[0051] 将第二待估像素点集合中的待估像素点按相邻种子点的个数进行降序排列，选取预定数量的待估像素点建立第三待估像素点集合；
[0052] 获取第三待估像素点集合中的待估像素点与相邻种子点之间最大的相位相似性，并将相似性大于第三阈值的待估像素点的相位因子值作为待估像素点的相位因子值；
[0053] 将相似性大于第三阈值的待估像素点添加为种子点对第一待估像素点集合中的待估像素点进行估计，直至第一待估像素点集合中的待估像素点的相邻种子点数为零。
[0054] 可选的，本步骤还可对第一待估像素点集合中的待估像素点的相位因子值进行确定，其过程包括：
[0055] 通过将第一阈值减小预定量确定第四待估像素点集合，并确定第四待估像素点集合中的待估像素点的相位因子值，直至第一阈值减为零；
[0056] 将第一待估像素点集合中未确定相位因子值的待估像素点的相位因子值确定为相邻种子点相位因子值的平均值。
[0057] 步骤14，将预定分辨率下的场图在最高分辨率下分别进行映射，获得在最高分辨率下的多幅场图，并将多幅场图合并后获得重构场图。
[0058] 其中的映射可以是以步骤13中获得的在预定分辨率下的场图为初始值，对最高分辨率下的场图进行选取，其过程可以包括：
[0059] 获取最高分辨率下每一个重构待估像素点的所有相位因子值与对应的预定分辨率下的像素点的相位因子值之间的相位相似性，将具有最高相似性的相位因子值确定为该重构待估像素点的相位因子值。
[0060] 步骤15，从重构场图中确定重构种子点，并基于重构种子点进行估计以获得重构待估像素点的相位因子值。
[0061] 其中，确定该重构待估像素点的相位因子值的过程可以包括：
[0062] 若重构待估像素点在所有的场图中均具有相同的相位因子值，则将重构待估像素点确定为新的种子点以及将相同的相位因子值确定为重构待估像素点的相位因子值；
[0063] 若重构待估像素点在所有的场图中具有不同的相位因子值，则将重构待估像素点重新设定为待估像素点以及将不同的相位因子值确定为相位因子候选值；
[0064] 以新的种子点为起点，并可采用如步骤13所述的估计方法在最高分辨率下对相位因子候选值进行估计，以确定重构待估像素点的相位因子值。
[0065] 步骤16，根据重构种子点和重构待估像素点的相位因子值得到两种预定成分的分离图像。
[0066] 在确定重构种子点和重构待估像素点的相位因子值后，利用重构场图进行平滑处理后即可用于得到两种预定成分的分离图像。
[0067] 采用本实施例提出的磁共振化学位移编码成像方法，通过在低分辨率下进行场图估计，可在低分辨率下将同时含有两种成分的区域识别为种子点，增加了种子点的数量和分布范围，利用自检验机制将多幅高分辨率的场图进行合并，排除了最高分辨率下由于相位因子值跳变造成的偏差，保证了最终种子点选择的正确性，从而能够对两者成分的信号进行准确的识别。
[0068] 下面通过具体的实施例对所述的磁共振化学位移编码成像方法进行详细说明：
[0069] 实施例一
[0070] 含有两种成分(本实施例以水和脂肪为例)的两点MR(Magnetic Resonance，磁共振)信号模型为：
[0071]
[0072] 其中，Sn表示在回波时间TEn下的信号强度，n＝1,2表示在本实施例中采集2个回波时间的MR信号；ρw和ρf分别表示水和脂肪的强度值；水的化学位移设定为0，fF表示脂肪相对于水的化学位移，本实施例中可设定为-3.5ppm；fB表示静磁场B0的局部场。定义参数c1、c2、d0、和d如下：
[0073]
[0074]
[0075]
[0076]
[0077] 其中，ΔTE＝TE1-TE2，d0表示在TE1时间内由于fB造成的水的相位，d表示在ΔTE时间内由于fB造成的相位累积，本实施例将d称之为相位因子phasor，其角度与fB成正比，范围是[-ππ]。由上式可知，d是关于fB的周期函数，具有相同d的不同fB在水脂图像重建过程中可得相同的结果，因此在本实施例中，没有估计实际的B0局部场fB，只估计得到相位因子d，其同样表征了B0局部场的不均匀性，亦可称之为场图。
[0078] 如果d已知，则由上式可解得d0：
[0079]
[0080] 上述的MR信号模型可写作矩阵的形式：
[0081]
[0082] 其 中 ，Re表示取复数实部，Im表示取复数虚部。
[0083] 根据变量投影方法(variable projection，VARPRO)，d可使用最小二乘法解得：
[0084] d＝argminerr(d)＝argmin||(I-AA+)S||2
[0085] 上式可以改写为以下形式，并记为式1：
[0086]
[0087] 其中，I表示4×4矩阵；上标‘+’表示取矩阵违逆，A+＝(ATA)-1AT；连接符 表示点乘；err(d)描述了phasor d与拟合误差之间的关系，表示为phasor-error谱，其定义为不同phasor值下对应的模型拟合误差，则相应的局部极小值即为该像素点对应的相位因子值，对于大多数的像素，其相位因子值为两个。
[0088] 将一块感兴趣区域(volume of interest，VOI)内所有像素点作为一个整体的phasor-error谱，则可确定下式，并记为式2：
[0089]
[0090] 其中， 和 均表示行向量，每一个元素代表一个像素点。
[0091] 为了从多个像素点中选取正确的相位因子值，本实施例提出了一种基于多分辨率最安全路径局部增长的自检验场图估计方法，实现了对相位因子值的正确选取，结合图2所示，所述方法包括：
[0092] 步骤21，选取初始种子点。
[0093] 本实施例将处于最高分辨率的二维原始图像进行低通滤波并根据预先定义的欠采样系数对该原始图像进行欠采样，得到不同的低分辨率下的二维原始图像。例如，可将欠采样系数确定为[33]、[55]、[88]和[1010]，则可得到四组低分辨率下的原始图像。对每个低分辨率下的图像中的每一个像素点对应最高分辨图像中的多个像素点，将该多个像素点作为整体并根据上述的式1和式2描述phasor-error谱，即可获得在低分辨下的每一个像素点的phasor-error谱，根据该phasor-error谱寻找使该像素点达到局部极小值的相位因子值，如果该像素点满足以下两个条件则被标记为种子点：
[0094] ①具有大于预定值的信噪比：m>th，其中m表示该像素点的幅值，th表示信噪比阈值。
[0095] ②具有唯一的相位因子值：以Amax和Amin分别表示phasor-error谱的最大值和最小值，以α表示拟合误差倍数，则可设定阈值为Amin+α*(Amax-Amin)，如果只有一个局部极小值使得拟合误差小于该阈值，则认为对应的像素点具有唯一的相位因子值；否则同时保留该像素点其它小于该阈值的局部解，则该像素点具有多个相位因子值。
[0096] 其中，像素点的幅值m可由多种计算方式得到，例如：
[0097] m＝max{abs(S1),abs(S2)}或m＝abs(S1)+abs(S2)等。
[0098] 对于th和α的选取既要排除SNR低和具有多个极小值解的种子点，又要保证得到足够多的种子点，一个候选取值可为：th取所有像素幅值的第98个百分位值的0.2倍，α取0.4。
[0099] 图3所示的为一个种子点的phasor-error谱，该种子点只有一个相位因子值。对于其余相位因子值多于一个的像素点(本实施例称之为待估像素点，其多个相位因子值合称为相位因子候选值)则需要从以下的步骤选取出正确的解。
[0100] 作为可选的，本步骤对上述四组低分辨图像进行种子点的选取。需要说明的是，上述的四组低分辨图像仅为举例，其它低分辨率若可适用于本发明，也应包含在本发明的保护范围以内，并以引用方式包含于此。
[0101] 步骤22，通过最安全路径局部增长算法确定待估像素点的相位因子值。
[0102] 在初始种子点的基础上，本实施例采用如下的局部增长算法完成待估像素点的相位因子估计，其过程可以包括：从待估像素点的相位因子候选值中选取与相邻种子点的相位因子值具有最大相似性的候选值作为当前待估像素点的相位因子值，即：
[0103] d＝argmax{D1,D2,...,DS}，其中Ds表示phasor相似性，定义为：
[0104]
[0105] 其中，K表示当前待估像素点的八邻域(若将当前待估像素点及相邻域划分为九宫格的格式，则除当前待估像素点所在区域以外的八个相邻的区域即为八邻域)内的种子点个数；mk表示第k个种子点的幅值，相应的相位因子值为dk，ds表示当前待估像素点的一个相位因子候选值；conj(.)表示取复共轭；angle(.)表示取角度。
[0106] 可选的，结合图4所示，本实施例提出的最安全路径局部增长算法如下：
[0107] 步骤221，从预定幅值图M中计算获得每个待估像素点的幅值；
[0108] 步骤222，设定幅值阈值thm，寻找幅值大于thm的待估像素点建立集合，该集合记为U；
[0109] 步骤223，计算集合U中所有待估像素点相邻种子点中两两之间的phasor相位差，如果最大相位差小于阈值thp，则保留相应的待估像素点，将保留的待估像素点建立集合，该集合记为V；如果最大相位差大于阈值thp，则剔除相应的待估像素点；
[0110] 步骤224，将集合V中所有待估像素点按照其相邻种子点的个数K进行降序排列；
[0111] 步骤225，选取该降序排列的前Q个待估像素点建立集合，该集合记为L；
[0112] 步骤226，对集合L中所有待估像素点分别计算其与相邻种子点之间最大的phasor相似性；
[0113] 步骤227，如果该相似性大于阈值thD，则选取相应的候选相位因子值作为待估像素点的相位因子值，并将其添加为种子点；否则，该像素点仍为待估像素点；
[0114] 步骤228，重复步骤223至步骤227，直至K＝0；
[0115] 步骤229，减小幅值阈值thm；
[0116] 步骤230，重复步骤222至步骤229，直至幅值阈值thm＝0；
[0117] 步骤231，对于步骤223和步骤227中未处理的待估像素点，其相位因子值选择相邻种子点相位因子值的平均值。
[0118] 上述最安全路径局部增长算法在每一个低分辨率下单独进行。
[0119] 幅值图M中每个像素点的幅值可定义为m＝max{abs(S1),abs(S2)}或m＝abs(S1)+abs(S2)等；幅值阈值thm的选取方法可有多种，例如依据所有像素幅值的百分位值进行选取；阈值thD和阈值thp可以根据候选场图的光滑性确定其数值，例如选取0.9和π/2。待估像素点相邻种子点可定义为其四邻域或八邻域内的种子点；Q值可选取集合V中待估像素点个数的1/4。
[0120] 步骤23，建立最高分辨率的场图。
[0121] 以步骤22中的最安全路径局部增长算法获得的低分辨率场图为初始值，对最高分辨率下的场图进行选取，本实施例中称之为场图映射，其实现过程可以包括：计算最高分辨率下每个像素点的所有相位因子候选值与其对应的低分辨率下的像素点的相位因子值之间的phasor相似性，具有最高相似性的候选值即为当前像素点的相位因子值。
[0122] 本实施例在多个不同的低分辨率分别进行种子点的选取和局部增长算法的估计，因此本步骤中也得到了多个最高分辨率的场图，然后利用自检验机制将上述多个场图进行合并得到一幅场图，其实现过程可以包括：若一个像素点在所有的场图中具有相同的相位因子值，则将该像素点选定为新的种子点；若一个像素点在所有的场图中具有不同的值，则将该像素点重新设定为待估像素点，该待估像素点的相位因子值为相位因子候选值中的一个。
[0123] 利用上述的最安全路径局部增长算法，以新的种子点为起点完成最高分辨率下的待估像素点的相位因子值选取，得到最终的场图。该步骤的原理与上述在低分辨下的局部增长算法相同，因此在本实施例中不再赘述。
[0124] 步骤24，水图和脂肪图的重建。
[0125] 将步骤23获得的最终的场图进行平滑处理，然后通过下式计算即可获得水图和脂肪图：
[0126]
[0127] 其中，W和F分别为计算获得的水图和脂肪图。
[0128] 为了验证本实施例提出的磁共振化学位移编码成像方法的可行性，下面分别在两种体组织实验进了测试，其中的测试对象分别为颈部和腹部。采集序列为FLASH，主磁场B0＝3T，成像参数为(颈部/腹部)：TR＝200ms/50ms，矩阵大小＝240×320/256×152，层厚＝3/6mm，翻转角＝30°/25°，TE＝[4.64 5.51]ms。实验数据的处理环境是具有Intel E5-
2650v2 CPU和64GB RAM的工作站，采用的数据处理软件是MATLAB。根据最终的试验结果可以确定，在体实验中得到的水图和脂肪图中无明显分离错误的组织。图7和图8所示的是根据颈部和腹部的其中一层的分离结果，可以看出在体组织实验中，得到的水图和脂肪图中无明显分离错误的组织。
[0129] 需要说明的是，本发明提出的磁共振化学位移编码成像方法不仅可以用于水脂分离成像，也可以用于分离其它化学位移成份的应用中，其中的区别仅仅在于含有两种成分的简化MR信号模型中输入对应的化学位移值。
[0130] 本实施例提出的磁共振化学位移编码成像方法不对TE施加任何的限制，能够达到更快的图像采集速度，使其适用于一些快速成像的医学应用，如动态成像和腹部屏气成像中；提出的最安全路径局部增长算法相比于已有的局部增长算法，可在低信噪比的区域也具有较高的鲁棒性；而利用多分辨率进行场图估计和局部增长，本实施例提出的方法可在低分辨率下将同时含有水和脂肪的区域识别为种子点，增加了种子点的数量和分布范围；低分辨下的场图估计和局部增长是相互独立的，即各个低分辨率之间的处理没有任何交叠，这就消除了不同分辨率之间可能存在的相位因子值跳变的影响；利用自检验机制将多幅高分辨率的场图进行合并，排除了最高分辨率下由于相位因子值跳变造成的偏差，保证了最终种子点选择的正确性。
[0131] 本发明还提出了一种磁共振化学位移编码成像装置，结合图5所示，所述装置包括：
[0132] 相位误差谱建立模块51，用于在预定分辨率下的采样图像中，基于两点磁共振信号模型建立phasor-error谱；
[0133] 种子点选取模块52，用于将phasor-error谱中具有唯一的相位因子值且使该谱达到局部极小值的像素点确定为初始种子点；
[0134] 估计模块53，用于根据初始种子点对待估像素点的相位因子值进行估计，获得在预定分辨率下的场图；
[0135] 场图重构模块54，用于将预定分辨率下的场图在最高分辨率下分别进行映射，获得在最高分辨率下的多幅场图，并将多幅场图合并后获得重构场图；
[0136] 相位因子确定模块55，用于从重构场图中确定重构种子点，并通过对重构种子点进行估计以获得重构待估像素点的相位因子值；
[0137] 成像模块56，用于根据重构种子点和重构待估像素点的相位因子值得到两种预定成分的分离图像。
[0138] 可选的，在种子点选取模块中52中包括：
[0139] 极小值确定子模块，用于获得所述预定分辨率图像中的每一个像素点的phasor-error谱，并根据所述phasor-error谱确定使所述像素点达到局部极小值的相位因子值；
[0140] 种子点确定子模块，用于若所述像素点的信噪比大于预定值且具有唯一的相位因子值，则将所述像素点确定为初始种子点。
[0141] 其中，该种子点确定子模块可用于若待估像素点包括有多个局部极小值，则将拟合误差小于预定值的局部极小值确定为待估像素点的多个相位因子值，并从待估像素点的多个相位因子值中选取与相邻种子点的相位因子值具有最大相似性的相位因子值作为待估像素点的相位因子值。
[0142] 可选的，在估计模块53中包括：
[0143] 第一集合建立子模块，用于选取幅值大于第一阈值的待估像素点建立第一待估像素点集合；
[0144] 第二集合建立子模块，用于获取所述第一待估像素点集合中每个待估像素点的相邻种子点中两两之间的相位差，选取最大相位差小于第二阈值的待估像素点建立第二待估像素点集合；
[0145] 第三集合建立子模块，用于将所述第二待估像素点集合中的待估像素点按所述相邻种子点的个数进行降序排列，选取预定数量的待估像素点建立第三待估像素点集合；
[0146] 第一相位因子确定子模块，用于获取所述第三待估像素点集合中的待估像素点与所述相邻种子点之间最大的相位相似性，并将所述相似性大于第三阈值的待估像素点的相位因子值作为待估像素点的相位因子值；
[0147] 第一估计子模块，用于将所述相似性大于第三阈值的待估像素点作为新的种子点对所述第一待估像素点集合中的待估像素点进行估计，直至所述第一待估像素点集合中的待估像素点的相邻种子点数为零；
[0148] 第四集合建立模块，用于通过将所述第一阈值减小预定量以确定第四待估像素点集合；
[0149] 第二相位因子确定子模块，用于确定所述第四待估像素点集合中的待估像素点的相位因子值，直至所述第一阈值减为零；
[0150] 第三相位因子确定子模块，用于将所述第一待估像素点集合中未确定相位因子值的待估像素点的相位因子值确定为相邻种子点的相位因子值的平均值。
[0151] 可选的，在场图重构模块54中包括：
[0152] 第四相位因子确定子模块，用于获取最高分辨率下每一个重构待估像素点的所有相位因子值与对应的预定分辨率下的待估像素点的相位因子值之间的相位相似性，将最高相似性的相位因子值确定为所述重构待估像素点的相位因子值。
[0153] 可选的，在相位因子确定模块55中包括：
[0154] 第五相位因子确定子模块，用于若所述重构待估像素点在所有的场图中均具有相同的相位因子值，则将所述重构待估像素点确定为新的种子点以及将所述相同的相位因子值确定为所述重构待估像素点的相位因子值；
[0155] 第六相位因子确定子模块，用于若所述重构待估像素点在所有的场图中具有不同的相位因子值，则将所述重构待估像素点重新设定为待估像素点以及将所述不同的相位因子值确定为相位因子候选值；
[0156] 第七相位因子确定子模块，用于以所述新的种子点为起点在最高分辨率下对所述相位因子候选值进行估计，以确定所述重构待估像素点的相位因子值。
[0157] 采用本实施例提出的磁共振化学位移编码成像装置，通过在低分辨率下进行场图估计，可在低分辨率下将同时含有两种成分的区域识别为种子点，增加了种子点的数量和分布范围，利用自检验机制将多幅高分辨率的场图进行合并，排除了最高分辨率下由于相位因子值跳变造成的偏差，保证了最终种子点选择的正确性，从而能够对两者成分的信号进行准确的识别。
[0158] 本发明还提出了一种设备，该设备包括处理器和存储有计算机可读指令的存储器，该处理器被配置来通过执行该计算机可读指令以执行如实施例一所述的方法。
[0159] 其中，图6是能够执行图4所示的指令流程并以图5的示例性磁共振化学位移编码成像装置的示例性处理器平台600的框图，该处理器平台600可以是服务器、个人计算机或任何其它类型的计算机设备。
[0160] 图6所示的处理器平台600包括处理器612。例如，处理器612可以由来自任何期望家庭或制造商的一个或者更多个微处理器或控制器来实现。
[0161] 处理器612包括本地存储器613(例如，高速缓冲存储器)，并且经由总线618与主存储器进行通信，主存储器包括易失性存储器614和非易失性存储器616。易失性存储器614可以由同步动态随机存取存储器(SDRAM)、动态随机存取存储器(DRAM)、RAM总线(RAMBUS)动态随机存取存储器(RDRAM)和/或任何其它类型的随机存取存储器装置来实现。非易失性存储器616可以由闪存和/或任何其它期望类型的存储设备来实现。对易失性存储器614和非易失性存储器616的访问是由存储器控制器来控制的。
[0162] 处理器平台600还可以包括一个或更多个大容量存储装置621，用于存储软件和数据。这样的大容量存储装置621的示例包括软盘驱动器、硬盘驱动器、光盘驱动器和数字多功能光盘(DVD)驱动器。图6所示的编码指令631可以被存储在大容量存储装置621、易失性存储器614、非易失性存储器616和/或可移动存储介质如CD或DVD上。
[0163] 本具体实施方式是对本发明的技术方案进行清楚、完整地描述，其中的实施例仅仅是本发明的一部分实施例，而并不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有经过创造性劳动的前提下所获得的所有其它实施方式都属于本发明的保护范围。
