(12)发明专利申请
(10)申请公布号 CN 110473210 A(43)申请公布日 2019.11.19
(21)申请号 201910715594.7(22)申请日 2019.08.05
(71)申请人 三峡大学
地址 443002 湖北省宜昌市大学路8号(72)发明人 夏平 师冬霞 任强 雷帮军
唐庭龙 邹耀斌 (74)专利代理机构 北京兴智翔达知识产权代理
有限公司 11768
代理人 蒋常雪(51)Int.Cl.
G06T 7/11(2017.01)G06T 7/143(2017.01)
权利要求书3页 说明书11页 附图4页
CN 110473210 A(54)发明名称
基于置信度传播的图像分割方法和装置(57)摘要
本发明实施例提供一种基于置信度传播的图像分割方法和装置，方法包括：对待分割图像进行小波变换，得到各尺度的特征向量；针对任一尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件；全局能量包括基于迭代分割结果和特征向量得到的标记场能量和特征场能量；将全局能量最小时对应的迭代分割结果作为该尺度的最终分割结果，将该尺度的最终分割结果作为下一尺度的初始分割结果，进行尺度间迭代分割，直至得到最细尺度的最终分割结果；将最细尺度的最终分割结果作为待分割图像的最终分割结果。本发明实施例提供的方法和装置，实现了全局最优的图像分割，提高了图像分割的准确性。
CN 110473210 A
权 利 要 求 书
1/3页
1.一种基于置信度传播的图像分割方法，其特征在于，包括：对待分割图像进行小波变换，得到各尺度的特征向量；针对任一所述尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件；其中，所述全局能量包括基于所述迭代分割结果和所述特征向量得到的标记场能量和特征场能量；
将所述全局能量最小时对应的所述迭代分割结果作为所述任一尺度的最终分割结果，将所述任一尺度的最终分割结果作为下一尺度的初始分割结果，进行尺度间迭代分割，直至得到最细尺度的最终分割结果；
将所述最细尺度的最终分割结果作为所述待分割图像的最终分割结果。2.根据权利要求1所述的基于置信度传播的图像分割方法，其特征在于，所述针对任一所述尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件，之前还包括：
基于最粗尺度的所述特征向量，对所述最粗尺度进行聚类，得到所述最粗尺度的初始分割结果。
3.根据权利要求1所述的基于置信度传播的图像分割方法，其特征在于，所述针对任一所述尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件，具体包括：
针对任一尺度，基于上一迭代分割结果，通过置信度传播算法，更新当前迭代分割结果；
对所述当前迭代分割结果进行平滑优化；基于所述当前迭代分割结果计算标记场能量；
基于所述当前迭代分割结果和所述任一尺度的特征向量，计算特征场能量；基于所述标记场能量和所述特征场能量，得到所述当前迭代分割结果对应的全局能量；
将所述当前迭代分割结果更新为所述上一迭代分割结果，并重复所述基于上一迭代分割结果，通过置信度传播算法，更新当前迭代分割结果，直至达到预设尺度内迭代条件。
4.根据权利要求3所述的基于置信度传播的图像分割方法，其特征在于，所述针对任一所述尺度，基于上一迭代分割结果，通过置信度传播算法，更新当前迭代分割结果，具体包括：
基于如下公式更新所述当前迭代分割结果：
2
CN 110473210 A
权 利 要 求 书
2/3页
式中，xp、xq∈label＝{1，2，...，L}为相邻节点的标记值，其中，label即标记；n为小波变换的尺度，t为当前迭代次数；相互作用的似然函数，
为可观察节点
与隐藏节点
之间的
为两个相邻的隐藏节点xp、xq之间相互作用的势函数，
为尺度n下第t次迭代时节点p向节点
N(p)\\q为不包括节点q的节点p的四邻域系统，q发送的消息，
为置信度。
5.根据权利要求4所述的基于置信度传播的图像分割方法，其特征在于，所述基于所述当前迭代分割结果计算所述标记场能量，具体包括：
基于如下公式计算所述标记场能量：
式中，为尺度n下第t次迭代时的标记场能量；Vc(x)为与基团c有关的势
函数；β为MLL模型势函数。
6.根据权利要求4所述的基于置信度传播的图像分割方法，其特征在于，所述基于所述当前迭代分割结果和所述任一尺度的特征向量，计算所述特征场能量，具体包括：
基于如下公式计算所述特征场能量：
式中，
xij为位置(i，j)的标记值；
为所述特征场能量，为尺度n下位置(i，j)的特征值，
表示在尺度n下的第t次迭代时位置(i，j)处所
对应的均值和方差。
7.根据权利要求3所述的基于置信度传播的图像分割方法，其特征在于，所述针对任一所述尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件，还包括：
基于所述当前迭代分割结果和所述任一尺度的特征向量，应用高斯模型估计当前区域均值和当前区域方差。
8.一种基于置信度传播的图像分割装置，其特征在于，包括：小波变换模块，用于对待分割图像进行小波变换，得到各尺度的特征向量；尺度内迭代模块，用于针对任一所述尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件；其中，所述全局能量包括基于所述迭代分割结果和所述特征向量得到的标记场能量和特征场能量；
3
CN 110473210 A
权 利 要 求 书
3/3页
尺度间迭代模块，用于将所述全局能量最小时对应的所述迭代分割结果作为所述任一尺度的最终分割结果，将所述任一尺度的最终分割结果作为下一尺度的初始分割结果，进行尺度间迭代分割，直至得到最细尺度的最终分割结果；
分割结果模块，用于将所述最细尺度的最终分割结果作为所述待分割图像的最终分割结果。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至7中任一项所述的基于置信度传播的图像分割方法的步骤。
10.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现如权利要求1至7中任一项所述的基于置信度传播的图像分割方法的步骤。
4
CN 110473210 A
说 明 书
基于置信度传播的图像分割方法和装置
1/11页
技术领域
[0001]本发明涉及图像分割技术领域，尤其涉及一种基于置信度传播的图像分割方法和装置。
背景技术
[0002]图像分割是一种重要的图像处理技术，是解决图像分析与识别的关键步骤。图像分割能够将原始图像中的目标从复杂的背景中提取出来，为后续的图像处理提供可靠依据。
[0003]马尔科夫随机场模型，是统计学方法中最常用的一种图模型，它将概率论与图论联系起来，并且有机地结合图像的时间和空间等信息，利用众多的变量独立关系，构造了一种基于联合概率分布的结构模型，为处理图像分割领域的不确定性问题提供了重要的渠道。
[0004]马尔科夫随机场模型因其可以表达图像像素间的空间信息，受到了广泛的关注。目前，马尔科夫随机场模型已应用于不同的图像分割领域，并存在针对不同的应用环境对马尔科夫随机场模型进行改进的多种方法。[0005]然而，将马尔科夫随机场模型应用在图像分割领域，存在如下问题：在进行能量函数最小化时，易陷入局部最小，导致搜寻效率低下、目标图像分割精确度不高，无法保证整幅图像分割中各目标区域内的一致性及不同目标区域间的差异性。
发明内容
[0006]本发明实施例提供一种基于置信度传播的图像分割方法和装置，用以解决现有的马尔科夫随机场图像分割结果的全局性能差的问题。[0007]第一方面，本发明实施例提供一种基于置信度传播的图像分割方法，包括：[0008]对待分割图像进行小波变换，得到各尺度的特征向量；[0009]针对任一所述尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件；其中，所述全局能量包括基于所述迭代分割结果和所述特征向量得到的标记场能量和特征场能量；
[0010]将所述全局能量最小时对应的所述迭代分割结果作为所述任一尺度的最终分割结果，将所述任一尺度的最终分割结果作为下一尺度的初始分割结果，进行尺度间迭代分割，直至得到最细尺度的最终分割结果；
[0011]将所述最细尺度的最终分割结果作为所述待分割图像的最终分割结果。[0012]第二方面，本发明实施例提供一种无线自组网通信节点，包括：[0013]小波变换模块，用于对待分割图像进行小波变换，得到各尺度的特征向量；[0014]尺度内迭代模块，用于针对任一所述尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件；其中，所述全局能量包括基于所述迭代分割结果和所述特征向量得到的标记场能量和特征场能量；
5
CN 110473210 A[0015]
说 明 书
2/11页
尺度间迭代模块，用于将所述全局能量最小时对应的所述迭代分割结果作为所述
任一尺度的最终分割结果，将所述任一尺度的最终分割结果作为下一尺度的初始分割结果，进行尺度间迭代分割，直至得到最细尺度的最终分割结果；[0016]分割结果模块，用于将所述最细尺度的最终分割结果作为所述待分割图像的最终分割结果。
[0017]第三方面，本发明实施例提供一种电子设备，包括处理器、通信接口、存储器和总线，其中，处理器，通信接口，存储器通过总线完成相互间的通信，处理器可以调用存储器中的逻辑指令，以执行如第一方面所提供的方法的步骤。[0018]第四方面，本发明实施例提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如第一方面所提供的方法的步骤。[0019]本发明实施例提供的一种基于置信度传播的图像分割方法和装置，在尺度内迭代中，将包含标记场能量和特征场能量的全局能量最小时对应的迭代分割结果作为最终分割结果，通过两种随机场的相互约束、共同作用来实现该尺度的分割；在尺度间迭代中，将大尺度的最终分割结果直接投影到相邻小尺度上作为初始分割结果，从而实现全局最优的图像分割，提高图像分割的准确性。附图说明
[0020]为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0021]图1为本发明实施例提供的基于置信度传播的图像分割方法的流程示意图；[0022]图2为本发明实施例提供的二阶邻域系统和原子图的示意图；
[0023]图3为本发明另一实施例提供的基于置信度传播的图像分割方法的流程示意图；[0024]图4为本发明实施例提供的小波变换示意图；
[0025]图5为本发明实施例提供的小波变换的状态空间与特征空间对应结构示意图；[0026]图6为本发明实施例提供的待分割图像；
[0027]图7为本发明实施例提供的待分割图像的最终分割结果；
[0028]图8为本发明实施例提供的基于置信度传播的图像分割装置的结构示意图；[0029]图9为本发明实施例提供的电子设备的结构示意图。
具体实施方式
[0030]为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。[0031]图1为本发明实施例提供的基于置信度传播的图像分割方法的流程示意图，如图1所示，基于置信度传播的图像分割方法，包括如下步骤：[0032]步骤110，对待分割图像进行小波变换，得到各尺度的特征向量。
6
CN 110473210 A[0033]
说 明 书
3/11页
具体地，小波变换(wavelet transform，WT)继承和发展了短时傅立叶变换局部化
的思想，同时又克服了窗口大小不随频率变化等缺点，能够提供一个随频率改变的“时间-频率”窗口，是进行信号时频分析和处理的理想工具。[0034]待分割图像即需要进行图像分割的图像。在预先设置好小波变换的参数，例如小波域分解层数后，对待分割图像进行小波变换，得到各个分解层号分别对应的特征向量，即各个尺度下的特征向量。[0035]步骤120，针对任一尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件；其中，全局能量包括基于迭代分割结果和特征向量得到的标记场能量和特征场能量。[0036]步骤130，将全局能量最小时对应的迭代分割结果作为该尺度的最终分割结果，将该尺度的最终分割结果作为下一尺度的初始分割结果，进行尺度间迭代分割，直至得到最细尺度的最终分割结果。[0037]具体地，在小波变换完成后，需要针对各个尺度分别进行尺度内和尺度间的分割迭代。针对任一尺度，首先获取上一尺度的最终分割结果作为该尺度的初始分割结果，并在初始分割结果的基础上，基于置信度传播算法进行尺度内的迭代分割，得到每次迭代下的迭代分割结果，并计算每次迭代下的迭代分割结果的标记场能量和特征场能量，并由此构成每次迭代分割结果的全局能量。此处，特征场通过描述同标记的小波系数之间的相关性反映每一像素观测特征属性，标记场通过考虑像素邻域标记间的相互约束、相互作用来反映待分割图像分割的区域一致性，特征场和标记场相互约束、共同作用来实现该尺度的分割。
[0038]在尺度内迭代时，应用置信度传播(Belief Propagation，BP)算法作为图像分割的迭代算法。置信度传播算法是一种在图模型上进行推断的消息传递算法，对于马尔科夫随机场中的每一个节点，通过消息传播，把该节点的概率分布状态传递给相邻的节点，从而影响相邻节点的概率分布状态。此外，预设尺度内迭代条件是针对尺度内迭代分割预先设定的迭代终止条件，预设尺度内迭代条件可以是迭代次数达到预设次数，或者迭代分割结果中的参数满足预设条件等，本发明实施例对此不作具体限定。[0039]在完成尺度内迭代分割后，即达到预设尺度内迭代条件后，从每次迭代下的迭代分割结果中，选取全局能量最小值对应的迭代分割结果作为该尺度的最终迭代结果，并作为下一尺度的初始分割结果参与下一尺度的尺度间迭代。[0040]步骤140，将最细尺度的最终分割结果作为待分割图像的最终分割结果。[0041]需要说明的是，尺度间迭代是由最粗尺度向最细尺度迭代的过程.完成尺度间的迭代，即可得到最细尺度的最终分割结果，该结果也是待分割图像的最终分割结果。[0042]本发明实施例提供的方法，在尺度内迭代中，将包含标记场能量和特征场能量的全局能量最小时对应的迭代分割结果作为最终分割结果，通过两种随机场的相互约束、共同作用来实现该尺度的分割；在尺度间迭代中，将大尺度的最终分割结果直接投影到相邻小尺度上作为初始分割结果，从而实现全局最优的图像分割，提高图像分割的准确性。[0043]基于上述实施例，该方法中，步骤120之前还包括：基于最粗尺度的特征向量，对最粗尺度进行聚类，得到最粗尺度的初始分割结果。[0044]具体地，在进行尺度间迭代时，通常将上一尺度的最终分割结果作为当前尺度的
7
CN 110473210 A
说 明 书
4/11页
初始分割结果。而最粗尺度是首个尺度，最粗尺度不存在对应的上一尺度，因此需要对最粗尺度进行聚类，将聚类结果作为初始分割结果。此处，聚类的方法可以是k均值聚类算法K-means，也可以是均值漂移聚类、基于密度的聚类(DBSCAN)等，本发明实施例对此不作具体限定。
[0045]基于上述任一实施例，该方法中，步骤120具体包括：[0046]步骤121，针对任一尺度，基于上一迭代分割结果，通过置信度传播算法，更新当前迭代分割结果。[0047]具体地，基于如下公式更新当前迭代分割结果：
[0048]
[0049]
[0050][0051]
式中，xp、xq∈label＝{1，2，...，L}为相邻节点的标记值，其中，label即标记；n为
为可观察节点
与隐藏节点
之
小波变换的尺度，t为当前迭代次数；间的相互作用的似然函数，
为两个相邻的隐藏节点xp、xq之间相互作用的势
为尺度n下第t次迭代时节点p
函数，N(p)\\q为不包括节点q的节点p的四邻域系统，向节点q发送的消息，
[0052]
为置信度。
即节点的小波系数，隐藏节点
即节点的标记
需要说明的是，可观察节点
值，当前迭代分割结果由当前迭代下所有节点的标记值构成。[0053]步骤122，对当前迭代分割结果进行平滑优化。[0054]具体地，基于置信度传播算法更新当前迭代分割结果后，还可以对当前迭代分割结果进行平滑，消除部分孤立点，以得到更好的分割效果。此处，平滑优化可以是采用3*3模板进行，也可以采用其他规格的模板进行。[0055]步骤123，基于当前迭代分割结果计算标记场能量。[0056]具体地，基于如下公式计算标记场能量：
[0057]
[0058]
8
CN 110473210 A[0059]
说 明 书
5/11页
式中，为尺度n下第t次迭代时的标记场能量；Vc(x)为与基团c有关
的势函数；β为MLL(multilevel logistic model，多级逻辑模型)模型势函数。此处，基团c
选取二阶邻域势团，图2为本发明实施例提供的二阶邻域系统和原子图的示意图，图2中斜线填充的方格为节点xs的邻域节点。[0060]步骤124，基于当前迭代分割结果和该尺度的特征向量，计算特征场能量。[0061]具体地，基于如下公式计算特征场能量：
[0062]
[0063]式中，为特征场能量，为尺度n下位置(i，j)的特征值，
xij为位置(i，j)的标记值；表示在尺度n下的第t次迭代时位置(i，j)处所
对应的均值和方差。[0064]步骤125，基于标记场能量和特征场能量，得到当前迭代分割结果对应的全局能量。
[0065]具体地，基于如下公式获取全局能量E(n)(t)：
[0066][0067]
式中，S为待分割图像的状态空间。[0068]步骤126，将当前迭代分割结果更新为上一迭代分割结果，并重复基于上一迭代分割结果，通过置信度传播算法，更新当前迭代分割结果，直至达到预设尺度内迭代条件。[0069]具体地，在完成当前迭代次数对应的迭代分割后，即可执行下一次迭代，此时当前迭代分割结果即上一次迭代对应的迭代分割结果，即上一迭代分割结果，并重新执行上述步骤121至步骤125，实现当前迭代次数对应的迭代分割，更新当前迭代分割结果，并计算当前迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件，结束尺度内迭代。[0070]基于上述任一实施例，该方法中，步骤120还包括：基于当前迭代分割结果和该尺度的特征向量，应用高斯模型估计当前区域均值和当前区域方差。[0071]具体地，基于如下公式获取当前区域均值和当前区域方差：
[0072]
[0073]
[0074]式中，和分别为当前区域均值和当前区域方差；(i，j)为像素的
9
CN 110473210 A
说 明 书
为尺度n下的像素(i，j)对应的特征值，
6/11页
坐标，S为待分割图像的状态空间，为尺
度下的像素(i，j)对应的标记值，η＝{(-1，-1)，(-1，0)，(-1，1)，(0，1)}为一阶平面区域，p为
[0075]
取值的概率，此概率仅与以(i，j)为中心的邻域η的像素状态有关。
由此得到的当前区域均值和当前区域方差，可以作为衡量当前迭代是否满足预设
尺度内迭代条件的参数条件，判断是否结束尺度内迭代。[0076]基于上述任一实施例，图3为本发明另一实施例提供的基于置信度传播的图像分割方法的流程示意图，如图3所示，该方法包括：[0077]301，设置小波变换层数为J-1，分类数为label＝(1，2，…，L)。[0078]302，设置预设尺度内迭代条件为P。[0079]303，对待分割图像I进行小波变换。[0080]此处，待分割图像I定义在大小为M×N的网格S上，做J-1层小波分解，每一小波尺度用对应的层号n来表示，1≤n≤J-1。[0081]小波域中，分辨率尺度n上构建的状态空间用S(n)描述，则所有分辨率状态空间构成集合S＝{S(0)，…S(n)，…，S(J-1)}；尺度n内状态空间表示为
表示尺度n上像素(i，j)的状态值；同理，构建各分辨
率特征集合w＝{w
(0)
，…，w
(n)
，…，w
(J-1)
}，尺度n内的特征场
[0082]
图4为本发明实施例提供的小波变换示意图，如图4所示，每一尺度由4个频带(LL、
图5为本发
LH、HL、HH)的小波系数组成，本发明实施例中，令
明实施例提供的小波变换的状态空间与特征空间对应结构示意图，如图5所示，状态空间S
和观测特征W空间一一对应。[0083]304，使用K-means聚类算法，对最粗尺度，即尺度J-1进行聚类，得到所述最粗尺度J-1的初始分割结果[0084]305，针对尺度J-1，基于置信度传播算法进行尺度内迭代分割：[0085]首先，使用高斯模型进行参数估计，得到当前区域均值和当前区域方差：
[0086]
[0087]
[0088]式中，和分别为当前区域均值和当前区域方差；(i，j)为像素的
为尺度下的像素(i，j)对应的标记
坐标，为尺度n下的像素(i，j)对应的特征值，
10
CN 110473210 A
说 明 书
7/11页
值，η＝{(-1，-1)，(-1，0)，(-1，1)，(0，1)}为一阶平面区域，p为取值的概率，此概率仅
与以(i，j)为中心的邻域η的像素状态有关。
[0089]基于当前区域均值和当前区域方差判断当前迭代分割结果是否满足预设尺度内迭代条件P，若满足，则将全局能量最小时对应的迭代分割结果作为该尺度的最终分割结果，并执行步骤306；[0090]否则，基于上一迭代分割结果，通过置信度传播算法，更新当前迭代分割结果：
[0091]
[0092]
[0093][0094]
式中，xp、xq∈label＝{1，2，...，L}为相邻节点的标记值，其中，label即标记；n为
为可观察节点
与隐藏节点
之
小波变换的尺度，t为当前迭代次数；间的相互作用的似然函数，
为两个相邻的隐藏节点xp、xq之间相互作用的势
为尺度n下第t次迭代时节点p
函数，N(p)\\q为不包括节点q的节点p的四邻域系统，向节点q发送的消息，
[0095][0096][0097]
为置信度。
此后，对当前迭代分割结果进行平滑优化。
再者，基于当前迭代分割结果计算标记场能量：
[0098][0099]
式中，为尺度n下第t次迭代时的标记场能量；Vc(x)为与基团c有关
的势函数；β为MLL(multilevellogistic model，多级逻辑模型)模型势函数。此处，基团c选
取二阶邻域势团。
[0100]基于当前迭代分割结果和该尺度的特征向量，计算特征场能量：
[0101]
11
CN 110473210 A[0102]
说 明 书
为特征场能量，
8/11页
式中，为尺度n下位置(i，j)的特征值，
xij为位置(i，j)的标记值；表示在尺度n下的第t次迭代时位置(i，j)处所
对应的均值和方差。
[0103]最终，基于标记场能量和特征场能量，得到当前迭代分割结果对应的全局能量E(n)(t)
：
[0104]
完成上述操作中，更新迭代次数，重新执行步骤305。
[0106]306，将下一尺度作为当前尺度，令n＝n-1，并判断n是否小于0，若n小于0，则结束上述流程，将尺度n+1下的最终分割结果作为图像的最终分割结果；否则，执行307。[0107]307，将尺度n+1下的最终分割结果作为尺度n的初始分割结果，使用高斯模型进行参数估计，得到当前区域均值和当前区域方差。
[0108]基于当前区域均值和当前区域方差判断当前迭代分割结果是否满足预设尺度内迭代条件P，若满足，则将全局能量最小时对应的迭代分割结果作为该尺度的最终分割结果，并执行步骤306；[0109]否则，基于上一迭代分割结果，通过置信度传播算法，更新当前迭代分割结果。[0110]此后，对当前迭代分割结果进行平滑优化。[0111]再者，基于当前迭代分割结果计算标记场能量，基于当前迭代分割结果和该尺度的特征向量，计算特征场能量。最终，基于标记场能量和特征场能量，得到当前迭代分割结果对应的全局能量。
[0112]完成上述操作中，更新迭代次数，重新执行步骤307。[0113]基于上述任一实施例，图6为本发明实施例提供的待分割图像，如图6所示，待分割图像为大小为512×512的含有结节的肺部CT图像。图7为本发明实施例提供的待分割图像的最终分割结果。对比图6和图7可知，上述实施例中提供的基于置信度传播的图像分割方法，可以有效提升分割结果的质量，在最终分割结果中不但可以清晰地分割出感兴趣区域，而且很好地保留了区域的细节部分。[0114]基于上述任一实施例，图8为本发明实施例提供的基于置信度传播的图像分割装置的结构示意图，如图8所示，该装置包括小波变换模块810、尺度内迭代模块820、尺度间迭代模块830和分割结果模块840；
[0115]小波变换模块810用于对待分割图像进行小波变换，得到各尺度的特征向量；[0116]尺度内迭代模块820用于针对任一所述尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件；其中，所述全局能量包括基于所述迭代分割结果和所述特征向量得到的标记场能量和特征场能量；
[0117]尺度间迭代模块830用于将所述全局能量最小时对应的所述迭代分割结果作为所述任一尺度的最终分割结果，将所述任一尺度的最终分割结果作为下一尺度的初始分割结果，进行尺度间迭代分割，直至得到最细尺度的最终分割结果；
12
[0105]
CN 110473210 A[0118]
说 明 书
9/11页
分割结果模块840用于将所述最细尺度的最终分割结果作为所述待分割图像的最
终分割结果。
[0119]本发明实施例提供的装置，在尺度内迭代中，将包含标记场能量和特征场能量的全局能量最小时对应的迭代分割结果作为最终分割结果，通过两种随机场的相互约束、共同作用来实现该尺度的分割；在尺度间迭代中，将大尺度的最终分割结果直接投影到相邻小尺度上作为初始分割结果，从而实现全局最优的图像分割，提高图像分割的准确性。[0120]基于上述任一实施例，该装置还包括：[0121]聚类模块，用于基于最粗尺度的所述特征向量，对所述最粗尺度进行聚类，得到所述最粗尺度的初始分割结果。[0122]基于上述任一实施例，该装置中，尺度内迭代模块820包括：[0123]置信度传播更新单元，用于针对任一尺度，基于上一迭代分割结果，通过置信度传播算法，更新当前迭代分割结果；[0124]平滑单元，用于对所述当前迭代分割结果进行平滑优化；[0125]标记场能量计算单元，用于基于所述当前迭代分割结果计算标记场能量；[0126]特征场能量计算单元，用于基于所述当前迭代分割结果和所述任一尺度的特征向量，计算特征场能量；
[0127]全局能量计算单元，用于基于所述标记场能量和所述特征场能量，得到所述当前迭代分割结果对应的全局能量；[0128]迭代判断单元，用于将所述当前迭代分割结果更新为所述上一迭代分割结果，并重复所述基于上一迭代分割结果，通过置信度传播算法，更新当前迭代分割结果，直至达到预设尺度内迭代条件。
[0129]基于上述任一实施例，置信度传播更新单元，具体用于：[0130]基于如下公式更新所述当前迭代分割结果：
[0131]
[0132]
[0133][0134]
式中，xp、xq∈label＝{1，2，...，L}为相邻节点的标记值，其中，label即标记；n为
为可观察节点
与隐藏节点
之
小波变换的尺度，t为当前迭代次数；间的相互作用的似然函数，
为两个相邻的隐藏节点xp、xq之间相互作用的势
为尺度n下第t次迭代时节点p
函数，N(p)\\q为不包括节点q的节点p的四邻域系统，
13
CN 110473210 A
说 明 书
为置信度。
10/11页
向节点q发送的消息，
[0135][0136][0137]
基于上述任一实施例，该装置中，标记场能量计算单元具体用于：基于如下公式计算所述标记场能量：
[0138][0139]
式中，为尺度n下第t次迭代时的标记场能量；Vc(x)为与基团c有关
的势函数；β为MLL模型势函数。
[0140]基于上述任一实施例，该装置中，特征场能量计算单元具体用于：[0141]基于如下公式计算所述特征场能量：
[0142]
[0143]式中，为所述特征场能量，为尺度n下位置(i，j)的特征
值，xij为位置(i，j)的标记值；表示在尺度n下的第t次迭代时位置(i，j)处
所对应的均值和方差。
[0144]基于上述任一实施例，该装置中，尺度内迭代模块820还包括：[0145]参数评估单元，用于基于所述当前迭代分割结果和所述任一尺度的特征向量，应用高斯模型估计当前区域均值和当前区域方差。
[0146]图9为本发明实施例提供的电子设备的实体结构示意图，如图9所示，该电子设备可以包括：处理器(processor)901、通信接口(Communications Interface)902、存储器(memory)903和通信总线904，其中，处理器901，通信接口902，存储器903通过通信总线904完成相互间的通信。处理器901可以调用存储在存储器903上并可在处理器901上运行的计算机程序，以执行上述各实施例提供的基于置信度传播的图像分割方法，例如包括：对待分割图像进行小波变换，得到各尺度的特征向量；针对任一所述尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件；其中，所述全局能量包括基于所述迭代分割结果和所述特征向量得到的标记场能量和特征场能量；将所述全局能量最小时对应的所述迭代分割结果作为所述任一尺度的最终分割结果，将所述任一尺度的最终分割结果作为下一尺度的初始分割结果，进行尺度间迭代分割，直至得到最细尺度的最终分割结果；将所述最细尺度的最终分割结果作为所述待分割图像的最终分割结果。[0147]此外，上述的存储器903中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分
14
CN 110473210 A
说 明 书
11/11页
可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-0nly Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0148]本发明实施例还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各实施例提供的基于置信度传播的图像分割方法，例如包括：对待分割图像进行小波变换，得到各尺度的特征向量；针对任一所述尺度，基于置信度传播算法进行尺度内迭代分割，并计算每一迭代分割结果对应的全局能量，直至达到预设尺度内迭代条件；其中，所述全局能量包括基于所述迭代分割结果和所述特征向量得到的标记场能量和特征场能量；将所述全局能量最小时对应的所述迭代分割结果作为所述任一尺度的最终分割结果，将所述任一尺度的最终分割结果作为下一尺度的初始分割结果，进行尺度间迭代分割，直至得到最细尺度的最终分割结果；将所述最细尺度的最终分割结果作为所述待分割图像的最终分割结果。
[0149]以上所描述的通信节点实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。[0150]通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。[0151]最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
15
CN 110473210 A
说 明 书 附 图
1/4页
图1
图2
16
CN 110473210 A
说 明 书 附 图
2/4页
图3
图4
17
CN 110473210 A
说 明 书 附 图
3/4页
图5
图6
18
CN 110473210 A
说 明 书 附 图
4/4页
图7
图8
图9
19
因篇幅问题不能全部显示，请点此查看更多更全内容