io/viewers/viewer3D-7-planes.cpp

This example shows the recognition of a simple digital plane with the COBA algorithm (class COBANaivePlaneComputer). Green points belong to the naive plane. There is no naive plane containing the green points and any one of the red points.

See also

How to recognize a plane ?

$ ./examples/io/viewers/viewer3D-7-planes

Piece of naive plane containing the four green points. There is no naive plane containing also any one of the red points.

#include <cstdlib>
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/io/viewers/PolyscopeViewer.h"
#include "DGtal/geometry/surfaces/COBANaivePlaneComputer.h"
 
using namespace std;
using namespace DGtal;

// Standard services - public :
template <typename Viewer3D, typename Domain, typename Predicate>
void
displayPredicate( Viewer3D & viewer,
                  const Domain & domain, const Predicate & pred )
{
  for ( typename Domain::ConstIterator itB = domain.begin(), itE = domain.end();
        itB != itE; ++itB )
    {
      if ( pred( *itB ) )
        viewer << *itB;
    }
}
 
int main( int argc, char** argv )
{
  using namespace Z3i;
 
  trace.beginBlock ( "Testing class COBANaivePlaneComputer" );
 
  unsigned int nbok = 0;
  unsigned int nb = 0;
 
  typedef COBANaivePlaneComputer<Z3, BigInteger> PlaneComputer;
  typedef PlaneComputer::Primitive Primitive;
  PlaneComputer plane;
 
  plane.init( 2, 100, 1, 1 );
  Point pt0( 0, 0, 0 );
  bool pt0_inside = plane.extend( pt0 );
  trace.info() << "(" << nbok << "/" << nb << ") Plane=" << plane
               << std::endl;
  Point pt1( 8, 1, 3 );
  bool pt1_inside = plane.extend( pt1 );
  ++nb, nbok += pt1_inside == true ? 1 : 0;
  trace.info() << "(" << nbok << "/" << nb << ") add " << pt1
               << " Plane=" << plane << std::endl;
  Point pt2( 2, 7, 1 );
  bool pt2_inside = plane.extend( pt2 );
  ++nb, nbok += pt2_inside == true ? 1 : 0;
  trace.info() << "(" << nbok << "/" << nb << ") add " << pt2
               << " Plane=" << plane << std::endl;
 
  Point pt3( 0, 5, 12 );
  bool pt3_inside = plane.extend( pt3 );
  ++nb, nbok += pt3_inside == false ? 1 : 0;
  trace.info() << "(" << nbok << "/" << nb << ") add " << pt3
               << " Plane=" << plane << std::endl;
 
  Point pt4( -5, -5, 10 );
  bool pt4_inside = plane.extend( pt4 );
  ++nb, nbok += pt4_inside == false ? 1 : 0;
  trace.info() << "(" << nbok << "/" << nb << ") add " << pt4
               << " Plane=" << plane << std::endl;
 
  Point pt5 = pt0 + pt1 + pt2 + Point( 0, 0, 1 );
  bool pt5_inside = plane.extend( pt5 );
  ++nb, nbok += pt5_inside == true ? 1 : 0;
  trace.info() << "(" << nbok << "/" << nb << ") add " << pt5
               << " Plane=" << plane << std::endl;
 
  Point pt6 = Point( 1, 0, 1 );
  bool pt6_inside = plane.extend( pt6 );
  ++nb, nbok += pt6_inside == true ? 1 : 0;
  trace.info() << "(" << nbok << "/" << nb << ") add " << pt5
               << " 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.axisWidth() < 1.0 ? 1 : 0;
  trace.info() << "(" << nbok << "/" << nb << ") axiswidth < 1 "
               << std::endl;
  ++nb, nbok += strip.diagonalWidth() < sqrt(3.0) ? 1 : 0;
  trace.info() << "(" << nbok << "/" << nb << ") axiswidth < sqrt(3) "
               << std::endl;
  trace.emphase() << ( nbok == nb ? "Passed." : "Error." ) << endl;
  trace.endBlock();
 
  PolyscopeViewer viewer;
  Color red( 255, 0, 0 );
  Color green( 0, 255, 0 );
  Color grey( 200, 200, 200 );
  Domain domain( Point( -5, -5, -5 ), Point( 12, 12, 12 ) );
  viewer << grey;
  viewer << ( pt0_inside ? green : red ) << pt0;
  viewer << ( pt1_inside ? green : red ) << pt1;
  viewer << ( pt2_inside ? green : red ) << pt2;
  viewer << ( pt3_inside ? green : red ) << pt3;
  viewer << ( pt4_inside ? green : red ) << pt4;
  viewer << ( pt5_inside ? green : red ) << pt5;
  viewer << ( pt6_inside ? green : red ) << pt6;
  viewer << grey;
  displayPredicate( viewer, domain, strip );
 
  viewer.show();
  return 0;
}
//                                                                           //