Merge pull request #4 from sgiavasis/0.9.2_cpac-tested

0.9.2 cpac tested

Author: stymy

Vagrantfile

@@ -20,19 +20,19 @@ $script = <<SCRIPT
 # qconf -aattr queue slots "2, [neuro=3]" main.q
 
 # install anaconda
-wget http://repo.continuum.io/miniconda/Miniconda-3.0.0-Linux-x86_64.sh -O miniconda.sh
+wget http://repo.continuum.io/miniconda/Miniconda-latest-Linux-x86_64.sh -O miniconda.sh
 chmod +x miniconda.sh
 ./miniconda.sh -b
 echo "export PATH=$HOME/miniconda/bin:\\$PATH" >> .bashrc
 
 # install nipype dependencies
-$HOME/anaconda/bin/conda update --yes conda
-$HOME/anaconda/bin/conda install --yes pip numpy scipy nose traits networkx
-$HOME/anaconda/bin/conda install --yes dateutil ipython-notebook matplotlib
-$HOME/anaconda/bin/pip install nibabel --use-mirrors
-$HOME/anaconda/bin/pip install https://github.com/RDFLib/rdflib/archive/master.zip
-$HOME/anaconda/bin/pip install https://github.com/satra/prov/archive/enh/rdf.zip
-$HOME/anaconda/bin/pip install https://github.com/nipy/nipype/archive/master.zip
+$HOME/miniconda/bin/conda update --yes conda
+$HOME/miniconda/bin/conda install --yes pip numpy scipy nose traits networkx
+$HOME/miniconda/bin/conda install --yes dateutil ipython-notebook matplotlib
+$HOME/miniconda/bin/pip install nibabel --use-mirrors
+$HOME/miniconda/bin/pip install https://github.com/RDFLib/rdflib/archive/master.zip
+$HOME/miniconda/bin/pip install https://github.com/satra/prov/archive/enh/rdf.zip
+$HOME/miniconda/bin/pip install https://github.com/nipy/nipype/archive/master.zip
 SCRIPT
 
 Vagrant.configure(VAGRANTFILE_API_VERSION) do |config|
@@ -42,7 +42,7 @@ Vagrant.configure(VAGRANTFILE_API_VERSION) do |config|
     engine_config.vm.box = "gridneuro"
     #engine_config.vm.box_url = "http://files.vagrantup.com/precise64.box"
     engine_config.vm.box_url = "https://dl.dropboxusercontent.com/u/363467/precise64_neuro.box"
-    engine_config.vm.network :forwarded_port, guest: 80, host: 8080
+    #engine_config.vm.network :forwarded_port, guest: 80, host: 8080
 
     #engine_config.vm.network :public_network, :bridge => 'en0: Wi-Fi (AirPort)'
     engine_config.vm.network :private_network, ip: "192.168.100.20"

nipype/info.py

@@ -104,6 +104,8 @@ def get_nipype_gitversion():
 NUMPY_MIN_VERSION = '1.3'
 SCIPY_MIN_VERSION = '0.7'
 TRAITS_MIN_VERSION = '4.0'
+DATEUTIL_MIN_VERSION = '1.0'
+NOSE_MIN_VERSION = '1.0'
 
 NAME                = 'nipype'
 MAINTAINER          = "nipype developers"
@@ -122,6 +124,7 @@ def get_nipype_gitversion():
 MICRO               = _version_micro
 ISRELEASE           = _version_extra == ''
 VERSION             = __version__
-REQUIRES            = ["nibabel (>=1.0)", "networkx (>=1.0)", "numpy (>=1.3)",
-                       "scipy (>=0.7)", "traits (>=4.0)"]
+REQUIRES            = ["nibabel>=1.0", "networkx>=1.0", "numpy>=1.3",
+                       "python-dateutil>1.0", "scipy>=0.7", "traits>=4.0",
+                       "nose>=1.0"]
 STATUS              = 'stable'

nipype/interfaces/afni/base.py

