BsplineShape.hpp

// This file is part of the imaging2 class library.
//
// University of Innsbruck, Infmath Imaging, 2009.
// http://infmath.uibk.ac.at
//
// All rights reserved.


#ifndef SHAPE_BSPLINESHAPE_H
#define SHAPE_BSPLINESHAPE_H

#include <shape/ShapeInterface.hpp>
#include <shape/DiscretizableShapeInterface.hpp>

#include <spline/PeriodicBspline.hpp>

namespace imaging
{
  class BsplineShape : public ShapeInterface, public DiscretizableShapeInterface<2>
  {
    PeriodicBspline< ublas::fixed_vector<float_t, 2> > _curve;
    
    class Discretizer;
    
  public:
    const static size_t SHAPE_DIMENSION = 2;

    BsplineShape() : _curve() {}
    
    BsplineShape(const BsplineShape & source) : _curve(source._curve) {}
    
    BsplineShape(const PeriodicBspline< ublas::fixed_vector<float_t, 2> > & curve) : _curve(curve) {}

    const BsplineShape & operator=(const BsplineShape & source)
    { _curve = source._curve; return *this; }
    
    virtual std::auto_ptr< BoundaryDiscretizer<2> > boundary_discretizer(size_t n_points) const;

    void exponential(const ublas::vector<float_t> & vector, ShapeInterface & shape) const;

    void logarithm(const ShapeInterface & shape, ublas::vector<float_t> & vector) const;

    size_t dimension() const { return 2 * _curve.n_coefficients(); }
    
    const PeriodicBspline< ublas::fixed_vector<float_t, 2> > & curve() const { return _curve; }
  };
  
  
  class BsplineShape::Discretizer : public BoundaryDiscretizer<2>
  {
    float_t _step_size;
    const BsplineShape & _spline_curve;

  public:
  
    Discretizer(const BsplineShape  & spline_curve, size_t n_points);
    
    void evaluate(size_t i, ublas::fixed_vector<float_t, 2> & point, ublas::fixed_vector<float_t, 2> & normal, float_t & curvature) const;
  };
}


#endif

---