Opencv2系列学习笔记4(灰度直方图)

和之前的一篇blog【http://blog.csdn.net/lu597203933/article/details/14104505】一样，这篇是opencv2的Mat格式。

一：一维直方图

Code：

[cpp]  view plain  copy

 print ?

1. int main()  
2. {  
3.     Mat image = imread("F:\\huangrong.jpg", 0);  
4.     if(!image.data)  
5.     {  
6.         cout << "fail to load image" << endl;  
7.         return 0;  
8.     }  
9.     MatND hist;       // 在cv中用CvHistogram *hist = cvCreateHist  
10.     int dims = 1;  
11.     float hranges[] = {0, 255};  
12.     const float *ranges[] = {hranges};   // 这里需要为const类型  
13.     int size = 256;  
14.     int channels = 0;  
15.     // 计算图像的直方图  
16.     calcHist(&image, 1, &channels, Mat(), hist, dims, &size, ranges);    // cv 中是cvCalcHist  
17.     int scale = 1;  
18.     Mat imageShow(size * scale, size, CV_8U, Scalar(0));  
19.     // 获取最大值和最小值  
20.     double minVal = 0;  
21.     double maxVal = 0;  
22.     minMaxLoc(hist,&minVal, &maxVal, 0, 0);  //  cv中用的是cvGetMinMaxHistValue  
23.     //显示直方图的图像  
24.     int hpt = saturate_cast<int>(0.9 * size);  
25.       
26.     for(int i = 0; i < 256; i++)  
27.     {  
28.         float value = hist.at<float>(i);           //   注意hist中是float类型    cv中用cvQueryHistValue_1D  
29.         int realValue = saturate_cast<int>(value * hpt/maxVal);  
30.         //line(imageShow, Point(i, size - 1), Point(i, size - realValue), Scalar(0));  
31.         rectangle(imageShow,Point(i*scale, size - 1), Point((i+1)*scale - 1, size - realValue), Scalar(255));  
32.     }  
33.     namedWindow("showImage");  
34.     imshow("showImage", imageShow);  
35.     waitKey(0);  
36.     return 0;  
37. }  

Explaination：

<1>代码中有注释将mat格式与cv格式进行了比对，大家注意

<2> float value = hist.at<float>(i);

       要注意直方图数据为浮点类型

<3>画出直方图其实有多种方式，可以使用line或rectangle都可以

Result：

 

二：二维直方图

Code：

[cpp]  view plain  copy

 print ?

1. int main()  
2. {  
3.     Mat image = imread("F:\\baboon.jpg", 1);  
4.     if(!image.data)  
5.     {  
6.         cout << "fail to load image" << endl;  
7.         return 0;  
8.     }  
9.     MatND hist;       // 在cv中用CvHistogram *hist = cvCreateHist  
10.     int dims = 2;  
11.     float r_hranges[] = {0, 255};  
12.     float g_hranges[] = {0, 255};  
13.     const float *ranges[] = {r_hranges, g_hranges};   // 这里需要为const类型  
14.     int size[2] = {256, 256};  
15.     int channels[] ={0,1};   //代表 r g通道 2代表b通道  
16.     // 计算图像的直方图  
17.     calcHist(&image, 1, channels, Mat(), hist, dims, size, ranges);    // cv 中是cvCalcHist  
18.     int scale = 1;  
19.     Mat imageShow(size[0] * scale, size[1]*scale, CV_8UC3, Scalar(255));  
20.     // 获取最大值和最小值  
21.     double minVal = 0;  
22.     double maxVal = 0;  
23.     minMaxLoc(hist,&minVal, &maxVal, 0, 0);  //  cv中用的是cvGetMinMaxHistValue  
24.     //显示直方图的图像  
25.     //int hpt = saturate_cast<int>(0.9 * size[]);  
26.       
27.     for(int i = 0; i < 256; i++)  
28.     {  
29.         for(int j = 0; j < 256; j++)  
30.         {  
31.         float value = hist.at<float>(i,j);           //   注意直方图的值是float类型    cv中用cvQueryHistValue_1D  
32.         int realValue = saturate_cast<int>(value * size[0]/maxVal);  
33.         //line(imageShow, Point(i, size - 1), Point(i, size - realValue), Scalar(0));  
34.         rectangle(imageShow,Point(i*scale, j*scale), Point( ((i+1)*scale - 1), (j+1)*scale - 1), Scalar(realValue,realValue, realValue));  
35.         }  
36.     }  
37.     namedWindow("baboon");  
38.     namedWindow("showImage");  
39.     imshow("showImage", imageShow);  
40.     imshow("baboon", image);  
41.     waitKey(0);  
42.     return 0;  
43.   
44. }  

Explaination：

需要注意的是取二维的通道数是通过变量channels来取得的，不需要将图片再进行split

Result：

另外推荐小魏的blog，关于直方图的描述写的也非常好：

链接：http://blog.csdn.net/xiaowei_cqu/article/details/8833799【Mat格式】

http://blog.csdn.net/xiaowei_cqu/article/details/7600666【cv】

作者：小村长  出处：http://blog.csdn.net/lu597203933 欢迎转载或分享，但请务必声明文章出处。 （新浪微博：小村长zack, 欢迎交流！）