Now possible with camino ProcStreamlines (convert.py) to collect created probability maps as output.

Author: mick-d

nipype/interfaces/camino/convert.py

@@ -7,10 +7,12 @@
 
 """
 import os
+import glob
 
 from nipype.interfaces.base import (CommandLineInputSpec, CommandLine, traits,
                                     TraitedSpec, File, StdOutCommandLine,
-                                    StdOutCommandLineInputSpec, isdefined)
+                                    OutputMultiPath, StdOutCommandLineInputSpec,
+                                    isdefined)
 from nipype.utils.filemanip import split_filename
 
 class Image2VoxelInputSpec(StdOutCommandLineInputSpec):
@@ -226,16 +228,21 @@ class ProcStreamlinesInputSpec(StdOutCommandLineInputSpec):
     outputtracts = traits.Bool(argstr='-outputtracts', desc="Output streamlines in raw binary format.")
 
     outputroot = File(exists=False, argstr='-outputroot %s',
-                    desc='root directory for output')
+                    desc='Prepended onto all output file names.')
 
     gzip = traits.Bool(argstr='-gzip', desc="save the output image in gzip format")
-    outputcp = traits.Bool(argstr='-outputcp', desc="output the connection probability map (Analyze image, float)")
-    outputsc = traits.Bool(argstr='-outputsc', desc="output the connection probability map (raw streamlines, int)")
-    outputacm = traits.Bool(argstr='-outputacm', desc="output all tracts in a single connection probability map (Analyze image)")
-    outputcbs = traits.Bool(argstr='-outputcbs', desc="outputs connectivity-based segmentation maps; requires target outputfile")
+    outputcp = traits.Bool(argstr='-outputcp', desc="output the connection probability map (Analyze image, float)",
+                           requires=['outputroot','seedfile'])
+    outputsc = traits.Bool(argstr='-outputsc', desc="output the connection probability map (raw streamlines, int)",
+                           requires=['outputroot','seedfile'])
+    outputacm = traits.Bool(argstr='-outputacm', desc="output all tracts in a single connection probability map (Analyze image)",
+                            requires=['outputroot','seedfile'])
+    outputcbs = traits.Bool(argstr='-outputcbs', desc="outputs connectivity-based segmentation maps; requires target outputfile",
+                            requires=['outputroot','targetfile','seedfile'])
 
 class ProcStreamlinesOutputSpec(TraitedSpec):
     proc = File(exists=True, desc='Processed Streamlines')
+    outputroot_files = OutputMultiPath(File(exists=True))
 
 class ProcStreamlines(StdOutCommandLine):
     """
@@ -256,9 +263,29 @@ class ProcStreamlines(StdOutCommandLine):
     input_spec=ProcStreamlinesInputSpec
     output_spec=ProcStreamlinesOutputSpec
 
+    def _run_interface(self, runtime):
+        outputroot = self.inputs.outputroot
+        if isdefined(outputroot):
+            outputroot = os.path.join('procstream_outfiles', outputroot)
+            base, filename, ext = split_filename(outputroot)
+            if not os.path.exists(base):
+                os.makedirs(base)
+            self.inputs.outputroot = outputroot
+            new_runtime = super(ProcStreamlines, self)._run_interface(runtime)
+            self.outputroot_files = self._get_ouputroot_files(outputroot)
+            return new_runtime
+        else:
+            new_runtime = super(ProcStreamlines, self)._run_interface(runtime)
+            return new_runtime
+
+    def _get_ouputroot_files(self, outputroot):
+        outputroot_files = glob.glob(os.path.join(os.getcwd(),outputroot+'*'))
+        return outputroot_files
+
     def _list_outputs(self):
         outputs = self.output_spec().get()
         outputs['proc'] = os.path.abspath(self._gen_outfilename())
+        outputs['outputroot_files'] = self.outputroot_files
         return outputs
 
     def _gen_outfilename(self):

nipype/interfaces/camino/tests/test_auto_ProcStreamlines.py

@@ -56,14 +56,18 @@ def test_ProcStreamlines_inputs():
     position=-1,
     ),
     outputacm=dict(argstr='-outputacm',
+    requires=['outputroot', 'seedfile'],
     ),
     outputcbs=dict(argstr='-outputcbs',
+    requires=['outputroot', 'targetfile', 'seedfile'],
     ),
     outputcp=dict(argstr='-outputcp',
+    requires=['outputroot', 'seedfile'],
     ),
     outputroot=dict(argstr='-outputroot %s',
     ),
     outputsc=dict(argstr='-outputsc',
+    requires=['outputroot', 'seedfile'],
     ),
     outputtracts=dict(argstr='-outputtracts',
     ),
@@ -103,7 +107,8 @@ def test_ProcStreamlines_inputs():
             yield assert_equal, getattr(inputs.traits()[key], metakey), value
 
 def test_ProcStreamlines_outputs():
-    output_map = dict(proc=dict(),
+    output_map = dict(outputroot_files=dict(),
+    proc=dict(),
     )
     outputs = ProcStreamlines.output_spec()