VC++直方图均衡化显示图像及直方图

#include <cv.h> #include <highgui.h> #include <cxcore.h> #include <stdio.h> #define HDIM 256 int main(int argc,char** argv) { IplImage *dst=0,*src=0; CvHistogram *histogram; CvHistogram *HistEqualization; IplImage *histOriginal;//原图像的直方图 IplImage *HistEquAfter;//均衡化后直方图 float HistogramRange1[2]={0,255}; float *HistogramRange[1]={&HistogramRange1[0]}; int n=HDIM; double nn0[HDIM]; int srcnn00[HDIM]; uchar T0[HDIM]; CvMat *T_mat0; int i=0; int sum = 0; //sum of pixels of the source image double val0=0; int val1=0; if((src=cvLoadImage("4.jpg",0)) == NULL) // force to gray image { printf("Read Image Error!"); return -1; } //原图像的直方图 histOriginal = cvCreateImage(cvSize(300,400),8,3); histOriginal->origin=1; histogram = cvCreateHist(1,&n,CV_HIST_ARRAY,HistogramRange,1); cvCalcHist(&src,histogram,0,0); // Compute intensity transformation sum = src->height * src->width; for( i = 0; i < n; i++ ) { val0 = val0 + cvGetReal1D (histogram->bins, i); nn0[i]=val0; T0[i] = (uchar) (255*nn0[i]/sum); // range is [0,255] srcnn00[i]=(int)(nn0[i]*255/sum); cvLine(histOriginal,cvPoint(i,0),cvPoint(i,(int)(cvQueryHistValue_1D(histogram,i)/3)),CV_RGB(127,127,127)); } // 直方图均衡化 dst = cvCloneImage(src); T_mat0 = cvCreateMatHeader(1,256,CV_8UC1); cvSetData(T_mat0,T0,src->width ); cvLUT( src, dst, T_mat0); //均衡化后的直方图 HistEquAfter = cvCreateImage(cvSize(300,400),8,3); HistEquAfter->origin=1; HistEqualization = cvCreateHist(1,&n,CV_HIST_ARRAY,HistogramRange,1); cvCalcHist(&dst,HistEqualization,0,0); for (i=0;i<n;i++) { cvLine(HistEquAfter,cvPoint(i,0),cvPoint(i,(int)(cvQueryHistValue_1D(HistEqualization,i)/3)),CV_RGB(127,127,127)); } cvNamedWindow( "source", 1 ); cvNamedWindow( "HistEquImage", 1 ); cvNamedWindow( "S_Histogram", 1 ); cvNamedWindow( "HistEqualization", 1 ); cvShowImage( "source", src );//原图像 cvShowImage( "HistEquImage", dst );//均衡化后图像 cvShowImage("S_Histogram",histOriginal);//原图像直方图 cvShowImage("HistEqualization",HistEquAfter);//均衡化后直方图 cvWaitKey(-1); cvReleaseImage( &src ); cvReleaseImage( &dst ); cvReleaseImage(&histOriginal); cvReleaseImage(&HistEquAfter); cvReleaseHist ( &histogram); cvReleaseHist(&HistEqualization); cvDestroyWindow("source"); cvDestroyWindow("HistEquImage"); cvDestroyWindow("S_Histogram"); cvDestroyWindow("HistEqualization"); return 0; }