opencv之图像直方图

#include<opencv2/opencv.hpp>
#include<iostream>
#include<vector>
using namespace std;
using namespace cv;
int main()
{
	
	 const char* imagename = "C://Users//huashuo111//Desktop//test2.bmp";
     //从文件中读入图像
     Mat img = imread(imagename,IMREAD_GRAYSCALE);
     //如果读入图像失败
     if(img.empty())
     {
         fprintf(stderr, "Can not load image %s\n", imagename);
         return-1;
     }
	// resize(img,img,Size(),0.5,0.5);
	int histsize=255;
	 float range[]={0,256};
	
	 const float* histrange={range};
	 bool uniform=true;
	 bool accumulate=false;
	 Mat b_hist,g_hist,r_hist;
	 vector<Mat> channels;
	 split(img,channels);//图像分成三通道
	
	 calcHist(&channels[0],1,0,Mat(),b_hist,1,&histsize,&histrange,uniform,accumulate);
	
	 //&channels[0]:图像地址（蓝色）
	 //1：图像数量
	 //0:计算直方图的通道数
	 //Mat():可选项，像素个数
	 //B_hist:输出直方图
	 //1:直方图维数
	 //&histsize:每一维直方图数组大小
	 //&histrange:用来进行统计的范围。比如
     //float rang1[] = {0, 20};
     //float rang2[] = {30, 40};
     //const float *rangs[] = {rang1, rang2};那么就是对0，20和30，40范围的值进行统计。
	 //uniform:每一个竖条的宽度是否相等。
	 //accumulate:默认0,不累加
	 int hist_w=512;
	 int hist_h=400;
	 int bin_w=cvRound((double)hist_w/histsize);
	 Mat histImage(hist_h,hist_w,CV_8UC3,Scalar(0,0,0));
	 normalize(b_hist,b_hist,0,histImage.rows,NORM_MINMAX,-1,Mat());//归一化到[0,400],400是最后显示直方图高度
	
	 for(int i=1;i<histsize;i++)
		 //图像中原点在左上角
	 {
		line(histImage,Point(bin_w*(i-1),hist_h-cvRound(b_hist.at<float>(i-1))),Point(bin_w*(i),hist_h-cvRound(b_hist.at<float>(i))),Scalar(255,0,0),1,8,0);
		
	 }
	imshow("直方图",histImage);
	imwrite("C://Users//huashuo111//Desktop//直方图.bmp",histImage);
     //此函数等待按键，按键盘任意键就返回
     waitKey();
     return 0;
}

原图：

蓝色通道直方图：