Measure region properties¶

import math
import matplotlib.pyplot as plt
import numpy as np

from skimage.draw import ellipse
from skimage.morphology import label
from skimage.measure import regionprops
from skimage.transform import rotate


image = np.zeros((600, 600))

rr, cc = ellipse(300, 350, 100, 220)
image[rr,cc] = 1

image = rotate(image, angle=15, order=0)

label_img = label(image)
regions = regionprops(label_img)

fig, ax = plt.subplots()
ax.imshow(image, cmap=plt.cm.gray)

for props in regions:
    y0, x0 = props.centroid
    orientation = props.orientation
    x1 = x0 + math.cos(orientation) * 0.5 * props.major_axis_length
    y1 = y0 - math.sin(orientation) * 0.5 * props.major_axis_length
    x2 = x0 - math.sin(orientation) * 0.5 * props.minor_axis_length
    y2 = y0 - math.cos(orientation) * 0.5 * props.minor_axis_length

    ax.plot((x0, x1), (y0, y1), '-r', linewidth=2.5)
    ax.plot((x0, x2), (y0, y2), '-r', linewidth=2.5)
    ax.plot(x0, y0, '.g', markersize=15)

    minr, minc, maxr, maxc = props.bbox
    bx = (minc, maxc, maxc, minc, minc)
    by = (minr, minr, maxr, maxr, minr)
    ax.plot(bx, by, '-b', linewidth=2.5)

ax.axis((0, 600, 600, 0))
plt.show()

aW1wb3J0IG1hdGgKaW1wb3J0IG1hdHBsb3RsaWIucHlwbG90IGFzIHBsdAppbXBvcnQgbnVtcHkgYXMgbnAKCmZyb20gc2tpbWFnZS5kcmF3IGltcG9ydCBlbGxpcHNlCmZyb20gc2tpbWFnZS5tb3JwaG9sb2d5IGltcG9ydCBsYWJlbApmcm9tIHNraW1hZ2UubWVhc3VyZSBpbXBvcnQgcmVnaW9ucHJvcHMKZnJvbSBza2ltYWdlLnRyYW5zZm9ybSBpbXBvcnQgcm90YXRlCgoKaW1hZ2UgPSBucC56ZXJvcygoNjAwLCA2MDApKQoKcnIsIGNjID0gZWxsaXBzZSgzMDAsIDM1MCwgMTAwLCAyMjApCmltYWdlW3JyLGNjXSA9IDEKCmltYWdlID0gcm90YXRlKGltYWdlLCBhbmdsZT0xNSwgb3JkZXI9MCkKCmxhYmVsX2ltZyA9IGxhYmVsKGltYWdlKQpyZWdpb25zID0gcmVnaW9ucHJvcHMobGFiZWxfaW1nKQoKZmlnLCBheCA9IHBsdC5zdWJwbG90cygpCmF4Lmltc2hvdyhpbWFnZSwgY21hcD1wbHQuY20uZ3JheSkKCmZvciBwcm9wcyBpbiByZWdpb25zOgogICAgeTAsIHgwID0gcHJvcHMuY2VudHJvaWQKICAgIG9yaWVudGF0aW9uID0gcHJvcHMub3JpZW50YXRpb24KICAgIHgxID0geDAgKyBtYXRoLmNvcyhvcmllbnRhdGlvbikgKiAwLjUgKiBwcm9wcy5tYWpvcl9heGlzX2xlbmd0aAogICAgeTEgPSB5MCAtIG1hdGguc2luKG9yaWVudGF0aW9uKSAqIDAuNSAqIHByb3BzLm1ham9yX2F4aXNfbGVuZ3RoCiAgICB4MiA9IHgwIC0gbWF0aC5zaW4ob3JpZW50YXRpb24pICogMC41ICogcHJvcHMubWlub3JfYXhpc19sZW5ndGgKICAgIHkyID0geTAgLSBtYXRoLmNvcyhvcmllbnRhdGlvbikgKiAwLjUgKiBwcm9wcy5taW5vcl9heGlzX2xlbmd0aAoKICAgIGF4LnBsb3QoKHgwLCB4MSksICh5MCwgeTEpLCAnLXInLCBsaW5ld2lkdGg9Mi41KQogICAgYXgucGxvdCgoeDAsIHgyKSwgKHkwLCB5MiksICctcicsIGxpbmV3aWR0aD0yLjUpCiAgICBheC5wbG90KHgwLCB5MCwgJy5nJywgbWFya2Vyc2l6ZT0xNSkKCiAgICBtaW5yLCBtaW5jLCBtYXhyLCBtYXhjID0gcHJvcHMuYmJveAogICAgYnggPSAobWluYywgbWF4YywgbWF4YywgbWluYywgbWluYykKICAgIGJ5ID0gKG1pbnIsIG1pbnIsIG1heHIsIG1heHIsIG1pbnIpCiAgICBheC5wbG90KGJ4LCBieSwgJy1iJywgbGluZXdpZHRoPTIuNSkKCmF4LmF4aXMoKDAsIDYwMCwgNjAwLCAwKSkKcGx0LnNob3coKQ==