Hello everybody,

first of all, thank you very much for the great work with fcon 1000 and the preprocessing scripts.
I would like to run the scripts on my own, similar dataset to compare it with some of your datasets.

the scripts worked fine on some of your datasets, so everything seems to be installed correctly.
Now when I try to run them on my data, an exception is being thrown as "global.1D" does not have enough data.

I checked the forum and already fixed the 3dskullstrip debian incompatibility and told it to use original voxel size
(in total added flags -no_use_edge -orig_vol to 3dskullstrip in preprocess_anat.sh).

Still, I get this error. My full protocol is written below.

Thank you so much for any suggestions!

Andy


preprocessing 12
--------------------------------------
!!!! PREPROCESSING ANATOMICAL SCAN!!!!
--------------------------------------
deobliquing 12 anatomical
++ 3drefit: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: RW Cox
++ Processing AFNI dataset mprage.nii.gz
 + loading and re-writing entire dataset mprage.nii.gz
++ 3drefit processed 1 datasets
Reorienting 12 anatomical
skull stripping 12 anatomical
++ 3dcalc: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: A cast of thousands
++ Output dataset ./mprage_brain.nii.gz
---------------------------------------
!!!! PREPROCESSING FUNCTIONAL SCAN !!!!
---------------------------------------
Dropping first TRs
++ 3dcalc: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: A cast of thousands
*+ WARNING:   If you are performing spatial transformations on an oblique dset,
  such as rest.nii.gz,
  or viewing/combining it with volumes of differing obliquity,
  you should consider running:
     3dWarp -deoblique
  on this and  other oblique datasets in the same session.
 See 3dWarp -help for details.
++ Oblique dataset:rest.nii.gz is 10.484265 degrees from plumb.
++ Output dataset ./rest_dr.nii.gz
Deobliquing 12
++ 3drefit: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: RW Cox
++ Processing AFNI dataset rest_dr.nii.gz
 + loading and re-writing entire dataset rest_dr.nii.gz
++ 3drefit processed 1 datasets
Reorienting 12
Motion correcting 12
++ 3dTstat: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: KR Hammett & RW Cox
++ Output dataset ./rest_ro_mean.nii.gz
++ 3dvolreg: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: RW Cox
++ Coarse del was 10, replaced with 3
++ Max displacement in automask = 0.25 (mm) at sub-brick 371
Skull stripping 12
++ 3dAutomask: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: Emperor Zhark
++ Loading dataset rest_mc.nii.gz
++ Forming automask
 + Fixed clip level = 279.744934
 + Used gradual clip level = 264.299774 .. 296.982086
 + Number voxels above clip level = 44079
 + Clustering voxels ...
 + Largest cluster has 43454 voxels
 + Clustering voxels ...
 + Largest cluster has 42786 voxels
 + Filled   461 voxels in small holes; now have 43247 voxels
 + Clustering voxels ...
 + Largest cluster has 43246 voxels
 + Clustering non-brain voxels ...
 + Clustering voxels ...
 + Largest cluster has 91922 voxels
 + Mask now has 43246 voxels
++ Dilating automask
 + Clustering voxels ...
 + Largest cluster has 85658 voxels
++ 49510 voxels in the mask [out of 135168: 36.63%]
++ first   9 x-planes are zero [from R]
++ last    8 x-planes are zero [from L]
++ first   0 y-planes are zero [from P]
++ last    4 y-planes are zero [from A]
++ first   0 z-planes are zero [from I]
++ last    0 z-planes are zero [from S]
++ Output dataset ./rest_mask.nii.gz
++ CPU time = 0.000000 sec
++ 3dcalc: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: A cast of thousands
++ Output dataset ./rest_ss.nii.gz
Getting example_func for registration for 12
++ 3dcalc: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: A cast of thousands
++ Output dataset ./example_func.nii.gz
Smoothing 12
Grand-mean scaling 12
Band-pass filtering 12
++ 3dFourier: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
Removing linear and quadratic trends for 12
++ 3dTstat: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: KR Hammett & RW Cox
++ Output dataset ./rest_filt_mean.nii.gz
++ 3dDetrend: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ 3dcalc: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: A cast of thousands
++ Output dataset ./rest_pp.nii.gz
Generating mask of preprocessed data for 12
------------------------------
!!!! RUNNING REGISTRATION !!!!
------------------------------
/bioimaging/andreash/topography/fcon_scripts/tissuepriors/3mm/
------------------------------
!!!! RUNNING SEGMENTATION !!!!
------------------------------
Segmenting brain for 12
Creating global mask
++ 3dcopy: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
Registering 12 csf to native (functional) space
Smoothing 12 csf
Registering 12 csf to standard space
Finding overlap between 12 csf and prior
Registering 12 csf back to native space
Threshold and binarize 12 csf probability map
Mask csf image by 12 functional
Registering 12 wm to native (functional) space
Smoothing 12 wm
Registering 12 wm to standard space
Finding overlap between 12 wm and prior
Registering 12 wm back to native space
Threshold and binarize 12 wm probability map
Mask wm image by 12 functional
--------------------------------------------
!!!! RUNNING NUISANCE SIGNAL REGRESSION !!!!
--------------------------------------------
Splitting up 12 motion parameters
Extracting global signal for 12
++ 3dmaskave: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
+++ 49510 voxels survive the mask
Extracting signal from csf for 12
++ 3dmaskave: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
+++ 319 voxels survive the mask
Extracting signal from white matter for 12
++ 3dmaskave: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
+++ 2113 voxels survive the mask
Modifying model file
Running feat model
Not enough data in /bioimaging/andreash/topography/12/func/nuisance/global.1D
++ 49510 voxels in mask
Running film to get residuals
Log directory is: /bioimaging/andreash/topography/12/func/nuisance/stats


An exception has been thrown
Unable to open /bioimaging/andreash/topography/12/func/nuisance/nuisance.mat
++ 3dTstat: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: KR Hammett & RW Cox
** FATAL ERROR: Can't open dataset /bioimaging/andreash/topography/12/func/nuisance/stats/res4d.nii.gz
** Program compile date = Jun 22 2011
++ 3dcalc: AFNI version=AFNI_2011_05_26_1457 (Jun 22 2011) [64-bit]
++ Authored by: A cast of thousands
** FATAL ERROR: can't open dataset /bioimaging/andreash/topography/12/func/nuisance/stats/res4d.nii.gz
** Program compile date = Jun 22 2011
** ERROR (nifti_image_read): failed to find header file for '/bioimaging/andreash/topography/12/func/rest_res'
** ERROR: nifti_image_open(/bioimaging/andreash/topography/12/func/rest_res): bad header info
ERROR: failed to open file /bioimaging/andreash/topography/12/func/rest_res
ERROR: Could not open image /bioimaging/andreash/topography/12/func/rest_res
Image Exception : #22 :: Failed to read volume /bioimaging/andreash/topography/12/func/rest_res
terminate called after throwing an instance of 'RBD_COMMON::BaseException'
/bioimaging/andreash/topography/fcon_scripts/5_nuisance.sh: line 95: 10368 Aborted                 flirt -ref ${reg_dir}/standard -in ${func_dir}/${rest}_res -out ${func_dir}/${rest}_res2standard -applyxfm -init ${reg_dir}/example_func2standard.mat -interp trilinear
+ + + + +
cat: rest: No such file or directory