Opencv2系列学习笔记6(直方图的应用)

前面我讲述了直方图的概念以及如何用opencv实现一维和二维的直方图。详见这两篇blog: Mat 格式：http://blog.csdn.net/lu597203933/article/details/16884409  

cv：http://blog.csdn.net/lu597203933/article/details/14104505 

这次主要想讲点直方图的应用，其中包括使用查找表修改图像的外观、直方图的均衡化、反投影直方图检测特定图像的内容、meanshift算法<均值漂移>跟踪物体和利用图像直方图检索相似图像<可靠性比较低>。

一：使用查找表修改图像的外观

查找表是简单的一对一(或者多对一)函数，定义了如何将像素值转换为新像素值。Opencv的cv::LUT对图像应用查找表以生成新图像。公式为: result[i] = lookup[image[i]].

<1>图像反转

Code：

Main.cpp

[cpp]  view plain  copy  print ?

1.        MatND hist = getHistogram(image);  
2. int i;  
3. for(i = 0; i < 256; i++)  
4. {  
5.     cout << "Value" << i << ": " << hist.at<float>(i) << endl;  
6. }  
7. /*int dim(256);  
8. Mat lut(1, &dim, CV_8U);  
9. for(i = 0; i < 256; i++)  
10. {  
11.     lut.at<uchar>(i) = 255 - i;   
12. }  
13. Mat result = applyLookUp(image, lut);  

Result：

 

<2>提高图像对比度

以下代码定义了一个提高图像对比度的查找表。方法就是通过拉伸直方图，能得到扩展后的对比度。

Code：

得到图像的直方图：

[cpp]  view plain  copy  print ?

1. MatND getHistogram(Mat &image)  //  得到图像的直方图  
2. {  
3.     MatND hist;  
4.     int channels[] = {0};  
5.     int dims = 1;  
6.     int histSize[] = {256};   
7.     float granges[] = {0, 255};  
8.     const float *ranges[] = {granges};  
9.     calcHist(&image, 1, channels, Mat(), hist, dims, histSize, ranges);  
10.     return hist;  
11. }  

将图像的直方图展示出来：

[cpp]  view plain  copy  print ?

1. Mat getHistogramImage(Mat &image, int scaleX = 1, int scaleY = 1)   //  将图像的直方图展示出来  
2. {  
3.     MatND hist = getHistogram(image);  
4.     Mat showImage(scaleY * 256, scaleX*256, CV_8U,Scalar(255));  
5.     int i;  
6.     double maxValue = 0;  
7.     minMaxLoc(hist, 0, &maxValue, 0, 0);  
8.     for(i = 0; i < 256; i++)  
9.     {  
10.         float value = hist.at<float>(i);  
11.         int intensity = saturate_cast<int>(256*scaleY - 256*scaleY * (value/maxValue));  
12.         rectangle(showImage, Point(i*scaleX, scaleY*256 - 1), Point((i+1)*scaleX-1, intensity), Scalar(0));  
13.     }  
14.     return showImage;  
15. }  

得到并应用查找表：

[cpp]  view plain  copy  print ?

1. Mat stretch(Mat &image, int minValue = 0)  
2. {  
3.     // 首先计算直方图  
4.     MatND hist = getHistogram(image);  
5.     int imin = 0,i = 0;  
6.     // 寻找直方图的左端  
7.     for(; i < 256; i ++)  
8.     {  
9.         float iminValue = hist.at<float>(i);  
10.         if(iminValue > minValue)  
11.         {  
12.             imin = i;  
13.             break;  
14.         }  
15.     }  
16.     //  寻找直方图的右端  
17.     int imax = 255;  
18.     for(i = 255; i >= 0 ; i--)  
19.     {  
20.         float imaxValue = hist.at<float>(i);  
21.         if(imaxValue > minValue){  
22.             imax = i;  
23.             break;  
24.         }  
25.     }  
26.     // 创建查找表  
27.     Mat lookUp(1, 256, CV_8U);  
28.     // 填充查找表  
29.     for(i = 0; i < 256; i++)  
30.     {  
31.         // 确保数值位于imin和imax之间  
32.         if(i < imin)  
33.             lookUp.at<uchar>(i) = 0;  
34.         else if(i > imax)  
35.             lookUp.at<uchar>(i) = 255;  
36.         else     // 线性映射  
37.             lookUp.at<uchar>(i) = saturate_cast<uchar>(255.0*(i - imin)/(imax - imin));   
38.     }  
39.     // 应用查找表  
40.     Mat result;  
41.     LUT(image, lookUp, result);  
42.     return result;  
43. }  
44. Mat applyLookUp(Mat &image, Mat& lut)  
45. {  
46.     Mat result;  
47.     LUT(image, lut, result);  
48.     return result;  
49. }  

Main.cpp

[cpp]  view plain  copy  print ?

1. int main()  
2. {  
3.     Mat image = imread("F:\\group.jpg", 0);  
4.     if(!image.data)  
5.     {  
6.         cout << "fail to load image" << endl;  
7.         return 0;  
8.     }  
9.     Mat showImage = getHistogramImage(image);  
10.     Mat result = stretch(image, 100);  
11.     Mat showImage2 = getHistogramImage(result);  
12.     namedWindow("image");  
13.     imshow("image", image);  
14.     namedWindow("showImage");  
15.     imshow("showImage", showImage);  
16.     namedWindow("showImage2");  
17.     imshow("showImage2", showImage2);  
18.     namedWindow("result");  
19.     imshow("result", result);  
20.     waitKey(0);  
21.     return 0;  
22. }  

Result：

 

二：直方图均衡化

直方图均衡化的目的也是提高图像的对比度，但它与通过查找表提高对比度方法不同。通过拉伸直方图，使它覆盖所有的取值范围。这个策略确实可以简单有效地提升图像质量。然而，大多数情况下，图像在视觉上的缺陷并非源于过窄的强度范围，而是由于某些颜色值出现的频率高于另外一些。事实上，一副高质量的图片应该均衡使用所有的像素值。

Code：

Main.cpp

[cpp]  view plain  copy  print ?

1. int main()  
2. {  
3.     Mat image = imread("F:\\group.jpg", 0);  
4.     if(!image.data)  
5.     {  
6.         cout << "fail to load image" << endl;  
7.         return 0;  
8.     }  
9.     Mat result;  
10.     equalizeHist(image, result);   // 直方图的均衡化====增加图像的对比度  
11.     namedWindow("image");  
12.     imshow("image", image);  
13.     namedWindow("result");  
14.     imshow("result", result);  
15.     Mat imageHist = getHistogramImage(image);  
16.     Mat resultHist = getHistogramImage(result);  
17.   
18.     namedWindow("imageHist");  
19.     imshow("imageHist", imageHist);  
20.     namedWindow("resultHist");  
21.     imshow("resultHist", resultHist);  
22.   
23.     waitKey(0);  
24.     return 0;  
25. }  

Explaination：

opencv2中提供的函数为equalizeHist();它的原理是这样的，在一副完全均衡的直方图中，所有的容器拥有同等数量的像素。公式为：

       lookup.at<uchar>(i)= static_cast<uchar>(255*p[i])；；其中p[i]代表强度值小于i的像素所占的比例。

Result：

三：反投影直方图检测特定图像的内容

如果一副图像的区域中显示的是一种独特的纹理或者是一个独特的物体，那么这个区域的直方图可以看做一个概率函数，它给出了某个像素属于该纹理或者物体的概率。直白的说，就是在一副图像中将与某特定感兴趣的区域全部找出来，使用的是概率。

Code：

Main.cpp：

[cpp]  view plain  copy  print ?

1. int main()  
2. {  
3.     Mat image = imread("F:\\clound.jpg", 0);  
4.     if(!image.data)  
5.     {  
6.         cout << "Fail to load image" << endl;  
7.         return 0;  
8.     }  
9.     Mat imageROI = image(Rect(80, 120, 40, 50));  
10.     MatND histImage = getHistogram(imageROI);  
11.     normalize(histImage, histImage, 1.0);   //归一化该直方图  
12.     int channels[] = {0};  
13.     Mat result;  
14.     float granges[] = {0, 255};  
15.     const float *ranges[] = {granges};  
16.     calcBackProject(&image, 1, channels, histImage, result, ranges, 255.0);   
17.     /*第二个参数：图像数量   第三个参数: 通道数， histImage: 归一化的直方图,  result: 概率映射图 
18.     ranges: 维度的范围 
19.     */  
20.     namedWindow("image", 0);  
21.     imshow("image", image);  
22.     namedWindow("result", 0);  
23.     imshow("result", result);  
24.     waitKey(0);  
25.     return 0;  
26. }  

