scipy.spatial.transform.

RotationSpline

class scipy.spatial.transform.RotationSpline(times, rotations)

以连续角速度和加速度插值旋转。

每个连续方向之间的旋转向量是时间的立方函数，并且保证角速度和加速度是连续的。这种插值类似于立方样条插值。

请参考 [1] 获取数学和实现细节。

参数：

timesarray_like, shape (N,)

已知旋转的时间。必须指定至少 2 个时间。

rotationsRotation 实例

用于执行插值的旋转。必须包含 N 个旋转。

参考文献

[1]

平滑姿态插值

示例

在您的浏览器中尝试！

>>> from scipy.spatial.transform import Rotation, RotationSpline
>>> import numpy as np

从欧拉角定义时间和旋转序列

>>> times = [0, 10, 20, 40]
>>> angles = [[-10, 20, 30], [0, 15, 40], [-30, 45, 30], [20, 45, 90]]
>>> rotations = Rotation.from_euler('XYZ', angles, degrees=True)

创建插值器对象

>>> spline = RotationSpline(times, rotations)

插值欧拉角、角速度和加速度

>>> angular_rate = np.rad2deg(spline(times, 1))
>>> angular_acceleration = np.rad2deg(spline(times, 2))
>>> times_plot = np.linspace(times[0], times[-1], 100)
>>> angles_plot = spline(times_plot).as_euler('XYZ', degrees=True)
>>> angular_rate_plot = np.rad2deg(spline(times_plot, 1))
>>> angular_acceleration_plot = np.rad2deg(spline(times_plot, 2))

在此图中，您可以看到欧拉角是连续且平滑的

>>> import matplotlib.pyplot as plt
>>> plt.plot(times_plot, angles_plot)
>>> plt.plot(times, angles, 'x')
>>> plt.title("Euler angles")
>>> plt.show()

角速度也很平滑

>>> plt.plot(times_plot, angular_rate_plot)
>>> plt.plot(times, angular_rate, 'x')
>>> plt.title("Angular rate")
>>> plt.show()

角加速度是连续的，但不是平滑的。另外请注意，角加速度不是分段线性函数，因为它不同于旋转向量的二阶导数（旋转向量的二阶导数是分段线性函数，就像立方样条中一样）。

>>> plt.plot(times_plot, angular_acceleration_plot)
>>> plt.plot(times, angular_acceleration, 'x')
>>> plt.title("Angular acceleration")
>>> plt.show()

返回在新标签页中打开

方法

__call__(times[, order])

计算插值值。