Hi,

 

I keep getting the same error every time I run fmriprep on subjects from the OASIS3 dataset.  There are 5-7 T1 scans over 3-4 time points for most subjects, and the registration is failing every time with the return code -9.  I can’t seem to find any indication of what that return code means, so I don’t know what to try to dix.  Does anybody have an idea?


Thanks,

Blaise

 

 

cat /n/scratch2/bdf1/OASIS/derivatives/fmriprep/sub-OAS30004/log/20190719-151129_2b54c614-de57-4bd6-9c95-9c721b09df7d/crash-20190719-153017-bdf1-t1_merge-59cde0ed-a554-4c50-8cc7-fd158f28f7ae.txt

Node: fmriprep_wf.single_subject_OAS30004_wf.anat_preproc_wf.anat_template_wf.t1_merge

Working directory: /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/t1_merge

 

Node inputs:

 

args = <undefined>

auto_detect_sensitivity = True

average_metric = <undefined>

environ = {'SUBJECTS_DIR': '/home/bdf1/software/freesurfer/subjects', 'OMP_NUM_THREADS': '4'}

fixed_timepoint = True

in_files = ['/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz', '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct1/sub-OAS30004_ses-d1101_run-02_T1w_corrected.nii.gz', '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct2/sub-OAS30004_ses-d2229_run-01_T1w_corrected.nii.gz', '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct3/sub-OAS30004_ses-d2229_run-02_T1w_corrected.nii.gz', '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct4/sub-OAS30004_ses-d2232_run-01_T1w_ras_corrected.nii.gz', '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct5/sub-OAS30004_ses-d2232_run-02_T1w_corrected.nii.gz', '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct6/sub-OAS30004_ses-d3457_T1w_ras_corrected.nii.gz']

in_intensity_scales = <undefined>

initial_timepoint = 1

initial_transforms = <undefined>

intensity_scaling = True

no_iteration = True

num_threads = 7

out_file = sub-OAS30004_ses-d1101_run-01_T1w_template.nii.gz

outlier_sensitivity = <undefined>

scaled_intensity_outputs = <undefined>

subjects_dir = /home/bdf1/software/freesurfer/subjects

subsample_threshold = 200

transform_outputs = True

 

Traceback (most recent call last):

  File "/home/bdf1/software/anaconda3/envs/fmriprep141sing/lib/python3.7/site-packages/nipype/pipeline/plugins/multiproc.py", line 69, in run_node

    result['result'] = node.run(updatehash=updatehash)

  File "/home/bdf1/software/anaconda3/envs/fmriprep141sing/lib/python3.7/site-packages/nipype/pipeline/engine/nodes.py", line 472, in run

    result = self._run_interface(execute=True)

  File "/home/bdf1/software/anaconda3/envs/fmriprep141sing/lib/python3.7/site-packages/nipype/pipeline/engine/nodes.py", line 563, in _run_interface

    return self._run_command(execute)

  File "/home/bdf1/software/anaconda3/envs/fmriprep141sing/lib/python3.7/site-packages/nipype/pipeline/engine/nodes.py", line 643, in _run_command

    result = self._interface.run(cwd=outdir)

  File "/home/bdf1/software/anaconda3/envs/fmriprep141sing/lib/python3.7/site-packages/nipype/interfaces/freesurfer/base.py", line 265, in run

    return super(FSCommandOpenMP, self).run(**inputs)

  File "/home/bdf1/software/anaconda3/envs/fmriprep141sing/lib/python3.7/site-packages/nipype/interfaces/freesurfer/base.py", line 144, in run

    return super(FSCommand, self).run(**inputs)

  File "/home/bdf1/software/anaconda3/envs/fmriprep141sing/lib/python3.7/site-packages/nipype/interfaces/base/core.py", line 375, in run

    runtime = self._run_interface(runtime)

  File "/home/bdf1/software/anaconda3/envs/fmriprep141sing/lib/python3.7/site-packages/nipype/interfaces/base/core.py", line 758, in _run_interface

    self.raise_exception(runtime)

  File "/home/bdf1/software/anaconda3/envs/fmriprep141sing/lib/python3.7/site-packages/nipype/interfaces/base/core.py", line 695, in raise_exception

    ).format(**runtime.dictcopy()))

RuntimeError: Command:

