Preset transformation matrices appear to be overwritten - assume inst…

…ead that one exists
jonclayden · Sep 24, 2024 · 8b61df2 · 8b61df2
1 parent b7e7476
commit 8b61df2
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/reg-lib/_reg_aladin.cpp b/src/reg-lib/_reg_aladin.cpp
@@ -322,10 +322,10 @@ void reg_aladin<T>::InitialiseRegistration()
     reg_thresholdImage<T>(this->FloatingPyramid[l],this->FloatingLowerThreshold, this->FloatingUpperThreshold);
   }
 
-  // Initialise the transformation
+  // Initialise the transformation - for R this is guaranteed to have been done
+#ifndef HAVE_R
   if (this->InputTransformName != NULL)
   {
-#ifndef HAVE_R
     if (FILE *aff = fopen(this->InputTransformName, "r")) {
       fclose(aff);
     }
@@ -338,7 +338,6 @@ void reg_aladin<T>::InitialiseRegistration()
       reg_exit();
     }
     reg_tool_ReadAffineFile(this->TransformationMatrix, this->InputTransformName);
-#endif
   }
   else  // No input affine transformation
   {
@@ -372,6 +371,7 @@ void reg_aladin<T>::InitialiseRegistration()
       this->TransformationMatrix->m[2][3] = floatingRealPosition[2] - referenceRealPosition[2];
     }
   }
+#endif
 }
 /* *************************************************************** */
 template<class T>