hofig4s4
Sequential thinning with a single hit-or-miss template

from morph import *

S=mmbinary([
    [0,0,0,0,0,0,0,0],
    [0,1,1,1,1,1,0,0],
    [0,1,1,1,1,1,1,0],
    [0,1,1,1,1,1,1,0],
    [0,0,0,0,1,1,1,0],
    [0,0,0,0,1,1,1,0],
    [0,0,0,0,1,1,1,0],
    [0,0,0,0,0,0,0,0]])

B1=mmbinary([[0,1,1],
             [0,1,1],
             [0,0,0]])

B2=mmbinary([[1,1,1],
             [0,1,1],
             [0,0,1]])

I=mmse2interval(mmimg2se(B1),mmimg2se(B2))

S1=mmsubm(S,mmsupgen(S,I))

S2=mmsubm(S1,mmsupgen(S1,I))

S3=mmsubm(S2,mmsupgen(S2,I))

S4=mmsubm(S3,mmsupgen(S3,I))

S5=mmsubm(S4,mmsupgen(S4,I))

fig_a = mmbshow(mmunion(B1,mmneg(B2)),B1)

fig_b = mmbshow(S>=0,S)

Warning: Converting input image from int32 to binary uint8.

fig_c = mmbshow(S>=0,S1)

Warning: Converting input image from int32 to binary uint8.

fig_d = mmbshow(S>=0,S2)

Warning: Converting input image from int32 to binary uint8.

fig_e = mmbshow(S>=0,S3)

Warning: Converting input image from int32 to binary uint8.

fig_f = mmbshow(S>=0,S4)

Warning: Converting input image from int32 to binary uint8.

fig_g = mmbshow(S>=0,S5)

Warning: Converting input image from int32 to binary uint8.

mmshow(fig_a)

mmshow(fig_b)

mmshow(fig_c)

mmshow(fig_d)

mmshow(fig_e)

mmshow(fig_f)

mmshow(fig_g)