【创龙TMS320C665X开发板试用】二值化

`在刚接触DSP图像算法时毫无头绪不知从哪里开始搞起，幸好看到jinyi7016的帖子，使我对这部分有了初步认识，在这里再次谢过。在本程序中参考了一部分jinyi7016的程序。由于例程中阈值的选取是靠人为设定一个值，因此我在阈值这部分做了一些改进。加入了自适应阈值选取（大津法取阈值）的算法。

int i,j,k,N;

double foregraysum ,backgraysum; //前景与背景

double n0=0,n1=0,u0=0,u1=0,w0=0,w1=0,max=-1,g=-1;

double histogram[256]={0};

unsigned char *thre***uf;

thre***uf=m->imgBuf;

N=m->info.biHeight*m->info.biWidth;

/***************************************灰度直方图*********************************************/

for(i=0;i
{

double temp=thre***uf[i];

histogram[(int)temp]++;

}

/****************************************计算阈值********************************************/

for(i=0;i<256;i++)

  {

  foregraysum=0;    //前景灰度总值

backgraysum=0;    //背景灰度总值

n0+=histogram[i]; //前景的像素数

n1=N-n0;    //背景的像素数

if( n1==0 ) {break;} //没有出现特征物体，跳出循环

w0=n0/N;    //前景像素点数占整幅图像的比例为w0

  w1=1-w0;    //背景像素点数占整幅图像的比例为w1

  for(j=0;j<=i;j++)

  {

foregraysum+=j*histogram[j]; //前景总灰度

}

u0=foregraysum/n0;    //前景平均灰度

for(k=i+1;k<256;k++)

  {

backgraysum+=k*histogram[k]; //背景总灰度

}

u1=backgraysum/n1;    //背景平均灰度

g=w0*w1*(u0-u1)*(u0-u1); //类间方差

if(max
  {

max=g;

T=i;

  }

  }

复制代码