code: def scipy_gvoronoi(labeled):
      from scipy import ndimage
  
      L1, L2 = ndimage.distance_transform_edt(
          labeled == 0, return_distances=False, return_indices=True
      )
      return labeled[L1, L2]
  
imports: import numpy
  from mahotas.segmentation import gvoronoi
label_desc: Use image node_id 1 for calling a global function with example usage: scipy_gvoronoi(labeled) and returns: labeled