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"


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

//STL
#include <vector>
#include <string>

using namespace DGtal;

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")
    ("input,i", po::value<std::string>(), "Input freeman chain file name.")
    ("info", "adds some info as comments at the beginning of the file.") 
    ("oneLine,o", " output the digital contour in one line like: X0 Y0 X1 Y1 ... XN YN");

  
  
  typedef FreemanChain<Z2i::Integer> FreemanChain;
  
  
  bool parseOK=true;
  po::variables_map vm;
  try{
    po::store(po::parse_command_line(argc, argv, general_opt), vm);  
  }catch(const std::exception& ex){
    parseOK=false;
    trace.info()<< "Error checking program options: "<< ex.what()<< std::endl;
  }

  po::notify(vm);    
  if(!parseOK||vm.count("help")||argc<=1 || (!(vm.count("input")) ) )
    {
      trace.info()<< "Transform freeman chain into a Sequence of Discrete Points. Result is given to std output. " <<std::endl << "Basic usage: "<<std::endl
                  << "\t freeman2sdp [input] > out.sdp  "<<std::endl
                  << general_opt << "\n";
      trace.info() << "Example:\n"
                   << "freeman2sdp -i ${DGtal}/tests/samples/contourS.fc > contourS.sdp \n";

      return 0;
    }
  
  bool oneline = vm.count("oneLine");
  if( vm.count("input") ){
    std::string fileName = vm["input"].as<std::string>();
    std::vector< FreemanChain > vectFcs =  PointListReader< Z2i::Point >:: getFreemanChainsFromFile<Z2i::Integer> (fileName); 
    
    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 ( vm.count("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;
    }
    
  }
  

}