Label image regions¶

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.patches as mpatches

from skimage import data
from skimage.filter import threshold_otsu
from skimage.segmentation import clear_border
from skimage.morphology import label, closing, square
from skimage.measure import regionprops
from skimage.color import label2rgb


image = data.coins()[50:-50, 50:-50]

# apply threshold
thresh = threshold_otsu(image)
bw = closing(image > thresh, square(3))

# remove artifacts connected to image border
cleared = bw.copy()
clear_border(cleared)

# label image regions
label_image = label(cleared)
borders = np.logical_xor(bw, cleared)
label_image[borders] = -1
image_label_overlay = label2rgb(label_image, image=image)

fig, ax = plt.subplots(ncols=1, nrows=1, figsize=(6, 6))
ax.imshow(image_label_overlay)

for region in regionprops(label_image):

    # skip small images
    if region.area < 100:
        continue

    # draw rectangle around segmented coins
    minr, minc, maxr, maxc = region.bbox
    rect = mpatches.Rectangle((minc, minr), maxc - minc, maxr - minr,
                              fill=False, edgecolor='red', linewidth=2)
    ax.add_patch(rect)

plt.show()

aW1wb3J0IG51bXB5IGFzIG5wCmltcG9ydCBtYXRwbG90bGliLnB5cGxvdCBhcyBwbHQKaW1wb3J0IG1hdHBsb3RsaWIucGF0Y2hlcyBhcyBtcGF0Y2hlcwoKZnJvbSBza2ltYWdlIGltcG9ydCBkYXRhCmZyb20gc2tpbWFnZS5maWx0ZXIgaW1wb3J0IHRocmVzaG9sZF9vdHN1CmZyb20gc2tpbWFnZS5zZWdtZW50YXRpb24gaW1wb3J0IGNsZWFyX2JvcmRlcgpmcm9tIHNraW1hZ2UubW9ycGhvbG9neSBpbXBvcnQgbGFiZWwsIGNsb3NpbmcsIHNxdWFyZQpmcm9tIHNraW1hZ2UubWVhc3VyZSBpbXBvcnQgcmVnaW9ucHJvcHMKZnJvbSBza2ltYWdlLmNvbG9yIGltcG9ydCBsYWJlbDJyZ2IKCgppbWFnZSA9IGRhdGEuY29pbnMoKVs1MDotNTAsIDUwOi01MF0KCiMgYXBwbHkgdGhyZXNob2xkCnRocmVzaCA9IHRocmVzaG9sZF9vdHN1KGltYWdlKQpidyA9IGNsb3NpbmcoaW1hZ2UgPiB0aHJlc2gsIHNxdWFyZSgzKSkKCiMgcmVtb3ZlIGFydGlmYWN0cyBjb25uZWN0ZWQgdG8gaW1hZ2UgYm9yZGVyCmNsZWFyZWQgPSBidy5jb3B5KCkKY2xlYXJfYm9yZGVyKGNsZWFyZWQpCgojIGxhYmVsIGltYWdlIHJlZ2lvbnMKbGFiZWxfaW1hZ2UgPSBsYWJlbChjbGVhcmVkKQpib3JkZXJzID0gbnAubG9naWNhbF94b3IoYncsIGNsZWFyZWQpCmxhYmVsX2ltYWdlW2JvcmRlcnNdID0gLTEKaW1hZ2VfbGFiZWxfb3ZlcmxheSA9IGxhYmVsMnJnYihsYWJlbF9pbWFnZSwgaW1hZ2U9aW1hZ2UpCgpmaWcsIGF4ID0gcGx0LnN1YnBsb3RzKG5jb2xzPTEsIG5yb3dzPTEsIGZpZ3NpemU9KDYsIDYpKQpheC5pbXNob3coaW1hZ2VfbGFiZWxfb3ZlcmxheSkKCmZvciByZWdpb24gaW4gcmVnaW9ucHJvcHMobGFiZWxfaW1hZ2UpOgoKICAgICMgc2tpcCBzbWFsbCBpbWFnZXMKICAgIGlmIHJlZ2lvbi5hcmVhIDwgMTAwOgogICAgICAgIGNvbnRpbnVlCgogICAgIyBkcmF3IHJlY3RhbmdsZSBhcm91bmQgc2VnbWVudGVkIGNvaW5zCiAgICBtaW5yLCBtaW5jLCBtYXhyLCBtYXhjID0gcmVnaW9uLmJib3gKICAgIHJlY3QgPSBtcGF0Y2hlcy5SZWN0YW5nbGUoKG1pbmMsIG1pbnIpLCBtYXhjIC0gbWluYywgbWF4ciAtIG1pbnIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGZpbGw9RmFsc2UsIGVkZ2Vjb2xvcj0ncmVkJywgbGluZXdpZHRoPTIpCiAgICBheC5hZGRfcGF0Y2gocmVjdCkKCnBsdC5zaG93KCk=