• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
1 /*M///////////////////////////////////////////////////////////////////////////////////////
2 //
3 //  IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING.
4 //
5 //  By downloading, copying, installing or using the software you agree to this license.
6 //  If you do not agree to this license, do not download, install,
7 //  copy or use the software.
8 //
9 //
10 //                           License Agreement
11 //                For Open Source Computer Vision Library
12 //
13 // Copyright (C) 2000-2008, Intel Corporation, all rights reserved.
14 // Copyright (C) 2009-2011, Willow Garage Inc., all rights reserved.
15 // Third party copyrights are property of their respective owners.
16 //
17 // Redistribution and use in source and binary forms, with or without modification,
18 // are permitted provided that the following conditions are met:
19 //
20 //   * Redistribution's of source code must retain the above copyright notice,
21 //     this list of conditions and the following disclaimer.
22 //
23 //   * Redistribution's in binary form must reproduce the above copyright notice,
24 //     this list of conditions and the following disclaimer in the documentation
25 //     and/or other materials provided with the distribution.
26 //
27 //   * The name of the copyright holders may not be used to endorse or promote products
28 //     derived from this software without specific prior written permission.
29 //
30 // This software is provided by the copyright holders and contributors "as is" and
31 // any express or implied warranties, including, but not limited to, the implied
32 // warranties of merchantability and fitness for a particular purpose are disclaimed.
33 // In no event shall the Intel Corporation or contributors be liable for any direct,
34 // indirect, incidental, special, exemplary, or consequential damages
35 // (including, but not limited to, procurement of substitute goods or services;
36 // loss of use, data, or profits; or business interruption) however caused
37 // and on any theory of liability, whether in contract, strict liability,
38 // or tort (including negligence or otherwise) arising in any way out of
39 // the use of this software, even if advised of the possibility of such damage.
40 //
41 //M*/
42 
43 #ifndef __OPENCV_VIDEOSTAB_GLOBAL_MOTION_HPP__
44 #define __OPENCV_VIDEOSTAB_GLOBAL_MOTION_HPP__
45 
46 #include <vector>
47 #include <fstream>
48 #include "opencv2/core.hpp"
49 #include "opencv2/features2d.hpp"
50 #include "opencv2/opencv_modules.hpp"
51 #include "opencv2/videostab/optical_flow.hpp"
52 #include "opencv2/videostab/motion_core.hpp"
53 #include "opencv2/videostab/outlier_rejection.hpp"
54 
55 #ifdef HAVE_OPENCV_CUDAIMGPROC
56 #  include "opencv2/cudaimgproc.hpp"
57 #endif
58 
59 namespace cv
60 {
61 namespace videostab
62 {
63 
64 //! @addtogroup videostab_motion
65 //! @{
66 
67 /** @brief Estimates best global motion between two 2D point clouds in the least-squares sense.
68 
69 @note Works in-place and changes input point arrays.
70 
71 @param points0 Source set of 2D points (32F).
72 @param points1 Destination set of 2D points (32F).
73 @param model Motion model (up to MM_AFFINE).
74 @param rmse Final root-mean-square error.
75 @return 3x3 2D transformation matrix (32F).
76  */
77 CV_EXPORTS Mat estimateGlobalMotionLeastSquares(
78         InputOutputArray points0, InputOutputArray points1, int model = MM_AFFINE,
79         float *rmse = 0);
80 
81 /** @brief Estimates best global motion between two 2D point clouds robustly (using RANSAC method).
82 
83 @param points0 Source set of 2D points (32F).
84 @param points1 Destination set of 2D points (32F).
85 @param model Motion model. See cv::videostab::MotionModel.
86 @param params RANSAC method parameters. See videostab::RansacParams.
87 @param rmse Final root-mean-square error.
88 @param ninliers Final number of inliers.
89  */
90 CV_EXPORTS Mat estimateGlobalMotionRansac(
91         InputArray points0, InputArray points1, int model = MM_AFFINE,
92         const RansacParams &params = RansacParams::default2dMotion(MM_AFFINE),
93         float *rmse = 0, int *ninliers = 0);
94 
95 /** @brief Base class for all global motion estimation methods.
96  */
97 class CV_EXPORTS MotionEstimatorBase
98 {
99 public:
~MotionEstimatorBase()100     virtual ~MotionEstimatorBase() {}
101 
102     /** @brief Sets motion model.
103 
104     @param val Motion model. See cv::videostab::MotionModel.
105      */
setMotionModel(MotionModel val)106     virtual void setMotionModel(MotionModel val) { motionModel_ = val; }
107 
108     /**
109     @return Motion model. See cv::videostab::MotionModel.
110     */
motionModel() const111     virtual MotionModel motionModel() const { return motionModel_; }
112 
113     /** @brief Estimates global motion between two 2D point clouds.
114 
115     @param points0 Source set of 2D points (32F).
116     @param points1 Destination set of 2D points (32F).
117     @param ok Indicates whether motion was estimated successfully.
118     @return 3x3 2D transformation matrix (32F).
119      */
120     virtual Mat estimate(InputArray points0, InputArray points1, bool *ok = 0) = 0;
121 
122 protected:
MotionEstimatorBase(MotionModel model)123     MotionEstimatorBase(MotionModel model) { setMotionModel(model); }
124 
125 private:
126     MotionModel motionModel_;
127 };
128 
129 /** @brief Describes a robust RANSAC-based global 2D motion estimation method which minimizes L2 error.
130  */
131 class CV_EXPORTS MotionEstimatorRansacL2 : public MotionEstimatorBase
132 {
133 public:
134     MotionEstimatorRansacL2(MotionModel model = MM_AFFINE);
135 
setRansacParams(const RansacParams & val)136     void setRansacParams(const RansacParams &val) { ransacParams_ = val; }
ransacParams() const137     RansacParams ransacParams() const { return ransacParams_; }
138 
setMinInlierRatio(float val)139     void setMinInlierRatio(float val) { minInlierRatio_ = val; }
minInlierRatio() const140     float minInlierRatio() const { return minInlierRatio_; }
141 
142     virtual Mat estimate(InputArray points0, InputArray points1, bool *ok = 0);
143 
144 private:
145     RansacParams ransacParams_;
146     float minInlierRatio_;
147 };
148 
149 /** @brief Describes a global 2D motion estimation method which minimizes L1 error.
150 
151 @note To be able to use this method you must build OpenCV with CLP library support. :
152  */
153 class CV_EXPORTS MotionEstimatorL1 : public MotionEstimatorBase
154 {
155 public:
156     MotionEstimatorL1(MotionModel model = MM_AFFINE);
157 
158     virtual Mat estimate(InputArray points0, InputArray points1, bool *ok = 0);
159 
160 private:
161     std::vector<double> obj_, collb_, colub_;
162     std::vector<double> elems_, rowlb_, rowub_;
163     std::vector<int> rows_, cols_;
164 
set(int row,int col,double coef)165     void set(int row, int col, double coef)
166     {
167         rows_.push_back(row);
168         cols_.push_back(col);
169         elems_.push_back(coef);
170     }
171 };
172 
173 /** @brief Base class for global 2D motion estimation methods which take frames as input.
174  */
175 class CV_EXPORTS ImageMotionEstimatorBase
176 {
177 public:
~ImageMotionEstimatorBase()178     virtual ~ImageMotionEstimatorBase() {}
179 
setMotionModel(MotionModel val)180     virtual void setMotionModel(MotionModel val) { motionModel_ = val; }
motionModel() const181     virtual MotionModel motionModel() const { return motionModel_; }
182 
183     virtual Mat estimate(const Mat &frame0, const Mat &frame1, bool *ok = 0) = 0;
184 
185 protected:
ImageMotionEstimatorBase(MotionModel model)186     ImageMotionEstimatorBase(MotionModel model) { setMotionModel(model); }
187 
188 private:
189     MotionModel motionModel_;
190 };
191 
192 class CV_EXPORTS FromFileMotionReader : public ImageMotionEstimatorBase
193 {
194 public:
195     FromFileMotionReader(const String &path);
196 
197     virtual Mat estimate(const Mat &frame0, const Mat &frame1, bool *ok = 0);
198 
199 private:
200     std::ifstream file_;
201 };
202 
203 class CV_EXPORTS ToFileMotionWriter : public ImageMotionEstimatorBase
204 {
205 public:
206     ToFileMotionWriter(const String &path, Ptr<ImageMotionEstimatorBase> estimator);
207 
setMotionModel(MotionModel val)208     virtual void setMotionModel(MotionModel val) { motionEstimator_->setMotionModel(val); }
motionModel() const209     virtual MotionModel motionModel() const { return motionEstimator_->motionModel(); }
210 
211     virtual Mat estimate(const Mat &frame0, const Mat &frame1, bool *ok = 0);
212 
213 private:
214     std::ofstream file_;
215     Ptr<ImageMotionEstimatorBase> motionEstimator_;
216 };
217 
218 /** @brief Describes a global 2D motion estimation method which uses keypoints detection and optical flow for
219 matching.
220  */
221 class CV_EXPORTS KeypointBasedMotionEstimator : public ImageMotionEstimatorBase
222 {
223 public:
224     KeypointBasedMotionEstimator(Ptr<MotionEstimatorBase> estimator);
225 
setMotionModel(MotionModel val)226     virtual void setMotionModel(MotionModel val) { motionEstimator_->setMotionModel(val); }
motionModel() const227     virtual MotionModel motionModel() const { return motionEstimator_->motionModel(); }
228 
setDetector(Ptr<FeatureDetector> val)229     void setDetector(Ptr<FeatureDetector> val) { detector_ = val; }
detector() const230     Ptr<FeatureDetector> detector() const { return detector_; }
231 
setOpticalFlowEstimator(Ptr<ISparseOptFlowEstimator> val)232     void setOpticalFlowEstimator(Ptr<ISparseOptFlowEstimator> val) { optFlowEstimator_ = val; }
opticalFlowEstimator() const233     Ptr<ISparseOptFlowEstimator> opticalFlowEstimator() const { return optFlowEstimator_; }
234 
setOutlierRejector(Ptr<IOutlierRejector> val)235     void setOutlierRejector(Ptr<IOutlierRejector> val) { outlierRejector_ = val; }
outlierRejector() const236     Ptr<IOutlierRejector> outlierRejector() const { return outlierRejector_; }
237 
238     virtual Mat estimate(const Mat &frame0, const Mat &frame1, bool *ok = 0);
239 
240 private:
241     Ptr<MotionEstimatorBase> motionEstimator_;
242     Ptr<FeatureDetector> detector_;
243     Ptr<ISparseOptFlowEstimator> optFlowEstimator_;
244     Ptr<IOutlierRejector> outlierRejector_;
245 
246     std::vector<uchar> status_;
247     std::vector<KeyPoint> keypointsPrev_;
248     std::vector<Point2f> pointsPrev_, points_;
249     std::vector<Point2f> pointsPrevGood_, pointsGood_;
250 };
251 
252 #if defined(HAVE_OPENCV_CUDAIMGPROC) && defined(HAVE_OPENCV_CUDAOPTFLOW)
253 
254 class CV_EXPORTS KeypointBasedMotionEstimatorGpu : public ImageMotionEstimatorBase
255 {
256 public:
257     KeypointBasedMotionEstimatorGpu(Ptr<MotionEstimatorBase> estimator);
258 
setMotionModel(MotionModel val)259     virtual void setMotionModel(MotionModel val) { motionEstimator_->setMotionModel(val); }
motionModel() const260     virtual MotionModel motionModel() const { return motionEstimator_->motionModel(); }
261 
setOutlierRejector(Ptr<IOutlierRejector> val)262     void setOutlierRejector(Ptr<IOutlierRejector> val) { outlierRejector_ = val; }
outlierRejector() const263     Ptr<IOutlierRejector> outlierRejector() const { return outlierRejector_; }
264 
265     virtual Mat estimate(const Mat &frame0, const Mat &frame1, bool *ok = 0);
266     Mat estimate(const cuda::GpuMat &frame0, const cuda::GpuMat &frame1, bool *ok = 0);
267 
268 private:
269     Ptr<MotionEstimatorBase> motionEstimator_;
270     Ptr<cuda::CornersDetector> detector_;
271     SparsePyrLkOptFlowEstimatorGpu optFlowEstimator_;
272     Ptr<IOutlierRejector> outlierRejector_;
273 
274     cuda::GpuMat frame0_, grayFrame0_, frame1_;
275     cuda::GpuMat pointsPrev_, points_;
276     cuda::GpuMat status_;
277 
278     Mat hostPointsPrev_, hostPoints_;
279     std::vector<Point2f> hostPointsPrevTmp_, hostPointsTmp_;
280     std::vector<uchar> rejectionStatus_;
281 };
282 
283 #endif // defined(HAVE_OPENCV_CUDAIMGPROC) && defined(HAVE_OPENCV_CUDAOPTFLOW)
284 
285 /** @brief Computes motion between two frames assuming that all the intermediate motions are known.
286 
287 @param from Source frame index.
288 @param to Destination frame index.
289 @param motions Pair-wise motions. motions[i] denotes motion from the frame i to the frame i+1
290 @return Motion from the frame from to the frame to.
291  */
292 CV_EXPORTS Mat getMotion(int from, int to, const std::vector<Mat> &motions);
293 
294 //! @}
295 
296 } // namespace videostab
297 } // namespace cv
298 
299 #endif
300