vpFeatureBuilderLine.cpp

/*
 * ViSP, open source Visual Servoing Platform software.
 * Copyright (C) 2005 - 2025 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 https://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.
 */

 
#include <visp3/core/vpDebug.h>
#include <visp3/core/vpMath.h>
#include <visp3/visual_features/vpFeatureBuilder.h>
 
BEGIN_VISP_NAMESPACE
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;
 
    // Move from (i,j) frame to (u,v) frame. rho unchanged but theta is changed
    // since  j <-> u, i <-> v so that sin(theta) = cos(thetap) and cos(theta) = sin(thetap)
    // which means theta = pi/2 - thetap
 
    // It would be more simple to output theta and rho in (u,v) frame
    // in vpMeLine.cpp by just inverting sin and cos in the atan2 of
    // vpMeLine::computeRhoTheta() to avoid the "hack" below.
    // Not done since I am not sure vpMeLine or this function are not used
    // for other purposes.
 
    thetap = M_PI / 2.0 - thetap;  // input in [-pi;pi] so output in [-pi/2;3pi/2]
    if (thetap > M_PI) {   // thetap in [-pi;pi]
      thetap -= 2.0 * M_PI;
    }
    vpPixelMeterConversion::convertLine(cam, rhop, thetap, rho, theta);
    s.buildFrom(rho, theta);
  }
  catch (...) {
    vpERROR_TRACE("Error caught");
    throw;
  }
}
#endif //#ifdef VISP_HAVE_MODULE_ME
END_VISP_NAMESPACE