Use ImageGeom::getCoordsf to get voxel center for ComputeShapes

* Fixes inconsistency between filters with where the voxel center is assumed to be Signed-off-by: Jared Duffey <[email protected]>
BlueQuartzSoftware · Nov 1, 2024 · 306a14d · 306a14d
1 parent 9bee9ce
commit 306a14d
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/src/Plugins/OrientationAnalysis/src/OrientationAnalysis/Filters/Algorithms/ComputeShapes.cpp b/src/Plugins/OrientationAnalysis/src/OrientationAnalysis/Filters/Algorithms/ComputeShapes.cpp
@@ -199,9 +199,10 @@ void ComputeShapes::find_moments()
       for(size_t k = 0; k < xPoints; k++)
       {
         int32_t gnum = featureIds[zStride + yStride + k];
-        x = float(k * modXRes) + (origin[0] * static_cast<float>(m_ScaleFactor));
-        y = float(j * modYRes) + (origin[1] * static_cast<float>(m_ScaleFactor));
-        z = float(i * modZRes) + (origin[2] * static_cast<float>(m_ScaleFactor));
+        FloatVec3 voxelCenter = imageGeom.getCoordsf(k, j, i);
+        x = voxelCenter[0] * static_cast<float>(m_ScaleFactor);
+        y = voxelCenter[1] * static_cast<float>(m_ScaleFactor);
+        z = voxelCenter[2] * static_cast<float>(m_ScaleFactor);
         x1 = x + (modXRes / 4.0f);
         x2 = x - (modXRes / 4.0f);
         y1 = y + (modYRes / 4.0f);