第13卷第4期 智能系统学报 Vol.13 No.4 2018年8月 CAAI Transactions on Intelligent Systems Aug.2018 D0:10.11992/tis.201701008 网络出版地址：http:/kns.cnki.net/kcms/detail/23.1538.TP.20170702.0425.022.html 自适应灰度加权的鲁棒模糊C均值图像分割 陆海青'，葛洪伟2 (1.江南大学物联网工程学院，江苏无锡214122,2.江南大学轻工过程先进控制教育部重点实验室，江苏无 锡214122) 摘要：针对传统模糊C均值(fuzzy C-meas,FCM算法以及结合空间信息的相关改进算法分割精度较低、对 噪声敏感的问题，提出一种自适应灰度加权的鲁棒模糊C均值图像分割算法。首先，通过定义像素间的局部灰 度相似性测度来反映各像素对局部邻域的影响程度，并根据邻域窗口中各像素的灰度差异，利用指数函数进一 步控制邻域像素的影响权重，实现像素灰度的自适应加权，从而提高像素灰度计算的准确性。其次，构造出一 种改进的距离测度代替传统的欧氏距离，用于计算各像素与聚类中心之间的相似距离，增强算法对噪声和异常 值的鲁棒性。最后，将提出的自适应灰度加权方法与改进的距离测度应用到FCM算法中，实现图像分割。实 验结果表明.该算法需根据图像噪声的强度适当地选取邻域窗口大小，在此条件下算法能够取得较优的分割效 果和运行效率，且对噪声具有较强的鲁棒性。 关键词：模糊C均值；图像分割；自适应灰度加权；空间信息；相似距离；抗噪性 中图分类号：TP391.4文献标志码：A文章编号：1673-4785(2018)04-0584-10 中文引用格式：陆海青，葛洪伟.自适应灰度加权的鲁棒模糊C均值图像分割.智能系统学报，2018,13(4)：584593。 英文引用格式：LU Haiqing,GE Hongwei.Adaptive gray-.weighted robust fuzzy C-means algorithm for image segmentation[J].. CAAI transactions on intelligent systems,2018,13(4):584-593. Adaptive gray-weighted robust fuzzy C-means algorithm for image segmentation LU Haiqing',GE Hongwei2 (1.School of Internet of Things,Jiangnan University,Wuxi214122,China;2.Ministry of Education Key Laboratory of Advanced Process Control for Light Industry,Jiangnan University,Wuxi214122,China) Abstract:The traditional fuzzy C-means(FCM)algorithm and its corresponding improved algorithm that is combined with spatial information have low segmentation accuracy and poor robustness to noise.To address these defects,we pro- pose a robust FCM image segmentation algorithm based on adaptive gray-weighting.First,we define a local grayscale similarity measure for pixels to reflect the influence of all pixels on the local neighborhood.Regarding the grayscale dif- ference between pixels in a neighborhood window,we utilize an exponential function to further control the influence weight of a neighborhood pixel and realize adaptive weighting of the pixel grayscale to improve its calculation accuracy Next,to strengthen the robustness of the algorithm to noise and outliers,we use an improved distance measure to re- place the traditional Euclidean distance and use it to calculate the similarity distance between the pixels and the cluster- ing center.Finally,we apply this new method based on adaptive gray weight and enhanced distance measurement to an FCM algorithm for image segmentation.Our experimental results show that,for the algorithm,the size of the neighbor- hood window must be properly selected on basis of the noise intensity of an image.Under this condition,an excellent segmentation effect and operational efficiency can be achieved,in addition to excellent robustness to noise. Keywords:fuzzy C-means;image segmentation;adaptive gray weight,spatial information;similarity distance;noise resistance 收稿日期：2017-01-11.网络出版日期：2017-07-02. 基金项目：江苏省普通高校研究生科研创新计划项目KYLX160781 图像分割是图像处理中十分重要的研究内 KYLX160782):江苏高校优势学科建设工程资助项 容，图像分割的质量将对后续的特征提取、图像 目(PAPD) 通信作者：葛洪伟.E-mail:ghw8601@163.com. 识别等工作产生直接的影响。模糊C均值(FCM)
DOI: 10.11992/tis.201701008 网络出版地址: http://kns.cnki.net/kcms/detail/23.1538.TP.20170702.0425.022.html 自适应灰度加权的鲁棒模糊 C 均值图像分割 陆海青1 ，葛洪伟1,2 （1. 江南大学 物联网工程学院，江苏 无锡 214122; 2. 江南大学 轻工过程先进控制教育部重点实验室，江苏 无 锡 214122） 摘 要：针对传统模糊 C 均值 (fuzzy C-means，FCM) 算法以及结合空间信息的相关改进算法分割精度较低、对 噪声敏感的问题，提出一种自适应灰度加权的鲁棒模糊 C 均值图像分割算法。首先，通过定义像素间的局部灰 度相似性测度来反映各像素对局部邻域的影响程度，并根据邻域窗口中各像素的灰度差异，利用指数函数进一 步控制邻域像素的影响权重，实现像素灰度的自适应加权，从而提高像素灰度计算的准确性。其次，构造出一 种改进的距离测度代替传统的欧氏距离，用于计算各像素与聚类中心之间的相似距离，增强算法对噪声和异常 值的鲁棒性。最后，将提出的自适应灰度加权方法与改进的距离测度应用到 FCM 算法中，实现图像分割。实 验结果表明，该算法需根据图像噪声的强度适当地选取邻域窗口大小，在此条件下算法能够取得较优的分割效 果和运行效率，且对噪声具有较强的鲁棒性。 关键词：模糊 C 均值；图像分割；自适应灰度加权；空间信息；相似距离；抗噪性 中图分类号：TP391.4 文献标志码：A 文章编号：1673−4785(2018)04−0584−10 中文引用格式：陆海青, 葛洪伟. 自适应灰度加权的鲁棒模糊 C 均值图像分割[J]. 智能系统学报, 2018, 13(4): 584–593. 英文引用格式：LU Haiqing, GE Hongwei. Adaptive gray-weighted robust fuzzy C-means algorithm for image segmentation[J]. CAAI transactions on intelligent systems, 2018, 13(4): 584–593. Adaptive gray-weighted robust fuzzy C-means algorithm for image segmentation LU Haiqing1 ，GE Hongwei1,2 (1. School of Internet of Things, Jiangnan University, Wuxi 214122, China; 2. Ministry of Education Key Laboratory of Advanced Process Control for Light Industry, Jiangnan University, Wuxi 214122, China) Abstract: The traditional fuzzy C-means (FCM) algorithm and its corresponding improved algorithm that is combined with spatial information have low segmentation accuracy and poor robustness to noise. To address these defects, we propose a robust FCM image segmentation algorithm based on adaptive gray-weighting. First, we define a local grayscale similarity measure for pixels to reflect the influence of all pixels on the local neighborhood. Regarding the grayscale difference between pixels in a neighborhood window, we utilize an exponential function to further control the influence weight of a neighborhood pixel and realize adaptive weighting of the pixel grayscale to improve its calculation accuracy Next, to strengthen the robustness of the algorithm to noise and outliers, we use an improved distance measure to replace the traditional Euclidean distance and use it to calculate the similarity distance between the pixels and the clustering center. Finally, we apply this new method based on adaptive gray weight and enhanced distance measurement to an FCM algorithm for image segmentation. Our experimental results show that, for the algorithm, the size of the neighborhood window must be properly selected on basis of the noise intensity of an image. Under this condition, an excellent segmentation effect and operational efficiency can be achieved, in addition to excellent robustness to noise. Keywords: fuzzy C-means; image segmentation; adaptive gray weight; spatial information; similarity distance; noise resistance 图像分割是图像处理中十分重要的研究内 容，图像分割的质量将对后续的特征提取、图像 识别等工作产生直接的影响。模糊 C 均值 (FCM) 收稿日期：2017−01−11. 网络出版日期：2017−07−02. 基金项目：江苏省普通高校研究生科研创新计划项目(KYLX16_0781， KYLX16_0782)；江苏高校优势学科建设工程资助项 目 (PAPD). 通信作者：葛洪伟. E-mail：ghw8601@163.com. 第 13 卷第 4 期 智 能 系 统 学 报 Vol.13 No.4 2018 年 8 月 CAAI Transactions on Intelligent Systems Aug. 2018
第4期 陆海青，等：自适应灰度加权的鲁棒模糊C均值图像分割 ·585· 算法是模糊聚类中最经典的方法，目前已广泛 到FCM算法中，在保证计算速度的条件下增强了 应用于医学、遥感影像分割2等方面。考虑到实 算法对噪声的鲁棒性。但该方法引入了两个额外 际应用中图像像素的归类具有不确定性，FCM算 参数，需进行人为设置，且参数选取的好坏将直 法利用模糊隶属度来衡量像素归属于某一类的程 接影响算法的分割效果。文献[16]提出了一种基 度，克服了传统硬聚类方法将像素归类一刀切的 于灰度与空间特征的FCM图像分割算法，从修改 缺陷。然而，FCM算法在应用于图像分割时仍存 像素空间隶属度的角度出发，利用像素的灰度相 在一些缺陷：计算复杂度高、对聚类初值和噪声 似度与空间邻域分布信息构造出新的空间隶属函 敏感、未考虑空间邻域信息。 数，提升了算法对噪声的抑制能力，具有良好的 为弥补这些不足，国内外诸多学者进行了大 聚类性能。除此之外，近年来涌现出一些结合图 量的研究工作以改进FCM算法的性能。针对 像局部熵.18、非局部均值去噪1920、中智学12 FCM算法计算效率较低的问题，现有改进方法主 等理论的先进方法。这些改进方法尽管在分割性 要从减少迭代次数和降低数据集维度两方面进行 能和抗噪性能方面取得了较大幅度的提高，但都 加速，如基于图像灰度级的FCM算法门，利用图 具有较高的计算复杂度，实时性较差。 像灰度级远小于像素个数的特性，采用一种线性 针对以上问题，本文采取一种自适应权重分 加权和图像代替原图像进行迭代计算，大大减少 配的策略，合理地为邻域窗口中的各像素赋予相 了参与计算的像素个数，进而提高了计算速度，但 应的权重；同时采用一种改进的距离测度代替原 算法的分割精度有所下降。针对FCM算法对初 始的欧氏距离，进而提出一种自适应灰度加权的 始聚类中心敏感的问题，Sikka等充分利用直方 鲁棒模糊C-均值算法(adaptive gray-weighted based 图的统计特性，通过检测直方图中各灰度区间内 robust fuzzy C-means algorithm,AGWRFCM)。对 的局部峰点来确定初始聚类中心和聚类数，并对 不同图像的分割结果表明，该算法能够在不同强 脑肿瘤图像进行分割，取得了比较准确的分割结 度的噪声条件下获得良好的分割性能。 果。鉴于智能优化算法具有良好的全局搜索能力， 因此也有学者利用智能优化算法对FCM进行初 1模糊C均值算法 始化0，在加速迭代收敛的同时避免了算法陷入 局部最优解。针对FCM算法中欧氏距离对噪声 考虑由n个p维样本组成的数据集X={x1,x2,…, 和异常点敏感的问题，国内外学者采用不同的距 x}∈RP,FCM算法旨在最小化目标函数Jm以完 离测度对欧氏距离进行改进，如核诱导距离， 成对样本数据的模糊划分，即 马氏距离等，以提升算法对噪声的鲁棒性。然 J(U.V) :-vall (1) 而距离测度中各参数（如核函数带宽等）的选取 k=1 i=1 对聚类分割性能具有较大的影响，且计算复杂度 式中：U=eR为隶属度矩阵，满足u=l 随图像规模的增大而增加。 针对FCM算法未考虑空间邻域信息的缺陷， i=1,2,…,n,V={w,2,…,}为聚类中心集合， 许多学者提出了结合空间信息的改进算法。Ahmed c∈[2，n为聚类数目，m∈[1，+o)为模糊指数，常取 等提出了一种基于空间邻域约束信息的FCM m=2。利用Lagrange乘子法对Jm进行迭代更新， 以实现目标函数的最小化，即 (fuzzy C-means with spatial constraints,FCM_S) 法，通过在FCM算法的目标函数中引入一空间惩 c,-v厂2m-) = (2) 罚项来调节邻域像素对中心像素的影响，提升了 FCM算法的抗噪性能，但每次迭代中均需计算图 像中各像素的邻域灰度，计算复杂度较高。为解 决该问题，Chen等对其进行了改进，通过融合 Vk= (3) 邻域像素的灰度均值与中值，提出了一种基于邻 域均值和中值的空间FCM(FCMS1,FCMS2)算 法，其分别对含有高斯噪声和椒盐噪声的图像具 反复迭代式(2(3)，直至FCM算法收敛。 有较好的分割效果。随后，Cai等提出了一种融 2 自适应灰度加权的鲁棒模糊C均 合局部空间与灰度差的快速广义FCM(fast gener- 值图像分割 alized fuzzy C-means,.FGFCM)算法，将各邻域窗 口中像素的局部空间距离与灰度差信息同时引入 现有的结合空间信息的FCM改进算法31主
算法[1]是模糊聚类中最经典的方法，目前已广泛 应用于医学、遥感影像分割[2-6]等方面。考虑到实 际应用中图像像素的归类具有不确定性，FCM 算 法利用模糊隶属度来衡量像素归属于某一类的程 度，克服了传统硬聚类方法将像素归类一刀切的 缺陷。然而，FCM 算法在应用于图像分割时仍存 在一些缺陷：计算复杂度高、对聚类初值和噪声 敏感、未考虑空间邻域信息。 为弥补这些不足，国内外诸多学者进行了大 量的研究工作以改进 FCM 算法的性能。针对 FCM 算法计算效率较低的问题，现有改进方法主 要从减少迭代次数和降低数据集维度两方面进行 加速，如基于图像灰度级的 FCM 算法[7] ，利用图 像灰度级远小于像素个数的特性，采用一种线性 加权和图像代替原图像进行迭代计算，大大减少 了参与计算的像素个数，进而提高了计算速度，但 算法的分割精度有所下降。针对 FCM 算法对初 始聚类中心敏感的问题，Sikka 等 [8]充分利用直方 图的统计特性，通过检测直方图中各灰度区间内 的局部峰点来确定初始聚类中心和聚类数，并对 脑肿瘤图像进行分割，取得了比较准确的分割结 果。鉴于智能优化算法具有良好的全局搜索能力， 因此也有学者利用智能优化算法对 FCM 进行初 始化[9-10] ，在加速迭代收敛的同时避免了算法陷入 局部最优解。针对 FCM 算法中欧氏距离对噪声 和异常点敏感的问题，国内外学者采用不同的距 离测度对欧氏距离进行改进，如核诱导距离[5, 11] 、 马氏距离[12]等，以提升算法对噪声的鲁棒性。然 而距离测度中各参数 (如核函数带宽等) 的选取 对聚类分割性能具有较大的影响，且计算复杂度 随图像规模的增大而增加。 针对 FCM 算法未考虑空间邻域信息的缺陷， 许多学者提出了结合空间信息的改进算法。Ahmed 等 [13]提出了一种基于空间邻域约束信息的 FCM (fuzzy C-means with spatial constraints，FCM_S) 算 法，通过在 FCM 算法的目标函数中引入一空间惩 罚项来调节邻域像素对中心像素的影响，提升了 FCM 算法的抗噪性能，但每次迭代中均需计算图 像中各像素的邻域灰度，计算复杂度较高。为解 决该问题，Chen 等 [14]对其进行了改进，通过融合 邻域像素的灰度均值与中值，提出了一种基于邻 域均值和中值的空间 FCM(FCM_S1，FCM_S2 ) 算 法，其分别对含有高斯噪声和椒盐噪声的图像具 有较好的分割效果。随后，Cai等 [15]提出了一种融 合局部空间与灰度差的快速广义 FCM(fast generalized fuzzy C-means，FGFCM) 算法，将各邻域窗 口中像素的局部空间距离与灰度差信息同时引入 到 FCM 算法中，在保证计算速度的条件下增强了 算法对噪声的鲁棒性。但该方法引入了两个额外 参数，需进行人为设置，且参数选取的好坏将直 接影响算法的分割效果。文献[16]提出了一种基 于灰度与空间特征的 FCM 图像分割算法，从修改 像素空间隶属度的角度出发，利用像素的灰度相 似度与空间邻域分布信息构造出新的空间隶属函 数，提升了算法对噪声的抑制能力，具有良好的 聚类性能。除此之外，近年来涌现出一些结合图 像局部熵[17-18] 、非局部均值去噪[19-20] 、中智学[21-22] 等理论的先进方法。这些改进方法尽管在分割性 能和抗噪性能方面取得了较大幅度的提高，但都 具有较高的计算复杂度，实时性较差。 针对以上问题，本文采取一种自适应权重分 配的策略，合理地为邻域窗口中的各像素赋予相 应的权重；同时采用一种改进的距离测度代替原 始的欧氏距离，进而提出一种自适应灰度加权的 鲁棒模糊 C-均值算法 (adaptive gray-weighted based robust fuzzy C-means algorithm，AGWRFCM)。对 不同图像的分割结果表明，该算法能够在不同强 度的噪声条件下获得良好的分割性能。 1 模糊 C 均值算法 X = {x1, x2,··· , xn} ∈ R n×p Jm 考虑由 n 个 p 维样本组成的数据集 ，FCM 算法旨在最小化目标函数 以完 成对样本数据的模糊划分，即 Jm(U,V) = ∑c k=1 ∑n i=1 u m ki∥xi −vk∥ 2 (1) U = {uki} ∈ R n×c ∑c k=1 uki = 1, ∀i = 1,2,··· ,n V = {v1, v2,··· , vc} c ∈ [2,n] m ∈ [1,+∞) m = 2 Jm 式中： 为隶属度矩阵，满足 ， 为聚类中心集合， 为聚类数目， 为模糊指数，常取 。利用 Lagrange 乘子法对 进行迭代更新， 以实现目标函数的最小化，即 uki = ∥xi −vk∥ −2/(m−1) ∑c j=1 xi −vj −2/(m−1) (2) vk = ∑n i=1 u m kixi ∑n i=1 u m ki (3) 反复迭代式 (2)~(3)，直至 FCM 算法收敛。 2 自适应灰度加权的鲁棒模糊 C 均 值图像分割 现有的结合空间信息的 FCM 改进算法[13-15]主 第 4 期 陆海青，等：自适应灰度加权的鲁棒模糊 C 均值图像分割 ·585·
·586· 智能系统学报 第13卷 要存在以下缺陷：1)这些算法并没有充分利用图 其次，对当前邻域中各像素的C取平均值，作 像中像素分布的特点，在引入空间信息时未能充 为该邻域的灰度分布度量，即 分考虑到邻域像素对中心像素灰度贡献的差异， 对像素灰度的计算不够准确：2)这些算法均引人 c (6) 了额外的空间参数，用于控制邻域像素对中心像 在此基础上定义邻域像素的影响因子y,将 素的约束程度，而空间参数需根据经验或大量实 各邻域点的C,与当前邻域进行比较，即 验试错来进行确定，这样会极大地增加算法的处 Yij=CJ-C.jEN (7) 理时间，并且若参数选取不佳将会直接影响算法 由式(7)可以看出，对于边缘点或噪声点， 的分割效果1：3)这些算法在计算像素与聚类中 y,的值较大，表明该点对当前邻域的影响较大；而 心之间的相似距离时均使用欧氏距离测度，难以 对于非边缘点与非噪声点，y的值较小，表明该点 准确地反映像素与聚类之间的关系，且对噪声和 对邻域的影响较小。因此y决定了局部邻域内各 异常点较敏感。 像素灰度的影响程度。 针对以上缺陷，本文主要从两方面对FCM算 在计算中心像素的灰度值时，对于y值较大 法进行改进：首先，为充分考虑邻域像素对中心 的点，应减小其对中心像素灰度计算的影响，反 像素影响的差异性，采用一种自适应灰度加权的 之应增大其影响程度。鉴于指数函数具有较快的 方法，根据局部窗口中各像素的灰度分布自适应 衰减速度，故采用指数函数进一步控制邻域像素 地调节邻域像素对中心像素的权重，使其在计算 的影响权重，具体形式为 中心像素的灰度时能够更充分地利用局部信息， wij exp(-Yij)=exp(-(Cj-C)) (8) 以提升像素灰度的计算精度；其次，为更好地反 式中：w为邻域像素相对于中心像素的权重。 映图像中各像素与聚类中心之间的相似性，用一 由式(8)可知，对于边缘点或噪声点，其对中 种改进的距离测度取代传统的欧氏距离测度，以 心像素的影响较大，w,的值更接近O,且随着y的 进一步增强算法对噪声的鲁棒性。 增加，权重下降幅度更大，从而进一步减小该邻 2.1自适应灰度加权 域点的影响程度；而对于非边缘点与非噪声点，w 对于一幅数字图像而言，其像素的灰度分布 能够取得较大的值，从而增大其对中心像素的影 通常具有以下特点：平坦区域内部的灰度分布较 响。因此w能够根据邻域像素的灰度差异自适应 为均匀，区域边缘或含噪区域的灰度差异较大， 地为各像素赋予合适的权重，进而达到提高中心 这在医学图像中表现尤为明显。由于人体组织结 像素灰度计算准确性的目的。 构的复杂性以及受医学成像中部分容积效应等因 最后，利用各邻域点的权重对邻域灰度作线 素的影响，使得获取的图像往往呈现出含噪较 性加权，将其作为该邻域中心像素的灰度值；同 多、对比度低、灰度分布不均匀、目标边界不连续 时为统一数据的量纲，将权重作归一化处理，得 等特性1。因此，在引入像素的空间信息时需考 到最终的待聚类图像，即 虑局部邻域内各像素对中心像素影响程度的差异 ∑wx 性，以进一步提高图像分割的准确性。 (9) 首先，对于图像中的任一邻域窗口N,定义局 部灰度相似性测度C,将当前邻域内的各像素与 jeN, 该邻域内所有像素的平均灰度进行比较，即 2.2距离测度的改进 C=xj-jEN 传统的欧氏距离无法解决算法对噪声敏感的 (4) 问题1。核诱导距离的实质是将像素映射到高 (5) 维特征空间中进行处理，以提升对高维数据的聚 式中：x为邻域像素，为邻域像素的平均灰度， 类性能。虽一定程度上弥补了欧氏距离的不足， N,表示以像素为中心的邻域，n,表示以像素为中 但其难以克服噪声对聚类性能的影响，对噪声的 心的邻域窗口大小，表示任一范数。 抑制能力仍不够强，因此核诱导距离也无法从根 由式(4)可以看出，对于灰度变化平缓的区 本上解决对噪声敏感的问题。 域，C,的值较小；而对于灰度发生突变的区域（如 为弥补这一不足，本文采用一种改进的距离 边缘、噪声等)，C,的值较大。因此C,反映了局部 测度，具体形式为 邻域内像素灰度分布的均匀程度。 dR(x,)=V1-R(x,) (10)
要存在以下缺陷：1) 这些算法并没有充分利用图 像中像素分布的特点，在引入空间信息时未能充 分考虑到邻域像素对中心像素灰度贡献的差异， 对像素灰度的计算不够准确；2) 这些算法均引入 了额外的空间参数，用于控制邻域像素对中心像 素的约束程度，而空间参数需根据经验或大量实 验试错来进行确定，这样会极大地增加算法的处 理时间，并且若参数选取不佳将会直接影响算法 的分割效果[19] ；3) 这些算法在计算像素与聚类中 心之间的相似距离时均使用欧氏距离测度，难以 准确地反映像素与聚类之间的关系，且对噪声和 异常点较敏感。 针对以上缺陷，本文主要从两方面对 FCM 算 法进行改进：首先，为充分考虑邻域像素对中心 像素影响的差异性，采用一种自适应灰度加权的 方法，根据局部窗口中各像素的灰度分布自适应 地调节邻域像素对中心像素的权重，使其在计算 中心像素的灰度时能够更充分地利用局部信息， 以提升像素灰度的计算精度；其次，为更好地反 映图像中各像素与聚类中心之间的相似性，用一 种改进的距离测度取代传统的欧氏距离测度，以 进一步增强算法对噪声的鲁棒性。 2.1 自适应灰度加权 对于一幅数字图像而言，其像素的灰度分布 通常具有以下特点：平坦区域内部的灰度分布较 为均匀，区域边缘或含噪区域的灰度差异较大， 这在医学图像中表现尤为明显。由于人体组织结 构的复杂性以及受医学成像中部分容积效应等因 素的影响，使得获取的图像往往呈现出含噪较 多、对比度低、灰度分布不均匀、目标边界不连续 等特性[18]。因此，在引入像素的空间信息时需考 虑局部邻域内各像素对中心像素影响程度的差异 性，以进一步提高图像分割的准确性。 Ni Cj 首先，对于图像中的任一邻域窗口 ，定义局 部灰度相似性测度 ，将当前邻域内的各像素与 该邻域内所有像素的平均灰度进行比较，即 Cj = xj − x¯ , j ∈ Ni (4) x¯ = 1 ni ∑ j∈Ni xj (5) xj x¯ Ni i ni i ∥·∥ 式中： 为邻域像素， 为邻域像素的平均灰度， 表示以像素 为中心的邻域， 表示以像素 为中 心的邻域窗口大小， 表示任一范数。 Cj Cj Cj 由式 (4) 可以看出，对于灰度变化平缓的区 域， 的值较小；而对于灰度发生突变的区域 (如 边缘、噪声等)， 的值较大。因此 反映了局部 邻域内像素灰度分布的均匀程度。 其次，对当前邻域中各像素的 Cj 取平均值，作 为该邻域的灰度分布度量，即 C¯= 1 ni ∑ j∈Ni Cj (6) γi j Cj 在此基础上定义邻域像素的影响因子 ，将 各邻域点的 与当前邻域进行比较，即 γi j = Cj −C¯, j ∈ Ni (7) γi j γi j γi j 由式 (7) 可以看出，对于边缘点或噪声点， 的值较大，表明该点对当前邻域的影响较大；而 对于非边缘点与非噪声点， 的值较小，表明该点 对邻域的影响较小。因此 决定了局部邻域内各 像素灰度的影响程度。 在计算中心像素的灰度值时，对于 γi j 值较大 的点，应减小其对中心像素灰度计算的影响，反 之应增大其影响程度。鉴于指数函数具有较快的 衰减速度，故采用指数函数进一步控制邻域像素 的影响权重，具体形式为 wi j = exp(−γi j) = exp(−(Cj −C¯)) (8) wi j 式中： 为邻域像素 j 相对于中心像素 i 的权重。 wi j γi jwi j wi j 由式 (8) 可知，对于边缘点或噪声点，其对中 心像素的影响较大， 的值更接近 0，且随着 的 增加，权重下降幅度更大，从而进一步减小该邻 域点的影响程度；而对于非边缘点与非噪声点， 能够取得较大的值，从而增大其对中心像素的影 响。因此 能够根据邻域像素的灰度差异自适应 地为各像素赋予合适的权重，进而达到提高中心 像素灰度计算准确性的目的。 x˜i 最后，利用各邻域点的权重对邻域灰度作线 性加权，将其作为该邻域中心像素的灰度值；同 时为统一数据的量纲，将权重作归一化处理，得 到最终的待聚类图像 ，即 x˜i = ∑ j∈Ni wi jxj ∑ j∈Ni wi j (9) 2.2 距离测度的改进 传统的欧氏距离无法解决算法对噪声敏感的 问题[23]。核诱导距离[11]的实质是将像素映射到高 维特征空间中进行处理，以提升对高维数据的聚 类性能。虽一定程度上弥补了欧氏距离的不足， 但其难以克服噪声对聚类性能的影响，对噪声的 抑制能力仍不够强，因此核诱导距离也无法从根 本上解决对噪声敏感的问题。 为弥补这一不足，本文采用一种改进的距离 测度，具体形式为 dR(xi , vk) = √ 1−R(xi , vk) (10) ·586· 智 能 系 统 学 报 第 13 卷
第4期 陆海青，等：自适应灰度加权的鲁棒模糊C均值图像分割 ·587· 式中： 行实验，将本文算法与FCMS,、FCMS,、En- R(xi.V)=exp(-Bllx:-velP) (11) FCMO、FGFCM'、FLICM(模糊局部信息C-均值， Zun-sr fuzzy local information C-means)2s1、文献 B=( [26](WFCM)几种算法进行对比，以测试算法的分 (12) n 割效果。实验中采用的测试图像如图1所示，几 (13) 种对比算法中的参数根据文献[7,14-15,25-26进 行设置。为获得良好的实验效果，本文实验中的 改进的距离测度建立在鲁棒统计理论的基础 相关参数设置如下：m=2,s=10-5,iter max=500, 上，因此对噪声或异常值具有较强的稳定性2。 邻域窗口大小ω=3（即3×3邻域），人工合成图像 该距离测度虽然在形式上与核诱导距离类似，但 中c=4,自然图像和医学图像中c=3。其中，迭代 其本质仍然是在原图像空间中进行处理，并未将 终止阈值ε为一较小数，其值通常根据人为经验 像素映射到高维特征空间中。下文实验中将验 进行选取。实验环境为MATLAB R2014a,Intel 证这一距离测度对噪声具有较强的鲁棒性，且均 (R)Core(TM)i3-4160 CPU @3.60 GHz,4.00 GB 优于欧式距离与核诱导距离。 内存。 2.3本文算法 根据以上分析，用dx,)代替式(1)中的 x-v,再用式(15)代替d(x,y)中的x,得到本 文算法的目标函数JAGWRFCM,即 (a)synthetic (b)synthetic, (c)eight 22a (14) 由Lagrange乘子法可得 [1-R(e,yg厂m-D lk 多n-E小 (15) (d)cameraman (e)MR, (f)MR, GRGE.Y) 图1实验中采用的测试图像 il Fig.1 Test images used in the experiment V&= (16) RE,n 3.1分割性能测试 =1 首先对人工合成图像(synthetic,244×244)进 AGWRFCM算法的具体步骤如下： 行分割实验（其中包含0、85、170、2554个灰度 1)设置聚类数目c,模糊加权指数m,迭代终 值)。用分割精度(segmentation accuracy,SA2作 止阈值e,最大迭代次数iter max,邻域窗口大小 为算法分割性能的客观评价指标，其定义为 w,当前迭代次数t=1: 2)初始化隶属度矩阵Uo和聚类中心Vo: SA=Anc (17) 3)根据式(10)(15)计算新的灰度图像元： j=1 4)根据式22)计算新的聚类中心V; 式中：c为聚类数，A为由分割算法得到的属于第 5)根据式(21)计算新的隶属度矩阵U; 类的像素点集，C:为标准图像中属于第类的像素 6)根据式(20)计算新的目标函数值J; 点集，SA表示正确分类的像素数占图像像素总数 7)若J+)-J 的比重。 iter_max,.则终止迭代，输出隶属度矩阵U和聚类 对该合成图像添加均值为0，方差为25的高 中心V;否则返回4)，继续下一次迭代。 斯噪声，采用7种算法对其进行分割，结果如图2 3实验结果与分析 所示。通过对比可以看出，EnFCM算法的分割结 果图中包含的噪声点最多，这是因为EnFCM 为验证本文算法的有效性，分别采用两组人 算法中未引入像素的空间邻域信息，对噪声的抑 工合成图像(synthetic,和syntheticz)、自然图像 制能力不强，因此该算法在含噪条件下的分割效 (eight和cameraman)和医学图像(MR,和MR2)进 果欠佳；FGFCM算法在EnFCM算法的基础上引
式中： R(xi , vk) = exp(−β∥xi −vk∥ 2 ) (11) β= ( ∑n i=1 ∥xi − x¯∥ 2 n ) −1 (12) x¯ = 1 n ∑n i=1 xi (13) 改进的距离测度建立在鲁棒统计理论的基础 上，因此对噪声或异常值具有较强的稳定性[24]。 该距离测度虽然在形式上与核诱导距离类似，但 其本质仍然是在原图像空间中进行处理，并未将 像素映射到高维特征空间中[24]。下文实验中将验 证这一距离测度对噪声具有较强的鲁棒性，且均 优于欧式距离与核诱导距离。 2.3 本文算法 d 2 R (xi , vk) ∥xi −vk∥ 2 d 2 R (xi , vk) xi JAGWRFCM 根据以上分析，用 代 替 式 (1) 中 的 ，再用式 (15) 代替 中的 ，得到本 文算法的目标函数 ，即 JAGWRFCM = ∑c k=1 ∑n i=1 u m kid 2 R (x˜i , vk) (14) 由 Lagrange 乘子法可得 uki = [1−R(x˜i , vk)]−1/(m−1) ∑c j=1 [1−R(x˜i , vj)]−1/(m−1) (15) vk = ∑n i=1 u m kiR(x˜i , vk) ˜xi ∑n i=1 u m kiR(x˜i , vk) (16) AGWRFCM 算法的具体步骤如下： c m ε ω t = 1 1) 设置聚类数目 ，模糊加权指数 ，迭代终 止阈值 ，最大迭代次数 iter_max，邻域窗口大小 ，当前迭代次数 ； U (0) V (0) 2) 初始化隶属度矩阵 和聚类中心 ； 3) 根据式 (10)~(15) 计算新的灰度图像 x˜i； 4) 根据式 (22) 计算新的聚类中心 V ； 5) 根据式 (21) 计算新的隶属度矩阵 U ； 6) 根据式 (20) 计算新的目标函数值 J ； J (t+1) − J (t) U V 7) 若 或当前迭代次数 iter_max，则终止迭代，输出隶属度矩阵 和聚类 中心 ；否则返回 4)，继续下一次迭代。 3 实验结果与分析 为验证本文算法的有效性，分别采用两组人 工合成图像 (synthetic1 和 synthetic2 )、自然图像 (eight 和 cameraman) 和医学图像 (MR1 和 MR2 ) 进 m = 2 ε = 10−5 iter_max = 500 ω = 3 c = 4 c = 3 ε 行实验，将本文算法与 FCM_S1 [14] 、FCM_S2 [14] 、EnFCM[7] 、FGFCM[15] 、FLICM(模糊局部信息 C-均值， fuzzy local information C-means)[ 2 5 ] 、文献 [26](WFCM) 几种算法进行对比，以测试算法的分 割效果。实验中采用的测试图像如图 1 所示，几 种对比算法中的参数根据文献[7, 14-15, 25-26]进 行设置。为获得良好的实验效果，本文实验中的 相关参数设置如下： ， ， ， 邻域窗口大小 (即 3×3邻域)，人工合成图像 中 ，自然图像和医学图像中 。其中，迭代 终止阈值 为一较小数，其值通常根据人为经验 进行选取。实验环境为 MATLAB R2014a，Intel (R) Core(TM) i3-4160 CPU @3.60 GHz，4.00 GB 内存。 3.1 分割性能测试 首先对人工合成图像 (synthetic1，244×244) 进 行分割实验 (其中包含 0、85、170、2 554 个灰度 值)。用分割精度 (segmentation accuracy，SA[25] ) 作 为算法分割性能的客观评价指标，其定义为 SA = ∑c i=1 Ai ∩Ci ∑c j=1 Cj (17) c Ai i Ci i SA 式中： 为聚类数， 为由分割算法得到的属于第 类的像素点集， 为标准图像中属于第 类的像素 点集， 表示正确分类的像素数占图像像素总数 的比重。 对该合成图像添加均值为 0，方差为 25 的高 斯噪声，采用 7 种算法对其进行分割，结果如图 2 所示。通过对比可以看出，EnFCM 算法的分割结 果图中包含的噪声点最多，这是因为 EnFCM 算法中未引入像素的空间邻域信息，对噪声的抑 制能力不强，因此该算法在含噪条件下的分割效 果欠佳；FGFCM 算法在 EnFCM 算法的基础上引 (e) MR1 (a) synthetic (c) eight 1 (f) MR2 (d) cameraman (b) synthetic2 图 1 实验中采用的测试图像 Fig. 1 Test images used in the experiment 第 4 期 陆海青，等：自适应灰度加权的鲁棒模糊 C 均值图像分割 ·587·
·588· 智能系统学报 第13卷 入了邻域像素的空间信息和灰度差，一定程度上 显下降，而AGWRFCM算法仍然能够取得较高的 提升了算法的抗噪性，因此其分割效果优于EFCM 分割精度。 算法，但图像中仍含有较多的噪声点，在含噪条 表17种算法对含噪合成图像的分割精度 件下的分割性能还不够高；FCMS、FCMS2算法 Table 1 Segmentation accuracy of seven algorithms on 结合了邻域像素的灰度均值和灰度中值，分割结 noisy synthetic image 果图中平坦区域较为平滑，但仍带有少量噪声， 分割算法 0=25 0-35 0=45 0=50 且图像边缘处存在大量的噪声点，因此分割结果 FCM_S1 0.9964 0.9939 0.9831 0.9769 不够理想；FLICM算法在目标函数中添加了包含 空间局部信息的模糊因子，更多地考虑了像素的 FCM_S2 0.9989 0.9955 0.9827 0.9768 空间信息，因此分割质量较高，视觉效果较好。 EnFCM 0.9954 0.9769 0.9379 0.9177 但其未能很好地保持边缘细节，导致图像中出现 FGFCM 0.9982 0.9925 0.9657 0.9483 了较多的伪边缘，产生了过平滑现象，因此该算 FLICM 0.9953 0.9955 0.9840 0.9788 法的分割效果仍不够理想；WFCM算法通过引入 WFCM 0.9972 0.9898 当前像素及其周围像素的灰度差来反映各像素对 0.97440.9653 分类的影响权重，其仍是对原图像进行分割处 AGWRFCM 0.99950.9991 0.98900.9847 理，并且目标函数中采用的是欧氏距离测度，对 其次，对自然图像(eight,.242×308)添加强度 噪声十分敏感，分割图像中仍包含大量噪声点， 为0.02的椒盐噪声，使用7种算法分别对其进行 因此WFCM算法难以有效地分割含噪图像；而 分割，结果如图3所示。 AGWRFCM算法是采用灰度加权后的新图像进 行分割，同时采用改进的距离测度进行计算，分 割结果十分接近原始图像，分割出的各区域连续 性强，边缘更清晰，很好地平衡了噪声抑制和边 缘保持之间的关系，不论是在视觉质量还是在分 (a)原图像 (b)噪声图像 (c)FCM S 割性能方面均优于其他几种算法。 (d)FCM S. (e)EnFCM (f)FGFCM (a)原图像 (b)噪声图像 (c)FCM S (g)FLICM (h)WFCM (i)AGWRFCM 图37种算法对含噪自然图像的分割结果 (d)FCM S, (e)EnFCM (①FGFCM Fig.3 Segmentation results of seven algorithms on noisy natural image 从图3中可以看出，FCMS1、FCM_S2、En FCM、FLICM以及WFCM算法的分割结果图中 均含有较多的噪声点，对噪声较为敏感，且 FLICM算法未能很好地将硬币中的纹理细节分 (g)FLICM (h)WFCM (①AGWRFCM 割出来，因此这些算法对硬币图像的分割效果均 图27种算法对含噪合成图像的分割结果 不佳：FGFCM算法虽能在一定程度上减弱噪声的 Fig.2 Segmentation results of seven algorithms on noisy 影响，但其对噪声的抑制能力有限：而AGWR synthetic image FCM算法可以很好地消除噪声的影响，且对硬币 表1中同时给出了各算法在噪声强度逐级递 中的纹理细节分割较为合理，这表明本文算法能 增（σ=25,35,45,50）时的SA值。从中可以看出，随 够获得较好的分割效果。 着噪声强度的增加，几种对比算法的分割精度明 为进一步评价本文算法的分割性能，采用划
入了邻域像素的空间信息和灰度差，一定程度上 提升了算法的抗噪性，因此其分割效果优于 EnFCM 算法，但图像中仍含有较多的噪声点，在含噪条 件下的分割性能还不够高；FCM_S1、FCM_S2 算法 结合了邻域像素的灰度均值和灰度中值，分割结 果图中平坦区域较为平滑，但仍带有少量噪声， 且图像边缘处存在大量的噪声点，因此分割结果 不够理想；FLICM 算法在目标函数中添加了包含 空间局部信息的模糊因子，更多地考虑了像素的 空间信息，因此分割质量较高，视觉效果较好。 但其未能很好地保持边缘细节，导致图像中出现 了较多的伪边缘，产生了过平滑现象，因此该算 法的分割效果仍不够理想；WFCM 算法通过引入 当前像素及其周围像素的灰度差来反映各像素对 分类的影响权重，其仍是对原图像进行分割处 理，并且目标函数中采用的是欧氏距离测度[26] ，对 噪声十分敏感，分割图像中仍包含大量噪声点， 因此 WFCM 算法难以有效地分割含噪图像；而 AGWRFCM 算法是采用灰度加权后的新图像进 行分割，同时采用改进的距离测度进行计算，分 割结果十分接近原始图像，分割出的各区域连续 性强，边缘更清晰，很好地平衡了噪声抑制和边 缘保持之间的关系，不论是在视觉质量还是在分 割性能方面均优于其他几种算法。 σ = 25,35,45,50 SA 表 1 中同时给出了各算法在噪声强度逐级递 增 ( ) 时的 值。从中可以看出，随 着噪声强度的增加，几种对比算法的分割精度明 显下降，而 AGWRFCM 算法仍然能够取得较高的 分割精度。 其次，对自然图像 (eight，242×308) 添加强度 为 0.02 的椒盐噪声，使用 7 种算法分别对其进行 分割，结果如图 3 所示。 从图 3 中可以看出，FCM_S1、FCM_S2、EnFCM、FLICM 以及 WFCM 算法的分割结果图中 均含有较多的噪声点，对噪声较为敏感， 且 FLICM 算法未能很好地将硬币中的纹理细节分 割出来，因此这些算法对硬币图像的分割效果均 不佳；FGFCM 算法虽能在一定程度上减弱噪声的 影响，但其对噪声的抑制能力有限；而 AGWRFCM 算法可以很好地消除噪声的影响，且对硬币 中的纹理细节分割较为合理，这表明本文算法能 够获得较好的分割效果。 为进一步评价本文算法的分割性能，采用划 表 1 7 种算法对含噪合成图像的分割精度 Table 1 Segmentation accuracy of seven algorithms on noisy synthetic image 分割算法 σ=25 σ=35 σ=45 σ=50 FCM_S1 0.996 4 0.993 9 0.983 1 0.976 9 FCM_S2 0.998 9 0.995 5 0.982 7 0.976 8 EnFCM 0.995 4 0.976 9 0.937 9 0.917 7 FGFCM 0.998 2 0.992 5 0.965 7 0.948 3 FLICM 0.995 3 0.995 5 0.984 0 0.978 8 WFCM 0.997 2 0.989 8 0.974 4 0.965 3 AGWRFCM 0.999 5 0.999 1 0.989 0 0.984 7 (a) 原图像 (b) 噪声图像 (c) FCM_S1 (d) FCM_S2 (e) EnFCM (f) FGFCM (g) FLICM (h) WFCM (i) AGWRFCM 图 2 7 种算法对含噪合成图像的分割结果 Fig. 2 Segmentation results of seven algorithms on noisy synthetic image (a) 原图像 (b) 噪声图像 (c) FCM_S1 (d) FCM_S2 (e) EnFCM (f) FGFCM (g) FLICM (h) WFCM (i) AGWRFCM 图 3 7 种算法对含噪自然图像的分割结果 Fig. 3 Segmentation results of seven algorithms on noisy natural image ·588· 智 能 系 统 学 报 第 13 卷
第4期 陆海青，等：自适应灰度加权的鲁棒模糊C均值图像分割 ·589· 分系数(partition coefficient,,Ve)Pm、划分嫡(parti- tion entropy,V)2以及聚类有效性指数(Xie-Beni index,V)3种评价指标对算法的聚类性能进行 定量评价，分别定义为 (a)原图像 (b)噪声图像 (c)FCM S. k=11 (18) (19) n 了x,- (d)FCM S, (e)EnFCM (f)FGFCM k=1e1 (20) n(minjatv,-ve) 由式(18)~(20)可知，V和V反映出划分矩 阵的模糊程度，V越大，计算出的划分矩阵模糊 性越小，像素的归类越明确，分割效果越好；同 理，V越小，算法的分割性能越好。而V则反映 (g)FLICM (h)WFCM (i)AGWRFCM 出各聚类之间的关联程度，V越小，同类中各像 图47种算法对含噪医学图像的分割结果 素之间的紧密性越高，不同类之间的关联性越 Fig.4 Segmentation results of seven algorithms on noisy 低，其聚类性能越好。 medical image 表2列出了7种算法对自然图像eight的划 表37种算法对含噪医学图像的聚类性能 分系数Vx、划分嫡V以及聚类有效性指数V。由 Table 3 Clustering performance of seven algorithms on 表2可以看出，与其余5种算法相比，AGWR- noisy medical image FCM算法具有更高的V值和更低的V值，这表明 分割算法 Vp Vpe V 由AGWRFCM算法计算得到的隶属度矩阵划分 FCM_S1 0.8873 0.0934 0.1041 性能更好，分割精度更高；同时AGWRFCM算法 FCM_S2 0.9477 0.0417 0.0895 具有更高的V值，表明类内像素的空间紧致性更 好，不同类像素的空间分离性更强，因此AGWR EnFCM 0.8668 0.1035 0.1547 FCM算法具有良好的聚类性能。 FGFCM 0.8750 0.1033 0.0804 表27种算法对含噪自然图像的聚类性能 FLICM 0.8333 0.1212 0.8129 Table 2 Clustering performance of seven algorithms on WFCM 0.9011 0.0956 0.1839 noisy natural image AGWRFCM 0.9504 0.0414 0.0136 分割算法 (平 Vpe Vih FCM_S 0.9074 0.0755 0.1110 由图4可以看出，在分割含噪医学MR图像 时，FCMS,、FCMS2、EnFCM、FGFCM、FLICM FCM_S2 0.9157 0.0713 0.1137 以及WFCM算法分割图中的背景噪声仍然较多， EnFCM 0.7875 0.1817 0.1017 抗噪性能较差，且FLICM算法分割出的脑组织内 FGFCM 0.7519 0.2064 0.0608 部结构细节信息丢失较严重，因此无法精确地将 FLICM 0.8912 0.0832 0.5548 MR脑组织分割出来；而AGWRFCM算法既能很 WFCM 0.8298 0.1425 0.2550 好地抑制噪声，又能保留更多的图像细节，且表3 AGWRFCM 0.9175 0.0688 0.0558 中的数据也充分反映出AGWRFCM算法具有更 优的聚类分割性能。因此，本文算法对医学图像 最后，对医学MR图像(MR1,256×256)2添加 也能够取得良好的分割效果。 强度为0.08的椒盐噪声，分割结果如图4所示， 3.2改进距离测度的抗噪性测试 各算法的划分系数V。c、划分熵V以及聚类有效性 由于本文中采用了一种改进的距离测度来计 指数V如表3所示。 算相似距离，因此为测试该距离测度对噪声的鲁
Vpc Vpe Vxb 分系数 (partition coefficient， ) [27] 、划分熵 (partition entropy， ) [27]以及聚类有效性指数 (Xie-Beni index， ) [28] 3 种评价指标对算法的聚类性能进行 定量评价，分别定义为 Vpc = ∑c k=1 ∑n i=1 u 2 ki n (18) Vpe = − ∑c k=1 ∑n i=1 uki loguki n (19) Vxb = ∑c k=1 ∑n i=1 u 2 ki∥xi −vk∥ 2 n(minj,k{ vj −vk 2 }) (20) Vpc Vpe Vpc Vpe Vxb Vxb 由式 (18)～(20) 可知， 和 反映出划分矩 阵的模糊程度， 越大，计算出的划分矩阵模糊 性越小，像素的归类越明确，分割效果越好；同 理， 越小，算法的分割性能越好。而 则反映 出各聚类之间的关联程度， 越小，同类中各像 素之间的紧密性越高，不同类之间的关联性越 低，其聚类性能越好。 Vpc Vpe Vxb Vpc Vpe Vxb 表 2 列出了 7 种算法对自然图像 eight 的划 分系数 、划分熵 以及聚类有效性指数 。由 表 2 可以看出，与其余 5 种算法相比，AGWRFCM 算法具有更高的 值和更低的 值，这表明 由 AGWRFCM 算法计算得到的隶属度矩阵划分 性能更好，分割精度更高；同时 AGWRFCM 算法 具有更高的 值，表明类内像素的空间紧致性更 好，不同类像素的空间分离性更强，因此 AGWRFCM 算法具有良好的聚类性能。 Vpc Vpe Vxb 最后，对医学 MR 图像 (MR1，256×256)[29]添加 强度为 0.08 的椒盐噪声，分割结果如图 4 所示， 各算法的划分系数 、划分熵 以及聚类有效性 指数 如表 3 所示。 由图 4 可以看出，在分割含噪医学 MR 图像 时，FCM_S1、FCM_S2、EnFCM、FGFCM、FLICM 以及 WFCM 算法分割图中的背景噪声仍然较多， 抗噪性能较差，且 FLICM 算法分割出的脑组织内 部结构细节信息丢失较严重，因此无法精确地将 MR 脑组织分割出来；而 AGWRFCM 算法既能很 好地抑制噪声，又能保留更多的图像细节，且表 3 中的数据也充分反映出 AGWRFCM 算法具有更 优的聚类分割性能。因此，本文算法对医学图像 也能够取得良好的分割效果。 3.2 改进距离测度的抗噪性测试 由于本文中采用了一种改进的距离测度来计 算相似距离，因此为测试该距离测度对噪声的鲁 表 2 7 种算法对含噪自然图像的聚类性能 Table 2 Clustering performance of seven algorithms on noisy natural image 分割算法 Vpc Vpe Vxb FCM_S1 0.907 4 0.075 5 0.111 0 FCM_S2 0.915 7 0.071 3 0.113 7 EnFCM 0.787 5 0.181 7 0.101 7 FGFCM 0.751 9 0.206 4 0.060 8 FLICM 0.891 2 0.083 2 0.554 8 WFCM 0.829 8 0.142 5 0.255 0 AGWRFCM 0.917 5 0.068 8 0.055 8 表 3 7 种算法对含噪医学图像的聚类性能 Table 3 Clustering performance of seven algorithms on noisy medical image 分割算法 Vpc Vpe Vxb FCM_S1 0.887 3 0.093 4 0.104 1 FCM_S2 0.947 7 0.041 7 0.089 5 EnFCM 0.866 8 0.103 5 0.154 7 FGFCM 0.875 0 0.103 3 0.080 4 FLICM 0.833 3 0.121 2 0.812 9 WFCM 0.901 1 0.095 6 0.183 9 AGWRFCM 0.950 4 0.041 4 0.013 6 (a) 原图像 (b) 噪声图像 (c) FCM_S1 (d) FCM_S2 (e) EnFCM (f) FGFCM (g) FLICM (h) WFCM (i) AGWRFCM 图 4 7 种算法对含噪医学图像的分割结果 Fig. 4 Segmentation results of seven algorithms on noisy medical image 第 4 期 陆海青，等：自适应灰度加权的鲁棒模糊 C 均值图像分割 ·589·
·590· 智能系统学报 第13卷 棒性，需将其与传统的欧氏距离和核诱导距离进 从表4和表5中可以看出，与传统的欧式距 行对比，即在本文算法中分别使用欧氏距离(AG- 离和核诱导距离相比，改进的距离测度在对不同 WFCM、核诱导距离(AGWKFCM)和改进的距离 水平的噪声图像进行分割时均获得了更高的PSNR 测度(AGWRFCM进行实验。 和MSSIM值，且随着噪声水平的提高，PSNR和 分别对人工合成图像(synthetic2,256像素× MSSIM值增加的幅度越来越大。这是由于改进 256像素)、自然图像(cameraman,256像素×256像 的距离测度能够克服噪声或野值数据对聚类性能 素)、医学MR图像(MR,386像素×331像素)2添 加不同强度的混合噪声，其中Mixed、Mixed2、 的影响，从而增强算法的抗噪能力。据上分析可 Mixed:、Mixed.4分别表示均值为0，方差为10、20、 知采用改进距离测度的FCM算法具有更强的抗 30、40的高斯噪声与椒盐噪声的叠加，并采用峰 噪性。 值信噪比(peak signal-to-noise ratio,.PSNR2)和平 图5同时给出了3种算法的PSNR和MSSIM 均结构相似度(mean structural similarity index, 曲线。从中可以更直观地看出，改进的距离测度 MSSIM30叫对上述3种方法的抗噪性能进行客观 相比于其余两种距离测度而言对噪声具有更强的 评价，以进一步验证本文算法对噪声的鲁棒性。 鲁棒性。因此本文算法在抗噪性能上具有一定的 实验结果如表4~5所示。 优越性。 表43种算法的峰值信噪比 32 ·-AGWFCM Table 4 PSNR of three algorithms dB AGWKFCM 30 eAGWRFCM 测试图像噪声类型AGWFCM AGWKFCM AGWRFCM 38 Mixed 30.846831.1736 31.4410 Mixed, 27.3314 27.5658 28.1008 synthetic2 Mixed; 24.3145 24.9942 25.3807 Mixed 21.8262 22.5281 22.6865 Mixed 21.6496 21.6897 21.7010 20 Mixed, 21.3147 21.3823 21.4049 5 1015202530354045 cameraman Mixed; 20.7327 20.7481 20.8024 (a)含噪合成图像的PSNR曲线 Mixed 20.0798 20.1066 20.1748 1.000 -AGWFCM Mixed122.4340 22.6024 22.6438 AGWKFCM 0.995 AGWRFCM Mixed.221.8828 22.0138 22.0376 MR2 Mixed; 20.8546 21.0278 21.0588 Mixeda 19.8566 19.9205 19.9824 表53种算法的平均结构相似度 Table 5 MSSIM of three algorithms 0.980 测试图像噪声类型AGWFCM AGWKFCM AGWRFCM 0.975 Mixed 0.9974 0.9976 0.9977 5 1015202530334045 Mixed2 0.9938 0.9942 0.9949 (b)含噪合成图像的MSSIM曲线 synthetic, Mixed3 0.9872 0.9892 0.9902 22.0 AGWECM Mixed 0.9764 0.9803 0.9811 21.8 AGWKFCM AGWRFCM Mixed 21.6 0.9368 0.9395 0.9405 21.4 Mixed, 0.9304 0.9343 0.9351 cameraman Mixed3 0.9189 0.9229 0.9238 Mixed 0.9039 0.9087 0.9100 20.8 20.6 Mixed 0.9585 0.9611 0.9615 20.4 Mixed2 0.9517 0.9543 0.9547 20.2 MR2 Mixed; 0.9372 0.9417 0.9419 20.0 5 1015202530354045 Mixed 0.9199 0.9239 0.9242 (C)含噪自然图像的PSNR曲线
棒性，需将其与传统的欧氏距离和核诱导距离进 行对比，即在本文算法中分别使用欧氏距离 (AGWFCM)、核诱导距离 (AGWKFCM) 和改进的距离 测度 (AGWRFCM) 进行实验。 分别对人工合成图像 (synthetic2，256 像素× 256 像素)、自然图像 (cameraman，256 像素×256 像 素)、医学 MR 图像 (MR2，386 像素×331 像素) [29]添 加不同强度的混合噪声，其中 Mixed1、Mixed2、 Mixed3、Mixed4 分别表示均值为 0，方差为 10、20、 30、40 的高斯噪声与椒盐噪声的叠加，并采用峰 值信噪比 (peak signal-to-noise ratio，PSNR[22] ) 和平 均结构相似度 (mean structural similarity index， MSSIM[30] ) 对上述 3 种方法的抗噪性能进行客观 评价，以进一步验证本文算法对噪声的鲁棒性。 实验结果如表 4～5 所示。 从表 4 和表 5 中可以看出，与传统的欧式距 离和核诱导距离相比，改进的距离测度在对不同 水平的噪声图像进行分割时均获得了更高的 PSNR 和 MSSIM 值，且随着噪声水平的提高，PSNR 和 MSSIM 值增加的幅度越来越大。这是由于改进 的距离测度能够克服噪声或野值数据对聚类性能 的影响，从而增强算法的抗噪能力。据上分析可 知采用改进距离测度的 FCM 算法具有更强的抗 噪性。 图 5 同时给出了 3 种算法的 PSNR 和 MSSIM 曲线。从中可以更直观地看出，改进的距离测度 相比于其余两种距离测度而言对噪声具有更强的 鲁棒性。因此本文算法在抗噪性能上具有一定的 优越性。 表 4 3 种算法的峰值信噪比 Table 4 PSNR of three algorithms dB 测试图像 噪声类型 AGWFCM AGWKFCM AGWRFCM synthetic2 Mixed1 30.846 8 31.173 6 31.441 0 Mixed2 27.331 4 27.565 8 28.100 8 Mixed3 24.314 5 24.994 2 25.380 7 Mixed4 21.826 2 22.528 1 22.686 5 cameraman Mixed1 21.649 6 21.689 7 21.701 0 Mixed2 21.314 7 21.382 3 21.404 9 Mixed3 20.732 7 20.748 1 20.802 4 Mixed4 20.079 8 20.106 6 20.174 8 MR2 Mixed1 22.434 0 22.602 4 22.643 8 Mixed2 21.882 8 22.013 8 22.037 6 Mixed3 20.854 6 21.027 8 21.058 8 Mixed4 19.856 6 19.920 5 19.982 4 表 5 3 种算法的平均结构相似度 Table 5 MSSIM of three algorithms 测试图像 噪声类型 AGWFCM AGWKFCM AGWRFCM synthetic2 Mixed1 0.997 4 0.997 6 0.997 7 Mixed2 0.993 8 0.994 2 0.994 9 Mixed3 0.987 2 0.989 2 0.990 2 Mixed4 0.976 4 0.980 3 0.981 1 cameraman Mixed1 0.936 8 0.939 5 0.940 5 Mixed2 0.930 4 0.934 3 0.935 1 Mixed3 0.918 9 0.922 9 0.923 8 Mixed4 0.903 9 0.908 7 0.910 0 MR2 Mixed1 0.958 5 0.961 1 0.961 5 Mixed2 0.951 7 0.954 3 0.954 7 Mixed3 0.937 2 0.941 7 0.941 9 Mixed4 0.919 9 0.923 9 0.924 2 AGWFCM AGWKFCM AGWRFCM 5 10 15 20 25 30 35 40 45 20 22 24 26 28 30 32 PSNR (a) 含噪合成图像的PSNR曲线 AGWFCM AGWKFCM AGWRFCM 5 10 15 20 25 30 35 40 45 0.975 0.980 0.985 0.990 0.995 1.000 MSSIM (b) 含噪合成图像的MSSIM曲线 AGWFCM AGWKFCM AGWRFCM 5 10 15 20 25 30 35 40 45 20.0 20.2 20.4 20.6 20.8 21.0 21.2 21.4 21.6 21.8 22.0 PSNR (c) 含噪自然图像的PSNR曲线 σ σ σ ·590· 智 能 系 统 学 报 第 13 卷
第4期 陆海青，等：自适应灰度加权的鲁棒模糊C均值图像分割 ·591· 0.950 AGWECM 3.3邻域窗口大小ω的选取 0.945 AGWKFCM 0.940 AGWRFCM 改进FCM算法涉及邻域窗口大小ω的选取.， 0.935 不同的窗口大小将会对图像分割的效果产生影 三0.930 响。因此有必要对不同噪声条件下ω的取值进行 分析。 0.915 表6给出了在不同窗口大小w下本文算法对 0.910 人工合成图像synthetic,.的分割精度（其中图像的 0.905 噪声方差σ从10~90依次递增)。通过对比可以 0.900 1015202530354045 看出，当σ≤50时，选取3×3窗口进行分割具有更 (d)含噪自然图像的MSSIM曲线 高的分割精度；而当σ>50时，选取5×5或7×7窗 230 ·-AGWFCM 口进行分割可以取得相对较高的分割精度，且随 22.5 -AGWKFCM AGWRFCM 着噪声水平的增加，邻域窗口越大，分割精度越 22.0 高。若此时仍选取3×3窗口，则分割精度将大幅 兰21.5 下降，这是因为图像中包含大量的噪声，3×3窗口 221.0 不足以抑制强度过高的噪声，因此分割性能较 20.5 差。然而若邻域窗口选取过大，算法的分割时间 20.0 将会越长，且容易过多地滤去图像中的纹理细节 953 等特征，从而产生过平滑现象。 101520 2530354045 综上分析可得，对于较低水平的噪声，可选 (e)含噪医学图像的PSNR曲线 取3×3窗口进行分割；对于较高水平的噪声，为更 0.97 AGWFCM 好地维持分割精度与分割时间之间的平衡，选取 AGWKFCM 0.96 AGWRFCM 5×5或7×7窗口进行分割较为合适。 3.4运行时间对比 0.95 为比较不同算法的运行效率，表7中列出了 7种分割算法对不同图像的处理时间，其中运行 0.93 时间取10次运行结果的平均值。 0.92 从表7中可以看出，本文算法的运行时间远 小于FLICM算法，与FCMS,、FCMS2以及WFCM 0.9 51015202530354045 算法相当，但不如EnFCM和FGFCM算法，这是 (①含噪医学图像的MSSIME曲线 因为EnFCM和FGFCM算法是基于灰度级的处 图53种算法的PSNR和MSSIM曲线 理，参与计算的像素数远小于图像的像素总数，因 Fig.5 PSNR and MSSIM curves of three algorithms 此运行效率高于本文算法，但本文算法在视觉效 果和分割质量方面均优于这两种算法。 表6本文算法在不同噪声强度与不同窗口大小下的分割精度对比 Table 6 SA comparison of proposed algorithm under different levels of noise and window sizes 0=10 0=20 0=30 0=40 0=50 0=60 0=70 0=80 0=90 3 0.9998 0.9990 0.9972 0.9934 0.9837 0.7191 0.7050 0.6949 0.6797 0.9988 0.9970 0.9941 0.9924 0.7386 0.7377 0.7349 0.7307 0.7269 7 0.9981 0.9951 0.9911 0.9904 0.7379 0.7375 0.7371 0.7370 0.7358 9 0.9973 0.9928 0.9876 0.9864 0.7372 0.7368 0.7365 0.7363 0.7360 11 0.9960 0.9899 0.9855 0.9845 0.7363 0.7361 0.7357 0.7350 0.7343
3.3 邻域窗口大小ω的选取 ω ω 改进 FCM 算法涉及邻域窗口大小 的选取， 不同的窗口大小将会对图像分割的效果产生影 响。因此有必要对不同噪声条件下 的取值进行 分析。 ω σ σ ⩽ 50 σ > 50 表 6 给出了在不同窗口大小 下本文算法对 人工合成图像 synthetic1 的分割精度 (其中图像的 噪声方差 从 10～90 依次递增)。通过对比可以 看出，当 时，选取 3×3 窗口进行分割具有更 高的分割精度；而当 时，选取 5×5 或 7×7 窗 口进行分割可以取得相对较高的分割精度，且随 着噪声水平的增加，邻域窗口越大，分割精度越 高。若此时仍选取 3×3 窗口，则分割精度将大幅 下降，这是因为图像中包含大量的噪声，3×3 窗口 不足以抑制强度过高的噪声，因此分割性能较 差。然而若邻域窗口选取过大，算法的分割时间 将会越长，且容易过多地滤去图像中的纹理细节 等特征，从而产生过平滑现象。 综上分析可得，对于较低水平的噪声，可选 取 3×3 窗口进行分割；对于较高水平的噪声，为更 好地维持分割精度与分割时间之间的平衡，选取 5×5 或 7×7 窗口进行分割较为合适。 3.4 运行时间对比 为比较不同算法的运行效率，表 7 中列出了 7 种分割算法对不同图像的处理时间，其中运行 时间取 10 次运行结果的平均值。 从表 7 中可以看出，本文算法的运行时间远 小于 FLICM 算法，与 FCM_S1、FCM_S2 以及 WFCM 算法相当，但不如 EnFCM 和 FGFCM 算法，这是 因为 EnFCM 和 FGFCM 算法是基于灰度级的处 理，参与计算的像素数远小于图像的像素总数，因 此运行效率高于本文算法，但本文算法在视觉效 果和分割质量方面均优于这两种算法。 表 6 本文算法在不同噪声强度与不同窗口大小下的分割精度对比 Table 6 SA comparison of proposed algorithm under different levels of noise and window sizes ω σ=10 σ=20 σ=30 σ=40 σ=50 σ=60 σ=70 σ=80 σ=90 3 0.999 8 0.999 0 0.997 2 0.993 4 0.983 7 0.719 1 0.705 0 0.694 9 0.679 7 5 0.998 8 0.997 0 0.994 1 0.992 4 0.738 6 0.737 7 0.734 9 0.730 7 0.726 9 7 0.998 1 0.995 1 0.991 1 0.990 4 0.737 9 0.737 5 0.737 1 0.737 0 0.735 8 9 0.997 3 0.992 8 0.987 6 0.986 4 0.737 2 0.736 8 0.736 5 0.736 3 0.736 0 11 0.996 0 0.989 9 0.985 5 0.984 5 0.736 3 0.736 1 0.735 7 0.735 0 0.734 3 AGWFCM AGWKFCM AGWRFCM 5 10 15 20 25 30 35 40 45 0.900 0.905 0.910 0.915 0.920 0.925 0.930 0.935 0.940 0.945 0.950 MSSIM (d) 含噪自然图像的MSSIM曲线 AGWFCM AGWKFCM AGWRFCM 5 10 15 20 25 30 35 40 45 19.5 20.0 20.5 21.0 21.5 22.0 22.5 23.0 PSNR (e) 含噪医学图像的PSNR曲线 AGWFCM AGWKFCM AGWRFCM 5 10 15 20 25 30 35 40 45 0.91 0.92 0.93 0.94 0.95 0.96 0.97 MSSIM (f) 含噪医学图像的MSSIM曲线 σ σ σ 图 5 3 种算法的 PSNR 和 MSSIM 曲线 Fig. 5 PSNR and MSSIM curves of three algorithms 第 4 期 陆海青，等：自适应灰度加权的鲁棒模糊 C 均值图像分割 ·591·
·592· 智能系统学报 第13卷 表77种算法的平均运行时间对比 算法).重庆邮电大学学报：自然科学版，2017,29(3) Table 7 Average running time comparison of seven al- 377-381 gorithms DING Xiaofeng,HE Kailin.A homogeneous pieces based 分割算法 FCM algorithm for image segmentation[J].Journal of synthetic eight MR Chongqing university of posts and telecommunications: FCM_S1 0.3679 1.5172 0.6477 natural science edition,2017,29(3):377-381. FCM_S2 0.3865 1.5808 0.4999 [5]XIANG Deliang,TANG Tao,HU Canbin,et al.A kernel EnFCM 0.0489 0.1548 0.0626 clustering algorithm with fuzzy factor:application to SAR FGFCM 0.0373 0.1761 0.0411 image segmentation[J].IEEE geoscience and remote sens- ing letters,2014,11(7):1290-1294. FLICM 3.6512 7.3905 4.2930 [6]LIU Guoying,ZHANG Yun,WANG Aimin.Incorporat- WFCM 0.4691 1.9057 1.0940 ing adaptive local information into fuzzy clustering for im- AGWRFCM 0.4953 1.9584 0.8051 age segmentation[J].IEEE transactions on image pro- cessing,2015,2411:3990-4000. 4结束语 [7]SZILAGYI L,BENYO Z,SZILAGYI S M,et al.MR brain image segmentation using an enhanced fuzzy C- 本文提出了一种自适应灰度加权的鲁棒模 means algorithm[Cl//Proceedings of the 25th Annual Inter- 糊C均值图像分割算法，该算法充分利用邻域像 national Conference of the IEEE Engineering in Medicine 素对中心像素影响的差异性，对局部邻域内的各 and Biology Society.Cancun,Mexico,2003:724-726 像素分配合适的权值，以提高像素灰度计算的准 [8]SIKKA K,SINHA N,SINGH P K,et al.A fully auto- 确性；同时采用一种改进的距离测度来计算像素 mated algorithm under modified FCM framework for im- 与聚类中心之间的相似距离，以提高算法对噪声 proved brain MR image segmentation[J].Magnetic reson- 的鲁棒性。大量仿真结果表明，所提算法相比于 ance imaging,2009,27(7):9941004. 经典FCM算法具有更优的分割效果和抗噪性能， [9]BENAICHOUCHE A N.OULHADJ H,SIARRY P.Im 且与其他分割算法相比也具有一定的竞争力。然 proved spatial fuzzy c-means clustering for image seg- 而该算法在进行图像分割时需根据图像噪声的强 mentation using PSO initialization,Mahalanobis distance and post-segmentation correction[J].Digital signal pro- 度选取合适的窗口大小，以取得更好的分割效 cessing,2013,23(23:1390-1400. 果，且算法的运行效率较低。因此下一步将致力 [10]孙立新，张栩之，邓先瑞，等.自适应果蝇算法优化模糊 于研究如何在噪声平滑与边缘保持之间取得更好 均值聚类算法图像分割[J】.控制工程，2016,23(4)： 的平衡，并进一步提升FCM算法的时间性能。 494-499. SUN Lixin,ZHANG Xuzhi,DENG Xianrui,et al.Image 参考文献： segmentation based on self-adaptive fruit fly algorithm [1]BEZDEK J C.Pattern recognition with fuzzy objective optimizing fuzzy means clustering algorithm[J].Control function algorithms[M].New York:Plenum Press,1981. engineering of China,2016,23(4):494-499 [2]唐文静，许兆新，张小峰.峰值检测FCM算法的医学图 [11]ZHANG Daoqiang,CHEN Songcan,PAN Zhisong,et al. 像分割.智能系统学报，2014,9(5)：584-589 Kernel-based fuzzy clustering incorporating spatial con- TANG Wenjing,XU Zhaoxin,ZHANG Xiaofeng.Medic- straints for image segmentation[Cl//Proceedings of the In- ternational Conference on Machine Learning and Cyber- al image segmentation based on FCM with peak detection netics.Xi'an,China,2003:2189-2192. [J].CAAI transactions on intelligent systems,2014,9(5): [12]王建英，孙德山，张永.基于马氏距离的FCM图像分割 584589 算法U.计算机工程与应用，2010,46(1)：147-149. [3]刘慧，张彩明，邓凯，等.改进局部自适应的快速FCM肺 WANG Jianying,SUN Deshan,ZHANG Yong.Mahalan- 结节分割方法).计算机辅助设计与图形学学报，2014， obis distance-based FCM image segmentation algorithm 26(10):1727-1736 [J].Computer engineering and applications,2010,46(1): LIU Hui,ZHANG Caiming,DENG Kai,et al.Research on 147-149. fast FCM pulmonary nodule segmentation algorithm using [13]AHMED M N,YAMANY S M,MOHAMED N,et al.A improved self-adaption[J].Journal of computer-aided modified fuzzy c-means algorithm for bias field estima- design and computer graphics,2014,26(10):1727-1736. tion and segmentation of MRI data[J].IEEE transactions [4]丁晓峰，何凯霖.一种基于一致性分片FCM的图像分割 on medical imaging,2002,21(3):193-199
4 结束语 本文提出了一种自适应灰度加权的鲁棒模 糊 C 均值图像分割算法，该算法充分利用邻域像 素对中心像素影响的差异性，对局部邻域内的各 像素分配合适的权值，以提高像素灰度计算的准 确性；同时采用一种改进的距离测度来计算像素 与聚类中心之间的相似距离，以提高算法对噪声 的鲁棒性。大量仿真结果表明，所提算法相比于 经典 FCM 算法具有更优的分割效果和抗噪性能， 且与其他分割算法相比也具有一定的竞争力。然 而该算法在进行图像分割时需根据图像噪声的强 度选取合适的窗口大小，以取得更好的分割效 果，且算法的运行效率较低。因此下一步将致力 于研究如何在噪声平滑与边缘保持之间取得更好 的平衡，并进一步提升 FCM 算法的时间性能。 参考文献： BEZDEK J C. Pattern recognition with fuzzy objective function algorithms[M]. New York: Plenum Press, 1981. [1] 唐文静, 许兆新, 张小峰. 峰值检测 FCM 算法的医学图 像分割[J]. 智能系统学报, 2014, 9(5): 584–589. TANG Wenjing, XU Zhaoxin, ZHANG Xiaofeng. Medical image segmentation based on FCM with peak detection [J]. CAAI transactions on intelligent systems, 2014, 9(5): 584–589. [2] 刘慧, 张彩明, 邓凯, 等. 改进局部自适应的快速 FCM 肺 结节分割方法[J]. 计算机辅助设计与图形学学报, 2014, 26(10): 1727–1736. LIU Hui, ZHANG Caiming, DENG Kai, et al. Research on fast FCM pulmonary nodule segmentation algorithm using improved self-adaption[J]. Journal of computer-aided design and computer graphics, 2014, 26(10): 1727–1736. [3] [4] 丁晓峰, 何凯霖. 一种基于一致性分片 FCM 的图像分割 算法[J]. 重庆邮电大学学报: 自然科学版, 2017, 29(3): 377–381. DING Xiaofeng, HE Kailin. A homogeneous pieces based FCM algorithm for image segmentation[J]. Journal of Chongqing university of posts and telecommunications: natural science edition, 2017, 29(3): 377–381. XIANG Deliang, TANG Tao, HU Canbin, et al. A kernel clustering algorithm with fuzzy factor: application to SAR image segmentation[J]. IEEE geoscience and remote sensing letters, 2014, 11(7): 1290–1294. [5] LIU Guoying, ZHANG Yun, WANG Aimin. Incorporating adaptive local information into fuzzy clustering for image segmentation[J]. IEEE transactions on image processing, 2015, 24(11): 3990–4000. [6] SZILAGYI L, BENYO Z, SZILAGYI S M, et al. MR brain image segmentation using an enhanced fuzzy Cmeans algorithm[C]//Proceedings of the 25th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. Cancun, Mexico, 2003: 724–726. [7] SIKKA K, SINHA N, SINGH P K, et al. A fully automated algorithm under modified FCM framework for improved brain MR image segmentation[J]. Magnetic resonance imaging, 2009, 27(7): 994–1004. [8] BENAICHOUCHE A N, OULHADJ H, SIARRY P. Improved spatial fuzzy c-means clustering for image segmentation using PSO initialization, Mahalanobis distance and post-segmentation correction[J]. Digital signal processing, 2013, 23(23): 1390–1400. [9] 孙立新, 张栩之, 邓先瑞, 等. 自适应果蝇算法优化模糊 均值聚类算法图像分割[J]. 控制工程, 2016, 23(4): 494–499. SUN Lixin, ZHANG Xuzhi, DENG Xianrui, et al. Image segmentation based on self-adaptive fruit fly algorithm optimizing fuzzy means clustering algorithm[J]. Control engineering of China, 2016, 23(4): 494–499. [10] ZHANG Daoqiang, CHEN Songcan, PAN Zhisong, et al. Kernel-based fuzzy clustering incorporating spatial constraints for image segmentation[C]//Proceedings of the International Conference on Machine Learning and Cybernetics. Xi’an, China, 2003: 2189–2192. [11] 王建英, 孙德山, 张永. 基于马氏距离的 FCM 图像分割 算法[J]. 计算机工程与应用, 2010, 46(1): 147–149. WANG Jianying, SUN Deshan, ZHANG Yong. Mahalanobis distance-based FCM image segmentation algorithm [J]. Computer engineering and applications, 2010, 46(1): 147–149. [12] AHMED M N, YAMANY S M, MOHAMED N, et al. A modified fuzzy c-means algorithm for bias field estimation and segmentation of MRI data[J]. IEEE transactions on medical imaging, 2002, 21(3): 193–199. [13] 表 7 7 种算法的平均运行时间对比 Table 7 Average running time comparison of seven algorithms s 分割算法 synthetic1 eight MR1 FCM_S1 0.367 9 1.517 2 0.647 7 FCM_S2 0.386 5 1.580 8 0.499 9 EnFCM 0.048 9 0.154 8 0.062 6 FGFCM 0.037 3 0.176 1 0.041 1 FLICM 3.651 2 7.390 5 4.293 0 WFCM 0.469 1 1.905 7 1.094 0 AGWRFCM 0.495 3 1.958 4 0.805 1 ·592· 智 能 系 统 学 报 第 13 卷
第4期 陆海青，等：自适应灰度加权的鲁棒模糊C均值图像分割 ·593· [14]CHEN Songcan,ZHANG Daoqiang.Robust image seg- 报，2016,28(4)：615-623 mentation using FCM with spatial constraints based on ZHAO Xuemei,LI Yu,ZHAO Quanhua.A fuzzy cluster- new kernel-induced distance measure[J].IEEE transac- ing image segmentation algorithm with double neighbor- tions on systems,man,and cybernetics.part B:cybernet- hood system combined with Markov Gaussian model[J] ics,2004,344):1907-1916. Journal of computer-aided design and computer graphics, [15]CAI Weiling,CHEN Songcan,ZHANG Daoqiang.Fast 2016,28(4):615-623 and robust fuzzy c-means clustering algorithms incorpor- [24]WU K L,YANG M S.Alternative c-means clustering al- ating local information for image segmentation[J].Pat- gorithms[J].Pattern recognition,2002,35(10):2267- tern recognition,2007,40(3):825-838 2278. [16]LI Ming,LI Yunsong.Fuzzy-c-means clustering based on [25]KRINIDIS S,CHATZIS V.A robust fuzzy local informa- the gray and spatial feature for image segmentation[C]// tion C-means clustering algorithm[J].IEEE transactions Proceedings of the 2006 International Conference on on image processing,2010,19(5):1328-1337 Computational Intelligence and Security.Guangzhou, [26]沙秀艳，何友，王贞俭.邻域灰度差加权的模糊C均值 China,2006:1641-1646. 聚类图像分割算法刀.火力与指挥控制，2008.33(12)： [17刀车娜，车翔玖，高占恒，等.基于局部嫡最小化的核磁共 34-36. 振脑图像二次分割算法).计算机研究与发展，2010， SHA Xiuyan,HE You,WANG Zhenjian.An image seg- 47(7):12941303 mentation algorithm of weighted with neighborhood gray CHE Na,CHE Xiangjiu,GAO Zhanheng,et al.Second- difference fuzzy C-means clustering[J].Fire control and ary segmentation algorithm for magnetic resonance brain command control,2008,33(12):34-36. image based on local entropy minimization[J].Journal of [27]CHUANG K S,TZENG H L,CHEN S,et al.Fuzzy c- computer research and development,2010,47(7):1294- means clustering with spatial information for image seg- 1303. mentation[J].Computerized medical imaging and graph- [18]张建伟，杨红，陈允杰，等.局部嫡驱动下的脑MR图像 ics.2006,30(1):9-15. 分割与偏移场恢复耦合模型[刀.中国图象图形学报 [28]XIE X L.BENI G.A validity measure for fuzzy cluster- 2013,18(8):1011-1018. ing[J].IEEE transactions on pattern analysis and machine ZHANG Jianwei,YANG Hong,CHEN Yunjie,et al. intelligence,1991.13(8):841-847. Brain MR image segmentation and bias correction model [29]KWAN RK S.EVANS A C.PIKE G B.MRI simulation- based on local entropy[J].Journal of image and graphics, based evaluation of image-processing and classification 2013,18(8)1011-1018. methods[J].IEEE transactions on medical imaging,1999, [19]WANG Jianzhong,KONG Jun,LU Yinghua,et al.A 18(11):1085-1097. modified FCM algorithm for MRI brain image segmenta- [30]WANG Zhou,BOVIK A C,SHEIKH H R,et al.Image tion using both local and non-local spatial constraints[J]. quality assessment:from error visibility to structural sim- Computerized medical imaging and graphics,2008,32(8): ilarity[J].IEEE transactions on image processing,2004, 685-698. 13(4)600-612. [20]MA Jingjing,TIAN Dayong,GONG Maoguo,et al. 作者简介： Fuzzy clustering with non-local information for image 陆海青，男，1992年生，硕士研究 segmentation[J].International journal of machine learn- 生，主要研究方向为图像处理和模式 ing and cybernetics,2014,5(6):845-859. 识别。 [21]GUO Yanhui,SENGUR A.NCM:Neutrosophic c-means clustering algorithm[J].Pattern recognition,2015,48(8): 2710-2724. [22]崔西希，吴成茂.核空间中智模糊聚类及图像分割应用 J.中国图象图形学报，2016,21(10)少：1316-1327. 葛洪伟，男，1967年生，教授，博 CUI Xixi,WU Chengmao.Neutrosophic C-means clus- 士生导师，博士，主要研究方向为人工 tering in kernel space and its application in image seg- 智能与模式识别、机器学习、图像处理 mentation[J].Journal of image and graphics,2016 与分析。主持和承担国家自然科学基 21(10:1316-1327. 金等国家级项目和省部级项目近 [23]赵雪梅，李玉，赵泉华.结合马尔可夫高斯模型的双邻 20项，获省部级科技进步奖多项。发 域模糊聚类分割算法J刀.计算机辅助设计与图形学学 表学术论文百余篇
CHEN Songcan, ZHANG Daoqiang. Robust image segmentation using FCM with spatial constraints based on new kernel-induced distance measure[J]. IEEE transactions on systems, man, and cybernetics. part B: cybernetics, 2004, 34(4): 1907–1916. [14] CAI Weiling, CHEN Songcan, ZHANG Daoqiang. Fast and robust fuzzy c-means clustering algorithms incorporating local information for image segmentation[J]. Pattern recognition, 2007, 40(3): 825–838. [15] LI Ming, LI Yunsong. Fuzzy-c-means clustering based on the gray and spatial feature for image segmentation[C]// Proceedings of the 2006 International Conference on Computational Intelligence and Security. Guangzhou, China, 2006: 1641–1646. [16] 车娜, 车翔玖, 高占恒, 等. 基于局部熵最小化的核磁共 振脑图像二次分割算法[J]. 计算机研究与发展, 2010, 47(7): 1294–1303. CHE Na, CHE Xiangjiu, GAO Zhanheng, et al. Secondary segmentation algorithm for magnetic resonance brain image based on local entropy minimization[J]. Journal of computer research and development, 2010, 47(7): 1294– 1303. [17] 张建伟, 杨红, 陈允杰, 等. 局部熵驱动下的脑 MR 图像 分割与偏移场恢复耦合模型[J]. 中国图象图形学报, 2013, 18(8): 1011–1018. ZHANG Jianwei, YANG Hong, CHEN Yunjie, et al. Brain MR image segmentation and bias correction model based on local entropy[J]. Journal of image and graphics, 2013, 18(8): 1011–1018. [18] WANG Jianzhong, KONG Jun, LU Yinghua, et al. A modified FCM algorithm for MRI brain image segmentation using both local and non-local spatial constraints[J]. Computerized medical imaging and graphics, 2008, 32(8): 685–698. [19] MA Jingjing, TIAN Dayong, GONG Maoguo, et al. Fuzzy clustering with non-local information for image segmentation[J]. International journal of machine learning and cybernetics, 2014, 5(6): 845–859. [20] GUO Yanhui, SENGUR A. NCM: Neutrosophic c-means clustering algorithm[J]. Pattern recognition, 2015, 48(8): 2710–2724. [21] 崔西希, 吴成茂. 核空间中智模糊聚类及图像分割应用 [J]. 中国图象图形学报, 2016, 21(10): 1316–1327. CUI Xixi, WU Chengmao. Neutrosophic C-means clustering in kernel space and its application in image segmentation[J]. Journal of image and graphics, 2016, 21(10): 1316–1327. [22] 赵雪梅, 李玉, 赵泉华. 结合马尔可夫高斯模型的双邻 域模糊聚类分割算法[J]. 计算机辅助设计与图形学学 [23] 报, 2016, 28(4): 615–623. ZHAO Xuemei, LI Yu, ZHAO Quanhua. A fuzzy clustering image segmentation algorithm with double neighborhood system combined with Markov Gaussian model[J]. Journal of computer-aided design and computer graphics, 2016, 28(4): 615–623. WU K L, YANG M S. Alternative c-means clustering algorithms[J]. Pattern recognition, 2002, 35(10): 2267– 2278. [24] KRINIDIS S, CHATZIS V. A robust fuzzy local information C-means clustering algorithm[J]. IEEE transactions on image processing, 2010, 19(5): 1328–1337. [25] 沙秀艳, 何友, 王贞俭. 邻域灰度差加权的模糊 C 均值 聚类图像分割算法[J]. 火力与指挥控制, 2008, 33(12): 34–36. SHA Xiuyan, HE You, WANG Zhenjian. An image segmentation algorithm of weighted with neighborhood gray difference fuzzy C-means clustering[J]. Fire control and command control, 2008, 33(12): 34–36. [26] CHUANG K S, TZENG H L, CHEN S, et al. Fuzzy cmeans clustering with spatial information for image segmentation[J]. Computerized medical imaging and graphics, 2006, 30(1): 9–15. [27] XIE X L, BENI G. A validity measure for fuzzy clustering[J]. IEEE transactions on pattern analysis and machine intelligence, 1991, 13(8): 841–847. [28] KWAN R K S, EVANS A C, PIKE G B. MRI simulationbased evaluation of image-processing and classification methods[J]. IEEE transactions on medical imaging, 1999, 18(11): 1085–1097. [29] WANG Zhou, BOVIK A C, SHEIKH H R, et al. Image quality assessment: from error visibility to structural similarity[J]. IEEE transactions on image processing, 2004, 13(4): 600–612. [30] 作者简介： 陆海青，男，1992 年生，硕士研究 生，主要研究方向为图像处理和模式 识别。 葛洪伟，男，1967 年生，教授，博 士生导师，博士，主要研究方向为人工 智能与模式识别、机器学习、图像处理 与分析。主持和承担国家自然科学基 金等国家级项目和省部级项目近 20 项，获省部级科技进步奖多项。发 表学术论文百余篇。 第 4 期 陆海青，等：自适应灰度加权的鲁棒模糊 C 均值图像分割 ·593·