@@ -144,7 +144,7 @@ def set_default_output_type(cls, outputtype):
         else:
             raise AttributeError('Invalid AFNI outputtype: %s' % outputtype)
 
-    def _overload_extension(self, value):
+    def _overload_extension(self, value, name=None):
         path, base, _ = split_filename(value)
         return os.path.join(path, base + Info.outputtype_to_ext(self.inputs.outputtype))
 

nipype/interfaces/afni/preprocess.py

@@ -560,7 +560,6 @@ class VolregInputSpec(AFNICommandInputSpec):
                    copyfile=False)
     out_file = File(name_template="%s_volreg", desc='output image file name',
                     argstr='-prefix %s', name_source="in_file")
-
     basefile = File(desc='base file for registration',
                     argstr='-base %s',
                     position=-6,
@@ -582,12 +581,18 @@ class VolregInputSpec(AFNICommandInputSpec):
                             argstr='-tshift 0')
     copyorigin = traits.Bool(desc='copy base file origin coords to output',
                              argstr='-twodup')
+    oned_matrix_save = File(name_template='%s.aff12.1D',
+                            desc='Save the matrix transformation',
+                            argstr='-1Dmatrix_save %s',
+                            keep_extension=True,
+                            name_source="in_file")
 
 
 class VolregOutputSpec(TraitedSpec):
     out_file = File(desc='registered file', exists=True)
     md1d_file = File(desc='max displacement info file', exists=True)
     oned_file = File(desc='movement parameters info file', exists=True)
+    oned_matrix_save = File(desc='matrix transformation from base to input', exists=True)
 
 
 class Volreg(AFNICommand):

nipype/interfaces/ants/__init__.py

@@ -11,7 +11,7 @@
 
 
 # Segmentation Programs
-from .segmentation import Atropos, N4BiasFieldCorrection
+from .segmentation import Atropos, LaplacianThickness, N4BiasFieldCorrection
 
 # Utility Programs
 from .utils import AverageAffineTransform, AverageImages, MultiplyImages, JacobianDeterminant

nipype/interfaces/ants/segmentation.py

@@ -160,6 +160,67 @@ def _list_outputs(self):
                 outputs['posteriors'].append(os.path.abspath(self.inputs.output_posteriors_name_template % (i + 1)))
         return outputs
 
