ParametricShapeArcLengthFunctor.h

 
#pragma once

 
#if defined(ParametricShapeArcLengthFunctor_RECURSES)
#error Recursive header files inclusion detected in ParametricShapeArcLengthFunctor.h
#else // defined(ParametricShapeArcLengthFunctor_RECURSES)
#define ParametricShapeArcLengthFunctor_RECURSES
 
#if !defined ParametricShapeArcLengthFunctor_h
#define ParametricShapeArcLengthFunctor_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
 
namespace DGtal
{

  // template class ParametricShapeArcLengthFunctor
  template <typename TParametricShape>
  class ParametricShapeArcLengthFunctor
  {
 
    // ----------------------- Standard services ------------------------------
  public:

    typedef TParametricShape ParametricShape;

    typedef typename TParametricShape::RealPoint RealPoint;

    typedef double Quantity;

    ParametricShapeArcLengthFunctor() = delete;
 

    ParametricShapeArcLengthFunctor(const ParametricShape &aShape): myShape(aShape) {};
 

    ~ParametricShapeArcLengthFunctor() = default;
 
 
    // ----------------------- Interface --------------------------------------
  public:
    ParametricShapeArcLengthFunctor & operator= ( const ParametricShapeArcLengthFunctor & other ) = delete;

    Quantity operator()(const RealPoint &aFirstPoint,const RealPoint &aSecondPoint) const
    {
      //determining nbSamples from the bounding box size of the shape
      RealPoint v = myShape.getUpperBound() - myShape.getLowerBound();
      double n = v.norm(RealPoint::L_infty);
      unsigned int nbSamples = (unsigned int) ceil( n*100 );
 
      //computes the angles
      double t = myShape.parameter( aFirstPoint );
      double t2 = myShape.parameter( aSecondPoint );
      return myShape.arclength( t, t2, nbSamples );
 
    }

    Quantity operator()() const
    {
      //determining nbSamples from the bounding box size of the shape
      RealPoint v = myShape.getUpperBound() - myShape.getLowerBound();
      double n = v.norm(RealPoint::L_infty);
      unsigned int nbSamples = (unsigned int) ceil( n*100 );
 
      return myShape.arclength( 0,2*M_PI, nbSamples );
    }
 
    // ------------------------- Private Data --------------------------------
  private:

    const ParametricShape& myShape;
 
    // ------------------------- Internals ------------------------------------
  private:
 
  }; // end of class ParametricShapeArcLengthFunctor
 
} // namespace DGtal
 
                                                                        //
 
#endif // !defined ParametricShapeArcLengthFunctor_h
 
#undef ParametricShapeArcLengthFunctor_RECURSES
#endif // else defined(ParametricShapeArcLengthFunctor_RECURSES)