Result：

四：使用meanshift算法<均值漂移>跟踪物体

Meanshift算法的具体介绍可以看我转载的一篇blog：http://blog.csdn.net/lu597203933/article/details/17042331  它就是以迭代的方式锁定概率函数的局部最大值

主要分为这三个步骤：首先需要得到目标的直方图，其次用该直方图对图片进行反投影，后对反投影的图像进行meanshift算法。其中迭代的次数和精度作为收敛的条件。

Code：

[cpp]  view plain  copy  print ?

1. /*计算一幅图像特定区域的hist*/  
2. MatND getHueHistogram(Mat &image, int minSaturation = 0)  
3. {  
4.     MatND hist;  
5.     Mat hsvColor;  
6.     cvtColor(image, hsvColor, CV_BGR2HSV);  
7.     vector<Mat> hsv;   
8.     split(hsvColor, hsv);  
9.     Mat mask;  
10.     if(minSaturation > 0)  
11.     {  
12.         threshold(hsv[1], mask, minSaturation, 255, THRESH_BINARY);  
13.         int channels[] = {0};  
14.         int histSize[] = {181};  
15.         float hRanges[] = {0, 180};  
16.         const float *ranges[] = {hRanges};  
17.         int dims = 1;  
18.         calcHist(&hsv[0], 1, channels, mask, hist, dims, histSize, ranges);  
19.     }  
20.     normalize(hist, hist, 1.0);  
21.     return hist;  
22. }  
23.   
24. /*通过上面特定区域的直方图得到图像的反投影直方图*/  
25. Mat getBackProject(Mat &image, MatND colorHist)  
26. {  
27.     int channels[] = {0};  
28.     Mat result;  
29.     float hRanges[] = {0, 180};  
30.     const float *ranges[] = {hRanges};  
31.     calcBackProject(&image, 1, channels, colorHist, result, ranges, 255);  
32.     int thre = 60;  
33.     threshold(result, result, 60, 255,THRESH_BINARY);  
34.     return result;  
35. }  
36. int main()  
37. {  
38.     // 载入第一幅图像  
39.     Mat image = imread("F:\\4.jpg",1);  
40.     // 定义需要跟踪的目标---感兴趣的区域  
41.     Mat imageROI = image(Rect(230, 320, 60, 35));  
42.     // 得到目标的直方图  ********  
43.     int minSat = 65;  
44.     MatND colorHist = getHueHistogram(imageROI, minSat);  
45.   
46.     //  载入需要跟踪的图片  
47.     Mat image2 = imread("F:\\10.jpg", 1);  
48.     // 转换颜色到hsv空间  
49.     Mat hsvColor;  
50.     cvtColor(image2, hsvColor, CV_BGR2HSV);  
51.     vector<Mat> hsv;  
52.     split(image2, hsv);  
53.     // 得到目标图像的反投影直方图  *********  
54.     Mat result = getBackProject(hsv[0], colorHist);  
55.     // 去除 低饱和区域  
56.     threshold(hsv[1], hsv[1], minSat, 255, THRESH_BINARY);  
57.     bitwise_and(result, hsv[1], result);  
58.   
59.     Rect rect(230, 320, 60, 35);  
60.     rectangle(image, rect, Scalar(0,0,255));  
61.     rectangle(image2, rect, Scalar(0,0,255));  
62.     // 定义迭代的次数以及迭代的精度  
63.     TermCriteria criteria(TermCriteria::MAX_ITER, 10, 0.01);  
64.     // meanshift 算法 跟踪目标 并且得到新目标的位置rect  
65.     meanShift(result, rect, criteria);  
66.     rectangle(image2, rect, cv::Scalar(0,255,0));  
67.     namedWindow("image");  
68.     imshow("image", image);  
69.     namedWindow("result");  
70.     imshow("result", result);  
71.     namedWindow("image2", 0);  
72.     imshow("image2", image2);  
73.     waitKey(0);  
74.     return 0;  
75.   
76. }  

