DGtal  1.0.0
testDistancePropagation.cpp
Go to the documentation of this file.
1 
30 #include <iostream>
32 #include <set>
33 #include "DGtal/base/Common.h"
34 #include "DGtal/kernel/CanonicEmbedder.h"
35 #include "DGtal/helpers/StdDefs.h"
36 #include "DGtal/graph/CUndirectedSimpleGraph.h"
37 #include "DGtal/graph/CGraphVisitor.h"
38 #include "DGtal/graph/GraphVisitorRange.h"
39 #include "DGtal/graph/DistanceBreadthFirstVisitor.h"
40 #include "DGtal/geometry/volumes/distance/LpMetric.h"
41 #include "DGtal/io/boards/Board2D.h"
42 #include "DGtal/io/Color.h"
43 #include "DGtal/io/colormaps/GradientColorMap.h"
44 #include "DGtal/shapes/Shapes.h"
46 
47 
48 using namespace std;
49 using namespace DGtal;
50 using namespace DGtal::concepts;
51 
53 // Functions for testing objects as graph.
55 
58 {
59  typedef Z2i::Space Space;
60  typedef Z2i::Point Point;
61  typedef Z2i::Domain Domain;
62  typedef Z2i::DigitalSet DigitalSet;
63  typedef Z2i::Object4_8 Object;
64 
65  BOOST_CONCEPT_ASSERT(( CUndirectedSimpleGraph<Z2i::Object4_8> ));
66 
67  trace.beginBlock( "Distance propagation in 2D object" );
68  Point p1( -41, -36 );
69  Point p2( 18, 18 );
70  Domain domain( p1, p2 );
71  Point c1( -2, -1 );
72  Point c2( -14, 5 );
73  Point c3( -30, -15 );
74  Point c4( -10, -20 );
75  Point c5( 12, -1 );
76  DigitalSet shape_set( domain );
77 
78  Shapes<Domain>::addNorm2Ball( shape_set, c1, 9 );
79  Shapes<Domain>::addNorm1Ball( shape_set, c2, 9 );
80  Shapes<Domain>::addNorm1Ball( shape_set, c3, 10 );
81  Shapes<Domain>::addNorm2Ball( shape_set, c4, 12 );
82  Shapes<Domain>::addNorm1Ball( shape_set, c5, 4 );
83 
84  Object obj(Z2i::dt4_8, shape_set);
85 
86 
87  GradientColorMap<int> cmap_grad( 0, 25);
88  cmap_grad.addColor( Color( 0, 0, 255 ) );
89  cmap_grad.addColor( Color( 0, 255, 0 ) );
90  cmap_grad.addColor( Color( 255, 0, 0 ) );
91 
92  Board2D board;
93  board << SetMode( domain.className(), "Paving" )
94  << domain
95  << SetMode( p1.className(), "Paving" );
96 
97  Image image = ImageFromSet<Image>::create(shape_set, 1);
98 
99  // Type definitions
100  typedef CanonicEmbedder<Space> VertexEmbedder;
102  typedef RealPoint::Coordinate Scalar;
103  typedef LpMetric<Space> Distance;
104  using DistanceToPoint = std::function<double(const Space::RealPoint &)>;
107 
108  BOOST_CONCEPT_ASSERT(( CGraphVisitor<Visitor> ));
109 
110 
111  VertexEmbedder embedder;
112  Distance distance(2.0);
113  DistanceToPoint distanceToPoint = std::bind(distance, embedder(c1), std::placeholders::_1);
114 
115  VertexFunctor vfunctor( embedder, distanceToPoint );
116  Visitor visitor( obj, vfunctor, c1 );
117 
118  while( ! visitor.finished() )
119  {
120  Scalar v = visitor.current().second;
121  image.setValue( visitor.current().first, v );
122  visitor.expand();
123  }
124 
125  string specificStyle = p1.className() + "/Paving";
126 
127  for ( DigitalSet::ConstIterator it = shape_set.begin();
128  it != shape_set.end();
129  ++it )
130  {
131  if( image(*it) == 0)
132  board << CustomStyle( specificStyle,
133  new CustomColors( Color::Black,
134  Color::Red ) );
135  else if( image(*it) > 0 )
136  board << CustomStyle( specificStyle,
137  new CustomColors( Color::Black,
138  cmap_grad( image(*it) ) ) );
139  else
140  board << CustomStyle( specificStyle,
141  new CustomColors( Color::Black,
142  cmap_grad( 0 ) ) );
143  board << *it;
144  }
145 
146  trace.info() << "- Output file testDistancePropagation.eps" << std::endl;
147  board.saveEPS("testDistancePropagation.eps");
148  trace.endBlock();
149 
150  trace.beginBlock( "Distance visitor as a range." );
151  typedef GraphVisitorRange<Visitor> VisitorRange;
152  VisitorRange range( new Visitor( obj, vfunctor, c1 ) );
153  Scalar d = -1.0;
154  unsigned int nb = 0;
155  unsigned int nbok = 0;
156  unsigned int nbperfect = 0;
157  for ( VisitorRange::NodeConstIterator it = range.beginNode(), itEnd = range.endNode();
158  it != itEnd; ++it )
159  { // Vertex is *it.first
160  Scalar next_d = (*it).second;
161  ++nb; nbok += (next_d >= d-0.75 ) ? 1 : 0;
162  nbperfect += (next_d >= d ) ? 1 : 0;
163  d = next_d;
164  }
165  trace.info() << "(" << nbok << "/" << nb
166  << ") number of vertices in approximate Euclidean distance ordering."<< std::endl;
167  trace.info() << "(" << nbperfect << "/" << nb
168  << ") number of vertices in perfect Euclidean distance ordering."<< std::endl;
169  trace.endBlock();
170  return nb == nbok;
171 }
172 
173 int main( int /*argc*/, char** /*argv*/ )
174 {
175  bool res = testDistancePropagation();
176  return res ? 0 : 1;
177 }
178 
179 
180 
181 
Aim: This class is useful to perform an exploration of a graph given a starting point or set (called ...
Aim: An object (or digital object) represents a set in some digital space associated with a digital t...
Definition: Object.h:119
void beginBlock(const std::string &keyword="")
int main(int, char **)
HyperRectDomain< Space > Domain
Aim: This class template may be used to (linearly) convert scalar values in a given range into a colo...
Aim: implements association bewteen points lying in a digital domain and values.
Definition: Image.h:69
MyDigitalSurface::ConstIterator ConstIterator
Aim: implements l_p metrics.
Definition: LpMetric.h:74
Trace trace
Definition: Common.h:144
Component Coordinate
Type for Point elements.
Definition: PointVector.h:617
Aim: Define utilities to convert a digital set into an image.
Definition: ImageFromSet.h:63
void addColor(const Color &color)
const Node & current() const
ImageSelector< Z2i::Domain, int >::Type Image
double endBlock()
Custom style class redefining the pen color and the fill color. You may use Board2D::Color::None for ...
Definition: Board2D.h:278
BreadthFirstVisitor< MyDigitalSurface > Visitor
Aim: Gathers several functions useful for concept checks.
Domain domain
bool testDistancePropagation()
PointVector< 3, double > RealPoint
Aim: Define a new Functor from the composition of two other functors.
void saveEPS(const char *filename, PageSize size=Board::BoundingBox, double margin=10.0) const
Definition: Board.cpp:805
std::string className() const
DGtal is the top-level namespace which contains all DGtal functions and types.
MyPointD Point
Definition: testClone2.cpp:383
Aim: A wrapper class around a STL associative container for storing sets of digital points within som...
Aim: Transforms a graph visitor into a single pass input range.
Aim: A trivial embedder for digital points, which corresponds to the canonic injection of Zn into Rn.
Aim: Represents the concept of local graph: each vertex has neighboring vertices, but we do not neces...
std::ostream & info()
Aim: This class is useful to perform a breadth-first exploration of a graph given a starting point or...
Modifier class in a Board2D stream. Useful to choose your own mode for a given class....
Definition: Board2D.h:247
Structure representing an RGB triple with alpha component.
Definition: Color.h:66
Aim: A utility class for constructing different shapes (balls, diamonds, and others).
Image image(domain)
Aim: This class specializes a 'Board' class so as to display DGtal objects more naturally (with <<)....
Definition: Board2D.h:70
Aim: Defines the concept of a visitor onto a graph, that is an object that traverses vertices of the ...