From f3bc1aede12fb65a555497d5d8e544768b63772c Mon Sep 17 00:00:00 2001
From: Alexey Spizhevoy <alexey.spizhevoy@itseez.com>
Date: Fri, 21 Sep 2012 10:46:00 +0400
Subject: [PATCH] removed obsolte global motion estimation method

---
 .../opencv2/videostab/global_motion.hpp       |  3 -
 modules/videostab/src/global_motion.cpp       | 66 -------------------
 2 files changed, 69 deletions(-)

diff --git a/modules/videostab/include/opencv2/videostab/global_motion.hpp b/modules/videostab/include/opencv2/videostab/global_motion.hpp
index 9616eae5c..c159c13fb 100644
--- a/modules/videostab/include/opencv2/videostab/global_motion.hpp
+++ b/modules/videostab/include/opencv2/videostab/global_motion.hpp
@@ -66,9 +66,6 @@ CV_EXPORTS Mat estimateGlobalMotionLeastSquares(
         InputOutputArray points0, InputOutputArray points1, int model = MM_AFFINE,
         float *rmse = 0);
 
-CV_EXPORTS Mat estimateGlobalMotionMedian(
-        InputArray points0, InputArray points1, int model, int size, int niters);
-
 CV_EXPORTS Mat estimateGlobalMotionRansac(
         InputArray points0, InputArray points1, int model = MM_AFFINE,
         const RansacParams &params = RansacParams::default2dMotion(MM_AFFINE),
diff --git a/modules/videostab/src/global_motion.cpp b/modules/videostab/src/global_motion.cpp
index 1eb2f0f46..d7ae1a55e 100644
--- a/modules/videostab/src/global_motion.cpp
+++ b/modules/videostab/src/global_motion.cpp
@@ -305,72 +305,6 @@ Mat estimateGlobalMotionLeastSquares(
 }
 
 
-Mat estimateGlobalMotionMedian(
-        InputArray points0, InputArray points1, int model, int size, int niters)
-{
-    // perform 'niters' iterations over points subsets ('size' elements each) estimating
-    // motions, after that select median motion parameters from the distribution
-
-    CV_Assert(model <= MM_AFFINE);
-    CV_Assert(points0.type() == points1.type());
-    const int npoints = points0.getMat().checkVector(2);
-    CV_Assert(points1.getMat().checkVector(2) == npoints);
-
-    const Point2f *points0_ = points0.getMat().ptr<Point2f>();
-    const Point2f *points1_ = points1.getMat().ptr<Point2f>();
-
-    // all estimated motions
-    vector<float> Ms[3][3];
-    for (int i = 0; i < 3; ++i)
-        for (int j = 0; j < 3; ++j)
-            Ms[i][j].resize(niters);
-
-    // current hypothesis
-    vector<int> indices(size);
-    vector<Point2f> subset0(size);
-    vector<Point2f> subset1(size);
-
-    RNG rng(0);
-
-    for (int iter = 0; iter < niters; ++iter)
-    {
-        for (int i = 0; i < size; ++i)
-        {
-            bool ok = false;
-            while (!ok)
-            {
-                ok = true;
-                indices[i] = static_cast<unsigned>(rng) % npoints;
-                for (int j = 0; j < i; ++j)
-                    if (indices[i] == indices[j])
-                        { ok = false; break; }
-            }
-        }
-        for (int i = 0; i < size; ++i)
-        {
-            subset0[i] = points0_[indices[i]];
-            subset1[i] = points1_[indices[i]];
-        }
-
-        Mat_<float> M = estimateGlobalMotionLeastSquares(subset0, subset1, model, 0);
-        for (int i = 0; i < 3; ++i)
-            for (int j = 0; j < 3; ++j)
-                Ms[i][j][iter] = M(i, j);
-    }
-
-    Mat_<float> medianM(3, 3);
-
-    for (int i = 0; i < 3; ++i)
-        for (int j = 0; j < 3; ++j)
-        {
-            nth_element(Ms[i][j].begin(), Ms[i][j].begin() + niters/2, Ms[i][j].end());
-            medianM(i, j) = Ms[i][j][niters/2];
-        }
-
-    return medianM;
-}
-
-
 Mat estimateGlobalMotionRansac(
         InputArray points0, InputArray points1, int model, const RansacParams &params,
         float *rmse, int *ninliers)