SamplingSchemeQSpaceRandomEstimation.cxx

Go to the documentation of this file.

#include "utl.h"
#include "itkSamplingSchemeQSpaceWriter.h"
#include "itkSamplingSchemeQSpace.h"

#include "SamplingSchemeQSpaceRandomEstimationCLP.h"

int 
main (int argc, char const* argv[])
{
  
  PARSE_ARGS;

  typedef itk::SamplingSchemeQSpace<double>  SamplingType;
  typedef SamplingType::Pointer  SamplingPointer;

  utlGlobalException(_NumberOfSamples.size()==0, "should give the numbers of samples in shells");

  SamplingPointer scheme = SamplingType::New();
  scheme->GenerateFromRandomPoints(_NumberOfSamples);

  typedef itk::SamplingSchemeQSpaceWriter<SamplingType> WriterType;
  WriterType::Pointer writer = WriterType::New();
  writer->SetSampling(scheme);
  writer->SetOrientationFile(_OutputOrientations);
  writer->SaveSingleShellOn();
  writer->SaveAllShellsInOneFileOff();
  writer->Update();
  
  return 0;
}