greedyAlphaThickDecomposition.cpp

Go to the documentation of this file.

 
 

#include <iostream>
#include "ConfigExamples.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/base/Common.h"
#include "DGtal/geometry/curves/AlphaThickSegmentComputer.h"
#include "DGtal/io/boards/Board2D.h"
#include "DGtal/geometry/curves/GreedySegmentation.h"
#include <DGtal/io/readers/GenericReader.h>
 
using namespace std;
using namespace DGtal;

 
int main()
{
  trace.beginBlock ( "Example of greedy alpha thick segment  decompotion" );
 
  typedef  std::vector<Z2i::RealPoint>::const_iterator ConstIterator;
  typedef  AlphaThickSegmentComputer<Z2i::RealPoint, ConstIterator > AlphaThickSegmentComputer2D;
 
  Board2D aBoard;
  std::string file = examplesPath + "samples/contourSnoisy.sdp";
  std::vector<Z2i::RealPoint> aContour = PointListReader<Z2i::RealPoint>::getPointsFromFile (file);
 
  typedef GreedySegmentation<AlphaThickSegmentComputer2D> DecompositionAT;
 
  // displaying contour
  aBoard << SetMode(aContour[0].className(), "Grid");
  std::vector<LibBoard::Point> poly;
  for (unsigned int i = 0; i< aContour.size(); i++)  poly.push_back(LibBoard::Point(aContour[i][0], aContour[i][1]));
  aBoard.setPenColor(DGtal::Color::Gray);
  aBoard.fillPolyline(poly);
 
  // Computing greedy Alpha Thick decomposition.
  aBoard << SetMode("AlphaThickSegment", "BoundingBox");

  DecompositionAT theDecomposition(aContour.begin(), aContour.end(), AlphaThickSegmentComputer2D(4));

  for ( DecompositionAT::SegmentComputerIterator
           it = theDecomposition.begin(),
           itEnd = theDecomposition.end();
         it != itEnd; ++it )
     {
       aBoard << CustomStyle( (*it).className(),
                              new CustomPenColor( Color::Blue ) );
       aBoard<< *it;
     }
 
  aBoard.saveEPS("greedyAlphaThickDecomposition.eps");
 
  trace.endBlock();
  return 0;
}
//                                                                           //