我正在尝试为来自单通道黑白视频流的每个帧计算直方图。最后,我将把这些直方图的投影反馈给OpenCV的camshift函数。
这是我的代码:
#include "opencv2/highgui/highgui.hpp"
#include "opencv2/imgproc/imgproc.hpp"
#include "opencv2/opencv.hpp"
#include "opencv2/core/core.hpp"
#include <iostream>
#include <vector>
using namespace std;
using namespace cv;
int main(int argc, char* argv[]) {
Mat frame, back, fore, hist, image, v_hist;
BackgroundSubtractorMOG2 bgs;
vector<vector<Point> > contours;
vector<Vec4i> hierarchy;
int bins = 8;
int histSize[] = { bins };
float range[] = { 0, 255 }; //values will range from 0 to 256
const float* histRange = { range };
bool uniform = true, accumulate = true;
VideoCapture cap(0);
if (!cap.isOpened()) {
cout << "Cannot open the video cam" << endl;
return -1;
}
namedWindow("histogram",CV_WINDOW_AUTOSIZE);
namedWindow("Foreground",CV_WINDOW_AUTOSIZE);
for (;;) {
if (!cap.read(frame)){
cout << "No frame" << endl;
break;
}
frame.copyTo(image);
//***start mask computation***
bgs.operator()(image, fore);
bgs.getBackgroundImage(back);
findContours(fore, contours, hierarchy, CV_RETR_EXTERNAL, CV_CHAIN_APPROX_SIMPLE);
int idx = 0;
for( ; idx >= 0; idx = hierarchy[idx][0] ){
Scalar color( 255, 255, 255 );
drawContours( fore, contours, idx, color, CV_FILLED, 8, hierarchy );
}
erode(fore, fore, Mat());
erode(fore, fore, Mat());
dilate(fore, fore, Mat());
dilate(fore, fore, Mat());
//****end mask computation****
calcHist( &fore, 1, 0, Mat(), v_hist, 1, histSize, &histRange, uniform, accumulate);
int hist_w = 512; int hist_h = 400;
int bin_w = cvRound( (double) hist_w/ (*histSize) );
Mat histImage( hist_h, hist_w, CV_8UC3, Scalar(0, 0, 0) );
normalize( v_hist, v_hist, 0, histImage.rows, NORM_MINMAX, -1, Mat() );
//draw histogram
for (int i = 1; i < *histSize; i++)
line ( histImage,
Point(bin_w*(i-1), hist_h - cvRound(v_hist.at<float>(i-1)) ),
Point( bin_w*(i), hist_h - cvRound(v_hist.at<float>(i)) ),
Scalar( 255, 255, 255), 2, 8, 0 );
if(!histImage.empty()) imshow("histogram", histImage);
if(!fore.empty()) imshow("Foreground", fore);
if (waitKey(30) >= 0) break;
}
return 0;
}
这是输出最初的样子(记住 - 一开始,前景蒙版全黑):
..但即使在前景提取后,直方图仍保持不变:
看起来这个直方图仅适用于第一个全黑图像,这就是为什么线从x =(x轴的长度)/ 8开始到y = 0的原因。这是有道理的,因为直方图有8个区间,并且这个原始图像中的所有值都是0.(FWIW:当只有2个区间时,线条从y = 0开始到y = 0(x轴的长度) / 2)
我过去曾使用OpenCV绘制线条和矩形,并且从来没有像我在这里做的那样做太多不同,所以我有理由相信这不是任何问题。绘图功能。如果是这种情况,似乎直方图没有随着新视频数据的出现而更新。
有没有人知道为什么会发生这种情况?
答案 0 :(得分:0)
虽然我还没有确切知道为什么绘制的直方图没有更新,但我已经做了一些测试,现在确信直方图在不同的帧之间发生了变化。
这是我在之前的程序中添加的代码:
//new declarations
Mat last_hist;
bool first = true;
//code is the same, leading up to...
normalize( v_hist, v_hist, 0, histImage.rows, NORM_MINMAX, -1, Mat() );
if (!first){
eq = std::equal( v_hist.begin<uchar>(), v_hist.end<uchar>(), last_hist.begin<uchar>() );
if (eq)
cout << "hists equal" << endl;
else
cout << "hists not equal" << endl;
}
else
first = false;
//draw histogram
for (int i = 1; i < *histSize; i++)
//...
if(!frame.empty())imshow("Frame", frame);
v_hist.copyTo(last_hist);
if (waitKey(30) >= 0) break;
//from here to end is the same
在equal()
来自how to check whether two matrixes are identical in OpenCV之后,在条件语句中使用STL normalize()
。
当我运行此代码时,在前景蒙版被充分打扰后,我会看到“hists not equal”打印几帧。