From f2e35db55e960abab4de1f6e1be11e9947de0ea9 Mon Sep 17 00:00:00 2001
From: Paul McCarthy <pauldmccarthy@gmail.com>
Date: Thu, 8 Oct 2020 13:34:53 +0100
Subject: [PATCH] BF: Missing normalisation step in decompose
---
fsl/transform/affine.py | 1 +
1 file changed, 1 insertion(+)
diff --git a/fsl/transform/affine.py b/fsl/transform/affine.py
index c7fb71f2a..57007df77 100644
--- a/fsl/transform/affine.py
+++ b/fsl/transform/affine.py
@@ -214,6 +214,7 @@ def decompose(xform, angles=True):
# The process of finding the scaling factors and shear parameters
# is interleaved. First, find s_x = |M'_1|.
sx = np.sqrt(np.dot(M1, M1))
+ M1 /= sx
# Then, compute an initial value for the xy shear factor,
# s_xy = M'_1 * M'_2. (this is too large by the y scaling factor).
--
GitLab