vpFeatureBuilderLine.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:
 * Conversion between tracker and visual feature line.
 *
 * Authors:
 * Eric Marchand
 *
 *****************************************************************************/

#include <visp3/core/vpMath.h>
#include <visp3/visual_features/vpFeatureBuilder.h>

void vpFeatureBuilder::create(vpFeatureLine &s, const vpLine &t)
{
  try {
    double A, B, C, D;
    s.setRhoTheta(t.getRho(), t.getTheta());

    if (fabs(t.cP[3]) > fabs(t.cP[7])) // |D1| > |D2|
    {
      A = t.cP[0];
      B = t.cP[1];
      C = t.cP[2];
      D = t.cP[3];
    } else {
      A = t.cP[4];
      B = t.cP[5];
      C = t.cP[6];
      D = t.cP[7];
    }

    s.setABCD(A, B, C, D);

  } catch (...) {
    vpERROR_TRACE("Error caught");
    throw;
  }
}

void vpFeatureBuilder::create(vpFeatureLine &s, const vpCylinder &t, int line)
{
  try {
    double a = t.getA();
    double b = t.getB();
    double c = t.getC();

    double x0 = t.getX();
    double y0 = t.getY();
    double z0 = t.getZ();

    double R = t.getR();

    double D =
        vpMath::sqr(x0) + vpMath::sqr(y0) + vpMath::sqr(z0) - vpMath::sqr(R) - vpMath::sqr(a * x0 + b * y0 + c * z0);

    double alpha1 = (1 - a * a) * x0 - a * b * y0 - a * c * z0;
    double beta1 = -a * b * x0 + (1 - b * b) * y0 - b * c * z0;
    double gamma1 = -a * c * x0 - b * c * y0 + (1 - c * c) * z0;

    D *= -1;

    if (D < 0) {
      alpha1 *= -1;
      beta1 *= -1;
      gamma1 *= -1;
      D *= -1;
    }

    s.setABCD(alpha1, beta1, gamma1, D);

    if (line == vpCylinder::line1) {

      s.setRhoTheta(t.getRho1(), t.getTheta1());

    } else {

      s.setRhoTheta(t.getRho2(), t.getTheta2());
    }
  } catch (...) {
    vpERROR_TRACE("Error caught");
    throw;
  }
}

#ifdef VISP_HAVE_MODULE_ME

void vpFeatureBuilder::create(vpFeatureLine &s, const vpCameraParameters &cam, const vpMeLine &t)
{
  try {
    double rhop = t.getRho();
    double thetap = t.getTheta();
    double rho;
    double theta;

    // Gives the rho and theta coordinates in the (u,v) coordinate system.
    if (thetap >= 0 && thetap < M_PI / 2) {
      thetap = M_PI / 2 - thetap;
    }

    else if (thetap >= M_PI / 2 && thetap < 3 * M_PI / 2) {
      thetap = 3 * M_PI / 2 + M_PI - thetap;
    }

    else if (thetap >= 3 * M_PI / 2 && thetap <= 2 * M_PI) {
      thetap = M_PI / 2 + 2 * M_PI - thetap;
    }

    // while (thetap > M_PI/2)  { thetap -= M_PI ; rhop *= -1 ; }
    // while (thetap < -M_PI/2) { thetap += M_PI ; rhop *= -1 ; }

    //  vpTRACE("pixel %f %f",rhop, thetap) ;
    vpPixelMeterConversion::convertLine(cam, rhop, thetap, rho, theta);

    while (theta > M_PI) {
      theta -= 2 * M_PI;
    }
    while (theta < -M_PI) {
      theta += 2 * M_PI;
    }
    //   vpTRACE("meter %f %f",rho, theta) ;
    /*

    while(theta < -M_PI)    theta += 2*M_PI ;
    while(theta >= M_PI)    theta -= 2*M_PI ;

    // If theta is between -90 and -180 get the equivalent
    // between 0 and 90
    if(theta <-M_PI/2)
    {
      theta += M_PI ;
      rho *= -1 ;
    }
    // If theta is between 90 and 180 get the equivalent
    // between 0 and -90
    if(theta >M_PI/2)
    {
      theta -= M_PI ;
      rho *= -1 ;
    }
    */
    s.buildFrom(rho, theta);
    //   vpTRACE("meter %f %f",rho, theta) ;

  } catch (...) {
    vpERROR_TRACE("Error caught");
    throw;
  }
}
#endif //#ifdef VISP_HAVE_MODULE_ME