基于边缘约束局部区域MRF 的图像分割方法

胡高珍，徐胜军，孟月波，刘光辉，冯 峰，段中兴

（西安建筑科技大学 信息与控制工程学院，西安 710055）

0 概述

图像分割是图像处理领域的重要研究方向。常用的图像分割方法主要有基于阈值的分割方法、基于聚类的图像分割方法和基于模型的分割方法。近年来，基于模型的分割方法［1］成为了研究热点，其中高斯混合模型（Gaussian Mixture Model，GMM）因其建模简单有效而受到广泛关注。而基于GMM 图像分割方法通常假设图像中的像素点相互独立，没有考虑图像像素间的空间约束关系，导致其分割准确性欠佳。为解决该问题，研究人员通过对马尔科夫随机场（Markov Random Field，MRF）模型引入图像的空间约束关系，建立的高斯马尔科夫随机场（Gaussian MRF，GMRF）可有效提升图像分割效果，并且在图像分割、图像恢复以及图像匹配等领域得到广泛应用。然而，MRF 模型经常使用Potts 类型先验知识建立图像的空间约束关系，这种简单的像素点对关系不能准确描述复杂图像的先验知识，且在分割过程中易产生过平滑现象。

为提高MRF 模型对复杂图像先验知识的描述能力，需引入更多邻域节点来表达图像的空间一致性约束。文献［2］针对传统GMM 难以自动获取类属数和对噪声敏感的问题，提出一种基于可变类空间约束GMM 的遥感图像分割方法。利用MRF 模型将GMM 的像素类属性建模为图像的先验概率，以提高GMM 的抗噪性。文献［3］提出一种局部空间自适应MRF 模型，通过在带权值的自适应空间GMM 中引入像素的局部空间信息，使得图像的噪声或者纹理突变具有一定的鲁棒性，但由于区域模型的平滑作用，将导致分割结果出现边缘带现象。文献［4］提出一种快速鲁棒的改进GMM（Modified GMM，MGMM）图像分割模型，该模型利用欧氏距离度量局部窗内像素间的相似性并引入图像局部空间信息。为保持图像的边缘结构，该模型设定了较小的MRF 模型先验参数β，但是对于自然图像而言，其全局非同态性使得具有较小先验惩罚的MRF 模型对复杂图像强度信息的突变较为敏感，不能得到更平滑的分割结果。文献［5］提出基于图像片权重MRF 分割模型，对邻域内不同的图像片根据相似度计算对应权重，同时通过K-L（Kullback-Leibler）距离引入熵的惩罚，从而解决过分割问题。文献［6］提出基于改进高斯混合模型的运动图像目标检测算法，将邻域像素值加权和作为像素重塑后的值，并利用像素空间距离建立滤波窗口的各像素权值，以有效抑制噪声并且更好地保留图像细节。文献［7］提出空间约束非对称GMM，采用K-L 距离计算邻域像素先验分布，引入了图像空间信息，并根据像素类内和类间先验知识和后验概率建立图像似然特征的非对称分布模型，有效拟制噪声干扰。

上述模型利用局部空间信息提高了分割模型对噪声或纹理突变干扰的鲁棒性，但由于这些模型没有考虑到图像的边缘信息，因此局部区域能量的最小化易导致分割结果过平滑，使得图像分割结果常出现边缘模糊、边界不清的问题。为解决该问题，研究人员将图像边缘信息融合到MRF 模型中，以提升图像分割的边缘效果。文献［8］基于隐马尔科夫随机场（Hidden MRF，HMRF）理论提出了HMRF-EMedge 模型，利用邻域像素相互作用引入空间信息，并使用图像边缘特征来保持分割边缘，但对于图像中边缘细小的区域，分割结果依然不理想。文献［9］提出一种结合形状信息的静态MST 区域划分和RHMRF-FCM 算法的图像分割方法，并定义一种MST 静态同质区域划分准则，通过借助MST 引入边界信息，可较好地表达图像不规则边界和形状信息，有效抑制几何噪声。文献［10］提出一种结合MRF能量和活动轮廓模型的方法，利用Garbo 纹理特征、DCE-MRI 时域特征和灰度特征构成的特征向量与聚类中心向量的距离构建一种模糊速度函数，并将其引入到活动轮廓模型中作为MRF 能量模型的边缘检测项。文献［11］针对多尺度马尔科夫随机场（MRMRF）模型非重叠区域在最优化过程中造成的块效应，提出一种具有边缘保持的MRMRF 模型图像分割方法，有效抑制MRF 能量模型在分割过程中向图像的边缘收敛，但这种Cauchy 分布对于图像弱边缘的惩罚较大，不利于保持图像的细节边缘特征。