Explaination：

<1>代码中我们使用HSV颜色空间中的色调分量以描述所要搜索的物体。

<2>还需要注意的是如果一个颜色的饱和度偏低，会导致色调信息变得不稳定以及不可靠。这是因为低饱和度的颜色中，红绿蓝三个分量几乎是相等的。

<3>meanshift：intcvMeanShift( const CvArr* prob_image, CvRect window, CvTermCriteria criteria,CvConnectedComp* comp );

参数

prob_image

目标直方图的反向投影(见 cvCalcBackProject).

window

初始搜索窗口

criteria

确定窗口搜索停止的准则

五：通过比较直方图检索相似图片

基本思想就是得到两幅图像的直方图，然后通过opencv提供的函数compareHist来得到它们的相似程度，返回的是个double值。

Code：

[cpp]  view plain  copy  print ?

1. int main()  
2. {  
3.     Mat imageSource = imread("F:\\test\\tongtong.jpg", 1);  
4.     colorReduce(imageSource, 32);  
5.     MatND sourceHist = getHistogram(imageSource);  
6.     stringstream ss;  
7.     string str;  
8.     string strBest = "";  
9.     double minDistance = 256*100*100*100;  
10.     for(int i = 1; i < 10; i++)  
11.     {  
12.         str = "F:\\test\\";  
13.         ss.clear();  
14.         ss << str;  
15.         ss << i;  
16.         ss << ".jpg";  
17.         ss >> str;  
18.         Mat imageDst = imread(str, 1);  
19.         colorReduce(imageDst, 1);  
20.         MatND dstHist = getHistogram(imageDst);  
21.         double distance = compareHist(sourceHist, dstHist, CV_COMP_INTERSECT);  
22.         if(distance < minDistance)  
23.         {  
24.             strBest = str;  
25.             minDistance = distance;  
26.         }  
27.     }  
28.     Mat best = imread(strBest,1);  
29.     namedWindow(strBest);  
30.     imshow(strBest,best);  
31.     waitKey(0);  
32.     return 0;  
33.   
34. }  

Explaination：

<1>代码中使用了降低颜色数。这是因为直方图的比较大多数都是基于逐个容器的，即比较直方图容器时并不考虑相邻容器的影像。因此，测量相似度之前减少颜色空间是很重要的。

<2>compareHist函数直接明了，只需提供两个直方图，函数边返回测量距离。通过一个标志参数可以指定测量方法。代码中的参数CV_COMP_INTERSECT表示交叉测量法。即简单比较每个直方图容器的值，并保留最小的一个。相似性测量值只是这些最小值的和。

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