vpFeatureMoment.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:
 * Base for all moment features
 */
 
#include <visp3/core/vpMath.h>
#include <visp3/core/vpMoment.h>
#include <visp3/core/vpMomentDatabase.h>
#include <visp3/visual_features/vpFeatureMoment.h>
#include <visp3/visual_features/vpFeatureMomentDatabase.h>
 
#include <visp3/core/vpException.h>
#include <visp3/visual_features/vpFeatureException.h>
 
#include <vector>
#include <visp3/core/vpDebug.h>
 
BEGIN_VISP_NAMESPACE
class vpBasicFeature;

void vpFeatureMoment::init()
{
  // feature dimension
  /*
   * The dimension of the visual feature is set according to the size of the
   * vpMoment associated to it. This partly explains why vpFeatureMomentBasic
   * cannot be used directly as a visual feature.
   */
  if (this->moment != nullptr)
    dim_s = static_cast<unsigned int>(this->moment->get().size());
  else
    dim_s = 0;
 
  nbParameters = 1;
 
  // memory allocation
  s.resize(dim_s);
  for (unsigned int i = 0; i < dim_s; i++)
    s[i] = 0;
 
  if (flags == nullptr)
    flags = new bool[nbParameters];
  for (unsigned int i = 0; i < nbParameters; i++)
    flags[i] = false;
}

int vpFeatureMoment::getDimension(unsigned int select) const
{
  int dim = 0;
 
  for (unsigned int i = 0; i < dim_s; ++i)
    if (vpBasicFeature::FEATURE_LINE[i] & select)
      dim++;
 
  return dim;
}

void vpFeatureMoment::print(unsigned int select) const
{
  for (unsigned int i = 0; i < dim_s; ++i) {
    if (vpBasicFeature::FEATURE_LINE[i] & select) {
      std::cout << s[i] << ",";
    }
  }
 
  std::cout << std::endl;
}

void vpFeatureMoment::display(const vpCameraParameters &cam, const vpImage<unsigned char> &I, const vpColor &color,
                              unsigned int thickness) const
{
  // visual representation of a moment doesn't often make sense
  (void)cam;
  (void)I;
  (void)color;
  (void)thickness;
}

void vpFeatureMoment::display(const vpCameraParameters &cam, const vpImage<vpRGBa> &I, const vpColor &color,
                              unsigned int thickness) const
{
  (void)cam;
  (void)I;
  (void)color;
  (void)thickness;
}

void vpFeatureMoment::update(double A_, double B_, double C_)
{
  this->A = A_;
  this->B = B_;
  this->C = C_;
 
  if (moment == nullptr) {
    bool found;
    this->moment = &(moments.get(momentName(), found));
    if (!found)
      throw vpException(vpException::notInitialized, "Moment not found for feature");
  }
  nbParameters = 1;
  if (this->moment != nullptr) {
    dim_s = static_cast<unsigned int>(this->moment->get().size());
 
    s.resize(dim_s);
 
    for (unsigned int i = 0; i < dim_s; i++)
      s[i] = this->moment->get()[i];
 
    if (flags == nullptr)
      flags = new bool[nbParameters];
    for (unsigned int i = 0; i < nbParameters; i++)
      flags[i] = false;
  }
  else
    dim_s = 0;
 
  compute_interaction();
}

vpMatrix vpFeatureMoment::interaction(unsigned int select)
{
  vpMatrix L(0, 0);
 
  for (unsigned int i = 0; i < dim_s; ++i) {
    if (vpBasicFeature::FEATURE_LINE[i] & select) {
      L.stack(interaction_matrices[i]);
    }
  }
 
  return L;
}

vpBasicFeature *vpFeatureMoment::duplicate() const
{
  vpFeatureMoment *feat = new vpMomentGenericFeature(moments, A, B, C, featureMomentsDataBase, moment);
  feat->interaction_matrices = interaction_matrices;
  feat->dim_s = dim_s;
  feat->nbParameters = nbParameters;
  // memory allocation
  feat->s.resize(dim_s);
  for (unsigned int i = 0; i < dim_s; i++)
    feat->s[i] = this->s[i];
 
  unsigned int n = static_cast<unsigned int>(nbParameters);
  feat->flags = new bool[n];
  for (unsigned int i = 0; i < n; i++)
    feat->flags[i] = flags[i];
 
  return feat;
}

void vpFeatureMoment::linkTo(vpFeatureMomentDatabase &featureMoments)
{
  m_name = name();
  this->featureMomentsDataBase = &featureMoments;
 
  featureMoments.add(*this, m_name);
}
 
void vpFeatureMoment::compute_interaction() { }

void vpFeatureMoment::printDependencies(std::ostream &os) const
{
  os << " WARNING : Falling back to base class version of "
    "printDependencies() in vpFeatureMoment. To prevent that, this has "
    "to be implemented in the derived classes!"
    << std::endl;
}
 
VISP_EXPORT std::ostream &operator<<(std::ostream &os, const vpFeatureMoment &featM)
{
  /*
   * - A static_cast is forced here since interaction() defined in vpBasicFeature()
   *   is not const. But introducing const in vpBasicFeature() can break a lot of
   *   client code.
   * - 6 corresponds to 6 velocities in standard interaction matrix
   */
  vpMatrix Lcomplete(static_cast<unsigned int>(featM.getDimension()), 6);
  Lcomplete = const_cast<vpFeatureMoment &>(featM).interaction(vpBasicFeature::FEATURE_ALL);
  Lcomplete.matlabPrint(os);
  return os;
}
 
END_VISP_NAMESPACE