DGtal 1.3.0
Loading...
Searching...
No Matches
volToOFF.cpp
Go to the documentation of this file.
1
55#include <iostream>
56#include "DGtal/base/Common.h"
57#include "DGtal/io/readers/VolReader.h"
58#include "DGtal/helpers/StdDefs.h"
59#include "DGtal/topology/helpers/Surfaces.h"
60#include "DGtal/topology/DigitalSurface.h"
61#include "DGtal/topology/SetOfSurfels.h"
62#include "DGtal/images/ImageSelector.h"
63#include "DGtal/images/imagesSetsUtils/SetFromImage.h"
65
67
68using namespace std;
69using namespace DGtal;
70using namespace Z3i;
71
73
74
75void usage( int, char** argv )
76{
77 std::cerr << "Usage: " << argv[ 0 ] << " <fileName.vol> <minT> <maxT> <int=0|ext=1>" << std::endl;
78 std::cerr << "\t - displays the boundary of the shape stored in vol file <fileName.vol> as an OFF geomview surface file. It is a kind of marching-cube surface, defined by duality with respect to the digital surface." << std::endl;
79 std::cerr << "\t - voxel v belongs to the shape iff its value I(v) follows minT <= I(v) <= maxT." << std::endl;
80 std::cerr << "\t - 0: interior adjacency, 1: exterior adjacency." << std::endl;
81}
82
83int main( int argc, char** argv )
84{
85 if ( argc < 5 )
86 {
87 usage( argc, argv );
88 return 1;
89 }
90 std::string inputFilename = argv[ 1 ];
91 unsigned int minThreshold = atoi( argv[ 2 ] );
92 unsigned int maxThreshold = atoi( argv[ 3 ] );
93 bool intAdjacency = atoi( argv[ 4 ] ) == 0;
94
96 trace.beginBlock( "Reading vol file into an image." );
98 Image image = VolReader<Image>::importVol(inputFilename);
99 DigitalSet set3d (image.domain());
101 minThreshold, maxThreshold);
102 trace.endBlock();
104
106 // Construct the Khalimsky space from the image domain
107 KSpace K;
108 bool space_ok = K.init( image.domain().lowerBound(),
109 image.domain().upperBound(), true );
110 if (!space_ok)
111 {
112 trace.error() << "Error in the Khamisky space construction."<<std::endl;
113 return 2;
114 }
116
118 typedef SurfelAdjacency<KSpace::dimension> MySurfelAdjacency;
119 MySurfelAdjacency surfAdj( intAdjacency ); // interior in all directions.
121
123 trace.beginBlock( "Extracting boundary by scanning the space. " );
125 typedef SetOfSurfels< KSpace, SurfelSet > MySetOfSurfels;
127
128 MySetOfSurfels theSetOfSurfels( K, surfAdj );
129 Surfaces<KSpace>::sMakeBoundary( theSetOfSurfels.surfelSet(),
130 K, set3d,
131 image.domain().lowerBound(),
132 image.domain().upperBound() );
133 MyDigitalSurface digSurf( theSetOfSurfels );
134 trace.info() << "Digital surface has " << digSurf.size() << " surfels."
135 << std::endl;
136 trace.endBlock();
138
140 trace.beginBlock( "Making OFF surface <marching-cube.off>. " );
141 ofstream out( "marching-cube.off" );
142 if ( out.good() )
143 digSurf.exportSurfaceAs3DOFF( out );
144 out.close();
145 trace.endBlock();
147
148 return 0;
149}
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...
void exportSurfaceAs3DOFF(std::ostream &out) const
Aim: implements association bewteen points lying in a digital domain and values.
Definition: Image.h:70
Aim: This class is a model of CCellularGridSpaceND. It represents the cubical grid as a cell complex,...
std::set< SCell > SurfelSet
Preferred type for defining a set of surfels (always signed cells).
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 connected surfels....
Definition: SetOfSurfels.h:74
static void sMakeBoundary(SCellSet &aBoundary, const KSpace &aKSpace, const PointPredicate &pp, const Point &aLowerBound, const Point &aUpperBound)
Aim: Represent adjacencies between surfel elements, telling if it follows an interior to exterior ord...
void beginBlock(const std::string &keyword="")
std::ostream & error()
std::ostream & info()
double endBlock()
DigitalSurface< MyDigitalSurfaceContainer > MyDigitalSurface
MyDigitalSurface::SurfelSet SurfelSet
DGtal is the top-level namespace which contains all DGtal functions and types.
Trace trace
Definition: Common.h:154
STL namespace.
ImageContainerBySTLVector< Domain, Value > Type
Definition: ImageSelector.h:78
Aim: Define utilities to convert a digital set into an image.
Definition: SetFromImage.h:64
static ImageContainer importVol(const std::string &filename, const Functor &aFunctor=Functor())
int main()
Definition: testBits.cpp:56
KSpace K