That worked! Thank you so much!! I assume that since I have the old @SSwarper, I am utilizing the MNI152_2009_template.nii.gz which has 1 volume instead of the MNI152_2009_template_SSW.nii.gz which has several?
Interestingly, I did try using the "old" syntax and @SSwarper wouldn't run. It seemed to be taking an awfully long time and produced several "junk" files. I have pasted the code I used and the terminal output below. If it is not an issue to use the "new" code, changing the name of my output files and then running afni_proc.py, I will continue doing this.
Thank you again for your assistance!
SSwarper code:
@SSwarper \
11987.anat+orig \
11987
SSwarper output (not working)
-------- freesurfer-Linux-centos6_x86_64-stable-pub-v6.0.0-2beb96c --------
Setting up environment for FreeSurfer/FS-FAST (and FSL)
FREESURFER_HOME /usr/local/freesurfer
FSFAST_HOME /usr/local/freesurfer/fsfast
FSF_OUTPUT_FORMAT nii.gz
SUBJECTS_DIR /home/bmiadmin/subjects
MNI_DIR /usr/local/freesurfer/mni
SUBJECTS_DIR= /home/bmiadmin/subjects
mitchelllab-linux:~/Tamara/PRR_Analysis/11987> @SSwarper \
? 11987.anat+orig \
? 11987
SUBJECTS_DIR= /home/bmiadmin/subjects
++ 3dUnifize: AFNI version=AFNI_18.1.09 (Apr 30 2018) [64-bit]
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as ./11987.anat+orig.BRIK,
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:./11987.anat+orig.BRIK is 0.142653 degrees from plumb.
+ Pre-processing: ADV................UWG
++ Output dataset ./anatU.11987.nii
++ ===== Elapsed = 23.3 sec
++ Notice SUMA_BrainWrap_ParseInput (SUMA_3dSkullStrip.c:1219 @14:36:45):
Setting fill_hole to 10
++ Notice 3dSkullStrip (SUMA_3dSkullStrip.c:1529 @14:36:45):
Loading dset, performing Spatial Normalization
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as anatU.11987.nii,
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:anatU.11987.nii is 0.142653 degrees from plumb.
3dSkullStrip:
Original resolution 1.000002, -1.000002, 1.000003
SpatNorm resolution 1.000002, 1.000002, 1.000003
++mri_brainormalize: normalizing
++mri_brainormalize: top clip above slice 254
++mri_brainormalize: bot clip below slice 84
++mri_brainormalize: Rescaling.
++mri_brainormalize: Eroding, clustering ...
+ Eroded 362572 voxels
+ Restored 300359 eroded voxels
+ Clustering voxels ...
+ Largest cluster has 3188439 voxels
+ Clustering voxels ...
+ Largest cluster has 2979588 voxels
++mri_brainormalize: Clipping extreme values.
+ scaling by fac=0.0524691
3dSkullStrip: -spatnorm: Expecting 6195700 voxels in in_vol dset (167 212 175)
3dSkullStrip: Edge threshold set to 47.128998. (minimum acceptable was 22.842131)
(92.000000 percentile =47.128998, 99.999001 percentile = 228.421310)
3dSkullStrip: Prepping volume...
SUMA_LoadPrepInVol: Volume 3213642.000000, radius 91.545380
3dSkullStrip: Beginning brain extraction...
3dSkullStrip: Have surface, OK for 1st entry.
3dSkullStrip: Stretching to avoid ventricles.
SO->Center is [1.055698 16.702696 -12.525011]
SUMA_StretchToFitLeCerveau:
In stage 1
SUMA_StretchToFitLeCerveau:
Stage1: satiated, low MaxExp
SUMA_StretchToFitLeCerveau:
reducing tightness, applying touchup with Stage2Type = 1
SUMA_StretchToFitLeCerveau:
Stage 2, type 1: 3650 troubled nodes, going for more...
SUMA_StretchToFitLeCerveau:
In stage 1
SUMA_StretchToFitLeCerveau:
Stage1: satiated, low MaxExp
SUMA_StretchToFitLeCerveau:
reducing tightness, applying touchup with Stage2Type = 1
SUMA_StretchToFitLeCerveau:
Stage 2, type 1: 0.100822 change in troubled nodes.
SUMA_StretchToFitLeCerveau:
Continuing with Stage.
SUMA_StretchToFitLeCerveau:
In stage 1
SUMA_StretchToFitLeCerveau:
Stage1: pastarea = 79507.312004
SUMA_StretchToFitLeCerveau:
Stage1: MaxExp = 0.683446, darea = 0.000000, going for more...
SUMA_StretchToFitLeCerveau:
In stage 1
SUMA_StretchToFitLeCerveau:
Stage1: satiated, small area differential.
SUMA_StretchToFitLeCerveau:
reducing tightness, applying touchup with Stage2Type = 1
SUMA_StretchToFitLeCerveau:
Stage 2, type 1: 0.086837 change in troubled nodes.
SUMA_StretchToFitLeCerveau:
Continuing with Stage.
SUMA_StretchToFitLeCerveau:
In stage 1
SUMA_StretchToFitLeCerveau:
Stage1: satiated, small area differential.
SUMA_StretchToFitLeCerveau:
reducing tightness, applying touchup with Stage2Type = 1
SUMA_StretchToFitLeCerveau:
Stage 2, type 1: 0.025359 change in troubled nodes.
SUMA_StretchToFitLeCerveau:
Continuing with Stage.
SUMA_StretchToFitLeCerveau:
In stage 1
SUMA_StretchToFitLeCerveau:
Stage1: satiated, small area differential.
SUMA_StretchToFitLeCerveau:
reducing tightness, applying touchup with Stage2Type = 1
SUMA_StretchToFitLeCerveau:
Stage 2, type 1: 0.039028 change in troubled nodes.
SUMA_StretchToFitLeCerveau:
Continuing with Stage.
SUMA_StretchToFitLeCerveau:
In stage 1
SUMA_StretchToFitLeCerveau:
Stage1: satiated, small area differential.
SUMA_StretchToFitLeCerveau:
reducing tightness, applying touchup with Stage2Type = 1
SUMA_StretchToFitLeCerveau:
Stage 2, type 1: 0.014250 change in troubled nodes.
SUMA_StretchToFitLeCerveau:
Continuing with Stage.
SUMA_StretchToFitLeCerveau:
In stage 1
SUMA_StretchToFitLeCerveau:
Stage1: satiated, small area differential.
SUMA_StretchToFitLeCerveau:
reducing tightness, applying touchup with Stage2Type = 1
SUMA_StretchToFitLeCerveau:
Stage 2, type 1: 0.000723 change in troubled nodes.
SUMA_StretchToFitLeCerveau:
Stage converged. Moving to new Stage or Type.
SUMA_StretchToFitLeCerveau:
In stage 1
SUMA_StretchToFitLeCerveau:
Stage1: satiated, small area differential.
SUMA_StretchToFitLeCerveau:
Stage 3, type 2: 1.000000 change in troubled nodes.
SUMA_StretchToFitLeCerveau:
Continuing with Stage.
3dSkullStrip: Checking for self intersection...
3dSkullStrip: No intersections found.
3dSkullStrip: Final smoothing of 72
3dSkullStrip: Touchup correction, pass 1 ...
3dSkullStrip: The beauty treatment smoothing.
3dSkullStrip: Locating voxels inside surface ...
3dSkullStrip: Setting output to orig_vol (iset = 0x2baece8) (osnd = 0x3d84c98)...
3dSkullStrip: Coercing...
3dSkullStrip: Writing masked volume ...
3dSkullStrip: Applying a bit of erosion and dilatation
++ OpenMP thread count = 8
++ 3dQwarp: AFNI version=AFNI_18.1.09 (Apr 30 2018) [64-bit]
++ Authored by: Zhark the (Hermite) Cubically Warped
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as anatU.11987.nii,
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:anatU.11987.nii is 0.142653 degrees from plumb.
++ Starting 3dAllineate (affine register) command:
3dAllineate -base /home/bmiadmin/abin/MNI152_2009_template.nii.gz[1] -source anatU.11987.nii -prefix AFN_vA68B827xfMNFbMsh9wpPA.nii -1Dmatrix_save AFN_vA68B827xfMNFbMsh9wpPA -cmass -final wsinc5 -float -master BASE -twobest 7 -zclip -fineblur 4.44
++ ###########################################################
++ 3dAllineate: AFNI version=AFNI_18.1.09 (Apr 30 2018) [64-bit]
++ Authored by: Zhark the Registrator
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as anatU.11987.nii,
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:anatU.11987.nii is 0.142653 degrees from plumb.
++ Source dataset: ./anatU.11987.nii
++ Base dataset: /home/bmiadmin/abin/MNI152_2009_template.nii.gz
++ Loading datasets
++ 6527249 voxels in source_automask+0
++ Zero-pad: ybot=0 ytop=1
++ Computing -autobox
++ 8007930 voxels [93.5%] in weight mask
++ Number of points for matching = 3763727
++ NOTE: base and source coordinate systems have different handedness
+ Orientations: base=Right handed (LPI); source=Left handed (ASR)
++ OpenMP thread count = 8
++ ======= Allineation of 1 sub-bricks using Hellinger metric =======
++ *** Coarse pass begins ***
++ *** Fine pass begins ***
+ Final fine fit Parameters:
x-shift=-1.8165 y-shift=-20.5583 z-shift= 5.4606
z-angle= 1.1994 x-angle= 9.2621 y-angle= 2.9879
x-scale= 0.9036 y-scale= 0.9894 z-scale= 0.8789
y/x-shear= 0.0144 z/x-shear=-0.0205 z/y-shear=-0.0852
++ Output dataset ./AFN_vA68B827xfMNFbMsh9wpPA.nii
++ Wrote -1Dmatrix_save AFN_vA68B827xfMNFbMsh9wpPA.aff12.1D
++ 3dAllineate: total CPU time = 0.0 sec Elapsed = 379.2
++ ###########################################################
++ # Please check results visually for alignment quality #
++ ###########################################################
++ 3dQwarp: replacing source dataset with 3dAllineate result junk.SSwarper.HCjcBWREg5A_TAL5_Allin.nii
+ 3dAllineate output files have been renamed
# mat44 3dAllineate matrix:
0.901998 0.026467 -0.046485 -1.816440
-0.007442 0.976627 0.158569 -20.558300
0.031974 -0.224018 0.853645 5.460660
++ Dataset zero-pad: xbot=24 xtop=23 ybot=18 ytop=27 zbot=24 ztop=18 voxels
++ +++++++++++ Begin warp optimization: base=/home/bmiadmin/abin/MNI152_2009_template.nii.gz source=./junk.SSwarper.HCjcBWREg5A_TAL5_Allin.nii
++ AFNI warpomatic: 240 x 274 x 235 volume ; autobbox = 30..209 31..243 30..204 [clock= 6m 21s 916ms]
lev=0 8..231 5..269 9..225: [first cost=-0.44210] .. done [cost:-0.44210==>-0.55226]
lev=1 patch=181x207x177 [clock= 8m 24s 617ms]........ done [cost:-0.55226==>-0.58763 ; 8 patches optimized, 0 skipped]
lev=2 patch=135x155x133 [clock= 13m 59s 902ms]..................... done [cost:-0.58763==>-0.65361 ; 21 patches optimized, 6 skipped]
lev=3 patch=101x117x99 [clock= 20m 15s 777ms]............... done [cost:-0.65361==>-0.71383 ; 39 patches optimized, 25 skipped]
lev=4 patch=77x87x75 [clock= 24m 58s 612ms]......... done [cost:-0.71383==>-0.77934 ; 80 patches optimized, 45 skipped]
lev=5 patch=57x65x57 [clock= 29m 25s 560ms]................. done [cost:-0.77934==>-0.83808 ; 160 patches optimized, 92 skipped]
++ ========== total number of parameters 'optimized' = 7497
++ Output dataset ./junk.SSwarper.HCjcBWREg5A_TAL5_AWARP.nii
++ Output dataset ./junk.SSwarper.HCjcBWREg5A_TAL5.nii
++ Output dataset ./junk.SSwarper.HCjcBWREg5A_TAL5_WARP.nii
++ Inverting warp *******............:..
++ Output dataset ./junk.SSwarper.HCjcBWREg5A_TAL5_WARPINV.nii
++ ===== clock time = 33m 44s 865ms
++ no -frac option: defaulting to -union
++ processing 1 input datasets...
++ padding all datasets by 1 (for dilations)
++ frac 0 over 1 volumes gives min count 0
++ voxel limits: 0 clipped, 2000189 survived, 6529832 were zero
++ writing result junk.SSwarper.HCjcBWREg5A_MASK.nii...
++ Output dataset ./junk.SSwarper.HCjcBWREg5A_MASK.nii
++ 3dcalc: AFNI version=AFNI_18.1.09 (Apr 30 2018) [64-bit]
++ Authored by: A cast of thousands
++ Output dataset ./junk.SSwarper.HCjcBWREg5A_TAL5mm.nii
++ 3dNwarpApply: AFNI version=AFNI_18.1.09 (Apr 30 2018) [64-bit]
++ Authored by: Zhark the Warped
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as anatS.11987.nii,
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:anatS.11987.nii is 0.142653 degrees from plumb.
++ -master dataset is 'anatS.11987.nii'
** AFNI converts NIFTI_datatype=2 (UINT8) in file junk.SSwarper.HCjcBWREg5A_TAL5mm.nii to FLOAT32
Warnings of this type will be muted for this session.
Set AFNI_NIFTI_TYPE_WARN to YES to see them all, NO to see none.
++ opened source dataset 'junk.SSwarper.HCjcBWREg5A_TAL5mm.nii'
++ Processing -nwarp
++ Warping:.Z
++ Output dataset ./junk.SSwarper.HCjcBWREg5A_TAL5ww.nii
++ total Elapsed time = 13.5 sec
++ no -frac option: defaulting to -union
++ processing 1 input datasets...
++ padding all datasets by 2 (for dilations)
++ frac 0 over 1 volumes gives min count 0
++ voxel limits: 0 clipped, 2117569 survived, 6412452 were zero
++ writing result junk.SSwarper.HCjcBWREg5A_MASK.nii...
++ Output dataset ./junk.SSwarper.HCjcBWREg5A_MASK.nii
++ 3dNwarpApply: AFNI version=AFNI_18.1.09 (Apr 30 2018) [64-bit]
++ Authored by: Zhark the Warped
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as anatS.11987.nii,
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:anatS.11987.nii is 0.142653 degrees from plumb.
++ -master dataset is 'anatS.11987.nii'
++ opened source dataset 'junk.SSwarper.HCjcBWREg5A_MASK.nii'
++ Processing -nwarp
++ Warping:.Z
++ Output dataset ./junk.SSwarper.HCjcBWREg5A_MASKO.nii
++ total Elapsed time = 11.0 sec
++ 3dcalc: AFNI version=AFNI_18.1.09 (Apr 30 2018) [64-bit]
++ Authored by: A cast of thousands
*+ WARNING: dataset 'b'=junk.SSwarper.HCjcBWREg5A_TAL5ww.nii has an obliquity difference of 0.142653 degrees with anatS.11987.nii
*+ WARNING: dataset 'c'=junk.SSwarper.HCjcBWREg5A_MASKO.nii has an obliquity difference of 0.142653 degrees with anatS.11987.nii
++ Output dataset ./anatSS.11987.nii
++ 3dAllineate: AFNI version=AFNI_18.1.09 (Apr 30 2018) [64-bit]
++ Authored by: Zhark the Registrator
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as anatSS.11987.nii,
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:anatSS.11987.nii is 0.142653 degrees from plumb.
** AFNI converts NIFTI_datatype=2 (UINT8) in file junk.SSwarper.HCjcBWREg5A_TAL5mm.nii to FLOAT32
Warnings of this type will be muted for this session.
Set AFNI_NIFTI_TYPE_WARN to YES to see them all, NO to see none.
++ Source dataset: ./anatSS.11987.nii
++ Base dataset: (not given)
++ Loading datasets
++ OpenMP thread count = 8
++ ========== Applying transformation to 1 sub-bricks ==========
++ wsinc5 interpolation setup:
+ taper function = Min sidelobe 3 term
+ taper cut point = 0.000
+ window radius = 5 voxels
+ window shape = Cubical
+ The above can be altered via the AFNI_WSINC5_* environment variables.
+ (To avoid this message, 'setenv AFNI_WSINC5_SILENT YES'.)
+ wsinc5 CUBE(5) mask has 1000 points
++ Output dataset ./junk.SSwarper.HCjcBWREg5A_AffSS.nii
++ 3dAllineate: total CPU time = 0.0 sec Elapsed = 3.5
++ ###########################################################
++ OpenMP thread count = 8
++ 3dQwarp: AFNI version=AFNI_18.1.09 (Apr 30 2018) [64-bit]
++ Authored by: Zhark the (Hermite) Cubically Warped
++ Dataset zero-pad: xbot=24 xtop=23 ybot=18 ytop=27 zbot=24 ztop=18 voxels
++ Weightizing the base image
++ +++++++++++ Begin warp optimization: base=/home/bmiadmin/abin/MNI152_2009_template.nii.gz source=./junk.SSwarper.HCjcBWREg5A_AffSS.nii
++ AFNI warpomatic: 240 x 274 x 235 volume ; autobbox = 43..197 38..228 25..188 [clock= 4s 379ms]
lev=3 patch=101x117x99 [clock= 5s 628ms][first cost=-1.60012] ..... done [cost:-1.60012==>-1.70164 ; 16 patches optimized, 11 skipped]
lev=4 patch=77x87x75 [clock= 2m 20s 979ms][first cost=-1.49978] ... done [cost:-1.70164==>-1.54957 ; 37 patches optimized, 27 skipped]
lev=5 patch=57x65x57 [clock= 4m 44s 493ms][first cost=-1.36680] .........:[cost=-1.42939]:....... done [cost:-1.54957==>-1.43467 ; 241 patches optimized, 191 skipped]
lev=6 patch=43x49x43 [clock= 10m 51s 733ms][first cost=-1.28611] .........:[cost=-1.34538]:............ done [cost:-1.43467==>-1.36036 ; 502 patches optimized, 394 skipped]
lev=7 patch=33x37x31 [clock= 16m 4s 742ms][first cost=-1.25315] .................:[cost=-1.30743]:................... done [cost:-1.36036==>-1.29359 ; 1261 patches optimized, 1159 skipped]
lev=8 patch=25x27x23 [clock= 21m 45s 351ms][first cost=-1.23734] .........................:[cost=-1.29184]:............................ done [cost:-1.29359==>-1.23711 ; 3047 patches optimized, 2803 skipped]
lev=9 patch=19x21x17 [clock= 27m 42s 299ms][first cost=-1.24843] .............................................. done [cost:-1.23711==>-1.30319 ; 3425 patches optimized, 3358 skipped]
lev=10 patch=13x15x13 [clock= 31m 28s 631ms][first cost=-1.26157] ........................................................
