X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=itk%2FclitkResampleImageWithOptionsFilter.txx;h=4c949cec33eff1ade7b4fe43cfac3f4bf313a1fe;hb=c24ade7228f4e5e3d79f74bb4ecbd82963cdbebe;hp=97dbb64faaa539ba3e0bd6d81af2d27816e023bb;hpb=cfb0abe4e63621400a2df8d1fe9617d52c48ffa3;p=clitk.git

diff --git a/itk/clitkResampleImageWithOptionsFilter.txx b/itk/clitkResampleImageWithOptionsFilter.txx
index 97dbb64..4c949ce 100644
--- a/itk/clitkResampleImageWithOptionsFilter.txx
+++ b/itk/clitkResampleImageWithOptionsFilter.txx
@@ -37,7 +37,7 @@
 //--------------------------------------------------------------------
 template <class InputImageType, class OutputImageType>
 clitk::ResampleImageWithOptionsFilter<InputImageType, OutputImageType>::
-ResampleImageWithOptionsFilter():itk::ImageToImageFilter<InputImageType, OutputImageType>() 
+ResampleImageWithOptionsFilter():itk::ImageToImageFilter<InputImageType, OutputImageType>()
 {
   static const unsigned int dim = InputImageType::ImageDimension;
   this->SetNumberOfRequiredInputs(1);
@@ -66,7 +66,7 @@ ResampleImageWithOptionsFilter():itk::ImageToImageFilter<InputImageType, OutputI
 template <class InputImageType, class OutputImageType>
 void
 clitk::ResampleImageWithOptionsFilter<InputImageType, OutputImageType>::
-SetInput(const InputImageType * image) 
+SetInput(const InputImageType * image)
 {
   // Process object is not const-correct so the const casting is required.
   this->SetNthInput(0, const_cast<InputImageType *>(image));
@@ -78,7 +78,7 @@ SetInput(const InputImageType * image)
 template <class InputImageType, class OutputImageType>
 void
 clitk::ResampleImageWithOptionsFilter<InputImageType, OutputImageType>::
-GenerateInputRequestedRegion() 
+GenerateInputRequestedRegion()
 {
   // call the superclass's implementation of this method
   Superclass::GenerateInputRequestedRegion();
@@ -99,7 +99,7 @@ GenerateInputRequestedRegion()
 template <class InputImageType, class OutputImageType>
 void
 clitk::ResampleImageWithOptionsFilter<InputImageType, OutputImageType>::
-GenerateOutputInformation() 
+GenerateOutputInformation()
 {
   static const unsigned int dim = InputImageType::ImageDimension;
 
@@ -124,8 +124,8 @@ GenerateOutputInformation()
   if (m_OutputIsoSpacing != -1) { // apply isoSpacing
     for(unsigned int i=0; i<dim; i++) {
       m_OutputSpacing[i] = m_OutputIsoSpacing;
-      // floor() is used to intentionally reduce the number of slices 
-      // because, from a clinical point of view, it's better to 
+      // floor() is used to intentionally reduce the number of slices
+      // because, from a clinical point of view, it's better to
       // remove data than to add data that privously didn't exist.
       if(inputSpacing[i]*m_OutputSpacing[i]<0)
         itkExceptionMacro( << "Input and output spacings don't have the same signs, can't cope with that" );
@@ -194,11 +194,11 @@ GenerateOutputInformation()
 
 //--------------------------------------------------------------------
 template <class InputImageType, class OutputImageType>
-void 
+void
 clitk::ResampleImageWithOptionsFilter<InputImageType, OutputImageType>::
-GenerateData() 
+GenerateData()
 {
-   
+
   // Get input pointer
   InputImagePointer input = dynamic_cast<InputImageType*>(itk::ProcessObject::GetInput(0));
   static const unsigned int dim = InputImageType::ImageDimension;
@@ -229,10 +229,9 @@ GenerateData()
   }
 
   // Compute origin based on image corner
-  typename FilterType::OriginPointType origin = input->GetOrigin();
   for(unsigned int i=0; i<OutputImageType::ImageDimension; i++) {
-    origin[i] -= 0.5 * input->GetSpacing()[i];
-    origin[i] += 0.5 * m_OutputSpacing[i];
+    m_OutputOrigin[i] -= 0.5 * input->GetSpacing()[i];
+    m_OutputOrigin[i] += 0.5 * m_OutputSpacing[i];
   }
 
   // Instance of the transform object to be passed to the resample
@@ -242,7 +241,7 @@ GenerateData()
   filter->SetOutputSpacing(m_OutputSpacing);
   filter->SetOutputOrigin(m_OutputOrigin);
   filter->SetDefaultPixelValue(m_DefaultPixelValue);
-  filter->SetNumberOfThreads(this->GetNumberOfThreads()); 
+  filter->SetNumberOfThreads(this->GetNumberOfThreads());
   filter->SetOutputDirection(m_OutputDirection); // <-- NEEDED if we want to keep orientation (in case of PermutAxes for example)
 
   // Select interpolator
@@ -317,9 +316,9 @@ GenerateData()
 
 //--------------------------------------------------------------------
 template<class InputImageType>
-typename InputImageType::Pointer 
-clitk::ResampleImageSpacing(typename InputImageType::Pointer input, 
-                            typename InputImageType::SpacingType spacing, 
+typename InputImageType::Pointer
+clitk::ResampleImageSpacing(typename InputImageType::Pointer input,
+                            typename InputImageType::SpacingType spacing,
                             int interpolationType)
 {
   typedef clitk::ResampleImageWithOptionsFilter<InputImageType> ResampleFilterType;