mri_robust_template --satit --fixtp --mov /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct1/sub-OAS30004_ses-d1101_run-02_T1w_corrected.nii.gz /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct2/sub-OAS30004_ses-d2229_run-01_T1w_corrected.nii.gz /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct3/sub-OAS30004_ses-d2229_run-02_T1w_corrected.nii.gz /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct4/sub-OAS30004_ses-d2232_run-01_T1w_ras_corrected.nii.gz /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct5/sub-OAS30004_ses-d2232_run-02_T1w_corrected.nii.gz /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct6/sub-OAS30004_ses-d3457_T1w_ras_corrected.nii.gz --inittp 1 --iscale --noit --template sub-OAS30004_ses-d1101_run-01_T1w_template.nii.gz --subsample 200 --lta /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/t1_merge/tp1.lta /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/t1_merge/tp2.lta /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/t1_merge/tp3.lta /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/t1_merge/tp4.lta /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/t1_merge/tp5.lta /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/t1_merge/tp6.lta /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/t1_merge/tp7.lta

Standard output:

$Id: mri_robust_template.cpp,v 1.54 2016/05/05 21:17:08 mreuter Exp $

 

--satit: Will estimate SAT iteratively!

--fixtp: Will map everything to init TP!

--mov: Using /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz as movable/source volume.

--mov: Using /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct1/sub-OAS30004_ses-d1101_run-02_T1w_corrected.nii.gz as movable/source volume.

--mov: Using /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct2/sub-OAS30004_ses-d2229_run-01_T1w_corrected.nii.gz as movable/source volume.

--mov: Using /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct3/sub-OAS30004_ses-d2229_run-02_T1w_corrected.nii.gz as movable/source volume.

--mov: Using /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct4/sub-OAS30004_ses-d2232_run-01_T1w_ras_corrected.nii.gz as movable/source volume.

--mov: Using /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct5/sub-OAS30004_ses-d2232_run-02_T1w_corrected.nii.gz as movable/source volume.

--mov: Using /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct6/sub-OAS30004_ses-d3457_T1w_ras_corrected.nii.gz as movable/source volume.

    Total: 7 input volumes

--inittp: Using TP 1 as target for initialization

--iscale: Enableing intensity scaling!

--noit: Will output only first template (no iterations)!

--template: Using sub-OAS30004_ses-d1101_run-01_T1w_template.nii.gz as template output volume.

--subsample: Will subsample if size is larger than 200 on all axes!

--lta: Will output LTA transforms

Setting iscale ...

reading source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz'...

converting source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz' to bspline ...

MRItoBSpline degree 3

reading source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct1/sub-OAS30004_ses-d1101_run-02_T1w_corrected.nii.gz'...

converting source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct1/sub-OAS30004_ses-d1101_run-02_T1w_corrected.nii.gz' to bspline ...

MRItoBSpline degree 3

reading source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct2/sub-OAS30004_ses-d2229_run-01_T1w_corrected.nii.gz'...

converting source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct2/sub-OAS30004_ses-d2229_run-01_T1w_corrected.nii.gz' to bspline ...

MRItoBSpline degree 3

reading source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct3/sub-OAS30004_ses-d2229_run-02_T1w_corrected.nii.gz'...

converting source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct3/sub-OAS30004_ses-d2229_run-02_T1w_corrected.nii.gz' to bspline ...

MRItoBSpline degree 3

reading source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct4/sub-OAS30004_ses-d2232_run-01_T1w_ras_corrected.nii.gz'...

converting source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct4/sub-OAS30004_ses-d2232_run-01_T1w_ras_corrected.nii.gz' to bspline ...

MRItoBSpline degree 3

reading source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct5/sub-OAS30004_ses-d2232_run-02_T1w_corrected.nii.gz'...

converting source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct5/sub-OAS30004_ses-d2232_run-02_T1w_corrected.nii.gz' to bspline ...

MRItoBSpline degree 3

reading source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct6/sub-OAS30004_ses-d3457_T1w_ras_corrected.nii.gz'...

converting source '/n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct6/sub-OAS30004_ses-d3457_T1w_ras_corrected.nii.gz' to bspline ...

MRItoBSpline degree 3

 

MultiRegistration::initializing Xforms (init 1 , maxres 0 , iterate 5 , epsit 0.01 ) :

 

[init] ========================= TP 5 to TP 1 ==============================

         Register TP 5 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct4/sub-OAS30004_ses-d2232_run-01_T1w_ras_corrected.nii.gz )

          to      TP 1 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz )

 

 

[init] ========================= TP 6 to TP 1 ==============================

        Register TP 6 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct5/sub-OAS30004_ses-d2232_run-02_T1w_corrected.nii.gz )

          to      TP 1 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz )

 

 

[init] ========================= TP 3 to TP 1 ==============================

         Register TP 3 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct2/sub-OAS30004_ses-d2229_run-01_T1w_corrected.nii.gz )

          to      TP 1 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz )

 

 

[init] ========================= TP 4 to TP 1 ==============================

         Register TP 4 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct3/sub-OAS30004_ses-d2229_run-02_T1w_corrected.nii.gz )

          to      TP 1 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz )

 

 

[init] ========================= TP 7 to TP 1 ==============================

         Register TP 7 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct6/sub-OAS30004_ses-d3457_T1w_ras_corrected.nii.gz )

          to      TP 1 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz )

 

 

[init] ========================= TP 2 to TP 1 ==============================

         Register TP 2 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct1/sub-OAS30004_ses-d1101_run-02_T1w_corrected.nii.gz )

          to      TP 1 ( /n/scratch2/bdf1/OASIS/work/fmriprep_wf/single_subject_OAS30004_wf/anat_preproc_wf/anat_template_wf/n4_correct/mapflow/_n4_correct0/sub-OAS30004_ses-d1101_run-01_T1w_corrected.nii.gz )

 

       -- Original : (0.999994, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (0.999996, 0.999996, 0.999996)mm and (176, 257, 257) voxels.

       -- Reslicing using cubic bspline

       -- Original : (1, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (1, 1, 1)mm and (177, 256, 256) voxels.

       -- Reslicing using cubic bspline

       -- Original : (0.999994, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (0.999996, 0.999996, 0.999996)mm and (176, 257, 257) voxels.

       -- Reslicing using cubic bspline

       -- Original : (1, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (1, 1, 1)mm and (177, 256, 256) voxels.

       -- Reslicing using cubic bspline

       -- Original : (0.999994, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (0.999996, 0.999996, 0.999996)mm and (176, 257, 257) voxels.

       -- Reslicing using cubic bspline

       -- Original : (0.999994, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (0.999996, 0.999996, 0.999996)mm and (176, 257, 257) voxels.

       -- Reslicing using cubic bspline

MRItoBSpline degree 3

MRItoBSpline degree 3

MRItoBSpline degree 3

MRItoBSpline degree 3

MRItoBSpline degree 3

MRItoBSpline degree 3

       -- Original : (0.999996, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (1, 1, 1)mm and (177, 256, 256) voxels.

       -- Reslicing using cubic bspline

MRItoBSpline degree 3

       -- Original : (0.999996, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (0.999996, 0.999996, 0.999996)mm and (176, 257, 257) voxels.

       -- Reslicing using cubic bspline

MRItoBSpline degree 3

       -- Original : (0.999996, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (0.999996, 0.999996, 0.999996)mm and (176, 257, 257) voxels.

       -- Reslicing using cubic bspline

MRItoBSpline degree 3

       -- Original : (0.999996, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (0.999996, 0.999996, 0.999996)mm and (176, 257, 257) voxels.

       -- Reslicing using cubic bspline

MRItoBSpline degree 3

       -- Original : (0.999996, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (0.999996, 0.999996, 0.999996)mm and (176, 257, 257) voxels.

       -- Reslicing using cubic bspline

       -- Original : (0.999996, 1, 1)mm and (176, 256, 256) voxels.

       -- Resampled: (1, 1, 1)mm and (177, 256, 256) voxels.

       -- Reslicing using cubic bspline

MRItoBSpline degree 3

MRItoBSpline degree 3

 

   - Max Resolution used: 2

     -- gpS ( 44 , 64 , 64 )

     -- gpT ( 44 , 64 , 64 )

   - running loop to estimate saturation parameter:

 

   - Max Resolution used: 2

     -- gpS ( 44 , 64 , 64 )

     -- gpT ( 44 , 64 , 64 )

   - running loop to estimate saturation parameter:

 

   - Max Resolution used: 2

     -- gpS ( 44 , 64 , 64 )

     -- gpT ( 44 , 64 , 64 )

   - running loop to estimate saturation parameter:

 

   - Max Resolution used: 2

     -- gpS ( 44 , 64 , 64 )

     -- gpT ( 44 , 64 , 64 )

   - running loop to estimate saturation parameter:

 

   - Max Resolution used: 2

     -- gpS ( 44 , 64 , 64 )

     -- gpT ( 44 , 64 , 64 )

   - running loop to estimate saturation parameter:

 

   - Max Resolution used: 2

     -- gpS ( 44 , 64 , 64 )

     -- gpT ( 44 , 64 , 64 )

Standard error:

 

Return code: -9

 

 

-------------------------------------------------------
Blaise Frederick
Associate Professor of Psychiatry/Biophysicist
Harvard Medical School/McLean Hospital
bbfrederick@mclean.harvard.edu

http://www.nirs-fmri.net