30-基于中值的Ostu图像分割方法

Halcom

30-基于中值的Ostu图像分割方法
百度网盘链接：https://pan.baidu.com/s/1Ciw1JIZC_joe-AAbzb2JnA
具体链接在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. im = imread('lena.tif');
   
4. figure(1), imshow(im)
   
5. 
   
6. [count, grayvalue] = imhist(im);
   
7. L = 256;
   
8. if(size(im,3)>1)
   
9.     im = rgb2gray(im);
   
10. end
    
11. [m,n] = size(im);
    
12. f = count./(m*n);
    
13. for i=1:L
    
14.     if(f(i)~=0)
    
15.         st = i;
    
16.         break;
    
17.     end
    
18. end
    
19. for i=L:-1:1
    
20.     if(f(i)~=0)
    
21.         ed = i;
    
22.         break;
    
23.     end
    
24. end
    
25. j=1;
    
26. for i=1:ed-st
    
27.     w0(j) = sum( f(1:i) );
    
28.     j=j+1;
    
29. end
    
30. w0 = w0+0.001;
    
31. w1 = 1-w0;
    
32. m0 = median(1:ed-st);
    
33. m1 = median(ed-st+1:L);
    
34. 
    
35. mad1=0;
    
36. for i=1:ed-st
    
37.     mad1 = mad1+f(i)./w0 * abs(i-m0);
    
38. end
    
39. mad2=0;
    
40. for i=ed-st+1:L
    
41.     mad2 = mad2+f(i)./w1 * abs(i-m1);
    
42. end
    
43. 
    
44. k=w0.*mad1 + w1.*mad2;
    
45. [value, index] = min(k);
    
46. thr = st + index;
    
47. 
    
48. bw = im>thr;
    
49. figure(2), imshow(bw)
    

复制代码