PSO的二维最大熵法阈值分割

Halcom

PSO的二维最大熵法阈值分割
链接：https://pan.baidu.com/s/1F3DToLeT1JuAGg7a2wgkhQ&shfl=sharepset 提取码：9p1k
具体链接在halcom.cn论坛，联系人QQ：3283892722
该论坛是一个学习交流平台，我会逐一的和大家分享学习。
欢迎大家录制视频，你可在论坛进行打赏分享。
视频专用播放器：http://halcom.cn/forum.php?mod=viewthread&tid=258&extra=page%3D1\

1. clc,clear,close all
   
2. warning off
   
3. format longG
   
4. % 问题背景
   
5. I =imread('car.bmp');
   
6. I = imresize( I, [256,256] );
   
7. if(size(I, 3)>1)
   
8.     im = rgb2gray(I);
   
9. else
   
10.     im = I;
    
11. end
    
12. h = fspecial('average', [2,2]);  % 2x2均值滤波
    
13. meanI = imfilter(im, h);         % 均值图像
    
14. figure(1), imshow(im, []);
    
15. %% PSO 参数
    
16. c1 = 1.4995;  
    
17. c2 = 1.4995;
    
18. Vmin = -5;
    
19. Vmax = 5;
    
20. maxiter = 10;  % 迭代次数
    
21. sizepop = 10;  % 种群数量
    
22. Lmax = 256;    % 最大灰阶
    
23. popmin = 1;  popmax = Lmax-1;   % 未知量个数，决策变量个数，阈值个数 x  取值范围
    
24. % 2个阈值，实现3分割
    
25. level = 3;
    
26. nvar = level-1;   % 未知量个数，决策变量个数，阈值个数
    
27. 
    
28. % 初始化种群
    
29. for i=1:sizepop
    
30.     pop(i,1) = popmin + fix( (popmax-popmin)*rand );
    
31.     pop(i,2) = popmin + fix( (popmax-popmin)*rand );
    
32.     pop(i, :) = sort(pop(i, :));
    
33.     fitness(i) = fun2(  im, meanI, pop(i,:) );
    
34.     V(i,:) = zeros(1, nvar);
    
35. end
    
36. % 记录一组最优值
    
37. [bestfitness,bestindex]=max(fitness);
    
38. zbest=pop(bestindex,:);   %全局最佳
    
39. gbest=pop;                %个体最佳
    
40. fitnessgbest=fitness;     %个体最佳适应度值
    
41. fitnesszbest=bestfitness; %全局最佳适应度值
    
42. % 迭代寻优
    
43. for i=1:maxiter
    
44.    disp(['当前迭代次数: ', num2str(i)])
    
45.     for j=1:sizepop
    
46.         % 速度更新
    
47.         V(j,:) = V(j,:) + c1*rand*(gbest(j,:) - pop(j,:)) + c2*rand*(zbest - pop(j,:));
    
48.         V(j, find(V(j,:)>Vmax)) = Vmax;
    
49.         V(j, find(V(j,:)<Vmin)) = Vmin;
    
50.         
    
51.         % 个体更新
    
52.         pop(j,:) = fix( pop(j,:) + 0.5 * V(j,:) );
    
53.         pop(j, find(pop(j,:)>popmax)) = popmax;
    
54.         pop(j, find(pop(j,:)<popmin)) = popmin;
    
55.         
    
56.         % 适应度更新
    
57.         pop(j, :) = sort(pop(j, :));
    
58.         fitness(j) = fun2(im, meanI, pop(j,:));
    
59.         
    
60.         % 比较  个体间比较
    
61.         if fitness(j)>fitnessgbest(j)
    
62.             fitnessgbest(j) = fitness(j);
    
63.             gbest(j,:) = pop(j,:);
    
64.         end
    
65.         if fitness(j)>bestfitness
    
66.             bestfitness = fitness(j);
    
67.             zbest =  pop(j,:);
    
68.         end
    
69.         
    
70.     end
    
71.     fitness_iter(i) = bestfitness;
    
72.    
    
73. end
    
74. %% 结果显示
    
75. disp('最优解')
    
76. zbest = sort(zbest);
    
77. disp(zbest)
    
78. fprintf('\n')
    
79. 
    
80. figure('color',[1,1,1])
    
81. plot(fitness_iter,'ro-','linewidth',2)
    
82. 
    
83. % 显示
    
84. Iout = imageGray(im, zbest);
    
85. figure('color',[1,1,1])
    
86. imshow(Iout,[])
    
87. colormap jet;
    
88. axis equal

复制代码

适应度函数如下：

游客，如果您要查看本帖隐藏内容请回复

参考：【32】一维最大熵法进行2分割和3分割（MATLAB）

wal1985

非常好的代码

余坤键

我来看看大佬的代码

yuyibintony

好资源，谢谢分享！

wwmwwmwwm

111111111111111111111