def constrained_voronoi(self, seed):
"""Voronoi parcellation of the field starting from the input seed
Parameters
----------
seed: int array of shape(p), the input seeds
Returns
-------
label: The resulting labelling of the data
Fixme
-----
deal with graphs with several ccs
"""
if np.size(self.field) == 0:
raise ValueError('No field has been defined so far')
seed = seed.astype(np.int)
weights = np.sqrt(np.sum((self.field[self.edges.T[0]] -
self.field[self.edges.T[1]]) ** 2, 1))
g = WeightedGraph(self.V, self.edges, weights)
label = g.voronoi_labelling(seed)
return label
