如何从Axis Studio导出的BVH文件获取头部/面部方向与旋转？

一、Axis Studio面部旋转的计算逻辑 #

Axis Studio是基于骨骼层级结构计算面部/头部旋转的：

• 头部骨骼的旋转是相对于父骨骼（颈部）的局部变换，XYZ欧拉顺序下，旋转按X→Y→Z的顺序依次叠加——先绕局部X轴旋转，再绕旋转后的局部Y轴，最后绕再次变换后的局部Z轴。
• 面部方向本质是头部骨骼的全局旋转方向，需要把局部欧拉角转换为全局旋转矩阵或四元数后，才能准确反映真实朝向；单纯的局部Y轴欧拉角受层级变换和欧拉顺序的万向锁影响，无法直接对应真实的面部左右转向。

二、XYZ欧拉顺序下Y轴旋转失效的核心原因 #

1. 万向锁限制：在XYZ顺序中，当X轴旋转接近±90°时，Y轴和Z轴会重合，导致Y轴旋转失去自由度，无法正确捕捉面部的左右转向。
2. 坐标系混淆：你直接使用的颈部到头部的局部Y轴欧拉角，是相对于颈部骨骼的局部坐标系，而非世界坐标系，自然无法直接对应你需要的面部全局朝向。

三、数据类型的选择建议 #

问题确实和导出的数据类型有关，推荐优先选择以下两种方案：

1. 四元数（Quaternion）：
   • 四元数不存在万向锁问题，能精准表示任意旋转，且插值更平滑。如果Axis Studio支持导出带四元数的骨骼数据（比如FBX、JSON格式），直接用四元数计算面部方向会更可靠。
   • 若只能导出BVH，可通过Python将欧拉角转换为四元数，再计算全局旋转来提取面部朝向，示例代码：

     from scipy.spatial.transform import Rotation as R
     # 读取BVH的局部XYZ欧拉角（单位：度）
     local_euler = [x_rot, y_rot, z_rot]
     # 转换为局部旋转四元数
     local_rot = R.from_euler('xyz', local_euler, degrees=True)
     # 结合颈部骨骼的全局旋转，计算头部的全局旋转
     parent_global_rot = R.from_quat(neck_global_quat)
     head_global_rot = parent_global_rot * local_rot
     # 提取面部前向向量（假设骨骼定义中Z轴为前向）
     forward_dir = head_global_rot.apply([0, 0, 1])
     

2. 旋转矩阵（Rotation Matrix）：
   • 旋转矩阵也能无歧义表示旋转，便于进行坐标系变换。可通过欧拉角转换得到全局旋转矩阵后，直接提取面部朝向向量。

如果必须保留BVH格式和XYZ欧拉顺序，建议不要直接使用局部欧拉角，而是将其转换为全局四元数或旋转矩阵后再提取面部方向——这样既不影响你其他基于XYZ顺序的计算，又能解决面部朝向的问题。

内容的提问来源于stack exchange，提问作者Claudio