pip config set global.index-url https://pypi.tuna.tsinghua.edu.cn/simple

pip install scipy

# coding=UTF-8
from scipy.spatial.transform import Rotation as R
 
# 定義歐拉角 XYZ順序
euler_angles = [90, 3, 3]

if (90 - abs(euler_angles[0])) < 0.5:
    euler_angles[0] = 89.5
euler_angles_zxy = [0, 0, 0]
euler_angles_zxy[0] = euler_angles[2]
euler_angles_zxy[1] = euler_angles[0]
euler_angles_zxy[2] = euler_angles[1]
euler_angles_zxy = [-x for x in euler_angles_zxy]
 
# 創(chuàng)建Rotation對象
r = R.from_euler('zxy', euler_angles_zxy, degrees=True)
 
# 獲取旋轉矩陣
rotation_matrix = r.as_dcm()
print('[')
for row in rotation_matrix:
    msg = ", ".join(map(str, row))
    msg = '[' + msg + '],'
    print(msg)
print(']')

# coding=UTF-8
import numpy as np
from scipy.spatial.transform import Rotation as R
 
# 假設我們有一個旋轉矩陣
#80 -7 45
rotation_matrix = np.array([
[0.9945219996629926, 0.10452647320585942, -0.00045671157999537665],
[-0.00045671157999537665, 0.008714576084760173, 0.9999619230641713],
[0.10452647320585942, -0.9944839227271639, 0.008714576084760173]
])
# 使用SciPy的Rotation對象轉換歐拉角
rotation = R.from_dcm(rotation_matrix)
euler_angles_zxy = rotation.as_euler('zxy', degrees=True)  # sco存儲使用'zxy'順序，如果需要度為單位，設置degrees=True

#左右手參考坐標系正負值替換
euler_angles_zxy = [-x for x in euler_angles_zxy]
euler_angles = [0, 0, 0]
euler_angles[0] = euler_angles_zxy[1]
euler_angles[1] = euler_angles_zxy[2]
euler_angles[2] = euler_angles_zxy[0]
print(euler_angles)