examplePlaneProbingParallelepipedEstimator.cpp

Go to the documentation of this file.

 
#include <iostream>
#include "ConfigExamples.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/base/Common.h"
#include "DGtal/geometry/surfaces/DigitalPlanePredicate.h"
#include "DGtal/geometry/surfaces/estimation/PlaneProbingParallelepipedEstimator.h"
 
using namespace std;
using namespace DGtal;
 
using DigitalPlane = DigitalPlanePredicate<Z3i::Space>;
using Point        = DigitalPlane::Vector;
using Integer      = Point::Coordinate;
using Vector       = DigitalPlane::Point;
using Estimator    = PlaneProbingParallelepipedEstimator<DigitalPlane, ProbingMode::R1>;
 
 
int main(int argc, const char** argv)
{
  const int height = (argc > 1) ? std::atoi(argv[1]) : 2;
  std::cout << "height=" << height << std::endl;
 
  Vector n(2, 6, 15);
 
  DigitalPlane plane(n, -height, n.norm1());
  Point o(0, 0, 0);
  Integer bound = n.norm1() + 1;
  std::array<Point, 3> m = { Point(1, 0, 0), Point(0, 1, 0), Point(0, 0, 1) };
  Estimator estimator(o, m, plane, bound);
 
  int it = 0;
  while (estimator.advance()) {
      it++;
 
      std::cout << "state=" << estimator.getState() << ", isInReverseState=" << estimator.isInReverseState() << std::endl;
 
      std::clog << "it = " << it << " "
          << estimator.m(0) << " " << estimator.m(1) << " " << estimator.m(2) << " "
          << estimator.getNormal() << std::endl;
  }
 
  ASSERT(estimator.getNormal() == n);
 
  return 0;
}
//                                                                           //