基于改进多目标和声搜索算法的交通图像分割方法

基于改进多目标和声搜索算法的交通图像分割方法
【技术领域】
[0001] 本发明属于图像处理技术领域，具体是一种基于改进多目标和声捜索算法的交通 图像分割方法。
【背景技术】
[0002] 图像分割结果的好坏直接地影响后续图像处理的质量，因此，图像分割是图像理 解领域中的重点和难点之一，该一艰难的任务给研究人员带来了巨大的挑战。图像分割的 目的是将一幅图像分割成若干个互不重叠的具有特定意义的区域，同一区域具有相似的特 性，不同区域差别较大。
[0003] 现有的图像分割方法主要可W划分为基于边缘的分割方法，基于区域的分割方 法，基于聚类的分割方法等。大多数聚类方法，其实质是目标函数的最优化问题。由于传统 的聚类方法存在易陷入局部最优、对初始聚类中屯、敏感等缺点，使得基于进化聚类的分割 方法成为当前图像分割方法中的热点，并广泛应用于各类图像的分割中。它很好的将进化 计算与聚类技术结合在一起，成功克服了传统聚类方法的不足，另外进化计算方法可W同 时优化多个目标，使得多目标进化聚类技术成为可能。鉴于进化计算的优点，研究合适的智 能优化算法，使其与聚类技术更好的结合，已成为一个备受关注的研究方向。
[0004] 和声捜索算法是用于解决全局最优问题的智能优化算法，具有控制参数少、简单 易懂、易于实现等优点，现已成功运用到各个领域。目前，基于多目标优化聚类的图像分割 方法大多只考虑两个分割标准，而每类图像都有各自的特点，因此存在应用面小，分割结果 不理想，分割正确率低等缺点。
[0005] 现有交通图像分割方法中，无论是多阔值分割方法，还是聚类分割方法，都需要人 为的设定分割目标的数目，而智能交通中基于视觉的智能车导航技术，需要计算机自动将 道路、车辆、障碍物等目标分离出来，显然传统的图像分割方法无法满足该一特定的要求。

【发明内容】

