数字图像的熵(C)

熵指的是体系的混乱的程度，对焦良好的图像的熵大于没有清晰对焦的图像，因此可以用熵作为一种对焦评价标准。熵越大，图像越清晰。

double LOG(double a)  
{  
   int N = 15;
   int k,nk;  
   double x,xx,y;  
   x = (a-1)/(a+1);  
   xx = x*x;  
   nk = 2*N+1;  
   y = 1.0/nk;  
   for(k=N;k>0;k--)  
   {  
   	nk = nk - 2;  
   	y = 1.0/nk+xx*y;  
   }  
   return 2.0*x*y;   
}  
double entropy_cal(U8* src,int size)
{
	int i,j,img_p;
	unsigned char img[size*size],img_gray[256];
	double H,Ph[256];
	H=0;
	img_p=0;
	for(i=0;i<256;i++)
	{
		img_gray[i]=0;
		Ph[i]=0;
	}
	for(j=pAppData->nHeight/2-size/2;j< pAppData->nHeight/2+size/2;j++)
	   for(i=pAppData->nWidth/2-size/2;i< pAppData->nWidth/2+size/2;i++)
	   {
		img[img_p++]=src[j*pAppData->nWidth*2+2*i];
	   }
	for(i=0;i<size*size;i++)
	{
		img_gray[img[i]]++;
	}
	for(i=0;i<256;i++)
	{
		if(img_gray[i])
		{
			Ph[i]=(double)img_gray[i]/(size*size);
			H=H-Ph[i]*LOG(Ph[i])/LOG(2);
		}
	}
	return H;
}