# Copy+pasteable commands from the AFNI Bootcamp's "afni11_roi.pdf"
# presentation


cd AFNI_data6/roi_demo

# Change grid/resolution of data set from (high res) anat to 
# (low res) EPI
3dresample                   \
    -master rall_vr+orig     \
    -prefix anat_roi_resam   \
    -inset anat_roi+orig     \
    -rmode NN


# Calculate average value per time point, dump to text file...
3dmaskave                      \
    -mask anat_roi_resam+orig  \
    -quiet                     \
    rall_vr+orig > epi_avg.1D

# ... and view dumped text file results, using either of these
more epi_avg.1D
1dplot -yaxis 1000:1200:2:1 epi_avg.1D


# Dump *values* of a dset that are within a mask region into a 
# text file; no coordinate locations saved in this case, but could 
# save the coordinates of each voxel, too.
3dmaskdump                                \
    -noijk                                \
    -mask anat_roi_resam+orig             \
    func_slim+orig'[2]' > Vrel-tstats.txt

# ... and view dumped text file results
more Vrel-tstats.txt


# Compute separate statistics for each ROI in a volume
3dROIstats                  \
    -mask 3rois+orig        \
    func_slim+orig'[0]'


# *** 3dclust example here, but use updated notation!! ***


# Find suprathreshold voxels, and make into clusters (kinda like GUI's
# Clusterize functionality and/or 3dclust)
3dmerge                       \
    -prefix func_roi          \
    -dxyz=1                   \
    -1clip 99.0               \
    -1clust_order 1 200       \
    func_slim+orig.'[0]'


# ============ "Getting around with spheres" example ==========
# Try this:

adwarp                       \
    -apar anat+tlrc          \
    -dpar func_slim+orig     \
    –dxyz 3

# 1) In afni GUI, select:
#         Underlay -> anat, 
#         Overlay  -> func_slim
# 2) Switch 'view' to Talairach
# 3) Overlay index   -> vrel_coef,
#    Threshold index -> vrel_tstat
# 4) Clusterize (default opts OK), and Save Tabl
1dcat Clust_table.1D'[4..6]' > Clust_PeakXYZ.1D
3dUndump -srad 7.5 -master func_slim+tlrc -orient RAI   \
      -prefix clust_spheres -xyz Clust_PeakXYZ.1D