erosion=imerode(binary,strel('disk',integerValue));
D = bwdist(erosion);
D = -D;
D(~erosion) = -Inf;
L = watershed(D);

distanceD = -cv2.distanceTransform(1-not_im_eros, distanceType=cv2.DIST_L2, maskSize=0)
distanceDD = np.zeros_like(distanceD)
distanceDD[distanceD == 0] = 255
distanceD_positive = -distanceD
distanceD_uint8 = (distanceD_positive/(distanceD_positive.max()/255.0)).astype('uint8')
cv2.imshow('distanceD', 255-distanceD_uint8)
Ld = mp.watershed(distanceDD, mask = distanceD_positive)