CN111402272B - 一种邻域加权模糊聚类图像分割方法 - Google Patents
一种邻域加权模糊聚类图像分割方法 Download PDFInfo
- Publication number
- CN111402272B CN111402272B CN202010197180.2A CN202010197180A CN111402272B CN 111402272 B CN111402272 B CN 111402272B CN 202010197180 A CN202010197180 A CN 202010197180A CN 111402272 B CN111402272 B CN 111402272B
- Authority
- CN
- China
- Prior art keywords
- pixel
- neighborhood
- objective function
- fuzzy
- formula
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开一种邻域加权模糊聚类图像分割方法,计算图像中每一个像素的邻域像素权重,并将邻域像素权重融入到模糊C均值聚类算法的目标函数中,得到目标函数JNWFCM,对目标函数JNWFCM进行优化得到使目标函数JNWFCM为局部最小值的每一个像素的模糊隶属度值,利用每一个像素的模糊隶属度值对图像进行分割。本发明的优点在于,能够有效改善FCM算法用于图像分割时对噪声的敏感性,明显提高分割的准确率。
Description
技术领域
本发明属于图像处理与模式识别领域,具体涉及一种邻域加权模糊聚类图像分割方法。
背景技术
图像分割是指根据图像的灰度、颜色或者纹理等特征将图像划分成若干个互不相交的区域,使得这些特征在每个区域内表现出一致性或相似性。图像分割是从图像处理到图像分析的关键步骤,也是图像理解的前提和基础。
模糊C均值聚类(Fuzzy C-means,FCM)算法是一种无监督的聚类方法,近年来被广泛地应用于图像分割领域。然而FCM算法用于图像分割时没有考虑像素的空间位置关系,因此FCM算法对噪声特别敏感,分割结果会产生许多孤立像素。为了克服FCM算法用于图像分割时对噪声的敏感性,提高分割的准确率,研究者对该问题进行了比较深入的研究,提出了一些方法,但是这些方法在不同程度上仍然受到噪声的影响。因此,需要探索新的方法,提高FCM用于图像分割时准确率。
发明内容
为了克服现有技术中存在的不足,本发明的目是提供一种邻域加权模糊聚类(NWFCM)图像分割方法,本发明所提供的图像分割方法能够提高图像分割的准确率和对噪声的鲁棒性。
为了实现上述目的,本发明所采用的技术方案是:一种邻域加权模糊聚类图像分割方法,计算图像中每一个像素的邻域像素权重,并将邻域像素权重融入到模糊C均值聚类算法的目标函数中,得到目标函数JNWFCM,对目标函数JNWFCM进行优化得到使目标函数JNWFCM为局部最小值的每一个像素的模糊隶属度值,利用每一个像素的模糊隶属度值对图像进行分割。
进一步的,具体包括以下步骤:
步骤一,设定聚类数目c、模糊加权指数m和迭代停止阈值εL;
步骤二,采用随机化的方法初始化聚类中心vk;
步骤三,计算每一个像素的邻域像素权重Wi;
步骤四,将步骤三计算得到的权重Wi={wir,r∈Ni}融入到FCM算法的目标函数中,得到目标函数JNWFCM为
式(5)中,c为聚类数目,uik为第i个像素属于第k类的模糊隶属度值,m为模糊加权指数,控制聚类结果模糊程度,vk为第k个类的中心,Ni为第i个像素的邻域,xr为像素灰度值;
步骤五,对步骤四得到的目标函数JNWFCM进行优化得到计算像素的模糊隶属度值的式(6)和计算聚类中心的式(7);
步骤六,对步骤五中的式(6)和式(7)进行迭代,直至(8)成立,则得到使目标函数JNWFCM为局部最小值的每一个像素的模糊隶属度值uik,否则继续迭代。
||vold-vnew||≤εL (8)
式(8)中,vnew为当前迭代的聚类中心,vold为前一步迭代的聚类中心;
步骤七,根据步骤六得到的uik,代入式(9)计算每个像素的分类标签,完成图像分割。
进一步的,所述步骤三中每一个像素的邻域像素权重Wi计算步骤如下:
首先,设X={x1,x2,...,xn}表示具有n个像素的图像,xi为第i个像素的灰度值,像素xi的邻域像素的灰度均值μi为
式(1)中,Ni为像素xi的邻域,NR为集合Ni的势,xr为像素xi的邻域像素,xr和μi之间的L1范数为
σir=|xr-μi| (2)
定义像素xi的邻域像素xr(r∈Ni)的权重为
式(3)中β为伸缩因子,这里取β=3;σi为σir(r∈Ni)的最大值,即对ξir进行归一化处理得权重wir:
记像素xi的邻域像素Ni的权重为Wi={wir,r∈Ni}。
进一步的,所述步骤五中采用拉格朗日乘子法进行优化计算。
进一步的,所述步骤三中的邻域为二阶邻域,在二阶邻域中,每个像素有8个邻点。
与现有技术相比,本发明至少具有以下有益效果,本发明定义了像素xi的邻域像素xr(r∈Ni)的权重wir,当像素xi的邻域像素xr属于噪声或边缘像素时,wir较小,反之较大;因此该方法能够检测邻域像素特征的各向异性;在此基础上,将权重wir融入FCM的目标函数中,对FCM算法引入了邻域约束,从而提供了一种邻域加权模糊聚类图像分割方法,本发明能够有效改善FCM算法用于图像分割时对噪声的敏感性,明显提高分割的准确率。
附图说明
图1为两个邻域像素及其权重;
图2为人工合成图像分割结果对比,其中,(a)为原始图像,(b)为加高斯噪声和椒盐噪声的混合噪声后的图像,(c)为FCM分割结果,(d)为FCM_S1分割结果,(e)为FLICM分割结果,(f)为NWFCM分割结果;
图3为真实图像分割结果对比,(a)为原始图像,(b)为加高斯噪声和椒盐噪声的混合噪声后的图像,(c)为FCM分割结果,(d)为FCM_S1分割结果,(e)为FLICM分割结果,(f)为NWFCM分割结果。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步描述。
计算图像中每一个像素的邻域像素权重,并将邻域像素权重融入到模糊C均值聚类(Fuzzy C-means,FCM)算法的目标函数中,得到目标函数JNWFCM,对目标函数JNWFCM进行优化得到使目标函数JNWFCM为局部最小值的每一个像素的模糊隶属度值,利用每一个像素的模糊隶属度值对图像进行分割;
具体包括以下步骤:
步骤一,设定聚类数目c、模糊加权指数m和迭代停止阈值εL;
步骤二,采用随机化的方法初始化聚类中心vk;
步骤三,计算每一个像素的邻域像素权重Wi;
首先,设X={x1,x2,...,xn}表示具有n个像素的图像,xi为第i个像素的灰度值,像素xi的邻域像素的灰度均值μi为
式(1)中,Ni为像素xi的邻域,NR为集合Ni的势,xr为像素xi的邻域像素,xr和μi之间的L1范数为
σir=|xr-μi| (2)
定义像素xi的邻域像素xr(r∈Ni)的权重为
式(3)中β为伸缩因子,这里取β=3;σi为σir(r∈Ni)的最大值,即对ξir进行归一化处理得权重wir:
记像素xi的邻域像素的权重为Wi={wir,r∈Ni}。
步骤四,将步骤三计算得到的权重Wi={wir,r∈Ni}融入到FCM算法的目标函数中,得到目标函数JNWFCM为
式(5)中,c为聚类数目,uik为第i个像素属于第k类的模糊隶属度值,m为模糊加权指数,控制聚类结果模糊程度,vk为第k个类的中心,Ni为第i个像素的邻域;
步骤五,采用拉格朗日乘子法对步骤四得到的目标函数JNWFCM进行优化得到计算像素的模糊隶属度值的式(6)和计算聚类中心的式(7);
步骤六,对步骤五中的式(6)和式(7)进行迭代,直至(8)成立,则得到使目标函数JNWFCM为局部最小值的每一个像素的模糊隶属度值uik,否则继续迭代;
||vold-vnew||≤εL (8)
式(8)中,vnew为当前迭代的聚类中心,vold为前一步迭代的聚类中心;
步骤七,根据步骤六得到的uik,代入式(9)计算每个像素的分类标签,完成图像分割。
如图1所示,从人工合成图像中选择的两个邻域及其权重,其中,(a)为选自边缘区域的邻域,(b)选自平滑区域的邻域,(a)和(b)的左边均为像素灰度值,右边均为相应的权重,图中边缘像素和噪声像素为粗体,从图1可以发现,当前邻域中,如果像素为边缘像素或者噪声像素,则该像素的相应的权重较小,因此该方法可以有效检测邻域像素的各项异性。
如图2和图3所示,其中,图2(a)和图3(a)均为原始图像,图2(b)和图3(b)均为加高斯噪声和椒盐噪声的混合噪声后的图像,图2(c)和图3(c)均为FCM分割结果,图2(d)和图3(d)均为FCM_S1分割结果,图2(e)和图3(e)均为FLICM分割结果,图2(f)和图3(f)均为NWFCM分割结果;从图中可以看出,FCM、FCM_S1、FLICM算法的分割结果不同程度受到噪声的影响,而本发明提供的NWFCM方法能够得到较好的分割效果。
表1为图1(a)人工合成图像加不同噪声后的分割准确率,从表1可以看出本发明提供的算法分割准确率最高。
表1四种方法分割准确率对比
Claims (3)
1.一种邻域加权模糊聚类图像分割方法,其特征在于,计算图像中每一个像素的邻域像素权重,并将邻域像素权重融入到模糊C均值聚类算法的目标函数中,得到目标函数JNWFCM,对目标函数JNWFCM进行优化得到使目标函数JNWFCM为局部最小值的每一个像素的模糊隶属度值,利用每一个像素的模糊隶属度值对图像进行分割;
具体包括以下步骤:
步骤一,设定聚类数目c、模糊加权指数m和迭代停止阈值εL;
步骤二,采用随机化的方法初始化聚类中心vk;
步骤三,计算每一个像素的邻域像素权重Wi;
步骤四,将步骤三计算得到的权重Wi={wir,r∈Ni}融入到FCM算法的目标函数中,得到目标函数JNWFCM为
式(5)中,c为聚类数目,uik为第i个像素属于第k类的模糊隶属度值,m为模糊加权指数,控制聚类结果模糊程度,vk为第k个类的中心,Ni为第i个像素的邻域,xr为像素灰度值;
步骤五,对步骤四得到的目标函数JNWFCM进行优化得到计算像素的模糊隶属度值的式(6)和计算聚类中心的式(7);
步骤六,对步骤五中的式(6)和式(7)进行迭代,直至式(8)成立,则得到使目标函数JNWFCM为局部最小值的每一个像素的模糊隶属度值uik,否则继续迭代;
||vold-vnew||≤εL (8)
式(8)中,vnew为当前迭代的聚类中心,vold为前一步迭代的聚类中心;
步骤七,根据步骤六得到的uik,代入式(9)计算每个像素的分类标签,完成图像分割;
所述步骤三中每一个像素的邻域像素权重Wi计算步骤如下:
首先,设X={x1,x2,...,xn}表示具有n个像素的图像,xi为第i个像素的灰度值,像素xi的邻域像素的灰度均值μi为
式(1)中,Ni为像素xi的邻域,NR为集合Ni的势,xr为像素xi的邻域像素,xr和μi之间的L1范数为
σir=|xr-μi| (2)
定义像素xi的邻域像素xr(r∈Ni)的权重为
式(3)中β为伸缩因子,这里取β=3;σi为σir(r∈Ni)的最大值,即对ξir进行归一化处理得权重wir:
记像素xi的邻域像素的权重为Wi={wir,r∈Ni}。
2.根据权利要求1所述的一种邻域加权模糊聚类图像分割方法,其特征在于,所述步骤五中采用拉格朗日乘子法进行优化计算。
3.根据权利要求1所述的一种邻域加权模糊聚类图像分割方法,其特征在于,所述步骤三中的邻域为二阶邻域,在二阶邻域中,每个像素有8个邻点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010197180.2A CN111402272B (zh) | 2020-03-19 | 2020-03-19 | 一种邻域加权模糊聚类图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010197180.2A CN111402272B (zh) | 2020-03-19 | 2020-03-19 | 一种邻域加权模糊聚类图像分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111402272A CN111402272A (zh) | 2020-07-10 |
CN111402272B true CN111402272B (zh) | 2023-09-29 |
Family
ID=71431033
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010197180.2A Active CN111402272B (zh) | 2020-03-19 | 2020-03-19 | 一种邻域加权模糊聚类图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111402272B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000073972A1 (en) * | 1999-05-28 | 2000-12-07 | Texas Tech University | Adaptive vector quantization/quantizer |
WO2011084142A1 (en) * | 2010-01-11 | 2011-07-14 | Tandent Vision Science, Inc. | Image segregation system with method for handling textures |
CN108062757A (zh) * | 2018-01-05 | 2018-05-22 | 北京航空航天大学 | 一种利用改进直觉模糊聚类算法提取红外目标的方法 |
CN108198193A (zh) * | 2018-01-16 | 2018-06-22 | 北京航空航天大学 | 一种利用改进直觉模糊聚类算法分割红外舰船图像的方法 |
CN109064470A (zh) * | 2018-08-28 | 2018-12-21 | 河南工业大学 | 一种基于自适应模糊聚类的图像分割方法及装置 |
CN109389608A (zh) * | 2018-10-19 | 2019-02-26 | 山东大学 | 以平面为聚类中心具有抗噪性的模糊聚类图像分割方法 |
CN110211126A (zh) * | 2019-06-12 | 2019-09-06 | 西安邮电大学 | 基于直觉模糊c均值聚类的图像分割方法 |
-
2020
- 2020-03-19 CN CN202010197180.2A patent/CN111402272B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000073972A1 (en) * | 1999-05-28 | 2000-12-07 | Texas Tech University | Adaptive vector quantization/quantizer |
WO2011084142A1 (en) * | 2010-01-11 | 2011-07-14 | Tandent Vision Science, Inc. | Image segregation system with method for handling textures |
CN108062757A (zh) * | 2018-01-05 | 2018-05-22 | 北京航空航天大学 | 一种利用改进直觉模糊聚类算法提取红外目标的方法 |
CN108198193A (zh) * | 2018-01-16 | 2018-06-22 | 北京航空航天大学 | 一种利用改进直觉模糊聚类算法分割红外舰船图像的方法 |
CN109064470A (zh) * | 2018-08-28 | 2018-12-21 | 河南工业大学 | 一种基于自适应模糊聚类的图像分割方法及装置 |
CN109389608A (zh) * | 2018-10-19 | 2019-02-26 | 山东大学 | 以平面为聚类中心具有抗噪性的模糊聚类图像分割方法 |
CN110211126A (zh) * | 2019-06-12 | 2019-09-06 | 西安邮电大学 | 基于直觉模糊c均值聚类的图像分割方法 |
Non-Patent Citations (3)
Title |
---|
基于分层的直觉模糊C均值聚类图像分割算法;池桂英等;《计算机工程与设计》(第12期);全文 * |
改进的基于邻域隶属度约束的FCM图像分割算法;仲崇峰等;《吉林大学学报(信息科学版)》(第06期);全文 * |
空间加权模糊C均值聚类图像分割算法;李小和等;《西安石油大学学报(自然科学版)》(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111402272A (zh) | 2020-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103886589B (zh) | 面向目标的自动化高精度边缘提取方法 | |
CN109448015B (zh) | 基于显著图融合的图像协同分割方法 | |
CN108198193B (zh) | 一种利用改进直觉模糊聚类算法分割红外舰船图像的方法 | |
CN109584282B (zh) | 一种基于sift特征与光流模型的非刚性图像配准方法 | |
CN103761726B (zh) | 基于fcm的分块自适应图像分割方法 | |
CN107742297B (zh) | 一种三维ct图像的局部三维最大类间方差分割方法 | |
Kuo et al. | Application of computer vision in the automatic identification and classification of woven fabric weave patterns | |
CN111709901B (zh) | 基于FCM聚类匹配+Wallis滤波的无重多/高光谱遥感图像匀色方法 | |
CN111553837A (zh) | 一种基于神经风格迁移的艺术文本图像生成方法 | |
CN110796667A (zh) | 基于改进小波聚类的彩色图像分割方法 | |
CN106600610B (zh) | 一种fcm图像分割方法及装置 | |
Zhao et al. | A novel Neutrosophic image segmentation based on improved fuzzy C-means algorithm (NIS-IFCM) | |
Lorette et al. | Fully unsupervised fuzzy clustering with entropy criterion | |
CN108846845B (zh) | 基于缩略图与分层模糊聚类的sar图像分割方法 | |
CN111402272B (zh) | 一种邻域加权模糊聚类图像分割方法 | |
Afshar et al. | Lung tumor area recognition in CT images based on Gustafson-Kessel clustering | |
Bahadure et al. | Performance analysis of image segmentation using watershed algorithm, fuzzy c-means of clustering algorithm and simulink design | |
CN110910392B (zh) | 一种基于多重空间信息加权的模糊聚类的图像分割方法 | |
CN106504260B (zh) | 一种fcm图像分割方法和系统 | |
CN113538445A (zh) | 一种基于加权鲁棒fcm聚类的图像分割方法及系统 | |
CN111079750B (zh) | 一种基于局部区域聚类的电力设备故障区域提取方法 | |
Li et al. | Adaptive spatially weighted fuzzy c-means clustering for image segmentation | |
Qin et al. | Pavement image segmentation based on fast FCM clustering with spatial information in internet of things | |
Liu et al. | Automatic Lung Parenchyma Segmentation of CT Images Based on Matrix Grey Incidence. | |
Hazarika et al. | A hierarchical spatial fuzzy C means algorithm for mammographic mass segmentation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |