exampleGridCurve3d.cpp

#include <iostream>

#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "ConfigExamples.h"

#include "DGtal/geometry/curves/GridCurve.h"

  #ifdef WITH_VISU3D_QGLVIEWER
#include <QtGui/qapplication.h>
#include "DGtal/io/viewers/Viewer3D.h"
#include "DGtal/io/DrawWithDisplay3DModifier.h"
   #endif


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

#include <boost/program_options/options_description.hpp>
#include <boost/program_options/parsers.hpp>
#include <boost/program_options/variables_map.hpp>

namespace po = boost::program_options;

int main( int argc, char** argv )
{
  // parse command line 
  po::options_description general_opt("Allowed options are");
  general_opt.add_options()
    ("help,h", "display this message")
    ("range,r",  po::value<string>()->default_value("gridcurve"), 
     " Either <gridcurve> (default), <scells>, <points>, <midpoints>, <arrows> " ); 
  
  po::variables_map vm;
  po::store(po::parse_command_line(argc, argv, general_opt), vm);  
  po::notify(vm);    
  if(vm.count("help"))
    {
      trace.info()<< "exampleGridCurve3d" << std::endl
                  << "Basic usage: "<<std::endl
                  << argv[0] << " " << std::endl
                  << general_opt << "\n";
      return 0;
    }
  
  //Parse options
  string type = vm["range"].as<string>(); 

  //curve
  string sinus = examplesPath + "samples/sinus.dat";
  
  // domain
  Point lowerBound = Point::diagonal( -100 );
  Point upperBound = Point::diagonal( 100 ); 

  K3 ks; ks.init( lowerBound, upperBound, true ); 
  GridCurve<K3> gc( ks ); 
  
  fstream inputStream;
  inputStream.open (sinus.c_str(), ios::in);

  gc.initFromVectorStream(inputStream);

  inputStream.close();

  bool flag = false; 
  #ifdef WITH_VISU3D_QGLVIEWER
  QApplication application(argc,argv);
  Viewer3D viewer;
  viewer.show();

  if (type == "gridcurve")
    {
      viewer  << gc;
    }
  else if (type == "scells")
    {
      viewer << gc.getSCellsRange(); 
    }
  else if (type == "points")
    {
      viewer << gc.getPointsRange(); 
    }
  else if (type == "midpoints")
    {
      viewer << gc.getMidPointsRange(); 
    }
  else if (type == "arrows")
    {
      viewer << gc.getArrowsRange(); 
    }
  else
    {
      trace.info() << "Display type not known. Use option -h" << std::endl; 
    }
  viewer << Viewer3D::updateDisplay;
  flag = application.exec();
   #endif
  
  return flag;
}
//                                                                           //