直方图均衡化

#include "highgui.h"
#include "cv.h"
#include <iostream>
#include <iomanip>
using namespace std;
#define HDIM 256 //直方图的维数
int main()
{
IplImage* src = cvLoadImage("lena.jpg",0);
IplImage* dst = NULL;
CvHistogram* hist =0;
int n = HDIM;
double nn[HDIM];
uchar T[HDIM];
CvMat *T_mat;

int x;
int sum;//图像中像素总和
double val = 0;
cvNamedWindow("source",1);
cvNamedWindow("result",1);
//计算直方图
hist = cvCreateHist(1,&n,CV_HIST_ARRAY,0,1);
cvCalcHist(&src,hist,0,0);
val = 0;
for (x=0;x<n;x++)
{
val = val + cvGetReal1D(hist->bins,x);
nn[x] = val;
}
//归一化直方图
sum = src->height*src->width;
for (x=0;x<n;x++)
{
T[x] = (uchar)(255*nn[x]/sum);
}
dst = cvCloneImage(src);
T_mat = cvCreateMatHeader(1,256,CV_8UC1);
cvSetData(T_mat,T,0);
//直接调用内部函数完成look-up-table的过程
cvLUT(src,dst,T_mat);
cvShowImage("source",src);
cvShowImage("result",dst);
cvWaitKey(0);
cvReleaseImage(&src);
cvReleaseImage(&dst);
cvDestroyAllWindows();
}