本发明公开了一种磁共振系统及其水脂分离成像方法、装置，包括：采集患者身体组织中的水和脂肪的磁共振图像；利用磁共振图像中的至少两幅反相位或至少两幅同相位的磁共振图像计算出附加相位；对附加相位进行去卷绕，所述去卷绕步骤包括：采用多项式拟合法对附加相位进行去卷绕；在多项式拟合法去卷绕的结果上采用区域增长法进行去卷绕处理；将去卷绕后的磁共振图像中的反相位的图像结合同相位的图像得出水图像和脂肪图像。本发明采用多项式拟合法和区域增长法相结合的方法对附加相位进行去卷绕，有效去除了多项式拟合法去卷绕产生的孤岛效应。
1.磁共振系统的水脂分离成像方法,其特征在于包括：
采集患者身体组织中的水和脂肪的磁共振图像，所述磁共振图像包括至少一幅同相位的磁共振图像和至少两幅反相位的磁共振图像，或所述磁共振图像包括至少一幅反相位的磁共振图像和至少两幅同相位的磁共振图像；
利用至少两幅反相位或至少两幅同相位的磁共振图像计算出附加相位；
对附加相位进行去卷绕，所述对附加相位进行去卷绕步骤包括：采用多项式拟合法对附加相位进行去卷绕；在多项式拟合法去卷绕的结果上采用区域增长法进行去卷绕处理；
将去卷绕后的磁共振图像中的反相位的图像结合同相位的图像得出水图像和脂肪图像。
2.如权利要求1所述的方法，其特征在于，采用区域增长法去卷绕步骤包括：
选取种子步骤，在未增长图像点中选择幅度最大且邻点为附加相位不卷绕的点作为种子，对于初始种子，若未找到符合要求的图像点，则选择未增长图像点中幅度次大且邻点为附加相位不卷绕的点作为种子；
种子幅度判断步骤，将新选取的种子的幅度与图像点中的最大幅度值进行比较，如果新选取的种子的幅度小于最大幅度值的1％或其它任一个可认为是远远小于的值时则结束增长，否则执行以下步骤；
种子相位调整步骤，将新选取的种子的相位与其最邻近的已增长图像点的相位比较，如果相位差大于2π则对该新选取的种子的相位进行加/减2π整数倍的相位调整，以使两者的相位差最小，否则将新选取的种子作为增长的种子；
增长去卷绕步骤，用于将增长的种子的相位与其周边相邻图像点的相位依次进行比较，若比较结果小于设定的第二相位阈值，则该相邻图像点被增长，并将该相邻图像点列为下一轮增长的种子之一；若比较结果大于或等于设定的第二相位阈值，则对该相邻图像点进行增/减2π整数倍后再次判断新的比较结果是否小于设定的第二相位阈值，如果是则将增/减2π整数倍后的相位作为该相邻图像点的新相位，并将该相邻图像点列为下一轮增长的种子之一，循环本步骤直到没有相邻图像点被增长；
判断是否所有图像点都已被增长，若是则结束增长，否则转向选取种子步骤。
3.如权利要求2所述的方法，其特征在于，在采用多项式拟合法和区域增长法进行去卷绕处理之前，所述对附加相位进行去卷绕步骤还包括：选取任一幅图像，对图像进行信号判断，识别出信号点；
在采用多项式拟合法对附加相位进行去卷绕步骤中仅对识别出的信号点的附加相位进行去卷绕，所述采用区域增长法去卷绕步骤中的图像点为信号点。
4.如权利要求3所述的方法，其特征在于，所述对图像进行信号判断包括：
幅度判断步骤，将图像中的像素的幅度和最终幅度阈值进行比较，如果像素的幅度大于或等于最终幅度阈值则认为该像素是信号点，如果像素的幅度小于最终幅度阈值则认为该像素不是信号点。
5.如权利要求4所述的方法，其特征在于，所述最终幅度阈值通过以下步骤选取：
计算初始幅度阈值的步骤，所述初始幅度阈值等于图像中的最大幅度值与最小幅度值的平均值；
像素分割步骤，将图像中所有像素的幅度值和初始幅度阈值进行比较，得到第一组像素集合和第二组像素集合，第一组像素集合为大于或等于初始幅度阈值的像素集合，第二组像素集合为小于初始幅度阈值的像素集合；
中间幅度阈值计算步骤，分别计算第一组像素集合的平均幅度值和第二组像素集合的平均幅度值，将第一组像素集合的平均幅度值和第二组像素集合的平均幅度值的平均值作为当前中间幅度阈值；
最终幅度阈值识别步骤，将该当前中间幅度阈值与其前一个中间幅度阈值进行比较，判断比较结果是否小于预先设定的参数，若是则该当前中间幅度阈值为最终幅度阈值，否则转向执行像素分割步骤，将图像中所有像素的幅度值和该当前中间幅度阈值进行比较。
6.如权利要求4或5所述的方法，其特征在于，所述对图像进行信号判断还包括：
相位判断步骤，将相邻像素的相位进行比较，如果两者的比较结果小于相位阈值，则认为相邻像素为信号点。
7.如权利要求6所述的方法，其特征在于，所述相位判断步骤包括：
选取初始种子步骤，选择图像中幅度最大的像素点为初始种子；
种子增长步骤，将增长的种子与其周边相邻像素点依次进行相位比较，如果相位比较结果小于第一相位阈值，则该相邻像素点被增长，并将该相邻像素点作为下一轮增长的种子之一，循环本步骤直到增长到没有种子，循环结束后种子的集合即为信号点的集合。
8.如权利要求7所述的方法，其特征在于，所述相位判断步骤还包括：
将种子增长步骤后得出的信号点个数与幅度判断步骤中得出的信号点个数进行比较，判断比较结果是否符合第一预定条件，若是，则以相位判断步骤中得出的信号点为最终结果，否则执行以下步骤；
按照设定步长将第一相位阈值进行递增后作为新的第一相位阈值，循环执行选取初始种子步骤和种子增长步骤直到相位判断步骤中得出的信号点个数与幅度判断步骤中得出的信号点个数的比较结果符合第一预定条件。
9.磁共振系统的水脂分离成像装置,其特征在于包括：
图像采集模块，用于采集患者身体组织中的水和脂肪的磁共振图像，所述磁共振图像包括至少一幅同相位的磁共振图像和至少两幅反相位的磁共振图像，或所述磁共振图像包括至少一幅反相位的磁共振图像和至少两幅同相位的磁共振图像；
附加相位计算模块，用于利用至少两幅反相位或至少两幅同相位的磁共振图像计算出附加相位；
去卷绕模块，用于对附加相位进行去卷绕，所述去卷绕模块包括拟合单元和第一区域增长单元，所述拟合单元用于采用多项式拟合法对附加相位进行去卷绕；所述第一区域增长单元用于在多项式拟合法去卷绕的结果上采用区域增长法进行去卷绕处理；
水脂分离模块，用于将去卷绕后的磁共振图像中的反相位的图像结合同相位的图像得出水图像和脂肪图像。
10.如权利要求9所述的装置，其特征在于，所述第一区域增长单元包括：
选取种子子单元，用于在未增长图像点中选择幅度最大且邻点为附加相位不卷绕的点作为种子，对于初始种子，若未找到符合要求的图像点，则选择未增长图像点中幅度次大且邻点为附加相位不卷绕的点作为种子；
种子幅度判断子单元，用于将新选取的种子的幅度与图像点中的最大幅度值进行比较，如果新选取的种子的幅度小于最大幅度值的1％或其它任一个可认为是远远小于的值时则结束增长；
种子相位调整子单元，用于当新选取的种子的幅度没有远远小于图像点中的最大幅度时将新选取的种子的相位与其最邻近的已增长图像点的相位比较，如果相位差大于2π则对该新选取的种子的相位进行加/减2π整数倍的相位调整，以使两者的相位差最小，否则将新选取的种子作为增长的种子；
增长去卷绕子单元，用于将增长的种子的相位与其周边相邻图像点的相位依次进行比较，若比较结果小于设定的第二相位阈值，则该相邻图像点被增长，并将该相邻图像点列为下一轮增长的种子之一；若比较结果大于或等于设定的第二相位阈值，则对该相邻图像点进行增/减2π整数倍后再次判断新的比较结果是否小于设定的第二相位阈值，如果是则将增/减2π整数倍后的相位作为该相邻图像点的新相位，并将该相邻图像点列为下一轮增长的种子之一，控制增长去卷绕子单元循环工作直到没有相邻图像点被增长；
信号点增长判断子单元，用于判断是否所有图像点都已被增长，若是则结束增长，否则控制选取种子子单元工作，继续在未增长图像点中选取种子。
11.如权利要求10所述的装置，其特征在于，所述去卷绕模块还包括：
信号点识别单元，用于选取任一幅图像，对图像进行信号判断，识别出信号点；
所述拟合单元用于在采用多项式拟合法对附加相位进行去卷绕时仅对信号点的附加相位进行去卷绕，所述第一区域增长单元中的图像点为信号点。
12.如权利要求11所述的装置，其特征在于，所述信号点识别单元包括：
幅度判断子单元，用于将图像中的像素的幅度和最终幅度阈值进行比较，如果像素的幅度大于或等于最终幅度阈值则认为该像素是信号点，如果像素的幅度小于最终幅度阈值则认为该像素不是信号点。
13.如权利要求12所述的装置，其特征在于，所述幅度判断子单元用于先计算初始幅度阈值，所述初始幅度阈值等于图像中的最大幅度值与最小幅度值的平均值；然后将图像中所有像素的幅度值和初始幅度阈值进行比较，得到第一组像素集合和第二组像素集合，第一组像素集合为大于或等于初始幅度阈值的像素集合，第二组像素集合为小于初始幅度阈值的像素集合；然后分别计算第一组像素集合的平均幅度值和第二组像素集合的平均幅度值，将第一组像素集合的平均幅度值和第二组像素集合的平均幅度值的平均值作为当前中间幅度阈值；将该当前中间幅度阈值与其前一个中间幅度阈值进行比较，判断比较结果是否小于预先设定的参数，若是则该当前中间幅度阈值为最终幅度阈值，否则像素分割子单元工作，将图像中所有像素的幅度值和该当前中间幅度阈值进行比较。
14.如权利要求12或13所述的装置，其特征在于，所述信号点识别单元还包括：
相位判断子单元，用于将相邻像素的相位进行比较，如果两者的比较结果小于相位阈值，则认为相邻像素为信号点。
15.如权利要求14所述的装置，其特征在于，所述相位判断子单元包括：
选取初始种子子单元，用于选择图像中幅度最大的像素点为初始种子；
种子增长子单元，用于将增长的种子与其周边相邻像素点依次进行相位比较，如果相位比较结果小于第一相位阈值，则该相邻像素点被增长，并将该相邻像素点作为下一轮增长的种子之一，控制种子增长子单元工作直到增长到没有种子，循环结束后种子的集合即为信号点的集合。
16.如权利要求15所述的装置，其特征在于，所述相位判断子单元还包括：
信号点个数判断子单元，用于将种子增长子单元得出的信号点个数与幅度判断子单元得出的信号点个数进行比较，判断比较结果是否符合第一预定条件，若是，则以种子增长子单元得出的信号点为最终结果；
第一相位阈值递增子单元，用于当种子增长子单元得出的信号点个数与幅度判断子单元得出的信号点个数不符合第一预定条件时，按照设定步长将第一相位阈值进行递增后作为新的第一相位阈值，控制选取初始种子子单元和种子增长子单元工作直到种子增长子单元得出的信号点个数与幅度判断子单元得出的信号点个数的比较结果符合第一预定条件。
17.一种磁共振系统，其特征在于包括如权利要求9至16中任一项所述的水脂分离成像装置。
磁共振系统及其水脂分离成像方法、装置
技术领域
[0001] 本发明涉及一种磁共振系统，尤其涉及磁共振系统的水脂分离成像方法、装置。
背景技术
[0002] 在磁共振成像中，人体组织中的脂肪中的氢质子和其它非脂肪中的氢质子所处的分子环境不同，使得它们的共振频率存在微小差异，称为化学位移。Dixon最早提出利用化学位移效应获得磁共振成像中的水图像和脂肪图像。Dixon方法的原理是假设人体中的每个体素只有水和脂肪两种成分，因而磁共振图像中的每个像素点由水的宏观磁化矢量和脂肪的宏观磁化矢量合成。根据化学位移原理，当水和脂肪中的氢质子同时受到射频脉冲激励后，它们的弛豫时间也不一样，在不同的回波时间采集信号，水和脂肪表现出不同的信号强度。因而在常规自旋回波序列基础上，通过调节不同的回波时间，可以在理论上精确地获得水和脂肪两个磁化矢量的任意夹角。
[0003] 而实际工程应用中由于系统噪声、涡流效应、主磁场不均匀性、Maxwell伴随场以及人体组织之间、组织与空气之间磁导率变化等都会导致磁场不均匀，磁场的不均匀性又会引起相位的严重畸变，因此如何消除这些误差的影响成为Dixon算法成败的关键。
[0004] 作为对Dixon方法的改进，一种基于不均匀场的相位误差校正的三点法被提出，该方法的原理是取得水和脂肪同相位的一副磁共振图像和反相位的两幅磁共振图像，根据反相位的两幅磁共振图像计算出磁场不均匀性导致的附加相位，然后与同相位的图像一起求得水图像和脂肪图像。
[0005] 该三点法的不足之处在于：由两幅反相位的图像求得的附加相位不能直接用来进行水脂分离，因为由磁场不均匀导致的附加相位可能超过一个域，例如超过±π，即发生相位卷绕现象，因此必须对相位进行去卷绕。为了对相位进行去卷绕，一种方案是假设磁场不均匀在空间上的变化是缓慢的，采用单一的多项式拟合的方法进行相位去卷绕，但这种方案求解不准确，容易导致孤岛效应。
发明内容
[0006] 本发明要解决的主要技术问题是，提供一种磁共振系统及其水脂分离成像方法、装置，更准确地进行相位去卷绕，从而更准确地生成分离的水图像和脂肪图像。
[0007] 根据本发明的一方面，提供一种磁共振系统的水脂分离成像方法，包括：
[0008] 采集患者身体组织中的水和脂肪的磁共振图像，所述磁共振图像包括至少一副同相位的磁共振图像和至少两幅反相位的磁共振图像，或所述磁共振图像包括至少一副反相位的磁共振图像和至少两幅同相位的磁共振图像；
[0009] 利用至少两幅反相位或至少两幅同相位的磁共振图像计算出附加相位；
[0010] 对附加相位进行去卷绕，所述去卷绕步骤包括：采用多项式拟合法对附加相位进行去卷绕；在多项式拟合法去卷绕的结果上采用区域增长法进行去卷绕处理；
[0011] 将去卷绕后的磁共振图像中的反相位的图像结合同相位的图像得出水图像和脂肪图像。
[0012] 根据本发明的另一方面，提供一种磁共振系统的水脂分离成像装置，包括：
[0013] 图像采集模块，用于采集患者身体组织中的水和脂肪的磁共振图像，所述磁共振图像包括至少一副同相位的磁共振图像和至少两幅反相位的磁共振图像，或所述磁共振图像包括至少一副反相位的磁共振图像和至少两幅同相位的磁共振图像；
[0014] 附加相位计算模块，用于利用至少两幅反相位或至少两幅同相位的磁共振图像计算出附加相位；
[0015] 去卷绕模块，用于对附加相位进行去卷绕，所述去卷绕模块包括拟合单元和第一区域增长单元，所述拟合单元用于采用多项式拟合法对附加相位进行去卷绕；所述第一区域增长单元用于在多项式拟合法去卷绕的结果上采用区域增长法进行去卷绕处理；
[0016] 水脂分离模块，用于将去卷绕后的磁共振图像中的反相位的图像结合同相位的图像得出水图像和脂肪图像。
[0017] 本发明还提供一种磁共振系统，包括上述的水脂分离成像装置。
附图说明
[0018] 图1为本发明一种实施例中磁共振系统的结构示意图；
[0019] 图2为本发明一种实施例中水脂分离成像装置的结构框图；
[0020] 图3为本发明一种实施例中第一区域增长单元的结构框图；
[0021] 图4为本发明一种实施例中去卷绕模块的结构框图；
[0022] 图5为本发明一种实施例中水脂分离成像的流程图；
[0023] 图6为本发明一种实施例中区域增长法去卷绕的流程图；
[0024] 图7为本发明另一种实施例中水脂分离成像的流程图；
[0025] 图8为本发明一种实施例中信号点识别的流程图。
具体实施方式
[0026] 下面通过具体实施方式结合附图对本发明作进一步详细说明。
[0027] 请参考图1，在一种实施例中，磁共振系统100包括磁体系统110、梯度磁场系统120、射频系统130和控制及处理系统140。磁体系统110包括磁体111、梯度磁场线圈112、发射线圈113和接收线圈114，磁体111可以采用永磁体、超导磁体或常导磁体，用于给待测物体(例如病人)提供一恒定的主磁场，梯度磁场线圈112用于在三维空间产生一梯度磁场，发射线圈113用于提供射频(RF)脉冲以激发待测物体内原子核的自旋，接收线圈114用于检测由待测物发出的回波信号。梯度磁场系统120和控制及处理系统140连接，用于在控制及处理系统的控制下驱动梯度磁场线圈112。射频系统130和控制及处理系统140连接，用于在控制及处理系统的控制下产生RF脉冲并经放大处理后施加给发射线圈113。
[0028] 控制及处理系统140既用于对各部分进行控制，也用于对回波信号进行处理。将接收线圈114检测到的回波信号传输到控制及处理系统140，控制及处理系统140包括水脂分离成像装置，根据回波信号得到磁共振图像，对磁共振图像进行处理后获得分离的水图像和脂肪图像。
[0029] 在本发明实施例中，水脂分离的原理是：采用反相位的两幅图像和同相位的一幅图像或同相位的两幅图像和反相位的一幅图像，如果不存在磁共振系统磁场不均匀导致的误差，则通过反相位或同相位的两幅图像相加减，即可得到水的图像和脂肪的图像。但由于磁共振系统磁场不均匀或其他原因，会在水宏观磁化矢量和脂肪的宏观磁化矢量上产生一附加相位，通过反相位或同相位的两幅图像可求得附加相位，对附加相位去卷绕后再将反相位的图像结合同相位的图像得出水图像和脂肪图像。本发明实施例在对附加相位去卷绕时包括两步：即先采用多项式拟合法对附加相位进行去卷绕，然后在多项式拟合法去卷绕的结果上采用区域增长法进行去卷绕处理。通过采用区域增长法进行去卷绕，有效去除了多项式拟合法去卷绕产生的孤岛效应。
[0030] 图2所示为一种实施例中水脂分离成像装置的结构框图，水脂分离成像装置200包括图像采集模块210、附加相位计算模块220、去卷绕模块230和水脂分离模块240。图像采集模块210用于采集患者身体组织中的水和脂肪的磁共振图像，例如采集至少一副同相位的磁共振图像和至少两幅反相位的磁共振图像，或采集至少一副反相位的磁共振图像和至少两幅同相位的磁共振图像。附加相位计算模块220用于利用至少两幅反相位或同相位的磁共振图像计算出附加相位。去卷绕模块230用于对附加相位进行去卷绕，所述去卷绕模块230包括拟合单元250和第一区域增长单元260，所述拟合单元250用于采用多项式拟合法对附加相位进行去卷绕；所述第一区域增长单元260用于在多项式拟合法去卷绕的结果上采用区域增长法进行去卷绕处理。水脂分离模块240用于将去卷绕后的磁共振图像中的反相位的图像结合同相位的图像得出水图像和脂肪图像，例如将去卷绕后的两幅反相位的图像结合一幅同相位的图像得出水图像和脂肪图像，或将去卷绕后的两幅同相位的图像结合一幅反相位的图像得出水图像和脂肪图像。
[0031] 在如图3所示实施例的一种具体实例中，所述第一区域增长单元260包括选取种子子单元261、种子幅度判断子单元262、种子相位调整子单元263、增长去卷绕子单元264和信号点增长判断子单元265。选取种子子单元261用于在确定的信号点中选取种子，并将选取的种子送到种子幅度判断子单元262进行判断，如果选取的种子符合要求，则将种子送到种子相位调整子单元263，根据预定规则对种子的相位进行调整，并将调整后的相位作为种子的相位。将调整相位后的种子送到增长去卷绕子单元264按照预定规则进行区域增长。在增长去卷绕子单元264增长完后信号点增长判断子单元265判断是否所有信号点都已被增长，若是则结束增长，否则控制选取种子子单元261工作，继续在未增长信号点中选取种子。
[0032] 在改进的实施例中，如图4所示，去卷绕模块230包括拟合单元250、第一区域增长单元260和信号点识别单元270，信号点识别单元270用于选取任一幅图像，对图像进行信号判断，先识别出信号点。所述拟合单元250用于在采用多项式拟合法对附加相位进行去卷绕时仅对信号点识别单元270识别出的信号点的附加相位进行去卷绕。
[0033] 其中，信号点识别单元270包括幅度判断子单元271，幅度判断子单元271用于将图像中的像素的幅度和最终幅度阈值进行比较，如果像素的幅度大于或等于最终幅度阈值则认为该像素是信号点，如果像素的幅度小于最终幅度阈值则认为该像素不是信号点。在另一实施例中，信号点识别单元270还包括相位判断子单元272，相位判断子单元272用于将相邻像素的相位进行两两比较，如果两者的比较结果小于相位阈值，则认为相邻像素为信号点。在一种实施例中，相位判断子单元272包括选取初始种子子单元272a、种子增长子单元272b、信号点个数判断子单元272c和第一相位阈值递增子单元272d。初始种子子单元272a用于根据设定条件选择初始种子，然后将获得的初始种子传输给种子增长子单元272b，种子增长子单元272b以初始种子为起始点，按照设定的第一相位阈值开始区域增长，被增长的像素点即认为是信号点。信号点个数判断子单元272c在种子增长子单元272b每完成循环得到信号点个数后，判断采用相位判断法得到的信号点个数是否符合要求，如果不符合要求则控制第一相位阈值递增子单元272d改变第一相位阈值，然后控制初始种子子单元272a和种子增长子单元272b再次工作。
[0034] 基于上述水脂分离成像装置，下面详细说明磁共振系统的水脂分离成像方法。
[0035] 在一种实施例中，如图5所示，水脂分离成像方法包括以下步骤：
[0036] 步骤S1，分别采集患者身体组织中的水和脂肪同相位和反相位的磁共振图像，同相位与反相位是指采集得到的信号中水信号和脂肪信号之间的角度，例如采集一副同相位磁共振图像和两幅反相位的磁共振图像，水和脂肪相位角度为(-π，0，π)。当然，也可以采集多幅同相位的磁共振图像和更多反相位的磁共振图像。由于磁共振系统的磁场不均匀或其他不良原因，将导致图像存在一个附加相位，附加相位用φ0表示，则采集的水和脂肪相对角度为(-π，0，π)三幅磁共振图像如下：
[0037]
[0038] 步骤S2，利用两幅反相位的磁共振图像计算出附加相位。由I1可以得到 ，用该值处理(1)式便得到：
[0039]
[0040] 由(2)式的两幅反相位的I0′以及I2′相除，即可得到附加相位，表达式如下：
[0041] ei2φ＝I′2*conj(I′0)----------------------------(3)
[0042] 步骤S4，对附加相位进行去卷绕。由于磁场不均匀导致的附加相位可能超过一个域，例如超过±π，即发生相位卷绕现象，因此必须对式(3)中的2φ角度去卷绕，所述去卷绕步骤包括：采用多项式拟合法对附加相位进行去卷绕；在多项式拟合法去卷绕的结果上采用区域增长法进行去卷绕处理。
[0043] 在采用多项式拟合法对附加相位进行去卷绕时，可采用已有的多项式进行拟合，多项式拟合法相位去卷绕的前提是假设图像中的信号点的相位是缓慢变化的，且其可以被多阶多项式拟合得到。本实施例中采用二次多项式拟合公式为：
[0044] 2φw＝ax+bx2+cy+dy2+exy+f -------------(4)
[0045] 上式中x、y是二维图像中的坐标，a、b、c、d、e、f是拟合系数。具体的实现过程是：
[0046] (1)先用atan函数得到卷绕的2φw相位：
[0047] 2φw＝atan(I3)--------------(5)
[0048] (2)对式(4)分别对x、y求导得到式(6)
[0049] -------------------(6)
[0050]
[0051] (3)用数据拟合式(6)前，先用阈值T0＝π/2剔除求导后数据中的尖峰值(目的是剔除2φw中的不连续值)，然后进行多项式拟合，拟合采用的是最小二乘法。其中f＝mean(2φw)(图像信号点相位的平均值)。
[0052] (4)得到拟合值后，计算of＝2φw-(ax+bx2+cy+dy2+exy+f)，然后根据该相位差of对信号点的卷绕相位加减2π的整数倍使of达到最小。
[0053] 采用多项式拟合法对附加相位进行去卷绕后得到的相位图有块状的不连续区域，即产生孤岛效应，因此在多项式拟合法对附加相位进行去卷绕的结果上，再采用区域增长法对附加相位进行去卷绕处理。
[0054] 区域增长法去卷绕的流程如图6所示，包括以下步骤：
[0055] 步骤S41，选取种子步骤。选取种子子单元261在未增长信号点中选择幅度最大且邻点为附加相位不卷绕的点作为种子。对于初始种子，若未找到符合要求的信号点，则选择未增长信号点中幅度次大且邻点为附加相位不卷绕的点作为种子。
[0056] 步骤S42，种子幅度判断步骤。种子幅度判断子单元262将新选取的种子的幅度与信号点中的最大幅度值进行比较，如果比较结果符合第二预定条件则结束增长，不符合第二预定条件则进行步骤S43。在一种具体实例中，计算新选取的种子的幅度与信号点中的最大幅度值的差值或比值，当新选取的种子的幅度远远小于信号点中的最大幅度值时，例如新选取的种子的幅度小于最大幅度值的1％或其它任一个可认为是远远小于的值时，认为新选取的种子不可能是需要的信号点，所以结束增长。当然，本领域技术人员可以理解，第二预定条件还可以是根据具体情况设定的其他条件。
[0057] 步骤S43，将新选取的种子的相位与其最邻近的已增长信号点的相位比较，如果相位比较结果不符合第三预定条件则执行步骤S44。否则执行步骤S45。本步骤中的相位比较可以是计算新选取的种子的相位与其最邻近的已增长信号点的相位的相位差或比值，优选是计算两者的相位差，第三预定条件也可根据比较的方式不同而根据具体情况和经验进行设定，例如当计算新选取的种子的相位与其最邻近的已增长信号点的相位差时，第三预定条件可以是判断相位差是否大于2π，当相位差大于2π时则执行步骤S44，对种子相位进行调整，否则直接用该种子进行区域增长。
[0058] 步骤S44，对种子相位进行调整。种子相位调整子单元263用于对该新选取的种子的相位进行加/减2π整数倍的相位调整。在对种子的相位进行相位调整时，具体加/减2π的倍数可根据相位差进行选定，目的是使新选取的种子的相位与其最邻近的已增长信号点的相位的相位差最小，例如相位差是5π，则需要对种子的相位进行加/减4π的调整，将调整后的相位作为新选取的种子的相位。如果相位比较结果符合第三预定条件则将新选取的种子作为增长的种子，在下面步骤中用该种子进行区域增长。
[0059] 步骤S45，增长去卷绕步骤。增长去卷绕子单元264用于将增长的种子的相位与其周边相邻信号点的相位依次进行比较，若种子的相位与某一相邻信号点的相位的比较结果小于设定的第二相位阈值，则该相邻信号点被增长，并将该相邻信号点列为下一轮增长的种子之一。本步骤中的相位比较可以是计算两者的相位差值，也可以是计算两者的相位比值，第二相位阈值也可根据比较的方式不同而根据具体情况和经验进行设定。在一种具体实例中，计算种子的相位与其周边相邻信号点的相位的相位差，第二相位阈值GrowT为一经验值，例如GrowT为π/4或π/8等，当相位差小于GrowT时，则该相邻信号点被增长。
[0060] 若增长的种子的相位与某一相邻信号点的相位的比较结果大于或等于设定的第二相位阈值，则对该相邻信号点的相位进行加/减2π整数倍，具体加/减2π的倍数可根据相位差进行选定，目的是使两者的相位差最小，然后再次将种子的相位与某一相邻信号点的加/减2π整数倍后的相位进行比较，判断新的比较结果是否小于设定的第二相位阈值，如果是则将增/减2π整数倍后的相位作为该相邻信号点的新相位，该相邻信号点被增长，并将该相邻信号点列为下一轮增长的种子之一。
[0061] 本步骤中已经用于区域增长(即用于比较)的种子在下一轮增长时将被排除，不再用于下一轮的比较和区域增长。
[0062] 当采用增长的种子比较完后，本步骤再循环一次，采用上次循环中得到的下一轮增长的种子进行相邻信号点的比较，符合条件的信号点被增长，并被列为更下一轮增长的种子之一。如此循环直到没有相邻信号点被增长，然后执行步骤S46。
[0063] 步骤S46，信号点增长判断子单元265判断是否所有信号点都已被增长，若是则结束增长，否则转向步骤S41，继续在未增长信号点中选取种子。
[0064] 在采用本实施例的方案对式(3)中的2φ角度进行去卷绕后，得到：
[0065] p＝coSθ，θ＝arg(eiφ)----------------------------(7)
[0066] 步骤S5，将去卷绕后的两幅反相位的图像结合一副同相位的图像得出水图像和脂肪图像。将I0′以及I2′相乘，并将I′和式(7)代入，得到水图像和脂肪图像如下。
[0067]
[0068] 采用多项式拟合去卷绕的前提是假设图像中信号点的相位是缓慢变化的，拟合后可对卷绕是4π或6π的相位去卷绕，但仍然存在2π的相位卷绕，这些存在相位卷绕的地方会导致最后的结果图像是一块一块的，导致孤岛效应，从而导致后续错误的水脂分离，使脂肪图像上有水信号，水图像上有脂肪信号，不利于医生的诊断。当将多项式拟合去卷绕结果代入到本发明实施例中的区域增长去卷绕方法中时，由于该区域增长去卷绕方法未假设图像中信号点的相位是缓慢变化的，可针对卷绕是2π的相位去卷绕，两种方法的结合正好去除孤岛效应，可提高后续水脂分离的正确性，为医生提供更准确的水图像和脂肪图像。
[0069] 在另一种实施例中，如图7所示，水脂分离成像方法包括以下步骤：
[0070] 步骤S1，分别采集患者身体组织中的水和脂肪同相位和反相位的磁共振图像。
[0071] 步骤S2，利用两幅反相位的磁共振图像计算出附加相位。
[0072] 步骤S3，选取任一幅图像，对图像进行信号判断，识别出信号点。
[0073] 步骤S4，对附加相位进行去卷绕。去卷绕时仅对识别出的信号点的附加相位进行去卷绕。
[0074] 步骤S5，将去卷绕后的两幅反相位的图像结合一副同相位的图像得出水图像和脂肪图像。
[0075] 对于整幅图像，有的像素点是由患者回波产生的信号点，这种信号点包括水的信号和脂肪的信号，由两者合成，而有的像素点是由干扰(例如患者移动、外界干扰或系统不稳)造成的噪声点。如果对每个像素点都进行去卷绕，必然用时较长，导致去卷绕处理时间长。本实施例中，首先在所有的像素点中识别出信号点，去卷绕时仅对识别出的信号点的附加相位进行去卷绕，一方面节约了处理时间，另一方面，减少了噪声点对后续去卷绕的影响，避免了这些噪声点在去卷绕过程的误差传递，所以本实施例还可有效抑制误差的传播，从而提高了水脂分离的成功率。
[0076] 在去卷绕时，可采用现有的去卷绕方法，也可采用本发明上述实施例中的方法，即先采用多项式拟合法对附加相位进行去卷绕，然后在多项式拟合法去卷绕的结果上采用区域增长法进行去卷绕处理。在采用多项式拟合法对附加相位进行去卷绕步骤中仅对识别出的信号点的附加相位进行去卷绕，然后将多项式拟合去卷绕的结果代入区域增长法去卷绕步骤中。在本实施例中，在采用区域增长法去卷绕步骤中也不需要对所有图像点进行增长处理，只需要对从所有的像素点中识别出的信号点进行增长处理。而上述实施例中，在采用区域增长法去卷绕步骤中处理的信号点是所有像素点(即图像点)，所以处理时间较长。
[0077] 本实施例在去卷绕之前先识别出信号点，所以保证识别出的信号点的准确性就变得非常重要。在本发明的具体实例中，可采用幅度判断法来识别信号点，也可采用相位判断法来识别信号点，还可采用两者结合的方法来识别信号点。
[0078] 幅度判断法是将图像中的像素的幅度和最终幅度阈值进行比较，如果像素的幅度大于或等于最终幅度阈值则认为该像素是信号点，如果像素的幅度小于最终幅度阈值则认为该像素不是信号点。
[0079] 相位判断法是将相邻像素的相位进行比较，如果两者的比较结果小于相位阈值，则认为相邻像素为信号点。
[0080] 下面详细说明将幅度判断法和相位判断法结合来进行信号点识别的实施例，其流程如图8所示，包括以下步骤：
[0081] 步骤S31，计算初始幅度阈值步骤。所述初始幅度阈值等于图像中的最大幅度值与最小幅度值的平均值，即初始幅度阈值估计值T＝(max(I)+min(I))/2，I是图像幅度数据，然后执行步骤S32。
[0082] 步骤S32，采用幅度阈值对像素进行分割步骤。将图像中所有像素的幅度值或灰度值和初始幅度阈值进行比较，得到第一组像素集合G1和第二组像素集合G2，第一组像素集合G1为大于或等于初始幅度阈值的像素集合，该集合中是信号点，从而可得到信号点的个数。第二组像素集合G2为小于初始幅度阈值的像素集合，为噪音点。当初次进行像素分割时，G1由所有灰度值大于T的像素组成，而G2由所有灰度值小于T的像素组成。
[0083] 步骤S33，中间幅度阈值计算步骤。分别计算第一组像素集合G1的平均幅度值U1和第二组像素集合G2的平均幅度值U2，将第一组像素集合的平均幅度值U1和第二组像素集合的平均幅度值U2的平均值作为当前中间幅度阈值，即新的中间幅度阈值Tnew＝(U1+U2)/2。
[0084] 步骤S34，将新计算出的中间幅度阈值与其前一个中间幅度阈值进行比较，例如比较两者的差值或比值，判断比较结果是否小于预先设定的参数，若是则该当前中间幅度阈值为最终幅度阈值，结束循环，如果后续进行相位判断法，则继续执行步骤S35。如果比较结果不小于预先设定的参数则转向执行步骤S32，在步骤S32中，将图像中所有像素的幅度值和该新计算出的中间幅度阈值进行比较，采用新计算出的中间幅度阈值进行像素分割，循环步骤S32至S34，直到新计算出的中间幅度阈值与其前一个中间幅度阈值的比较结果小于预先设定的参数。分割后得到的第一组像素集合G1中的像素点的个数就是采用幅度判断法得到的信号点的个数。
[0085] 步骤S35，选取初始种子步骤。选取初始种子子单元272a选择图像中幅度最大的像素点为初始种子。由于信号点的幅度一般比噪音的幅度大，所以先假设幅度最大的像素点为信号点，并以它为初始种子进行区域增长。
[0086] 步骤S36，种子增长步骤。种子增长子单元272b用于将增长的种子的相位与其周边相邻像素点的相位依次进行比较，例如比较两者的差值或比值，如果相位比较结果小于第一相位阈值，则该相邻像素点被增长，并将该相邻像素点作为下一轮增长的种子之一，循环增长直到没有相邻像素点被增长，即没有新的增长种子，循环结束后种子的集合即为信号点的集合。本步骤中初始的第一相位阈值可根据比较的方式和经验进行设定，例如根据磁共振系统的性能设定。
[0087] 步骤S37，信号点个数判断子单元272c将步骤S36中增长结束后得出的信号点个数与采用幅度判断法得出的信号点个数进行比较，例如比较两者的差值或比值，判断比较结果是否符合第一预定条件，若是，则采用相位判断法最后得出的信号点为最终结果，用于后续的附加相位去卷绕。否则执行步骤S38。在本步骤的一种具体实例中，判断步骤S36中增长结束后得出的信号点个数与采用幅度判断法得出的信号点个数的比值是否符合第一预定条件，第一预定条件为经验值。
[0088] 步骤S38，第一相位阈值递增子单元272d用于对第一相位阈值进行重新计算，按照设定步长将第一相位阈值进行递增，设定步长为经验值。递增后的值作为新的第一相位阈值，然后转向执行步骤S35，循环直到得出的信号点个数与幅度判断子单元得出的信号点个数的比较结果符合第一预定条件。
[0089] 单独采用幅度判断法或相位判断法来识别信号点的处理流程可分别参照上述相应的步骤。
[0090] 通过将幅度判断法和相位判断法结合来识别信号点可修正由于幅度信号判断导致的图像边缘信号判断失误问题，得到很好的效果。一方面，由于射频线圈的边缘效应，边缘磁场不均匀，导致边缘点的图像信号幅度小，在采用幅度判断法过程中容易被误判为噪声，而相位判断法是基于相位进行判断，不考虑幅度，因此会将符合相位条件的边缘点判断为信号点。另一方面，在磁共振检测过程中，由于外环境的影响，图像上可能出现亮度较大的噪声，如果按照幅度判断法，噪声将会被判断成信号点。但此类噪声的相位通常比较杂乱，相邻点相位差较大，因此在后续的相位判断法中，大量噪声会被排除掉，不会被识别成信号点。因此，相位判断法正好弥补了幅度判断法的不足。
[0091] 另外，本实施例在相位判断法识别信号点之后还对相位判断法的识别结果进行校正，因在图像编码过程中可能造成相位错误，导致采用相位判断法造成误判，因此用幅度判断法的结果对相位判断法的识别结果进行校正，有效避免了采用相位判断法造成误判。
[0092] 上述实施例中，采用区域增长法时选取的初始种子，除了上述提到的获取方案，初始种子还可采用其他获取方法和约束条件。
[0093] 上述实施例可应用于任何磁场强度的磁共振系统中，尤其对于低场磁共振系统，效果更明显。
[0094] 上述实施例中的相位判断法识别信号点还可应用于其他的图像识别系统，区域增长去卷绕方法还可应用于其他领域的角度去卷绕处理。
[0095] 在另外的实施例中，如果分别采集患者身体组织中水和脂肪的一副反相位磁共振图像和两幅同相位的磁共振图像，水和脂肪相位角度为(0，π，2π)，则公式(1)变为：
[0096]
[0097] 公式(2)变为：
[0098]
[0099] 公式(3)变为：
[0100] ei2φ＝I′2-------------------------(3′)
[0101] 以上内容是结合具体的实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。
