python opencv 分水岭算法对图像进行分割

import numpy as np
import cv2
from matplotlib import pyplot as plt

img = cv2.imread('coin.jpg')
gray = cv2.cvtColor(img,cv2.COLOR_BGR2GRAY)
ret, thresh = cv2.threshold(gray,0,255,cv2.THRESH_BINARY_INV+cv2.THRESH_OTSU)

print(thresh.shape)
# cv2.imshow(thresh)
# cv2.waitKeyai(0)
# cv2.destroyAllWindows()
plt.imshow(thresh,'gray')

(312, 252)





<matplotlib.image.AxesImage at 0x2119d576438>

# noise removal
kernel = np.ones((3,3),np.uint8)
opening = cv2.morphologyEx(thresh,cv2.MORPH_OPEN,kernel, iterations = 2)

# sure background area
sure_bg = cv2.dilate(opening,kernel,iterations=3)

# Finding sure foreground area
dist_transform = cv2.distanceTransform(opening,cv2.DIST_L2,5)
ret, sure_fg = cv2.threshold(dist_transform,0.7*dist_transform.max(),255,0)
print(sure_fg.dtype)
# Finding unknown region
sure_fg = np.uint8(sure_fg)
print(sure_fg.dtype)
unknown = cv2.subtract(sure_bg,sure_fg)

plt.imshow(sure_fg,'gray')

float32
uint8





<matplotlib.image.AxesImage at 0x2119fc4e668>

# Marker labelling
ret, markers = cv2.connectedComponents(sure_fg)
print(markers.shape)
# Add one to all labels so that sure background is not 0, but 1
markers = markers+1

# Now, mark the region of unknown with zero
markers[unknown==255]=0
plt.imshow(sure_fg)

(312, 252)





<matplotlib.image.AxesImage at 0x2119fc95f98>

markers = cv2.watershed(img,markers)
img[markers == -1] = [255,0,0]

plt.imshow(img)

<matplotlib.image.AxesImage at 0x2119fd9b208>

plt.imshow(markers)

<matplotlib.image.AxesImage at 0x211a0dc4cf8>