from morph import *
from Numeric import shape
C=mmreadgray('club.tif')
A=mmsubm(C,C)
A[shape(C)[0]/2:,shape(C)[1]/2]=1
A[shape(C)[0]/2,shape(C)[1]/2:]=1
E=mmsupcanon(mmintersec(A,C),mmendpoints())
A=mmdil(A)
A=mmintersec(A,C)
B=mmsedisk(10)
D=mmcdil(A,C,B)
Ed=mmgradm(mmdil(E,B))
mmshow(A)
mmshow(mmseshow(B))
mmshow(C)
mmshow(D)
mmshow(C,D,A,Ed)