def construct_from_parameters_add_tree(): parameters = Skeleton.define_parameters('2017-01-12_FD0156-2', (11.24, 11.24, 32)) skel = Skeleton(parameters=parameters) nodes = skel.define_nodes([40000, 40100, 40200], [45000, 45100, 45200], [1000, 1100, 1200], [1, 2, 3]) edges = [(1, 3), (2, 3)] skel.add_tree(nodes, edges) skel.write_nml('testdata/PA_gen.nml')
nodes = skel.define_nodes( position_x=positions[:, 0].tolist(), position_y=positions[:, 1].tolist(), position_z=positions[:, 2].tolist(), id=list(range(min_id, max_id + 1)) ) edges = np.array([ [1, 2], [1, 3], [1, 4], [1, 5], [2, 3], [2, 4], [3, 5], [4, 6] ]) + skel.max_node_id() if zi not in skel.group_ids: skel.add_group(id=zi, name='plane{:03d}:X'.format(zi)) skel.add_tree( nodes=nodes, edges=edges, name='patch{:02d}.{:02d}:X'.format(xi, yi), group_id=zi) min_id = max_id + 1 #skel.write_nml(path_nml_out)