为有效表达复杂图像局部统计特征并保持图像分割区域的边缘结构，本文提出一种基于边缘约束局部区域MRF（Edge Constrained Local Region MRF，ECLRMRF）的分割模型。该模型先利用图像局部区域内像素的欧氏距离度量像素的相似度，构建局部空间约束高斯混合模型（Local Spatial Constrained Gaussian Mixture Model，LSCGMM），以提取图像的局部区域似然特征。为约束图像分割区域的边缘，利用Canny 边缘检测算子提取图像边缘特征和图像分割区域的边缘结构，并在MRF 模型框架下将图像的局部区域似然特征和边缘特征相融合，建立具有边缘约束的局部区域MRF 分割能量模型，再采用Gibbs 算法优化MRF 模型并得到最终分割结果。

1 基于GMRF 的图像分割模型

假设S={1,2,…,N}表示一个格点集合，其中N表示集合S中的格点总数。令X={xs|s∊S}表示观测图像，其中xs是指图像中第s个像素的强度值。令Y={ys|s∊S}表示分割图像的标签场，ys∊{1,2,…,K}表示第s个像素的标签，其中K为图像的类别总数，图像中每一个像素xs对应唯一的标签ys。

在基于模型的图像分割方法中，GMM 是描述图像像素分配标签过程的理想模型之一。假定分割过程是一个高斯过程，即图像像素的强度值与其对应标签的概率分布可用GMM 概率密度函数描述为：

其中，π(ys=k)表示给一个像素分配为第k类图像标签的先验概率，且满足约束条件0 ≤π(ys=k) ≤1 和（ys=k）=1；k为分割类别且k∊{1,2,…,K}；Θ={θ1,θ2,…,θK}是GMM参数。p(xs|θk)表示GMM第k类高斯分布的概率密度函数，计算公式如下：

其中，θk=(μk,σk)分别为第k类高斯分布的均值和方差。

传统GMM 通常假定图像像素是独立同分布的，即不考虑邻域像素点的空间约束关系时，由式（1）可得基于GMM 的对数似然图像分割模型为：

为了在基于GMM 图像分割模型中引入邻域像素的相关性，常利用MRF 建模图像先验概率π(ys=k)刻画邻域像素的空间约束关系。假设N(s)={Ns|∀s∊S}为定义在图像域X上的一个邻域系统，Ns为第s个节点的邻域像素集合，且s∉Ns。根据Hammersley-Clifford 定理［15］，基于MRF 描述的图像空间Gibbs 先验分布可表示为：

其中，Z(β)为归一化常数，T为温度参数，β为模型先验参数且是一个Gibbs 先验能量项。根据贝叶斯定理并结合式（1）和式（4），可得到GMRF 后验概率密度联合分布为：

其中，P(Y|X,β,Θ)为标签场Y在给定观察条件X,β,Θ下的后验概率分布，P(X|Y,Θ)为观察图像场X关于标签场Y的似然概率分布，Θ={θ1,θ2,…,θK}为GMM 参数，P(Y|β)为标签场Y的先验概率分布，β为先验参数，k为类别数。

GMRF 模型同时考虑了图像强度特征和空间约束先验知识，能够更好地对图像进行分割。但是常规Pairwise MRF 模型中的先验参数β使空间邻域内像素点趋于取相同的类标签，对于含有丰富特征的图像易产生细节特征的过平滑问题，且不能很好地描述边缘，这是由邻域内每个像素对中心像素的惩罚相同的假设条件不合理造成的［16］。

2 局部区域一致性MRF 的分割模型

