ostu法图像分割

Halcom

主函数：

1. global ImageGray
   
2. hist = imhist(ImageGray)./(size(ImageGray,1)*size(ImageGray,2));  % 直方图
   
3. for i=1:256
   
4.     fit(i) = image_segmentation_ostu(hist,i);  % Ostu
   
5. end
   
6. [maxfit,thre] = max(fit);
   
7. bw = ImageGray > thre;   % 二值化操作

复制代码

子函数：

1. function y = image_segmentation_ostu( hist, thre )
   
2. % 最大类间方差法
   
3. Lmax = 256;
   
4. 
   
5. y = 0;
   
6. p1 = 0; p2=0; p3 =0; p4 =0; p5=0;
   
7. 
   
8. for i=1:Lmax
   
9.     p1 = p1+i*hist(i);
   
10.     if i<thre
    
11.         p2 = p2 + hist(i);
    
12.     else
    
13.         p3 = p3 + hist(i);
    
14.     end
    
15. end
    
16. 
    
17. for i=1:Lmax
    
18.     if i<thre
    
19.         p4 = p4 + i*hist(i)/p2;
    
20.     else
    
21.         p5 = p5 + i*hist(i)/p3;
    
22.     end
    
23. end
    
24. 
    
25. y = p2*(p4-p1)*(p4-p1);
    
26. y = y+p3*(p5-p1)*(p5-p1);

复制代码