vpTemplateTrackerWarpRT.cpp

/****************************************************************************
 *
 * ViSP, open source Visual Servoing Platform software.
 * Copyright (C) 2005 - 2019 by Inria. All rights reserved.
 *
 * This software is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 * See the file LICENSE.txt at the root directory of this source
 * distribution for additional information about the GNU GPL.
 *
 * For using ViSP with software that can not be combined with the GNU
 * GPL, please contact Inria about acquiring a ViSP Professional
 * Edition License.
 *
 * See http://visp.inria.fr for more information.
 *
 * This software was developed at:
 * Inria Rennes - Bretagne Atlantique
 * Campus Universitaire de Beaulieu
 * 35042 Rennes Cedex
 * France
 *
 * If you have questions regarding the use of this file, please contact
 * Inria at visp@inria.fr
 *
 * This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
 * WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
 *
 * Description:
 * Template tracker.
 *
 * Authors:
 * Amaury Dame
 * Aurelien Yol
 * Fabien Spindler
 *
 *****************************************************************************/
#include <visp3/tt/vpTemplateTrackerWarpRT.h>

vpTemplateTrackerWarpRT::vpTemplateTrackerWarpRT()
{
  nbParam = 3;
}

void vpTemplateTrackerWarpRT::getParamPyramidDown(const vpColVector &p, vpColVector &p_down)
{
  p_down[0] = p[0];
  p_down[1] = p[1] / 2.;
  p_down[2] = p[2] / 2.;
}

void vpTemplateTrackerWarpRT::getParamPyramidUp(const vpColVector &p, vpColVector &p_up)
{
  p_up[0] = p[0];
  p_up[1] = p[1] * 2.;
  p_up[2] = p[2] * 2.;
}

void vpTemplateTrackerWarpRT::getdW0(const int &v, const int &u, const double &dv, const double &du, double *dIdW)
{
  dIdW[0] = -v * du + u * dv;
  dIdW[1] = du;
  dIdW[2] = dv;
}

void vpTemplateTrackerWarpRT::getdWdp0(const int &v, const int &u, double *dIdW)
{
  dIdW[0] = -v;
  dIdW[1] = 1.;
  dIdW[2] = 0;

  dIdW[3] = u;
  dIdW[4] = 0;
  dIdW[5] = 1.;
}

void vpTemplateTrackerWarpRT::warpX(const int &v1, const int &u1, double &v2, double &u2, const vpColVector &p)
{
  double c = cos(p[0]);
  double s = sin(p[0]);

  u2 = (c * u1) - (s * v1) + p[1];
  v2 = (s * u1) + (c * v1) + p[2];
}

void vpTemplateTrackerWarpRT::warpX(const vpColVector &X1, vpColVector &X2, const vpColVector &p)
{
  double c = cos(p[0]);
  double s = sin(p[0]);

  X2[0] = (c * X1[0]) - (s * X1[1]) + p[1];
  X2[1] = (s * X1[0]) + (c * X1[1]) + p[2];
}

void vpTemplateTrackerWarpRT::dWarp(const vpColVector &X, const vpColVector &, const vpColVector &p,
                                    vpMatrix &dM)
{
  double u = X[0];
  double v = X[1];
  double c = cos(p[0]);
  double s = sin(p[0]);

  dM[0][0] = - s * u - c * v;
  dM[0][1] = 1;
  dM[0][2] = 0;

  dM[1][0] = c * u - s * v;
  dM[1][1] = 0;
  dM[1][2] = 1;
}

void vpTemplateTrackerWarpRT::dWarpCompo(const vpColVector &, const vpColVector &,
                                         const vpColVector &p, const double *dwdp0, vpMatrix &dM)
{
  double c = cos(p[0]);
  double s = sin(p[0]);

  for (unsigned int i = 0; i < nbParam; i++) {
    dM[0][i] = (c * dwdp0[i]) - (s * dwdp0[i + nbParam]);
    dM[1][i] = (s * dwdp0[i]) + (c * dwdp0[i + nbParam]);
  }
}

void vpTemplateTrackerWarpRT::warpXInv(const vpColVector &X1, vpColVector &X2, const vpColVector &p)
{
  double c = cos(p[0]);
  double s = sin(p[0]);

  X2[0] = (c * X1[0]) - (s * X1[1]) + p[1];
  X2[1] = (s * X1[0]) + (c * X1[1]) + p[2];
}

void vpTemplateTrackerWarpRT::getParamInverse(const vpColVector &p, vpColVector &p_inv) const
{
  double c = cos(p[0]);
  double s = sin(p[0]);
  double u = p[1];
  double v = p[2];

  p_inv[0] = atan2(-s, c);
  p_inv[1] = -(c * u + s * v);
  p_inv[2] =  s * u - c * v;
}

void vpTemplateTrackerWarpRT::pRondp(const vpColVector &p1, const vpColVector &p2, vpColVector &p12) const
{
  double c1 = cos(p1[0]);
  double s1 = sin(p1[0]);
  double c2 = cos(p2[0]);
  double s2 = sin(p2[0]);
  double u1 = p1[1];
  double v1 = p1[2];
  double u2 = p2[1];
  double v2 = p2[2];

  p12[0] = atan2(s1 * c2 + c1 * s2, c1 * c2 - s1 * s2);
  p12[1] = c1 * u2 - s1 * v2 + u1;
  p12[2] = s1 * u2 + c1 * v2 + v1;
}