3

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

根据此处和此处的 SO 答案，我尝试使用

sqrt(位移_in_x^2 + 位移_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 并且没有构建错误。但程序不会运行并抛出以下内容：

构建错误

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

c0nsole

我不明白我在哪里犯了错误。谁能帮帮我。这里有类似的 SO 问题，答案是我使用的代码。我仍然遇到上述问题

你可以在这里或这里找到我的输入视频

4

2 回答 2

0

你确定线路`

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

被允许？Point2f不是cv::Mat类文档中描述的类型之一。尝试使用cv::Vec2fas 模板参数。然后，您可以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 函数。在调试模式下，在您的程序中一步一步地进行，以便我们知道是哪一个。

于 2013-04-15T06:44:47.427 回答
0

代码行

std::swap(prevgray, gray);

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

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

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

gray.copyTo(prevgray);

为prevgray分配自己的内存。否则prevgray的数据指针指向gray的内容。

于 2013-08-06T17:06:55.497 回答