geometry/curves/exampleAlphaThickSegmentTgtCover.cpp

This example shows the tangential cover of an noisy contour from the alpha-thick segment primitive with alpha = 4.

See also

Alpha-thick Segment Recognition

Tangential cover with alpha-thick segments with alpha = 4

 
#include <iostream>
#include "ConfigExamples.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/base/Common.h"
#include "DGtal/geometry/curves/AlphaThickSegmentComputer.h"
#include "DGtal/geometry/curves/SaturatedSegmentation.h"
#include "DGtal/base/Circulator.h"
#include "DGtal/geometry/curves/SegmentComputerUtils.h"
 
 
#include "DGtal/io/boards/Board2D.h"
 
 
using namespace std;
using namespace DGtal;
 
 
int main()
{
  trace.beginBlock ( "Example of AlphaThickSegment tangential cover obtained from SaturatedSegmentation" );
 
  typedef std::vector<Z2i::RealPoint>::const_iterator RAConstIterator;
  typedef Circulator<RAConstIterator> ConstCirculator;
  typedef AlphaThickSegmentComputer<Z2i::RealPoint, ConstCirculator> AlphaThickSegmentComputer2D;
  typedef SaturatedSegmentation<AlphaThickSegmentComputer2D> AlphaSegmentation;
 
 
  Board2D aBoard, aBoard2;
  std::string file = examplesPath + "samples/contourSnoisy.sdp";
  std::vector<Z2i::RealPoint> aContour = PointListReader<Z2i::RealPoint>::getPointsFromFile (file);
 
  // displaying contour
  aBoard << SetMode(aContour[0].className(), "Grid"); 
  aBoard2 << 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(200, 200, 200));
  aBoard.fillPolyline(poly);
  aBoard2.setPenColor(DGtal::Color(200, 200, 200));
  aBoard2.fillPolyline(poly);
 
  Circulator<RAConstIterator> circulator (aContour.begin(), aContour.begin(), aContour.end());
 
  AlphaThickSegmentComputer2D computer(4, functions::Hull2D::EuclideanThickness);
  AlphaSegmentation segmentator(circulator, circulator, computer);
  
  AlphaSegmentation::SegmentComputerIterator i = segmentator.begin();
  AlphaSegmentation::SegmentComputerIterator end = segmentator.end();
 
  for ( ; i != end; ++i) {
    AlphaThickSegmentComputer2D current(*i);
    aBoard << SetMode(current.className(), "BoundingBox");
    aBoard << current;
  }
  aBoard.saveEPS("exampleAlphaThickSegmentTgtCover.eps");
  trace.endBlock();
 
  trace.beginBlock ( "Example of AlphaThickSegment tangential cover obtained for a single point" );
 
  // Displaying Tangential cover associated to a single point:
  
  unsigned int index = 80;
  firstMaximalSegment(computer, circulator+index, circulator, circulator);
  AlphaThickSegmentComputer2D first (computer);
  lastMaximalSegment(computer, circulator+index, circulator, circulator);
  AlphaThickSegmentComputer2D last (computer);
  
  aBoard2.setPenColor(DGtal::Color::Blue);
  aBoard2.fillCircle((*(circulator+index))[0], (*(circulator+index))[1], 1.0);
  while(first.end() != last.end()){
    aBoard2 << SetMode(first.className(), "BoundingBox");
    aBoard2 << first;
    nextMaximalSegment(first, circulator);
  }
  aBoard2 << last;
  
  stringstream name;
  name << "exampleAlphaThickSegmentTgtCoverPoint" << index<< ".eps";
  aBoard2.saveEPS(name.str().c_str());
  trace.endBlock();
  
 
 
 
  return 0;
}
//                                                                           //