为有效保持图像分割的边缘，并提高基于GMRF的图像分割方法在分割过程中对噪声、纹理等突变信号的鲁棒性，本文提出一种局部区域一致性MRF 分割模型。在GMRF 模型中，图像的标签过程缺乏像素间的空间信息，易受噪声或图像非同态性的干扰。为解决上述问题，本文提出模型在图像观察场X中利用K-L距离引入邻域内像素空间关系，建立局部空间约束高斯混合模型（Local Spatial Constrained GMM，LSCGMM）。在标签场Y中，针对常规MRF 的点对交互结构不能有效表达复杂图像的局部先验知识问题，本文提出一种基于区域的部分二阶Potts 先验模型，并对局部区域内标签取不同值时进行不一致性惩罚。在MRF 框架下，将图像的局部空间约束高斯似然特征和局部区域一致性先验知识相结合来建立图像的局部区域一致性MRF分割模型。

2.1 局部空间约束高斯混合模型

常规GMM 模型的独立性假设忽略了图像的空间先验关系，为解决此问题，建立图像局部区域内像素与标签带空间约束的似然关系并提出一种LSCGMM 模型。令xs为图像局部区域内中心像素点，定义其邻域Ns内的像素点为xr，即xr={xr|r∊Ns} 。在局部区域内，倾向于取相同标签的相邻像素强度值相近，反之则相邻像素强度值相差较大。基于此假设，在距离度量方式中引入局部区域内像素的空间约束关系，获得中心像素xs与邻域像素xr的空间相似性权重，并将该权重引入到GMM 中，得到LSCGMM 似然模型。

其中，wr=为一个归一化函数，对于区域内所有像素，当w(yr)=1 时，wr为区域内像素的个数。包含空间信息的权重函数w(yr)可表示为：

其中，δ2是一个高斯核函数的方差。显然，在邻域系统Ns中，中心像素点xs与其邻域像素点xr之间距离越大，则相似度越低、赋予像素权重越小。通过引入邻域系统像素间相似性，在分割过程中引入像素的空间相关性信息，能够将像素类的分布根据邻域空间信息自适应调整，约束图像区域内具有相似像素强度值的标签一致性。从式（6）、式（7）可知，LSCGMM 模型利用像素间的距离描述图像的局部空间关系，通过空间相似性权重约束邻域内像素的标签一致性，使得同质区域像素具有相同类标签，避免了噪声对分割结果的干扰。

2.2 部分二阶区域Potts 先验模型

常规Pairwise MRF 模型的简单点对结构限制了MRF 模型的表达能力，使其无法表达复杂自然图像的区域和全局先验知识［17］，而高阶MRF（Higherorder MRF，HMRF）虽然具有较强的表达能力，能够表达更加复杂的图像先验知识和统计信息，但是由于高阶势能项学习和推理较高的复杂度限制了HMRF 模型在图像处理领域中的进一步应用。为解决此问题，本文提出一种部分二阶区域Potts 先验模型。该模型首先对图像进行区域划分，建立大小为w×w的局部区域集合W，假定第s个局部区域ws的标签场是一个MRF 分布，局部区域内标签场的先验使用部分二阶区域Potts 先验模型描述，根据MRF 与Gibbs 分布的等价性，则建立的部分二阶区域Potts 先验模型可描述为：

其中，Z(β)是归一化常量，β为Potts 模型先验参数，控制局部区域分割的平滑度。Gd为局部区域内中心节点在水平、垂直、对角和反对角等8 个方向上的惩罚因子，对于水平和垂直方向的邻域节点，令Gd=1，对于对角和反对角方向的邻域节点，由于其距离中心节点较远，因此令为delta 函数，可表示为：

