大津法问题

wx_o19nT8nV

不知道怎么自己原理理解错了还是怎么，大津法返回来的阈值是一个定值，程序如下

1. void get_theos()
   
2. {
   
3.   const uint16 hui_du = 256;                 //灰度级数
   
4.   const unsigned long long pixelsum = ROW * COL;        //像素总数
   
5.   double zhifangtu[256]={0};                    //直方图数组
   
6.   double pixelcount[256]={0};                    //各个灰度级数的概率
   
7.   uint16 i=0,j=0,x=0,y=0,t=0;
   
8.   double p0=0,p1=0,u1=0,u0=0,u00=0,u11=0;
   
9.   double threshold=0;
   
10.   double fangcha=0.0,fangchamax=0.0;
    
11. 
    
12.   
    
13.   for(i=0;i<ROW;i++)
    
14.     for(j=0;j<COL;j++)
    
15.       zhifangtu[image[i][j]]++;            //算出直方图
    
16.   /*********测试直方图************/
    
17. // for(t=0;t<hui_du;t++)
    
18.   //SCI_Up_acc_gyr((short) zhifangtu[t],(short) 0);
    
19. ///uart_putchar(UART2,zhifangtu[t]);
    
20.    
    
21.   for(x=0;x<hui_du;x++)
    
22.   {
    
23.      pixelcount[x]=zhifangtu[x]/pixelsum;
    
24. //    if((0<pixelcount[x])&&(pixelcount[x]<1))
    
25.    //   y++;
    
26.   }
    
27.   
    
28.   for(x=0;x<hui_du;x++)
    
29.   {
    
30.     p0=p1=u0=u1=fangcha=u00=u11=0.0;
    
31.     for(y=0;y<hui_du;y++)
    
32.     {
    
33.       if(y<=x)
    
34.       {
    
35.         u00+=y*pixelcount[y];
    
36.         p0+=pixelcount[y];
    
37.       }
    
38.       else
    
39.       {
    
40.         u11+=y*pixelcount[y];
    
41.       }
    
42.     }
    
43.     p1=1-p0;
    
44.     u0=u00/p0;
    
45.     u1=u11/p1;
    
46.     fangcha=p0*p1*(u0-u1)*(u0-u1);
    
47.     if(fangcha>fangchamax)
    
48.     {
    
49.       fangchamax=fangcha;
    
50.       threshold=i;
    
51.     }
    
52.   }
    
53. // SCI_Up_acc_gyr((short) threshold,(short) y);
    
54. }
    

复制代码

哞哞哞哒

目前我也不太理解大津法

哞哞哞哒

http://www.doc88.com/p-666157374161.html
我照着敲了一遍，还没用

wx_o19nT8nV

哞哞哞哒 发表于 2018-1-28 20:24
http://www.doc88.com/p-666157374161.html
我照着敲了一遍，还没用

看这个程序大致框架和思路是没有问题的，你可以试试。好了可以交流下

哞哞哞哒

返回来的值为0，不知道错哪了，我换了双峰法试了一下，值是会变了，但效果不怎么好。晚上再调试看看错哪了