geometry/curves/estimation/exampleLMST2D.cpp

This example shows the basic usage of the Lambda maximal segment tangent estimation in 2D.

See also

Digital straight lines and segments

 
#include <iostream>
#include <iterator>
#include <vector>
#include "DGtal/base/Common.h"
#include "DGtal/base/BasicTypes.h"
#include "DGtal/helpers/StdDefs.h"
 
#include "DGtal/geometry/curves/ArithmeticalDSSComputer.h"
#include "DGtal/geometry/curves/SaturatedSegmentation.h"
#include "DGtal/geometry/curves/estimation/LambdaMST2D.h"
 
 
using namespace std;
using namespace DGtal;
using namespace Z2i;
 
 
int main()
{
  typedef vector < Point > Container;
  typedef Container::const_iterator ConstIterator;
  typedef ArithmeticalDSSComputer < ConstIterator, int, 8 > SegmentComputer;
  typedef SaturatedSegmentation<SegmentComputer> Segmentation;
  LambdaMST2D < Segmentation > lmst;
  
  // Input points
  Container contour;
  contour.push_back(Point(18,18));
  contour.push_back(Point(17,19));
  contour.push_back(Point(16,20));
  contour.push_back(Point(15,21));
  contour.push_back(Point(14,22));
  contour.push_back(Point(13,23));
  contour.push_back(Point(12,24));
  contour.push_back(Point(11,25));
  contour.push_back(Point(10,26));
  contour.push_back(Point(9,27));
  contour.push_back(Point(8,28));
  
  // Initialization of tangential cover
  Segmentation segmenter ( contour.begin(), contour.end(), SegmentComputer() );
  lmst.attach ( segmenter );
  
  for ( ConstIterator it = contour.begin(); it != contour.end(); ++it )
    lmst.eval ( *it );
  
  lmst.init ( contour.begin(), contour.end() );
  std::vector < RealVector > tangent;
  lmst.eval < back_insert_iterator< vector < RealVector > > > ( contour.begin(), contour.end(),  back_inserter ( tangent ) );
  
  return 0;
}
//                                                                           //