vpFeatureDisplay.cpp

/*
 * ViSP, open source Visual Servoing Platform software.
 * Copyright (C) 2005 - 2024 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:
 * Interface with the image for feature display.
 *
*****************************************************************************/
 
#include <visp3/core/vpFeatureDisplay.h>
 
// Meter/pixel conversion
#include <visp3/core/vpMeterPixelConversion.h>
 
// display
#include <visp3/core/vpDisplay.h>
 
// math
#include <visp3/core/vpMath.h>
 
#include <visp3/core/vpImagePoint.h>
 
BEGIN_VISP_NAMESPACE
void vpFeatureDisplay::displayPoint(double x, double y, const vpCameraParameters &cam, const vpImage<unsigned char> &I,
                                    const vpColor &color, unsigned int thickness)
{
  vpImagePoint ip; // pixel coordinates in float
  vpMeterPixelConversion::convertPoint(cam, x, y, ip);
  vpDisplay::displayCross(I, ip, 15, color, thickness);
}
 
void vpFeatureDisplay::displayLine(double rho, double theta, const vpCameraParameters &cam,
                                   const vpImage<unsigned char> &I, const vpColor &color, unsigned int thickness)
{
  //  --comment: x times cos(theta) plus y times sin(theta) minus rho equals 0
  double rhop, thetap;
  vpMeterPixelConversion::convertLine(cam, rho, theta, rhop, thetap);
 
  //  --comment: u times cos thetap plus v times sin thetap minus rhop equals 0
  double co = cos(thetap);
  double si = sin(thetap);
  double c = -rhop;
 
  double a = si;
  double b = co;
  vpImagePoint ip1, ip2;
 
  if (fabs(a) < fabs(b)) {
    ip1.set_ij(0, (-c) / b);
    double h = I.getHeight() - 1;
    ip2.set_ij(h, (-c - (a * h)) / b);
    vpDisplay::displayLine(I, ip1, ip2, color, thickness);
  }
  else {
    ip1.set_ij((-c) / a, 0);
    double w = I.getWidth() - 1;
    ip2.set_ij((-c - (b * w)) / a, w);
    vpDisplay::displayLine(I, ip1, ip2, color, thickness);
  }
}
 
void vpFeatureDisplay::displayCylinder(double rho1, double theta1, double rho2, double theta2,
                                       const vpCameraParameters &cam, const vpImage<unsigned char> &I,
                                       const vpColor &color, unsigned int thickness)
{
  displayLine(rho1, theta1, cam, I, color, thickness);
  displayLine(rho2, theta2, cam, I, color, thickness);
}
 
void vpFeatureDisplay::displayEllipse(double x, double y, double n20, double n11, double n02,
                                      const vpCameraParameters &cam, const vpImage<unsigned char> &I,
                                      const vpColor &color, unsigned int thickness)
{
  vpImagePoint center;
  double n20_p, n11_p, n02_p;
  const unsigned int index_0 = 0;
  const unsigned int index_1 = 1;
  const unsigned int index_2 = 2;
  const unsigned int index_3 = 3;
  const unsigned int index_4 = 4;
  vpCircle circle;
  circle.p[index_0] = x;
  circle.p[index_1] = y;
  circle.p[index_2] = n20;
  circle.p[index_3] = n11;
  circle.p[index_4] = n02;
 
  vpMeterPixelConversion::convertEllipse(cam, circle, center, n20_p, n11_p, n02_p);
  vpDisplay::displayEllipse(I, center, n20_p, n11_p, n02_p, true, color, thickness);
}
 
void vpFeatureDisplay::displayPoint(double x, double y, const vpCameraParameters &cam, const vpImage<vpRGBa> &I,
                                    const vpColor &color, unsigned int thickness)
{
  vpImagePoint ip; // pixel coordinates in float
  vpMeterPixelConversion::convertPoint(cam, x, y, ip);
 
  vpDisplay::displayCross(I, ip, 15, color, thickness);
}
 
void vpFeatureDisplay::displayLine(double rho, double theta, const vpCameraParameters &cam, const vpImage<vpRGBa> &I,
                                   const vpColor &color, unsigned int thickness)
{
  // --comment: x times cos of theta plus y times sin of theta minus rho equals 0
  double rhop, thetap;
  vpMeterPixelConversion::convertLine(cam, rho, theta, rhop, thetap);
 
  // --comment: u times cos of thetap plus v times sin of thetap minus rhop equals 0
  double co = cos(thetap);
  double si = sin(thetap);
  double c = -rhop;
 
  double a = si;
  double b = co;
  vpImagePoint ip1, ip2;
 
  if (fabs(a) < fabs(b)) {
    ip1.set_ij(0, (-c) / b);
    double h = I.getHeight() - 1;
    ip2.set_ij(h, (-c - (a * h)) / b);
    vpDisplay::displayLine(I, ip1, ip2, color, thickness);
  }
  else {
    ip1.set_ij((-c) / a, 0);
    double w = I.getWidth() - 1;
    ip2.set_ij((-c - (b * w)) / a, w);
    vpDisplay::displayLine(I, ip1, ip2, color, thickness);
  }
}
 
void vpFeatureDisplay::displayCylinder(double rho1, double theta1, double rho2, double theta2,
                                       const vpCameraParameters &cam, const vpImage<vpRGBa> &I, const vpColor &color,
                                       unsigned int thickness)
{
  displayLine(rho1, theta1, cam, I, color, thickness);
  displayLine(rho2, theta2, cam, I, color, thickness);
}
 
void vpFeatureDisplay::displayEllipse(double x, double y, double n20, double n11, double n02,
                                      const vpCameraParameters &cam, const vpImage<vpRGBa> &I, const vpColor &color,
                                      unsigned int thickness)
{
  vpImagePoint center;
  double n20_p, n11_p, n02_p;
  const unsigned int index_0 = 0;
  const unsigned int index_1 = 1;
  const unsigned int index_2 = 2;
  const unsigned int index_3 = 3;
  const unsigned int index_4 = 4;
  vpCircle circle;
  circle.p[index_0] = x;
  circle.p[index_1] = y;
  circle.p[index_2] = n20;
  circle.p[index_3] = n11;
  circle.p[index_4] = n02;
 
  vpMeterPixelConversion::convertEllipse(cam, circle, center, n20_p, n11_p, n02_p);
  vpDisplay::displayEllipse(I, center, n20_p, n11_p, n02_p, true, color, thickness);
}
END_VISP_NAMESPACE