X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=filters%2FclitkImageResampleGenericFilter.cxx;h=aa84b8e4fba1eafcced7237585e9a7a043954e71;hb=1bbbd127f84b5f5949b2cbb82badb9193b3f59e8;hp=da510c6d1e50b712a3480fba1be7e92b668770c8;hpb=9bcf64b7f7c5a7f0a39843df634c05b5d9e4e237;p=clitk.git

diff --git a/filters/clitkImageResampleGenericFilter.cxx b/filters/clitkImageResampleGenericFilter.cxx
index da510c6..aa84b8e 100644
--- a/filters/clitkImageResampleGenericFilter.cxx
+++ b/filters/clitkImageResampleGenericFilter.cxx
@@ -1,21 +1,53 @@
+/*=========================================================================
+  Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
+
+  Authors belong to:
+  - 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
+
+  This software is distributed WITHOUT ANY WARRANTY; without even
+  the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+  PURPOSE.  See the copyright notices for more information.
+
+  It is distributed under dual licence
+
+  - BSD        See included LICENSE.txt file
+  - CeCILL-B   http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
+======================================================================-====*/
 #ifndef CLITKIMAGERESAMPLEGENERICFILTER2_CXX
 #define CLITKIMAGERESAMPLEGENERICFILTER2_CXX
-
 /**
  -------------------------------------------------------------------
  * @file   clitkImageResampleGenericFilter.cxx
  * @author David Sarrut <David.Sarrut@creatis.insa-lyon.fr>
  * @date   23 Feb 2008 08:37:53
 
- * @brief  
+ * @brief
 
  -------------------------------------------------------------------*/
 
 #include "clitkImageResampleGenericFilter.h"
 
+// itk include
+#include "itkImage.h"
+#include "itkImageFileReader.h"
+#include "itkImageSeriesReader.h"
+#include "itkImageFileWriter.h"
+#include "itkRecursiveGaussianImageFilter.h"
+#include "itkResampleImageFilter.h"
+#include "itkAffineTransform.h"
+#include "itkNearestNeighborInterpolateImageFunction.h"
+#include "itkWindowedSincInterpolateImageFunction.h"
+#include "itkLinearInterpolateImageFunction.h"
+#include "itkBSplineInterpolateImageFunction.h"
+#include "itkBSplineInterpolateImageFunctionWithLUT.h"
+#include "itkCommand.h"
+
 //--------------------------------------------------------------------
 clitk::ImageResampleGenericFilter::ImageResampleGenericFilter():
