1. 引言
1.1. LBF模型的重要性与应用
局部二值拟合(Local Binary Fitting, LBF)模型作为图像分割领域的重要工具,由Li等人提出,因其能够有效处理图像灰度不均匀的特性而受到广泛关注;LBF模型通过引入局部区域信息,该模型通过引入高斯核函数,在局部区域内拟合轮廓内外两侧的灰度均值,克服了传统全局模型(如CV模型)在分割非均匀图像时的局限性,处理非匀质图像时表现出更高的鲁棒性与准确性,从而在演化过程中更准确地捕捉感兴趣区域(the Region of Interest, ROI)的边界特征;LBF模型在医疗影像处理、普通图像切割、人工智能及工业质检等多领域表现出优异性能;以医学图像为例,该模型能有效区分边界模糊且对比度较低的目标区域,这对临床诊断具有重要价值;具体实际应用充分证明:LBF模型在提升分割准确度和系统稳定性方面优势明显。
然而,尽管LBF模型在理论上具备诸多优点,其实际应用仍面临一些挑战。特别是在复杂场景下,当图像包含强噪声或对比度极低的目标时,LBF模型的表现往往受到限制,难以捕捉到ROI;这些问题促使研究者不断探索改进的方法,以进一步提升其性能。因此,基于改进LBF模型来提升分割效果不仅是当前图像分割研究的核心之一,也是推动相关领域技术发展的关键驱动力。
1.2. LBF模型初始轮廓敏感性问题
LBF模型虽然在灰度不均匀图像分割中表现出色,但其对初始轮廓位置的高度依赖性成为制约其性能的重要因素。具体而言,LBF模型的演化过程严重依赖于初始轮廓的设定,当初始轮廓远离目标边界时,模型可能无法准确收敛到真实边界位置。这种现象的成因在于,LBF模型的能量函数设计主要基于局部区域信息,而缺乏对全局图像特征的充分利用。因此,当初始轮廓偏离目标较远时,局部信息难以提供足够的引导力,导致演化曲线无法正确逼近目标边界。
在实际应用中,初始轮廓敏感性问题带来了诸多困扰。例如,在医学影像分割中,由于图像通常包含复杂的背景和噪声,手动设定精确的初始轮廓不仅耗时耗力,而且容易引入人为误差。此外,在自动驾驶和环境感知等领域,目标物体的位置和形状变化多样,难以预先确定合适的初始轮廓。这些问题使得LBF模型在实际场景中的应用受到限制,亟需通过改进算法设计来降低其对初始轮廓的依赖性。
1.3. 研究目的与意义
T. Chan和L. Vese [1]提出最小化Mumford-Shah函数的水平集算法;T. Wang等[2]提出了由局部高斯分布拟合能量驱动的主动轮廓模型;Li等[3]提出了一个具有局部强度拟合项的能量泛函数和一个辅助全局强度拟合项的模型;Chunming Li和Chiu_Yen [4]等提出最小化区域–可扩展拟合能量用于图像分割;Zhang等[5]结合改进CV模型,提出了一种融合局部和全局图像信息的活动轮廓模型。Wang等[6]结合CV与LBF模型的优点,提出了局部和全局灰度拟合模型,该模型能够分割灰度不均匀图像,且对初始轮廓的敏感性降低,但其权重系数受初始轮廓位置影响,限制了模型的应用;任鸽等[7]研究结合全局和局部信息的“两阶段”活动轮廓模型,王顺凤等[8]利用CV模型对图像目标进行定位,再利用LBF模型进一步精确化第一阶段的分割结果,这一类方法改进了CV模型和LBF模型的分割效果,但是分割精度还有待改善;Liu等[9]等构造新的局部和全局能量项,Fang等[10]基于迭代次数设计自适应函数,实现全局与局部能量权重系数的自动选取。
本研究旨在通过结合全局信息对LBF模型进行改进,以降低其对初始轮廓的敏感性并提高演化速度。具体而言,我们将在LBF模型的基础上融入全局图像特征,通过设计新的能量函数和演化方程,使模型能够在演化过程中充分利用全局和局部信息的优势[8];改进后的模型对于提升图像分割效果具有重要意义。首先,通过引入全局信息,模型能够在初始轮廓设定不够精确的情况下仍能准确收敛到目标边界,从而减少人工干预的需求。其次,全局信息的引入可以提高模型对噪声和复杂背景的鲁棒性,使其在低对比度图像分割中表现更加稳定。最后,该改进方法对弱边界图像进行分割能取得较好的效果,而且在时间效率上也有所提高。
2. 相关工作
2.1. 水平集算法在图像分割中的应用原理
图像分割问题存在多种数学建模方式。本文将重点阐述介绍其中一种较为实用的模型方法M-S模型[1] (由Mumford以及Shah于1989年提出)。该模型将观察到的电子图像定义在区域
上二维函数
来表示,并且假设
,其中
是区域
边界形成的一条曲线,
是一个正整数。图像分割问题可转化为求解最优解:寻找一个光滑函数连续可微函数
,使下面函数
(1)
满足特定泛函的极值条件;当中
表示曲线
的长度,
是预设的常数。该优化问题仍具挑战性;为此,可将图像分割任务转化为寻求分段常值函数的最优解:通过找到一个分段单值函数,即在
上
为常值
,使下面函数
(2)
达到最小值。
2.2. CV模型
2.2.1. CV模型理论
2001年,Chan与Vese共同建立了CV分割模型[2],该算法基于目标与背景存在显著差异的假设,采用水平集技术来最小化能量函数;其核心能量函数表达式如下所示:
(3)
式(3)中,
表示图像区域,
表示水平集函数,
表示输入图像,
分别表示演化曲线内部和外部的平均灰度值,
和
分别表示Dirac、Heaviside函数,
和
均设为正的权重系数,利基于梯度下降法,可推导出CV模型的演化方程:
(4)
其中:
(5)
(6)
式(6)中,
为一个全局量,
代表演化曲线的内部灰度平均值,
代表演化曲线的外部灰度平均值;在实际求解时,Heaviside函数
通常用一个函数
来代替:
(7)
同样的,对Dirac函数也通常用一个函数来代替:
(8)
2.2.2. CV模型实验结果与分析
在图1中三组图像分别展示了不同灰度特征的样本:首列为均匀分布,次列呈现非均匀特性,末列则显示目标与背景灰度相似。实验结果表明,CV算法在均匀灰度样本上分割精度较高,但在后两种情况下表现欠佳。CV模型作为依赖区域特征的全局性算法,通过分析像素灰度统计特性,在处理灰度分布均匀的图像时能获得理想的分割效果。然而当面对灰度不均的图像时,该模型不仅收敛速度迟缓,分割精度也会显著下降。
Figure 1. Experimental results of the CV model
图1. CV模型实验结果
2.3. LBF模型
2.3.1. LBF模型理论
与CV模型不同LBF (Local Binary Fitting)模型通过引入局部图像信息来处理强度不均匀性,能够更好地适应图像中的灰度变化,其实现方式采用高斯核函数。设
表示图像域
中的待分割图像,轮廓曲线
由一个Lipschitz函数的零水平集表示,水平集函数
定义如下:
(9)
在方程(9)中定义了水平集函数
后,LBF模型中的能量泛函定义为:
(10)
其中,
表示该模型的能量泛函,表达式如下:
(11)
式(11)中,
,
和
为均是正的权重系数,
是与点
无关的空间像素,
为输入图像,
表示标准差为
的高斯核函数,
和
则是近似局部的两个拟合函数分别在轮廓内外的灰度平均值。
表示
的梯度。
和
为函数,其中
为带有小常数
的Heaviside函数,而Dirac函数
为
的导数。两个函数
和
写为:
(12)
(13)
当保持
值固定时,基于变分法最小化能量函数式得到:
(14)
(15)
分析表明,局部拟合函数
和
会依据空间位置差异产生变化,这使得LBF模型在处理灰度分布不均的图像时表现优异。该模型仅聚焦于局部特征,未纳入全局信息考量;鉴于图像各区域间存在显著的灰度差异,分割效果很大程度上取决于初始轮廓线的位置和尺寸;相较于C-V方法,LBF算法基于非均匀强度假设;然而每次迭代仍需执行四次卷积操作,运算开销较大。同时其分割效果易受初始位置影响,且难以适应光照变化。
2.3.2. LBF模型实验结果与分析
在图2中,前两幅呈现灰度分布不均的特征,第三幅则是边界不清晰的非均匀灰度图像。LBF算法对前两种图像的分割效果良好,但对第三种图像的处理效果欠佳。研究数据显示,LBF模型因其出色的局部特征提取能力,有效解决了传统区域强度不均匀方法的不足。该模型不仅运算速度快、结果精确,在处理灰度分布不均的图像时表现优异。然而,当面对强度差异极大的图像时,其分割性能会明显下降。
Figure 2. LBF experimental model
图2. LBF实验模型
3. 基于全局信息改进LBF模型的实现
3.1. G-LBF模型(Global-Local Binary Fitting)
基于上述不足,本研究提出了一种结合全局信息和局部信息的改进LBF模型[7] [10]。该模型通过融入全局信息,不仅能够提供更准确的ROI,还能有效降低对初始轮廓位置的依赖性。同时,通过优化能量函数和演化方程,本研究旨在进一步提高模型的演化速度,从而提升整体分割效率。该模型即G-LBF模型,将CV模型和LBF模型能量项相加,是一种基于局部区域与全局区域的活动轮廓模型。该模型的能量泛函是基于CV和LBF两个模型局部能量泛函的组合,公式如下:
(16)
其中,
是正常量,
。将
固定,然后对
极小化,得到它们的表达式;最后利用变分原理及梯度下降流得到该模型的水平集函数演化方程:
(17)
其中:
(18)
(19)
在G-LBF模型中,通过
值来调节图像全局和局部的权重比例。该模型从两个维度利用全局信息:首先,借助整体统计建模提升对图像宏观结构的认知;其次,通过构建全局能量函数来改进演化过程的优化效果。以水平集活动轮廓模型为例,该方法通过整合全局特征构建能量泛函的整体项,使算法在迭代时能同时考量目标与背景的宏观差异,显著提升边界定位精度。此外,将全局特征与局部特征相融合,可有效克服单一信息源的局限性,从而优化模型的分割效果。
3.2. 关键参数θ分析
在能量函数(式(16))中,参数
扮演着平衡全局信息与局部信息贡献的关键角色:
当
时,模型退化为经典的Chan-Vese (CV)模型,完全依赖图像全局统计信息;此时模型对灰度均匀但边界模糊的图像有效,但对灰度不均匀区域敏感。
当
时,模型退化为经典的局部二值拟合(LBF)模型,完全依赖局部高斯拟合信息;此时模型能有效处理灰度不均匀,但对噪声和初始轮廓位置更敏感。
当
时,模型是CV与LBF的线性组合,其驱动力(式(17))是全局力
与局部力
的加权和;
的值决定了模型在“全局一致性”与“局部适应性”之间的权衡。
为了直观展示
对分割结果的影响,我们选取了具有不同特征的典型图像进行实验,固定其他参数(
,
,
,核尺度
),仅改变
值。
对于存在明显灰度不均匀的图像,
应设置为一个较小的正值(如0.1~0.4),让局部力主导,同时借助少量全局力提升稳定性;对于高噪声图像,应适当增大
值(如0.4~0.7),增强全局力的作用以抑制噪声,但需权衡细节保持能力;对于弱边界图像,需要综合局部与全局信息,通常设置
在0.5左右能取得较好平衡。本模型设置
。
3.3. 本模型算法设计
改进的模型算法具体步骤见下表1。
Table 1. Algorithmic steps of the G-LBF model
表1. G-LBF模型算法步骤
Step 1. 输入图像I,选择一个初始轮廓位置并初始化水平集函数
; |
Step 2. 设置算法所需参数,计算图像的全局灰度均值
和
,以及局部灰度均值
和
; |
Step 3. 计算全局能量项
和局部能量项
的值; |
Step 4. 迭代计算逐步更新水平集函数
; |
Step 5. 依据分割结果判断是否满足停止条件,满足则停止分割;否则返回Step 2。 |
4. 实验结果与分析
本文采用新提出的模型,其与LBF方法进行性能比较,重点考察初始轮廓差异对最终分割结果与ROI差距的影响;本文的评价算法评估体系主要包含两个指标:算法收敛所需的迭代次数和从轮廓初始化到分割结束所消耗CPU的处理时间;
针对一张存在轻微灰度不均匀的人工合成图像进行分析,图3显示光照等因素导致目标边缘产生阴影效应;下面使用本文模型与LBF模型通过设置不同初始轮廓进行分割实验,以评估本文模型对于初始轮廓的敏感程度;不同初始轮廓的位置如图3的第一列中红色曲线表示,每行代表不同的初始轮廓;图3的第二列为LBF模型迭代停止的分割结果;图3的第三列为本文模型迭代停止的分割结果。
Figure 3. Comparison of segmentation results with different initial contours
图3. 不同初始轮廓分割结果对比
Table 2. Iteration counts and CPU processing time (seconds) of different models
表2. 不同模型迭代次数及CPU处理时间(秒)
初始轮廓 |
LBF模型迭代次数 |
LBF模型CPU
处理时间 |
本文模型迭代次数 |
本文模型CPU
处理时间 |
轮廓1 |
30 |
3.031037 |
10 |
0.467690 |
轮廓2 |
30 |
2.995676 |
10 |
0.446348 |
轮廓3 |
30 |
2.960719 |
10 |
0.498609 |
为检验本章算法对初始轮廓的适应能力,采用人工合成图像与LBF模型进行对比分析。如图3所示,实验设置了三种差异化的初始轮廓:完全包覆目标的大矩形轮廓、局部覆盖目标的小矩形轮廓,以及远离目标区域的小矩形轮廓。实验数据表明,初始轮廓的差异会显著影响LBF模型的性能表现,当仅能覆盖目标局部区域时,该模型容易与ROI相差甚远;而本文模型基本与ROI密切吻合;从表2展示不同模型迭代次数与运算耗时对比可见,本方案的计算效率明显优于LBF方法。
5. 结束
笔者提出了一种融合全局信息和局部信息活动轮廓模型对灰度不均的图像有较好的分割效果,增强演化曲线对目标边界的捕捉能力和对LBF模型进行改进,在能量函数定义及演化方程中融入全局信息,使模型在演化过程中能充分利用图像的宏观特征;实验结果表明,改进后的模型显著克服了对初始轮廓的敏感性,即使在初始轮廓设置不理想的情况下,仍能有效收敛到目标边界;同时,模型的演化速度得到明显提升,大幅地缩短了分割所需时间。