Ellipse2D.h

#pragma once

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

#define Ellipse2D_RECURSES

#if !defined Ellipse2D_h

#define Ellipse2D_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/shapes/parametric/StarShaped2D.h"
#include <cmath>

namespace DGtal
{

  // template class Ellipse2D
  template <typename TSpace>
  class Ellipse2D:  public StarShaped2D<TSpace>
  {
    // ----------------------- Standard services ------------------------------
  public:

    typedef TSpace Space;
    typedef typename Space::Point Point;
    typedef typename Space::RealPoint RealPoint2D;
    typedef typename Space::RealVector RealVector2D;
      
    ~Ellipse2D();
    
    Ellipse2D( const double x0, const double y0, 
         const double a0, const double a1, const double theta);

    Ellipse2D(const RealPoint2D &aPoint,
        const double a0, const double a1, const double theta);

    Ellipse2D(const Point &aPoint,
        const double a0, const double a1, const double theta);

    
  // ------------- Implementation of 'StarShaped' services ------------------
  public:

    RealPoint2D getLowerBound() const
    {
      return RealPoint2D(myCenter[0] - myAxis1, myCenter[1] - myAxis1);
    }

    RealPoint2D getUpperBound() const
    {
      return RealPoint2D(myCenter[0] + myAxis1, myCenter[1] + myAxis1);
    }

    RealPoint2D center() const
    {
      return myCenter;
    }
   
    double parameter( const RealPoint2D & p ) const;


    RealPoint2D x( const double t ) const;

    RealVector2D xp( const double t ) const;

    RealVector2D xpp( const double t ) const;
    

    // ------------------------- data ----------------------------
  private:

    RealPoint2D myCenter;
    
    double myAxis1;


    double myAxis2;
    
    double myTheta;

    // ----------------------- Interface --------------------------------------
  public:

    void selfDisplay ( std::ostream & out ) const;

    bool isValid() const;


    // ------------------------- Hidden services ------------------------------
  protected:

    Ellipse2D();

  private:

    //  Ellipse2D ( const Ellipse2D & other );

    Ellipse2D & operator= ( const Ellipse2D & other );

    // ------------------------- Internals ------------------------------------
  private:

  }; // end of class Ellipse2D


  template <typename T>
  std::ostream&
  operator<< ( std::ostream & out, const Ellipse2D<T> & object );

} // namespace DGtal


// Includes inline functions.
#include "DGtal/shapes/parametric/Ellipse2D.ih"

//                                                                           //

#endif // !defined Ellipse2D_h

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