In [2]:

    

from cremi import Annotations, Volume
from cremi.io import CremiFile
import numpy as np


    

In [12]:

    




    

In [4]:

    

import mahotas as mh


    

In [5]:

    

seg = np.zeros((125,1250,1250), dtype=np.uint64)
for z in range(125):
    
    i = mh.imread('/tmp/creminew/test/Aplus/raw/0/0/0/neuroproof/000000000_000000000_000000'+str(z).zfill(3)+'_neuroproof.png')
    
    rhoana_single = np.zeros((i.shape[0], i.shape[1]), dtype=np.uint64)
    rhoana_single[:, :] = i[:,:,0]*256*256 + i[:,:,1]*256 + i[:,:,2]
    
    seg[z] = rhoana_single


    

In [16]:

    

segB = np.zeros((125,1250,1250), dtype=np.uint64)
for z in range(125):
    
    i = mh.imread('/tmp/creminewB/test/Bplus/raw/0/0/0/neuroproof/000000000_000000000_000000'+str(z).zfill(3)+'_neuroproof.png')
    
    rhoana_single = np.zeros((i.shape[0], i.shape[1]), dtype=np.uint64)
    rhoana_single[:, :] = i[:,:,0]*256*256 + i[:,:,1]*256 + i[:,:,2]
    
    segB[z] = rhoana_single


    

In [19]:

    

segC = np.zeros((125,1250,1250), dtype=np.uint64)
for z in range(125):
    
    i = mh.imread('/tmp/creminewC/test/Cplus/raw/0/0/0/neuroproof/000000000_000000000_000000'+str(z).zfill(3)+'_neuroproof.png')
    
    rhoana_single = np.zeros((i.shape[0], i.shape[1]), dtype=np.uint64)
    rhoana_single[:, :] = i[:,:,0]*256*256 + i[:,:,1]*256 + i[:,:,2]
    
    segC[z] = rhoana_single


    

In [13]:

    

file = CremiFile("/home/d/data/cremiA+.hdf", "w")

file.write_neuron_ids(Volume(seg, resolution=(40.0,4.0,4.0)))

file.close()


    

In [18]:

    

file = CremiFile("/home/d/data/cremiB+.hdf", "w")

file.write_neuron_ids(Volume(segB, resolution=(40.0,4.0,4.0)))

file.close()


    

In [20]:

    

file = CremiFile("/home/d/data/cremiC+.hdf", "w")

file.write_neuron_ids(Volume(segC, resolution=(40.0,4.0,4.0)))

file.close()


    

In [15]:

    

print 'a'


    

    




a

In [ ]: