Correction of an orientation error in the TODs estimated with animaTO…

…DEstimator
Inria-Empenn · Apr 29, 2024 · 7023720 · 7023720
1 parent 185203f
commit 7023720
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/Anima/diffusion/odf/tod_estimator/animaTODEstimatorImageFilter.hxx b/Anima/diffusion/odf/tod_estimator/animaTODEstimatorImageFilter.hxx
@@ -84,6 +84,14 @@ namespace anima
             for (int i = 0; i < nbPoints; i++)
             {
                 DirType dir = GetFiberDirection(i, fiber);
+                if(output->GetDirection()[0][0] > 0)
+                    dir[0] *= -1;
+                if(output->GetDirection()[1][1] < 0)
+                    dir[1] *= -1;
+                if(output->GetDirection()[2][2] < 0)
+                    dir[2] *= -1;
+
+
                 PointType point = GetCenterVoxel(i, fiber);
                 itk::Index<3> index;