DGtal 1.4.0
Loading...
Searching...
No Matches
Typedefs | Functions
fullConvexityShortestPaths3D.cpp File Reference
DGtal Examples
#include <iostream>
#include <queue>
#include "DGtal/base/Common.h"
#include "DGtal/io/viewers/Viewer3D.h"
#include "DGtal/io/DrawWithDisplay3DModifier.h"
#include "DGtal/io/Color.h"
#include "DGtal/io/colormaps/SimpleDistanceColorMap.h"
#include "DGtal/shapes/Shapes.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/helpers/Shortcuts.h"
#include "DGtal/images/ImageContainerBySTLVector.h"
#include "DGtal/geometry/volumes/TangencyComputer.h"
#include "ConfigExamples.h"
Include dependency graph for fullConvexityShortestPaths3D.cpp:

Go to the source code of this file.

Typedefs

typedef Z3i::Space Space
 
typedef Z3i::KSpace KSpace
 
typedef Z3i::Domain Domain
 
typedef Z3i::SCell SCell
 
typedef Shortcuts< KSpace > SH3
 
typedef Space::Point Point
 
typedef Space::RealPoint RealPoint
 
typedef Space::Vector Vector
 

Functions

int reaction (void *viewer, DGtal::int32_t name, void *data)
 
int main (int argc, char **argv)
 

Detailed Description

This program is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.

Author
Jacques-Olivier Lachaud (jacqu.nosp@m.es-o.nosp@m.livie.nosp@m.r.la.nosp@m.chaud.nosp@m.@uni.nosp@m.v-sav.nosp@m.oie..nosp@m.fr ) Laboratory of Mathematics (CNRS, UMR 5127), University of Savoie, France
Date
2021/06/20

An example file named fullConvexityShortestPaths3D

This file is part of the DGtal library.

Definition in file fullConvexityShortestPaths3D.cpp.

Typedef Documentation

◆ Domain

typedef Z3i::Domain Domain

Definition at line 82 of file fullConvexityShortestPaths3D.cpp.

◆ KSpace

typedef Z3i::KSpace KSpace

Definition at line 81 of file fullConvexityShortestPaths3D.cpp.

◆ Point

typedef Space::Point Point

Definition at line 85 of file fullConvexityShortestPaths3D.cpp.

◆ RealPoint

typedef Space::RealPoint RealPoint

Definition at line 86 of file fullConvexityShortestPaths3D.cpp.

◆ SCell

typedef Z3i::SCell SCell

Definition at line 83 of file fullConvexityShortestPaths3D.cpp.

◆ SH3

typedef Shortcuts< KSpace > SH3

Definition at line 84 of file fullConvexityShortestPaths3D.cpp.

◆ Space

typedef Z3i::Space Space

Definition at line 80 of file fullConvexityShortestPaths3D.cpp.

◆ Vector

typedef Space::Vector Vector

Definition at line 87 of file fullConvexityShortestPaths3D.cpp.

Function Documentation

◆ main()

int main ( int argc,
char ** argv )

[Tangency3D-shortest-paths]

[Tangency3D-shortest-paths]

[Tangency3D-shortest-path]

[Tangency3D-shortest-path]

Definition at line 98 of file fullConvexityShortestPaths3D.cpp.

