Skip to content

Commit

Permalink
FIX: Drop safeguard errors; manually verified weird references
Browse files Browse the repository at this point in the history
  • Loading branch information
effigies committed Oct 10, 2024
1 parent 23c4636 commit 1607996
Showing 1 changed file with 2 additions and 9 deletions.
11 changes: 2 additions & 9 deletions fmriprep/utils/transforms.py
Original file line number Diff line number Diff line change
Expand Up @@ -49,7 +49,8 @@ def load_ants_h5(filename: Path) -> nt.base.TransformBase:
# * Always return a nitransforms TransformBase
# * Construct warp affine from fixed parameters
#
# This should be upstreamed into nitransforms
# This has been upstreamed into nitransforms.
# Future versions should switch to using nitransforms directly.
h = h5py.File(filename)
xform = ITKCompositeH5.from_h5obj(h)

Expand Down Expand Up @@ -78,14 +79,6 @@ def load_ants_h5(filename: Path) -> nt.base.TransformBase:
spacing = transform2['TransformFixedParameters'][6:9]
direction = transform2['TransformFixedParameters'][9:].reshape((3, 3))

# We are not yet confident that we handle anisotropic spacing
# or direction cosine ordering correctly.
# If we confirm or fix, we can remove these checks.
if not np.allclose(spacing, spacing[0]):
raise ValueError(f'Anisotropic spacing: {spacing}')
if not np.allclose(direction, direction.T):
raise ValueError(f'Asymmetric direction matrix: {direction}')

# ITK uses LPS affines
lps_affine = compose_affine(T=origin, R=direction, Z=spacing)
ras_affine = np.diag([-1, -1, 1, 1]) @ lps_affine
Expand Down

0 comments on commit 1607996

Please sign in to comment.