DGtal  0.9.2
imageGridCurveEstimator.cpp
1 
31 #include <iostream>
33 #include <fstream>
34 #include <algorithm>
36 
38 #include "DGtal/base/Common.h"
39 #include "DGtal/helpers/StdDefs.h"
40 #include "ConfigExamples.h"
42 
44 #include "DGtal/base/BasicFunctors.h"
45 #include "DGtal/kernel/BasicPointPredicates.h"
46 #include "DGtal/io/readers/PGMReader.h"
47 #include "DGtal/images/ImageContainerBySTLVector.h"
49 
51 #include "DGtal/topology/helpers/Surfaces.h"
53 
55 #include "DGtal/geometry/curves/estimation/DSSLengthEstimator.h"
57 
58 //display
59 #include "DGtal/io/boards/Board2D.h"
60 
61 //segmentation
62 #include "DGtal/geometry/curves/GreedySegmentation.h"
63 
65 
66 using namespace DGtal;
67 using namespace DGtal::functors;
68 
69 int main()
70 {
71  //image import
73  std::string filename = examplesPath + "samples/contourS.pgm";
74  Image image = DGtal::PGMReader<Image>::importPGM(filename);
75 
77  //predicate from the image
79  Binarizer b(1, 135);
80  PointFunctorPredicate<Image,Binarizer> predicate(image, b);
82 
84  Z2i::KSpace ks; //Khalimsky space
85  ks.init( image.domain().lowerBound(), image.domain().upperBound(), true );
86  SurfelAdjacency<2> sAdj( true ); //adjacency
88 
90  //extraction of all the contours
91  std::vector< std::vector< Z2i::SCell > > contours;
93  ::extractAll2DSCellContours( contours, ks, sAdj, predicate );
95 
96  if (contours.size() > 0)
97  {
98 
100  //init grid curve from the first retrieved contour
101  Z2i::Curve c;
102  c.initFromSCellsVector( contours.at(1) );
104 
106  //range of points
108  Range r = c.getPointsRange();
110 
112  //length estimation based on a DSS segmentation
114  DSSlength.init(1, r.c(), r.c());
115  trace.info() << "Length: " << DSSlength.eval() << std::endl;
117 
118  //DSS segmentation display
120  typedef GreedySegmentation<SegmentComputer> Segmentation;
121 
122  Segmentation theSegmentation( r.c(), r.c(), SegmentComputer() );
123  Segmentation::SegmentComputerIterator i = theSegmentation.begin();
124  Segmentation::SegmentComputerIterator end = theSegmentation.end();
125 
126  DGtal::Board2D aBoard;
127  aBoard << SetMode("PointVector", "Grid");
128  for ( ; i != end; ++i) {
129  SegmentComputer::Primitive dss = i->primitive();
130  aBoard << SetMode(dss.className(), "Points") << dss;
131  aBoard << SetMode(dss.className(), "BoundingBox") << dss;
132  }
133  aBoard.saveEPS("DisplayDSSSegmentationTuto3.eps");
134 
135  } else trace.info() << "no contour" << std::endl;
136 
137  return 0;
138 
139 }
140 
PointsRange getPointsRange() const
Definition: GridCurve.h:426
Aim: implements association bewteen points lying in a digital domain and values.
Definition: Image.h:69
Trace trace
Definition: Common.h:130
Aim: A utility class for constructing surfaces (i.e. set of (n-1)-cells).
Definition: Surfaces.h:78
Aim: Represent adjacencies between surfel elements, telling if it follows an interior to exterior ord...
functors namespace gathers all DGtal functors.
Aim: This class is a wrapper around ArithmeticalDSS that is devoted to the dynamic recognition of dig...
Aim: Computes the greedy segmentation of a range given by a pair of ConstIterators. The last element of a given segment is the first one one of the next segment.
Aim: model of CBidirectionalRangeFromPoint that adapts any range of elements bounded by two iterators...
static ImageContainer importPGM(const std::string &aFilename, const Functor &aFunctor=Functor(), bool topbotomOrder=true)
Quantity eval() const
Aim: A small functor with an operator () that compares one value to an interval.
GreedySegmentation::SegmentComputerIterator begin() const
bool init(const Point &lower, const Point &upper, bool isClosed)
void init(const double h, const ConstIterator &itb, const ConstIterator &ite)
Aim: a model of CGlobalCurveEstimator that segments the digital curve into DSS and computes the lengt...
DGtal is the top-level namespace which contains all DGtal functions and types.
Aim: model of CConstBidirectionalRange that adapts any range of elements bounded by two iterators [it...
std::ostream & info()
void saveEPS(const char *filename, PageSize size=Board::BoundingBox, double margin=10.0) const
Definition: Board.cpp:805
Modifier class in a Board2D stream. Useful to choose your own mode for a given class. Realizes the concept CDrawableWithBoard2D.
Definition: Board2D.h:247
const Domain & domain() const
Definition: Image.h:192
Aim: describes, in a cellular space of dimension n, a closed or open sequence of signed d-cells (or d...
Definition: GridCurve.h:172
bool initFromSCellsVector(const std::vector< SCell > &aVectorOfSCells)
Aim: This class is a model of CCellularGridSpaceND. It represents the cubical grid as a cell complex...
Aim: The predicate returns true when the predicate returns true for the value assigned to a given poi...
Aim: This class specializes a 'Board' class so as to display DGtal objects more naturally (with <<)...
Definition: Board2D.h:70