颜色直方图

////////////////////////////////////////////////

//////////////////////////////////////////////////这是一个计算颜色直方图的程序，可以计算出一幅多通道图像的颜色直方图，但是目前还无法显示。

#include <StdAfx.h>
#include <cv.h>
#include <highgui.h>
#include <cxcore.h>

void main( int argc, char** argv )
{
 IplImage *src = 0;
 IplImage *histImg = 0;
 CvHistogram *hist = 0;
 int hdims[] = {256,256,256}; // 划分hist的个数，
// float hranges_arr[] = {0,256};
// float *hranges = hranges_arr;
 float r_ranges[] = {0,255};
 float g_ranges[] = {0,255};
 float b_ranges[] = {0,255};
 float* hranges[] = {r_ranges,g_ranges,b_ranges};
 int  aa=255;
 int fdims[3]={aa,aa,aa};
 float *pranges[]={r_ranges,g_ranges,b_ranges};

 int i;

 // 显示原图
 src = cvLoadImage("C://test.bmp",3);
 cvNamedWindow( "Histogram", 1 );

 //创建直方图 c
 hist = cvCreateHist(3,hdims,CV_HIST_ARRAY,hranges,1);
 //计算直方图
 IplImage *pold_B=cvCreateImage(cvGetSize(src),8,1);
 IplImage *pold_G=cvCreateImage(cvGetSize(src),8,1);
 IplImage *pold_R=cvCreateImage(cvGetSize(src),8,1);
 IplImage *pall[]={pold_R,pold_G,pold_B};
 cvSplit(src,pold_B,pold_G,pold_R,NULL);

 cvCalcHist( pall, hist, 0, 0 );

// int r,g,b;
// for(i=0;i<=255;i++)
// cvGetHistValue_3D(hist,r,g,b);

 cvShowImage( "Histogram", src );
 cvWaitKey(0);

}
 

 

///////////////////////////////////////////////////////////////

以上为创建3通道直方图的代码，下面为单通道直方图的计算及相关的显示，和原图的显示代码

//////////////////////////////////////////////////////////////////////////
// 另外的显示图象的方法
//////////////////////////////////////////////////////////////////////////
CPaintDC dc(GetDlgItem(IDC_STATIC_KEY)); //用此句,得到picture控件的DC,图像将被绘制在控件上 
    CRect rcclient;
    GetDlgItem(IDC_STATIC_KEY)->GetClientRect(&rcclient);
    CDC memdc;
    memdc.CreateCompatibleDC(&dc); 
    CBitmap bitmap;

 hbmp = (HBITMAP)::LoadImage(AfxGetInstanceHandle(),
        "c://test.bmp", IMAGE_BITMAP, 0, 0, LR_CREATEDIBSECTION|LR_LOADFROMFILE);
 
 if( hbmp == NULL )
        return ;

 bitmap.Attach( hbmp );

//    bitmap.CreateCompatibleBitmap(&dc, rcclient.Width(), rcclient.Height());
    memdc.CreateCompatibleDC(&dc);
 memdc.SelectObject( &bitmap );
  //  CWnd::DefWindowProc(WM_PAINT, (WPARAM)memdc.m_hDC , 0);
 memdc.BitBlt( 0, 0, rcclient.Width(), rcclient.Height(), &memdc,
        rcclient.left, rcclient.top, SRCCOPY);
/*    CDC maskdc;

    maskdc.CreateCompatibleDC(&dc);

    CBitmap maskbitmap;

    maskbitmap.CreateBitmap(rcclient.Width(), rcclient.Height(), 1, 1, NULL);

    maskdc.SelectObject( &maskbitmap );

    maskdc.BitBlt( 0, 0, rcclient.Width(), rcclient.Height(), &memdc,

        rcclient.left, rcclient.top, SRCCOPY);
*/
  //////////////////////////////////////////////////////////////////////////
  /// 以下为显示图片和直方图，及其它
  //////////////////////////////////////////////////////////////////////////
