一种基于区域生长和蚁群聚类的图像分割方法技术

本发明专利技术涉及一种基于区域生长和蚁群聚类的图像分割方法，技术特征在于：首先，由于考虑的前提是从灰度值最高的像素点开始生长，因此去除图像中的噪声很重要。因此采用四邻域的中值滤波去滤除噪声；接着，选择具有灰度值最大的像素点作为种子点进行区域生长。然后利用区域生长后提到的空间信息和灰度信息定义一种新的引导函数，用到蚁群算法中在区域之间进行聚类合并，获得最终的分割结果。本发明专利技术提出的一种基于区域生长和蚁群聚类的图像分割方法，具有三个明显的优点：第一，克服了区域生长得不到有意义区域的不足；第二，大大提高了蚁群聚类算法的搜索时间。第三，新的引导函数的定义可准确有效引导蚁群聚类，提高图像分割的准确率。

【技术实现步骤摘要】

本专利技术涉及，属于计算机视觉、图像理解以及模式识别和人工智能等领域。
技术介绍
图像分割就是指把图像分成各具特性的区域并提取出感兴趣目标的技术和过程。图像在分割后的处理，如特征提取、目标识别等都依赖图像分割的质量，所以图像分割一直是计算机视觉和模式识别领域的研究热点。 目前已经有多种方法被人们提出应用于图像分割领域，传如阈值法、边缘检测法、数学形态学法、基于区域处理方法等，这些方法针对不同图像都取得了很好的效果。但是对于不同应用目的的以及不同图像特性，上述方法又表现出很大局限性。例如阈值法，这种方法较高的计算效率，但是对噪声敏感，会误将噪声作为目标来处理；边缘检测算子存在边界不连续或边界不准确的问题；数学形态学方法在一定程度上降低了噪声对图像的影响，但是开、闭、腐蚀、等运算会导致图像的过度平滑，从而导致图像变形及细节丢失。 越来越多的学者开始将模糊理论、马尔科夫模型、遗传算法理论、分形理论和小波理论等研究成果运用于图像分割的研究，取得了很大进展。尽管人们在图像分割方面做了许多研究工作，由于尚无通用的分割理论，现已提出了大量算法都是针对具体问题的，并没有一种适合于所有图像的通用的分割方法。蚁群算法是一种具有离散性、并行性、鲁棒性和模糊聚类能力的进化方法。它的这些特点对于离散的数字图像非常适用。而现有的蚁群算法在图像分割的应用中其收敛时间过长，易陷入局部最优。
技术实现思路
要解决的技术问题 为了避免现有技术的不足之处，本专利技术提出，该方法首先用区域生长法对图像做初始分割，然后利用初始分割后的空间信息和灰度信息定义了一种新的引导函数，再利用蚁群算法搜索最优解的能力，在区域之间进行聚类合并，获得最终的分割结果。 技术方案 本专利技术的基本思想是利用具有离散性、并行性、正反馈性、鲁棒性和模糊聚类能力的蚁群算法做后继的分割。 本专利技术的技术特征在于步骤如下 1)对于M1×N的图像采用四邻域的中值滤波去滤除噪声； 2)采用区域生长的方法对图像进行预处理 2.1、在没有被标记的像素中选择灰度最大点为种子点； 2.2、根据式|Iseed-I|＜λ|Imax-Imin|，取λ＝0.3，进行区域生长，将满足此公式的邻接像素加入到种子区域，并对它们进行标记；其中I表示像素的灰度值，Iseed表示种子点的灰度值，Imax与Imin分别表示图像中的最大灰度值与最小灰度值，λ是可调节的参数； 2.3、搜索未被标记的像素，进行上述的两步，当不再有未被标记的像素时，输出被划分的区域； 3)提取出灰度信息和空间信息特征 所述的提取灰度信息特征区域生长后计算区域Ri和Rj之间的平均灰度差gij＝|Ii-Ij|，其中其中I是区域生长后得到小区域的平均灰度值，M是区域生长后得到小区域数(即R1，R2，…，RM)，h代表每个小区域中的灰度值，有k个灰度级，即h1，h2，…，hk； 所述的提取空间信息特征根据得到空间特征；其中dij表示Xi到Xj之间的加权欧氏距离，即某区域Xi到不同食物源的距离，m是蚂蚁的维数，Pk为加权因子，加权因子可以根据数据各分量对聚类的影响程度来设定； 4)利用蚁群算法进行聚类合并和分割 4.1将区域生长后的图像区域数M看作每个蚂蚁，初始化N，r，ε0，α，β，τs(0)＝0，P0； 4.2根据式计算各路径上的信息量τij，其中r表示聚类半径；计算引导函数式ηij′，根据式计算Xi归并到Xj的概率pij(t)，其中，S是蚂蚁Xs下一步可以选择的路径集合即S＝{Xs|dsj≤r，s＝1，2，…，N}；α，β为调节因子，起到防止所有蚂蚁沿相同路径得到相同结果所产生的停滞搜索的作用； 4.3根据判断原则1Pij(t)≥P0成立，则Xi归并到Xj邻域内，继续执行下一步；根据判断原则2Pij(t)＜P0，则转到步骤3中计算i+1点到j点的灰度信息特征gi+1，j和di+1，j，根据步骤4.2计算概率Pi＝1，j(t)，按照判断原则1或原则2决定执行步骤； 4.4根据Xk∈Cj计算聚类中心Cj，其中Cj表示所有归并到Xj邻域的数据集合，J为Cj类中元素的个数，即Cj＝{Xk|dkj≤r，k＝1，2，…，J}； 4.5计算各类的类间距离，当类间距小于阈值ε时，将两类合并为一类，更新聚类中心； 4.6如果还有待分类区域，则返回步骤3，否则输出聚类结果，得到分割的图像。 有益效果 本专利技术提出的，具有三个明显的优点第一，克服了区域生长得不到有意义区域的不足；第二，大大提高了蚁群聚类算法的搜索时间。第三，新的引导函数的定义可准确有效引导蚁群聚类，提高图像分割的准确率。 附图说明 图1基于区域生长和蚁群聚类的图像分割流程图 图2基于区域生长和蚁群聚类的图像分割结果1 a原图像 b区域生长后的图像 c基本蚁群算法得到的图像分割结果 d本专利技术算法得到的图像分割结果 图3基于区域生长和蚁群聚类的图像分割结果2 a原图像 b区域生长后的图像 c基本蚁群算法得到的图像分割结果 d本专利技术算法得到的图像分割结果 具体实施例方式 现结合附图对本专利技术作进一步描述 首先，由于考虑的前提是从灰度值最高的像素点开始生长，因此去除图像中的噪声很重要。因此采用四邻域的中值滤波去滤除噪声；接着，选择具有灰度值最大的像素点作为种子点进行区域生长。然后利用区域生长后提到的空间信息和灰度信息定义一种新的引导函数，用到蚁群算法中在区域之间进行聚类合并，获得最终的分割结果。 由于考虑的前提是从灰度值最高的像素点开始生长，因此去除图像中的噪声很重要。考虑到需要去除的只是一些灰度值不正常(偏高或者偏低)的噪声点，采用四邻域的中值滤波，结果显示对不正常的像素点有很好的效果。 用区域生长方法对图像进行初始划分。该方法利用图像区域或像素之间的连续性与邻接性来进行处理。根据事前定义的规则将像素或子区域聚合成更大的区域。它的基本思路是从一个或多个种子点出发，不断地加入满足相似性规则的邻居点来生长出图像区域。本专利技术的区域生长过程如下 1)种子点的选择上，每次选择具有灰度值最大的像素点作为种子点进行区域生长。 2)在空间上采用八邻域连通方案对邻接的相似像素进行搜索。 3)在相似性准则的选取上，定义如下的公式用于选择邻近的像素 |Iseed-I|＜λ|Imax-Imin| 其中I表示像素的灰度值，Iseed表示种子点的灰度值，Imax与Imin分别表示图像中的最大灰度值与最小灰度值，λ是可调节的参数，用来控制像素之间的相似度门限，将满足此公式的邻接像素加入到种子区域。 4)生长的过程中当没有像素满足加入某个种子区域的条件时，区域生长中止。 在实现上，程序递归地调用该算法直到所有的像素都被划分区域。当区域生长完成时，输出的是一系列空间上连续的种子区域。本专利技术在区域生长完后将一些过于琐碎的区域(总像素少于10)简单地并入到与其邻接的相似度最接近的一个区域中，因为这些区域的数量比较多，这样做可以避免复杂的计算量，同时也影响不到图像中的主要信息，而且可以使蚁群聚类的问题求解规模减小。 本专利技术用区域生长后的灰度信息和空间信息来改变基本蚁群聚类算法中的引导函数本文档来自技高网...

