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

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

一:一维直方图

Code:

[cpp]  view plain  copy
 print ? 在CODE上查看代码片 派生到我的代码片
  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 ? 在CODE上查看代码片 派生到我的代码片
  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/7600666cv】

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

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值