[BUG] ANTsRegistrator resamples MRI modalities with nearest-neighbor interpolation by default
@LMZimmer is already working on this.
Since Sep 18, 2026.
Assessment
This issue has not been assessed yet.
Description
Version: brainles_preprocessing 0.6.13 (same in 0.6.10)
ANTsRegistrator.transform defaults to interpolator="nearestNeighbor", and ANTsRegistrator.register calls self.transform(...) without passing an interpolator. Modality.register / Modality.transform also never pass one (brainles_preprocessing/modality.py:265, :398, :411).
As a result, every resampling step in AtlasCentricPreprocessor.run (coregistration to the center modality, atlas registration incl. the combined coregistration+atlas transform, atlas correction) resamples the MRI modalities with nearest neighbor.
There currently seems to be no supported way to change this from the user side: passing ANTsRegistrator(transformation_params={"interpolator":"linear"}) fails with "TypeError: apply_transforms() got multiple values for keyword argument 'interpolator'" because transform forwards interpolator=interpolator, **transform_kwargs to ants.apply_transforms (ANTs.py:217). The only workaround is subclassing ANTsRegistrator and overriding the default of transform.
Suggested fix: default interpolator to "linear" in ANTsRegistrator.transform (and inverse_transform), and let Modality.register / Modality.transform pass an interpolator.
Reproducibility
import tempfile
from pathlib import Path
import ants
import numpy as np
from brainles_preprocessing.registration import ANTsRegistrator
tmp = Path(tempfile.mkdtemp())
img = tmp / "img.nii.gz"
ants.image_write(ants.from_numpy(np.random.rand(16, 16, 16).astype(np.float32)), str(img))
tx = tmp / "identity.mat"
ants.write_transform(ants.new_ants_transform(dimension=3, transform_type="AffineTransform"), str(tx))
registrator = ANTsRegistrator(transformation_params={"interpolator": "linear"})
registrator.transform(
fixed_image_path=img,
moving_image_path=img,
transformed_image_path=tmp / "out.nii.gz",
matrix_path=tx,
log_file_path=tmp / "log.txt",
)
- Dominant language
- C
- Stars
- 39
- Forks
- 11
- PR merge metrics
- No merged PRs in 30d
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
More from BrainLesion/preprocessing
-
bug
Difficulty 2/5 1-3 hours Newbie friendliness 68/100
BrainLesion/preprocessing#169 · 7 comments ·
-
bug
Difficulty 3/5 1-2 days Newbie friendliness 35/100
BrainLesion/preprocessing#190 ·
-
bug
BrainLesion/preprocessing#188 · 1 reaction · 2 assignees ·
-
enhancement
Difficulty 4/5 3-5 days Newbie friendliness 35/100
BrainLesion/preprocessing#185 ·
-
bug
Difficulty 3/5 1-2 days Newbie friendliness 58/100
BrainLesion/preprocessing#177 · 3 comments ·
All issues in BrainLesion/preprocessing
Similar issues
-
Difficulty 2/5 1-3 hours Newbie friendliness 65/100
duckdb/duckdb-wasm#2258 ·
-
Difficulty 1/5 Under an hour Newbie friendliness 90/100
rofl0r/microsocks#110 ·
-
Difficulty 2/5 1-3 hours Newbie friendliness 65/100
espressif/esp-aws-iot#261 ·
-
os:linux
Difficulty 2/5 1-3 hours Newbie friendliness 70/100
mpv-player/mpv#18510 · 2 comments ·
-
Difficulty 1/5 Under an hour Newbie friendliness 85/100