exampleGeometricalDCA.cpp

#include <iostream>

#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "ConfigExamples.h"

#include "DGtal/geometry/curves/GeometricalDCA.h"


using namespace std;
using namespace DGtal;
using namespace Z2i; 


int main( int argc, char** argv )
{
  trace.beginBlock ( "Example for GeometricalDCA" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << endl;
  
  std::string filename = examplesPath + "samples/DCA.dat";
  ifstream instream; // input stream
  instream.open (filename.c_str(), ifstream::in);
  
  Curve c; //grid curve
  c.initFromVectorStream(instream);


  trace.beginBlock("Simple example");

  Curve::IncidentPointsRange r = c.getIncidentPointsRange(); 
  Curve::IncidentPointsRange::ConstIterator itEnd (r.end()); 

  GeometricalDCA<Curve::IncidentPointsRange::ConstIterator> s; 
  //extension
  s.init( r.begin() );
  while ( ( s.end() != itEnd )
        &&( s.extendForward() ) ) {}

  trace.info() << s << endl;  

  trace.endBlock();


  return 0;
}
//                                                                           //