imageGridCurveEstimator.cpp

#include <iostream>
#include <fstream>
#include <algorithm>

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

#include "DGtal/base/BasicFunctors.h"
#include "DGtal/kernel/BasicPointPredicates.h"
#include "DGtal/io/readers/PNMReader.h"
#include "DGtal/images/ImageContainerBySTLVector.h"

#include "DGtal/topology/helpers/Surfaces.h"

#include "DGtal/geometry/curves/estimation/DSSLengthEstimator.h"

//display
#include "DGtal/io/boards/Board2D.h"

//segmentation
#include "DGtal/geometry/curves/GreedySegmentation.h"


int main()
{
  //image import
  typedef DGtal::ImageContainerBySTLVector< Z2i::Domain, int> Image;
  std::string filename =  examplesPath + "samples/contourS.pgm";
  Image image = DGtal::PNMReader<Image>::importPGM(filename); 

  //predicate from the image
  typedef IntervalThresholder<Image::Value> Binarizer; 
  Binarizer b(1, 135); 
  PointFunctorPredicate<Image,Binarizer> predicate(image, b); 

  Z2i::KSpace ks;                                            //Khalimsky space 
  ks.init( image.domain().lowerBound(), image.domain().upperBound(), true );
  SurfelAdjacency<2> sAdj( true );                           //adjacency

  //extraction of all the contours
  std::vector< std::vector< Z2i::SCell > > contours;
  Surfaces<Z2i::KSpace>
    ::extractAll2DSCellContours( contours, ks, sAdj, predicate );

  if (contours.size() > 0)
  {
    
    //init grid curve from the first retrieved contour
    Z2i::Curve c;
    c.initFromSCellsVector( contours.at(1) );  

    //range of points
    typedef Z2i::Curve::PointsRange Range; 
    Range r = c.getPointsRange(); 

    //length estimation based on a DSS segmentation
    DSSLengthEstimator< Range::ConstIterator > DSSlength;
    DSSlength.init(1, r.begin(), r.end(), c.isClosed());
    double length = DSSlength.eval();
    trace.info() << "Length: " << length << endl; 
    
    //DSS segmentation display
    typedef Z2i::Curve::PointsRange::ConstCirculator ConstCirculator; 
    typedef ArithmeticalDSS<ConstCirculator,int,4> SegmentComputer;
    typedef GreedySegmentation<SegmentComputer> Segmentation;

    Segmentation theSegmentation( r.c(), r.c(), SegmentComputer() );
    Segmentation::SegmentComputerIterator i = theSegmentation.begin();
    Segmentation::SegmentComputerIterator end = theSegmentation.end();
    
    DGtal::Board2D aBoard;
    aBoard << SetMode("PointVector", "Grid");
    for ( ; i != end; ++i) {
      aBoard << SetMode(i->className(), "Points") << *i; 
      aBoard << SetMode(i->className(), "BoundingBox") << *i; 
    } 
    aBoard.saveEPS("DisplayDSSSegmentationTuto3.eps");
  
  } else trace.info() << "no contour" << endl; 
  
  return 0;

}