CN103971363B - 基于分层模糊c均值的图像分割方法 - Google Patents

基于分层模糊c均值的图像分割方法 Download PDF

Info

Publication number
CN103971363B
CN103971363B CN201410127317.1A CN201410127317A CN103971363B CN 103971363 B CN103971363 B CN 103971363B CN 201410127317 A CN201410127317 A CN 201410127317A CN 103971363 B CN103971363 B CN 103971363B
Authority
CN
China
Prior art keywords
degree
membership
image
fuzzy
pixel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410127317.1A
Other languages
English (en)
Other versions
CN103971363A (zh
Inventor
张辉
陈北京
郑钰辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangzhou Kunchuan Network Technology Co ltd
Guangzhou Yuanbo Information Technology Co.,Ltd.
Original Assignee
Nanjing University of Information Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Information Science and Technology filed Critical Nanjing University of Information Science and Technology
Priority to CN201410127317.1A priority Critical patent/CN103971363B/zh
Publication of CN103971363A publication Critical patent/CN103971363A/zh
Application granted granted Critical
Publication of CN103971363B publication Critical patent/CN103971363B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种基于分层模糊c均值的图像分割方法，通过将均值模板同时应用于隶属度和距离函数这两项上，从而获得更好的图像分割结果，有效地利用了图像中的空间上下文信息，获得了更好的图像分割质量，同时算法的运行时间和计算量都较小，通过将距离函数采用分层模糊c模型，即将距离函数看作是子模糊c均值模型，其比传统的欧式距离函数具有更好的抗噪性本发明的算法对于图像噪声具有更好的鲁棒性。

Description

基于分层模糊c均值的图像分割方法
技术领域
本发明属于图像处理领域，特别涉及一种基于分层模糊c均值的图像分割方法。
背景技术
图像分割是图像处理中的重要研究课题，它决定着最终的图像分析和图像理解的结果和质量。正是由于图像分割的重要性，国内外许多学者对其展开了大量的研究，提出了各种各样不同的分割算法，但是这些算法大都是针对具体研究对象的，至今还没有通用的分割理论提出，因此人们仍然在不断地探索新的分割算法和分割理论。
模糊C均值算法是基于目标函数的聚类算法中理论最完善、应用最广泛的算法。由于成功地将模糊概念引入到图像像素的隶属度，模糊C均值算法能够保留更多的原始图像信息。虽然模糊C均值算法由于自身的优势得到了广泛的应用和发展，但是该算法在进行图像分割时需要事先确定分类数、对噪声敏感、易陷入局部极值。
Ahmed等通过引入空间约束修改经典的模糊c均值(Fuzzy C-means，FCM)算法的目标函数，取得了较好的图像分割结果，该算法称为MFCM算法。陈松灿等提出了关于MFCM算法的变化形式:FCM_S算法，并且加快了MFCM算法的执行速度。FLICM考虑了像素的局部空间信息而HMRF-FCM则考虑了像素空间关系的马尔科夫随机场(Markov Random Field,MRF)模型。然而，这些已有的研究算法具有计算复杂度高，计算时间长，依然容易受到图像噪声影响等弊端。
发明内容
本发明为了解决现有的算法具有计算复杂度高，计算时间长，依然容易受到图像噪声影响等弊端的问题，提出了一种基于分层模糊c均值的图像分割方法。
为解决上述问题，本发明采用的方法是：一种基于分层模糊c均值的图像分割方法，包括以下步骤：
(1)、将图像数据分到J个类中，在每个类j中，数据再被分到k个子类中，确定第一层的类数J，第二层的类数K，初始化隶属度子隶属度和初始化聚类中心μjk；
(2)、设置循环计数l＝0；
(3)、通过公式计算新的聚类中心
(4)、通过公式计算更新隶属度
(5)、通过公式计算更新子隶属度
(6)、如果目标函数收敛，则结束迭代估值，否则l＝l+1，重复步骤(3)-(5)。
作为本发明的一种优选，所述目标函数的推导过程包括如下步骤：
(1)、令yi表示图像中第i个点的像素值，其中i＝(1,2,…,N),N为图像的像素总个数，j(j＝1,2,…,K)表示像素点i所对应的类，模糊c均值的目标函数可以表示为：其中，uij为隶属度，表示像素yi属于第j类的可能性程度，m为权重指数，dij是点yi到聚类中心μj的距离测度，称为距离函数；
(2)、定义距离函数为一“子”模糊模型，基于标准FCM中的笛卡尔距离将“子”模糊模型表示为：其中是子距离函数，vijk是子隶属度；
(3)、将步骤(2)中的公式带入步骤(1)中的公式里，得到HFCM目标函数：
(4)、将步骤(3)中的HFCM目标函数结合均值模板的分层模糊c均值得到目标函数为：其中权重因子为高斯函数，δ＝(窗口大小-1)/4。
作为本发明的一种优选，所述的公式 是通过对目标函数求偏导得到。
作为本发明的一种优选，所述的初始化隶属度子隶属度和初始化聚类中心μjk通过k-means算法进行预处理得到。
有益效果：
1、本发明将均值模板同时应用于隶属度和距离函数这两项上，从而获得更好的图像分割结果，有效地利用了图像中的空间上下文信息，获得了更好的图像分割质量，同时算法的运行时间和计算量都较小。
2、本发明的距离函数采用分层模糊c模型，即将距离函数看作是子模糊c均值模型，其比传统的欧式距离函数具有更好的抗噪性。因此，我们的算法对于图像噪声具有更好的鲁棒性。
附图说明
图1本发明不同的方法分类三个数据点集的分类结果图。
图2采用本发明采用一幅人造的黑，灰，白三色图像进行分割的实验图。
图3采用本发明不同的方法对高斯噪声图像分割结果图。
图4采用本发明不同的方法对斑点噪声图像的分割结果图。
图5为伯克利数据库的原始图像集。
图6为采用本发明的算法得到的图像分割结果图。
图7为不同的方法对噪声图像的分割结果，误差率对比表格。
图8为基于伯克利数据集的不同方法的图像分割结果数据表。
具体实施方式
下面结合附图和具体实施方式，进一步阐明本发明。
现有的模糊c均值的图像分割方法主要采用以下的办法，令yi表示图像中第i个点的像素值，其中i＝(1,2,…,N),N为图像的像素总个数。j(j＝1,2,…,K)表示像素点i所对应的类。那么，模糊c均值(fuzzy c-means，FCM)的目标函数可以表示为：
其中，uij为隶属度，表示像素yi属于第j类的可能性程度。m为权重指数，dij是点yi到聚类中心μj的距离测度，称为距离函数。标准FCM中采用笛卡尔距离作为距离函数，表示为：
dij＝||yi-μj||2 (2)
FCM算法由最小化目标函数Jm得到以下的迭代公式：
其中
这里，我们介绍分层模糊c均值算法(Hierarchical Fuzzy c-Means,HFCM)。我们的算法简单，直观且有效。我们认为距离函数为一“子”模糊模型。因此，我们的方法具有很强的通用性，可以扩展到各种距离函数，如l1,lp,l∞范数和核距离函数等。为了明确阐述我们的算法，我们基于标准FCM中的笛卡尔距离来介绍我们的算法。因此，“子”模糊模型可以表示为：
这里，是子距离函数，vijk是子隶属度。
事实上，我们的HFCM可以看作是两层的FCM模型：在第一层中，数据被分到J个类中；而在第二层中，在每个类j中，数据再被分到k个子类中。将公式(4)代入(1)得到HFCM的目标函数为：
同样，我们可以使用经典的FCM迭代算法来对我们的HFCM做参数估计，得到：
我们考虑均值模板作为空间约束项，并且同时将其应用于隶属度和距离函数上。这样做的好处是考虑到了空间的上下文信息，因此对于噪声具有更强的鲁棒性。结合了均值模板的分层模糊c均值(GHFCM)的目标函数为：
其中权重因子为高斯函数，δ＝(窗口大小-1)/4。权重因子选择为高斯函数，像素点距离模板中心点越远，其权重因子越小。我们同样使用均值模板作用于隶属度上，并采用经典的FCM迭代算法对我们的GHFCM做参数估计，得到：
根据上述步骤计算推导出本发明采用的方法基于的目标函数，以及通过本发明的方法基于的目标函数推导出隶属度uij，子隶属度vijk和聚类中心μjk等参数的公式。然后结合这些公式完成图像的分割，具体过程如下：
[1]确定第一层的类数J，第二层的类数(子聚类数)K，初始化隶属度子隶属度和初始化聚类中心μjk，其中第一层的类数J根据待分割图像的颜色来确定，比如待分割图像有3种颜色，则J就取3，K的数值一般取2-3，一般来说，k值越大，相应的模型也越复杂，其解决问题的能力也越好。由于我们应用于图像处理，其噪声程度影响有限，因此我们取2-3，对于这个应用(图像分割)足够了。
[2]设置循环计数l＝0。
[3]通过公式(10)计算新的聚类中心
[4]通过公式(8)计算更新隶属度具体的说通过计算隶属度，通过更新隶属度。
[5]通过公式(9)计算更新子隶属度具体的说通过计算子隶属度，通过更新子隶属度。
[6]如果目标函数(公式(7))收敛，则结束迭代估值。否则l＝l+1，重复步骤[3]-[5]。
我们采用了三个实验来验证我们的算法的有效性，并且和现有的MFCM算法，FCM_S算法，FLICM算法，HMRF-FCM算法做比较，以验证我们的算法GHFCM的鲁棒性。
第一个实验用来验证HFCM算法的有效性。我们取1800个模拟数据点集，这些点由三类高斯分布组成，每组600个点。高斯分布的中心分别为(0,0),(3,2),(0,5)，高斯分布的方差分别为diag(1/2,1/2),diag(1/8,1/8),diag(1/2,1/2)。然后这些点集被1800个点的离群点(噪声)所干扰，这些离群点满足[-6 6]上的均匀分布。被噪声点干扰的数据点集如图1的a部分所示，三类数据点分别用绿色，红色，蓝色表示，离群点用黑色表示。采用FCM和HFCM进行分类的结果分别由图1中b部分和图1中c部分所示。从图1中b部分，我们可以看出FCM的分类结果并不好，部分红色和绿色点集混合在了一起，并且其分类误差率(MCR)也较大，为5.33％。而从图1中c部分中我们可以看出HFCM获得了较好的分类结果，绿色，红色，蓝色都得到了正确的分类，并且其误差率也较小，为0.72％。这说明了我们的方法的有效性。
第二个是对一幅人造的黑，灰，白三色图像进行分割的实验。此模拟图像如图2中a部分所示。受到高斯噪声和斑点噪声的污染图像分别如图2中b部分和c部分所示。不同的方法对高斯噪声图像和斑点噪声图像的分割结果分别如图3和图4所示。其中图3中的(a)部分采用MFCM方法,其分类误差率MCR＝29.7％；(b)部分采用FCM_S方法,其分类误差率MCR＝13.11％；(c)部分采用FLICM算法,其分类误差率MCR＝9.39％；(d)部分采用HMRF-FCM算法,其分类误差率MCR＝11.3％；(e)部分采用本发明的GHFCM,其分类误差率MCR＝3.66％，可以看出我们的算法GHFCM获得了最好的分割结果以及最小的误差率。
各种方法对于不同噪声强度的图像分割结果如图7所示。图4中的(a)部分采用MFCM方法,其分类误差率MCR＝15.05％；(b)部分采用FCM_S方法,其分类误差率MCR＝5.48％；(c)部分采用FLICM算法,其分类误差率MCR＝9.05％％；(d)部分采用HMRF-FCM算法,其分类误差率MCR＝7.94％；(e)部分采用本发明的GHFCM,其分类误差率MCR＝2.36％，同样可以看出我们的GHFCM具有最小的误差率，以及最强的抗噪声性和鲁棒性。
第三个实验采用著名的美国伯克利大学图像数据库来验证我们的算法。部分实验图像如图5所示。采用我们的算法GHFCM得到的图像分割结果如图6所示。我们采用PRI值(Probabilistic Rand Index)来评估实验结果，PRI越大，表示分割结果越好。相关的实验结果如图8中的表可以看出，和现有的方法相比，我们的方法具有最优的图像分割结果(最大PRI值)。
本发明方案所公开的技术手段不仅限于上述技术手段所公开的技术手段，还包括由以上技术特征任意组合所组成的技术方案。以上所述是本发明的具体实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。

Claims (3)

1.一种基于分层模糊c均值的图像分割方法，其特征在于：包括以下步骤：
(1)、令yi表示图像中第i个点的像素值，其中i＝1,2,…,N,N为图像的像素总个数，将图像数据分到J个类中，j＝1,2,…,J表示像素点i所对应的类，在每个类j中，数据再被分到k个子类中，k＝1,2,…,K表示像素点i所对应的子类，确定第一层的类数J，第二层的类数K，初始化隶属度子隶属度和初始化聚类中心μjk；
(2)、设置循环计数l＝0；
(3)、通过公式计算新的聚类中心
(4)、通过公式计算更新隶属度
(5)、通过公式计算更新子隶属度
(6)、如果目标函数收敛，则结束迭代估值，否则l＝l+1，重复步骤(3)-(5)；所述目标函数的推导过程包括如下步骤：
(1)、令yi表示图像中第i个点的像素值，其中i＝1,2,…,N,N为图像的像素总个数，j＝1,2,…,J表示像素点i所对应的类，模糊c均值的目标函数表示为：其中，uij为隶属度，表示像素yi属于第j类的可能性程度，m为权重指数，dij是点yi到聚类中心μj的距离测度，称为距离函数；
(2)、定义距离函数为一“子”模糊模型，基于标准FCM中的笛卡尔距离将“子”模糊模型表示为：其中是子距离函数，vijk是子隶属度；
(3)、将步骤(2)中的公式带入步骤(1)中的公式里，得到HFCM目标函数：
(4)、将步骤(3)中的HFCM目标函数结合均值模板得到分层模糊c均值的目标函数为：其中权重因子为高斯函数，δ＝(窗口大小-1)/4。
2.根据权利要求1所述的基于分层模糊c均值的图像分割方法，其特征在于：所述的公式 是通过对目标函数求偏导得到。
3.根据权利要求1所述的基于分层模糊c均值的图像分割方法，其特征在于：所述的初始化隶属度子隶属度和初始化聚类中心μjk通过k-means算法进行预处理得到。
CN201410127317.1A 2014-03-31 2014-03-31 基于分层模糊c均值的图像分割方法 Active CN103971363B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410127317.1A CN103971363B (zh) 2014-03-31 2014-03-31 基于分层模糊c均值的图像分割方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410127317.1A CN103971363B (zh) 2014-03-31 2014-03-31 基于分层模糊c均值的图像分割方法

Publications (2)

Publication Number Publication Date
CN103971363A CN103971363A (zh) 2014-08-06
CN103971363B true CN103971363B (zh) 2017-09-05

Family

ID=51240812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410127317.1A Active CN103971363B (zh) 2014-03-31 2014-03-31 基于分层模糊c均值的图像分割方法

Country Status (1)

Country Link
CN (1) CN103971363B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105976373A (zh) * 2016-05-05 2016-09-28 江南大学 基于邻域信息熵的核模糊c均值图像分割算法
CN106570871B (zh) * 2016-11-02 2019-10-18 天津大学 模糊c均值颈动脉超声图像内中膜厚度测量方法及系统
CN106780515B (zh) * 2017-01-04 2019-09-27 南京审计大学 一种脑部核磁共振图像中脑胶质瘤区域提取方法
CN109300137B (zh) * 2018-09-20 2021-05-07 北京航空航天大学 一种多曲面估计区间二型模糊聚类磁共振脑图像分割方法
CN113688960A (zh) * 2021-10-27 2021-11-23 南昌工程学院 基于灰狼优化ghfcm的居民电力数据聚类方法及装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101980298A (zh) * 2010-11-04 2011-02-23 西安电子科技大学 基于多智能体遗传聚类算法的图像分割方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007001236A1 (en) * 2005-06-28 2007-01-04 Agency For Science, Technology And Research Wavelet shrinkage pre-filtering of mr images for brain tissue segmentation

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101980298A (zh) * 2010-11-04 2011-02-23 西安电子科技大学 基于多智能体遗传聚类算法的图像分割方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A Fast Kernel-based Clustering Algorithm with Application in MRI Image Segmentation;Liang Liao;《2009 IEEE International Conference on Intelligent Computing and Intelligent Systems》;20091120;第4卷;第405-410页 *
基于FCM聚类的彩色图像分割算法研究;朱小雨;《中国优秀硕士学位论文全文数据库信息科技辑》;20120515(第05期);摘要，第23-32页 *

Also Published As

Publication number Publication date
CN103971363A (zh) 2014-08-06

Similar Documents

Publication Publication Date Title
CN103971363B (zh) 基于分层模糊c均值的图像分割方法
CN107644429B (zh) 一种基于强目标约束视频显著性的视频分割方法
CN104182772B (zh) 一种基于深度学习的手势识别方法
CN108875935B (zh) 基于生成对抗网络的自然图像目标材质视觉特征映射方法
CN108399428B (zh) 一种基于迹比准则的三元组损失函数设计方法
CN102930246B (zh) 一种基于点云片段分割的室内场景识别方法
CN108932536A (zh) 基于深度神经网络的人脸姿态重建方法
US7536064B2 (en) Image comparison by metric embeddings
CN105654453B (zh) 一种鲁棒性的fcm图像分割方法
US20100098331A1 (en) System and method for segmenting foreground and background in a video
CN107153816A (zh) 一种用于鲁棒人脸识别的数据增强方法
CN103824049A (zh) 一种基于级联神经网络的人脸关键点检测方法
CN108615229B (zh) 基于曲率点聚类及决策树的碰撞检测优化方法
JPH06150000A (ja) 画像クラスタリング装置
CN109829449A (zh) 一种基于超像素时空上下文的rgb-d室内场景标注方法
CN110516537B (zh) 一种基于自步学习的人脸年龄估计方法
CN102542593A (zh) 一种基于视频解译的交互式视频风格化渲染方法
CN109816714B (zh) 一种基于三维卷积神经网络的点云物体类型识别方法
Sandhu et al. A nonrigid kernel-based framework for 2D-3D pose estimation and 2D image segmentation
CN112614070B (zh) 一种基于DefogNet的单幅图像去雾方法
WO2021068061A1 (en) System and method for generating 3d models from specification documents
CN111507334A (zh) 一种基于关键点的实例分割方法
Zhao et al. Deep Adaptive Log‐Demons: Diffeomorphic Image Registration with Very Large Deformations
CN107146251A (zh) 一种三维人脸模型的对称性分析方法
CN116682021A (zh) 一种高分辨率遥感影像建筑物矢量轮廓数据提取方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20190329

Address after: 510000 Room 3301-3302, No. 1 Jinsui Road, Tianhe District, Guangzhou City, Guangdong Province

Patentee after: GUANGZHOU FENGSHEN NETWORK TECHNOLOGY Co.,Ltd.

Address before: 510000 2414-2416 of the main building 371, five mountain road, Tianhe District, Guangzhou, Guangdong.

Patentee before: GUANGDONG GAOHANG INTELLECTUAL PROPERTY OPERATION Co.,Ltd.

Effective date of registration: 20190329

Address after: 510000 2414-2416 of the main building 371, five mountain road, Tianhe District, Guangzhou, Guangdong.

Patentee after: GUANGDONG GAOHANG INTELLECTUAL PROPERTY OPERATION Co.,Ltd.

Address before: 210000 No. 219 Ningliu Road, Pukou District, Nanjing City, Jiangsu Province

Patentee before: Nanjing University of Information Science and Technology

TR01 Transfer of patent right
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201012

Address after: Room 405, no.602-1, Guangshan 2nd Road, Tianhe District, Guangzhou City, Guangdong Province 510000

Patentee after: Guangzhou Fengshen Intelligent Technology Co.,Ltd.

Address before: 510000 Room 3301-3302, No. 1 Jinsui Road, Tianhe District, Guangzhou City, Guangdong Province

Patentee before: GUANGZHOU FENGSHEN NETWORK TECHNOLOGY Co.,Ltd.

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: Room 3301, No.1 Jinsui Road, Tianhe District, Guangzhou City, Guangdong Province, 510000 (for office purposes only)

Patentee after: Guangzhou Yuanbo Information Technology Co.,Ltd.

Country or region after: China

Address before: 510000 Room 405, no.602-1, Guangshan 2nd Road, Tianhe District, Guangzhou City, Guangdong Province

Patentee before: Guangzhou Fengshen Intelligent Technology Co.,Ltd.

Country or region before: China

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240515

Address after: Room 499, 4th Floor, No. 89 Yanling Road, Tianhe District, Guangzhou City, Guangdong Province 510000. Self made No. 134 (for office only)

Patentee after: Guangzhou Kunchuan Network Technology Co.,Ltd.

Country or region after: China

Address before: Room 3301, No.1 Jinsui Road, Tianhe District, Guangzhou City, Guangdong Province, 510000 (for office purposes only)

Patentee before: Guangzhou Yuanbo Information Technology Co.,Ltd.

Country or region before: China