-  ImageToImageGenericFilter<Self>("ImageResample") {
+  ImageToImageGenericFilter<Self>("ImageResample")
+{
   mApplyGaussianFilterBefore = false;
   mDefaultPixelValue = 0.0;
   mInterpolatorName = "NN";
@@ -29,18 +61,17 @@ clitk::ImageResampleGenericFilter::ImageResampleGenericFilter():
 
 //--------------------------------------------------------------------
 template<unsigned int Dim>
-void clitk::ImageResampleGenericFilter::InitializeImageTypeWithDim() {      
-  ADD_IMAGE_TYPE(Dim, char);
-  ADD_IMAGE_TYPE(Dim, short);
-  ADD_IMAGE_TYPE(Dim, int);
-  ADD_IMAGE_TYPE(Dim, float);
+void clitk::ImageResampleGenericFilter::InitializeImageTypeWithDim()
+{
+  ADD_DEFAULT_IMAGE_TYPES(Dim);
 }
 //--------------------------------------------------------------------
 
 
 //--------------------------------------------------------------------
 template<class InputImageType>
-void clitk::ImageResampleGenericFilter::UpdateWithInputImageType() {
+void clitk::ImageResampleGenericFilter::UpdateWithInputImageType()
+{
 
   // Some typedefs
   typedef typename InputImageType::SizeType    SizeType;
@@ -78,7 +109,7 @@ void clitk::ImageResampleGenericFilter::UpdateWithInputImageType() {
   // Create Image Filter
   typedef itk::ResampleImageFilter<InputImageType,InputImageType> 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, InputImageType::ImageDimension> TransformType;
@@ -102,36 +133,38 @@ void clitk::ImageResampleGenericFilter::UpdateWithInputImageType() {
 
   // Select interpolator
   if (mInterpolatorName == "nn") {
-    typedef itk::NearestNeighborInterpolateImageFunction<InputImageType, double> InterpolatorType;     
+    typedef itk::NearestNeighborInterpolateImageFunction<InputImageType, double> InterpolatorType;
     typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
     filter->SetInterpolator(interpolator);
-  }
-  else { 
+  } else {
     if (mInterpolatorName == "linear") {
-      typedef itk::LinearInterpolateImageFunction<InputImageType, double> InterpolatorType;     
+      typedef itk::LinearInterpolateImageFunction<InputImageType, double> InterpolatorType;
       typename InterpolatorType::Pointer interpolator =  InterpolatorType::New();
       filter->SetInterpolator(interpolator);
-    }
-    else {
-      if (mInterpolatorName == "bspline") {
-	typedef itk::BSplineInterpolateImageFunction<InputImageType, double> InterpolatorType; 
-	typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
-	interpolator->SetSplineOrder(mBSplineOrder);
-	filter->SetInterpolator(interpolator);
-      }
-      else {
-	if (mInterpolatorName == "blut") {
-	  typedef itk::BSplineInterpolateImageFunctionWithLUT<InputImageType, double> InterpolatorType; 
-	  typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
-	  interpolator->SetSplineOrder(mBSplineOrder);
-	  interpolator->SetLUTSamplingFactor(mSamplingFactors[0]);
-	  filter->SetInterpolator(interpolator);
-	}
-	else {
-	  std::cerr << "Sorry, I do not know the interpolator '" << mInterpolatorName 
-		    << "'. Known interpolators are :  nn, linear, bspline, blut" << std::endl;
-	  exit(0);
-	}
+    } else {
+      if (mInterpolatorName == "windowed sinc") {
+        typedef itk::WindowedSincInterpolateImageFunction<InputImageType, 4> InterpolatorType;
+        typename InterpolatorType::Pointer interpolator =  InterpolatorType::New();
+        filter->SetInterpolator(interpolator);
+      } else {
+        if (mInterpolatorName == "bspline") {
+          typedef itk::BSplineInterpolateImageFunction<InputImageType, double> InterpolatorType;
+          typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
+          interpolator->SetSplineOrder(mBSplineOrder);
+          filter->SetInterpolator(interpolator); 
+        } else {
+          if (mInterpolatorName == "blut") {
+            typedef itk::BSplineInterpolateImageFunctionWithLUT<InputImageType, double> InterpolatorType;
+            typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
+            interpolator->SetSplineOrder(mBSplineOrder);
+            interpolator->SetLUTSamplingFactor(mSamplingFactors[0]);
+            filter->SetInterpolator(interpolator);
+          } else {
+            std::cerr << "Sorry, I do not know the interpolator '" << mInterpolatorName
+              << "'. Known interpolators are :  nn, linear, bspline, blut" << std::endl;
+            exit(0);
+          }
+        }
       }
     }
   }
@@ -153,23 +186,21 @@ void clitk::ImageResampleGenericFilter::UpdateWithInputImageType() {
       else gaussianFilters[i]->SetInput(gaussianFilters[i-1]->GetOutput());
     }
     filter->SetInput(gaussianFilters[InputImageType::ImageDimension-1]->GetOutput());
-  }
-  else {
+  } else {
     filter->SetInput(input);
   }
 
   // Go !
-  try { 
+  try {
     filter->Update();
-  }
-  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;
     exit(0);
   }
 
   // Get result
-   typename InputImageType::Pointer outputImage = filter->GetOutput();
+  typename InputImageType::Pointer outputImage = filter->GetOutput();
 
   // Write/save results
   this->SetNextOutput<InputImageType>(outputImage);
@@ -179,27 +210,31 @@ void clitk::ImageResampleGenericFilter::UpdateWithInputImageType() {
 
 
 //--------------------------------------------------------------------
-void clitk::ImageResampleGenericFilter::SetOutputSize(const std::vector<int> & size) {
+void clitk::ImageResampleGenericFilter::SetOutputSize(const std::vector<int> & size)
+{
   mOutputSize.resize(size.size());
   std::copy(size.begin(), size.end(), mOutputSize.begin());
 }
 //--------------------------------------------------------------------
 
 //--------------------------------------------------------------------
-void clitk::ImageResampleGenericFilter::SetOutputSpacing(const std::vector<double> & spacing) {
+void clitk::ImageResampleGenericFilter::SetOutputSpacing(const std::vector<double> & spacing)
+{
   mOutputSpacing.resize(spacing.size());
   std::copy(spacing.begin(), spacing.end(), mOutputSpacing.begin());
 }
 //--------------------------------------------------------------------
 
 //--------------------------------------------------------------------
-void clitk::ImageResampleGenericFilter::SetInterpolationName(const std::string & inter) {
+void clitk::ImageResampleGenericFilter::SetInterpolationName(const std::string & inter)
+{
   mInterpolatorName = inter;
 }
 //--------------------------------------------------------------------
 
 //--------------------------------------------------------------------
-void clitk::ImageResampleGenericFilter::SetGaussianSigma(const std::vector<double> & sigma) {
+void clitk::ImageResampleGenericFilter::SetGaussianSigma(const std::vector<double> & sigma)
+{
   mApplyGaussianFilterBefore = true;
   mSigma.resize(sigma.size());
   std::copy(sigma.begin(), sigma.end(), mSigma.begin());