從一個SO答案中,使用光流分割需要密集光流。我修改了this code。它使用calcOpticalFlowFarneback函數計算兩個連續幀之間的密集流場。OpenCV C++基於光流的分割代碼拋出異常
SQRT來計算位移的大小(displacement_in_x^2 + displacement_in_y^2)
該函數返回位移CV_32FC2
這是我的代碼: -
#include "stdafx.h"
#include "cv.h"
#include "highgui.h"
#include <iostream>
using namespace cv;
using namespace std;
void drawOptFlowMap(const Mat& flow, Mat& cflowmap, int step,
double scale, const Scalar& color)
{
for(int y = 0; y < cflowmap.rows; y += step)
for(int x = 0; x < cflowmap.cols; x += step)
{
const Point2f& fxy = flow.at<Point2f>(y, x);
line(cflowmap, Point(x,y), Point(cvRound(x+fxy.x), cvRound(y+fxy.y)),
color);
circle(cflowmap, Point(x,y), 2, color, -1);
}
}
int main()
{
VideoCapture cap("vip.avi");
if(!cap.isOpened())
return -1;
Mat prevgray, gray, flow, cflow, frame;
double col = cap.get(CV_CAP_PROP_FRAME_WIDTH);
double row = cap.get(CV_CAP_PROP_FRAME_HEIGHT);
Mat flow_img = Mat::zeros(row,col,CV_8UC1);
namedWindow("flow", 1);
namedWindow("segmented",1);
for(;;)
{
cap >> frame;
cvtColor(frame, gray, CV_BGR2GRAY);
if(prevgray.data)
{
calcOpticalFlowFarneback(prevgray, gray, flow, 0.5, 3, 15, 15, 5, 1.2, 0);
cvtColor(prevgray, cflow, CV_GRAY2BGR);
drawOptFlowMap(flow, cflow, 8, 1.5, CV_RGB(0, 255, 0));
imshow("flow", cflow);
}
//This part I added to compute Magnitude of flow
for(int i=0;i<gray.rows;i++)
for(int j=0;j<gray.cols;j++)
flow_img.at<uchar>(i,j)= cvRound( sqrt (pow (flow.at<Vec2f>(i,j)[0] ,2)+ pow (flow.at<Vec2f>(i,j)[1],2)));
//threshold image
threshold(flow_img,flow_img,10,255,THRESH_BINARY); //some arbitrary thresh value 10
imshow("segmented",flow_img);
if(waitKey(30)>=0)
break;
std::swap(prevgray, gray);
}
return 0;
}
我正在使用Microsoft Visual C++ 2010,並且沒有生成錯誤。但是,該方案不會運行,並拋出了以下內容:
在控制檯窗口中,我得到:
我不明白,我所犯的錯誤。任何人都可以請幫助我。有類似的SO問題here,答案是我使用的代碼。儘管如此,我得到上述
哪一行導致異常? Visual Studio應該通過標記來顯示它,如果沒有通過調試代碼檢查它。 – cyriel
它沒有標記任何東西...一個窗口彈出顯示圖像(更新在我的qn)...當我調試它並沒有顯示任何問題...代碼運行顯示流向量的輸出幾秒鐘,然後退出顯示上面.. –
上傳某處thie AVI文件 - 我會檢查,所以我可以自己檢查。 – cyriel