itk/clitkGenericVectorInterpolator.txx

   1 #ifndef __clitkGenericVectorInterpolator_txx
   2 #define __clitkGenericVectorInterpolator_txx
   3 #include "clitkGenericVectorInterpolator.h"
   4
   5
   6 namespace clitk
   7 {
   8
   9   //=========================================================================================================================
  10   //constructor
  11   template <class args_info_type, class ImageType, class TCoordRep>
  12   GenericVectorInterpolator<args_info_type, ImageType, TCoordRep>::GenericVectorInterpolator()
  13   {
  14     m_Interpolator=NULL;
  15     m_Verbose=false;
  16   }
  17
  18
  19   //=========================================================================================================================
  20   //Get the pointer
  21   template <class args_info_type, class ImageType, class TCoordRep>
  22   typename GenericVectorInterpolator<args_info_type, ImageType, TCoordRep>::InterpolatorPointer
  23   GenericVectorInterpolator<args_info_type, ImageType, TCoordRep>::GetInterpolatorPointer()
  24   {
  25     //============================================================================
  26     // We retrieve the type of interpolation from the command line
  27     //============================================================================
  28     typename InterpolatorType::Pointer interpolator;
  29
  30     switch ( m_ArgsInfo.interpVF_arg )
  31       {
  32       case 0:
  33
  34         interpolator= itk::VectorNearestNeighborInterpolateImageFunction< ImageType,TCoordRep >::New();
  35         if (m_Verbose) std::cout<<"Using nearestneighbor interpolation..."<<std::endl;
  36         break;
  37
  38       case 1:
  39
  40         interpolator = itk::VectorLinearInterpolateImageFunction< ImageType,TCoordRep >::New();
  41         if (m_Verbose) std::cout<<"Using linear interpolation..."<<std::endl;
  42         break;
  43
  44       case 2:
  45         {
  46           typename clitk::VectorBSplineInterpolateImageFunction< ImageType,TCoordRep >::Pointer m =clitk::VectorBSplineInterpolateImageFunction< ImageType,TCoordRep >::New();
  47           m->SetSplineOrder(m_ArgsInfo.interpVFOrder_arg);
  48           interpolator=m;
  49           if (m_Verbose) std::cout<<"Using Bspline interpolation..."<<std::endl;
  50           break;
  51         }
  52
  53       case 3:
  54         {
  55           typename clitk::VectorBSplineInterpolateImageFunctionWithLUT< ImageType,TCoordRep >::Pointer m =clitk::VectorBSplineInterpolateImageFunctionWithLUT< ImageType,TCoordRep >::New();
  56           m->SetSplineOrder(m_ArgsInfo.interpVFOrder_arg);
  57           m->SetLUTSamplingFactor(m_ArgsInfo.interpVFSF_arg);
  58           interpolator=m;
  59           if (m_Verbose) std::cout<<"Using BLUT interpolation..."<<std::endl;
  60           break;
  61         }
  62
  63       }//end of switch
  64
  65
  66     //============================================================================
  67     //return the pointer
  68     return interpolator;
  69   }
  70
  71 }
  72
  73 #endif