In [1]:

    

%run common.ipynb


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

import glob
path = "/Users/arve/Dokumenter/TFY4500/well z-scan/slide--S00/chamber--U00--V00/field--X01--Y01/"
files = glob.glob(path + "*.ome.tif")


    

In [3]:

    

img = pylab.imread(files[0])
for f in files:
    tmp = pylab.imread(f)
    img = numpy.maximum(tmp, img)


    

In [4]:

    

imshow(img)


    

    
Out[4]:





<matplotlib.image.AxesImage at 0x110d93290>






    

In [5]:

    

imsave(path+"max.tif", img)


    

In [6]:

    

imsave(path+"max-gray.tif", img[:,:,1])


    

In [7]:

    

img = pylab.imread(files[0])
for i, f in enumerate(files):
    if i%3 == 0:
        tmp = pylab.imread(f)
        img = numpy.maximum(tmp, img)
imshow(img)


    

    
Out[7]:





<matplotlib.image.AxesImage at 0x110e64810>






    

In [8]:

    

imsave(path+"max3.tif", img)


    

In [10]:

    

img = smooth(pylab.imread(files[0])[:,:,1])
for f in files:
    tmp = smooth(pylab.imread(f)[:,:,1])
    img = numpy.maximum(tmp, img)


    

In [11]:

    

imsave(path+"1smooth max.tif", img)