One option is to combine the spheres with white matter with 3dcalc. You can generate white matter masks with FreeSurfer or even 3dSeg. Mask out the data with something like this:
3dcalc -a masks+orig. -b wm_mask+orig -expr 'a*not(b)' -prefix masks_no_wm