[发明专利]一种基于局部灰度聚类特征的图像水平集分割方法

摘要

本发明提供一种基于局部灰度聚类特征的图像水平集分割方法，该方法为读取待分割图像，采用正交基函数的线性加权和拟合偏移场，并初始化各基函数的权重值，初始化图像的水平集函数集合，建立图像水平集分割的能量泛函，根据待分割图像设定水平集分割控制参数，分别更新聚簇中心集合、图像水平集函数集合、基函数权重列向量，直至满足迭代终止条件，得到迭代的能量泛函，根据当前更新后的图像水平集函数集合构造图像的隶属度函数，即待分割图像的分割结果，并根据更新后的基函数权重列向量与基函数列向量得到待分割图像的偏移场估计，该方法克服了弱边界、图像噪声以及灰度不一致性对图像分割精度的负面影响，起到校正图像灰度的作用。

主权项

1.一种基于局部灰度聚类特征的图像水平集分割方法，其特征在于，包括以下步骤：步骤1：读取待分割图像I；步骤2：采用M个正交基函数的线性加权和拟合偏移场b，并初始化各基函数的权重值；所述偏移场b的公式如下所示：b＝wTG；其中，w＝(w1，w2，...，wM)T为基函数权重列向量，G＝(g1，g2，...，gM)T为M个基函数构成的列向量，g1，g2，...，gM为两两正交的基函数，w1，w2，...，wM为初始化的各基函数的权重值；步骤3：初始化图像的水平集函数集合：根据待分割图像I待划分的区域数N，设定初始的k个零水平集分割曲线，根据设定的零水平集分割曲线生成高一维的水平集函数集合Φ＝(φ1，φ2，...，φk)；步骤4：建立图像水平集分割的能量泛函F(中，w，c)，根据待分割图像设定水平集分割控制参数：图像数据项权重系数λn、水平集函数正则项系数μ、水平集分割曲线平滑项系数v；所述能量泛函如下所示：其中，为图像数据项，n∈{1，2，...，N}，C为函数值为1的常值函数，*为卷积操作符号，K为归一化的非负偶函数，c＝(c₁，c₂，...，c_N)为聚簇中心集合，M_n(Φ(x))为第n个子区域Ω_n的隶属度函数，Φ(x)为x像素点处k个水平集函数值构成的集合；为水平集函数正则项，为φ_i(x)的梯度，为水平集分割曲线平滑项，δ(φ_i(x))为φ_i(x)的Dirac函数，为φ_i(x)的Heaviside函数；步骤5：设定最大迭代次数J和图像水平集分割的能量泛函F(Φ，w，c)的收敛阈值；步骤6：分别更新第j次迭代的聚簇中心集合cj、图像水平集函数集合Φj、基函数权重列向量wj，得到第j次迭代的能量泛函Fj(Φ，w，c)；步骤7：判断是否达到迭代终止条件：若当前迭代次数达到最大迭代次数或者Fj(Φ，w，c)与Fj‑1(Φ，w，c)的差值小于图像水平集分割的能量泛函F(Φ，w，c)的收敛阈值，则执行步骤8，否则，迭代次数加1，返回步骤6；步骤8：根据当前更新后的图像水平集函数集合构造图像的隶属度函数即待分割图像I的分割结果，并根据更新后的基函数权重列向量w^j与基函数列向量G得到待分割图像I的偏移场估计。