homotopicThinning3D.cpp File Reference

#include <iostream>
#include <queue>
#include "DGtal/base/Common.h"
#include "DGtal/io/viewers/PolyscopeViewer.h"
#include "DGtal/io/Color.h"
#include "DGtal/shapes/Shapes.h"
#include "DGtal/helpers/StdDefs.h"

Include dependency graph for homotopicThinning3D.cpp:

Go to the source code of this file.

Functions
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

Bertrand Kerautret (kerau.nosp@m.tre@.nosp@m.loria.nosp@m..fr ) LORIA (CNRS, UMR 7503), University of Nancy, France

Date

2011/01/04

An example file named Viewer.

This file is part of the DGtal library.

Definition in file homotopicThinning3D.cpp.

Function Documentation

main()

int main	(	int	argc,
		char **	argv )

Definition at line 63 of file homotopicThinning3D.cpp.

{
 
  trace.beginBlock ( "Example simple example of 3DViewer" );
  
  PolyscopeViewer<> viewer;
  
  // Domain cretation from two bounding points.
  Point c( 0, 0, 0 );
  Point p1( -50, -50, -50 );
  Point p2( 50, 50, 50 );
  Domain domain( p1, p2 );
  
  trace.warning() << "Constructing a ring DigitalSet  ... ";
  DigitalSet shape_set( domain );
  for (Domain::ConstIterator it = domain.begin(); it != domain.end(); ++it )
    {
      if ( ((*it - c ).norm() <= 25) && ((*it - c ).norm() >= 18)
     && ( (((*it)[0] <= 3)&& ((*it)[0] >= -3))|| (((*it)[1] <= 3)&& ((*it)[1] >= -3)))){
  shape_set.insertNew( *it );
      }
    }
  trace.warning() << "  [Done]" << std::endl;
  
  trace.beginBlock ( "Thinning" );
  Object18_6 shape( dt18_6,  shape_set );
  int nb_simple=0; 
  DigitalSet::Iterator it, itE;
  do 
    {
      DigitalSet & S = shape.pointSet();
      std::queue<DigitalSet::Iterator> Q;
      it = S.begin(); 
      itE = S.end();
#ifdef DGTAL_WITH_OPENMP
      std::vector<DigitalSet::Iterator> v( S.size() );
      std::vector<uint8_t> b( v.size() );
      for ( size_t i = 0; it != itE; ++it, ++i )
        v[ i ] = it;
#pragma omp parallel for schedule(dynamic)
      for ( size_t i = 0; i < v.size(); ++i )
        b[ i ] = shape.isSimple( *(v[ i ]) );
 
      for ( size_t i = 0; i < v.size(); ++i )
        if ( b[ i ] ) Q.push( v[ i ] ); 
#else
      for ( ; it != itE; ++it )
        if ( shape.isSimple( *it ) )
          Q.push( it );
#endif
      nb_simple = 0;
      while ( ! Q.empty() )
        {
          DigitalSet::Iterator itt = Q.front();
          Q.pop();
          if ( shape.isSimple( *itt ) )
            {
              S.erase( *itt );
              ++nb_simple;
            }
        }
    }
  while ( nb_simple != 0 );
  DigitalSet & S = shape.pointSet();
  trace.endBlock();
 
  // Display by using two different list to manage OpenGL transparency.
 
  viewer << Color(25,25,255, 255);
  viewer << S ; 
 
  viewer << Color(250, 0,0, 25);
  viewer << shape_set;
 
  
  trace.endBlock();
  viewer.show();  
  return 0;
 
}

References DGtal::DigitalSetByAssociativeContainer< TDomain, TContainer >::begin(), domain, DGtal::Z3i::dt18_6, DGtal::DigitalSetByAssociativeContainer< TDomain, TContainer >::end(), DGtal::DigitalSetByAssociativeContainer< TDomain, TContainer >::erase(), DGtal::DigitalSetByAssociativeContainer< TDomain, TContainer >::insertNew(), DGtal::Object< TDigitalTopology, TDigitalSet >::isSimple(), DGtal::Object< TDigitalTopology, TDigitalSet >::pointSet(), DGtal::PolyscopeViewer< Space, KSpace >::show(), DGtal::DigitalSetByAssociativeContainer< TDomain, TContainer >::size(), and DGtal::trace.