一维最大熵法进行2分割和3分割（MATLAB）

Halcom

一维最大熵法进行2分割和3分割（MATLAB）
链接：https://pan.baidu.com/s/1M15vV4l8R96hahprmhXDEA&shfl=sharepset   提取码：hobi
具体链接在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. % 2个阈值，实现3分类分割
   
4. I = imread('car.bmp');
   
5. % I = imread('screen.jpg');
   
6. figure('color',[1,1,1]), imshow(I,[])
   
7. t1 = 128;  % 初始化最佳阈值
   
8. t2 = 180;  % 初始化最佳阈值
   
9. Lmax = 256;  % 最大灰阶, uin8
   
10. k=1;
    
11. for t1=0:Lmax-1
    
12.     for t2 = t1+1:Lmax
    
13.         pop(k,:) = [t1,t2];
    
14.         H1(k) = H1_segFun2(I, pop(k,:) );
    
15.         k=k+1;
    
16.     end
    
17. end
    
18. [Hmax, bestT]=max(H1);
    
19. 
    
20. % 显示
    
21. Iout = imageGray(I, pop(bestT,:) );
    
22. figure('color',[1,1,1])
    
23. imshow(Iout,[])
    
24. colormap jet;
    
25. axis equal

复制代码

适应度函数H1_segFun2如下：

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

显示函数imageGray如下：

1. function imgOut = imageGray(img, thre)
   
2. % 图像按照不同的灰阶值显示
   
3. limits = [0, thre, 255];
   
4. imageSize = size(img);
   
5. imgOut(:,:,:)=img*0;
   
6. 
   
7. k=1;
   
8. for i=1:imageSize(1,1)
   
9.     for j=1:imageSize(1,2)
   
10.         while(k<size(limits, 2))
    
11.             if(img(i,j)>=limits(1,k) && img(i,j)<=limits(1,k+1))
    
12.                 imgOut(i,j,1) = limits(1, k);  
    
13.             end
    
14.             k=k+1;
    
15.         end
    
16.         k=1;
    
17.     end
    
18. end

复制代码

我成功失败了

感谢管理员分享的资源！

yuyibintony

好东西，谢谢分享！

842013760

非常好得代码