initial image centralization (BLUT) is optional

- with --centre flag

diff --git a/registration/clitkBLUTDIR.ggo b/registration/clitkBLUTDIR.ggo
index ed3efce63f5e4c8d4ae1e80b76eaeca1f519fee9..9db5646e25c8ea98b08dde443362c839eb7257d2 100644 (file)
--- a/registration/clitkBLUTDIR.ggo
+++ b/registration/clitkBLUTDIR.ggo
@@ -34,6 +34,7 @@ option "after"                        -       "Difference image after "               string          no
 section "Transform (Note that only one of --control, --spacing is required. The other will be adjusted to fit the region and allow exact representation. SamplingFactor will be set accordingly"
 
 option "initMatrix"            -       "Prior rigid/affine transform matrix from reference to target space"    string          no      
+option "centre"   - "Centre images before starting registration (ignored if initMatrix was given)"  flag    off
 option "initCoeff"             -       "Initial coefficient image"                                             string          no
 option "order"                 -       "Spline Order FFD"                                              int       no    multiple        default="3"
 option "control"               -       "Internal control points for each dimension"                    int       no    multiple        

diff --git a/registration/clitkBLUTDIRGenericFilter.cxx b/registration/clitkBLUTDIRGenericFilter.cxx
index 9f808e7e1ffa24b8126457acf3ea72894e44704e..06045c52defc3a3774844c2bd30334cb4add04a1 100644 (file)
--- a/registration/clitkBLUTDIRGenericFilter.cxx
+++ b/registration/clitkBLUTDIRGenericFilter.cxx
@@ -496,9 +496,9 @@ namespace clitk
         itk::Vector<double,3> finalTranslation = clitk::GetTranslationPartMatrix3D(rigidTransformMatrix);
         rigidTransform->SetTranslation(finalTranslation);
       }
-      else
+      else if (m_ArgsInfo.centre_flag)
       {
-        if(m_Verbose) std::cout<<"No itinial matrix given. Centering all images..."<<std::endl;
+        if(m_Verbose) std::cout<<"No itinial matrix given and \"centre\" flag switched on. Centering all images..."<<std::endl;
         
         rigidTransform=RigidTransformType::New();