]> Creatis software - clitk.git/blob - registration/clitkGenericVectorInterpolator.txx
Moved from repository clitk to clitk.private/tests_dav
[clitk.git] / 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