+class LaplacianThicknessInputSpec(ANTSCommandInputSpec):
+    input_wm = File(argstr='%s', mandatory=True, copyfile=True,
+                       desc=('white matter segmentation image'),
+                       position=1)
+    input_gm = File(argstr='%s', mandatory=True, copyfile=True,
+                       desc=('gray matter segmentation image'),
+                       position=2)
+    output_image = File(desc='name of output file', argstr='%s', position=3,
+                        genfile=True, hash_files=False)
+    smooth_param = traits.Float(argstr='smoothparam=%d', desc='', position=4)
+    prior_thickness = traits.Float(argstr='priorthickval=%d', desc='',
+                                   position=5)
+    dT = traits.Float(argstr='dT=%d', desc='', position=6)
+    sulcus_prior = traits.Bool(argstr='use-sulcus-prior', desc='', position=7)
+    opt_tolerance = traits.Float(argstr='optional-laplacian-tolerance=%d', 
+                                 desc='', position=8)
+
+
+class LaplacianThicknessOutputSpec(TraitedSpec):
+    output_image = File(exists=True, desc='Cortical thickness')
+
+
+class LaplacianThickness(ANTSCommand):
+    """Calculates the cortical thickness from an anatomical image
+
+    Examples
+    --------
+
+    >>> from nipype.interfaces.ants import LaplacianThickness
+    >>> cort_thick = LaplacianThickness()
+    >>> cort_thick.inputs.input_wm = 'white_matter.nii.gz'
+    >>> cort_thick.inputs.input_gm = 'gray_matter.nii.gz'
+    >>> cort_thick.inputs.output_image = 'output_thickness.nii.gz'
+    >>> cort_thick.cmdline
+    'LaplacianThickness white_matter.nii.gz gray_matter.nii.gz output_thickness.nii.gz'
+
+    """
+
+    _cmd = 'LaplacianThickness'
+    input_spec = LaplacianThicknessInputSpec
+    output_spec = LaplacianThicknessOutputSpec
+
+    def _gen_filename(self, name):
+        if name == 'output_image':
+            output = self.inputs.output_image
+            if not isdefined(output):
+                _, name, ext = split_filename(self.inputs.input_wm)
+                output = name + '_thickness' + ext
+            return output
+        return None
+
+    def _list_outputs(self):
+        outputs = self._outputs().get()
+        _, name, ext = split_filename(os.path.abspath(self.inputs.input_wm))
+        outputs['output_image'] = os.path.join(os.getcwd(),
+                                               ''.join((name,
+                                                        self.inputs.output_image,
+                                                        ext)))
+        return outputs
+
+
 class N4BiasFieldCorrectionInputSpec(ANTSCommandInputSpec):
     dimension = traits.Enum(3, 2, argstr='--image-dimension %d',
                             usedefault=True,
@@ -228,4 +289,4 @@ def _gen_filename(self, name):
     def _list_outputs(self):
         outputs = self._outputs().get()
         outputs['output_image'] = os.path.abspath(self._gen_filename('output_image'))
-        return outputs
+        return outputs

nipype/interfaces/base.py

@@ -731,6 +731,7 @@ def _get_trait_desc(self, inputs, name, spec):
         desc = spec.desc
         xor = spec.xor
         requires = spec.requires
+        argstr = spec.argstr
 
         manhelpstr = ['\t%s' % name]
 
@@ -753,6 +754,17 @@ def _get_trait_desc(self, inputs, name, spec):
                                    initial_indent='\t\t',
                                    subsequent_indent='\t\t')
 
+        if argstr:
+            pos = spec.position
+            if pos is not None:
+                manhelpstr += wrap('flag: %s, position: %s' % (argstr, pos), 70,
+                                   initial_indent='\t\t',
+                                   subsequent_indent='\t\t')
+            else:
+                manhelpstr += wrap('flag: %s' % argstr, 70,
+                                   initial_indent='\t\t',
+                                   subsequent_indent='\t\t')
+
         if xor:
             line = '%s' % ', '.join(xor)
             manhelpstr += wrap(line, 70,
@@ -1480,13 +1492,13 @@ def _filename_from_source(self, name):
             _, _, ext = split_filename(retval)
             if trait_spec.keep_extension and ext:
                 return retval
-            return self._overload_extension(retval)
+            return self._overload_extension(retval, name)
         return retval
 
     def _gen_filename(self, name):
         raise NotImplementedError
 
-    def _overload_extension(self, value):
+    def _overload_extension(self, value, name=None):
         return value
 
     def _list_outputs(self):

nipype/interfaces/camino/__init__.py

@@ -13,3 +13,4 @@
                   ComputeTensorTrace,  ComputeEigensystem, DTMetric)
 from .calib import (SFPICOCalibData, SFLUTGen)
 from .odf import (QBallMX, LinRecon, SFPeaks, MESD)
+from .utils import ImageStats

nipype/interfaces/camino/tests/test_auto_ImageStats.py

