基于变分法的细胞演化建模
上QQ阅读APP看书,第一时间看更新

3.3.2 图像分割的Mumford-Shah模型和Chan-Vese模型

图像分割的目的是将图像中的灰度同质的区域分离出来,并通过各个同质区域的边界来表达.通常基于几何活动轮廓线模型分割图像的方法仅依赖于包围待分割目标或在目标内部的轮廓线所在位置的图像局部信息,难于综合图像区域的全局信息.如果轮廓线位于边缘尖锐度较弱的平滑边缘,则可能越过边缘,出现“冒顶”现象,且不会返回到正确位置.

这里介绍用于图像分割的Mumford-Shah模型和Chan-Vese模型(简称C-V模型).Mumford-Shah模型可以直接对图像进行运算,仅基于图像数据的驱动来完成分割,而不依赖于图像局部梯度信息.Mumford-Shah模型的能量函数包含了对图像区域、边界的描述,并且可以通过优化能量函数,一次获得图像(包括带噪声的)边界、区域以及平滑图像.

1)Mumford-Shah模型

Mumford和Shah[53]提出了一个关于图像分割的目标函数,并且通过泛函优化的方法进行图像分割.设图像Ix, y)的定义域为Ω,当前考察的图像边界C将图像Ix, y)划分为若干近似同值(质)区域,并得到分割图像Iox, y).Mumford-Shah图像分割模型就是通过寻找真正的图像边界Co,将图像Ix, y)划分为若干同值区域,并且所得分割图像I的误差比所有分割图像和原图像的误差还小,即最小化如下能量方程:

其中,

式中,λμ为正权数;dσ=dxdy为关于面积积分;|C|为曲线C的长度;最小点函数argmin参见式(3-13).当FMSIoC)最小时,所得边界Co将图像划分为若干个平滑区域,并且保留了尖锐的边界Co

2)简化的Mumford-Shah图像分割模型

对于复杂图像边界,长度|C|难以计算,所以Ambrosio等[54-55]提出了改进的能量函数,针对原模型使用各向同性线性扩散做恢复,Teboul等[56]提出利用非线性保边正则化的能量函数.

为方便计算,Chan和Vese提出一种利用水平集方法求解简化的Mumford-Shah模型的图像分割方法,推动了Mumford-Shah模型的应用.Chan和Vese[57-58]提出的简化Mumford-Shah模型原理如下:

假设图像中每个同值区域的灰度值都是常数,对于区域Di⊂(Ω\C),有IDi)=ci, ci是常数,则最小化能量泛函FMS的目的就是寻找最优化分割线(边界)Co,使得分割图像Io和原图像I之间的差异最小.此即为C-V模型原理

这里首先引入泛函来理解C-V模型.设原图像Ix, y)被活动轮廓线C划分为目标Ωo=inside (C)和背景Ωb=outside (C)两个区域,模型的拟合能量泛函如下:

式中,C是任意闭合活动轮廓线;cocb是依赖于C的两个常数.当闭合活动轮廓线C没有位于两个同值区域的边界Co时,FC)不能达到最小值.

据此,Chan和Vese提出了如下的图像分割能量函数:

式中,|C|是闭合轮廓线C的长度;μλoλbμ≥0, λo>0, λb>0)是各个能量项的权重系数;F的首项是平滑项.最终分割轮廓线C的位置以及未知数cocb经最优化式(3-22)得到,即

由于此模型利用了图像全图信息,因此通过最优化能量泛函,可以得到全局最优的图像分割结果.

3)基于水平集求解C-V模型

C-V模型的完整推导也需利用变分法,这里我们简化了理论推导,但其实质还是变分分析过程.

Φo是根据初始轮廓线Co构造的符号距离函数(signed distance function, SDF),即{Co|Φox, y)=0}.并设Φ为内正外负型的SDF,即

Φo[inside (C)]>0, Φo[outside (C)]<0

可以证明,以水平集函数表达的轮廓线C的长度为

上式中的Ω是水平集函数的定义域,赫维赛德(Heaviside)函数Hz)表示如下:

是狄拉克(Dirac)函数.以水平集函数ϕ表达为

下面关于cocbϕFϕ, cocb)的极小值.由于Fϕ, cocb)达到极小值时,一定是Fϕ, cocb)关于每一个变量的偏导数为零的时刻,因此,将ϕ固定,对能量函数Fϕ, cocb)关于co求偏导,得

整理得

同样对Fϕ, cocb)关于cb求导数,得

整理得

现固定cocb,对Fϕ, cocb)关于ϕ求导数,此时有

首先考虑

利用泰勒公式对上式进行展开,有

故有

记div为向量场的散度.由φ的任意性可知

为了计算式(3-17)至(3-23)的解,引入时间变量

式中,

从式(3-28)可看出,偏微分方程右边所涉及的图像函数Ix, y)范围是全图像,而不像测地线活动轮廓模型,仅利用轮廓线C所在位置的图像数据.式中的cocb也在图像定义域内,具有全局特征.因此,Chan和Vese强调该方法的一个特点就是全局优化,仅使用一条初始闭合轮廓线,就可把带有内部空洞的目标的内部检测出来,不需要为检测内部空洞而另做特别处理.

4)图像分割方法的数值解法

这里采用显式有限差分数值解法而不是C-V模型的隐式解法,对分割做介绍,以期式(3-28)的数值解法更好的稳定性.具体如下:

h为二维网格步长,(xi, yj)=(ih, jh), ij(1≤i, jM)为网格坐标,则是关于hϕx, y)在网格上的近似,有n≥0, ϕ0=ϕ0

引入如下差分格式记号:

根据Osher和Sethian[52]求解水平集的“熵守恒”差分方法,得到如下数值解表达式:

其中

Ki, j是水平集函数在(i, j)处的曲率,通过下式定义:

也可以由中心差分近似计算出来.

在实际算法应用中,为避免不必要的迭代计算,给出了如下迭代收敛判断条件:

式中,M是满足的网格点数目.当QQr时,迭代结束,Qr是门限.