]> Creatis software - clitk.git/blob - tools/clitkVFInterpolate.cxx
Update Landmark coordinates with transformation matrix
[clitk.git] / tools / clitkVFInterpolate.cxx
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 CLITKVFRESAMPLE_CXX
19 #define CLITKVFRESAMPLE_CXX
20
21 // clitk
22 #include "clitkVFInterpolate_ggo.h"
23 #include "clitkIO.h"
24 #include "clitkVFInterpolateGenericFilter.h"
25
26 //--------------------------------------------------------------------
27 int main(int argc, char * argv[])
28 {
29
30     // Init command line
31     GGO(clitkVFInterpolate, args_info);
32     CLITK_INIT;
33
34     // Read input image header to check image dimension
35     itk::ImageIOBase::Pointer header = clitk::readImageHeader(args_info.input1_arg);
36     unsigned int dim = header->GetNumberOfDimensions();
37     std::string pixelTypeName = header->GetComponentTypeAsString(header->GetComponentType());
38
39     // Print image info if verbose
40     if (args_info.verbose_flag) {
41         std::cout << "Input image <" << args_info.input1_arg << "> is ";
42         clitk::printImageHeader(header, std::cout);
43         std::cout << std::endl;
44     }
45
46     // Get input size/spacing
47     std::vector<int>    inputSize;
48     std::vector<double> inputSpacing;
49     inputSize.resize(dim);
50     inputSpacing.resize(dim);
51     for (unsigned int i=0; i<dim; i++) {
52         inputSpacing[i] = header->GetSpacing(i);
53         inputSize[i] = header->GetDimensions(i);
54     }
55
56     if (args_info.verbose_flag) {
57         std::cout << "Output image will be : " << std::endl;
58     }
59
60     // Create a filter
61     clitk::VFInterpolateGenericFilter::Pointer filter = clitk::VFInterpolateGenericFilter::New();
62     filter->SetInputFilename(args_info.input1_arg);
63     filter->SetInputFilename2(args_info.input2_arg);
64     filter->SetInterpolationName(args_info.interp_arg);
65     filter->SetDistance(args_info.distance_arg);
66 //     filter->SetBSplineOrder(args_info.order_arg);
67 //     filter->SetBLUTSampling(args_info.sampling_arg);
68     filter->SetOutputFilename(args_info.output_arg);
69
70     // Go !
71     filter->Update();
72
73     // this is the end my friend
74     return 0;
75 }// end main
76 //--------------------------------------------------------------------
77
78 #endif  /* end #define CLITKVFRESAMPLE_CXX */