获取图像像素点的概率分布图

来源：互联网发布：apache 双向认证编辑：程序博客网时间：2024/05/07 10:44
获取图像像素点的概率分布，不细讲，看代码
#include<cv.h>  #include<highgui.h>  #define cvGetHistValue_1D( hist, idx0 ) ((float*)(cvPtr1D( (hist)->bins, (idx0), 0 )) int main(){int hist_size[] = { 16, 16, 16 };float range[] = { 0, 255 };float *ranges[] = { range };int i, j, bin_w;float max_value, min_value;int min_idx, max_idx;char *name[3] = { "b", "g", "r" };double mean = 0, variance = 0;IplImage* img = cvLoadImage("E:\\VS\\picture\\24.jpg", 1);IplImage *pImage = NULL;IplImage *pImg[3];pImg[0] = cvCreateImage(cvGetSize(img), 8, 1);pImg[1] = cvCreateImage(cvGetSize(img), 8, 1);pImg[2] = cvCreateImage(cvGetSize(img), 8, 1);cvSplit(img, pImg[0], pImg[1], pImg[2], NULL);pImage = cvCloneImage(img);CvRect rect = cvRect(0, 0, 500, 600);//cvSetImageROI(pImage, rect);//创建一个图像用来存放直方图  IplImage *histImage[3];CvHistogram *hist[3];for (j = 0; j<3; j++){histImage[j] = cvCreateImage(cvSize(320, 200), 8, 3);//cvZero(histImage[j]);  hist[j] = cvCreateHist(1, &hist_size[j], CV_HIST_ARRAY, ranges, 1);//计算直方图并作用到hist变量中  cvCalcHist(&pImg[j], hist[j], 0, NULL);//得到直方图的最值及标号  cvGetMinMaxHistValue(hist[j], &min_value, &max_value, &min_idx, &max_idx);//缩放其最大值和最小值让其融入图像  cvScale(hist[j]->bins, hist[j]->bins, ((double)histImage[j]->height) / max_value, 0);//设置所有的直方图的数值为255  cvSet(histImage[j], cvScalarAll(255), 0);//建一个比例因子  沿宽度释放  bin_w = cvRound((double)histImage[j]->width / hist_size[j]);mean = 0;for (i = 0; i<hist_size[j]; i++){CvScalar scalar = cvScalarAll(0);scalar.val[j] = (i * 255 / hist_size[j]);cvRectangle(histImage[j], cvPoint(i*bin_w, histImage[j]->height),cvPoint((i + 1)*bin_w, histImage[j]->height - cvRound(cvGetReal1D(hist[j]->bins, i))), scalar, -1, 8, 0);float *bins = cvGetHistValue_1D(hist[j], i));//***//增加均值  mean += bins[0];//std::cout<<bins[0]<<" "<<bins[1]<<std::endl;  //printf("%d %d\n",bins[0],bins[1]);  }cvNamedWindow(name[j], 0);cvShowImage(name[j], histImage[j]);mean /= hist_size[j];//根据均值计算变化量  for (i = 0; i<hist_size[j]; i++){float* bins = cvGetHistValue_1D(hist[j], i));//***variance += pow((bins[0] - mean), 2);}variance /= hist_size[j];printf("histgram Mean:%f\n", mean);}//创建窗口  cvNamedWindow("Original", 0);cvShowImage("Original", pImage);cvWaitKey(0);cvReleaseImage(&img);for (j = 0; j<3; j++){cvDestroyWindow(name[j]);cvReleaseImage(&pImg[j]);cvReleaseImage(&histImage[j]);cvReleaseHist(&hist[j]);}cvDestroyWindow("Original");return 0;}
0 0