testContourHelper.cpp

Go to the documentation of this file.

 
#include <iostream>
#include "DGtal/base/Common.h"
#include "ConfigTest.h"
#include "DGtal/geometry/helpers/ContourHelper.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/kernel/PointVector.h"
 
using namespace std;
using namespace DGtal;
 
// Functions for testing class ContourHelper.
 
bool testContourHelper()
{
  unsigned int nbok = 0;
  unsigned int nb = 0;
  
  trace.beginBlock ( "Test ContourHelper:  getgetBarycenter and isCounterClockWise" );
  std::vector<Z2i::Point> aContour; 
  aContour.push_back(Z2i::Point(0,0));
  aContour.push_back(Z2i::Point(10,0));
  aContour.push_back(Z2i::Point(10,4));
  aContour.push_back(Z2i::Point(0,4));
  Z2i::RealPoint midPoint = ContourHelper::getBarycenter(aContour);   
 
  std::vector<Z2i::Point> aContour2; 
  aContour2.push_back(Z2i::Point(0,0));
  aContour2.push_back(Z2i::Point(0,10));
  aContour2.push_back(Z2i::Point(5,10));
 
  nbok += midPoint==Z2i::Point(5,2) ? 1 : 0; 
  nb++;
  nbok += ContourHelper::isCounterClockWise(aContour) ? 1 : 0; 
  nb++;
  nbok += ContourHelper::isCounterClockWise(aContour2) ? 0 : 1; 
  nb++;
  trace.endBlock();
 
  trace.beginBlock("Test ContourHelper:  pixels2pixels8c ");
 
  std::vector<Z2i::Point> aContour3;
  aContour3.push_back(Z2i::Point(0,0));
  aContour3.push_back(Z2i::Point(1,0));
  aContour3.push_back(Z2i::Point(1,1));
  aContour3.push_back(Z2i::Point(2,1));
  aContour3.push_back(Z2i::Point(3,1));
  aContour3.push_back(Z2i::Point(4,1));
  aContour3.push_back(Z2i::Point(4,2));
  aContour3.push_back(Z2i::Point(5,2));
  aContour3.push_back(Z2i::Point(6,2));
  std::vector<Z2i::Point> res, res2;
    
  ContourHelper::pixels2pixels8C(aContour3.begin(),
                                 aContour3.end(), std::back_inserter(res));
  trace.info() << "Apply pixels2pixels8c (from std::vector contour)init size : " << res.size() << " (should be: " << aContour3.size() -2  << " ) "<< std::endl;  
  trace.info() << "Third point : " << res[2] << " (should be: " << Z2i::Point(2,1) << " ) "<< std::endl;  
  nb++;
  nbok += res.size() ==  aContour3.size() -2 || Z2i::Point(2,1)==res[2];
 
  typedef GridCurve< KhalimskySpaceND<2>> GridCurve2d ;
  GridCurve2d gc;
  gc.initFromPointsVector(aContour3);
  ContourHelper::pixels2pixels8C(gc.getPointsRange().begin(),
                                 gc.getPointsRange().end(), std::back_inserter(res2));
  trace.info() << "Apply pixels2pixels8c (from GridCurve contour) init size : " << res2.size() << " (should be: " << aContour3.size() -2  << " ) "<< std::endl;  
  trace.info() << "Third point : " << res2[2] << " (should be: " << Z2i::Point(2,1) << " ) "<< std::endl;  
  nb++;
  nbok += res2.size() ==  aContour3.size() -2 || Z2i::Point(2,1)==res2[2];
  
  trace.info() << "(" << nbok << "/" << nb << ") "<< std::endl;
  
  for(auto p: res)
    trace.info() << p << std::endl;
  return nbok == nb;
}
 
// Standard services - public :
 
int main( int argc, char** argv )
{
  trace.beginBlock ( "Testing class ContourHelper" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << endl;
 
  bool res = testContourHelper(); // && ... other tests
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //