70{
71 if ( argc < 4 )
72 {
73 usage( argc, argv );
74 return 1;
75 }
76 std::string inputFilename = argv[ 1 ];
77 unsigned int minThreshold = atoi( argv[ 2 ] );
78 unsigned int maxThreshold = atoi( argv[ 3 ] );
79 unsigned int idxP = (argc <= 4) ? 0 : atoi( argv[ 4 ] );
80
87 minThreshold, maxThreshold);
90
91
93 trace.
beginBlock(
"Construct the Khalimsky space from the image domain." );
95 bool space_ok = ks.
init( image.domain().lowerBound(),
96 image.domain().upperBound(), true );
97 if (!space_ok)
98 {
99 trace.
error() <<
"Error in the Khamisky space construction."<<std::endl;
100 return 2;
101 }
104
107 MySurfelAdjacency surfAdj( true );
109
120
122 for ( idxP = idxP % digSurf.size(); idxP != 0; --idxP ) ++it;
123 bel = *it;
125
127 trace.
beginBlock(
"Extracting boundary by distance tracking from an initial bel." );
135 MyDistanceVisitor;
136 typedef MyDistanceVisitor::Node MyNode;
137 typedef MyDistanceVisitor::Scalar MySize;
138
139 SCellEmbedder embedder( ks );
140 Distance distance(2.0);
141 DistanceToPoint distanceToPoint = std::bind( distance, embedder( bel ), std::placeholders::_1 );
142 VertexFunctor vfunctor( embedder, distanceToPoint );
143 MyDistanceVisitor visitor( digSurf, vfunctor, bel );
144
145 unsigned long nbSurfels = 0;
146 MyNode node;
147 while ( ! visitor.finished() )
148 {
149 node = visitor.current();
150 ++nbSurfels;
151 visitor.expand();
152 }
153 MySize maxDist = node.second;
156
159 QApplication application(argc,argv);
161 viewer.show();
163 MyDistanceVisitor visitor2( digSurf, vfunctor, bel );
166 << bel;
167 visitor2.expand();
168 std::vector< MyDistanceVisitor::Node > layer;
169 while ( ! visitor2.finished() )
170 {
171 MyNode n = visitor2.current();
172 Color c = hueShade( n.second );
174 << n.first;
175 visitor2.expand();
176 }
177 viewer << Viewer3D<>::updateDisplay;
178 trace.
info() <<
"nb surfels = " << nbSurfels << std::endl;
180 return application.exec();
182}
Structure representing an RGB triple with alpha component.
Aim: A model of CDigitalSurfaceContainer which defines the digital surface as the boundary of a given...
Aim: A wrapper class around a STL associative container for storing sets of digital points within som...
Aim: Represents a set of n-1-cells in a nD space, together with adjacency relation between these cell...
DigitalSurfaceContainer::SurfelConstIterator ConstIterator
Aim: This class is useful to perform an exploration of a graph given a starting point or set (called ...
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.
Aim: This class is a model of CCellularGridSpaceND. It represents the cubical grid as a cell complex,...
bool init(const Point &lower, const Point &upper, bool isClosed)
Specifies the upper and lower bounds for the maximal cells in this space.
Aim: A model of CDigitalSurfaceContainer which defines the digital surface as the boundary of an impl...
Aim: implements l_p metrics.
Aim: Implements basic operations that will be used in Point and Vector classes.
Component Coordinate
Type for Point elements.
static SCell findABel(const KSpace &K, const PointPredicate &pp, unsigned int nbtries=1000)
Aim: Represent adjacencies between surfel elements, telling if it follows an interior to exterior ord...
void beginBlock(const std::string &keyword="")
Aim: Define a new Functor from the composition of two other functors.
DigitalSurface< MyDigitalSurfaceContainer > MyDigitalSurface
DigitalSetBoundary< KSpace, DigitalSet > MyDigitalSurfaceContainer
Space::RealPoint RealPoint
Aim: A trivial embedder for signed cell, which corresponds to the canonic injection of cell centroids...
ImageContainerBySTLVector< Domain, Value > Type
Aim: Define utilities to convert a digital set into an image.
Modifier class in a Display3D stream. Useful to choose your own mode for a given class....
Represents a signed cell in a cellular grid space by its Khalimsky coordinates and a boolean value.
std::string className() const
Return the style name used for drawing this object.
static ImageContainer importVol(const std::string &filename, const Functor &aFunctor=Functor())
