ShapeGeometricFunctors.h

#pragma once

#if defined(ShapeGeometricFunctors_RECURSES)
#error Recursive header files inclusion detected in ShapeGeometricFunctors.h
#else // defined(ShapeGeometricFunctors_RECURSES)

#define ShapeGeometricFunctors_RECURSES

#if !defined ShapeGeometricFunctors_h

#define ShapeGeometricFunctors_h

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

namespace DGtal
{
  namespace functors
  {
    
  namespace ShapeGeometricFunctors 
  {

    template <typename TShape>
    struct ShapePositionFunctor {
      typedef TShape Shape;
      typedef typename Shape::RealPoint RealPoint;
      typedef RealPoint Argument;
      typedef RealPoint Quantity;
      typedef Quantity Value;
     
      ShapePositionFunctor( ConstAlias<Shape> aShape = 0 ) : myShape( aShape ) {}
      
      void attach( ConstAlias<Shape> aShape )
      {
        myShape = aShape;
      }

      Quantity operator()( const RealPoint & p ) const
      {
        return p;
      }

    private:
      CountedConstPtrOrConstPtr<Shape> myShape;
    };

    template <typename TShape>
    struct ShapeNormalVectorFunctor {
      typedef TShape Shape;
      typedef typename Shape::RealPoint RealPoint;
      typedef typename Shape::RealVector RealVector;
      typedef RealPoint Argument;
      typedef typename RealVector::Component Scalar;
      typedef RealVector Quantity;
      typedef Quantity Value;

      ShapeNormalVectorFunctor( ConstAlias<Shape> aShape = 0 ) : myShape( aShape ) {}
      
      void attach( ConstAlias<Shape> aShape )
      {
        myShape = aShape;
      }

      Quantity operator()( const RealPoint & p ) const
      {
        RealVector v = myShape->gradient( p );
        Scalar norm = v.norm();
        return ( norm != 0 ) ? v / norm : v;
      }

    private:
      CountedConstPtrOrConstPtr<Shape> myShape;
    };


    template <typename TShape>
    struct ShapeMeanCurvatureFunctor {
      typedef TShape Shape;
      typedef typename Shape::RealPoint RealPoint;
      typedef typename Shape::RealVector RealVector;
      typedef typename RealVector::Component Scalar;
      typedef RealPoint Argument;
      typedef Scalar Quantity;
      typedef Quantity Value;
      
      ShapeMeanCurvatureFunctor( ConstAlias<Shape> aShape = 0 ) : myShape( aShape ) {}
      
      void attach( ConstAlias<Shape> aShape )
      {
        myShape = aShape;
      }

      Quantity operator()( const RealPoint & p ) const
      {
        return myShape->meanCurvature( p );
      }

    private:
      CountedConstPtrOrConstPtr<Shape> myShape;
    };
    
  } // namespace ShapeGeometricFunctors
  }
} // namespace DGtal


// Includes inline functions.
#include "DGtal/shapes/ShapeGeometricFunctors.ih"

//                                                                           //

#endif // !defined ShapeGeometricFunctors_h

#undef ShapeGeometricFunctors_RECURSES
#endif // else defined(ShapeGeometricFunctors_RECURSES)