BUG: Boundary condition was not propagated to FFTPad

Author: phcerdan

include/itkFFTPadPositiveIndexImageFilter.h

@@ -91,8 +91,17 @@ class FFTPadPositiveIndexImageFilter:
   itkGetConstMacro(SizeGreatestPrimeFactor, SizeValueType);
   itkSetMacro(SizeGreatestPrimeFactor, SizeValueType);
   /** Set/get the boundary condition. */
-  itkSetMacro(BoundaryCondition, BoundaryConditionPointerType);
   itkGetConstMacro(BoundaryCondition, BoundaryConditionPointerType);
+  virtual void SetBoundaryCondition(const BoundaryConditionPointerType boundaryCondition)
+    {
+    if ( this->m_BoundaryCondition != boundaryCondition )
+      {
+      this->m_BoundaryCondition = boundaryCondition;
+      this->m_FFTPadFilter->SetBoundaryCondition(this->m_BoundaryCondition);
+      this->m_FFTPadFilter->Modified();
+      this->Modified();
+      }
+    }
 
 protected:
   FFTPadPositiveIndexImageFilter();

include/itkFFTPadPositiveIndexImageFilter.hxx

@@ -53,14 +53,20 @@ FFTPadPositiveIndexImageFilter< TInputImage, TOutputImage >
   const typename OutputImageType::RegionType & outputRequestedRegion =
     outputPtr->GetRequestedRegion();
 
+  typename OutputImageType::RegionType shiftedOutputRequestedRegion;
+  typename OutputImageType::RegionType::IndexType shiftedRequestedIndex =
+    outputRequestedRegion.GetIndex() - this->m_ChangeInfoFilter->GetOutputOffset();
+  shiftedOutputRequestedRegion.SetIndex(shiftedRequestedIndex);
+  shiftedOutputRequestedRegion.SetSize(outputRequestedRegion.GetSize());
+
   // Ask the boundary condition for the input requested region.
   if ( !m_BoundaryCondition )
     {
     itkExceptionMacro( << "Boundary condition is ITK_NULLPTR so no request region can be generated.");
     }
   typename InputImageType::RegionType inputRequestedRegion =
     m_BoundaryCondition->GetInputRequestedRegion( inputLargestPossibleRegion,
-      outputRequestedRegion );
+      shiftedOutputRequestedRegion );
 
   inputPtr->SetRequestedRegion( inputRequestedRegion );
 }