TwoStepLocalLengthEstimator.ih

#include <cstdlib>

// IMPLEMENTATION of inline methods.

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



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



// Interface - public :



template <typename T>
inline
void
DGtal::TwoStepLocalLengthEstimator<T>::init(const double h, 
                                            const ConstIterator& itb, 
                                            const ConstIterator& ite, 
                                            const bool& )
{
  myH = h;
  myBeginIt = itb;
  myEndIt = ite;
  myIsInitBefore = true;
}

template <typename T>
inline
typename DGtal::TwoStepLocalLengthEstimator<T>::Quantity
DGtal::TwoStepLocalLengthEstimator<T>::eval() const
{
  ASSERT(myH > 0);
  ASSERT(myIsInitBefore);
  ASSERT(myBeginIt != myEndIt);

  unsigned int nbDirect=0;
  unsigned int nbDiag=0;

  ConstIterator i = myBeginIt, ii = myBeginIt;
  ++ii;
 
  while (( i != myEndIt) && (ii != myEndIt))
    {
      if (((*i).second).dot((*ii).second) == 0)
        nbDiag++;
      else
        nbDirect+=2;

      ++i; 
      if (i == myEndIt) 
        continue;
      ++i;

      ++ii; 
      if (ii == myEndIt) 
        continue;
      ++ii;
    }   
  
  if ( i!= myEndIt)
    nbDirect ++;
  
  return (nbDiag*myWeightDiagonal + nbDirect*myWeightDirect)*myH;
}







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

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



// Implementation of inline functions                                        //

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

//                                                                           //