Question

从一个SO答案中，需要使用光流密集光流进行分割。我修改了this code 。它使用calcOpticalFlowFarneback函数来计算两个连续帧之间的密集流场。

根据SO答案here和here，我尝试使用

计算位移的大小

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，没有构建错误。但该计划不会运行并抛出以下内容：

build error

在控制台窗口中，我得到：

c0nsole

我不明白我犯了什么错误。任何人都可以帮助我。有类似的SO问题here，答案是我使用的代码。我仍然遇到上述问题

您可以找到我的输入视频 here或here

Answer 1

你确定这行`

const Point2f& fxy = flow.at<Point2f>(y, x);

是允许的吗？ Point2f不是cv::Mat类文档中描述的类型之一。请尝试使用cv::Vec2f作为模板参数。然后，您可以使用结果填充Point2f变量。

下一个问题：在计算流量的8位幅度的行中，使用saturate_cast运算符会更安全：

flow_img.at<uchar>(i,j)= cv::saturate_cast<uchar>(cvRound(  sqrt ( pow (flow.at<Vec2f>(i,j)[0] ,2 )+ pow ( flow.at<Vec2f>(i,j)[1],2) )));

不确定是否会导致您的程序崩溃。

错误告诉您的是您正在使用错误的参数调用OpenCV函数。在调试模式下，在程序中逐步进行，以便我们知道哪一个。

Answer 2

代码行

std::swap(prevgray, gray);

可能会导致问题，因为std :: swap的行为是

c(a); a=b; b=c;

以这种方式使用=运算符，它不应用输入图像的深层副本。更好地使用

gray.copyTo(prevgray);

为prevgray分配自己的内存。另外，prevgray的数据指针指向灰色内容。

OpenCV C ++基于光流的分段代码抛出异常

2 个答案: