def load_particles_from_spatiocyte(filename, index=0, ps=None): if not os.path.isfile(filename): return ps try: reader = SpatiocyteLogReader(filename) if ps is None: header = reader.getHeader() col_size = header['aColSize'] row_size = header['aRowSize'] layer_size = header['aLayerSize'] lspecies = header['latticeSpecies'] voxel_radius = header['aVoxelRadius'] ps = LatticeParticleSpace(col_size, row_size, layer_size, lspecies, voxel_radius) species = reader.skipSpeciesTo(index) molecules = species['Molecules']; #debug("index : %d" % index) for sp in molecules: sid = sp['index'] for coord in sp['Coords']: ps.add_particle(LatticeParticle(sid, coord)) finally: reader.close() return ps
def main(): ''' main function ''' if (len(sys.argv) != 2): print "USAGE: python %s filename" % sys.argv[0] return name = sys.argv[1] if (not os.path.exists(name)): print "%s does not exist!" % name return reader = SpatiocyteLogReader(name) header = reader.header print '[header]\t', header indexSize = reader.getIndexSize() index = indexSize - 1 print "LastSeek : %d, FooterSeek : %d" % (reader.tell(), reader.footerSeek) print "Index : %d" % index species = reader.skipSpeciesTo(index) for spiece in species['Molecules']: print spiece reader.close()
def load_particles_from_spatiocyte(filename, index=0, ps=None): if not os.path.isfile(filename): return ps try: reader = SpatiocyteLogReader(filename) if ps is None: header = reader.getHeader() col_size = header['aColSize'] row_size = header['aRowSize'] layer_size = header['aLayerSize'] lspecies = header['latticeSpecies'] voxel_radius = header['aVoxelRadius'] ps = LatticeParticleSpace(col_size, row_size, layer_size, lspecies, voxel_radius) species = reader.skipSpeciesTo(index) molecules = species['Molecules'] #debug("index : %d" % index) for sp in molecules: sid = sp['index'] for coord in sp['Coords']: ps.add_particle(LatticeParticle(sid, coord)) finally: reader.close() return ps
def main(): ''' main function ''' if (len(sys.argv) != 2): print "USAGE: python %s filename" % sys.argv[0] return name = sys.argv[1] if (not os.path.exists(name)): print "%s does not exist!" % name return reader = SpatiocyteLogReader(name) header = reader.header print '[header]\t',header indexSize = reader.getIndexSize() index = indexSize - 1 print "LastSeek : %d, FooterSeek : %d" % (reader.tell(), reader.footerSeek) print "Index : %d" % index species = reader.skipSpeciesTo(index) for spiece in species['Molecules']: print spiece reader.close()