彩色图像直方图均衡处理

1、彩色图像直方图均衡处理彩色图像直方图均衡处理，在灰度图像直方图均衡处理的基础上实现彩色直方图均衡处理。三、彩色直方图均衡#includestdafx.h#include#includeint_tmain(intargc,_TCHAR*argv)IplImage*src=cvLoadImage(1.jpg);cvNamedWindow(，原图,1);cvShowImage(原图”,src);IplImage*equlimage=cvCreateImage(cvGetSize(src),8,3);分别均衡化每个信道IplImage*redlmage=cvCreateImage(cvGetSize(

2、src),8,1);IplImage*greenImage=cvCreateImage(cvGetSize(src),8,1);IplImage*blueImage=cvCreateImage(cvGetSize(src),8,1);cvSplit(src,blueImage,greenImage,redImage,NULL);cvEqualizeHist(redImage,redImage);cvEqualizeHist(greenImage,greenImage);cvEqualizeHist(blueImage,blueImage);均衡后的图像cvMerge(blueImage,gre

3、enImage,redImage,NULL,equlimage);cvMerge(blueImage,NULL,NULL,NULL,equlimage);cvMerge(NULL,greenImage,NULL,NULL,equlimage);cvMerge(NULL,NULL,redImage,NULL,equlimage);cvNamedWindow(均衡图,1);cvShowImage(均衡图”,equlimage);/*cvNamedWindow(均衡-红,1);cvShowImage(均衡-红”,redImage);cvNamedWindow(均衡-绿,1);cvShowImage(

4、均衡-绿”,greenImage);cvNamedWindow(均衡-蓝,1);cvShowImage(均衡-蓝”,blueImage);*/IplImage*gray_plane=cvCreateImage(cvGetSize(src),8,1);IplImage*equl_plane=cvCreateImage(cvGetSize(src),8,1);cvCvtColor(src,gray_plane,CV_BGR2GRAY);cvCvtColor(equlimage,equl_plane,CV_BGR2GRAY);inthist_size=256;/直方图尺寸inthist_height

5、=256;floatrange=0,255;/灰度级的范围float*ranges=range;创建一维直方图，统计图像在0,255像素的均匀分布CvHistogram*gray_hist=cvCreateHist(1,&hist_size,CV_HIST_ARRAY,ranges,1);CvHistogram*equl_hist=cvCreateHist(1,&hist_size,CV_HIST_ARRAY,ranges,1);计算灰度图像的一维直方图cvCalcHist(&gray_plane,gray_hist,0,0);cvCalcHist(&equl_plane,equl_hist,

6、0,0);归一化直方图归一化直方图cvNormalizeHist(gray_hist,1.0);cvNormalizeHist(equl_hist,1.0);intscale=1;创建一张一维直方图的“图像，横坐标为灰度级，纵坐标为像素个数(*scale)图像位深度为8位整形，每个元素通道号为3IplImage*hist_image=cvCreateImage(cvSize(hist_size*scale*2,hist_height),8,3);初始化，图像清零cvZero(hist_image);统计直方图中的最大直方块floatmax_value=0;floatmax_value2=0;c

7、vGetMinMaxHistValue(gray_hist,O,&max_value,0,0);cvGetMinMaxHistValue(equl_hist,O,&max_value2,0,0);分别将每个直方块的值绘制到图中for(inti=O;ihist_size;i+)floatbin_val=cvQueryHistValue_1D(gray_hist,i);floatbin_val2=cvQueryHistValue_1D(equl_hist,i);intintensity=cvRound(bin_val*hist_height/max_value);intintensity2=cvRound(bin_val2*hist_height/max_value2);cvRectangle(hist_image,cvPoint(i*scale,hist_height-1),cvPoint(i+1)*scale-1,hist_height-intensity),CV_RGB(255,255,255);cvRectangle(hist_image,cvPoint(i+hist_height)*scale,hist_height-1),cvPoint(i+hist_height+1)*scale-1,hist_height-intensity2),CV_RGB(255,255,255)