Thank you both Gang and Daniel, I think I'll stick with the surface-based approach.
One tiny question: Does it make sense to mapping the group analysis results that are totally from volume world to surface rendering, just for showing the cortex contrasts?
Like, select the statistical results from method 1 as the overlay, and open suma, have the spec and sv as the group average, does it make sense?