volDTGranulo.cpp File Reference

An example of generating a grid curve from a parametric shape and estimating its length. More...

#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/PolyscopeViewer.h"
#include "DGtal/io/colormaps/GradientColorMap.h"
#include <boost/algorithm/minmax_element.hpp>

Include dependency graph for volDTGranulo.cpp:

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Detailed Description

An example of generating a grid curve from a parametric shape and estimating its length.

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

Tristan Roussillon (trist.nosp@m.an.r.nosp@m.oussi.nosp@m.llon.nosp@m.@liri.nosp@m.s.cn.nosp@m.rs.fr)

Date

2010/10/17

This file is part of the DGtal library.

Definition in file volDTGranulo.cpp.

Function Documentation

main()

int main	(	int	argc,
		char **	argv )

Definition at line 57 of file volDTGranulo.cpp.

{
 
  //Loarding the image
  typedef ImageContainerBySTLVector<Z3i::Domain, unsigned int> Image;
  Image image = VolReader<Image>::importVol(argv[1]);
  trace.info() << image << std::endl;
 
  //Viewer
  PolyscopeViewer<> viewer;
 
  for(Image::Domain::ConstIterator it = image.domain().begin(),
        itend = image.domain().end(); it != itend; ++it)
    if (image(*it) != 0)
      viewer << *it;
 
  trace.info() << "viewer launched..."<<std::endl;
  viewer.show();
 
 
  //DT
  Z3i::L2Metric l2;
  typedef functors::SimpleThresholdForegroundPredicate<Image> Predicate;
  Predicate binaryshape(image, 0);
 
  typedef DistanceTransformation< Z3i::Space, Predicate, Z3i::L2Metric> DT;
 
  DT distancemap(image.domain(), binaryshape, l2);
 
  //Viewer
  PolyscopeViewer<> viewer2;
  DT::Value maxDT = (*boost::first_max_element(distancemap.constRange().begin(),
                                       distancemap.constRange().end()));
 trace.info() << "we display the dt map"<<std::endl;
  int cpt=0;
  viewer2 << DGtal::ClippingPlane(1,0,0,-10.1);
 
  for(DT::Domain::ConstIterator it = distancemap.domain().begin(),
  itend = distancemap.domain().end(); it != itend;
      ++it)
    if (distancemap(*it) > 0)
      {
        viewer2 << DGtal::WithQuantity(*it, "value", distancemap(*it));
        cpt++;
      }
  trace.info() << "Got "<<cpt<<" points."<<std::endl;
  trace.info() << "viewer2  launched..."<<std::endl;
  viewer2.show();
 
  //Granulo
  Image imageGranulo ( image.domain() );
  for(Image::Range::Iterator it = imageGranulo.range().begin(), itend= imageGranulo.range().end();
      it != itend; ++it)
    *it = 0;
 
 
  trace.info() << "Computing the granulometry"<<std::endl;
  cpt=0;
  for(Image::Domain::ConstIterator it = imageGranulo.domain().begin(),
        itend = imageGranulo.domain().end(); it != itend; ++it)
    {
      if (distancemap(*it) > 0)
        {
          //Construct the sphere with radius from distancemap(*it)
          Z3i::RealPoint center = *it;
          unsigned int radius = (unsigned int)distancemap(*it);
          ImplicitBall<Z3i::Space> ball(center,radius);
          Z3i::Point low(ball.getLowerBound() - Z3i::RealPoint::diagonal(1.0), functors::Floor<>());
          Z3i::Point up (ball.getUpperBound() + Z3i::RealPoint::diagonal(1.0), functors::Ceil<>());
          Z3i::Domain dom(low,up);
 
          for(Z3i::Domain::ConstIterator itball = dom.begin(), itendball= dom.end();
              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;
  VolWriter<Image, functors::Cast<unsigned char> >::exportVol("granulo.vol", imageGranulo);
  trace.info() << "Save OK"<< std::endl;
 
 
  //Viewer
  PolyscopeViewer<> viewer3;
  Image::Value maxG = (*boost::first_max_element(imageGranulo.constRange().begin(),
                                         imageGranulo.constRange().end()));
 
  viewer3 << DGtal::ClippingPlane(1,0,0,-10.1);
  cpt=0;
  for(Image::Domain::ConstIterator it = imageGranulo.domain().begin(),
  itend = imageGranulo.domain().end(); it != itend;
      ++it)
    if (imageGranulo(*it) > 0)
      {
        viewer3 << DGtal::WithQuantity(*it, "value", imageGranulo(*it));
        cpt++;
      }
  trace.info() << "Got "<<cpt<<" points."<<std::endl;
  trace.info() << "viewer3  launched..."<<std::endl;
  viewer3.show();
 return 0;
 
}

References DGtal::HyperRectDomain< TSpace >::begin(), DGtal::SimpleRandomAccessConstRangeFromPoint< TConstIterator, DistanceFunctor >::begin(), DGtal::SimpleRandomAccessRangeFromPoint< TConstIterator, TIterator, DistanceFunctor >::begin(), DGtal::Image< TImageContainer >::constRange(), DGtal::PointVector< dim, double >::diagonal(), DGtal::Image< TImageContainer >::domain(), DGtal::HyperRectDomain< TSpace >::end(), DGtal::SimpleRandomAccessConstRangeFromPoint< TConstIterator, DistanceFunctor >::end(), DGtal::SimpleRandomAccessRangeFromPoint< TConstIterator, TIterator, DistanceFunctor >::end(), DGtal::ImplicitBall< TSpace >::getLowerBound(), DGtal::ImplicitBall< TSpace >::getUpperBound(), image(), DGtal::VolReader< TImageContainer, TFunctor >::importVol(), DGtal::HyperRectDomain< TSpace >::isInside(), DGtal::Image< TImageContainer >::range(), DGtal::Image< TImageContainer >::setValue(), DGtal::PolyscopeViewer< Space, KSpace >::show(), and DGtal::trace.