99{
100 trace.info() << "Usage: " << argv[ 0 ] << " <input.vol> <m> <M> <opt>" << std::endl;
101 trace.info() << "\tComputes shortest paths to a source point" << std::endl;
102 trace.info() << "\t- input.vol: choose your favorite shape" << std::endl;
103 trace.info() << "\t- m [==0], M [==255]: used to threshold input vol image" << std::endl;
104 trace.info() << "\t- opt >= sqrt(3): secure shortest paths, 0: fast" << std::endl;
105 string inputFilename = examplesPath + "samples/Al.100.vol";
106 std::string fn= argc > 1 ? argv[ 1 ] : inputFilename; //< vol filename
107 int m = argc > 2 ? atoi( argv[ 2 ] ) : 0; //< low for thresholding
108 int M = argc > 3 ? atoi( argv[ 3 ] ) : 255; //< up for thresholding
109 double opt = argc > 4 ? atof( argv[ 4 ] ) : sqrt(3.0); //< exact (sqrt(3)) or inexact (0) computations
110
111 QApplication application(argc,argv);
112 Viewer3D<> viewer;
113 viewer.setWindowTitle("fullConvexityShortestPaths3D");
114 viewer.show();
115
116 // Set up shortcuts parameters.
117 auto params = SH3::defaultParameters();
118 params( "thresholdMin", m )( "thresholdMax", M );
119 params( "surfaceComponents" , "All" );
120
121 // Domain creation from two bounding points.
122 trace.info() << "Building set or importing vol ... ";
123 KSpace K;
124 auto bimage = SH3::makeBinaryImage( fn, params );
125 K = SH3::getKSpace( bimage );
126 trace.info() << " [Done]" << std::endl;
127
128 // Compute surface
129 auto surface = SH3::makeDigitalSurface( bimage, K, params );
130
131 // Compute interior boundary points
132 // They are less immediate interior points than surfels.
133 std::vector< Point > points;
134 std::map< SCell, int > surfel2idx;
135 std::map< Point, int > point2idx;
136 int idx = 0;
137 for ( auto s : (*surface) )
138 {
139 // get inside point on the border of the shape.
140 Dimension k = K.sOrthDir( s );
141 auto voxel = K.sIncident( s, k, K.sDirect( s, k ) );
142 Point p = K.sCoords( voxel );
143 auto it = point2idx.find( p );
144 if ( it == point2idx.end() )
145 {
146 points.push_back( p );
147 surfel2idx[ s ] = idx;
148 point2idx [ p ] = idx++;
149 }
150 else
151 surfel2idx[ s ] = it->second;
152 }
153 trace.info() << "Shape has " << points.size() << " interior boundary points"
154 << std::endl;
155
156 // Select a starting point.
157 typedef Viewer3D<> MViewer3D;
158 DGtal::int32_t selected_surfels[ 2 ] = { 0, 0 };
159 auto surfels = SH3::getSurfelRange ( surface );
160 for ( int i = 0; i < 2; i++ )
161 {
162 MViewer3D viewerCore( K );
163 viewerCore.show();
164 Color colSurfel( 200, 200, 255, 255 );
165 Color colStart( 255, 0, 0, 255 );
166 DGtal::int32_t name = 0;
167 viewerCore << SetMode3D( surfels[ 0 ].className(), "Basic");
168 viewerCore.setFillColor( colSurfel );
169 for ( auto && s : surfels ) viewerCore << SetName3D( name++ ) << s;
170 viewerCore << SetSelectCallback3D( reaction, &selected_surfels[ i ],
171 0, surfels.size() - 1 );
172 viewerCore << MViewer3D::updateDisplay;
173 application.exec();
174 }
175
176 // Get selected surfel/point
177 const auto s0 = surfels[ selected_surfels[ 0 ] ];
178 Dimension k0 = K.sOrthDir( s0 );
179 auto voxel0 = K.sIncident( s0, k0, K.sDirect( s0, k0 ) );
180 Point p0 = K.sCoords( voxel0 );
181 auto start0 = point2idx[ p0 ];
182 std::cout << "Start0 index is " << start0 << std::endl;
183 const auto s1 = surfels[ selected_surfels[ 1 ] ];
184 Dimension k1 = K.sOrthDir( s1 );
185 auto voxel1 = K.sIncident( s1, k1, K.sDirect( s1, k1 ) );
186 Point p1 = K.sCoords( voxel1 );
187 auto start1 = point2idx[ p1 ];
188 std::cout << "Start1 index is " << start1 << std::endl;
189
190 // (I) Extracts shortest paths to a target
191
192 typedef TangencyComputer< KSpace >::Index Index;
194 TangencyComputer< KSpace > TC( K );
195 TC.init( points.cbegin(), points.cend() );
196 auto SP = TC.makeShortestPaths( opt );
197 SP.init( start0 ); //< set source
198 double last_distance = 0.0;
199 while ( ! SP.finished() )
200 {
201 last_distance = std::get<2>( SP.current() );
202 SP.expand();
203 }
204 std::cout << "Max distance is " << last_distance << std::endl;
206
207 {
208 const int nb_repetitions = 10;
209 const double period = last_distance / nb_repetitions;
210 SimpleDistanceColorMap< double > cmap( 0.0, period, false );
211 MViewer3D viewerCore;
212 viewerCore.show();
213 Color colSurfel( 200, 200, 255, 128 );
214 Color colStart( 255, 0, 0, 128 );
215
216 viewerCore.setUseGLPointForBalls(true);
217 for ( auto i = 0; i < points.size(); ++i )
218 {
219 const double d_s = SP.distance( i );
220 Color c_s = cmap( fmod( d_s, period ) );
221 viewerCore.setFillColor( c_s );
222 viewerCore.addBall( RealPoint( points[ i ][ 0 ],
223 points[ i ][ 1 ],
224 points[ i ][ 2 ] ), 12.0 );
225 }
226
227 // JOL: Left if you wish to display it as a surface, and to display paths.
228
229 // auto surfels = SH3::getSurfelRange ( surface );
230 // viewerCore << SetMode3D( surfels[ 0 ].className(), "Basic");
231 // for ( auto && s : surfels )
232 // {
233 // const double d_s = SP.distance( surfel2idx[ s ] );
234 // Color c_s = cmap( fmod( d_s, period ) );
235 // viewerCore.setFillColor( c_s );
236 // viewerCore << s;
237 // }
238 // viewerCore.setFillColor( colStart );
239 // viewerCore.setLineColor( Color::Green );
240 // viewerCore << s;
241 // for ( Index i = 0; i < SP.size(); i++ ) {
242 // Point p1 = SP.point( i );
243 // Point p2 = SP.point( SP.ancestor( i ) );
244 // viewerCore.addLine( p1, p2, 1.0 );
245 // }
246 viewerCore << MViewer3D::updateDisplay;
247 application.exec();
248 }
249
250 // (II) Extracts a shortest path between two points.
251
253 auto SP0 = TC.makeShortestPaths( opt );
254 auto SP1 = TC.makeShortestPaths( opt );
255 SP0.init( start0 ); //< source point
256 SP1.init( start1 ); //< target point
257 std::vector< Index > Q; //< the output shortest path
258 while ( ! SP0.finished() && ! SP1.finished() )
259 {
260 auto n0 = SP0.current();
261 auto n1 = SP1.current();
262 auto p0 = std::get<0>( n0 );
263 auto p1 = std::get<0>( n1 );
264 SP0.expand();
265 SP1.expand();
266 if ( SP0.isVisited( p1 ) )
267 {
268 auto c0 = SP0.pathToSource( p1 );
269 auto c1 = SP1.pathToSource( p1 );
270 std::copy(c0.rbegin(), c0.rend(), std::back_inserter(Q));
271 Q.pop_back();
272 std::copy(c1.begin(), c1.end(), std::back_inserter(Q));
273 break;
274 }
275 }
276 // Q is empty if there is no path.
278
279 // Display computed distances and shortest path
280 {
281 const int nb_repetitions = 10;
282 const double period = last_distance / nb_repetitions;
283 SimpleDistanceColorMap< double > cmap( 0.0, period, false );
284 MViewer3D viewerCore;
285 viewerCore.show();
286 Color colSurfel( 200, 200, 255, 128 );
287 Color colStart( 255, 0, 0, 128 );
288 viewerCore.setUseGLPointForBalls(true);
289 for ( auto i = 0; i < points.size(); ++i )
290 {
291 const double d_s0 = SP0.isVisited( i ) ? SP0.distance( i ) : SP0.infinity();
292 const double d_s1 = SP1.isVisited( i ) ? SP1.distance( i ) : SP1.infinity();
293 const double d_s = std::min( d_s0, d_s1 );
294 Color c_s = ( d_s != SP0.infinity() )
295 ? cmap( fmod( d_s, period ) )
296 : Color::Black;
297 viewerCore.setFillColor( c_s );
298 viewerCore.addBall( RealPoint( points[ i ][ 0 ],
299 points[ i ][ 1 ],
300 points[ i ][ 2 ] ), 12 );
301 }
302 viewerCore.setLineColor( Color::Green );
303 for ( auto i = 1; i < Q.size(); i++ )
304 {
305 Point p1 = TC.point( Q[ i-1 ] );
306 Point p2 = TC.point( Q[ i ] );
307 viewerCore.addLine( p1, p2, 18.0 );
308 }
309 viewerCore << MViewer3D::updateDisplay;
310 application.exec();
311 }
312
313 // (III) Extracts multiple shortest paths between many sources and two targets.
314
315 std::vector< Index > sources;
316 std::vector< Index > dests;
317 for ( int i = 0; i < 20; i++ )
318 sources.push_back( rand() % TC.size() );
319 dests.push_back( start0 );
320 dests.push_back( start1 );
321 auto paths = TC.shortestPaths( sources, dests, opt );
322
323 // Display them.
324 {
325 MViewer3D viewerCore;
326 viewerCore.show();
327 Color colSurfel( 200, 200, 255, 128 );
328 Color colStart( 255, 0, 0, 128 );
329 viewerCore.setUseGLPointForBalls(true);
330 for ( auto i = 0; i < points.size(); ++i )
331 {
332 viewerCore.setFillColor( Color( 150, 150, 150, 255 ) );
333 viewerCore.addBall( RealPoint( points[ i ][ 0 ],
334 points[ i ][ 1 ],
335 points[ i ][ 2 ] ), 12 );
336 }
337 viewerCore.setLineColor( Color::Green );
338 for ( auto path : paths )
339 {
340 for ( auto i = 1; i < path.size(); i++ )
341 {
342 Point p1 = TC.point( path[ i-1 ] );
343 Point p2 = TC.point( path[ i ] );
344 viewerCore.addLine( p1, p2, 18.0 );
345 }
346 trace.info() << "length=" << TC.length( path ) << std::endl;
347 }
348 viewerCore << MViewer3D::updateDisplay;
349 application.exec();
350 }
351
352 return 0;
353}
DGtal::Color
Structure representing an RGB triple with alpha component.
Definition Color.h:68
DGtal::Color::Green
static const Color Green
Definition Color.h:417
DGtal::KhalimskySpaceND
Aim: This class is a model of CCellularGridSpaceND. It represents the cubical grid as a cell complex,...
Definition KhalimskySpaceND.h:394
DGtal::KhalimskySpaceND::sOrthDir
Dimension sOrthDir(const SCell &s) const
Given a signed surfel [s], returns its orthogonal direction (ie, the coordinate where the surfel is c...
DGtal::KhalimskySpaceND::sCoords
Point sCoords(const SCell &c) const
Return its digital coordinates.
DGtal::KhalimskySpaceND::sDirect
bool sDirect(const SCell &p, Dimension k) const
Return 'true' if the direct orientation of [p] along [k] is in the positive coordinate direction.
DGtal::KhalimskySpaceND::sIncident
SCell sIncident(const SCell &c, Dimension k, bool up) const
Return the forward or backward signed cell incident to [c] along axis [k], depending on [up].
DGtal::Shortcuts::getKSpace
static KSpace getKSpace(const Point &low, const Point &up, Parameters params=parametersKSpace())
Definition Shortcuts.h:332
DGtal::Shortcuts::getSurfelRange
static SurfelRange getSurfelRange(CountedPtr< ::DGtal::DigitalSurface< TDigitalSurfaceContainer > > surface, const Parameters &params=parametersDigitalSurface())
Definition Shortcuts.h:1547
DGtal::Shortcuts::makeDigitalSurface
static CountedPtr< DigitalSurface > makeDigitalSurface(CountedPtr< TPointPredicate > bimage, const KSpace &K, const Parameters &params=parametersDigitalSurface())
Definition Shortcuts.h:1209
DGtal::Shortcuts::defaultParameters
static Parameters defaultParameters()
Definition Shortcuts.h:203
DGtal::Shortcuts::makeBinaryImage
static CountedPtr< BinaryImage > makeBinaryImage(Domain shapeDomain)
Definition Shortcuts.h:561
DGtal::SimpleDistanceColorMap
Aim: simple blue to red colormap for distance information for instance.
Definition SimpleDistanceColorMap.h:65
DGtal::TangencyComputer
Aim: A class that computes tangency to a given digital set. It provides services to compute all the c...
Definition TangencyComputer.h:74
DGtal::Trace::info
std::ostream & info()
DGtal::Viewer3D::show
virtual void show()
Overload QWidget method in order to add a call to updateList() method (to ensure that the lists are w...
surface
CountedPtr< SH3::DigitalSurface > surface
Definition dgtalCalculus-geodesic.cpp:66
reaction
int reaction(void *viewer, DGtal::int32_t name, void *data)
Definition fullConvexityShortestPaths3D.cpp:90
RealPoint
Space::RealPoint RealPoint
Definition fullConvexityShortestPaths3D.cpp:86
Index
SMesh::Index Index
Definition fullConvexitySphereGeodesics.cpp:117
DGtal::Dimension
DGtal::uint32_t Dimension
Definition Common.h:136
DGtal::trace
Trace trace
Definition Common.h:153
DGtal::int32_t
boost::int32_t int32_t
signed 32-bit integer.
Definition BasicTypes.h:72
DGtal::SetMode3D
Modifier class in a Display3D stream. Useful to choose your own mode for a given class....
Definition DrawWithDisplay3DModifier.h:74
Point
MyPointD Point
Definition testClone2.cpp:383
K
KSpace K
Definition testCubicalComplex.cpp:62

References DGtal::Color::Black, DGtal::Shortcuts< TKSpace >::defaultParameters(), DGtal::Shortcuts< TKSpace >::getKSpace(), DGtal::Shortcuts< TKSpace >::getSurfelRange(), DGtal::Color::Green, DGtal::Trace::info(), K, DGtal::Shortcuts< TKSpace >::makeBinaryImage(), DGtal::Shortcuts< TKSpace >::makeDigitalSurface(), reaction(), DGtal::KhalimskySpaceND< dim, TInteger >::sCoords(), DGtal::KhalimskySpaceND< dim, TInteger >::sDirect(), DGtal::Viewer3D< TSpace, TKSpace >::show(), DGtal::KhalimskySpaceND< dim, TInteger >::sIncident(), DGtal::KhalimskySpaceND< dim, TInteger >::sOrthDir(), surface, and DGtal::trace.

◆ reaction()

int reaction ( void * viewer,
DGtal::int32_t name,
void * data )
Examples
geometry/volumes/fullConvexityShortestPaths3D.cpp.

Definition at line 90 of file fullConvexityShortestPaths3D.cpp.

91{
92 DGtal::int32_t* selected = (DGtal::int32_t*) data;
93 *selected = name;
94 std::cout << "Selected surfel=" << *selected << std::endl;
95 return 0;
96}

Referenced by main().