基于空间信息的模糊C-均值噪声图像分割算法.pdf

年 无线电工程 第 卷 第 期：引用格式：李力，陈息坤基于空间信息的模糊均值噪声图像分割算法无线电工程，（）：，（）：基于空间信息的模糊均值噪声图像分割算法李力，陈息坤（广东科学技术职业学院 物联网工程学院，广东 广州；上海大学 机电工程与自动化学院，上海）摘要：针对传统模糊均值（，）聚类算法对噪声鲁棒性差的问题，提出一种基于空间信息的模糊均值噪声图像分割算法。将区域级信息加入目标函数中，并用核度量方法代替传统欧氏距离，计算区域级空间信息与聚类中心的距离，提高算法对噪声的鲁棒性；用原始图像与区域级空间信息的绝对差的倒数和其本身约束原始图像和区域信息项，实现约束项参数的自适应选择；利用连通分量滤波，消除聚类结果中出现的过分割现象，提高分割精度。含噪合成图像和彩色图像实验表明，所提算法在模糊分割系数、模糊分割熵、分割精确度、平均交互比和归一化互信息等方面均优于其他几种聚类算法。关键词：噪声图像分割；模糊均值聚类；区域级信息约束；核度量方法；连通分量滤波中图分类号：；文献标志码：开放科学（资源服务）标识码（）：文 章 编 号：（），（，；，）：（），：；收稿日期：基金项目：广东省普通高校工程技术中心“智能装备制造工程技术研究中心”（）：（）引言图像分割是根据像素之间的相似度（包括亮度、颜色和纹理等）将数字图像分为若干个像素组，从而简化后续分析的过程。近年来，针对不同的应用场合，出现了多种图像分割算法，如阈值法、区域划分法、边缘检测法和深度学习方法等。在真实环境中，由于采集和转换过程中的噪声干扰、图像的质量下降，用上述方法难以产生有效的分割信号与信息处理 结果。比如面对噪声时，阈值法出现目标和背景混叠的现象，区域划分法产生过分割和欠分割，边缘检测法分割的结果边缘贴合度差，难以满足实际需求，基于深度学习的方法需要大量的噪声图像的学习样本训练模型，复杂度和成本较高。模糊均值（，）聚类算法是一种无监督训练算法，因其不需要提前训练样本、算法结构简单和聚类效果好等优点被广泛应用于图像分割领域。由于传统聚类算法未考虑图像的空间信息，因此对噪声较为敏感，为提高传统算法的鲁棒性和有效性，已经提出许多变体算法以实现更好的性能。目前对于算法的改进大致可分为以下种：基于局部空间信息的算法，该方法考虑了中心像素以外的其他邻域像素，提高了算法的鲁棒性，如模糊局部信息均值算法（）、基于空间约束和偏差校正的模糊均值算法（）和基于空间约束和隶属度链接的模糊均值算法（）；基于直方图的方法，该方法将聚类过程由原始图像转换为在直方图上进行聚类，提高了算法的时间效率，如增强型模糊均值算法（）、快速广义模糊均值算法（）和快速超像素模糊均值算法（）；基于区域级信息的算法，该方法通过考虑更多的图像冗余信息以提高算法的鲁棒性，如基于区域信息和均值隶属度链接的模糊子空间聚类算法（）、区域级信息模糊均值算法（）和自适应区域级模糊均值算法（）；基于特征加权的算法，该方法主要通过考虑图像的多个特征以加强对彩色图像的分割性能，如特征和聚类加权模糊均值算法（）、聚类加权和组特征加权的模糊均值算法（）。虽然上述算法提高了对噪声的鲁棒性，但其局部空间信息或区域级信息与聚类中心的距离采用欧式距离进行计算。尽管此度量方法在计算上很简单，但是，在分割已被噪声、离群值或其他图像伪影损坏的图像时，使用欧氏距离可能会产生不稳健的结果，因此需要有一种新的度量方法代替欧氏距离，提高算法的鲁棒性。核度量方法由此引入，基于核函数的模糊均值算法（）在的基础上用核度量方法代替欧氏距离，分割性能显著提高。该算法虽然能对噪声图像进行有效分割，但其对局部空间的约束参数无法进行自适应选择，且聚类结果出现过分割现象，影响分割精度。本文提出一种基于空间信息的模糊均值噪声图像分割算法。首先，将区域级信息加入目标函数中，并用核度量方法代替传统欧氏距离，计算区域级空间信息与聚类中心的距离，提高算法对噪声的鲁棒性；其次，用原始图像与区域级空间信息绝对差的倒数和它本身来约束原始图像和区域信息项，实现约束项参数自适应选择；最后，利用连通分量滤波，消除聚类结果中出现的过分割现象，提高分割精度。相关工作 聚类算法传统聚类算法通过反复迭代，使目标函数达到最小，将图像中的个像素指定为类：，（）式中：为第个聚类簇的中心，为第个像素属于第个聚类簇的隶属度，为隶属度权重指数，通常设置为。通过拉格朗日数乘法，式（）通过式（）和式（）最小化：()，（）。（）区域级信息传统没有考虑局部空间信息，对噪声比较敏感，为了增强算法对噪声的鲁棒性，将局部空间信息引入目标函数中。但在高噪声图像或混合噪声图像中，中心像素的相邻像素可能大都是噪声像素，因此，考虑其他局部像素块比以中心像素为中心的像素块更合理，区域级信息由此引入目标函数中，区域级信息表示为：（）（），()，（）（）（），()，（）式中：为原始图像的第个像素，为落在区域（）的第个像素，（）和（）分别为以像素和像素为中心的局部像素块，大小为（），为高斯标准差，一般设置为，为平滑度参数，一般设置为。局部像素块与区域级信息的关系如图所示。信号与信息处理 年 无线电工程 第 卷 第 期 图图像中的局部像素块与区域级信息的关系 本文算法传统算法对噪声图像较为敏感，难以得到理想的分割结果。本文考虑了图像的区域级信息，并用核度量方法代替传统欧氏距离，计算区域级空间信息与聚类中心的距离，然后用原始图像与区域级空间信息绝对差的倒数和它本身来约束原始图像和区域信息项，最后用连通分量滤波，消除聚类过程中出现的过分割现象。基于核度量的区域级信息提取为了提高的抗噪声性能，许多变体算法被提出，其目标函数可以概括为：，（）式中：为模糊度因子，为约束项，不同的算法具有不同的模糊度因子。在基于区域级信息的模糊均值算法（）中，模糊度因子为：（）（），()。（）观察区域级信息约束项，所在中前面部分为原始图像的区域级信息，则可表示为，其中为图像的区域级信息，因此可以看作区域级信息与聚类中心的欧氏距离计算的结果。尽管此度量方法在计算上很简单，但是使用欧氏距离可能会导致对因噪声、离群值和其他成像伪影而损坏的图像进行分割时出现非稳健的结果。为提高鲁棒性，可将核度量方法引入目标函数代替欧氏距离。特征空间中的核可以表示为以下函数：（），（），（）式中：（）表示隐式非线性映射，（），（）为内积运算。高斯径向基函数核（）是一种常用的核算法，其表达式为：（，）（），（）式中：为向量的维数，为内核带宽，一般设置为，显然（，）。（）和（）在特征空间的内积为（）（）（，），通过核替换有：（）（）（，）（，）（，）。（）通过这种方法，得到了原始数据空间中一类新的非欧几里德距离测度，由于（，），因此模糊因子可以用下式表示：()()。（）综上所述，基于核度量的区域级信息提取步骤如下：首先，选取图像的一个像素为中心像素，记为，其周围 的邻域范围为局部区域，记为（）；其次，在以为中心的区域内搜寻一个与灰度值相同（或相差不大）的像素，记为，其局部区域为（），则由上述个局部区域（）和（）构成的区域即为以像素为中心的区域级信息；最后，利用式（）和式（）计算区域级信息，并将其引入目标函数中，遍历整幅图像，直至所有像素完成聚类。连通分量滤波消除过分割现象聚类算法将每个像素视为独立的样本，因此在聚类过程中，容易产生过分割现象，即分割结果中包含大量孤立的小区域，如图所示。图使用连通分量滤波前后的对比 信号与信息处理 由图（）可以看出，聚类结果包含很多无用的小区域，从而影响分割精度，利用连通分量滤波可以有效地去除这些小区域。连通分量滤波技术，首先计算所有连通区域的面积，然后由大到小对这些连通分量进行排序，生成面积决策图，接着求出各个连通区域的区间长度，最大区间对应一个区域，该区域的面积被认为是阈值，最后用得到的消除孤立区域。文献提出可利用面积密度平衡算法与连通分量滤波相结合来快速寻找阈值。用表示第个数据点的归一化面积，其中，表示第个区域的面积，表示在半径下附近的的个数，可用公式表示为：，（），其他，（）式中：为由聚类算法得到的连通区域的个数。为的映射：。（）将归一化的值作为阈值进行连通分量滤波消除聚类结果中出现的过分割现象。如图（）所示，经过连通分量滤波后，图（）中孤立的小区域被滤除，进而提高了分割精度。本文算法目标函数为了提高传统算法对噪声的鲁棒性，提出一种自适应区域级信息约束方法，目标函数为：（）()()，（）式中：为迭代次数，和为和区域级信息的约束，当图像受噪声影响较大时，需要选择较小的和较大的，因此可以在一定程度上保留了原始图像的同时，还能有效地约束区域级信息。对于噪声图像，像素受噪声干扰越大，原始像素与区域级像素的强度差越大，反之越小，因此用原始像素与区域级像素的强度差的倒数和它本身来自适应约束和区域级信息。因此令，同时为了防止和趋于无穷大和，在等式中加入一个常数和，为一个极小的数，一般设置为，经过多次实验，得到较好的分割结果。由此，最终，。使用拉格朗日数乘法，得到和：（）（），（）（）（）（）（）。（）通过式（）和式（）更新和，使目标函数最小化，在输出的隶属度中，将每个像素指定为隶属度最大的聚类簇中，完成图像分割。本文算法步骤如算法所示。算法：基于区域信息和连通分量滤波的算法输入：图像，聚类数目，误差参数，最大迭代次数，局部像素块尺寸，搜索半径输出：隶属度矩阵，聚类中心计算原始图像的区域级信息初始化隶属度矩阵、聚类中心利用式（）计算区域级信息与聚类中心的距离 重复执行用式（）更新聚类中心用式（）更新隶属度矩阵用式（）更新目标函数（）直到（）（）或返回隶属度矩阵，聚类中心实验结果与分析为了验证本文算法的有效性，采用合成图像和彩色图像进行测试，并将实验结果与、和算法进行对比。为了说明算法对不同类型的噪声都具有较强的鲁棒性，实验选用等量的高斯白噪声、椒盐噪声和均匀分布乘性噪声的混合噪声污染图像。指标采用模糊分割系数、模糊分割熵、分割精准度、平均交并比和归一化互信息：，（）（），（），（），（）信号与信息处理 年 无线电工程 第 卷 第 期 式中：为分割结果中第个聚类簇的像素集合，为参考图像中第个聚类簇的像素集合。、和越大，分割结果越好，越小分割性能越好。对于图像和：（，）（）（），（）式中：（，）表示和的互信息，（）和（）分别表示和的熵，越大，分割结果越好。对几种算法的参数设置如下：对于所有算法，设置隶属度指数 ，最小误差，最大迭代次数 ；对于算法，、和分别设置为、和；对于和算法，设置局部窗口，搜索窗口，平滑度参数；对于本文算法的参数敏感性将在 节讨论。合成图像分割实验原始图像大小为 ，在实验过程中分别加入、和的混合噪声，聚类数目 。图为种算法对含混合噪声合成图像的分割结果，种算法对含混合噪声合成图像的定量指标结果如表所示。首先，讨论分割的视觉效果。视觉效果最差的为和算法，因为算法没考虑图像的任何空间信息，对噪声较为敏感，而算法是基于特征的聚类算法，即在算法的输入端输入图像的特征而非图像本身，而在面对噪声时，原始图像的特征被噪声破坏，难以表征图像信息，因此的结果视觉效果较差。和算法虽然考虑了图像的空间信息，但度量像素与聚类中心的距离仍然采用原始的欧氏距离，因此结果中含有一部分噪声斑点，且存在误分割现象，视觉效果较差。本文算法的视觉效果最好，这主要归因于区域级信息和核度量方法。其次，分析各个算法的性能指标。由表可以看出，总体上，随着噪声密度的增加，所有算法的性能均有所下降。其中，和算法的性能指标最低，如在噪声密度为时，和算法的指标分别为 和，表明在面对高密度噪声时，这种算法无法正确执行。算法的性能指标较好，接近本文算法，个别指标甚至超过本文算法，如在噪声密度为时，算法的为，超过本文算法。本文算法在大多数指标上性能均优于其他几种比较算法，稳定性较好。图 种算法对含混合噪声合成图像的分割结果 表 种算法对含不同混合噪声的合成图像的性能指标 测试算法分割指标及其混合噪声密度 本文算法 彩色图像分割实验 节的合成图像实验说明本文算法是一种鲁棒性较强的图像分割算法，本小节采用真实的彩色图像验证算法的有效性。为了实验的严谨性，本文测试了来自伯努利数据库（）的所有图像，并给出了所有算法的平均性能指标，如表所信号与信息处理 示。同时为了展示视觉效果，本文选取几幅具有代表性的图像，将所有算法的分割结果展示在图中。由图可以看出，算法的分割结果中充满噪声，视觉效果较差。算法存在严重的误分割现象，几乎无法辨别图像中的目标和背景，这表明，算法虽然考虑了图像的多个特征，但面对噪声时，其稳健性逐渐降低。算法对、的分割结果较好，但在分割时存在误分割，将飞机和天空误分为一类，和的结果中存在少许噪声斑点。算法是最近提出的一种基于非局部空间信息的变体算法，其性能较好，因此可以分割大多的彩色图像，但也有例外，如的分割结果中目标比较模糊，的分割结果含有少量噪声。相比前几种算法，本文算法不仅能有效地分割目标，同时可以较好地抑制噪声，视觉效果最好。在性能指标方面，本文给出了所有算法在数据库中的平均性能指标。与合成图像类似，随着噪声的增加，各个指标逐渐降低，但相较于其他算法，本文算法下降得较为缓慢，即本文算法具有更好的稳定性。这主要归因于本文算法考虑了图像的区域信息，同时又以核函数为距离度量，在面对噪声时稳定性更好。表 种算法对含不同混合噪声彩色图像的性能指标 测试算法分割指标及其混合噪声密度 本文算法 图 种算法对含混合噪声彩色图像的分割结果 信号与信息处理 年 无线电工程 第 卷 第 期 参数敏感性分析讨论参数对算法性能的影响，即局部窗口尺寸和搜索窗口尺寸的选择。首先，确定的值，为局部窗口尺寸，决定邻域范围的大小。如果太小，则无法有效滤除噪声，而过大的导致空间运算（即区域信息的提取过程）过分平滑图像，使得图像的细节和边缘信息丢失，进而影响分割精度，因此需要选择合适的值。如图所示，随着增大，算法的性能也相应提升，当增加到时，其性能逐渐趋于稳定，因此本文设置 。其次，确定的值，为搜索窗口的尺寸，当较小时，算法执行速度较快，但鲁棒性较差，当较大时，可以有效滤除噪声，但算法的时间效率也相对降低，因此选择合适的值是至关重要的。如图所示，随着的增加，算法的性能也在逐渐提升，但当增加到之后，其各个性能增加缓慢，综合时间效率和分割性能考虑，本文设置。图不同对算法性能的影响 图不同对算法性能的影响 时间复杂度比较时间复杂度是评价算法性能的重要指标，表列出了种算法的时间复杂度，其中为像素总个数，为特征数，为迭代步长。由表可知，算法的复杂度最低，因为原始算法并未考虑任何空间信息，因此最为简单。复杂度较低的为算法，该算法只考虑图像的特征而没有进行空间运算，因此复杂度较低，本文算法和其他几种算法都是基于区域级信息的算法，因此复杂度较高，为（），但同时鲁棒性也较好。表 种算法的复杂度比较 算法计算步骤复杂度（）（）（）（）本文算法（）结束语本文提出了一种基于空间信息的模糊均值噪声图像分割算法。针对传统聚类算法对噪声缺乏鲁棒性的问题，将原始图像的区域级信息加入的目标函数之中，并用核度量方法代替欧氏距离，计算区域级信息与聚类中心之间的距离，提高了算法对噪声的鲁棒性，而且将原始图像与区域级空间信息绝对差的倒数和它本身作为约束项，来自适应约束原始图像和与区域级信息，最后再利用连通分量滤波对聚类结果中出现的孤立小区域进行去除，提高分割精度。通过合成图像和彩色图像分割实验对比，表明本文算法在模糊分割系数、模糊分割熵、分割精确度、平均交互比和归一化互信息等方面均优于、和算法。但本文算法输入参数较多，如何减少输入参数将是下一步的研究目标。?参考文献，：，：王小鹏，王庆圣，焦建军，等快速自适应非局部空间加权与隶属度连接的模糊均值噪声图像分割算法电子与信息学报，（）：信号与信息处理 杜吉生，张华忠基于高分影像的冰灾遥感监测方法研究无线电工程，（）：王莹，朱家明，徐婷宜，等融合深度残差结构的脑肿瘤分割无线电工程，（）：张淑军，彭中，李辉：基于和自注意力机制的医学图像分割方法电子学报，（）：，（）：，（）：，（）：，：，：，：，：，：，（）：，（）：，：，（）：，：，：，：，（）：，（）：，（）：，：，（）：作者简介李力女，（），副教授。主要研究方向：控制工程技术、计算机应用技术等。陈息坤男，（），博士，教授，博士生导师。主要研究方向：新型储能变换与控制技术、脉冲功率技术等。信号与信息处理