遗传算法的图象分割matlab程序

frogfish

用matlab写的遗传算法分割程序，用的是Sheffield大学gatbs工具箱
能够运行得到真确结果！

a10.m

1. O=imread('Lenna.bmp');
   
2. figure(1);imshow(O);
   
3. 
   
4. NIND = 40; % Number of individuals
   
5. MAXGEN =50; % Maximum no. of generations
   
6. PRECI = 8; % Precision of variables
   
7. GGAP = 0.9; % Generation gap
   
8. % Build field descriptor
   
9. FieldD = [8;1;256;1;0;1;1];
   
10. % Initialise population
    
11. Chrom = crtbp(NIND, PRECI);
    
12. gen = 0;% Counter
    
13. phen=bs2rv(Chrom,FieldD);% Evaluate initial population
    
14. ObjV = target(O,phen);
    
15. while gen < MAXGEN,% Generational loop
    
16. % Assign fitness values to entire population
    
17. FitnV = ranking(-ObjV);
    
18. SelCh = select('sus', Chrom, FitnV, GGAP);% Select individuals for breeding
    
19. SelCh = recombin('xovsp',SelCh,0.7);% Recombine individuals (crossover)
    
20. SelCh = mut(SelCh);% Apply mutation
    
21. phenSel = bs2rv(SelCh,FieldD);% Evaluate offspring, call objective function
    
22. ObjVSel=target(O,phenSel);
    
23. [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);% Reinsert offspring into population
    
24. gen = gen+1;% Increment counter
    
25. end
    
26. [Y,I]=max(ObjV);
    
27. M=bs2rv(Chrom(I,:),FieldD)
    
28. [m,n]=size(O);
    
29. for i=1:m
    
30.     for j=1:n
    
31.         if O(i,j)>M
    
32.             O(i,j)=256;
    
33.         end
    
34.     end
    
35. end
    
36. figure(2)
    
37. imshow(O);

复制代码

target.m

1. function f=target(T,M)
   
2. dT=double(T);
   
3. [U,V]=size(T); W=length(M);f=zeros(W,1);
   
4. for k=1:W
   
5. s1=0;I=0;s2=0;J=0;
   
6. for i=1:U
   
7.     for j=1:V
   
8.         if T(i,j)<=M(k)
   
9.             s1=s1+dT(i,j);I=I+1;
   
10.         end
    
11.         if T(i,j)>M(k)
    
12.             s2=s2+dT(i,j);J=J+1;
    
13.         end
    
14.     end
    
15. end
    
16. p1=s1/I;p2=s2/J;
    
17. f(k)=I*J*(p1-p2)*(p1-p2)/(256*256*256*256);
    
18. end

复制代码

本程序由研学hwh4103

amanda-mat

能不能说一下具体要怎么运行？还需要什么初始文件？图像分割程序，那是不是要图像的原始文件啊？谢谢！

hjd2003

谢谢楼主分享！学习学习！！！！！