hofig6s15
Flat zones

from morph import *

from Numeric import zeros

f=mmbinary(zeros((100,100)))

f[50,50]=1

d=mmdist(mmneg(f),mmsebox(),'EUCLIDEAN');

g1=mmthreshad(d,0,45);

Warning: Converting input image from int32 to uint16.
Warning: Converting input image from int32 to uint16.

g2=mmthreshad(d,0,30);

Warning: Converting input image from int32 to uint16.
Warning: Converting input image from int32 to uint16.

g3=mmdil(f,mmsesum(mmimg2se([[0,1,0],
                             [0,1,0],
                             [1,1,1]]),17));

Warning: Converting input image from int32 to binary uint8.

g=mmaddm(mmgray(mmsymdif(g1,g3),'uint8',1),
         mmgray(g2,'uint8',2));

g_flat = mmlabelflat(g,mmsebox())

mmshow(g)

mmlblshow(g_flat)

f=mmreadgray('pcb-gray.tif');

f=f[:60,-60:];

f_flat=mmlabelflat(f,mmsebox());

mmshow(f)

mmlblshow(f_flat)