由式（8）可知，本文所提先验模型在局部区域ws内中心节点标签ys与邻域节点标签yr不相等时，这个不一 致惩罚 为。由此可知，该不一致性惩罚值随局部区域内不一致标签数量的增加而呈线性增长，即区域内较多像素点的类标签趋于一致时，具有较小的惩罚值，且避免了常规的区域同一性惩罚造成的过平滑现象。仅当中心节点标签ys与邻域节点标签yr全部相等时，对局部区域内标签的不一致性惩罚为0。与传统基于四邻域结构的Pairwise MRF 模型相比，部分二阶区域结构通过对八邻域结构引入更多的图像边缘线过程特征，有助于描述更为复杂的图像边缘特征。虽然本文所提先验模型引入的较大局部区域包含了图像更多的局部先验知识，但是并未显著增加本文所提先验模型的计算复杂度。

2.3 基于局部区域一致性约束MRF 分割模型

在常规的MRF 分割框架下，根据贝叶斯定理，基于LSCGMM 似然特征和部分二阶区域Potts 先验构造图像的MRF 分割能量模型为：

其中，令EData(X,Y)=-lnP(X|Y,Θ)表示图像观察场与标签场的似然特征能量项，根据提出的LSCGMM 似然模型式（6），则EData(X,Y)可定义为：

令ESmooth(Y,β)=-lnP(Y|β)表示图像标签场的空间先验能量约束项，根据所提部分二阶区域Potts 先验模型式（8），则ESmooth(Y,β)可定义为：

根据式（12）可知，在局部区域ws内，EData(X,Y)利用局部像素空间相关性信息，可有效约束更相似的邻域像素取相同的标签值，对于图像的噪声或者纹理突变信息的干扰具有鲁棒性。ESmooth(Y,β)对于局部区域内标签的不一致惩罚的大小随着邻接标签取值不一致个数的增大而呈线性增长趋势，这种类别不一致性惩罚避免了局部区域内异质像素相同的惩罚值造成的过平滑现象，且能够有效保证图像局部区域一致性约束。因此，本文提出的局部区域一致性MRF 分割模型不仅能对噪声图像或纹理信号的干扰具有鲁棒性，而且可有效提升图像分割结果的边缘效果，从而得到更平滑的分割结果。

3 基于边缘约束局部区域MRF 的分割模型

虽然本文建立的局部区域一致性约束MRF 模型对噪声或纹理突变信息干扰具有较强的鲁棒性，但是由于图像分割对象的边缘呈高度尖峰态分布，而局部区域一致性模型常约束区域内像素倾向于取相同的标签，因此本文所提模型难以描述图像强烈突变的高维复杂性边缘特征，易导致边缘区域的过平滑，使得对含有丰富细节特征的图像产生较大的误分割。图像边缘是图像的基本特征之一，通常携带有图像的大部分信息，这种边缘信息能够有效约束图像边缘分割的准确性。为解决局部区域一致性约束MRF 模型对图像边缘的过平滑问题，本文提出利用图像边缘特征建立分割对象的一致性边缘约束，并融合局部区域一致性模型表示的区域特征，建立一种ECLRMRF 分割模型。

3.1 边缘能量约束先验模型

为建立ECLRMRF 模型，基于Canny 边缘检测算子提取图像的边缘特征，提出一种边缘能量约束先验模型，并融合局部区域一致性MRF 图像分割模型，建立所提ECLRMRF 模型的边缘能量项，且在分割过程中对图像目标分割的边界进行一致性边缘约束。本文所提边缘约束先验模型使用Canny 边缘检测算子提取图像的边缘特征，得到图像的边缘信息B={b1,b2,…,bN}，令bs表示边缘图像中第s个节点，则图像在该点的边缘特征bs可表示为：

其中，bs=1 表明第s个节点是边缘节点，bs=0 表明第s个节点不是边缘节点。

本文定义边缘先验约束能量项为：

其中，γ表示边缘先验参数，其大小衡量边缘约束的重要性。Ns表示第s个节点的邻域节点集合，二阶边缘势函数φsr(xs,xr)可表示为：

式（15）表明，对于邻接节点对(s,r)，如果2 个邻域节点均是边缘上的点或者不是边缘上的点，则边缘约束项的惩罚为0，以保证分割结果不过于平滑，有利于保持图像的边缘约束，反之则给该点增加一个边缘惩罚γ。因此，边缘特征的引入增加了对图像目标边界的有效约束，在能量模型优化过程中，通过边缘能量项的最小化约束图像分割结果向边缘收敛，能够有效提高图像分割边缘的准确度。

