如何将MLKit输出的3D人体关键点转换为符合规范的BVH文件

直接给核心结论：
BVH文件里存储的每帧骨骼旋转值，既不是相邻两帧的姿态差值，也不是当前帧和首帧的姿态差值，参照基准是两个：

1. 父骨骼的局部坐标系
2. HIERARCHY段提前定义好的骨骼初始绑定姿态（REST姿态，通常是T pose）

你之前纠结的三个参照选项里，正确的是「当前骨骼当前姿态 vs HIERARCHY定义的初始骨骼姿态」，且所有计算都必须在对应父骨骼的局部空间内完成，没有帧间增量的概念——BVH每帧存储的是独立的绝对局部姿态，跳帧读取也能正确还原姿态，不需要依赖前一帧的数据。

你当前用的算法逻辑从根上不成立：

for all landmark position:
  u = vector representing the bone (landmark to next landmark in hierarchy)
  v = vector representing previous bone in hierarchy
  q = quaternion between u and v
  return euler from q

这段逻辑算的是相邻两根骨骼的夹角，既没有引入绑定姿态的基准向量，也没有做坐标空间转换，输出的结果根本不是BVH要求的局部旋转参数，效果不符合预期是必然的。

前提：你已经把MLKit输出的坐标系转换为BVH通用的右手坐标系（Y轴向上，和你HIERARCHY定义的坐标系轴向完全对齐），且已经完成了MLKit骨骼拓扑到通用人形骨骼拓扑的映射。

• 第一步：固定HIERARCHY段的静态内容
  • 按你用的通用人形拓扑写好骨骼父子层级
  • 给每个骨骼填写固定的OFFSET值：这个值就是绑定姿态下，当前骨骼原点相对于父骨骼原点的向量（比如T pose下，左大臂到左小臂的OFFSET就是沿大臂轴向指向手肘的固定长度向量），OFFSET是静态值，所有动画帧都不会改变。
  • 给每个节点定义好CHANNEL顺序，常规人形BVH的根节点（髋部）通道一般是Xposition Yposition Zposition Zrotation Xrotation Yrotation，其他子骨骼只有三个旋转通道，通道顺序一旦定了，后面欧拉角解算必须严格匹配这个顺序。
• 第二步：逐帧逐节点计算运动参数
  • 根节点处理：三个平移值直接填根节点（髋部）转换完坐标系后的世界坐标即可；根节点的旋转计算逻辑和其他子骨骼一致，参照空间是世界空间。
  • 非根节点旋转计算：
    1. 计算当前帧下该骨骼的世界空间朝向向量：bone_vec_world = 当前骨骼的子关节世界坐标 - 当前骨骼自身的世界坐标，做归一化
    2. 取出该骨骼的绑定姿态朝向向量：就是HIERARCHY里该骨骼的OFFSET值，归一化得到bone_vec_bind
    3. 把上述两个向量都转换到当前骨骼的父骨骼局部坐标系下：你需要提前从根节点开始，逐节点维护父节点到世界空间的变换矩阵，求逆之后就能把世界空间向量转到父节点局部空间
    4. 计算从bone_vec_bind旋转到当前帧bone_vec_world的四元数q_local。注意：仅靠骨骼两个端点的坐标只能算出两个轴向的旋转，骨骼绕自身长轴的扭转角（比如大臂绕自己轴线转手腕的角度）是无法通过两个端点计算的，MLKit输出的关键点本身不包含这个自由度的信息，普通场景直接把扭转角设为0即可，如果需要更高精度，可以通过手肘/膝盖的侧面朝向辅助点补算扭转值。
    5. 严格按照你之前在CHANNEL里定义的旋转顺序，把四元数q_local转换为对应顺序的欧拉角，单位为度。
• 第三步：后处理避免跳变
  逐通道对欧拉角做连续化处理：如果相邻两帧同个通道的角度差超过180度，就对当前帧角度加/减360度，把差值修正到±180度区间内，避免播放时出现骨骼转圈的跳变问题。

• 旋转顺序绝对不能错：BVH的旋转是按CHANNEL里写的顺序从右往左叠加变换的，比如通道写的是Zrotation Xrotation Yrotation，变换矩阵的乘法顺序就是R_y * R_x * R_z，转欧拉角的时候必须匹配这个顺序，否则姿态会完全扭曲。
• 不要在世界空间直接算旋转：BVH的变换是层级叠加的，子骨骼的旋转永远是相对于父骨骼局部空间的，跳过坐标空间转换直接在世界空间算出来的旋转，叠加上父节点的旋转之后会出现双倍偏移。
• 不要试图用帧间差值算旋转：BVH不存储增量，存的是每帧的绝对局部姿态，用帧间差值算出来的参数会出现误差累计，播放几帧之后姿态就会完全错位。

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