def mesh(self, leaf_key, L_shape, Lw_shape, length, s_base, s_top, incline=1, flipx = False):
""" Compute mesh for a leaf element.
- shape is a x,y,s,r tuple descriibing leaf shape
- L_shape is the length of the scaled shape
- Lw_shape is the width of the scaled shape
- length is the total visible length to be meshed
- s_base and s_top are relative proportion (on length) of the element to represent
"""
shape = self.get_leaf(leaf_key)
shape = incline_leaf(shape, incline)
if flipx:
shape = (-shape[0],)+shape[1:] # to position leaves along tiller emitted
leaf_mesh = fitting.mesh4(shape, L_shape, length, s_base, s_top, Lw_shape, twist=self.twist)
if leaf_mesh:
pts, ind = leaf_mesh
if len(ind) < 1:
#raise AdelError('ERROR less than 1 triangles') # mesh4 filters triangle < 1e-6
mesh = None
else:
mesh = fitting.plantgl_shape(pts, ind)
else:
if length > 0:
print 'ERROR No mesh', s_base, s_top, length
pass
mesh = None
return mesh
def leaf_to_mesh_new(leaf, lmax, l, rmax, twist=False, nb_twist=1.,
nb_waves=8):
pts, ind = fitting._mesh(leaf, lmax, l, rmax,
functor=fitting.leaf_to_mesh_new,
twist=twist, nb_twist=nb_twist, nb_waves=nb_waves)
mesh = fitting.plantgl_shape(pts, ind)
return mesh
def mesh(self, leaf_key, L_shape, Lw_shape, length, s_base, s_top, incline=1, flipx = False):
""" Compute mesh for a leaf element.
- shape is a x,y,s,r tuple descriibing leaf shape
- L_shape is the length of the scaled shape
- Lw_shape is the width of the scaled shape
- length is the total visible length to be meshed
- s_base and s_top are relative proportion (on length) of the element to represent
"""
shape = self.get_leaf(leaf_key)
shape = incline_leaf(shape, incline)
if flipx:
shape = (-shape[0],)+shape[1:] # to position leaves along tiller emitted
leaf_mesh = fitting.mesh4(shape, L_shape, length, s_base, s_top, Lw_shape, twist=self.twist)
if leaf_mesh:
pts, ind = leaf_mesh
if len(ind) < 1:
#raise AdelError('ERROR less than 1 triangles') # mesh4 filters triangle < 1e-6
mesh = None
else:
mesh = fitting.plantgl_shape(pts, ind)
else:
if length > 0:
print 'ERROR No mesh', s_base, s_top, length
pass
mesh = None
return mesh
def leaf_to_mesh_cicloid_twist(leaf,
lmax,
l,
rmax,
twist=True,
nb_twist=2,
nb_waves=8):
pts, ind = fitting.mesh3(leaf, lmax, l, rmax)
mesh = fitting.plantgl_shape(pts, ind)
return mesh
def test7(leaf=leaf, scene=None):
if scene is None:
scene = pgl.Scene()
Viewer.display(scene)
x, y, s, r = leaf
spline_leaf, leaf_surface = fitting.fit_leaf(x, y, s, r)
pts, ind = fitting.mesh(spline_leaf, 30, 7, 7, 1)
fitting.write_smf('leaf_full.smf', pts, ind)
Viewer.display(fitting.plantgl_shape(pts, ind))
def leaf_to_mesh_new(leaf,
lmax,
l,
rmax,
twist=False,
nb_twist=1.,
nb_waves=8):
pts, ind = fitting._mesh(leaf,
lmax,
l,
rmax,
functor=fitting.leaf_to_mesh_new,
twist=twist,
nb_twist=nb_twist,
nb_waves=nb_waves)
mesh = fitting.plantgl_shape(pts, ind)
return mesh
def test8(leaf=leaf, scene=None):
global translation, zt
if scene is None:
scene = pgl.Scene()
Viewer.display(scene)
x, y, s, r = leaf
leaf_new, leaf_surface = fitting.fit2(x, y, s, r)
pts, ind = fitting.mesh2(leaf_new, 7, 7, 1)
#pts2, ind2 = fitting.qslim(13, pts, ind)
#mesh = fitting.plantgl_shape(pts2, ind2)
#sc=pgl.SurfComputer(pgl.Discretizer())
#mesh.apply(sc)
#scale_z = leaf_surface*7 / (sc.surface)
#mesh_final = mesh.transform(pgl.Scaling((1,1,scale_z)))
#mesh_final = mesh
scene += pgl.Translated(translation, fitting.plantgl_shape(pts, ind))
#scene += pgl.Translated(translation+(0,yt/3.,0), mesh_final)
Viewer.update()
def leaf_to_mesh(leaf, lmax, l, rmax, antisens):
pts, ind = fitting.mesh3(leaf, lmax, l, rmax, antisens)
mesh = fitting.plantgl_shape(pts, ind)
return mesh
def leaf_element(leaf, lmax, l, s_base, s_top, rmax):
pts, ind = fitting.mesh4(leaf, lmax, l, s_base, s_top, rmax)
mesh = fitting.plantgl_shape(pts, ind)
return mesh
def leaf_to_mesh(leaf, lmax, l, rmax, antisens):
pts, ind = fitting.mesh3(leaf, lmax, l, rmax, antisens)
mesh = fitting.plantgl_shape(pts, ind)
return mesh
def leaf_element(leaf, lmax, l, s_base, s_top, rmax):
pts, ind = fitting.mesh4(leaf, lmax, l, s_base, s_top, rmax)
mesh = fitting.plantgl_shape(pts, ind)
return mesh
def leaf_to_mesh_cicloid_twist(leaf, lmax, l, rmax,
twist=True, nb_twist=2, nb_waves=8):
pts, ind = fitting.mesh3(leaf, lmax, l, rmax)
mesh = fitting.plantgl_shape(pts, ind)
return mesh