3.2 边缘约束局部区域MRF 模型

基于MRF 框架，在能量函数中融合LSCGMM模型表达的似然特征、部分二阶区域Potts 模型建立的区域先验知识和边缘能量先验约束模型描述的图像边缘特征，建立基于边缘约束的局部区域MRF（ECLRMRF）模型可表示为：

本文提出的模型由以下3 个部分组成：局部空间约束高斯混合似然能量项EData(X,Y)、部分二阶区域Potts 先验项ESmooth(Y)和边缘能量约束项EEdge(B)。似然能量项EData(X,Y)基于相似性测量的方式引入了融合邻域空间信息的高斯混合似然特征，保持区域内像素与分配标签的似然一致性。先验项ESmooth(Y)根据局部区域内和中心标签值不一致数量呈正比的线性惩罚建立局部区域先验的一致性平滑约束，并且引入更多的先验空间约束，有效提高分割过程的鲁棒性。边缘能量约束项EEdge(B)通过引入图像边缘特征，增加对目标分割边界的约束，提高图像分割边缘的准确度。

最终根据MAP 准则，本文所提ECLRMRF 模型的最优估计标签Y*为：

3.3 基于ECLRMRF 模型的Gibbs 分割算法

对本文所提ECLRMRF 模型求解最优估计标签Y*时，首先需要估计ECLRMRF 模型的参数(Θ,β,γ)。部分二阶区域Potts 先验模型参数β通常采用人工试错法进行估计，选择分割结果最好的值作为分割图像的最优参数值，选择的原则是对噪声强或者纹理较复杂的图像选择较大的β值，反之则选择较小的β值。边缘先验约束参数γ由人工给定。局部空间约束GMM 参 数Θ*={μ1,μ2,…,μK；σ1,σ2,…,σK}，μk,σk分别是基于LSCGMM 的似然模型的第k类分布的均值和方差。由于LSCGMM 是一个含有未知量Y的参数模型，对于该模型，期望最大化（Expectation Maximization，EM）算法是一种有效的估计算法。EM 算法由期望步（E 步）和最大化步（M 步）组成。E 步根据在前一次迭代估计的参数Θ(t-1)计算隐含变量Y的期望，M 步是对似然函数用最大似然法求解新参数Θ(t)的估计值。EM 算法求解LSCGMM 模型参数的迭代公式为：

本文提出基于ECLRMRF 分割模型的Gibbs 采样算法步骤为：

步骤1给定输入图像X，对待分割图像进行参数初始化：人工给定分类数K，部分二阶区域Potts 先验模型参数β和边缘先验参数γ，Gibbs 采样算法的温度参数初值T0=4℃，局部区域大小给定为3×3。

步骤2LSCGMM 的参数初始化由Kmeans 算法估计。

步骤3采用EM 算法根据式（18）和式（19）估计具有LSCGMM 参 数Θ*={μK,σK}K=1,2,…，k，并根据式（8）计算权重w(yn)。

步骤4由步骤1 给定的参数β、步骤3 估计的LSCGMM 参 数Θ*={μK,σK}K=1,2,…，k以 及权重w(yn)，按照最大似然法根据式（20）对图像进行初始化分割。

步骤5根据ECLRMRF 模型式（16），采用Gibbs 采样算法进行迭代估计给像素xs分配一个标签ys时的接受率，Gibbs 采样算法接受这个分配的标签的概率如式（21）所示：

步骤6算法迭代运行，如果|E(t)(Y|X,B,Θ,β,γ)-E(t-1)(Y|X,B,Θ,β,γ)|≤ε（ε是一个很小的值）或达到最大迭代次数，则算法终止，并输出最优分割结果Y*={y1,y2,…,yN}，否则减小Gibbs 采样算法温度T=0.95T(t-1)，并返回步骤5 继续迭代。

4 实验结果与分析

