CircleFrom3Points.h

#pragma once

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

#define CircleFrom3Points_RECURSES

#if !defined CircleFrom3Points_h

#define CircleFrom3Points_h

// Inclusions
#include <iostream>

#include "DGtal/base/Common.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/kernel/NumberTraits.h"
#include "DGtal/io/Color.h"

namespace DGtal
{


  // template class CircleFrom3Points
  template <typename TPoint>
  class CircleFrom3Points
  {

    BOOST_STATIC_ASSERT(( TPoint::dimension == 2 ));
    
    // ----------------------- associated types ------------------------------
  public:

    typedef typename TPoint::Coordinate Coordinate;
    typedef Coordinate Distance; //to promote
    typedef TPoint Point;
    typedef TPoint Vector;

    // ----------------------- Standard services ------------------------------
  public:

    CircleFrom3Points();


    CircleFrom3Points(const Point& aFirstPoint, const Point& aSecondPoint, const Point& aThirdPoint);

    void init(const Point& aFirstPoint, const Point& aSecondPoint, const Point& aThirdPoint);

    CircleFrom3Points ( const CircleFrom3Points & other );

    CircleFrom3Points & operator= ( const CircleFrom3Points & other );


    ~CircleFrom3Points();

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

    Distance signedDistance(const Point& aP) const;

    //------------------ accessors -------------------------------
    bool isValid() const;
  
    void getParameters(double& cx, double& cy, double& rr) const;

    double getCurvature() const;
    
    const Point & p() const
    {
      return myP;
    };
    
    const Point & q() const
    {
      return myQ;
    };
    
    const Point & r() const
    {
      return myR;
    };


    //------------------ display -------------------------------

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

    
    // ------------------------- Protected Datas ------------------------------
  private:
    // ------------------------- Private Datas --------------------------------
  private:
    //the three points that uniquely define the circle
    Point myP;
    Point myQ;
    Point myR;
    // ------------------------- Hidden services ------------------------------
  protected:


  private:



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




  }; // end of class CircleFrom3Points


  template <typename TPoint>
  inline
  std::ostream&
  operator<< ( std::ostream & out, 
        const CircleFrom3Points<TPoint> & object )
  {
    object.selfDisplay( out );
    return out;
  }


} // namespace DGtal


// Includes inline functions.
#include "DGtal/shapes/fromPoints/CircleFrom3Points.ih"

//                                                                           //

#endif // !defined CircleFrom3Points_h

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