Merge pull request #1227 from effigies/enh/catch_svd_failure

effigies · web-flow · commit 977d044ba056 · 2024-09-05T18:15:13.000-04:00
ENH: Catch SVD failure and raise informative HeaderDataError
diff --git a/nibabel/nifti1.py b/nibabel/nifti1.py
@@ -1097,7 +1097,10 @@ def set_qform(self, affine, code=None, strip_shears=True):
         # (a subtle requirement of the NIFTI format qform transform)
         # Transform below is polar decomposition, returning the closest
         # orthogonal matrix PR, to input R
-        P, S, Qs = npl.svd(R)
+        try:
+            P, S, Qs = npl.svd(R)
+        except np.linalg.LinAlgError as e:
+            raise HeaderDataError(f'Could not decompose affine:\n{affine}') from e
         PR = np.dot(P, Qs)
         if not strip_shears and not np.allclose(PR, R):
             raise HeaderDataError('Shears in affine and `strip_shears` is False')