/* CPaintDC dc(this);// device context for painting
 //GetDlgItem(IDC_STATIC_KEY)->GetClientRect(&m_ClientRect);
 GetDlgItem(IDC_STATIC_KEY)->GetWindowRect(&m_ClientRect);
 ScreenToClient(&m_ClientRect);
// CClientDC dc(GetDlgItem(IDC_STATIC_KEY));

 CDC *memDC = new CDC;
// memDC = new CDC;
  memDC->CreateCompatibleDC(&dc);
  memDC->SelectObject(bm);
 dc.StretchBlt(m_ClientRect.left,m_ClientRect.top,m_ClientRect.Width(),m_ClientRect.Height(),
   memDC,0,0,bmMetric.bmWidth,bmMetric.bmHeight,SRCCOPY);
*/
  /* if(bm->m_hObject!=NULL)
 {
  bm->GetBitmap(&bmMetric);

 MessageBox("bmfeikong");
 }
*/ //DeleteDC(memDC);

/* CRect m_ZhiFangTu;
 CString out;
 CBitmap bm,*pbm;
 BITMAP bmMetric;
 IplImage *src= 0;
    CvHistogram *hist = 0;
   
    int hdims = 50;     // 划分HIST的个数，越高越精确
    float hranges_arr[] = {0,255};
    float* hranges = hranges_arr;
    float max_val;
    int i;
 

// m_pWnd=GetDlgItem(IDC_STATIC_KEY);
 GetDlgItem(IDC_STATIC_KEY)->GetClientRect(&m_ClientRect);
 CClientDC dc1(GetDlgItem(IDC_STATIC_KEY));

 GetDlgItem(IDC_ZHIFANGTU)->GetClientRect(&m_ZhiFangTu);
 CClientDC dc2(GetDlgItem(IDC_ZHIFANGTU));

 CDC memDC;
 if(m_FileName !="")
 {
// m_FileName = "c://test.bmp";
 HBITMAP hbmp = (HBITMAP)::LoadImage(AfxGetInstanceHandle(),
        m_FileName, IMAGE_BITMAP, 0, 0, LR_CREATEDIBSECTION|LR_LOADFROMFILE);
    if( hbmp == NULL )
   MessageBox("加载位图失败！");

    bm.Attach( hbmp );
 bm.GetBitmap(&bmMetric);
 
  memDC.CreateCompatibleDC(&dc1);
  pbm=memDC.SelectObject(&bm);

  dc1.StretchBlt(0,0,m_ClientRect.right,m_ClientRect.bottom,
   &memDC,0,0,bmMetric.bmWidth,bmMetric.bmHeight,SRCCOPY);
  memDC.SelectObject(pbm);
  bm.DeleteObject();
  memDC.DeleteDC();
  
 //以上显示图片，下面计算颜色直方图，并且在框中显示
 // 首先 得到框的位置

 // 计算直方图
 src = cvLoadImage(m_FileName, 0);
    hist = cvCreateHist( 1, &hdims, CV_HIST_ARRAY, &hranges, 1 );  // 计算直方图
    cvCalcHist( &src, hist, 0, 0 ); // 计算直方图
    cvGetMinMaxHistValue( hist, 0, &max_val, 0, 0 );  // 只找最大值
    cvConvertScale( hist->bins, hist->bins, max_val ? 255. / max_val : 0., 0 ); // 缩放 bin 到区间 [0,255]
   
 for( i = 0; i < hdims; i++ )
    {
        double val = (cvGetReal1D(hist->bins,i)*m_ZhiFangTu.Height()/255 );
  dc2.MoveTo(i*m_ZhiFangTu.Width()/50,m_ZhiFangTu.bottom);
  dc2.LineTo(i*m_ZhiFangTu.Width()/50,m_ZhiFangTu.Height()-val);
 }

 }
 else
  MessageBox("m_FileName 为空");

*/

//////////////////////////////////////////////////////////////////////////
// 以上为显示图片和直方图
//////////////////////////////////////////////////////////////////////////