ODFNormalization.cxx

Go to the documentation of this file.

#include "utlITK.h"

#include "itkMultiplyByConstantVectorImageFilter.h"
#include "ODFNormalizationCLP.h"

#include "itkNormalizeODFImageFilter.h"

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

  typedef double FloatType;

  typedef itk::VectorImage<FloatType,3> InputImageType;
  typedef itk::VectorImage<FloatType,3> OutputImageType;
  
  InputImageType::Pointer input = InputImageType::New();
  itk::ReadImage(_InputFile, input);


  typedef itk::NormalizeODFImageFilter<InputImageType, OutputImageType>  FilterType;
  FilterType::Pointer filter = FilterType::New();
  
  filter->SetInput(input);

  if (_Type=="SH")
    filter->SetODFType(FilterType::FunctorType::SH);
  else if (_Type=="SAMPLE")
    filter->SetODFType(FilterType::FunctorType::SAMPLE);

  filter->Update();
  
  OutputImageType::Pointer output = filter->GetOutput();
  itk::SaveImage(output, _OutputFile);

  return 0;
}