Thank you so much! I am using an in-house template which may be why it's in ORIG space. After converting my normalized functional to ORIG space, I am able to overlay my results on top of the templates (left image). However, if I use my template as an overlay on top of my normalized anatomical image, the template becomes undersampled (right image). I believe it's due to the template having a bigger matrix than my images, so the voxels become lost? Is there a way to alleviate this? I've tried 3dresample with various rmodes

Also, would it be possible to overlay both my functional and atlas over my anatomical scans?

Thank you for your help