@@ -0,0 +1,46 @@
+# AUTO-GENERATED by tools/checkspecs.py - DO NOT EDIT
+from nipype.testing import assert_equal
+from nipype.interfaces.camino.utils import ImageStats
+
+def test_ImageStats_inputs():
+    input_map = dict(args=dict(argstr='%s',
+    ),
+    environ=dict(nohash=True,
+    usedefault=True,
+    ),
+    ignore_exception=dict(nohash=True,
+    usedefault=True,
+    ),
+    in_files=dict(argstr='-images %s',
+    mandatory=True,
+    position=-1,
+    ),
+    out_type=dict(argstr='-outputdatatype %s',
+    usedefault=True,
+    ),
+    output_root=dict(argstr='-outputroot %s',
+    mandatory=True,
+    ),
+    stat=dict(argstr='-stat %s',
+    mandatory=True,
+    units='NA',
+    ),
+    terminal_output=dict(mandatory=True,
+    nohash=True,
+    ),
+    )
+    inputs = ImageStats.input_spec()
+
+    for key, metadata in input_map.items():
+        for metakey, value in metadata.items():
+            yield assert_equal, getattr(inputs.traits()[key], metakey), value
+
+def test_ImageStats_outputs():
+    output_map = dict(out_file=dict(),
+    )
+    outputs = ImageStats.output_spec()
+
+    for key, metadata in output_map.items():
+        for metakey, value in metadata.items():
+            yield assert_equal, getattr(outputs.traits()[key], metakey), value
+

nipype/interfaces/camino/utils.py

@@ -0,0 +1,63 @@
+"""
+    Change directory to provide relative paths for doctests
+    >>> import os
+    >>> filepath = os.path.dirname( os.path.realpath( __file__ ) )
+    >>> datadir = os.path.realpath(os.path.join(filepath, '../../testing/data'))
+    >>> os.chdir(datadir)
+
+"""
+from nipype.interfaces.base import (traits, TraitedSpec, File,
+                                    CommandLine, CommandLineInputSpec, isdefined,
+                                    InputMultiPath)
+from nipype.utils.filemanip import split_filename
+import os
+
+class ImageStatsInputSpec(CommandLineInputSpec):
+    in_files = InputMultiPath(File(exists=True), argstr='-images %s', mandatory=True,
+                              position=-1, desc=('List of images to process. They must '
+                                                 'be in the same space and have the same '
+                                                 'dimensions.'))
+    stat = traits.Enum("min", "max", "mean", "median", "sum", "std", "var",
+                       argstr='-stat %s', units='NA', mandatory=True,
+                       desc="The statistic to compute.")
+
+    out_type = traits.Enum("float", "char", "short", "int", "long", "double",
+                           argstr='-outputdatatype %s', usedefault=True,
+                           desc=('A Camino data type string, default is "float". '
+                                 'Type must be signed.'))
+    output_root = File(argstr='-outputroot %s', mandatory=True,
+                       desc=('Filename root prepended onto the names of the output '
+                             ' files. The extension will be determined from the input.'))
+
+class ImageStatsOutputSpec(TraitedSpec):
+    out_file = File(exists=True, desc='Path of the file computed with the statistic chosen')
+
+class ImageStats(CommandLine):
+    """
+    This program computes voxelwise statistics on a series of 3D images. The images
+    must be in the same space; the operation is performed voxelwise and one output
+    is produced per voxel.
+
+    Examples
+    --------
+
+    >>> import nipype.interfaces.camino as cam
+    >>> imstats = cam.ImageStats()
+    >>> imstats.inputs.in_files = ['im1.nii','im2.nii','im3.nii']
+    >>> imstats.inputs.stat = 'max'
+    >>> imstats.run()                  # doctest: +SKIP
+    """
+    _cmd = 'imagestats'
+    input_spec = ImageStatsInputSpec
+    output_spec = ImageStatsOutputSpec
+
+    def _list_outputs(self):
+        outputs = self.output_spec().get()
+        outputs['out_file'] = os.path.abspath(self._gen_outfilename())
+        return outputs
+
+    def _gen_outfilename(self):
+        output_root = self.inputs.output_root
+        first_file = self.inputs.in_files[0]
+        _, _ , ext = split_filename(first_file)
+        return output_root + ext

nipype/interfaces/fsl/base.py

@@ -239,7 +239,7 @@ def _gen_fname(self, basename, cwd=None, suffix=None, change_ext=True,
                                 use_ext=False, newpath=cwd)
         return fname
 
-    def _overload_extension(self, value):
+    def _overload_extension(self, value, name=None):
         return value + Info.output_type_to_ext(self.inputs.output_type)