freeman2sdp.cpp

 
#include <iostream>
 
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
 
//image
#include "DGtal/io/readers/PointListReader.h"
 
 
//contour
#include "DGtal/geometry/curves/FreemanChain.h"
 
 
//STL
#include <vector>
#include <string>
 
 
#include "CLI11.hpp"
 
 
using namespace DGtal;
 
int main( int argc, char** argv )
{
 
 
// parse command line using CLI ----------------------------------------------
  CLI::App app;
  std::string inputFileName;
  bool oneline {false};
  bool info {false};
  
  app.description("Transform freeman chain into a Sequence of Discrete Points. Result is given to std output.\n Example:\n freeman2sdp ${DGtal}/tests/samples/contourS.fc > contourS.sdp \n");
  app.add_option("-i,--input,1", inputFileName, "Input freeman chain file name." )
    ->required()
    ->check(CLI::ExistingFile);
  app.add_flag("-o,--oneLine", oneline, "output the digital contour in one line like: X0 Y0 X1 Y1 ... XN YN");
  app.add_flag("--info", info, "adds some info as comments at the beginning of the file.");
 
 
  app.get_formatter()->column_width(40);
  CLI11_PARSE(app, argc, argv);
  // END parse command line using CLI ----------------------------------------------
 
 
  typedef FreemanChain<Z2i::Integer> FreemanChain;
 
  std::vector< FreemanChain > vectFcs =  PointListReader< Z2i::Point >:: getFreemanChainsFromFile<Z2i::Integer> (inputFileName); 
    
  for(unsigned int i=0; i< vectFcs.size(); i++){
    bool isClosed = vectFcs.at(i).isClosed(); 
    std::cout << "# grid curve " << i+1 << "/" << vectFcs.size()
              << ( (isClosed)?" closed":" open" ) << std::endl;
    if (info)
      std::cout << "# SDP contour" << i+1<< "/" << vectFcs.size() << " "
                << "# size=" << vectFcs.at(i).size() << std::endl;
    std::vector<Z2i::Point> vectPts; 
    FreemanChain::getContourPoints( vectFcs.at(i), vectPts ); 
    for(unsigned int k=0; k < vectPts.size(); k++){
      std::cout << vectPts.at(k)[0] << " "<< vectPts.at(k)[1] ;
      if(!oneline)
      {
        std::cout << std::endl;
      }
      else
      {
        std::cout << " "; 
      }
          
    }
    std::cout << std::endl;
  }
  
}