Thank you Daniel. I'm familiar with the -omask option. However, is there an existing function to use results of the "omask" option of whereami in breaking down 3dclust output (peak voxel, center of mass coordinates, etc) for each cluster? I suppose this might be overkill in terms of reporting results since mapping between significant clusters and ROIs on any atlas are ultimately approximations to begin with.