VC++实现一维图像直方图计算与显示

int  getHistogram1DImage(IplImage * pImgGray)
{
     int histSize=255; //直方图针数
 float range0[]={0,256}; //第0维数值变化范围
 float * ranges[]={range0}; //第1维数值变化范围
 int i,binW;    //下标号
 float maxValue,minValue; //直方图数值的max和min
 int minIdx=0,maxIdx=0; //对应上述值时的下标号
 double mean=0,variance=0; //均值和方差
 
//     CvRect rect=cvRect(0,0,256,256);
//  cvSetImageROI(pImgGray,rect);    //设置ROI
 
   IplImage * histImg=cvCreateImage(cvGetSize(pImgGray),IPL_DEPTH_8U,1);//用于存放直方图
   CvHistogram * hist=cvCreateHist(1,&histSize,CV_HIST_ARRAY,ranges,1);//创建一张一维直方图的“图”，横坐标为灰度级，纵坐标为像素个数（*scale）  
   cvCalcHist(&pImgGray,hist,0,NULL);//计算直方图
   cvGetMinMaxHistValue(hist,&minValue,&maxValue,&minIdx,&maxIdx);//得到最大最小值及其标号，最大值即在某处标号处像素的最大个数
   //缩放最大值最小值以溶入图像
   cvScale(hist->bins,hist->bins,((double)histImg->height)/maxValue,0);
   cvSet(histImg,cvScalar(255),0);
   binW=cvRound((double)histImg->width/histSize);//　对一个double型的数进行四舍五入，并返回一个整型数！
 //把直方图画到图像中
   for (i=0;i<histSize;i++)
{
   cvRectangle(histImg,cvPoint(i*binW,histImg->height), cvPoint((i+1)*binW,histImg->height-cvRound(cvGetReal1D(hist->bins, i))),
   //过对角线上的两个顶点绘制简单、指定粗细或者带填充的矩形,cvGetReal1D返回单通道数组的指定元素
   cvScalarAll(0),-1,8,0);
   float * bins=cvGetHistValue_1D(hist,i);//像素i的概率 即灰度是i的像素个数
   mean+=bins[0];
 }