Shapes¶

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

from skimage.draw import (line, polygon, circle,
                          circle_perimeter,
                          ellipse, ellipse_perimeter,
                          bezier_curve)


fig, (ax1, ax2) = plt.subplots(ncols=2, nrows=1, figsize=(10, 6))


img = np.zeros((500, 500, 3), dtype=np.double)

# draw line
rr, cc = line(120, 123, 20, 400)
img[rr, cc, 0] = 255

# fill polygon
poly = np.array((
    (300, 300),
    (480, 320),
    (380, 430),
    (220, 590),
    (300, 300),
))
rr, cc = polygon(poly[:, 0], poly[:, 1], img.shape)
img[rr, cc, 1] = 1

# fill circle
rr, cc = circle(200, 200, 100, img.shape)
img[rr, cc, :] = (1, 1, 0)

# fill ellipse
rr, cc = ellipse(300, 300, 100, 200, img.shape)
img[rr, cc, 2] = 1

# circle
rr, cc = circle_perimeter(120, 400, 15)
img[rr, cc, :] = (1, 0, 0)

# Bezier curve
rr, cc = bezier_curve(70, 100, 10, 10, 150, 100, 1)
img[rr, cc, :] = (1, 0, 0)

# ellipses
rr, cc = ellipse_perimeter(120, 400, 60, 20, orientation=math.pi / 4.)
img[rr, cc, :] = (1, 0, 1)
rr, cc = ellipse_perimeter(120, 400, 60, 20, orientation=-math.pi / 4.)
img[rr, cc, :] = (0, 0, 1)
rr, cc = ellipse_perimeter(120, 400, 60, 20, orientation=math.pi / 2.)
img[rr, cc, :] = (1, 1, 1)

ax1.imshow(img)
ax1.set_title('No anti-aliasing')
ax1.axis('off')


from skimage.draw import line_aa, circle_perimeter_aa


img = np.zeros((100, 100), dtype=np.double)

# anti-aliased line
rr, cc, val = line_aa(12, 12, 20, 50)
img[rr, cc] = val

# anti-aliased circle
rr, cc, val = circle_perimeter_aa(60, 40, 30)
img[rr, cc] = val


ax2.imshow(img, cmap=plt.cm.gray, interpolation='nearest')
ax2.set_title('Anti-aliasing')
ax2.axis('off')

plt.show()

aW1wb3J0IG1hdGgKaW1wb3J0IG51bXB5IGFzIG5wCmltcG9ydCBtYXRwbG90bGliLnB5cGxvdCBhcyBwbHQKCmZyb20gc2tpbWFnZS5kcmF3IGltcG9ydCAobGluZSwgcG9seWdvbiwgY2lyY2xlLAogICAgICAgICAgICAgICAgICAgICAgICAgIGNpcmNsZV9wZXJpbWV0ZXIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgZWxsaXBzZSwgZWxsaXBzZV9wZXJpbWV0ZXIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgYmV6aWVyX2N1cnZlKQoKCmZpZywgKGF4MSwgYXgyKSA9IHBsdC5zdWJwbG90cyhuY29scz0yLCBucm93cz0xLCBmaWdzaXplPSgxMCwgNikpCgoKaW1nID0gbnAuemVyb3MoKDUwMCwgNTAwLCAzKSwgZHR5cGU9bnAuZG91YmxlKQoKIyBkcmF3IGxpbmUKcnIsIGNjID0gbGluZSgxMjAsIDEyMywgMjAsIDQwMCkKaW1nW3JyLCBjYywgMF0gPSAyNTUKCiMgZmlsbCBwb2x5Z29uCnBvbHkgPSBucC5hcnJheSgoCiAgICAoMzAwLCAzMDApLAogICAgKDQ4MCwgMzIwKSwKICAgICgzODAsIDQzMCksCiAgICAoMjIwLCA1OTApLAogICAgKDMwMCwgMzAwKSwKKSkKcnIsIGNjID0gcG9seWdvbihwb2x5WzosIDBdLCBwb2x5WzosIDFdLCBpbWcuc2hhcGUpCmltZ1tyciwgY2MsIDFdID0gMQoKIyBmaWxsIGNpcmNsZQpyciwgY2MgPSBjaXJjbGUoMjAwLCAyMDAsIDEwMCwgaW1nLnNoYXBlKQppbWdbcnIsIGNjLCA6XSA9ICgxLCAxLCAwKQoKIyBmaWxsIGVsbGlwc2UKcnIsIGNjID0gZWxsaXBzZSgzMDAsIDMwMCwgMTAwLCAyMDAsIGltZy5zaGFwZSkKaW1nW3JyLCBjYywgMl0gPSAxCgojIGNpcmNsZQpyciwgY2MgPSBjaXJjbGVfcGVyaW1ldGVyKDEyMCwgNDAwLCAxNSkKaW1nW3JyLCBjYywgOl0gPSAoMSwgMCwgMCkKCiMgQmV6aWVyIGN1cnZlCnJyLCBjYyA9IGJlemllcl9jdXJ2ZSg3MCwgMTAwLCAxMCwgMTAsIDE1MCwgMTAwLCAxKQppbWdbcnIsIGNjLCA6XSA9ICgxLCAwLCAwKQoKIyBlbGxpcHNlcwpyciwgY2MgPSBlbGxpcHNlX3BlcmltZXRlcigxMjAsIDQwMCwgNjAsIDIwLCBvcmllbnRhdGlvbj1tYXRoLnBpIC8gNC4pCmltZ1tyciwgY2MsIDpdID0gKDEsIDAsIDEpCnJyLCBjYyA9IGVsbGlwc2VfcGVyaW1ldGVyKDEyMCwgNDAwLCA2MCwgMjAsIG9yaWVudGF0aW9uPS1tYXRoLnBpIC8gNC4pCmltZ1tyciwgY2MsIDpdID0gKDAsIDAsIDEpCnJyLCBjYyA9IGVsbGlwc2VfcGVyaW1ldGVyKDEyMCwgNDAwLCA2MCwgMjAsIG9yaWVudGF0aW9uPW1hdGgucGkgLyAyLikKaW1nW3JyLCBjYywgOl0gPSAoMSwgMSwgMSkKCmF4MS5pbXNob3coaW1nKQpheDEuc2V0X3RpdGxlKCdObyBhbnRpLWFsaWFzaW5nJykKYXgxLmF4aXMoJ29mZicpCgoKZnJvbSBza2ltYWdlLmRyYXcgaW1wb3J0IGxpbmVfYWEsIGNpcmNsZV9wZXJpbWV0ZXJfYWEKCgppbWcgPSBucC56ZXJvcygoMTAwLCAxMDApLCBkdHlwZT1ucC5kb3VibGUpCgojIGFudGktYWxpYXNlZCBsaW5lCnJyLCBjYywgdmFsID0gbGluZV9hYSgxMiwgMTIsIDIwLCA1MCkKaW1nW3JyLCBjY10gPSB2YWwKCiMgYW50aS1hbGlhc2VkIGNpcmNsZQpyciwgY2MsIHZhbCA9IGNpcmNsZV9wZXJpbWV0ZXJfYWEoNjAsIDQwLCAzMCkKaW1nW3JyLCBjY10gPSB2YWwKCgpheDIuaW1zaG93KGltZywgY21hcD1wbHQuY20uZ3JheSwgaW50ZXJwb2xhdGlvbj0nbmVhcmVzdCcpCmF4Mi5zZXRfdGl0bGUoJ0FudGktYWxpYXNpbmcnKQpheDIuYXhpcygnb2ZmJykKCnBsdC5zaG93KCk=