openCV：三通道直方图计算

#include <iostream>  
#include<opencv2/opencv.hpp>
#include "math.h"  

using namespace std;
using namespace cv;

int main(int argc, char** argv)
{
	Mat src, dst;//初始化一个操作对象
	src = imread("C:/Users/JAY/Desktop/Others/lena.jpg");
	if (!src.data)//判断图片是否加载进来
	{
		cout << "不能加载图片" << endl;
		return -1;
	}
	//namedWindow("加载的图片", WINDOW_AUTOSIZE);
	imshow("加载的图片", src);//""内命名一致，才能显示在一个窗口

	//分通道显示
	vector<Mat> bgr_planes;
	split(src, bgr_planes);//把多通道图像分为多个单通道图像
	//imshow("单通道图像", bgr_planes[0]);

	//计算直方图
	int histSize = 256;
	float range[] = { 0,256 };
	const float *histRange = { range };
	//bool uniform = true;
	//bool accumulate = false;
	Mat b_hist, g_hist, r_hist;
	calcHist(&bgr_planes[0], 1, 0, Mat(), b_hist, 1, &histSize, &histRange, true, false);//计算直方图公式
	calcHist(&bgr_planes[1], 1, 0, Mat(), g_hist, 1, &histSize, &histRange, true, false);
	calcHist(&bgr_planes[2], 1, 0, Mat(), r_hist, 1, &histSize, &histRange, true, false);

	//归一化
	int hist_h = 400;
	int hist_w = 512;
	int bin_v = hist_w / histSize;
	Mat histImage(hist_w, hist_h, CV_8UC3, Scalar(0, 0, 0));
	normalize(b_hist, b_hist, 0, hist_h, NORM_MINMAX, -1, Mat());//归一化
	normalize(g_hist, g_hist, 0, hist_h, NORM_MINMAX, -1, Mat());//归一化
	normalize(r_hist, r_hist, 0, hist_h, NORM_MINMAX, -1, Mat());//归一化
	
	//绘制直方图
	for (int i = 1; i < histSize; i++)
	{	
		line(histImage, Point((i - 1)*bin_v, hist_h - cvRound(b_hist.at<float>(i - 1))), Point((i)*bin_v, hist_h - cvRound(b_hist.at<float>(i))), Scalar(255, 0, 0), 1, LINE_AA);
		line(histImage, Point((i - 1)*bin_v, hist_h - cvRound(g_hist.at<float>(i - 1))), Point((i)*bin_v, hist_h - cvRound(g_hist.at<float>(i))), Scalar(0, 255, 0), 1, LINE_AA);
		line(histImage, Point((i - 1)*bin_v, hist_h - cvRound(r_hist.at<float>(i - 1))), Point((i)*bin_v, hist_h - cvRound(r_hist.at<float>(i))), Scalar(0, 0, 255), 1, LINE_AA);
	}
	imshow("绘制直方图", histImage);

	waitKey(0);
	return 0;
}