opencv 绘制直方图学习过程中的疑问记录

一、 BGR与RGB

opencv中默认的通道存储顺序的是BGR，而不是熟知的RGB顺序，即BLUE,GREEN,RED，是按照此顺序排列的。例如在颜色表示函数scalar(a,b,c)中，scalar(255,0,0)表示纯蓝色，scalar(0,255,0)表示纯绿色。

二、 二维直方图distHist的意义

文章来源http://blog.csdn.net/xiaowei_cqu/article/details/7600666

首先来看一维直方图，一维直方图很好理解，横轴是bins，即每个特征空间子区段的数目，比如直方图的高度就是统计出来的数目，比如灰度图像的直方图，横坐标是各个“强度”的灰度值，纵坐标就是满足灰度值“强度”的像素点个数。

然后，高维直方图可以理解为图像在每个维度上灰度级分布的直方图。常见的是二维直方图。如红-蓝直方图的两个分量分别表示红光图像的灰度值和蓝光图像灰度值的函数。其图像坐标（Dr,Db）处表示在红光图像中具有灰度级Dr同时在蓝光图像中具有灰度级Db的像素个数。满足条件的像素个数的多少可以用灰度来表示，例如（binValue/maxValue）*255，则，像素个数越多，在图像上显示的颜色越白。

三、 rectangle()函数中的减法

以下图绘制一维灰度直方图为例，代码第30行计算了要绘制的高度intensity，第33行绘制直方图的矩形时，为什么要用矩形长条的高度减去intensity呢？因为histimage是一个矩阵，矩阵的元素排列顺序是从左上角开始的，Y轴方向向下，矩阵左上角是（0,0）位置，而左下角是（0,255）位置，绘制的直方图是X-Y轴的，Y轴方向向上，以左下角为（0,0）位置，左上角才是（0,255），说到底，是用矩阵绘制一个直方图，而不是直接生成一个直方图。

int main( )  
{  
    IplImage * src= cvLoadImage("baboon.jpg");  
    IplImage* gray_plane = cvCreateImage(cvGetSize(src),8,1);  
    cvCvtColor(src,gray_plane,CV_BGR2GRAY);  
  
    int hist_size = 256;    //直方图尺寸  
    int hist_height = 256;  
    float range[] = {0,255};  //灰度级的范围  
    float* ranges[]={range};  
    //创建一维直方图，统计图像在[0 255]像素的均匀分布  
    CvHistogram* gray_hist = cvCreateHist(1,&hist_size,CV_HIST_ARRAY,ranges,1);  
    //计算灰度图像的一维直方图  
    cvCalcHist(&gray_plane,gray_hist,0,0);  
    //归一化直方图  
    cvNormalizeHist(gray_hist,1.0);  
  
    int scale = 2;  
    //创建一张一维直方图的“图”，横坐标为灰度级，纵坐标为像素个数（*scale）  
    IplImage* hist_image = cvCreateImage(cvSize(hist_size*scale,hist_height),8,3);  
    cvZero(hist_image);  
    //统计直方图中的最大直方块  
    float max_value = 0;  
    cvGetMinMaxHistValue(gray_hist, 0,&max_value,0,0);  
      
    //分别将每个直方块的值绘制到图中  
    for(int i=0;i<hist_size;i++)  
    {  
        float bin_val = cvQueryHistValue_1D(gray_hist,i); //像素i的概率  
        int intensity = cvRound(bin_val*hist_height/max_value);  //要绘制的高度  
        cvRectangle(hist_image,  
            cvPoint(i*scale,hist_height-1),  
            cvPoint((i+1)*scale - 1, hist_height - intensity),  
            CV_RGB(255,255,255));    
    }  
    cvNamedWindow( "GraySource", 1 );  
    cvShowImage("GraySource",gray_plane);  
    cvNamedWindow( "H-S Histogram", 1 );  
    cvShowImage( "H-S Histogram", hist_image );  
  
    cvWaitKey(0);  
}