import numpy as np from nibabel import load, save from nipy.neurospin.image import Image, asNifti1Image, transform_image from nipy.neurospin.registration import IconicRegistration, Affine, GridTransform from os.path import join datadir = join('d:\\', 'home', 'Alexis', 'My Dropbox', 'example_images') I = load(join(datadir, 'b0.nii')) J = load(join(datadir, 'b1000.nii')) # Create a registration instance print('Creating matcher instance') R = IconicRegistration(Image(I), Image(J)) R.set_source_fov(fixed_npoints=64**3) R.similarity = 'cr' R.interp = 'tri' T = R.optimize(Affine())