/usr/bin/nipy_4d_realign is in python-nipy 0.3.0+git262-gbb838d7-1.
This file is owned by root:root, with mode 0o755.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 | #! /usr/bin/python
# emacs: -*- mode: python; py-indent-offset: 4; indent-tabs-mode: nil -*-
# vi: set ft=python sts=4 ts=4 sw=4 et:
"""This is a wrapper of SpaceTimeRealign
Based on:
Alexis Roche (2011) A Four-Dimensional Registration Algorithm With Application to Joint Correction of Motion and Slice Timing in fMRI. IEEE Trans. Med. Imaging 30(8): 1546-1554
"""
import os
import os.path as op
import nipy.algorithms.registration as reg
import nipy.externals.argparse as argparse
parser = argparse.ArgumentParser()
parser.add_argument('TR', type=float, metavar='Float', help="""The TR of the measurement""")
parser.add_argument('input', type=str, metavar='File',
help="""Path to a nifti file, or to a folder containing nifti files. If a path to a folder is provided, the order of motion correction will be np.sort(list_of_files). The outputs will be '*_mc.par' (containing 3 translation and three rotation parameters) and '*_mc.nii.gz' containing the motion corrected data (unless 'apply' is set to False)""")
parser.add_argument('--slice_order', type=str, metavar='String',
help="""The order of slice aquisition {'ascending', 'descending' (default), or the name of a function from `nipy.algorithms.slicetiming.timefuncs`}""", default='descending')
parser.add_argument('--slice_dim', type=int, metavar='Int', help="""Integer
denoting the axis in `images` that is the slice axis. In a 4D image, this will
often be axis = 2 (default).""", default=2)
parser.add_argument('--slice_dir', type=int, metavar='Int', help=""" 1 if the
slices were acquired slice 0 first (default), slice -1 last, or -1 if acquire slice -1 first, slice 0 last.""", default=1)
parser.add_argument('--make_figure', type=bool, metavar='Bool',
help="""Whether to generate a '.png' figure with the motion parameters across runs. {True, False}. Default: False """, default=False)
parser.add_argument('--save_path', type=str, metavar='String',
help="""Full path to a file-system location for the output files. Defaults to the same location as the input files""",
default='none')
parser.add_argument('--save_params', type=bool, metavar='Bool',
help="""Whether to save the motion corrections parameters (3 rotations, 3 translations). {True, False}. Default: False. NOTE: The rotations are not Euler angles, but a rotation vector. Use `nipy.algorithms.registration.to_matrix44` to convert to a 4-by-4 affine matrix""", default=False)
# parse the command line
args = parser.parse_args()
if __name__ == '__main__':
if args.save_path == 'none':
save_path = op.split(args.input)[0]
else:
save_path = args.save_path
xform = reg.space_time_realign(args.input, float(args.TR),
slice_order=args.slice_order,
slice_dim=int(args.slice_dim),
slice_dir=int(args.slice_dir),
apply=True, # We always apply the xform in the cli
make_figure=args.make_figure,
out_name=save_path)
if args.save_params:
f = file(op.join(save_path, 'mc.par'), 'w')
for x in xform:
euler_rot = reg.aff2euler(x.as_affine())
for r in euler_rot:
f.write('%s\t'%r)
for t in x.translation[:-1]:
f.write('%s\t'%t)
f.write('%s\n'%x.translation[-1])
f.close()
|