Matlab双峰法求阈值的代码

marazone

从其他论坛上找到的跟大家分享一下,记住要把图像文件放在m文件那个目录，smart.jpg是我电脑里面的图片，你们需要自己修改，开头和末尾的tic toc是求运行时间的

1. tic
   
2. 
   
3. clear;clc;clear all
   
4. p=imread('smart.jpg');
   
5. if ndims(p) == 3
   
6.     p = rgb2gray(p);  %转化为灰度图像
   
7. end
   
8. subplot(1,3,1)
   
9. imshow(p)                                               
   
10. title('原来的灰度图像')
    
11. [count,x]=imhist(p);
    
12. subplot(1,3,2);
    
13. plot(x,count)
    
14. title('分布直方图')
    
15. 
    
16. zl=min(min(p));              %图像的最小灰度值
    
17. zk=max(max(p));              %最大图像的灰度值
    
18. zl=double(zl);
    
19. zk=double(zk);
    
20. T(1)=(zl+zk)/2;              
    
21. 
    
22. T(1)=round(T(1));            %设定初值为最值的平均值
    
23. 
    
24. for k=1:6
    
25.      n1=sum(count(1:T(k)));
    
26.      n2=sum(count(T(k)+1:256));
    
27.      sum1=0;sum2=0;
    
28.      for i=1:T(k)
    
29.          sum1=sum1+count(i)*i;
    
30.      end
    
31.      for i=T(k)+1:256
    
32.          sum2=sum2+count(i)*i;
    
33.      end
    
34.      ua1=sum1/n1;   
    
35.      ua2=sum2/n2;  
    
36.     T(k+1)=(ua1+ua2)/2;
    
37.     T(k+1)=round(T(k+1));     
    
38.     if T(k+1)==T(k)
    
39.         break;
    
40.     end
    
41. end
    
42. if T(k+1)==T(k);
    
43.     TT=T(k+1);
    
44.     disp('The threshold is')
    
45.     TT
    
46. else
    
47.     TT=0;
    
48.     disp('error');   
    
49. end
    
50. %得到了阈值TT
    
51. [m,n]=size(p);
    
52. R=zeros(m,n);
    
53. for i=1:m
    
54.     for j=1:n
    
55.         if p(i,j)<TT
    
56.            R(i,j)=0;
    
57.         else
    
58.            R(i,j)=255;
    
59.         end
    
60.     end
    
61. end
    
62. R=uint8(R);
    
63. subplot(1,3,3);
    
64. imshow(R);
    
65. title('转化后的图像');
    
66. 
    
67. toc

复制代码

sleep_ball

为什么要用matlab求阈值

marazone

sleep_ball 发表于 2012-3-2 20:56
为什么要用matlab求阈值

matlab方便观看效果，如果用C语言处理的话，看到二值化前后的图片没有那么方便，matlab只是工具，最后还是得用C语言实现

abc李瑞华

王谦623

佩服！