【技术保护点】
一种基于区域生长和蚁群聚类的图像分割方法，其特征在于步骤如下：　　　　１）对于Ｍ１×Ｎ的图像采用四邻域的中值滤波去滤除噪声；　　　　２）采用区域生长的方法对图像进行预处理：　　　　２．１、在没有被标记的像素中选择灰度最大点为种子点；　　　　２．２、根据式｜Ｉ↓［ｓｅｅｄ］－Ｉ｜＜λ｜Ｉ↓［ｍａｘ］－Ｉ↓［ｍｉｎ］｜，取λ＝０．３，进行区域生长，将满足此公式的邻接像素加入到种子区域，并对它们进行标记；其中Ｉ表示像素的灰度值，Ｉ↓［ｓｅｅｄ］表示种子点的灰度值，Ｉ↓［ｍａｘ］与Ｉ↓［ｍｉｎ］分别表示图像中的最大灰度值与最小灰度值，λ是可调节的参数；　　　　２．３、搜索未被标记的像素，进行上述的两步，当不再有未被标记的像素时，输出被划分的区域；　　　　３）提取出灰度信息和空间信息特征：　　　　所述的提取灰度信息特征：区域生长后计算区域Ｒ↓［ｉ］和Ｒ↓［ｊ］之间的平均灰度差ｇ↓［ｉｊ］＝｜Ｉ↓［ｉ］－Ｉ↓［ｊ］｜，其中Ｉ＝１／ｋ＊ｈ↓［ｉ］；其中：Ｉ是区域生长后得到小区域的平均灰度值，Ｍ是区域生长后得到小区域数（即Ｒ↓［１］，Ｒ↓［２］，…，Ｒ↓［Ｍ］），ｈ代表每个小区域中的灰度值，有ｋ个灰度级，即ｈ↓［１］，ｈ↓［２］，…，ｈ↓［ｋ］；　　　　所述的提取空间信息特征：根据＊＊＊得到空间特征；其中ｄ↓［ｉｊ］表示Ｘ↓［ｉ］到Ｘ↓［ｊ］之间的加权欧氏距离，即某区域Ｘ↓［ｉ］到不同食物源的距离，ｍ是蚂蚁的维数，Ｐ↓［ｋ］为加权因子，加权因子可以根据数据各分量对聚类的影响程度来设定；　　　　４）利用蚁群算法进行聚类合并和分割：　　　　４．１将区域生长后的图像区域数Ｍ看作每个蚂蚁，初始化Ｎ，ｒ，ε↓［０］，α，β，τ↓［ｓ］（０）＝０，Ｐ↓［０］；　　　　４．２根据式＊＊＊计算各路径上的信息量τ↓［ｉｊ］，其中ｒ表示聚类半径；计算引导函数式η′↓［ｉｊ］，η′↓［ｉｊ］＝ｒ／ｄ↓［ｉｊ］＋ｇ↓［ｉｊ］＝ｒ／＊＊＊＋｜Ｉ↓［ｉ］－Ｉ↓［ｊ］｜；根据式ｐ↓［ｉｊ］（ｔ）＝τ↓［ｉｊ］↑［α］（ｔ）η′↓［ｉｊ］↑［β］（ｔ）／＊τ↓［ｓｊ］↑［α］（ｔ）η′↓［ｓｊ］↑［β］（ｔ）计算Ｘ↓［ｉ］归并到Ｘ↓［ｊ］的概率ｐ↓［ｉｊ］（ｔ），其中，ｓ是蚂蚁Ｘ↓［ｓ］下一步可以选择的路径集合即Ｓ＝｛Ｘ↓［ｓ］｜ｄ↓［ｓｊ］≤ｒ，ｓ＝１，２，…，Ｎ｝；α，β为调节因子，起到防止所有蚂蚁沿相同路径得到相同结果所产生的停滞搜索的作用；　　　　４...

【技术特征摘要】

【专利技术属性】
技术研发人员：郭雷，杨卫莉，赵天云，肖谷初，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：87[中国|西安]