SimpleTk.Euler3DTransform和scipy.spatial.transform.Rotation.from_euler之间的区别？

Question

使用以下两个库函数：

• SimpleITK.Euler3DTransform
• scipy.spatial.transform.Rotation.from_euler

要从Euler角度创建简单旋转矩阵，请执行以下操作：

import numpy as np
import SimpleITK as sitk
from scipy.spatial.transform import Rotation
from math import pi

euler_angles = [pi / 10, pi / 18, pi / 36]

sitk_matrix = sitk.Euler3DTransform((0, 0, 0), *euler_angles).GetMatrix()
sitk_matrix = np.array(sitk_matrix).reshape((3,3))
print(np.array_str(sitk_matrix, precision=3, suppress_small=True))

order = 'XYZ' # Different results for any order in ['XYZ','XZY','YZX','YXZ','ZXY','ZYX','xyz','xzy','yzx','yxz','zxy','zyx']
scipy_matrix = Rotation.from_euler(order, euler_angles).as_matrix()
print(np.array_str(scipy_matrix, precision=3, suppress_small=True))

我得到两个不同的结果：

[[ 0.976 -0.083  0.2  ]
 [ 0.139  0.947 -0.288]
 [-0.165  0.309  0.937]]
[[ 0.981 -0.086  0.174]
 [ 0.136  0.943 -0.304]
 [-0.138  0.322  0.937]]

为什么?？如何使用scipy计算与SimpleTk相同的矩阵

Answer 1

什么是您的“订单”字符串。当我使用order='xyz'运行代码时，对于SimpleTk和scipy的旋转，我得到了相同的结果