Hi Bin,
If you're lucky enough to have multiple anatomical scans from a single subject (from a single session), you can run 3dUniformize to do an intensity normalization on the images (makes for nice gray-white matter contrast), register the images with 3dvolreg, and then average them with 3dMean. Do a Talairach transformation and the end result is a beautiful anatomical underlay for displaying your group's fMRI data. Here's an example:
[
afni.nimh.nih.gov]
(see Slide #5)
Good Luck,
Peggy