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)