1

我正在尝试使用 Irrlicht 作为图形引擎和物理 ODE 在 C++ 中编写 3d 模拟。然后我使用一个函数将 ODE 四元数转换为 Irrlicht Euler 角。为了做到这一点，我正在使用这段代码。

void QuaternionToEuler(const dQuaternion quaternion, vector3df &euler)
{
    dReal w,x,y,z;

    w = quaternion[0];
    x = quaternion[1];
    y = quaternion[2];
    z = quaternion[3];

    double sqw = w*w;    
    double sqx = x*x;    
    double sqy = y*y;    
    double sqz = z*z; 

    euler.Z = (irr::f32) (atan2(2.0 * (x*y + z*w),(sqx - sqy - sqz + sqw)) * (180.0f/irr::core::PI));
    euler.X = (irr::f32) (atan2(2.0 * (y*z + x*w),(-sqx - sqy + sqz + sqw)) * (180.0f/irr::core::PI));          
    euler.Y = (irr::f32) (asin(-2.0 * (x*z - y*w)) * (180.0f/irr::core::PI));

}

它可以很好地绘制正确的位置和旋转，但问题与asin指令有关。它只返回0..90-范围内的值0..-90，我需要从0..360度数中获取一个范围。至少我需要0..360在我打电话时的范围内进行轮换node->getRotation().Y。

4

2 回答 2

6

欧拉角（任何类型）具有奇点。在您使用的那些特定欧拉角的情况下（看起来像 Tait-Bryan 角或其一些变体），奇点位于正负 90 度的俯仰角 (Y)。这是欧拉角的固有限制，也是它们很少在任何严肃环境中使用的主要原因之一（飞机动力学除外，因为所有飞机的速度矢量（可能不是水平的）俯仰能力非常有限，所以它们很少靠近那个奇点）。

这也意味着您的计算实际上只是两个等效解决方案之一。对于给定的四元数，欧拉角有两种解法，代表相同的旋转，一种在奇点的一侧，另一种反映第一个。由于两种解决方案是等效的，您只需选择最简单的一侧，即俯仰角在 -90 到 90 度之间的那个。

此外，您的代码需要处理接近奇点以避免获得 NaN。换句话说，您必须检查您是否接近（小公差）奇异点（俯仰角为 -90 和 90 度），如果是，请使用替代公式（只能计算一个最接近的角度）旋转）。

如果您有任何方法可以完全避免使用欧拉角，我强烈建议您这样做，几乎任何旋转表示都比欧拉角更可取。Irrlicht 本机使用矩阵，还支持通过轴角表示设置/获取旋转，这更容易使用（并且更容易从四元数获得，并且没有奇点）。

于 2013-01-21T21:41:58.203 回答
1

想想地球的地球。其上的每个点只能使用纬度（在范围内[-90, 90]）和经度（在范围内[-180, 180]）来定义。因此，可以使用这些角度来指定球体上的每个点。现在球体上的一个点指定一个向量，球体上的所有点指定所有可能的向量。因此，正如本文所指出的，您使用的公式将生成所有可能的方向。

希望这可以帮助。

于 2013-01-21T21:38:49.000 回答