图像处理2：图像直方图（c++）

计算图像的直方图

要注意opencv中起始点是左上角是（0,0）为了符合常见的直观，要进行上下镜像变换。

void imagehist(Mat src)//, Mat dst
{
	Mat EH = Mat::zeros(src.size(), CV_8UC1);
	double hist[3][256];
	for (int i = 0; i < 3; i++)
	{
		for (int j = 0; j < 256; j++)
		{
			hist[i][j] = 0;
		}
	}

	for (int k = 0; k < 3; k++)
	{
		for (int i = 0; i < src.rows; i++)
		{
			for (int j = 0; j < src.cols; j++)
			{
				hist[k][src.at<Vec3b>(i, j)[k]]++;
			}
		}
	}
	char color[3] = {'B','G','R'};
	for (int k = 0; k < 3; k++)
	{
		Mat frame(257, 256, CV_8U, Scalar(0));
		for (int i = 0; i < 256; i++)
		{
			double max_value = max_double(hist[k], 256); //
			//cout << "Value of str is : " << endl;
			int pointend = cvRound(hist[k][i] / max_value * 256 * 0.9);
			line(frame, Point(i, 0+256), Point(i, 256-pointend), Scalar(255, 255, 255), 2);//
			//line(frame, (i,0),(i,hist[0][i]/ src.rows/ src.cols), Scalar(0, 0, 255), 2);
		}
		char color_j = color[k];
        imshow("灰度直方图", frame);
		waitKey(0);
	}
}

void CMFCworkDlg::OnBnClickedButton1()
{
	// TODO: 在此添加控件通知处理程序代码
	Mat image = imread("ema.jpg", 1);

	imagehist(image);
	//equalizeHist(image);

}