【图像分割】粒子群算法OTSU图像分割【含Matlab源码 3688期】.zip

function [G_best1,G_best2,D]=PSO(Pso_no,Max_iteration,lb,ub,I,image_c) L_a=max(max(image_c)); L_b=min(min(image_c)); V_max=4; Leader_score1=-inf; %change this to -inf for maximization problems Leader_score2=-inf; %统计图像整体灰度平均值u，背景灰度平均值u11，背景概率w11 [ua11,w11,pp1,nd1,u1]=Otsu_q(I); %统计图像整体灰度平均值u，背景灰度平均值u11，背景概率w11 [ua22,w22,pp2,nd2,u2]=Otsu_q(image_c); %%%%粒子群位置及速度的初始化 Positions1=intialize_position_rand(Pso_no,ub,lb);%粒子位置初始化 Positions2=intialize_position_rand(Pso_no,L_a,L_b);%粒子位置初始化 V_pso2=intialize_position_rand(Pso_no,2-V_max,V_max-2);%粒子位置初始化 V_pso1=intialize_position_rand(Pso_no,-V_max,V_max);%粒子位置初始化 D=[];t=0; P_best1=zeros(1,Pso_no); P_best2=zeros(1,Pso_no); d111=zeros(1,Pso_no); d222=zeros(1,Pso_no); while t<Max_iteration for i=1:Pso_no %粒子位置重置 if Positions1(i)>ub %粒子位置超出上界设置为255 Positions1(i)=255; end if Positions1(i)<lb %粒子位置超出下界设置为0 Positions1(i)=0; end % Return back the search agents that go beyond the boundaries of the search space if Positions2(i)>ub %粒子位置超出上界设置为255 Positions2(i)=255; end if Positions2(i)<lb %粒子位置超出下界设置为0 Positions2(i)=0; end %粒子速度重置 if V_pso1(i)>V_max %粒子速度超出上界设置为V_max V_pso1(i)=V_max; end if V_pso1(i)<-V_max %粒子速度超出下界设置为-V_max V_pso1(i)=-V_max; end if V_pso2(i)>V_max-2 %粒子速度超出上界设置为V_max-2 V_pso2(i)=V_max-2; end if V_pso2(i)<2-V_max %粒子速度超出下界设置为2-V_max V_pso2(i)=2-V_max; end end %由最大类间方差公式计算每个位置（灰度值）的类间方差 [Max_jubu1,d1,d11]=Otsu_h(ua11,w11,Positions1,pp1,nd1,u1);%P_best1是得带过程中每个粒子的最优位置 [Max_jubu2,d2,d22]=Otsu_h(ua22,w22,Positions2,pp2,nd2,u2);%d22是类间方差组 %全体最优位置 if d1>Leader_score1 Leader_score1=d1; G_best1=Max_jubu1;%G_best1是迭代中全局最优位置 end if d2>Leader_score2 Leader_score2=d2; G_best2=Max_jubu2;%G_best2是迭代中全局最优位置 end %P_best1是得带过程中每个粒子的最优位置 for i=1:Pso_no if d11(i)>d111(i) d111(i)=d11(i); P_best1(i) =Positions1(i); end if d22(i)>d222(i) d222(i)=d22(i); P_best2(i) =Positions2(i); end end %给各个系数赋值 w=0.5; r1=rand; r2=rand; c1=2; c2=2; for i=1:Pso_no V_pso1(i)=w*V_pso1(i)+c1*r1*(P_best1(i)-Positions1(i))+c2*r2*(G_best1-Positions1(i));%改变粒子的速度 V_pso2(i)=w*V_pso2(i)+c1*r1*(P_best2(i)-Positions2(i))+c2*r2*(G_best2-Positions2(i)); Positions1(i)=Positions1(i)+V_pso1(i);%改变粒子的位置 Positions2(i)=Positions2(i)+V_pso2(i); end Positions1=double(int32(Positions1)); Positions2=double(int32(Positions2)); D=[D,G_best1]; t=t+1; end end