一种基于磷虾优化算法的摇床矿带分带图像分割法
【专利摘要】本发明涉及基于磷虾优化算法的摇床矿带分带图像分割法，属图像处理领域。首先将彩色的摇床矿带分带图像转换为灰度图像；再对磷虾进行初始化，后计算磷虾的适应度函数值，再对磷虾进行排列，随着迭代次数的变化对惯性权重进行更新，计算磷虾的运动向量，觅食运动向量和物理扩散运动向量，更新磷虾的位置，计算磷虾个体的适应度函数值，在满足一定的迭代次数后，通过适应度函数的最大化，找出适应度函数值最优的磷虾，该磷虾对应的位置即为摇床矿带分带图像的最佳阈值，根据最佳阈值对摇床矿带分带图像进行分割。本发明随着迭代次数的变化对惯性权重进行更新，算法搜索到最佳阈值的速度更快和精度更高，适合于摇床矿带分带图像的分割。
【专利说明】
一种基于磷虾优化算法的摇床矿带分带图像分割法
技术领域
[0001] 本发明涉及一种基于磷虾优化算法的摇床矿带分带图像分割法，属于图像处理技 术领域。
【背景技术】
[0002] 目前，在我国摇床选矿中，大部分都使用肉眼来观察摇床矿带分带的变化情况和 人工的方法来分割摇床矿带分带，这种方法存在浪费劳力、实时性差、分带精度低和金属回 收率低等缺点。为了克服其缺点，提出了基于数字图像处理技术的摇床矿带分带图像分割， 其不需要人工干预，直接采用图像分割技术对摇床矿带分带图像进行实时分割，其中图像 分割法是关键技术。目前图像分割方法有很多种，考虑到摇床矿带分带图像需要实时分割 和摇床床面上形成的精矿、次精矿、中矿和尾矿的分带在颜色和灰度上存在一定的不同，因 此采用阈值分割算法对摇床矿带分带图像进行分割。
[0003] 传统的阈值分割方法，对单阈值分割非常有效，但是对于多阈值分割存在计算复 杂度高和计算时间长的缺点，为了解决此问题，提出基于智能群体优化算法的阈值分割法
[0004] 磷虾优化算法是一种新的智能群体优化算法，相关的文献中提出其性能优于其它 的许多优化算法。但基本的磷虾优化算法仍然存在容易陷入局部最优，收敛速度慢的问题。
[0005] 因此，需要对基本的磷虾优化算法进行改进，并将其应用于摇床矿带分带图像的 分割。
【发明内容】
[0006] 本发明提供了一种基于磷虾优化算法的摇床矿带分带图像分割法，目的是克服肉 眼观察摇床床面上形成的矿带分带图像和人工分割摇床矿带分带存在的缺点，提出一种新 的基于磷虾优化算法的摇床矿带分带图像分割法，针对基本的磷虾优化算法应用于摇床矿 带分带图像分割时产生的收敛速度慢和计算时间长问题，提出一种改进的磷虾优化算法， 用于摇床矿带分带图像分割，提高收敛速度和收敛精度。
[0007] 为了达到上述目的，本发明首先对彩色的摇床矿带分带图像进行灰度化处理，将 彩色的图像转换为灰度图像;对磷虾位置进行初始化，让磷虾随机地分布在整个搜索空间 中，以Kapur熵作为适应度函数，计算所有磷虾的适应度函数值，并按其对磷虾进行排序，每 个磷奸个体计算自己的运动向量、觅食运动向量和物理扩散运动向量，对自己的位置进行 更新。通过多次迭代，在满足最大迭代次数时，找出适应度函数值最优即全局最优的磷虾个 体，磷虾所处的位置为摇床矿带分带图像的最佳阈值，最后利用最佳阈值对矿带分带进行 分割。该算法在搜索最佳阈值的过程中，根据迭代次数对惯性权重《"进行更新，使磷虾算 法的收敛速度更快和精度更高。
[0008] 改进后的磷虾优化算法惯性权重ω n更新方式如下所示：
[0009]
[0010]其中ωη为惯性权重，t为迭代次数，iter_max为最大迭代次数，CO max为惯性权重的 最大值，《min为惯性权重的最小值；
[0011] 实现本发明的技术方案的具体方法和步骤如下：
[0012] (1)在摇床选矿的过程中，实时拍摄摇床床面上形成的矿带分带图像，然后将图像 输入计算机中，对图像进行灰度化处理，将彩色的摇床矿带分带图像转换为灰度图像；
[0013] ⑵磷虾个体的初始化:设置磷虾算法的参数，觅食速度Vf，最大扩散速度D max，最 大移动速度Nmax，最大迭代次数iterjnax和磷虾的个数M;产生M个磷虾个体均匀地分布在摇 床矿带分带图像[0，255 ]的灰度直方图空间中；
[0014] (3)计算磷虾的适应度函数值，并根据适应度函数值，对磷虾个体进行从最好到最 差的排列；
[0015] (4)利用公： ?惯性权重《"进行更新，式中：ωη为惯性 权重，t为迭代次数，iter_max为最大迭代次数，ω max为惯性权重的最大值，ω min为惯性权重 的最小值；
[0016] (5)对每个磷虾个体i，执行下列步骤：
[0017] a)计算磷虾个体i的运动向量iV产+<￥) + ftvV,°w，式中：化为运动向 量，Nmax为最大移动速度为磷虾个体i的受周围邻居吸引产生的运动向量，a,为当 前最优磷虾个体吸引磷虾个体i产生的运动向量，ω η为惯性权重，汉f为磷虾个体i的上一 次运动向量；
[0018] b)计算觅食运动向量巧=心爲+，其中：汊=汊为食物的吸 引力，爲tenS到目前为止适应度函数值最大的磷虾个体i的影响力，F 1为磷虾个体i的觅食 运动向量;Vf为觅食速度，觅食运动向量的惯性权重，为磷虾个体i的上一次觅食 运动向量；
[0019] c)计算物理扩散运动向量，式中=D1为物理扩散运动向 量，Dmax为最大扩散速度，t为当前的迭代次数，iterjnax为最大迭代次数，δ为随机方向向 量； Jy
[0020] d)对磷虾的位置进行更新，更新公式为：? + Δ/) = Α^) + Δ(-γ，式中ddt+A at t)为磷虾个体i在t+Δ t时刻的位置，Xi(t)为磷虾个体i在t时刻的位置，△ t为时间间隔；
[0021] (6)计算磷虾个体的适应度函数值，找出当前适应度函数值最优的磷虾个体；
[0022] (7)不断循环(4)、（5)、（6)三个步骤，直到满足最大迭代次数tlter_ max，通过适应度 函数的最大化，找出经过多次迭代后适应度函数值最优的磷虾个体，该磷虾所处的位置即 为摇床矿带分带图像的最佳阈值，最后根据求出的最佳阈值对摇床矿带分带图像进行分 割。
[0023]本发明中所述步骤(2)中的参数，最大迭代次数i ter_max范围为[30，50 ]，磷4下的 个数M的范围为[25,50]。
[0024]本发明与现有技术相比具有下列优点：
[0025] 1、本发明提出基于磷虾优化算法的摇床矿带分带图像分割法，可以有效地克服肉 眼观察摇床矿带分带情况和人工分割摇床矿带分带存在的问题，节约劳动力，提高摇床矿 带分带分割的实时性，提高矿物的回收率和利用率；
[0026] 2、本发明提出的磷虾优化算法，提出随着迭代次数的变化对惯性权重进行更新的 方法，提高算法搜索全局最优值即最佳阈值的速度和精度，算法能够快速和准确地找到全 局最优即摇床矿带分带图像的最佳阈值。
【附图说明】
[0027] 图1为本发明基于磷虾优化算法的摇床矿带图像分割法的流程图。
【具体实施方式】
[0028] 实施例1:参见图1，以从云南锡业集团大屯选矿厂拍摄的摇床床面上形成的锡矿 矿带分带图像为例，利用VC++软件编写程序，对锡矿的精矿、中矿和尾矿矿带进行分割，采 用的方法和具体步骤如下：
[0029] (1)对彩色的摇床矿带分带图像进行灰度化处理，将彩色的摇床矿带分带图像转 换为灰度图像；
[0030] (2)磷虾个体的初始化:对参数进行设置，觅食速度Vf = 0.02,最大扩散速度Dmax = 0.005，最大移动速度Nmax = 0.01，最大迭代次数iter_max = 30和磷虾的个数M= 25;最大值 惯性权重ω max = 〇 . 9，惯性权重的最小值ω min = 〇. 1;产生M个磷虾个体均匀地分布在摇床矿 带分带图像[0,255]的灰度直方图空间中；
[0031 ] ⑶利用Kapur熵的公式H( ti，t2) =Ηι+Η2+Η3计算磷奸的适应度函数值，并根据适应 度函数值，对磷虾个体进行最好到最差的排列，25个磷虾个体的适应度函数值为：
[0032] H= [12.3020060145,12.2795029028,12.0847884278,11.6167119173， 11.5683895895,11.0988169532,10.8340270715,10.7686284683,10.5340270715,0,0,0,
[0033] (4)利用公?
f惯性权重《"进行更新，式中：ωη为惯性 权重，t为迭代次数，iter_max为最大迭代次数，ω max为惯性权重的最大值，ω min为惯性权重 的最小值；
[0034] (5)对每个磷虾个体i，执行下列步骤：
[0035] a)计算磷虾个体i的运动向量= + <_，+气汉f，式中：化为运动向 量，Nmax为最大移动速度，为磷虾个体i的受周围邻居吸引产生的运动向量，a,为当 前最优磷虾个体吸引磷虾个体i产生的运动向量，ω η为惯性权重，为磷虾个体i的上一 次运动向量；
[0036] b)计算觅食运动向量朽=匕病+巧/^，其中：爲=爲~ +及^,/?/°°(/为食物的吸 引力，#产"为到目前为止适应度函数值最大的磷虾个体i的影响力，F 1为磷虾个体i的觅食 运动向量;Vf为觅食速度，COf为觅食运动向量的惯性权重，F '为磷奸个体i的上一次觅食 运动向量；
[0037] c)计算物理扩散运动向1
式中：Di为物理扩散运动向 量，Dmax为最大扩散速度，t为当前的迭代次数，iterjnax为最大迭代次数，δ为随机方向向 量；
[0038] d)对磷虾的位置进行更新，更新公式为：
式中=X1U+Δ t)为磷虾个体i在t+Δ t时刻的位置，Xi(t)为磷虾个体i在t时刻的位置，△ t为时间间隔；
[0039] (6)计算磷虾个体的适应度函数值，找出当前适应度函数值最优的磷虾个体；
[0040] (7)循环（4)、（5)、（6)三个步骤30次后，搜索到适应度函数值的最大值 12.5319036419，该磷虾所处的位置[107，161 ]为摇床矿带分带图像的最佳阈值，根据求出 的最佳阈值将摇床矿带分割为精矿、中矿和尾矿。
[0041 ]实施例2:参见图1，利用VC++软件对摇床床面上形成的锡矿的精矿、次精矿、中矿 和尾矿矿带进行分割，采用的方法和步骤与实施例1相同，其中最大迭代次数为40:
[0042] (1)对彩色的摇床矿带分带图像进行灰度化处理，将彩色的摇床矿带分带图像转 换为灰度图像；
[0043] (2)磷虾个体的初始化:对参数进行设置，觅食速度Vf = 0.02,最大扩散速度Dmax = 0.005，最大移动速度Nmax = 0.01，最大迭代次数iterjnax = 30和磷虾的个数M= 25;最大值 惯性权重ω max = 〇 . 9，惯性权重的最小值ω min = 〇. 1;产生M个磷虾个体均匀地分布在摇床矿 带分带图像[0,255]的灰度直方图空间中；
[0044] (3)利用Kapur熵的公式H(ti，t2, t3) =Ηι+Η2+Η3+Η4计算磷奸的适应度函数值，并根 据适应度函数值，对磷虾个体进行最好到最差的排列，25个磷虾个体的适应度函数值为：
[0045] H= [14.4145962493,14.8931997116,13.5053071374,13.7716234406， 11.0845701850,12.5734712465,15.1428814259
0];
[0046] (4)利用公式 （寸惯性权重~进行更新，式中：ωη为惯性 权重，t为迭代次数，iter_max为最大迭代次数，ω max为惯性权重的最大值，ω min为惯性权重 的最小值；
[0047] (5)对每个磷虾个体i，执行下列步骤：
[0048] a)计算磷虾个体i的运动向量焉=#1^(<；^+?严，+ ?"7/；'式中41为运动向 量，Nmax为最大移动速度，为磷虾个体i的受周围邻居吸引产生的运动向量，"为当 前最优磷虾个体吸引磷虾个体i产生的运动向量，ω n为惯性权重，祝f为磷虾个体i的上一 次运动向量；
[0049] b)计算觅食运动向量巧:=匕扃+巧，其中：爲=汊广'爲为食物的吸 引力，爲^为到目前为止适应度函数值最大的磷虾个体i的影响力，F1为磷虾个体i的觅食 运动向量;V f为觅食速度，觅食运动向量的惯性权重，为磷虾个体i的上一次觅食 运动向量；
[0050] c)计算物理扩散运动向量
1式中：Di为物理扩散运动向 量，Dmax为最大扩散速度，t为当前的迭代次数，iterjnax为最大迭代次数，δ为随机方向向 量；
[0051] d)对磷虾的位置进行更新，更新公式为：
式中=X1U+Δ t)为磷虾个体i在t+Δ t时刻的位置，Xi(t)为磷虾个体i在t时刻的位置，△ t为时间间隔；
[0052] (6)计算磷虾个体的适应度函数值，找出当前适应度函数值最优的磷虾个体；
[0053] (7)循环(4)、（5)、（6)三个步骤40次后，通过适应度函数的最大化，搜索到适应度 函数值的最大值15.2997820216，该磷虾所处的位置[117，150，182]为摇床矿带分带图像的 最佳阈值，根据求出的最佳阈值将摇床矿带分割为精矿、次精矿、中矿和尾矿。
[0054]实施例3:参见图1，对锡矿的精矿、次精矿、中矿和尾矿矿带进行分割，采用的方法 和步骤与实施例1相同，最大迭代次数为40，其中利用Kapur熵H(ti，t2，t3) =?+?+?+?,来 计算磷虾个体的适应度函数值，通过迭代步骤(4)、（5)、（6)40次后，通过使适应度函数值最 大化，搜索到适应度函数的最大值15.3772825987,搜索出矿带图像的最佳阈值[118,150， 181]，根据搜索到的最佳阈值对摇床矿带分带进行分割，从而把锡矿矿带分带分割为精矿、 次精矿、中矿和尾矿四个矿带。
[0055]利用Kapur熵计算磷虾个体的适应度函数值，迭代次数为40次时，25个磷虾的适应 度函数值分别为：
[0056] H= [15.3772825987,15. 1725749681,14.8666985095,13.8595522641 , 13.6070001036,12.7124563267,12.3188929863,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0]〇
[0057] 实施例4:为了验证改进的磷虾优化算法的性能，本实施例中选取了四幅摇床矿带 分带图像，利用基本的磷虾优化算法和本发明中改进的磷虾优化算法对四幅图像进行分 害J，两种算法采用的参数与实施例1中的相同，本发明中提出的磷虾优化算法随着迭代次数 的变化对惯性权重ω "进行更新，提高了算法的搜索速度和搜索精度，能够快速和准确地搜 索到摇床矿带分带图像的最佳适应度函数值和最佳阈值。
[0058]表1:实验比较结果
[0060]上面结合附图对本发明的【具体实施方式】作了详细说明，但是本发明并不限于上述 实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前 提下作出各种变化。
【主权项】
1. 一种基于憐邮优化算法的摇床矿带分带图像分割法，其特征在于:包括如下步骤： (1)在摇床选矿的过程中，实时拍摄摇床床面上形成的矿带分带图像，将彩色的摇床矿 带分带图像转换为灰度图像； (2 )憐邮个体的初始化:设置憐邮算法的参数，觅食速度Vf，最大扩散速度，最大移动 速度Ν"3χ，最大迭代次数i ter_max和憐邮的个数Μ;产生[0，255 ]之间均匀分布的憐邮个体， 让Μ个憐邮个体均匀地分布在摇床矿带分带图像的灰度直方图空间中； (3) 计算憐邮的适应度函数值，并根据适应度函数值，对憐邮个体进行从最好到最差的 排列； (4) 利用公式巧惯性权重ωη进行更新，式中：ωη为惯性权重， t为迭代次数，iter_max为最大迭代次数，ω max为惯性权重的最大值，ω min为惯性权重的最 小值； (5) 对每个憐邮个体i，执行下列步骤： a)计算憐邮个体i的运动向量式中:Ni为运动向量，Nmax 为最大移动速度，为憐邮个体i的受周围邻居吸引产生的运动向量，苗为当前最优 憐邮个体吸引憐邮个体i产生的运动向量，ωη为惯性权重，为憐邮个体i的上一次运动 向量； b )计算觅食运动向量Fi = Vf扣+ ω ，其中：wd+pibest，心DDd为食物的吸引力， 0产％到目前为止适应度函数值最大的憐邮个体i的影响力，F功憐邮个体i的觅食运动向 量;Vf为觅食速度，Wf为觅食运动向量的惯性权重，Fi°id为憐邮个体i的上一次觅食运动向 量； C)计算物理扩散运动向量，式中：Di为物理扩散运动向量， 为最大扩散速度，t为当前的迭代次数，iter_max为最大迭代次数，δ为随机方向向量； d)对憐邮的位置进行更新，更新公式为式中：Xi(t+At)为 憐邮个体i在t+At时刻的位置，Xi(t)为憐邮个体i在t时刻的位置，At为时间间隔； (6) 计算憐邮个体的适应度函数值，找出当前适应度函数值最优的憐邮个体； (7) 不断循环(4)、（5)、（6)Ξ个步骤，直到满足最大迭代次数titer_"ax，通过适应度函数 的最大化，找出经过多次迭代后适应度函数值最优的憐邮个体，该憐邮所处的位置即为摇 床矿带分带图像的最佳阔值，最后根据求出的最佳阔值对摇床矿带分带图像进行分割。2. 根据权利要求1所述的基于憐邮优化算法的摇床矿带分带图像分割方法，其特征在 于:最大迭代次数iterjnax范围为[30,50 ]，憐邮的个数Μ的范围为[25,50 ]。3. 根据权利要求1所述的基于憐邮优化算法的摇床矿带分带图像分割法，其特征在于： 步骤(4)中惯性权重的最大值ω max为0.9，惯性权重的最小值ω min为0.1。
【文档编号】G06Q50/02GK105844628SQ201610159915
【公开日】2016年8月10日
【申请日】2016年3月21日
【发明人】和丽芳, 郭思哲, 黄宋魏, 童雄, 郝鹏宇, 司绪, 张元元
【申请人】昆明理工大学