Determining cortical functional areas is an important goal for neurosciences and clinical neurosurgery. This paper presents a method for connectivity-based parcellation of the entire human cortical surface, exploiting the idea that each cortex region has a specific connection profile. The connectivity matrix of the cortex is computed using analytical Q-ball-based tractography. The parcellation is achieved independently for each subject and applied to the subset of the cortical surface endowed with enough connections to estimate safely a connectivity profile, namely the top of the cortical gyri. The key point of the method lies in a twofold reduction of the connectivity matrix dimension. First, parcellation amounts to iterating the clustering of Vorono?i patches of the cortical surface into parcels endowed with homogeneous profiles. The parcels without intersection with the patch boundaries are selected for the final parcellation. Before clustering a patch, the complete profiles are col...