发明名称 基于免疫密母聚类的图像分割方法
摘要 本发明公开了一种基于免疫密母聚类的图像分割方法,主要解决已有聚类技术稳定性差、易陷入局部极值问题。其实现步骤为:1)输入一幅待分割图像;2)提取图像纹理特征;3)分水岭预分割图像,产生聚类数据;4)初始化第一个种群a;5)计算第一个种群a中个体适应度值,得到第二个种群b;6)对第一个种群a和第二个种群b同时进行进化操作,对第一个种群a依次进行交叉操作、变异操作、个体学以及选择操作,对第二个种群b依次进行克隆操作、变异操作及克隆选择操作;7)从精英种群c中取最大适应度值所对应的个体进行标记作为最终的图像分割结果。本发明具有稳定性高、分割结果区域一致性好和保留信息完整的优点,可用于SAR图像目标识别。
申请公布号 CN102663751A 申请公布日期 2012.09.12
申请号 CN201210099828.8 申请日期 2012.04.07
申请人 西安电子科技大学 发明人 马文萍;焦李成;李聪玲;黄媛媛;马晶晶
分类号 G06T7/00(2006.01)I;G06N3/12(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于免疫密母聚类的图像分割方法,包括如下步骤:(1)输入一幅待分割图像;(2)提取待分割图像的纹理特征,用灰度共生矩阵方法提取图像的前12维纹理特征向量,用小波分解方法提取图像的后10维纹理特征向量,得到N个D维的纹理特征向量,其中N为图像的像素个数,D为用灰度共生矩阵方法与小波分解方法提取待分割图像的纹理特征向量维数之和;(3)对待分割图像进行分水岭预分割,得到N1个不重叠的区域块,对每一个区域块所包含像素点的纹理特征向量取平均值,得到N1个D维的纹理特征向量,作为聚类的输入数据样本,N1为分水岭分割块数;(4)将第一个种群a分为两部分进行初始化:4a)根据最小生成树MST原理生成第一部分种群a1;4b)随机产生第二部分种群a2,每个个体基因位上的值是1到K之间的随机数,K为待分割图像的分类数;(5)对第一个种群a中的个体按照如下公式计算适应度值:<maths num="0001"><![CDATA[<math><mrow><mi>fitness</mi><mo>=</mo><mfrac><mn>1</mn><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>RN</mi><mi>i</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>PN</mi><mi>i</mi></msub></munderover><mi>dis</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>ijk</mi></msub><mo>,</mo><msub><mi>m</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><mn>1</mn></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>dis</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>ijk</mi></msub><mo>,</mo><msub><mi>m</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>g</mi><mo>=</mo><mn>1</mn></mrow><mi>D</mi></munderover><msup><mrow><mo>(</mo><msub><mi>p</mi><msub><mi>ijk</mi><mi>g</mi></msub></msub><mo>-</mo><msub><mi>m</mi><msub><mi>i</mi><mi>g</mi></msub></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>,</mo></mrow></math>]]></maths>RN<sub>i</sub>是第i类包含的区域块数,PN<sub>i</sub>是第i类中第j个区域块包含的像素点数,p<sub>ijk</sub>为第i类中第j区域块的第k个像素点,m<sub>i</sub>是第i类的聚类中心,<img file="FDA0000151123190000013.GIF" wi="73" he="48" />为第i类中第j区域块的第k个像素点的第g个特征,<img file="FDA0000151123190000014.GIF" wi="55" he="49" />为第i类聚类中心的第g个特征,D为用灰度共生矩阵方法与小波分解方法提取待分割图像的纹理特征维数之和;(6)计算出第一个种群a中所有个体的适应度值后,将适应度值由高到低排序,取前d个个体的聚类中心作为第二个种群b的抗体,将前d个个体对应的适应度值作为第二个种群b中抗体的亲和度值,d为第二个种群b中抗体的个数;(7)进化第一个种群a7a)用基于同类之间交叉方法对第一个种群a中的个体进行交叉操作,重组第一个种群a中的个体信息;7b)用单点变异方法对交叉操作后的第一个种群a中的个体进行变异操作,以提高第一个种群a的多样性;7c)用个体学习方法对变异操作后的第一个种群a中的个体进行局部搜索,使第一个种群a不会陷入局部最优解;7d)用精英联赛机制对个体学习操作后的第一个种群a中的个体进行选择操作,并把适应度值最高的个体保存到精英种群c中;(8)进化第二个种群b8a)对第二个种群b中的抗体进行比例克隆操作;8b)用均匀变异操作对比例克隆操作后的第二个种群b中的抗体进行变异操作,以提高第二个种群b的多样性;8c)用克隆选择方法对均匀变异操作后的第二个种群b中的抗体进行选择操作,并用选择操作之后第二个种群b中抗体取代选择操作后的第一个种群a中适应度值最低的d个个体;(9)如果进化种群a和b执行的次数达到G<sub>max</sub>次,则执行步骤(10),否则执行步骤(5),G<sub>max</sub>为进化两个种群a和b的最大迭代次数;(10)从精英种群c中,选择适应度最大值对应的个体作为最优个体;(11)在最优个体中,由每个块的类标找到每个像素点的类标,从闭区间[0,255]中任意选择一个整数作为所标记像素的灰度值,得到分割结果。
地址 710071 陕西省西安市太白南路2号