registration/clitkGenericVectorInterpolator.txx

   1 /*=========================================================================
   2   Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
   3
   4   Authors belong to:
   5   - University of LYON              http://www.universite-lyon.fr/
   6   - Léon Bérard cancer center       http://www.centreleonberard.fr
   7   - CREATIS CNRS laboratory         http://www.creatis.insa-lyon.fr
   8
   9   This software is distributed WITHOUT ANY WARRANTY; without even
  10   the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  11   PURPOSE.  See the copyright notices for more information.
  12
  13   It is distributed under dual licence
  14
  15   - BSD        See included LICENSE.txt file
  16   - CeCILL-B   http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
  17 ===========================================================================**/
  18 #ifndef __clitkGenericVectorInterpolator_txx
  19 #define __clitkGenericVectorInterpolator_txx
  20 #include "clitkGenericVectorInterpolator.h"
  21
  22
  23 namespace clitk
  24 {
  25
  26   //=========================================================================================================================
  27   //constructor
  28   template <class args_info_type, class ImageType, class TCoordRep>
  29   GenericVectorInterpolator<args_info_type, ImageType, TCoordRep>::GenericVectorInterpolator()
  30   {
  31     m_Interpolator=NULL;
  32     m_Verbose=false;
  33   }
  34
  35
  36   //=========================================================================================================================
  37   //Get the pointer
  38   template <class args_info_type, class ImageType, class TCoordRep>
  39   typename GenericVectorInterpolator<args_info_type, ImageType, TCoordRep>::InterpolatorPointer
  40   GenericVectorInterpolator<args_info_type, ImageType, TCoordRep>::GetInterpolatorPointer()
  41   {
  42     //============================================================================
  43     // We retrieve the type of interpolation from the command line
  44     //============================================================================
  45     typename InterpolatorType::Pointer interpolator;
  46
  47     switch ( m_ArgsInfo.interpVF_arg )
  48       {
  49       case 0:
  50
  51         interpolator= itk::VectorNearestNeighborInterpolateImageFunction< ImageType,TCoordRep >::New();
  52         if (m_Verbose) std::cout<<"Using nearestneighbor interpolation..."<<std::endl;
  53         break;
  54
  55       case 1:
  56
  57         interpolator = itk::VectorLinearInterpolateImageFunction< ImageType,TCoordRep >::New();
  58         if (m_Verbose) std::cout<<"Using linear interpolation..."<<std::endl;
  59         break;
  60
  61       case 2:
  62         {
  63           typename clitk::VectorBSplineInterpolateImageFunction< ImageType,TCoordRep >::Pointer m =clitk::VectorBSplineInterpolateImageFunction< ImageType,TCoordRep >::New();
  64           m->SetSplineOrder(m_ArgsInfo.interpVFOrder_arg);
  65           interpolator=m;
  66           if (m_Verbose) std::cout<<"Using Bspline interpolation..."<<std::endl;
  67           break;
  68         }
  69
  70       case 3:
  71         {
  72           typename clitk::VectorBSplineInterpolateImageFunctionWithLUT< ImageType,TCoordRep >::Pointer m =clitk::VectorBSplineInterpolateImageFunctionWithLUT< ImageType,TCoordRep >::New();
  73           m->SetSplineOrder(m_ArgsInfo.interpVFOrder_arg);
  74           m->SetLUTSamplingFactor(m_ArgsInfo.interpVFSF_arg);
  75           interpolator=m;
  76           if (m_Verbose) std::cout<<"Using BLUT interpolation..."<<std::endl;
  77           break;
  78         }
  79
  80       }//end of switch
  81
  82
  83     //============================================================================
  84     //return the pointer
  85     return interpolator;
  86   }
  87
  88 }
  89
  90 #endif