Fig. 2.5:
.
from morph import *
B=mmsedisk(10)
F=mmpad(mmreadgray('numbers.tif'),B)
E=mmclose(F,B)
O=mmcloseth(F,B)
mmshow(F)
mmshow(E)
mmshow(O)