X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;ds=sidebyside;f=registration%2FclitkBLUTDIRGenericFilter.cxx;fp=registration%2FclitkBLUTDIRGenericFilter.cxx;h=9640282ec800d2cf1d3a4aa622b3f244f6248c8d;hb=2d83e4d660b0fb679977954c0ffcb8a34e8eb455;hp=3e6c6ec18b36d8569b9d37e2e0e10a179e69f2c2;hpb=dded8b0644eb54aecd495c15cef179001aefb234;p=clitk.git

diff --git a/registration/clitkBLUTDIRGenericFilter.cxx b/registration/clitkBLUTDIRGenericFilter.cxx
index 3e6c6ec..9640282 100755
--- a/registration/clitkBLUTDIRGenericFilter.cxx
+++ b/registration/clitkBLUTDIRGenericFilter.cxx
@@ -340,13 +340,23 @@ namespace clitk
         }
         if (m_Verbose)std::cout <<"Reference image mask was read..." <<std::endl;
 
+        // Resample mask
+        typedef itk::ResampleImageFilter<ImageMaskType, ImageMaskType> ResamplerType;
+        typename ResamplerType::Pointer resampler = ResamplerType::New();
+        typedef itk::NearestNeighborInterpolateImageFunction<ImageMaskType, TCoordRep> InterpolatorType;
+        typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
+        resampler->SetOutputParametersFromImage(fixedImage);
+        resampler->SetInterpolator(interpolator);
+        resampler->SetInput(maskReader->GetOutput());
+        resampler->Update();
+
         // Set the image to the spatialObject
-        fixedMask->SetImage( maskReader->GetOutput() );
+        fixedMask->SetImage(resampler->GetOutput());
 
         // Find the bounding box of the "inside" label
         typedef itk::LabelGeometryImageFilter<ImageMaskType> GeometryImageFilterType;
         typename GeometryImageFilterType::Pointer geometryImageFilter=GeometryImageFilterType::New();
-        geometryImageFilter->SetInput(maskReader->GetOutput());
+        geometryImageFilter->SetInput(resampler->GetOutput());
         geometryImageFilter->Update();
         typename GeometryImageFilterType::BoundingBoxType boundingBox = geometryImageFilter->GetBoundingBox(1);