volDTGranulo-template.cpp

 
#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"
 
using namespace DGtal;
 
int main(int argc, char ** argv)
{
 
  //Load a vol image and store it into a ImageContainerBySTLVector
  //
  // ...
  //  
 
 
  //Visualization 2
  //
 
  //Create a new viewer instance
  QApplication application(argc,argv);
  Viewer3D<> viewer;
  viewer.show();
 
  // - iterate over the loaded image domain and send strictly positive
  //values to the viewer
  //
  // Do not forget a  viewer2 << DGtal::Viewer3D<>::updateDisplay;
  // when you're done.
 
 
  //Distance transformation
  
  //Create a Point Predicate from the input image and a thresholder
  // (see SimpleThresholdForegroundPredicate)
 
  // ...
 
 
  // Compute the distance map for the L2 metric
 
  // ...
 
 
 
  //Visualization 2
  //
 
  // Create a new instance of Viewer3D
 
  //Create the colormap from the maximal distance value
  //and a GradientColorMap
  DT::Value maxDT = (*boost::first_max_element(distancemap.constRange().begin(), 
                                       distancemap.constRange().end()));
  GradientColorMap<DT::Value> gradient( 0, maxDT);
  gradient.addColor(DGtal::Color::Blue);
  gradient.addColor(DGtal::Color::Green);
  gradient.addColor(DGtal::Color::Yellow);
  gradient.addColor(DGtal::Color::Red);
 
  
 
  //Granulometry
 
 
  //Create an empty image with the same size as the input one
  
  //...
 
 
  //Compute the granulometry function
  
  //...
 
  
  //Export the granulometry map to a vol
  VolWriter<Image, CastFunctor<unsigned char> >::exportVol("granulo.vol", imageGranulo);     
  
 
 
  //Visualization 3
  //
  //  Create a visualization of the granulometry function
  
  return 0;
}