1 /*=========================================================================
2 Program: vv http://www.creatis.insa-lyon.fr/rio/vv
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
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.
13 It is distributed under dual licence
15 - BSD See included LICENSE.txt file
16 - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
17 ===========================================================================**/
19 #ifndef CLITKRESAMPLEIMAGEGENERICFILTER_TXX
20 #define CLITKRESAMPLEIMAGEGENERICFILTER_TXX
23 #include "clitkResampleImageWithOptionsFilter.h"
25 //--------------------------------------------------------------------
26 // Update with the number of dimensions and the pixeltype
27 //--------------------------------------------------------------------
28 template<class InputImageType>
30 clitk::ResampleImageGenericFilter::UpdateWithInputImageType()
34 typename InputImageType::Pointer input = this->template GetInput<InputImageType>(0);
37 typedef typename InputImageType::PixelType PixelType;
38 typedef InputImageType OutputImageType; // to change to float is user ask it (?)
41 typedef clitk::ResampleImageWithOptionsFilter<InputImageType, OutputImageType> ResampleImageFilterType;
42 typename ResampleImageFilterType::Pointer filter = ResampleImageFilterType::New();
43 filter->SetInput(input);
46 filter->SetVerboseOptions(mArgsInfo.verbose_flag);
49 static const unsigned int dim = OutputImageType::ImageDimension;
50 typename OutputImageType::SpacingType spacing;
51 typename OutputImageType::SizeType size;
52 typename OutputImageType::PointType origin;
53 typename OutputImageType::DirectionType direction;
55 if (mArgsInfo.like_given) {
56 itk::ImageIOBase::Pointer header = clitk::readImageHeader(mArgsInfo.like_arg);
58 for(unsigned int i=0; i<dim; i++){
59 spacing[i] = header->GetSpacing(i);
60 size[i] = header->GetDimensions(i);
61 origin[i] = header->GetOrigin(i);
63 for(unsigned int i=0; i<dim; i++) {
64 for(unsigned int j=0;j<dim;j++) {
65 direction(i,j) = header->GetDirection(i)[j];
68 filter->SetOutputSpacing(spacing);
69 filter->SetOutputSize(size);
70 filter->SetOutputOrigin(origin);
71 filter->SetOutputDirection(direction);
74 std::cerr << "*** Warning : I could not read '" << mArgsInfo.like_arg << "' ***" << std::endl;
79 if (mArgsInfo.spacing_given == 1) {
80 filter->SetOutputIsoSpacing(mArgsInfo.spacing_arg[0]);
82 else if ((mArgsInfo.spacing_given != 0) && (mArgsInfo.size_given != 0)) {
83 std::cerr << "Error: use spacing or size, not both." << std::endl;
86 else if (mArgsInfo.spacing_given) {
87 if ((mArgsInfo.spacing_given != 0) && (mArgsInfo.spacing_given != dim)) {
88 std::cerr << "Error: spacing should have one or " << dim << " values." << std::endl;
91 for(unsigned int i=0; i<dim; i++)
92 spacing[i] = mArgsInfo.spacing_arg[i];
93 filter->SetOutputSpacing(spacing);
95 else if (mArgsInfo.size_given) {
96 if ((mArgsInfo.size_given != 0) && (mArgsInfo.size_given != dim)) {
97 std::cerr << "Error: size should have " << dim << " values." << std::endl;
100 for(unsigned int i=0; i<dim; i++)
101 size[i] = mArgsInfo.size_arg[i];
102 filter->SetOutputSize(size);
104 itk::ImageIOBase::Pointer header = clitk::readImageHeader(mArgsInfo.input_arg);
105 for(unsigned int i=0; i<dim; i++){
106 origin[i] = header->GetOrigin(i);
108 for(unsigned int i=0; i<dim; i++) {
109 for(unsigned int j=0;j<dim;j++) {
110 direction(i,j) = header->GetDirection(i)[j];
113 filter->SetOutputOrigin(origin);
114 filter->SetOutputDirection(direction);
117 // Set temporal dimension
118 filter->SetLastDimensionIsTime(mArgsInfo.time_flag);
121 filter->SetGaussianFilteringEnabled(mArgsInfo.autogauss_flag);
122 if (mArgsInfo.gauss_given != 0) {
123 typename ResampleImageFilterType::GaussianSigmaType g;
124 for(unsigned int i=0; i<dim; i++) {
125 g[i] = mArgsInfo.gauss_arg[i];
127 filter->SetGaussianSigma(g);
131 std::string interp = std::string(mArgsInfo.interp_arg);
132 if (interp == "nn") {
133 filter->SetInterpolationType(ResampleImageFilterType::NearestNeighbor);
135 if (interp == "linear") {
136 filter->SetInterpolationType(ResampleImageFilterType::Linear);
138 if (interp == "bspline") {
139 filter->SetInterpolationType(ResampleImageFilterType::BSpline);
141 if (interp == "blut") {
142 filter->SetInterpolationType(ResampleImageFilterType::B_LUT);
144 if (interp == "windowed sinc") {
145 filter->SetInterpolationType(ResampleImageFilterType::WSINC);
147 std::cerr << "Error. I do not know interpolation '" << mArgsInfo.interp_arg
148 << "'. Choose among: nn, linear, bspline, blut, windowed sinc" << std::endl;
156 // Set default pixel value
157 filter->SetDefaultPixelValue(mArgsInfo.default_arg);
160 if (mArgsInfo.thread_given) {
161 filter->SetNumberOfThreads(mArgsInfo.thread_arg);
166 typename OutputImageType::Pointer outputImage = filter->GetOutput();
167 this->template SetNextOutput<OutputImageType>(outputImage);
169 //--------------------------------------------------------------------
171 #endif /* end #define CLITKRESAMPLEIMAGEGENERICFILTER_TXX */