opencv 2.3.1 绘制H-S直方图

LOVE青青小城

于 2012-07-11 14:44:02 发布

阅读量647

点赞数

分类专栏：图像视频分析

本文链接：https://blog.csdn.net/LQQXC/article/details/7735478

版权

图像视频分析专栏收录该内容

3 篇文章 0 订阅

订阅专栏

#include <cv.h>
#include <highgui.h>
#include <iostream>
using namespace std;
using namespace cv;
 
 
int main( int argc, char** argv )
{
	Mat src=imread("D:/lena.jpg");
	Mat hsv;
	cvtColor(src,hsv,CV_BGR2HSV);

	// let's quantize the hue to 30 levels
	// and the saturation to 32 levels

	int hbins = 30, sbins = 32;

	int histSize[] = {hbins, sbins};

	// hue varies from 0 to 179, see cvtColor

	float hranges[] = { 0, 180 };

	// saturation varies from 0 (black-gray-white) to

	// 255 (pure spectrum color)

	float sranges[] = { 0, 256 };

	const float* ranges[] = { hranges, sranges };

	MatND hist;

	// we compute the histogram from the 0-th and 1-st channels

	int channels[] = {0, 1};

	calcHist( &hsv, 1, channels, Mat(), hist, 2, histSize, ranges,true, false );
	normalize( hist,hist,1);

	double maxVal=0;

	minMaxLoc(hist, 0, &maxVal, 0, 0);

	int scale = 10;

	Mat histImg = Mat::zeros(sbins*scale, hbins*10, CV_8UC3);

	for( int h = 0; h < hbins; h++ )

	for( int s = 0; s < sbins; s++ )

	{


	float binVal = hist.at<float>(h, s);


	int intensity = cvRound(binVal*255/maxVal);


	rectangle( histImg, Point(h*scale, s*scale),Point( (h+1)*scale - 1, (s+1)*scale - 1),Scalar::all(intensity),CV_FILLED );

	}

	namedWindow( "Source", 1 );

	imshow( "Source", src );

	namedWindow( "H-S Histogram", 1 );

	imshow( "H-S Histogram", histImg );

	waitKey();

}