FPLengthEstimator.ih

#include <cstdlib>

// IMPLEMENTATION of inline methods.

// ----------------------- Standard services ------------------------------

template <typename T>
inline
DGtal::FPLengthEstimator<T>::FPLengthEstimator()
{
  myIsInitBefore = false;
}

template <typename T>
inline
DGtal::FPLengthEstimator<T>::~FPLengthEstimator()
{
}



// Interface - public :



template <typename T>
inline
void
DGtal::FPLengthEstimator<T>::init(const double h, 
     const ConstIterator& itb, const ConstIterator& ite, const bool& isClosed)
{
  myH = h;
  myIsInitBefore = true;
  myRep.clear();

  FaithfulPolygon fp( itb, ite, isClosed ); 
  myRep.resize( fp.size() );
  fp.copyFP( myRep.begin() );

  if ( (isClosed)&&(myRep.size()!=0) )
    myRep.push_back(myRep.front());

}

template <typename T>
inline
typename DGtal::FPLengthEstimator<T>::Quantity
DGtal::FPLengthEstimator<T>::eval() const
{
  ASSERT(myH > 0);
  ASSERT(myIsInitBefore);

  Quantity val = 0;
  
  if (myRep.size() > 1) {

    typename std::vector<Point>::const_iterator i = myRep.begin();
    typename std::vector<Point>::const_iterator j = i;
    ++j; 
    for ( ; j != myRep.end(); ++i, ++j) {
      Vector v( *j - *i ); 
      val += v.norm(Vector::L_2);
    }

  }
  
  return val*myH;
}





template <typename T>
inline
void
DGtal::FPLengthEstimator<T>::selfDisplay ( std::ostream & out ) const
{
  out << "[FPLengthEstimator]";
  if (myIsInitBefore) 
    out <<" myH="<<myH;
  else
    out<< " not initialized";
}

template <typename T>
inline
bool
DGtal::FPLengthEstimator<T>::isValid() const
{
    return myIsInitBefore;
}



// Implementation of inline functions                                        //

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

//                                                                           //