optical_flow.cpp

/*M///////////////////////////////////////////////////////////////////////////////////////
//
//  IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING.
//
//  By downloading, copying, installing or using the software you agree to this license.
//  If you do not agree to this license, do not download, install,
//  copy or use the software.
//
//
//                           License Agreement
//                For Open Source Computer Vision Library
//
// Copyright (C) 2000-2008, Intel Corporation, all rights reserved.
// Copyright (C) 2009-2011, Willow Garage Inc., all rights reserved.
// Third party copyrights are property of their respective owners.
//
// Redistribution and use in source and binary forms, with or without modification,
// are permitted provided that the following conditions are met:
//
//   * Redistribution's of source code must retain the above copyright notice,
//     this list of conditions and the following disclaimer.
//
//   * Redistribution's in binary form must reproduce the above copyright notice,
//     this list of conditions and the following disclaimer in the documentation
//     and/or other materials provided with the distribution.
//
//   * The name of the copyright holders may not be used to endorse or promote products
//     derived from this software without specific prior written permission.
//
// This software is provided by the copyright holders and contributors "as is" and
// any express or implied warranties, including, but not limited to, the implied
// warranties of merchantability and fitness for a particular purpose are disclaimed.
// In no event shall the Intel Corporation or contributors be liable for any direct,
// indirect, incidental, special, exemplary, or consequential damages
// (including, but not limited to, procurement of substitute goods or services;
// loss of use, data, or profits; or business interruption) however caused
// and on any theory of liability, whether in contract, strict liability,
// or tort (including negligence or otherwise) arising in any way out of
// the use of this software, even if advised of the possibility of such damage.
//
//M*/

#include "precomp.hpp"
#include "opencv2/video.hpp"
#include "opencv2/videostab/optical_flow.hpp"
#include "opencv2/videostab/ring_buffer.hpp"

#ifdef HAVE_OPENCV_CUDAARITHM
  #include "opencv2/cudaarithm.hpp"
#endif

namespace cv
{
namespace videostab
{

void SparsePyrLkOptFlowEstimator::run(
        InputArray frame0, InputArray frame1, InputArray points0, InputOutputArray points1,
        OutputArray status, OutputArray errors)
{
    calcOpticalFlowPyrLK(frame0, frame1, points0, points1, status, errors, winSize_, maxLevel_);
}


#ifdef HAVE_OPENCV_CUDAOPTFLOW

SparsePyrLkOptFlowEstimatorGpu::SparsePyrLkOptFlowEstimatorGpu()
{
    CV_Assert(cuda::getCudaEnabledDeviceCount() > 0);
    optFlowEstimator_ = cuda::SparsePyrLKOpticalFlow::create();
}


void SparsePyrLkOptFlowEstimatorGpu::run(
        InputArray frame0, InputArray frame1, InputArray points0, InputOutputArray points1,
        OutputArray status, OutputArray errors)
{
    frame0_.upload(frame0.getMat());
    frame1_.upload(frame1.getMat());
    points0_.upload(points0.getMat());

    if (errors.needed())
    {
        run(frame0_, frame1_, points0_, points1_, status_, errors_);
        errors_.download(errors.getMatRef());
    }
    else
        run(frame0_, frame1_, points0_, points1_, status_);

    points1_.download(points1.getMatRef());
    status_.download(status.getMatRef());
}


void SparsePyrLkOptFlowEstimatorGpu::run(
        const cuda::GpuMat &frame0, const cuda::GpuMat &frame1, const cuda::GpuMat &points0,
        cuda::GpuMat &points1, cuda::GpuMat &status, cuda::GpuMat &errors)
{
    optFlowEstimator_->setWinSize(winSize_);
    optFlowEstimator_->setMaxLevel(maxLevel_);
    optFlowEstimator_->calc(frame0, frame1, points0, points1, status, errors);
}


void SparsePyrLkOptFlowEstimatorGpu::run(
        const cuda::GpuMat &frame0, const cuda::GpuMat &frame1, const cuda::GpuMat &points0,
        cuda::GpuMat &points1, cuda::GpuMat &status)
{
    optFlowEstimator_->setWinSize(winSize_);
    optFlowEstimator_->setMaxLevel(maxLevel_);
    optFlowEstimator_->calc(frame0, frame1, points0, points1, status);
}


DensePyrLkOptFlowEstimatorGpu::DensePyrLkOptFlowEstimatorGpu()
{
    CV_Assert(cuda::getCudaEnabledDeviceCount() > 0);
    optFlowEstimator_ = cuda::DensePyrLKOpticalFlow::create();
}


void DensePyrLkOptFlowEstimatorGpu::run(
        InputArray frame0, InputArray frame1, InputOutputArray flowX, InputOutputArray flowY,
        OutputArray errors)
{
    frame0_.upload(frame0.getMat());
    frame1_.upload(frame1.getMat());

    optFlowEstimator_->setWinSize(winSize_);
    optFlowEstimator_->setMaxLevel(maxLevel_);

    if (errors.needed())
    {
        CV_Error(Error::StsNotImplemented, "DensePyrLkOptFlowEstimatorGpu doesn't support errors calculation");
    }
    else
    {
        cuda::GpuMat flow;
        optFlowEstimator_->calc(frame0_, frame1_, flow);

        cuda::GpuMat flows[2];
        cuda::split(flow, flows);

        flowX_ = flows[0];
        flowY_ = flows[1];
    }

    flowX_.download(flowX.getMatRef());
    flowY_.download(flowY.getMatRef());
}

#endif // HAVE_OPENCV_CUDAOPTFLOW

} // namespace videostab
} // namespace cv