DM_r_convert.cpp

//  =================================================================
//   converting points into an R data frame objects
//  =================================================================
 
#include <iostream>
#include <vector>
 
#include "DM/Datamanager.hpp"
#include "DM/RConverter.hpp"
 
using namespace std;
 
 
void DM_r_convert(const char *filename)
{
#ifdef DM_WITH_R_SUPPORT
  DM::RConverter::initEmbeddedR();
 
  //open the odm
  DM::DatamanagerHandle dm = DM::IDatamanager::load(filename,true/*readOnly*/,false/*threadSafety*/);
 
  //if the dm wasn't opened sucessful exit function
  if (!dm)
  {
    cout << "Unable to open ODM '" << filename << "'" << endl;
    return;
  }
 
  //create a layout for conversion
  DM::AddInfoLayoutFactoryHandle f = DM::IAddInfoLayoutFactory::New();
  f->addColumn(DM::ColumnSemantic::Amplitude);        //pre-defined attribute "Amplitude"
  f->addColumn(DM::ColumnSemantic::EchoWidth);        //pre-defined attribute "EchoWidth"
  DM::AddInfoLayoutHandle layout = f->getLayout();    
 
 
  DM::PointIndexHandle pi = dm->getPointIndex();
 
  //convert each leaf into a data frame
  std::vector<void *> dfs;
  for(auto it = pi->beginLeaf(); it != pi->endLeaf(); ++it)
  {
    dfs.push_back( DM::RConverter::asDataFrame( it->beginPoint(), it->endPoint(), layout, true) );
  }
 
  //create a single data frame 
  __int64 idx = 0;
  void *df = DM::RConverter::createDataFrame( pi->sizePoint(), layout, true );
  for(auto it = pi->beginLeaf(); it != pi->endLeaf(); ++it)
  {
    DM::RConverter::fillDataFrame(df, idx, it->beginPoint(), it->endPoint());
    idx += it->sizePoint();
  }
#endif //#ifdef DM_WITH_R_SUPPORT
}