PointListReader.ih

// IMPLEMENTATION of inline methods.

#include <cstdlib>
#include <sstream>
#include <fstream>



// Implementation of inline methods                                          //




template<typename TPoint>
inline
std::vector<TPoint> 
DGtal::PointListReader<TPoint>::getPointsFromFile (const std::string &filename,  std::vector<unsigned int> aVectPosition)
{
  ifstream infile;
  infile.open (filename.c_str(), ifstream::in);
  return PointListReader<TPoint>::getPointsFromInputStream(infile, aVectPosition);
}



template<typename TPoint>
inline
std::vector<TPoint> 
DGtal::PointListReader<TPoint>::getPointsFromInputStream (std::istream &in,  std::vector<unsigned int>  aVectPosition)
{
  if(aVectPosition.size()==0){
    for(unsigned int i=0; i<TPoint::dimension; i++){
      aVectPosition.push_back(i);
    }
  }
  vector<TPoint> vectResult;
  string str;
  getline(in, str );
  while ( in.good() ){
    if ( ( str != "" ) && ( str[ 0 ] != '#' ) ){
      istringstream in_str( str );
      unsigned int idx = 0;
      typename TPoint::Component val;
      unsigned int nbFound=0;
      TPoint p;
      while ( in_str.good()&& (nbFound<TPoint::dimension)){
  bool isOK = (in_str >> val);
  for(unsigned int j=0; j< TPoint::dimension; j++){
    if (isOK && (idx == aVectPosition.at(j)) ){
      nbFound++;
      p[j]=val;
    }
  }
  ++idx;
      }      
      if(nbFound==TPoint::dimension){
  vectResult.push_back(p);
      }
    }
    getline(in, str );
  }
  return vectResult;

}




template<typename TPoint>
inline
std::vector< std::vector<TPoint> > 
DGtal::PointListReader<TPoint>::getPolygonsFromFile(const std::string &filename){
  ifstream infile;
  infile.open (filename.c_str(), ifstream::in);
  return PointListReader<TPoint>::getPolygonsFromInputStream(infile);
}


template<typename TPoint>
inline
std::vector< std::vector<TPoint> > 
DGtal::PointListReader<TPoint>::getPolygonsFromInputStream(std::istream & in){
  vector< vector< TPoint > > vectResult;
  string str;
  getline(in, str );
  while ( in.good() ){
    if ( ( str != "" ) && ( str[ 0 ] != '#' ) ){
      vector <TPoint> aContour;
      istringstream in_str( str );
      double x=0.0, y=0.0;
      while ( in_str.good()){
  bool isOK = (in_str >> x);
  isOK = isOK &&  (in_str >>y);
  if(isOK){
    TPoint p;
    p[0]=x;
    p[1]=y;
    aContour.push_back(p);
  }
      }
      vectResult.push_back(aContour);
    }
    getline(in, str );
  }
  return vectResult;

}


  
  



template<typename TPoint>
template<typename TInteger>
inline
std::vector<  FreemanChain< TInteger>  >
DGtal::PointListReader<TPoint>::getFreemanChainsFromFile (const std::string &filename){
  std::vector< FreemanChain<  TInteger> >   vectResult;
  ifstream infile;
  infile.open (filename.c_str(), ifstream::in);
  string str;
  getline(infile, str );
  while ( infile.good() ){
    if ( ( str != "" ) && ( str[ 0 ] != '#' ) ){
      istringstream in_str( str );
      int x0, y0;
      string fcChain;
      bool isOK = (in_str >> x0) && (in_str >> y0) && (in_str >> fcChain);
      FreemanChain< TInteger> fc(fcChain, x0, y0);
      if(isOK){
  vectResult.push_back(fc);
      }else{
  cerr << "Ignoring entry invalid FreemanChain" << endl;
      }
    }
    getline(infile, str );
  }  
  return vectResult;
}
  
  


//                                                                           //