为验证本文所提模型的可行性和有效性，以人工合成图像和自然图像作为实验对象，分别采用传统高斯马尔科夫随机场（CGMRF）模型［20］、HMRFEM 算法［8］、局部区域高斯马尔科夫随机场（LRGMRF）模型［21］、RLSI-HMRF 模型［22］以及本文所提ECLRMRF 模型进行分割实验。在实验中估计GMM 参数时，由于EM 算法对随机初值敏感，因而上述模型中的GMM 初值都采用K-means 算法进行初始化，不同图像Potts 先验模型参数β和边缘先验约束参数γ由试错法给定。根据大量实验验证，参数β取值区间（0，6］、参数γ取值区间为［0，1.0］可得到较好的分割结果，实验设定γ=0.5。对于噪声弱、较平坦的图像，平滑参数β和边缘先验约束参数γ可取较小值，而对于噪声强、纹理丰富的图像，可增加β和γ的取值，分类数K根据图像由人工指定，且邻域窗口大小设置为3×3。

4.1 人工合成图像分割实验

人工合成图像如图1（a）所示，图像大小为128 像素×128 像素、灰度级为256。图像分割为4 类，各类均值分别为0、85、170 和255。为评估本文所提模型对噪声图像的鲁棒性，人工叠加均值为零、方差为500 的高斯噪声，实验中Potts 先验模型参数β设为1.2，对比模型GMM 参数分别用EM 算法进行估计。

5 种不同模型对人工合成图像的分割图像对比如图1 所示。从图1 可以看出：CGMRF 模型对噪声比较敏感且存在较多的误分点，使得分割效果最差；LRGMRF 模型由于利用了较多的图像空间信息，因而对噪声具有一定的鲁棒性，但仍然存在一些误分割斑点且较大的区域平滑导致边缘的误分割；HMRF-EM 模型受噪声影响出现了误分割现象，且分割结果出现了明显的边缘带；RLSI-HMRF 模型由于建立了远距离的空间关系，分割结果更为平滑但缺乏边缘约束，导致分割结果出现边缘带；与上述4 种模型相比，在本文所提ECLRMRF 模型中，局部区域一致性约束对图像噪声干扰具有更强的鲁棒性，同时边缘特征约束图像分割结果向边缘收敛，有效提高了图像分割边缘的准确度。

图1 5 种模型对合成图像的分割图像对比Fig.1 Comparison of five models for synthetic image segmentation

为更加准确地表明本文所提模型的有效性，对人工合成噪声图像的5 种模型分割结果进行量化评估，并采用分类正确率（Correct Classification Ratio，CCR）［23］作为量化评估指标，且CCR 可定义为：

其中，gtk表示第k个分割的真实结果，segk表示第k个分割的模型分割结果，其数值越接近于1，表示分割结果越好。

表1 给出了5 种模型在合成图像分割过程中的运行时间、迭代次数以及CCR。从表1 可以看出，在相同的噪声水平下，本文所提ECLRMRF 模型的迭代次数低于CGMRF 模型、LRGMRF 模型和RLSIHMRF 模型，且具有较高的CCR 值。这说明在5 种对比模型中，本文所提ECLRMRF 模型的分割结果准确性最优，并且与CGMRF、LRGMRF、HMRFEM、RLSI-HMRF 模型相比，其分割正确率分别提高了约5.23、1.32、41.92 和0.45 个百分点。由于每次迭代计算的区域信息量增多，因此和HMRF-EM 模型、CGMRF 模型、RLSI-HMRF 模型相比，本文模型的迭代时间有所增加，但是迭代时间仍小于LRGMRF 模型。因此，在分割准确度和算法运行时间方面而言，本文所提模型具有较高的综合性能。

表1 5 种模型对合成图像的分割结果对比Table 1 Comparison of segmentation results of five models for synthetic images

4.2 自然图像分割实验

