testWindingNumbersShape.cpp

Go to the documentation of this file.

#include <iostream>
#include "DGtal/base/Common.h"
#include "ConfigTest.h"
#include "DGtalCatch.h"
#include "DGtal/helpers/StdDefs.h"
 
#include <DGtal/shapes/WindingNumbersShape.h>
#include <DGtal/shapes/GaussDigitizer.h>

 
using namespace std;
using namespace DGtal;
using namespace Z3i;

// Functions for testing class WindingNumbersShape.
 
TEST_CASE( "Testing WindingNumbersShape" )
{
  using WNShape = WindingNumbersShape<Z3i::Space>;
 
  SECTION("Checking concept")
    {
      BOOST_CONCEPT_ASSERT( (DGtal::concepts::CEuclideanOrientedShape<WNShape>) );
    }
 
  SECTION("Simple test of the API")
  {
    Eigen::MatrixXd points(4,3);
    points << 0,0,0,
              0,1,0,
              1,0,0,
              1,1,1;
    Eigen::MatrixXd normals(4,3);
    normals << 0,0,-1,
               0,0,-1,
               0,0,-1,
               0,0,1;
 
    WNShape wnshape(points,normals);
 
    RealPoint p(-2.0,-2.0,-2.0);
    DGtal::Orientation ori = wnshape.orientation(p);
    REQUIRE( ori == DGtal::OUTSIDE);
 
    p = RealPoint(2.0,2.0,2.0);
    ori = wnshape.orientation(p);
    REQUIRE( ori == DGtal::OUTSIDE);
 
    RealPoint q= RealPoint(.2,.2,.2);
    auto ori2 = wnshape.orientation(q);
    REQUIRE( ori2 == DGtal::INSIDE);
  }
 
 SECTION("Tesing with the GaussDigitizer")
 {
     Eigen::MatrixXd points(4,3);
     points << 0,0,0,
               0,1,0,
               1,0,0,
               1,1,1;
     Eigen::MatrixXd normals(4,3);
     normals << 0,0,-1,
                0,0,-1,
                0,0,-1,
                0,0,1;
 
     WNShape wnshape(points,normals);
     GaussDigitizer<Z3i::Space, WNShape> gauss;
     gauss.attach(wnshape);
     gauss.init(Z3i::RealPoint(0,0,0),Z3i::RealPoint(1.5,1.5,1.5), 0.5);
     auto cpt=0;
     for(auto p: gauss.getDomain())
         if (gauss.orientation(p) == DGtal::INSIDE)
             ++cpt;
     REQUIRE( cpt == 8);
 }
};