Hi,
I want to define a 2D square region on a SUMA flatmap and obtain the indices / values of the corresponding AFNI voxels that map to that region.
A couple of years ago, I was using Surf2VolCoord to get this mapping between AFNI volume coordinates and flatmap surface coordinates in two steps:
Surf2VolCoord -i_fs /data/surf/$1/SUMA/lh.smoothwm.asc -grid_parent Processed/pb02.SUBJ.rall.scale_aligned_cropped_detrend+orig -sv "$1"_SurfVol_stripped_Alnd_Exp+orig > Processed/lhsurfNodeMapping.txt
Surf2VolCoord -i_fs /data/surf/$1/SUMA/lh.full.flat.patch.3d.asc -grid_parent Processed/pb02.SUBJ.rall.scale_aligned_cropped_detrend+orig > Processed/lhflatModel.txt
cut -d' ' -f5,6 Processed/lhflatModel.txt > Processed/lhflatVerts.txt
cut -d' ' -f15,16,17 Processed/lhsurfNodeMapping.txt > Processed/lhvoxelPositions.txt
paste Processed/lhvoxelPositions.txt Processed/lhflatVerts.txt > lhVoxel2Flat.txt
However, these commands written verbatim no longer work; when running either of the Surf2VolCoords commands I get the following error:
-- Error Surf2VolCoord_demo (SUMA_Surf2VolCoord_demo.c:447):
Have -1 dims for xyz, and 3 for surf.
Do you have any suggestions on what could have changed about Surf2VolCoords in the meantime or perhaps another way to get the mapping?
Thanks,
MCI