from morph import *
f=mmreadgray('chromosomes.tif')
L=30
r0 = mmopen(f,mmseline(L,0))
r45 = mmopen(f,mmseline(L,45))
r90 = mmopen(f,mmseline(L,90))
rm45 = mmopen(f,mmseline(L,-45))
r0_rec = mminfrec(r0,f)
r45_rec = mminfrec(r45,f)
r90_rec = mminfrec(r90,f)
rm45_rec = mminfrec(rm45,f)
g = mmintersec(r0_rec,r45_rec,r90_rec,rm45_rec)
fig_a = f
fig_b = g
# mmshow([fig3_11_3_1a, fig3_11_3_1b]);
fig_c =mmunion(mmgray(r0),mmgray(r0_rec,'uint8',140))
fig_d =mmunion(mmgray(r90),mmgray(r90_rec,'uint8',140))
fig_e =mmunion(mmgray(r45),mmgray(r45_rec,'uint8',140))
fig_f =mmunion(mmgray(rm45),mmgray(rm45_rec,'uint8',140))
mmshow(fig_a)
mmshow(fig_b)
mmshow(fig_c)
mmshow(fig_d)
mmshow(fig_e)
mmshow(fig_f)
