Changed ().I for np.linalg.pinv to address Singular mat error

spinalcordtoolbox · Jun 24, 2020 · a456fa5 · a456fa5
1 parent 111289a
commit a456fa5
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/spinalcordtoolbox/centerline/nurbs.py b/spinalcordtoolbox/centerline/nurbs.py
@@ -673,9 +673,9 @@ def reconstructGlobalApproximation(self, P_x, P_y, P_z, p, n, w):
             Tz.append(somme)
         Tz = np.matrix(Tz)
 
-        P_xb = (R.T * W * R).I * Tx.T
-        P_yb = (R.T * W * R).I * Ty.T
-        P_zb = (R.T * W * R).I * Tz.T
+        P_xb = np.linalg.pinv(R.T * W * R) * Tx.T
+        P_yb = np.linalg.pinv(R.T * W * R) * Ty.T
+        P_zb = np.linalg.pinv(R.T * W * R) * Tz.T
 
         # Modification of first and last control points
         P_xb[0], P_yb[0], P_zb[0] = P_x[0], P_y[0], P_z[0]

diff --git a/unit_testing/test_centerline.py b/unit_testing/test_centerline.py
@@ -20,7 +20,7 @@
 from sct_utils import init_sct
 
 init_sct(log_level=2)  # Set logger in debug mode
-VERBOSE = 2  # Set to 2 to save images, 0 otherwise
+VERBOSE = 0  # Set to 2 to save images, 0 otherwise
 
 
 # Generate a list of fake centerlines: (dummy_segmentation(params), dict of expected results)