mirror of
https://github.com/wassname/scikit-image.git
synced 2026-06-28 23:42:02 +08:00
Johannes Schönberger
177 lines
6.4 KiB
Python
177 lines
6.4 KiB
Python
import numpy as np
|
|
from skimage.morphology import skeletonize, medial_axis
|
|
import numpy.testing
|
|
from skimage import draw
|
|
from scipy.ndimage import correlate
|
|
from skimage.io import imread
|
|
from skimage import data_dir
|
|
import os.path
|
|
|
|
|
|
class TestSkeletonize():
|
|
def test_skeletonize_no_foreground(self):
|
|
im = np.zeros((5, 5))
|
|
result = skeletonize(im)
|
|
numpy.testing.assert_array_equal(result, np.zeros((5, 5)))
|
|
|
|
def test_skeletonize_wrong_dim1(self):
|
|
im = np.zeros((5))
|
|
numpy.testing.assert_raises(ValueError, skeletonize, im)
|
|
|
|
def test_skeletonize_wrong_dim2(self):
|
|
im = np.zeros((5, 5, 5))
|
|
numpy.testing.assert_raises(ValueError, skeletonize, im)
|
|
|
|
def test_skeletonize_not_binary(self):
|
|
im = np.zeros((5, 5))
|
|
im[0, 0] = 1
|
|
im[0, 1] = 2
|
|
numpy.testing.assert_raises(ValueError, skeletonize, im)
|
|
|
|
def test_skeletonize_unexpected_value(self):
|
|
im = np.zeros((5, 5))
|
|
im[0, 0] = 2
|
|
numpy.testing.assert_raises(ValueError, skeletonize, im)
|
|
|
|
def test_skeletonize_all_foreground(self):
|
|
im = np.ones((3, 4))
|
|
skeletonize(im)
|
|
|
|
def test_skeletonize_single_point(self):
|
|
im = np.zeros((5, 5), np.uint8)
|
|
im[3, 3] = 1
|
|
result = skeletonize(im)
|
|
numpy.testing.assert_array_equal(result, im)
|
|
|
|
def test_skeletonize_already_thinned(self):
|
|
im = np.zeros((5, 5), np.uint8)
|
|
im[3, 1:-1] = 1
|
|
im[2, -1] = 1
|
|
im[4, 0] = 1
|
|
result = skeletonize(im)
|
|
numpy.testing.assert_array_equal(result, im)
|
|
|
|
def test_skeletonize_output(self):
|
|
im = imread(os.path.join(data_dir, "bw_text.png"), as_grey=True)
|
|
|
|
# make black the foreground
|
|
im = (im == 0)
|
|
result = skeletonize(im)
|
|
|
|
expected = np.load(os.path.join(data_dir, "bw_text_skeleton.npy"))
|
|
numpy.testing.assert_array_equal(result, expected)
|
|
|
|
def test_skeletonize_num_neighbours(self):
|
|
# an empty image
|
|
image = np.zeros((300, 300))
|
|
|
|
# foreground object 1
|
|
image[10:-10, 10:100] = 1
|
|
image[-100:-10, 10:-10] = 1
|
|
image[10:-10, -100:-10] = 1
|
|
|
|
# foreground object 2
|
|
rs, cs = draw.line(250, 150, 10, 280)
|
|
for i in range(10):
|
|
image[rs + i, cs] = 1
|
|
rs, cs = draw.line(10, 150, 250, 280)
|
|
for i in range(20):
|
|
image[rs + i, cs] = 1
|
|
|
|
# foreground object 3
|
|
ir, ic = np.indices(image.shape)
|
|
circle1 = (ic - 135)**2 + (ir - 150)**2 < 30**2
|
|
circle2 = (ic - 135)**2 + (ir - 150)**2 < 20**2
|
|
image[circle1] = 1
|
|
image[circle2] = 0
|
|
result = skeletonize(image)
|
|
|
|
# there should never be a 2x2 block of foreground pixels in a skeleton
|
|
mask = np.array([[1, 1],
|
|
[1, 1]], np.uint8)
|
|
blocks = correlate(result, mask, mode='constant')
|
|
assert not numpy.any(blocks == 4)
|
|
|
|
def test_lut_fix(self):
|
|
im = np.zeros((6, 6), np.uint8)
|
|
im[1, 2] = 1
|
|
im[2, 2] = 1
|
|
im[2, 3] = 1
|
|
im[3, 3] = 1
|
|
im[3, 4] = 1
|
|
im[4, 4] = 1
|
|
im[4, 5] = 1
|
|
result = skeletonize(im)
|
|
expected = np.array([[0, 0, 0, 0, 0, 0],
|
|
[0, 0, 1, 0, 0, 0],
|
|
[0, 0, 0, 1, 0, 0],
|
|
[0, 0, 0, 0, 1, 0],
|
|
[0, 0, 0, 0, 0, 1],
|
|
[0, 0, 0, 0, 0, 0]], dtype=np.uint8)
|
|
assert np.all(result == expected)
|
|
|
|
|
|
class TestMedialAxis():
|
|
def test_00_00_zeros(self):
|
|
'''Test skeletonize on an array of all zeros'''
|
|
result = medial_axis(np.zeros((10, 10), bool))
|
|
assert np.all(result == False)
|
|
|
|
def test_00_01_zeros_masked(self):
|
|
'''Test skeletonize on an array that is completely masked'''
|
|
result = medial_axis(np.zeros((10, 10), bool),
|
|
np.zeros((10, 10), bool))
|
|
assert np.all(result == False)
|
|
|
|
def test_01_01_rectangle(self):
|
|
'''Test skeletonize on a rectangle'''
|
|
image = np.zeros((9, 15), bool)
|
|
image[1:-1, 1:-1] = True
|
|
#
|
|
# The result should be four diagonals from the
|
|
# corners, meeting in a horizontal line
|
|
#
|
|
expected = np.array([[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
|
|
[0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0],
|
|
[0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0],
|
|
[0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0],
|
|
[0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0],
|
|
[0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0],
|
|
[0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0],
|
|
[0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0],
|
|
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]],
|
|
bool)
|
|
result = medial_axis(image)
|
|
assert np.all(result == expected)
|
|
result, distance = medial_axis(image, return_distance=True)
|
|
assert distance.max() == 4
|
|
|
|
def test_01_02_hole(self):
|
|
'''Test skeletonize on a rectangle with a hole in the middle'''
|
|
image = np.zeros((9, 15), bool)
|
|
image[1:-1, 1:-1] = True
|
|
image[4, 4:-4] = False
|
|
expected = np.array([[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
|
|
[0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0],
|
|
[0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0],
|
|
[0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0],
|
|
[0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0],
|
|
[0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0],
|
|
[0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0],
|
|
[0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0],
|
|
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]],
|
|
bool)
|
|
result = medial_axis(image)
|
|
assert np.all(result == expected)
|
|
|
|
def test_narrow_image(self):
|
|
"""Test skeletonize on a 1-pixel thin strip"""
|
|
image = np.zeros((1, 5), bool)
|
|
image[:, 1:-1] = True
|
|
result = medial_axis(image)
|
|
assert np.all(result == image)
|
|
|
|
|
|
if __name__ == '__main__':
|
|
np.testing.run_module_suite()
|