from morph import *
from Numeric import sum
f1=mmneg(mmbinary(mmreadgray('toner1.tif'),128));
f2=mmneg(mmbinary(mmreadgray('toner2.tif'),128));
f3=mmneg(mmbinary(mmreadgray('toner3.tif'),128));
mmshow(f1);
mmshow(f2);
mmshow(f3);
g1=mmopentransf(f1);
processing r= 0
processing r= 1
processing r= 2
processing r= 3
processing r= 4
processing r= 5
processing r= 6
processing r= 7
processing r= 8
processing r= 9
processing r= 10
processing r= 11
processing r= 12
processing r= 13
processing r= 14
processing r= 15
processing r= 16
processing r= 17
processing r= 18
processing r= 19
g2=mmopentransf(f2);
processing r= 0
processing r= 1
processing r= 2
processing r= 3
processing r= 4
processing r= 5
processing r= 6
processing r= 7
processing r= 8
processing r= 9
processing r= 10
processing r= 11
processing r= 12
processing r= 13
processing r= 14
processing r= 15
processing r= 16
processing r= 17
processing r= 18
processing r= 19
processing r= 20
processing r= 21
processing r= 22
processing r= 23
processing r= 24
processing r= 25
processing r= 26
processing r= 27
processing r= 28
processing r= 29
processing r= 30
processing r= 31
g3=mmopentransf(f3);
processing r= 0
processing r= 1
processing r= 2
processing r= 3
processing r= 4
processing r= 5
processing r= 6
processing r= 7
processing r= 8
processing r= 9
processing r= 10
processing r= 11
processing r= 12
processing r= 13
processing r= 14
processing r= 15
processing r= 16
processing r= 17
processing r= 18
processing r= 19
processing r= 20
processing r= 21
processing r= 22
processing r= 23
processing r= 24
processing r= 25
print mmstats(g1,'mean')
2.58168587022
print mmstats(g2,'mean')
4.12531635335
print mmstats(g3,'mean')
5.03731958763
h1=mmhistogram(g1);
h1=h1[1:];
h2=mmhistogram(g2);
h2=h2[1:];
h3=mmhistogram(g3);
h3=h3[1:];
mmplot([[(h1+0.)/sum(h1)]],[['style','histeps']])
[Errno 22] Invalid argument
mmplot([[(h2+0.)/sum(h2)]],[['style','histeps']])
[Errno 22] Invalid argument
mmplot([[(h3+0.)/sum(h3)]],[['style','histeps']])
[Errno 22] Invalid argument
dmax=uint8(mmunion(mmstats(g1,'max'),mmstats(g2,'max'),mmstats(g3,'max')));
Warning: Converting input image from double to uint16.
Warning: Converting input image from double to uint16.
Warning: Converting input image from double to uint16.
g1[0,0]=dmax[0];
g2[0,0]=dmax[0];
g3[0,0]=dmax[0]; # trick for displaying all three at same enhancement
mmshow(g1);
mmshow(g2);
mmshow(g3);
