DSSLengthEstimator.ih

#include <cstdlib>

// IMPLEMENTATION of inline methods.

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

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

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



// Interface - public :



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

  //segments into DSS
  DSSComputer computer; 
  deprecated::GreedyDecomposition<DSSComputer> decomposition ( itb, ite, computer, isClosed ); 

  //computes the resulting polygonal representation
  typename deprecated::GreedyDecomposition<DSSComputer>::SegmentIterator segIt = decomposition.begin(); 
  for ( ; segIt != decomposition.end(); ++segIt ) {
    myRep.push_back( *segIt.getBack() );
  } 
  //last point
  ConstIterator i( segIt.getFront() ); 
  myRep.push_back(*--i); 
}

template <typename T>
inline
typename DGtal::DSSLengthEstimator<T>::Quantity
DGtal::DSSLengthEstimator<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::DSSLengthEstimator<T>::selfDisplay ( std::ostream & out ) const
{
  out << "[DSSLengthEstimator]";
  if (myIsInitBefore) 
    out <<" myH="<<myH;
  else
    out<< " not initialized";
}

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



// Implementation of inline functions                                        //

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

//                                                                           //