为进一步验证本文所提ECLRMRF 模型的有效性，实验对5 幅自然图像进行分割，而实验图像取自Berkeley 数据集BSDS500［24］。如图2（a）所示，从左到右图像ID 分别为80099、135069、253036、28083、55067，分割类别k分别设定为2、2、3、4、6，实验中Potts 先验模型参数β分别设定为2.1、1.4、2.3、3.2、5.3。从图2 可以看出，CGMRF 模型的分割结果较差，这是由于Pairwise MRF 模型简单的点对交互结构不能充分描述图像的空间先验特征，因此易受噪声或者纹理突变信号的干扰而出现误分割。例如，第1 幅图中的“水面”有较多的分割斑点噪声，第3 幅图中的“草地”“树”和 第4幅图中的“山”“草地”“树”这样纹理丰富的区域都出现了较为明显的误分割。与CGMRF 模型相比，LRGMRF 模型可利用更大的局部统计信息包含更多的空间先验特征，因而受噪声或者纹理突变信号的干扰较小，在噪声较弱或纹理平坦区域有较好的分割结果，但是由于该模型采用的局部区域模型的平滑作用，造成图像的细节特征缺失，且在图像边缘处易出现边缘带现象。第1 幅图中“水牛”与“水面”之间出现的边缘细节缺失，且“水面”区域出现分割斑点。采用HMRF-EM 模型的分割结果在图像纹理丰富区域和特征相近区域效果较差，易出现误分割现象，如第2 幅图中“尾巴”，第3 幅图中的“天空”和第4 幅图中“草地”和“树”的误分割。RLSI-HMRF 模型利用鲁棒的Pn Potts 模型表达了更复杂的图像局部区域特征，因此图像分割结果背景平滑且轮廓清晰，但在纹理复杂的区域以及相似区域仍然会出现少量斑点以及误分割情况，比如第1 幅图中“水面”区域，第2 幅图中“尾巴”区域，第3 幅图中“草地”和“树”以及第4 幅图中“山”和“树”都出现误分割现象。本文所提ECLRMRF 模型与其他4 种模型相比，不仅对噪声或者纹理突变信息具有较强的鲁棒性，且在纹理细节丰富的区域取得了更为平滑的分割结果，如第1 幅图中“水面”、第2 幅图的“翅膀”和“尾巴”、第3 幅图中的“草地”和“树”以及第4 幅图中的“山”均取得了良好的分割结果，而且同时考虑了图像的边缘特征约束，可较好地保持图像的边缘信息，如第5 幅图中“山”之间更为清晰的分割边界。

图2 5 种模型对自然图像的分割效果Fig.2 Segmentation effect of five models on natural images

对于自然图像分割结果常用概率随机索引（Probabilistic Rand Index，PRI）指标衡量［25］，且该指标可定义为：

其中，Sseg表示图像分割结果，Sgt表示标准分割结果。csr=1 表示像素s和r属于同一类别，csr=0 表示像素s和r属于不同类别。PRI∊[0,1]，且其值越大，说明分割准确度越高。

表2 给出了5 幅图像在不同模型下分割结果的PRI 值。从表2 可以看出，本文所提ECLRMRF 模型的PRI 值均高于其他模型，这说明本文模型的分割性能优于其他对比模型。

表2 5 种模型对自然图像分割结果的PRI 值Table 2 PRI values of five models for natural images segmentation

为进一步评价本文所提分割算法的性能，表3给出了实验图像在5 种不同模型下分割过程迭代次数和运行时间的对比。从表3 可以看出，与其他4 种模型相比，本文所提ECLRMRF 模型具有相对较少的迭代次数和运行时间。

表3 5 种模型的迭代次数和运行时间对比结果Table 3 Comparison results of iteration times and running time of five models

5 结束语

本文提出一种基于边缘约束局部区域MRF（ECLRMRF）模型的图像分割方法。通过对局部空间约束高斯混合模型（LSCGMM）引入邻接像素似然特征的空间约束相关性，提升图像的似然特征模型表达能力。部分二阶区域Potts 先验模型利用较大的区域信息显著提升局部区域的先验知识表达能力，使得Canny 边缘特征约束分割结果向图像分割边缘收敛。实验结果表明，本文所提模型不仅对图像噪声或者纹理突变信息的干扰具有较强的鲁棒性，而且能够有效保持图像分割结果的边缘信息，验证了该模型的可行性和有效性。Gibbs 采样算法在采样过程中由于计算信息增多而使得分割算法耗时较长，因此下一步将在保持分割结果准确度的前提下引入群体智能优化算法，以进一步提高计算效率。