Understanding brain structure and function entails the inclusion of anatomical and functional information in a common space, in order to study how these different informations relate to each other in a population of subjects. In this paper, we revisit the parcellation model and explicitly combine anatomical features, i.e. a segmentation of the cortex into gyri, with a functional information under the form of several cortical maps, which are used to further subdivide the gyri into functionally consistent regions. A probabilistic model is introduced, and the parcellation model is estimated using a Variational Bayes approach. The number of regions in the model is validated based on cross-validation. It is found that about 250 patches of cortex can be delineated both in the left and right hemisphere based on this procedure.