|
56 | 56 | import nipype.interfaces.cmtk as cmtk |
57 | 57 | import nipype.interfaces.dipy as dipy |
58 | 58 | import inspect |
59 | | -import os.path as op # system functions |
| 59 | +import os, os.path as op # system functions |
60 | 60 | from nipype.workflows.dmri.fsl.dti import create_eddy_correct_pipeline |
61 | 61 | from nipype.workflows.dmri.camino.connectivity_mapping import select_aparc_annot |
62 | 62 | from nipype.utils.misc import package_check |
63 | 63 | import warnings |
64 | 64 | from nipype.workflows.dmri.connectivity.nx import create_networkx_pipeline, create_cmats_to_csv_pipeline |
| 65 | +from nipype.workflows.smri.freesurfer import create_tessellation_flow |
65 | 66 |
|
66 | 67 | try: |
67 | 68 | package_check('cmp') |
|
82 | 83 | fs.FSCommand.set_default_subjects_dir(subjects_dir) |
83 | 84 | fsl.FSLCommand.set_default_output_type('NIFTI') |
84 | 85 |
|
| 86 | +fs_dir = os.environ['FREESURFER_HOME'] |
| 87 | +lookup_file = op.join(fs_dir,'FreeSurferColorLUT.txt') |
| 88 | + |
85 | 89 | """ |
86 | 90 | This needs to point to the fdt folder you can find after extracting |
87 | 91 |
|
|
328 | 332 |
|
329 | 333 | CFFConverter = pe.Node(interface=cmtk.CFFConverter(), name="CFFConverter") |
330 | 334 | CFFConverter.inputs.script_files = op.abspath(inspect.getfile(inspect.currentframe())) |
331 | | -giftiSurfaces = pe.Node(interface=util.Merge(8), name="GiftiSurfaces") |
| 335 | +giftiSurfaces = pe.Node(interface=util.Merge(9), name="GiftiSurfaces") |
332 | 336 | giftiLabels = pe.Node(interface=util.Merge(2), name="GiftiLabels") |
333 | 337 | niftiVolumes = pe.Node(interface=util.Merge(3), name="NiftiVolumes") |
334 | 338 | fiberDataArrays = pe.Node(interface=util.Merge(4), name="FiberDataArrays") |
|
344 | 348 | NxStatsCFFConverter = pe.Node(interface=cmtk.CFFConverter(), name="NxStatsCFFConverter") |
345 | 349 | NxStatsCFFConverter.inputs.script_files = op.abspath(inspect.getfile(inspect.currentframe())) |
346 | 350 |
|
| 351 | +tessflow = create_tessellation_flow(name='tessflow', out_format='gii') |
| 352 | +tessflow.inputs.inputspec.lookup_file = lookup_file |
| 353 | + |
347 | 354 | """ |
348 | 355 | Connecting the workflow |
349 | 356 | ======================= |
|
371 | 378 | mapping.connect([(inputnode, FreeSurferSourceRH,[("subjects_dir","subjects_dir")])]) |
372 | 379 | mapping.connect([(inputnode, FreeSurferSourceRH,[("subject_id","subject_id")])]) |
373 | 380 |
|
| 381 | +mapping.connect([(inputnode, tessflow,[("subjects_dir","inputspec.subjects_dir")])]) |
| 382 | +mapping.connect([(inputnode, tessflow,[("subject_id","inputspec.subject_id")])]) |
| 383 | + |
374 | 384 | mapping.connect([(inputnode, parcellate,[("subjects_dir","subjects_dir")])]) |
375 | 385 | mapping.connect([(inputnode, parcellate,[("subject_id","subject_id")])]) |
376 | 386 | mapping.connect([(parcellate, mri_convert_ROI_scale500,[('roi_file','in_file')])]) |
|
516 | 526 | mapping.connect([(mris_convertRHinflated, giftiSurfaces,[("converted","in6")])]) |
517 | 527 | mapping.connect([(mris_convertLHsphere, giftiSurfaces,[("converted","in7")])]) |
518 | 528 | mapping.connect([(mris_convertRHsphere, giftiSurfaces,[("converted","in8")])]) |
| 529 | +mapping.connect([(tessflow, giftiSurfaces,[("outputspec.meshes","in9")])]) |
519 | 530 |
|
520 | 531 | mapping.connect([(mris_convertLHlabels, giftiLabels,[("converted","in1")])]) |
521 | 532 | mapping.connect([(mris_convertRHlabels, giftiLabels,[("converted","in2")])]) |
|
0 commit comments