/*=========================================================================
Program: vv http://www.creatis.insa-lyon.fr/rio/vv
/*=========================================================================
Program: vv http://www.creatis.insa-lyon.fr/rio/vv
- University of LYON http://www.universite-lyon.fr/
- Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
- CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
- University of LYON http://www.universite-lyon.fr/
- Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
- CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
//--------------------------------------------------------------------
clitk::VFResampleGenericFilter::VFResampleGenericFilter():
//--------------------------------------------------------------------
clitk::VFResampleGenericFilter::VFResampleGenericFilter():
InitializeImageType<2>();
InitializeImageType<3>();
InitializeImageType<4>();
InitializeImageType<2>();
InitializeImageType<3>();
InitializeImageType<4>();
ADD_IMAGE_TYPE(Dim, float);
}
//--------------------------------------------------------------------
ADD_IMAGE_TYPE(Dim, float);
}
//--------------------------------------------------------------------
std::cerr << "Error, only one components ? Use clitkImageResample instead." << std::endl;
exit(0);
}
typedef typename ImageType::PixelType PixelType;
if (mNbOfComponents == 2) Update_WithDimAndPixelTypeAndComponent<ImageType::ImageDimension,PixelType,2>();
std::cerr << "Error, only one components ? Use clitkImageResample instead." << std::endl;
exit(0);
}
typedef typename ImageType::PixelType PixelType;
if (mNbOfComponents == 2) Update_WithDimAndPixelTypeAndComponent<ImageType::ImageDimension,PixelType,2>();
- if (mNbOfComponents == 3) Update_WithDimAndPixelTypeAndComponent<ImageType::ImageDimension,PixelType,3>();
- if (mNbOfComponents == 4) Update_WithDimAndPixelTypeAndComponent<ImageType::ImageDimension,PixelType,4>();
+ if (mNbOfComponents == 3) Update_WithDimAndPixelTypeAndComponent<ImageType::ImageDimension,PixelType,3>();
+ if (mNbOfComponents == 4) Update_WithDimAndPixelTypeAndComponent<ImageType::ImageDimension,PixelType,4>();
}
//--------------------------------------------------------------------
//--------------------------------------------------------------------
template<unsigned int Dim, class PixelType, unsigned int DimCompo>
}
//--------------------------------------------------------------------
//--------------------------------------------------------------------
template<unsigned int Dim, class PixelType, unsigned int DimCompo>
// Reading input
typedef itk::Vector<PixelType, DimCompo> DisplacementType;
typedef itk::Image< DisplacementType, Dim > ImageType;
// Reading input
typedef itk::Vector<PixelType, DimCompo> DisplacementType;
typedef itk::Image< DisplacementType, Dim > ImageType;
-typename ImageType::Pointer
-clitk::VFResampleGenericFilter::ComputeImage(typename ImageType::Pointer inputImage) {
+typename ImageType::Pointer
+clitk::VFResampleGenericFilter::ComputeImage(typename ImageType::Pointer inputImage)
+{
// Create Image Filter
typedef itk::VectorResampleImageFilter<ImageType,ImageType> FilterType;
typename FilterType::Pointer filter = FilterType::New();
// Create Image Filter
typedef itk::VectorResampleImageFilter<ImageType,ImageType> FilterType;
typename FilterType::Pointer filter = FilterType::New();
// Instance of the transform object to be passed to the resample
// filter. By default, identity transform is applied
typedef itk::AffineTransform<double, ImageType::ImageDimension> TransformType;
// Instance of the transform object to be passed to the resample
// filter. By default, identity transform is applied
typedef itk::AffineTransform<double, ImageType::ImageDimension> TransformType;
typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
filter->SetInterpolator(interpolator);
typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
filter->SetInterpolator(interpolator);
typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
filter->SetInterpolator(interpolator);
typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
filter->SetInterpolator(interpolator);
else gaussianFilters[i]->SetInput(gaussianFilters[i-1]->GetOutput());
}
filter->SetInput(gaussianFilters[ImageType::ImageDimension-1]->GetOutput());
else gaussianFilters[i]->SetInput(gaussianFilters[i-1]->GetOutput());
}
filter->SetInput(gaussianFilters[ImageType::ImageDimension-1]->GetOutput());
- }
- catch( itk::ExceptionObject & err ) {
- std::cerr << "Error while filtering " << mInputFilenames[0].c_str()
- << " " << err << std::endl;
+ } catch( itk::ExceptionObject & err ) {
+ std::cerr << "Error while filtering " << mInputFilenames[0].c_str()
+ << " " << err << std::endl;
}
//--------------------------------------------------------------------
//--------------------------------------------------------------------
}
//--------------------------------------------------------------------
//--------------------------------------------------------------------
mOutputSize.resize(size.size());
std::copy(size.begin(), size.end(), mOutputSize.begin());
}
//--------------------------------------------------------------------
//--------------------------------------------------------------------
mOutputSize.resize(size.size());
std::copy(size.begin(), size.end(), mOutputSize.begin());
}
//--------------------------------------------------------------------
//--------------------------------------------------------------------
mOutputSpacing.resize(spacing.size());
std::copy(spacing.begin(), spacing.end(), mOutputSpacing.begin());
}
//--------------------------------------------------------------------
//--------------------------------------------------------------------
mOutputSpacing.resize(spacing.size());
std::copy(spacing.begin(), spacing.end(), mOutputSpacing.begin());
}
//--------------------------------------------------------------------
//--------------------------------------------------------------------
mInterpolatorName = inter;
}
//--------------------------------------------------------------------
//--------------------------------------------------------------------
mInterpolatorName = inter;
}
//--------------------------------------------------------------------
//--------------------------------------------------------------------
mApplyGaussianFilterBefore = true;
mSigma.resize(sigma.size());
std::copy(sigma.begin(), sigma.end(), mSigma.begin());
mApplyGaussianFilterBefore = true;
mSigma.resize(sigma.size());
std::copy(sigma.begin(), sigma.end(), mSigma.begin());