#include <iostream>
#include <fstream>
#include <algorithm>
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/images/ImageContainerBySTLVector.h"
#include "DGtal/io/readers/VolReader.h"
#include "DGtal/io/writers/VolWriter.h"
#include "DGtal/images/SimpleThresholdForegroundPredicate.h"
#include "DGtal/geometry/volumes/distance/DistanceTransformation.h"
#include "DGtal/shapes/implicit/ImplicitBall.h"
#include "DGtal/base/BasicFunctors.h"
#include "DGtal/io/viewers/Viewer3D.h"
#include "DGtal/io/colormaps/GradientColorMap.h"
int main(int argc, char ** argv)
{
QApplication application(argc,argv);
itend = image.
domain().end(); it != itend; ++it)
if (image(*it) != 0)
viewer << *it;
viewer << DGtal::Viewer3D<>::updateDisplay;
trace.
info() <<
"viewer launched..."<<std::endl;
bool res = application.exec();
Predicate binaryshape(image, 0);
DT distancemap(image.
domain(), binaryshape, l2);
QApplication application2(argc,argv);
DT::Value maxDT = (*std::max_element(distancemap.constRange().begin(),
distancemap.constRange().end()));
trace.
info() <<
"we display the dt map"<<std::endl;
int cpt=0;
for(DT::Domain::ConstIterator it = distancemap.domain().begin(),
itend = distancemap.domain().end(); it != itend;
++it)
if (distancemap(*it) > 0)
{
DT::Value val= distancemap( *it );
cpt++;
}
trace.
info() <<
"Got "<<cpt<<
" points."<<std::endl;
viewer2 << DGtal::Viewer3D<>::updateDisplay;
trace.
info() <<
"viewer2 launched..."<<std::endl;
res = res && application2.exec();
Image imageGranulo ( image.
domain() );
it != itend; ++it)
*it = 0;
trace.
info() <<
"Computing the granulometry"<<std::endl;
cpt=0;
itend = imageGranulo.domain().end(); it != itend; ++it)
{
if (distancemap(*it) > 0)
{
unsigned int radius = (unsigned int)distancemap(*it);
itball != itendball; itball++)
if (imageGranulo.domain().isInside(*itball) &&
( ball(*itball) > 0) &&
(imageGranulo(*itball) < radius))
imageGranulo.setValue(*itball, radius);
cpt++;
}
}
trace.
info() <<
"Granulometry ok nbBalls="<<cpt<< std::endl;
QApplication application3(argc,argv);
Image::Value maxG = (*std::max_element(imageGranulo.constRange().begin(),
imageGranulo.constRange().end()));
viewer3 << DGtal::ClippingPlane(1,0,0,-10.1);
cpt=0;
itend = imageGranulo.domain().end(); it != itend;
++it)
if (imageGranulo(*it) > 0)
{
cpt++;
}
trace.
info() <<
"Got "<<cpt<<
" points."<<std::endl;
viewer3 << DGtal::Viewer3D<>::updateDisplay;
trace.
info() <<
"viewer3 launched..."<<std::endl;
return application3.exec();
}