def test_syn_registration():
with nbtmp.InTemporaryDirectory() as tmpdir:
warped_moving, mapping = syn_registration(subset_b0,
subset_t2,
moving_affine=hardi_affine,
static_affine=MNI_T2_affine,
step_length=0.1,
metric='CC',
dim=3,
level_iters=[5, 5, 5],
sigma_diff=2.0,
radius=1,
prealign=None)
npt.assert_equal(warped_moving.shape, subset_t2.shape)
mapping_fname = op.join(tmpdir, 'mapping.nii.gz')
write_mapping(mapping, mapping_fname)
file_mapping = read_mapping(mapping_fname, subset_b0_img,
subset_t2_img)
# Test that it has the same effect on the data:
warped_from_file = file_mapping.transform(subset_b0)
npt.assert_equal(warped_from_file, warped_moving)
# Test that it is, attribute by attribute, identical:
for k in mapping.__dict__:
assert (np.all(
mapping.__getattribute__(k) == file_mapping.__getattribute__(
k)))
def test_syn_registration():
with nbtmp.InTemporaryDirectory() as tmpdir:
warped_moving, mapping = syn_registration(subset_b0,
subset_t2,
moving_affine=hardi_affine,
static_affine=MNI_T2_affine,
step_length=0.1,
metric='CC',
dim=3,
level_iters=[10, 10, 5],
sigma_diff=2.0,
prealign=None)
npt.assert_equal(warped_moving.shape, subset_t2.shape)
mapping_fname = op.join(tmpdir, 'mapping.nii.gz')
write_mapping(mapping, mapping_fname)
file_mapping = read_mapping(mapping_fname,
subset_b0_img,
subset_t2_img)
# Test that it has the same effect on the data:
warped_from_file = file_mapping.transform(subset_b0)
npt.assert_equal(warped_from_file, warped_moving)
# Test that it is, attribute by attribute, identical:
for k in mapping.__dict__:
assert (np.all(mapping.__getattribute__(k) ==
file_mapping.__getattribute__(k)))