[0006] 本发明的目的在于克服上述已有方法的不足，提出了一种基于改进多目标和声捜 索算法的交通图像分割方法，本发明同时优化=个目标函数，并将动态聚类的算法融入到 多目标和声捜索算法中，提高了图像分割精度，实现了自动确定分割类别数，为后续交通图 像处理提供了较好的信息基础。
[0007] 基于改进多目标和声捜索算法的交通图像分割方法，包括如下步骤：
[000引 （1)输入待分割灰度图像I，统计图像I的灰度直方图GH=化1，1 = 0, 1，...，255}，hi为灰度级1的像素点出现的次数，设置迭代计数器t，初始值为1 ;
[0009] (2)初始化和声记忆库；
[0010] (2a)随机初始化和声记忆库歷中的和声Xi所表示的聚类中屯、的分割类别数K1， 1 = 1，...，歷5，歷8为和声记忆库大小，町￡怔。化，1^]，町为整数，1(。化为图像1最小分割 类别数，Km"为最大分割类别数；
[0011] (2b)根据图像I中像素灰度值范围，随机生成分割类别数为Ki的聚类中屯、Ci，。 ={Ck=化+帅-Lb)*r，k= 1，. ..，KJ，却为第k类的聚类中心U3、Ub分别是图像I中所 有像素值中的最小值、最大值，r是随机产生的一个0到1之间的实数；
[001引 （2c)将每个聚类中屯、。随机编排到对应和声X冲，Xi是一个KmJ隹向量；
[0013] (2d)找出歷中每个和声Xi中所有值不等于NOTcenter的分量作为该和声所表示 的聚类中屯、。，评估聚类中屯、。的S个适应度函数值f1，f2，f3，并将该S个适应度函数值作 为对应和声Xi的适应度值；
[0014] (3)生成一个新解：
[00巧](3a)创作的新和声X胃=(xr，xr，...，xr)中的每一个变量xr都是经过^下； 种机理产生，即保留和声记忆库中分量，音调微调，随机选择音调，n=Km"，该过程描述为：
[0016]
[0017] 其中ri、r2均为随机产生的0到1之间的实数，HMCR、PAR、bw分别为和声忆库保 留概率、音调调节概率、音调调节步长；
[0018] (3b)对新产生的和声XD6W执行变异操作，生成变异后的解向量X-ut;
[0019] (3c)对解向量X-ut执行随机操作，生成解向量Xand.
[0020] (3d)按照步骤（2d)的方式，评估新生成的解向量^^嘴立个适应度函数值fI，f2， fs;
[002U (4)重复步骤（3)，直到生成歷S个新解；
[002引 （5)将新生成的歷S个解与和声记忆库歷中的解组合在一起，形成一个具有 2X歷S个解的组合和声记忆库，并根据各个解的适应度函数值，采用NSGA-II中的快速非 支配排序策略，对组合和声记忆库中所有的解进行排序；
[002引 做删除歷中原有的和声，根据NSGA-II中拥挤比较操作，从组合和声记忆库中选 取最优的HMS个解保存到和声记忆库HM中，形成新的和声记忆库HM;
[0024] (7)判断当前迭代次数t是否大于最大迭代次数tm"，如果满足该条件，则执行步 骤巧），否则返回步骤（3)，t=t+1 ;
[0025] (8)从迭代结束后得到的和声记忆库中找出每个和声中所有值不等于NOTcenter 的分量作为该和声所表示的聚类中屯、，计算其聚类有效性指标PBM值，选取PBM值最大的对 应的聚类中屯、作为最优聚类中屯、，其中PBM表示为：
[00測其中El对于同一幅图像是一个常数，Dc= {maxd(c1，Cj),i,j= 1，. . .，K}是各聚 类中屯、间欧式距离最大值，K为分割类别数，d(x。Ck)为图像I中第i个像素点Xi到第k类 的聚类中屯、Ck的欧式距离；
[0029] (9)计算图像I各像素的灰度值与最优聚类中屯、的欧式距离，把该像素划分到距 离它的欧氏距离最小的聚类中屯、的类别中，得到每一个像素所属的类别，从而得到灰度图 像的分类结果。
[0030] 本发明的实现还在于：步骤（2a)中的所述的和声记库大小歷S在80至Ij120之间 取值效果较好。歷S取值过小，不易于保持和声记忆库中解的多样性，从而容易过早成熟， HMS取值过大会造成计算量大，导致计算时间过长，HMS的取值在一定程度上影响着算法的 收敛速度和收敛精度，通过大量实验仿真得出歷S在80至Ij120之间取值，算法能够同时获 得较好的收敛速度和收敛精度。
[0031] 本发明的实现还在于：步骤（2c)中每个聚类中屯、。随机编排到对应和声Xi中，其 过程包括；初始化一个Kmax维的空向量V1=(，，...，），将聚类中屯、C海类的中屯、Ck，k= 1，..，Ki，随机放到Vi中任意一空维，直到聚类中屯、Ci中的所有类中屯、被放完，对于Vi中仍 存在的空维用NOTcenter替代，NOTcenter是一个负常数，表示不是类中屯、，最后将Vi赋给 Xi〇
[0032] 本发明的实现还在于：步骤（3a)中的音调微调采用一种随迭代次数t增加而逐渐 减小的音调调节步长bw，表不为；
[0033]
[0034] 其中bWm"，bWmi。为音调调节步长的上下界，9为一个正常数。
[0035] 在本发明中，采用动态变化的bw使算法在迭代早期有较好的全局捜索能力，而迭 代后期有较好的局部捜索能力。
[0036] 本发明的实现还在于：步骤（3b)中对新产生的和声XD6W执行变异操作，其过程包 括：
[0037] (3b. 1)对和声xnew中变量坤"'进行变异操作，j=1，. .，n，若X胃的值为 NOTcenter，则xf''保持不变，即=x^e'v，否则变异，变异可表示为；
[00%]
[0039] 其中r是随机产生的一个0到1之间的实数，Pmut为变异概率，W〇.Ub~U))是均 kx 值为0,标准差为的正态分布随机产生的一个实数，kx是一个正系数；kx
[0040] (3b.。若变异后变量xf的值超出范围[化化]，则将与xf"接近的边界值赋给 X?*;
[0041] (3b. 3)对和声xnew中每个变量依次按步骤（3b. 1) (3b. 2)执行，得到经过变异操作 后的解向量X?t。
[0042] 在本发明中，对新产生的和声中的每 个变量依次执行变异操作，能够增强算法逃 避局部最优的能力，同时使变异后产生的解所表示的聚类中屯、的分割类别数保持动态变 化，有利于捜寻到最佳分割类别数。
[0043] 本发明的实现还在于：步骤（3c)中对解向量rut进行随机操作，操作过程包括； 统计解向量X-ut中值不为NOTcenter的分量的个数K，K表示分割类别数，若K>Kmi。，X-ut保 持不变，即r^=xmut，否则随机生成分割类别数K，KG[Kmi"，KmJ，并在[化，化]范围内随 机生成分割类别数为K的聚类中屯、C，按步骤（2c)中所述的随机编码方式将C随机编排到 Xrand中。
[0044] 在本发明中，对变异后的解执行随机操作，有利于确保最后产生的解所表示的聚 类中屯、的分割类别数大于等于Kmi。，同时能够增强解的多样性。
[0045] 本发明的实现还在于：步骤（3d)中所述的S个目标函数值表示为：
[0049] 其中K为分割类别数，N为图像I像素点总数，Uki为图像I中第i个像素点X i属 于第k类的模糊隶属度，屯1。={mind(c。Cj),i,j= 1，. . .，K,i声j}是各聚类中屯、间欧式 距离最小值，d(X。Ck)为图像I中第i个像素点X剧第k类的聚类中屯、ck的欧式距离，up。 是第q类的聚类中屯、c。属于第p类的聚类中屯、cp的模糊隶属度，uki，yp。分别可表示为：
[0化2] 在本发明中，选用S个目标函数来评估分割结果的质量，第一目标函数代表聚类 紧凑性，第二个目标函数代表致密分离性，第=个目标函数代表模糊分离性，该=个评价指 标有利于算法朝着最优聚类中屯、方向进行捜索。
[0化3] 本发明与现有技术相比具有W下优点：
[0054] 1.本发明采用了S个互补的目标函数，增加了评价的多样性，提高了图像分割的 准确率，具有更好的抗噪能力。
[0055] 2.本发明将动态聚类的思想融入到多目标和声捜索算法中，动态变化分割类别数 有利于最佳分割类别数的捜寻，实现了自动确定分割类别数，解决了需要人为确定分割类 别数该一难题。
[0056] 3.本发明将动态变化的音调调节步长、变异操作、随机操作引入到多目标和声捜 索算法中，提高了算法的全局和局部寻优能力，增强了算法逃离局部最优的能力，提高了算 法的收敛速度和收敛精度，为最优聚类中屯、的捜索提供了有力保障。
【附图说明】
[0化7] 图1是本发明的实现流程图；
[005引图2中的（a)、化）、（C)和（d)是本发明与现有技术在一幅有道路、多车辆的交通 图像上的分割实验结果对比图；
[0059]图3中的（a)、化）、（C)和（d)是本发明与现有技术在一幅有道路、车辆、障碍物的 交通图像上的分割实验结果对比图；
[0060]图4中的（a)、化）、（C)和（d)是本发明与现有技术在一幅道路交通图像上的分割 实验结果对比图；
[0061]图5中的（a)、化）、（C)和（d)是本发明与现有技术在一幅有道路、单车辆的交通 图像上的分割实验结果对比图；
[006引图6中的（a)、化）、（C)、（d)和（e)是本发明与现有技术在一幅高斯噪声为 (0, 0.01)的交通标志图像上的分割实验结果对比图；
[0063]图7中的（a)、化）、（C)和（d)是本发明与现有技术在一幅高斯噪声为化0. 015) 的交通标志图像上的分割实验结果对比图。
【具体实施方式】
[0064] 下面结合附图对本发明做进一步详细描述。
[0065] 实施例1
[0066] 本发明提出了一种基于改进多目标和声捜索算法的交通图像分割方法，属于图像 处理技术领域，主要设及到交通图像分割技术。本例仿真是在主频2. 00G监的Intel(时 Core(TM)2DuoCPUT5870、内存1.96GB的硬件环境和MATLABR2010a的软件环境下进行的。
[0067] 参照图1，本发明的具体实现步骤如下：
[0068] 步骤1，输入一幅待分割的灰度图像I，该灰度图像用采样像素8位的非线性尺度 保存，有256级灰度。
[0069] 在本实施例中，输入一幅有道路、多车辆的交通图像Ii，参见图2(a)，其大小为 250X214,统计图像I的灰度直方图GH=化1，1 = 0, 1，...，255}，hi为灰度级1的像素点 出现的次数，设置迭代计数器t，初始值为1。
[0070] 步骤2,初始化和声记忆库；
[007U (2a)随机初始化和声记忆库歷中的和声Xi所表示的聚类中屯、的分割类别数K1， 1 = 1，...，歷5，歷8为和声记忆库大小，町￡怔。化，1^]，町为整数，1(。化为图像1最小分割 类别数，Km"为最大分割类别数。Kmi。通常为2，Km"的取值需要根据待分割图像中可能存在 的最大分割类别数来确定，对于交通图像而言，其最大分割类别数通常小于8 ;
[0072] 姊）根据图像I中像素灰度值范围，随机生成分割类别数为Ki的聚类中屯、C1，。 ={Ck=化+帅-Lb)*r，k= 1,. . . ,K;}，Ck为第k类的聚类中心化、Ub分别是图像I中 所有像素值中的最小值、最大值，r是随机产生的一个0到1之间的实数；
[007引 （2c)将每个聚类中屯、。随机编排到对应和声X冲，Xi是一个Km。,维向量，其过 程包括：初始化一个Km"维的空向量Vi=(，，...，），将聚类中屯、Ci每类的中屯、Ck，k= 1，..，Ki，随机放到Vi中任意一空维，直到聚类中屯、Ci中的所有类中屯、被放完，对于Vi中仍 存在的空维用NOTcenter替代，NOTcenter是一个负常数，表示不是类中屯、，最后将Vi赋给 Xi0
[0074] (2d)找出歷中每个和声Xi中所有值不等于NOTcenter的分量作为该和声所表示 的聚类中屯、评估聚类中屯、。的S个适应度函数值fs，并将该S个适应度函数值作 为对应和声Xi的适应度值，=个目标函数值表示为；
[007引其中K为分割类别数，N为图像I像素点总数，咕为图像I中第i个像素点Xi属 于第k类的模糊隶属度，屯1。={mind(c。Cj),i,j= 1，. . .，K,i声j}是各聚类中屯、间欧式 距离最小值，d(X。Ck)为图像I中第i个像素点X剧第k类的聚类中屯、ck的欧式距离，up。 是第q类的聚类中屯、c。属于第p类的聚类中屯、cp的模糊隶属度，uki，yp。分别可表示为：
[00川在本实施例中，HMS取值100,Kmi。取值2,Km"取值8,NOTcenter设置为-100。[0082]步骤3,生成一个新解：
[008引（3a)创作的新和声X胃=(xr，xr，...，xr)中的每一个变量x:?都是经过^下； 种机理产生，即保留和声记忆库中分量，音调微调，随机选择音调，n=Km"，该过程描述为：
[0084]
[00财其中均为随机产生的0到1之间的实数，HMCR、PAR、bw分别为和声忆库保 留概率、音调调节概率、音调调节步长，其中bw可表示为：
[0086]
[0087] 其中bWm"，bWmi。为音调调节步长的上下界，CP为一个正常数。9取较小值时，算法 收敛速度慢，但更容易逃离局部最优，9取较大值时，收敛速度快，但易陷入局部最优，权衡 算法收敛速度和逃离局部最优的能力，取值2。
[008引 （3b)对新产生的和声XD6W执行变异操作，生成变异后的解向量Xuut，其过程包括；
[0089] (3b. 1)对和声xnew中变量xfw执行变异操作，j二1，. .，n，若xfw的值为 NOTcenter，则X；*?保持不变，即X;"" = ，否则变异，变异可表示为；
[0090]
[OOW] 其中r是随机产生的一个0到1之间的实数，Pmut为变异概率，N(0，Ub-促均 kx IIK-Th 值为0,标准差为的正态分布随机产生的一个实数，kx是一个正系数；kx
[0092] (3b. 2)若变异后变量xijD"的值超出范围[化，化]，则将与X；*"接近的边界值赋给 X?* ；
[0093] (3b. 3)对和声xnew中每个变量依次按步骤（3b. 1 ) (3b. 2)执行，得到经过变异操作 后的解向量x?t;
[0094] (3c)对解向量r"销行随机操作，生成解向量Xand,其过程包括操作过程包括；统 计解向量rut中值不为NOTcenter的分量的个数K，K表示分割类别数，若K>Kmi。，X-ut保 持不变，即r^=xmut，否则随机生成分割类别数K，KG怔并在[化，化]范围内随 机生成分割类别数为K的聚类中屯、C，按步骤（2c)中所述的随机编码方式将C随机编码到 Xrand中；
[009引 （3d)按照步骤（2d)的方式，评估新生成的解向量^^的立个适应度函数值fI，f2， fg。
[0096]在本实施例中，HMCR=0. 95,PAR=0. 45,bWm化=1，bWmax= 0.IX(化-化），pmut =1/Kmax，kx= 20。
[0097] 步骤4,重复按照步骤3,直到生成歷S个新解。
[009引步骤5,将新生成的歷S个解与和声记忆库歷中的解组合在一起，形成一个具有 2X歷S个解的组合和声记忆库，并根据各个解的适应度函数值，采用NSGA-II中的快速非 支配排序策略，对组合和声记忆库中所有的解进行排序。
[0099] 步骤6,删除歷中原有的和声，根据NSGA-II中拥挤比较操作，从组合和声记忆库 中选取最优的歷S个解保存到和声记忆库歷中，形成新的和声记忆库歷。
[0100] 步骤7,判断当前迭代次数t是否大于最大迭代次数tm"，如果满足该条件，则执行 步骤8,否则返回步骤3,t=t+1。
[0101] 在本实施例中，最大迭代次数tm"= 150。
[0102] 步骤8,从迭代结束后得到的和声记忆库中找出每个和声中所有值不等于 NOTcenter的分量作为该和声所表示的聚类中屯、，计算其聚类有效性指标PBM值，选取PBM 值最大的对应的聚类中屯、作为最优聚类中屯、，其中PBM表示为：
[010引其中El对于同一幅图像是一个常数，Dc= {maxd(c1，Cj),i,j= 1，. . .，K}是各聚 类中屯、间欧式距离最大值，K为分割类别数，d(x。Ck)为图像I中第i个像素点Xi到第k类 的聚类中屯、Ck的欧式距离。
[0106] 步骤9,计算图像I各像素的灰度值与最优聚类中屯、的欧式距离，把该像素划分到 距离它的欧氏距离最小的聚类中屯、的类别中，得到每一个像素所属的类别，从而得到灰度 图像的分类结果。还原的分类结果参见图2化）。
[0107] 本发明的对比实验为K-means、模糊C-均值图像分割方法，由于K-means、模糊 C-均值图像分割方法都需要人为设定分割类别数，为便于比较，将其分割类别数均设置为 本发明方法最后分割结果所得的分割类别数。图2(a)为交通图像Ii的原始灰度图像，图 2化）、图2 (C)、图2 (d)分别为采用本发明方法、K-means、模糊C-均值方法得到的分割结 果图，从图2(b)可W看出，本发明方法自动将交通图像Ii划分成了立类，车辆为一类，道路 为一类，道路上的分道线为一类，该发明方法很好的将车道和车辆分割出来。从图2化）、图 2 (C)、图2 (d)可W看出，本发明很好的将车牌上的数字及字母分割出来了，而K-means、模 糊C-均值分割方法均存在分割不完整情况，同时，本发明克服了该两种方法均需人为设置 分割类别数的缺点。
[010引实施例2
[0109] 基于改进多目标和声捜索算法的交通图像分割方法同实施例1，在实施例2中， 对比实验为K-means、模糊C-均值图像分割方法，K-means、模糊C-均值图像分割方法的 分割类别数K均设置为3。实施例2中输入有道路、车辆、障碍物的交通图像I,，其大小为 276X223像素，灰度级为256。义用本发明得到最终分类结果图参见3化）。
[0110] 图3(a)为交通图像12的原始灰度图像，图3化）、图3(c)、图3(d)分别为采用本 发明方法、K-means、模糊C-均值方法得到的分割结果图。从图3(b)可W看出，本发明自动 将图像12划分成=类，车辆及障碍物为一类，道路为一类，道路两边的树木为一类，很好的 将障碍物、车辆从道路上分割出来。从图3 (a)、图3化）、图3 (C)、图3 (d)可W看出，本发明 方法有更好的抗噪能力，原灰度图3(a)中路面灰尘较多，K-means、模糊C-均值方法均存在 较多的将该些灰尘误分为障碍物，而本发明方法明显地减少了误分区域，具有更高的分割 准确率。
[0111] 实施例3
[0112] 基于改进多目标和声捜索算法的交通图像分割方法同实施例1，在实施例3中，对 比实验为K-means、模糊C-均值图像分割方法，K-means、模糊C-均值图像分割方法的分割 类别数K均设置为3。实施例3中输入有道路的交通图像l3,其大小为318X245像素，灰 度级为256。采用本发明得到最终分类结果图参见4化）。
[0113] 图4(a)为交通图像Is的原始灰度图像，图4化）、图4(c)、图4(d)分别为采用本 发明方法、K-means、模糊C-均值方法得到的分割结果图。从图4(b)可W看出，本发明自 动将图像13划分成=类，道路为一类，道路边线及分道线为一类，道路两旁的树木为一类， 很好的将道路从图像中分割出来。从图4化）、图4(c)、图4(d)可W看出，本发明得到了与 K-means、模糊C-均值图像分割方法相似的结果。
[0114] 实施例4
[0115] 基于改进多目标和声捜索算法的交通图像分割方法同实施例1，在实施例4中，对 比实验为K-means、模糊C-均值图像分割方法，K-means、模糊C-均值图像分割方法的分割 类别数K均设置为2。实施例4中输入有道路、单车辆的交通图像14,其大小为278X244像 素，灰度级为256。采用本发明得到最终分类结果图参见5化）。
[0116] 图5(a)为交通图像L的原始灰度图像，图5化）、图5(c)、图5(d)分别为采用本发 明方法、K-means、模糊C-均值方法得到的分割结果图。从图4(b)可W看出，本发明自动 将图像L划分成二类，车辆为一类，分道线为一类，很好的将车辆从图像中分割出来。从图 5化）、图5 (C)、图5 (d)可W看出，本发明得到了与K-means、模糊C-均值图像分割方法相似 的结果。
[0117] 实施例5
[0118] 基于改进多目标和声捜索算法的交通图像分割方法同实施例1，在实施例5中，对 比实验为K-means、模糊C-均值图像分割方法，K-means、模糊C-均值图像分割方法的分割 类别数K均设置为2。实施例5中输入高斯噪声为（0, 0.01)的交通标志图像Ig，其大小为 499X495像素，灰度级为256。采用本发明得到最终分类结果图参见6(c)。
[0119] 图6(a)为无噪声交通标志图像Is的原始灰度图像，图6化）是高斯噪声为 (0, 0.01)的交通标志图像le的原始灰度图像，图6(c)、图6(d)、图6(e)分别为采用本发明 方法、K-means、模糊C-均值方法得到的分割结果图。从图6(c)可W看出，本发明自动将图 像I。划分成二类，圆中的矩形为一类，圆中除去矩形的部分为一类，符合分割要求，体现了 本发明方法良好的抗噪能力。从图6(c)、图6(d)、图6(e)及表1可W看出，本发明方法误 分像素的个数要少于其他两种方法，具有更强的抗噪能力及更高分割准确率。
[0120] 表种算法性能比较
[0121]
[0122] 实施例6
[0123] 基于改进多目标和声捜索算法的交通图像分割方法同实施例1，在实施例6中，对 比实验为K-means、模糊C-均值图像分割方法，K-means、模糊C-均值图像分割方法的分割 类别数K均设置为2。实施例6中输入高斯噪声为化0.015)的交通标志图像le，其大小 为499X495像素，灰度级为256。采用本发明得到最终分类结果图参见7化）。
[0124] 图7(a)为高斯噪声为化0.015)交通标志图像le的原始灰度图像，图7化）、图 7(c)、图7(d)分别为采用本发明方法、K-means、模糊C-均值方法得到的分割结果图。从图 7(b)可W看出，本发明自动将图像le划分成二类，圆中的 矩形为一类，圆中除去矩形的部分 为一类，符合分割要求，体现了本发明方法良好的抗噪能力。从图7化）、图7 (C)、图7 (d)及 表1可W看出，本发明方法误分像素的个数要少于其他两种方法，具有更强的抗噪能力更 高分割准确率。
[0125] 综上，本发明提出的基于改进多目标和声捜索算法的交通图像分割方法，主要解 决交通图像需要人为确定分割类别数的难题，同时克服了传统聚类算法分割中存在的抗噪 能力弱、对初始聚类中屯、敏感等问题。本发明方法针对灰度图像像素值进行处理，根据灰度 范围来初始化和声记忆库，并将动态聚类的思想融入到多目标和声捜索算法中；本发明还 在于选取=个互补的目标函数同时优化，克服了评价的单一性；本发明的关键在于在多目 标和声捜索算法中，引入动态音调调节步长、变异操作、随机操作，增强了算法的全局和局 部捜索能力，全面提升了算法的收敛速度和收敛精度，为最优聚类中屯、的捜寻提供了有力 保障。
[01%] 本发明能够自动确定交通图像中的分割类别数，分割精度高，抗噪能力强，为后续 交通图像处理工作提供良好的信息基础。
【主权项】
1. 一种基于改进多目标和声搜索算法的交通图像分割方法，其特征在于包括有如下步 骤： (1) 输入待分割灰度图像I，统计图像I的灰度直方图GH = Ih1, 1 = 0, 1，...，255} 为灰度级1的像素点出现的次数，设置迭代计数器t，初始值为1 ; (2) 初始化和声记忆库： (2a)随机初始化和声记忆库HM中的和声Xi所表示的聚类中心的分割类别数Ki, i = 1，...，HMS，HMS为和声记忆库大小，Ki e [K min，Kmax]，&为整数，K min为图像I最小分割类别 数，Kmax为最大分割类别数； (2b)根据图像I中像素灰度值范围，随机生成分割类别数为Ki的聚类中心C i，Ci= {c k =Lb+(Ub-Lb)*r，k = 1，. . .，Kj，ck为第k类的聚类中心，Lb、Ub分别是图像I中所有像 素值中的最小值、最大值，r是随机产生的一个0到1之间的实数； (2c)将每个聚类中心Ci随机编排到对应和声X i中，X i是一个K max维向量； (2d)找出HM中每个和声Xi*所有值不等于NOTcenter的分量作为该和声所表示的聚 类中心Ci，评估聚类中心Ci的三个适应度函数值f i，f2，f3，并将这三个适应度函数值作为对 应和声Xi的适应度值； (3) 生成一个新解： (3a)即兴创作一个和声又_=(\1~4广，...4，），新创作的和声中的每一个变量乂广 都是经过以下三种机理产生，即保留和声记忆库中分量，音调微调，随机选择音调，n = Kmax， 该过程描述为：其中A、1*2均为随机产生的0到1之间的实数，HMCR、PAR、bw分别为和声忆库保留概 率、音调调节概率、音调调节步长； (3b)对新产生的和声Xn?执行变异操作，生成变异后的解向量Xmut; (3c)对解向量Xmut执行随机操作，生成解向量X ?d; (3d)按照步骤（2d)的方式，评估新生成的解向量Xrand的三个适应度函数值f Pfyf3; (4) 重复步骤（3)，直到生成HMS个新解； (5) 将新生成的HMS个解与和声记忆库HM中的解组合在一起，形成一个具有2 X HMS个 解的组合和声记忆库，并根据各个解的适应度函数值，采用NSGA-II中的快速非支配排序 策略，对组合和声记忆库中所有的解进行排序； (6) 删除HM中原有的和声，根据NSGA-II中拥挤比较操作，从组合和声记忆库中选取最 优的HMS个解保存到和声记忆库HM中，形成新的和声记忆库HM ; (7) 判断当前迭代次数t是否大于最大迭代次数tmax，如果满足该条件，则执行步骤 (8)，否则返回步骤（3)，t = t+l; (8) 从迭代结束后得到的和声记忆库中找出每个和声中所有值不等于NOTcenter的分 量作为该和声所表示的聚类中心，计算其聚类有效性指标PBM值，选取PBM值最大的对应的 聚类中心作为最优聚类中心； (9)计算图像I各像素的灰度值与最优聚类中心的欧式距离，把该像素划分到距离它 的欧氏距离最小的聚类中心的类别中，得到每一个像素所属的类别，从而得到灰度图像的 分类结果。2. 根据权利要求1所述的基于改进多目标和声搜索算法的交通图像分割方法，其特征 在于：步骤（2a)中的所述的和声记库大小HMS在80到120之间取值。3. 根据权利要求1所述的基于改进多目标和声搜索算法的交通图像分割方法，其特征 在于：步骤（2c)中每个聚类中心Ci随机编排到对应和声X，其过程包括：初始化一个K max 维的空向量Vi=(，，...，），将聚类中心C i每类的中心c k，k = 1，..，Ki,随机放到Vi中任意 一空维，直到聚类中心(^中的所有类中心被放完，对于V i中仍存在的空维用NOTcenter替 代，NOTcenter是一个负常数，表示不是类中心，最后将Vi赋给X i。4. 根据权利要求1所述的基于改进多目标和声搜索算法的交通图像分割方法，其特征 在于：步骤（3a)中的音调微调采用一种随迭代次数t增加而逐渐减小的音调调节步长bw， 表示为："imx 其中bwmax，bwmin为音调调节步长的上下界，φ为一个正常数。5. 根据权利要求1所述的基于改进多目标和声搜索算法的交通图像分割方法，其特征 在于：步骤（3b)中对新产生的和声Χη?执行变异操作，其过程包括： (3b. 1)对和声Xnew中变量行变异操作，j = 1，· ·，n，若χΓ的值为NOTcenter，则 保持不变，即χΓ =ΧΓ，否则变异，变异可表示为：其中r是随机产生的一个0到1之间的实数，pmut为变异概率，N(0.Ub^b)是均值为 KX I Jh - I b 〇,标准差为一^的正态分布随机产生的一个实数，kx是一个正系数； KX (3b.2)若变异后变量χΓ的值超出范围[Lb，Ub]，则将与Xmt接近的边界值赋给ΧΓ"; (3b. 3)对和声Xn?中每个变量依次按步骤（3b. I) (3b. 2)执行，得到经过变异操作后的 解向量Xmut。6. 根据权利要求1所述的基于改进多目标和声搜索算法的交通图像分割方法，其特征 在于：步骤（3c)中对解向量Xmut进行随机操作，操作过程包括：统计解向量X mut中值不为 NOTcenter的分量的个数K，K表示分割类别数，若K彡Kmin，X mut保持不变，即X _= X mut，否 则随机生成分割类别数K，Ke [Kmin，Kmax]，并在[Lb，Ub]范围内随机生成分割类别数为K的 聚类中心C，按步骤（2c)中所述的随机编码方式将C随机编排到X rand中。7. 根据权利要求1所述的基于改进多目标和声搜索算法的交通图像分割方法，其特征 在于：步骤（3d)中所述的三个目标函数值表示为：其中K为分割类别数，N为图像I像素点总数，Uki为图像I中第i个像素点X i属于第 k类的模糊隶属度，dmin= {min d(c i, Cj), i, j = 1，. . .，K, i辛j}是各聚类中心间欧式距离 最小值，d(Xi，ck)为图像I中第i个像素点X i到第k类的聚类中心c k的欧式距离，u M是第 q类的聚类中心c,属于第p类的聚类中心c 5的模糊隶属度，u ki，μ M分别可表示为：
【专利摘要】本发明公开了基于改进多目标和声搜索算法的交通图像分割方法，主要解决现有交通图像分割技术中评价指标单一、抗噪声能力弱、分割准确率不高、需人为确定图像分割类别数的问题。其实现步骤主要包括：读入灰度图像，统计灰度直方图；初始化和声记忆库；产生一群新解，个数等于和声记忆库大小；更新和声记忆库；判断是否满足终止条件；根据PBM评价指标，从和声记忆库中选取最优聚类中心；根据最优聚类中心，对图像中各像素分类，得到分割结果。本发明与现有交通图像分割技术相比，评价指标多元化，分割准确率高，抗噪声能力强，能自动确定图像中的分割类别数，可用于对交通图像的分割。
【IPC分类】G06T7/00
【公开号】CN104881871
【申请号】CN201510266684
【发明人】袁小芳, 戴香山, 向永忠, 王耀南
【申请人】湖南大学, 袁小芳
【公开日】2015年9月2日
【申请日】2015年5月22日