Hi Clara,

So what you might want is to take the ROI coordinates from
+tlrc space, warp them to +orig space, and then drop a 3x3x3
box down at each coordinate (for a set number of ROIs). How
does that sound?

I have not done this, but I suspect it can be done as follows:

1. Given a list of coordinates in text file coords.tlrc.1D,
transform the coordinates to +orig space for this subject:

Vecwarp -apar ${subj}_anat+tlrc -backward \
-input coords.tlrc.1D -output coords.orig.1D

2. Add ROI values at each of those coordinates (assuming that
vals.1D has the values that match those coordinates):

1dcat coords.orig.1D vals.1D > $subj.coords.vals.1D

3. Turn this into a dataset matching the subject's stats dset:

3dUndump -master ${subj}_stats+orig -prefix $subj.one.vox \
-xyz $subj.coords.vals.1D

It is possible that you may want '-orient RAI' here...

4. Then use 3dcalc with differential sub-scripts to blow the
single voxels in 3x3x3 voxel blocks, or perhaps little
spheres.

Maybe this will do what you want.

- rick