I want to draw a $3$-dimensional bezier path with thickness=5 instead of the default value 2, using the code

[[(0,0,0),(1,0,0),(0,1,0),(0,1,1)]]}

curve = bezier3d(path, thickness=5, color='blue')}

recommended in the manual. It seems however, that the thickness does not change.

What do you suggest?

edit retag close merge delete

Sort by » oldest newest most voted There exists at least two options to control thickness for lines and curves in 3D: thickness and radius. If the latter is given, lines are rendered as tubes. Both options work for functions like line3d or parametric_plot3d. However, although these functions are used in the implementation of bezier3d, the radius option has no effect. I haven’t found a simple workaround. But one can always adapt the code of bezier3d. Consider, for example, the following code:

from sage.misc.decorators import options, rename_keyword
@rename_keyword(alpha='opacity')
@options(opacity=1, color="blue", aspect_ratio=[1,1,1], thickness=2, radius=None, plot_points='automatic')
def my_bezier3d(path, **options):
from sage.plot.plot3d import parametric_plot3d as P3D
from sage.modules.free_module_element import vector
from sage.symbolic.ring import SR

plot_options = {'color': options['color'], 'aspect_ratio': options['aspect_ratio'],
'thickness': options['thickness'], 'opacity': options['opacity']}
line_options = {'color': options['color'], 'thickness': options['thickness'],
'opacity': options['opacity']}
if options['plot_points']=='automatic': plot_options.update({'plot_points': 75})
p0 = vector(path[-1])
t = SR.var('t')
if len(path) > 2:
B = (1-t)**3*vector(path)+3*t*(1-t)**2*vector(path)+3*t**2*(1-t)*vector(path[-2])+t**3*p0
G = P3D.parametric_plot3d(list(B), (0, 1), **plot_options)
else:
G = line3d([path, p0], **line_options)

for curve in path[1:]:
if len(curve) > 1:
p1 = vector(curve)
p2 = vector(curve[-2])
p3 = vector(curve[-1])
B = (1-t)**3*p0+3*t*(1-t)**2*p1+3*t**2*(1-t)*p2+t**3*p3
G += P3D.parametric_plot3d(list(B), (0, 1), **plot_options)
else:
G += line3d([p0,curve], **line_options)
p0 = curve[-1]
return G


This code defines my_bezier3d, which has exactly the same syntax as bezier3d, but it admits two additional keywords: radius and plot_points. So, we can now use it:

path_1 = [[(0,0,0),(1,0,0),(1,1,0)], [(0,1,0),(0,0,1),(1,1,1)], [(0,1,1)]]
path_2 = [[(1,0,1),(0.1,0.9,0.7),(0.1,0.9,0.3),(1,0,0)]]
path_3 = [[(0,1,0),(0.7,0.9,0.2),(1,0.8,0.4),(0.6,0.4,0.6),(0.5,0,1)], [(0,0,1),(0,0,0.5),(0.5,0.2,0)]]
curve_1 = my_bezier3d(path_1)
curve_3 = my_bezier3d(path_3, color='red', radius=0.01, plot_points=25)
show(curve_1+curve_2+curve_3, viewer='threejs')


See the resulting graphic in this SageMath Cell.

more