粒子群图像分割

pso 图像分割本科毕业论文里的东西,无论是彩色还是灰度都可以

clc;clear all;close all;
t1=clock;
I=imread（‘a.jpg‘）;
%J=rgb2gray（I）;%如果是彩色图像要加此句，但去掉下面一句
 J=I;
[ab]=size（J）;
figure（1）
imshow（J）
[px]=imhist（J256）;
figure（2）
plot（p）;
title（‘灰度直方图‘）
xlabel（‘灰度值‘）
ylabel（‘像素点个数‘）
L=x‘;
LP=p‘/（a*b）;

n=256
c1=2;
c2=2;
wmax=0.9;
wmin=0.4;
G=100;
M=10;
X=min（L）+fix（（max（L）-min（L））*rand（1M））;
V=min（L）+（max（L）-min（L））*rand（1M）;


m=0;
for i=1:1:n
    m=m+L（i）*LP（i）;
end
pbest=zeros（M2）;
gbest1=0;
gbest2=0;
GG=0;
t2=clock;
for k=1:1:G
    w（k）=wmax-（wmax-wmin）*k/G;
    for i=1:1:M
        t=length（find（X（i）>=L））;
        r=0;
        s=0;
        for j=1:1:t
            r=r+LP（j）;
            s=s+L（j）*LP（j）;
        end
        W0（i）=r;
        W1（i）=1-r;
        U0（i）=s/r;
        U1（i）=（m-s）/（1-r）;
    end
    for i0=1:1:M
        BB（i0）=W0（i0）*W1（i0）*（（U1（i0）-U0（i0））^2）;
    end
  
   for i=1:1:M
       if pbest（i2）           pbest（i2）=BB（i）;            
           pbest（i1）=X（i）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      21034  2012-11-30 10:51  图像分割1\a.jpg

     文件       4170  2012-05-18 11:14  图像分割1\pic1.fig

     文件      65155  2012-05-18 11:14  图像分割1\pic1.jpg

     文件       4284  2012-05-18 11:20  图像分割1\pic2.fig

     文件      74224  2012-05-18 11:20  图像分割1\pic2.jpg

     文件       4497  2012-05-18 11:30  图像分割1\pic3.fig

     文件      62023  2012-05-18 11:30  图像分割1\pic3.jpg

     文件       9096  2012-05-18 10:47  图像分割1\picture1.jpg

     文件     114592  2012-05-18 10:49  图像分割1\picture2.jpg

     文件     165880  2012-05-18 10:59  图像分割1\picture3.jpg

     文件       1593  2012-11-30 14:49  图像分割1\pso.m

     目录          0  2012-12-10 21:15  图像分割1

----------- ---------  ---------- -----  ----

               526548                    12