The code snippet below shows a simple scenario for geometry analysis of the reconstructed cytoskelton.
from cytoskeleton_analyser.inout import Paths # simulation-specific paths
from cytoskeleton_analyser.position.drivers import process # driver fuction
from sim_specs import cell, data_in, data_out, rinds
# Set general parameters of the input simulation. (a shortlist
# of config settings).
params = {
'edge_len_fine': 0.008, # (μm) length of filament edge.
'iscoarse': True, # use coarse-grained representation
'coarsegraining': 25, # coarse-graining factor
'use_final': False, # if True, limit to final snapshot only
'cell': cell, # cell object
'slice_limits': {'bottom': 0., 'top': 0.8}, # z-pos. of slicing planes
}
# Specify features to process. If this is None,
# all features are included.
features = [
'RadialMass',
]
# Decide if this is a new analysis, or we merely want to import one.
new = True
# This dict will become populated only if analysis
# is imported ('new' = False).
# Otherwise, the results are stored on a disc.
result = {}
if __name__ == '__main__':
for ci in rinds:
paths = Paths(data_in, data_out, params['cell'], ci)
result[str(ci)] = process(new, paths, ci, params, True, features)