c++ - 基于 OpenCV C++ 光流的分割代码抛出异常

Question

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

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

计算位移的大小

sqrt(displacement_in_x^2 + displacement_in_y^2)

the function returns displacement as 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答案是我使用的代码。我还是遇到了上面提到的问题

你可以找到我输入的视频 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 函数。在 Debug模式下逐步执行您的程序，以便我们可以知道是哪一个。