我尝试在我的橙色虚线上显示以下点[control = [0.20,0.40,0.60,0.80]]。这些点表示橙色虚线的每个分段的一段内的位置。根据下图,我做了一个注释,它在旁边显示了它们,但我需要它们在上面作为我做的红色x,但对于所有的小节。如何在橙色虚线上突出显示所有这些点?

import numpy as np
import scipy.special
import matplotlib.pyplot as plt
def calc_bezier_path(control_points, n_points=100):
"""
Compute bezier path (trajectory) given control points.
:param control_points: (numpy array)
:param n_points: (int) number of points in the trajectory
:return: (numpy array)
"""
traj = []
for t in np.linspace(0, 1, n_points):
traj.append(bezier(t, control_points))
return np.array(traj)
def bernstein_poly(n, i, t):
"""
Bernstein polynom.
:param n: (int) polynom degree
:param i: (int)
:param t: (float)
:return: (float)
"""
return scipy.special.comb(n, i) * t ** i * (1 - t) ** (n - i)
def bezier(t, control_points):
"""
Return one point on the bezier curve.
:param t: (float) number in [0, 1]
:param control_points: (numpy array)
:return: (numpy array) Coordinates of the point
"""
n = len(control_points) - 1
return np.sum([bernstein_poly(n, i, t) * control_points[i] for i in range(n + 1)], axis=0)
def line_bezier(visx, visy, control, mod="nothing"):
vis = np.column_stack((visx,visy))
path_x, path_y = np.array([]),np.array([])
setting = {"nothing":[len(vis)-2, 1, 1], "start":[len(vis)-1, 0, 0], "end":[len(vis)-1, 1, 0], "both":[len(vis), 0, -1]}
epoch = setting[mod][0]
start = setting[mod][1]
end = setting[mod][2]
if len(vis) > 2:
current_control = vis[0]
for x in range(epoch):
if x != (epoch-1):
for y in control:
if y == control[0]:
mid_control = [(vis[x+start,0]+(vis[x+(start+1),0]-vis[x+start,0])*y), (vis[x+start,1]+(vis[x+(start+1),1]-vis[x+start,1])*y)]
plt.annotate("*", mid_control)
bezier_line = calc_bezier_path(np.array([current_control,vis[x+start], mid_control]))
path_x = np.append(path_x, bezier_line.T[0])
path_y = np.append(path_y, bezier_line.T[1])
current_control = mid_control
else:
mid_control = [(vis[x+start,0]+(vis[x+(start+1),0]-vis[x+start,0])*y), (vis[x+start,1]+(vis[x+(start+1),1]-vis[x+start,1])*y)]
plt.annotate("*", mid_control)
bezier_line = calc_bezier_path(np.array([current_control, mid_control]))
path_x = np.append(path_x, bezier_line.T[0])
path_y = np.append(path_y, bezier_line.T[1])
current_control = mid_control
else:
if mod == "end" or mod == "both":
bezier_line = calc_bezier_path(np.array([current_control, vis[x+(end+1)]]))
else:
bezier_line = calc_bezier_path(np.array([current_control, vis[x+end], vis[x+(end+1)]]))
path_x = np.append(path_x, bezier_line.T[0])
path_y = np.append(path_y, bezier_line.T[1])
else:
path_x, path_y = visx, visy
return path_x, path_y
visx, visy = [1,2,10,15,20,25,21], [0,5,1,4,2,3,3]
control = [0.20,0.40,0.60,0.80]
path_x, path_y = line_bezier(visx, visy, control,mod="end")
plt.plot(path_x, path_y)
plt.plot(visx, visy, "--o")
for xy in range(len(visx)):
plt.annotate(f"P{xy}", [visx[xy], visy[xy]])
plt.xlabel('X')
plt.ylabel('Y')
plt.show()发布于 2020-06-23 05:49:02
更改:
for y in control:
if y == control[0]:
mid_control = [(vis[x+start,0]+(vis[x+(start+1),0]-vis[x+start,0])*y), (vis[x+start,1]+(vis[x+(start+1),1]-vis[x+start,1])*y)]
plt.annotate("*", mid_control)
bezier_line = calc_bezier_path(np.array([current_control,vis[x+start], mid_control]))
path_x = np.append(path_x, bezier_line.T[0])
path_y = np.append(path_y, bezier_line.T[1])
current_control = mid_control
else:
mid_control = [(vis[x+start,0]+(vis[x+(start+1),0]-vis[x+start,0])*y), (vis[x+start,1]+(vis[x+(start+1),1]-vis[x+start,1])*y)]
plt.annotate("*", mid_control)
bezier_line = calc_bezier_path(np.array([current_control, mid_control]))
path_x = np.append(path_x, bezier_line.T[0])
path_y = np.append(path_y, bezier_line.T[1])
current_control = mid_control至:
for y in control:
mid_control = [(vis[x+start,0]+(vis[x+(start+1),0]-vis[x+start,0])*y), (vis[x+start,1]+(vis[x+(start+1),1]-vis[x+start,1])*y)]
mid_control2 = [(vis[x+start,0]+(vis[x+(start+1),0]-vis[x+start,0])*y-0.1), (vis[x+start,1]+(vis[x+(start+1),1]-vis[x+start,1])*y-0.1)]
if y == control[0]:
plt.annotate("*", mid_control2)
bezier_line = calc_bezier_path(np.array([current_control,vis[x+start], mid_control]))
path_x = np.append(path_x, bezier_line.T[0])
path_y = np.append(path_y, bezier_line.T[1])
current_control = mid_control
else:
plt.annotate("*", mid_control2)
bezier_line = calc_bezier_path(np.array([current_control, mid_control]))
path_x = np.append(path_x, bezier_line.T[0])
path_y = np.append(path_y, bezier_line.T[1])
current_control = mid_control总而言之:
import numpy as np
import scipy.special
import matplotlib.pyplot as plt
def calc_bezier_path(control_points, n_points=100):
"""
Compute bezier path (trajectory) given control points.
:param control_points: (numpy array)
:param n_points: (int) number of points in the trajectory
:return: (numpy array)
"""
traj = []
for t in np.linspace(0, 1, n_points):
traj.append(bezier(t, control_points))
return np.array(traj)
def bernstein_poly(n, i, t):
"""
Bernstein polynom.
:param n: (int) polynom degree
:param i: (int)
:param t: (float)
:return: (float)
"""
return scipy.special.comb(n, i) * t ** i * (1 - t) ** (n - i)
def bezier(t, control_points):
"""
Return one point on the bezier curve.
:param t: (float) number in [0, 1]
:param control_points: (numpy array)
:return: (numpy array) Coordinates of the point
"""
n = len(control_points) - 1
return np.sum([bernstein_poly(n, i, t) * control_points[i] for i in range(n + 1)], axis=0)
def line_bezier(visx, visy, control, mod="nothing"):
vis = np.column_stack((visx,visy))
path_x, path_y = np.array([]),np.array([])
setting = {"nothing":[len(vis)-2, 1, 1], "start":[len(vis)-1, 0, 0], "end":[len(vis)-1, 1, 0], "both":[len(vis), 0, -1]}
epoch = setting[mod][0]
start = setting[mod][1]
end = setting[mod][2]
if len(vis) > 2:
current_control = vis[0]
for x in range(epoch):
if x != (epoch-1):
for y in control:
mid_control = [(vis[x+start,0]+(vis[x+(start+1),0]-vis[x+start,0])*y), (vis[x+start,1]+(vis[x+(start+1),1]-vis[x+start,1])*y)]
mid_control2 = [(vis[x+start,0]+(vis[x+(start+1),0]-vis[x+start,0])*y-0.1), (vis[x+start,1]+(vis[x+(start+1),1]-vis[x+start,1])*y-0.1)]
if y == control[0]:
plt.annotate("*", mid_control2)
bezier_line = calc_bezier_path(np.array([current_control,vis[x+start], mid_control]))
path_x = np.append(path_x, bezier_line.T[0])
path_y = np.append(path_y, bezier_line.T[1])
current_control = mid_control
else:
plt.annotate("*", mid_control2)
bezier_line = calc_bezier_path(np.array([current_control, mid_control]))
path_x = np.append(path_x, bezier_line.T[0])
path_y = np.append(path_y, bezier_line.T[1])
current_control = mid_control
else:
if mod == "end" or mod == "both":
bezier_line = calc_bezier_path(np.array([current_control, vis[x+(end+1)]]))
else:
bezier_line = calc_bezier_path(np.array([current_control, vis[x+end], vis[x+(end+1)]]))
path_x = np.append(path_x, bezier_line.T[0])
path_y = np.append(path_y, bezier_line.T[1])
else:
path_x, path_y = visx, visy
return path_x, path_y
visx, visy = [1,2,10,15,20,25,21], [0,5,1,4,2,3,3]
control = [0.20,0.40,0.60,0.80]
path_x, path_y = line_bezier(visx, visy, control,mod="end")
plt.plot(path_x, path_y)
plt.plot(visx, visy, "--o")
for xy in range(len(visx)):
plt.annotate(f"P{xy}", [visx[xy], visy[xy]])
plt.xlabel('X')
plt.ylabel('Y')
plt.show()输出:

您看,在您的if else语句中,每个块中都有mid_control = [...],这不是必需的,您可以在if else语句之外定义mid_control = [...]。
在每个条件块(if y == control[0]和else)中,您使用了两次变量mid_control,一次用于"*",一次用于calc_bezier_path()。
正如您所知道的,calc_bezier_path()的绘图是准确的,而"*"的绘图是向右移动的。
如果我修改mid_control,"*"图将被纠正,但calc_bezier_path()图将被破坏。
因此,我创建了另一个变量,将mid_control赋给它,并为"*"图修改了一个变量,为calc_bezier_path()图保留了另一个变量。
https://stackoverflow.com/questions/62523906
复制相似问题