viewer3D-7-stdplane.cpp File Reference

#include <cstdlib>
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/io/viewers/PolyscopeViewer.h"
#include "DGtal/geometry/surfaces/COBAGenericStandardPlaneComputer.h"

Include dependency graph for viewer3D-7-stdplane.cpp:

Go to the source code of this file.

Functions
template<typename Viewer3D, typename Domain, typename Predicate>
void	displayPredicate (Viewer3D &viewer, const Domain &domain, const Predicate &pred)
template<typename Viewer3D, typename InputIterator>
void	displayRange (Viewer3D &viewer, InputIterator it, InputIterator itE)
template<typename Domain>
std::vector< typename Domain::Point >	pointsInStandardPlane (const Domain &domain, typename Domain::Integer a, typename Domain::Integer b, typename Domain::Integer c, typename Domain::Integer mu)
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

Jacques-Olivier Lachaud (jacqu.nosp@m.es-o.nosp@m.livie.nosp@m.r.la.nosp@m.chaud.nosp@m.@uni.nosp@m.v-sav.nosp@m.oie..nosp@m.fr ) Laboratory of Mathematics (CNRS, UMR 5127), University of Savoie, France

Date

2012/03/05

Functions for testing class COBANaivePlaneComputer.

This file is part of the DGtal library.

Definition in file viewer3D-7-stdplane.cpp.

Function Documentation

displayPredicate()

template<typename Viewer3D, typename Domain, typename Predicate>

void displayPredicate	(	Viewer3D &	viewer,
		const Domain &	domain,
		const Predicate &	pred )

Definition at line 72 of file viewer3D-7-stdplane.cpp.

{
  for ( typename Domain::ConstIterator itB = domain.begin(), itE = domain.end();
        itB != itE; ++itB )
    {
      if ( pred( *itB ) )
        viewer << *itB;
    }
}

References domain.

Referenced by main().

displayRange()

template<typename Viewer3D, typename InputIterator>

void displayRange	(	Viewer3D &	viewer,
		InputIterator	it,
		InputIterator	itE )

Definition at line 87 of file viewer3D-7-stdplane.cpp.

{
  for ( ; it != itE; ++it )
    viewer << *it;
}

Referenced by main().

main()

int main	(	int	argc,
		char **	argv )

Definition at line 122 of file viewer3D-7-stdplane.cpp.

{
  using namespace Z3i;
  unsigned int nb = 0;
  unsigned int nbok = 0;
 
  unsigned int diameter = argc > 1 ? atoi( argv[ 1 ] ) : 10;
  int a = argc > 2 ? atoi( argv[ 2 ] ) : 2;
  int b = argc > 3 ? atoi( argv[ 3 ] ) : 3;
  int c = argc > 4 ? atoi( argv[ 4 ] ) : 5;
  int mu = argc > 5 ? atoi( argv[ 5 ] ) : 0;
  trace.beginBlock ( "Testing class COBAGenericStandardPlaneComputer" );
  trace.info() << "Recognizing plane "
               << mu << " <= " << a << " * x + "
               << b << " * y + " << c << " * z < "
               << (mu+abs(a)+abs(b)+abs(c)) << std::endl;
  Domain domain1( Point( -diameter, -diameter, -diameter ),
                  Point(  diameter,  diameter,  diameter ) );
 
  typedef COBAGenericStandardPlaneComputer<Z3, BigInteger> PlaneComputer;
  typedef PlaneComputer::Primitive Primitive;
  PlaneComputer plane;
  plane.init( 2*diameter, 1, 1 );
 
  std::vector<Point> recognized = pointsInStandardPlane( domain1,
                                                         a, b , c, mu );
  ++nb, nbok += plane.extend( recognized.begin(), recognized.end() ) ? 1 : 0;
  trace.info() << "(" << nbok << "/" << nb
               << ") All points are recognized." << std::endl;
  trace.info() << " - Plane=" << plane
               << std::endl;
  Primitive strip = plane.primitive();
  trace.info() << "strip=" << strip
               << " axis=" << strip.mainAxis()
               << " axiswidth=" << strip.axisWidth()
               << " diag=" << strip.mainDiagonal()
               << " diagwidth=" << strip.diagonalWidth()
               << std::endl;
  ++nb, nbok += ( strip.diagonalWidth() < sqrt(3.0) ) ? 1 : 0;
  trace.info() << "(" << nbok << "/" << nb
               << ") Diagonal width < sqrt(3)." << std::endl;
  trace.emphase() << ( (nb == nbok) ? "Passed." : "Error." ) << endl;
  trace.endBlock();
 
  PolyscopeViewer viewer;
  Color red( 255, 0, 0 );
  Color green( 0, 255, 0 );
  Color grey( 200, 200, 200 );
  Domain domain2( Point( -2*diameter, -2*diameter, -2*diameter ),
                  Point(  2*diameter,  2*diameter,  2*diameter ) );
  viewer << red;
  for ( std::vector<Point>::const_iterator it = recognized.begin(),
          itE = recognized.end(); it != itE; ++it )
    if ( ! strip( *it ) ) viewer << *it;
 
  viewer << green; 
  displayRange( viewer, plane.begin(), plane.end() );
 
  viewer << grey;
  displayPredicate( viewer, domain2, strip );
 
  trace.info() << "- Points in green have been recognized as belonging to this standard plane." << std::endl;
  trace.info() << "- Points in grey belongs also to the parallel strip of the recognized standard plane." << std::endl;
  trace.info() << "- Points in red belongs to the parallel strip of the recognized standard plane but not to the input standard plane: NONE should be red." << std::endl;
 
  viewer.show();
  return 0;
}

References displayPredicate(), displayRange(), DGtal::COBAGenericStandardPlaneComputer< TSpace, TInternalInteger >::init(), pointsInStandardPlane(), DGtal::PolyscopeViewer< Space, KSpace >::show(), and DGtal::trace.

pointsInStandardPlane()

template<typename Domain>

std::vector< typename Domain::Point > pointsInStandardPlane	(	const Domain &	domain,
		typename Domain::Integer	a,
		typename Domain::Integer	b,
		typename Domain::Integer	c,
		typename Domain::Integer	mu )

Definition at line 98 of file viewer3D-7-stdplane.cpp.

{
  typedef typename Domain::Integer Integer;
  typedef typename Domain::Point Point;
  typedef typename Domain::ConstIterator ConstIterator;
  std::vector<Point> pts;
  Integer mup = mu + abs(a) + abs(b) + abs(c);
  for ( ConstIterator it = domain.begin(), itE = domain.end();
        it != itE; ++it )
    {
      Point p = *it;
      Integer r = a * p[ 0 ] + b * p[ 1 ] + c * p[ 2 ];
      if ( ( mu <= r ) && ( r < mup ) )
        pts.push_back( p );
    }
  return pts;
}

References domain.

Referenced by main().