hofig3s13
Reconstruction using geodesic dilation

from morph import *

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

M=mmintersec(A,0)

Warning: Converting input image from int32 to binary uint8.

M[2,3]=1

M1=mmcdil(M,A,mmsebox(),1)

M2=mmcdil(M,A,mmsebox(),2)

M3=mmcdil(M,A,mmsebox(),3)

M4=mmcdil(M,A,mmsebox(),4)

M5=mmcdil(M,A,mmsebox(),5)

mmshow(mmbshow(A,M))

mmshow(mmbshow(A,M1))

mmshow(mmbshow(A,M2))

mmshow(mmbshow(A,M3))

mmshow(mmbshow(A,M4))

mmshow(mmbshow(A,M5))