直方图均衡化

直方图均衡化的步骤如下

参考 http://wenku.baidu.com/link?url=yHCFrzOZTiqWhYTa0TH-8kpxyQry7D2dZQSombvAKrOrIP8LcL_VXCEKb3Pi9EaQdgkYdIPpPsTgvj5ZzQEcGWxm2LKvZg9cKmvHF8gSCNy

以上两张PPT讲的非常清楚了。

下面贴出代码

#include <stdio.h>
#include <tchar.h>
#include "highgui.h"
#include "cv.h"
using namespace cv;

void DrawHistImg(int hist[])
{
	printf("Hist image :\n");
	for (int i = 0;i<255;i++)
	{
		printf("pixel %3d:",i);
		//太长显示难看，因此除以50
		for(int j = 0;j<hist[i]/50;j++)
		{
			printf("-");
		}
		printf("\n");
	}
}

int _tmain(int argc, _TCHAR* argv[])
{
	//读图
	Mat img = imread("E:/【照片图片】/image/007.jpg",0);
	imshow("Source Image",img);
	//waitKey(0);

	int hist[255];
	memset(hist,0,sizeof(int)*255);
	//统计直方图
	long sumPixel = 0;
	for (int i = 0;i<img.rows ;i++)
		for (int j = 0;j<img.cols ;j++)
	{
		int tmp = img.at<uchar>(i,j);
		sumPixel ++;
		hist[tmp]++;
	}
	//绘制直方图
	DrawHistImg(hist);

	double d_hist[255];
	double sum_d_hist[255];
	sum_d_hist[0] = 0;
	for (int i = 0;i<255;i++)
	{
		d_hist[i] = hist[i]/(double)sumPixel;
		if (i>0)
		{
			sum_d_hist[i] = sum_d_hist[i-1]+d_hist[i];
		}
	}

	int new_hist[255];
	for (int i = 0;i<255;i++)
	{
		new_hist[i] = (int)(255*sum_d_hist[i]+0.5);
	}

	for (int i = 0;i<img.rows ;i++)
		for (int j = 0;j<img.cols ;j++)
		{
			int tmp = img.at<uchar>(i,j);
			img.at<uchar>(i,j) = new_hist[tmp];
		}
	DrawHistImg(new_hist);
		
	imshow("HistEqualize Image",img);
	waitKey(0);

	return 0;
}

效果如下

直方图显示效果不大好，仅供参照