直方图均衡化—图像增强

void HistogramEqualization(unsigned char *pGrey,int nWidth, int nHeight)
{
	unsigned char* pData,pInResult[256];
	int			   pHist[256],i,j,nMulti;
	float		   pIntensity[256];
	pData = (unsigned char*)malloc(nWidth*nHeight*sizeof(unsigned char));
	memcpy(pData, pGrey, nWidth*nHeight*sizeof(unsigned char));
	memset(pHist, 0, sizeof(int)*256);
	nMulti = nWidth * nHeight;

	//统计原始输入图像各灰度级的像素数目
	for (i=0; i<nMulti; i++)
	{
		pHist[pData[i]]++;
	}
	//计算原始图像直方图，即各灰度级的概率密度,并计算累积分布函数
	memset(pIntensity, 0, sizeof(float)*256);
	pIntensity[0] = (float)pHist[0] / (float)nMulti;
	for (i=1; i<256; i++)
	{
		pIntensity[i] = (float)pHist[i] / (float)nMulti;
		pIntensity[i] += pIntensity[i-1];
	}
	//计算最后的输出灰度级

	for (i=0; i<256; i++)
	{
		pInResult[i] = (int)(255*pIntensity[i] + 0.5);
	}

	//修改原图像的灰度级，获得输出图像，其直方图为近似均匀分布
	for(i=0; i<nMulti; i++)
	{
		pGrey[i] = pInResult[pData[i]];
	}

	free(pData);

}