The larger picture is that I am interested in exploring how two patient groups differ in global-based connectivity (3dTCorrMap), and want to evaluate each voxel's connectivity with "all other voxels," ensuring that voxels identified belong to gray matter, independently based on Freesurfer gray matter maps.
One approach (the original post) I considered was using the gray matter map during preprocessing.