python - 无法从欧拉角创建 quat,偏航角超过 90 glm
问题描述
这里有一些背景:我正在开发一个游戏引擎,我最近添加了一个物理引擎(在本例中为 physx)。问题是我的变换类使用欧拉角进行旋转,而物理引擎的变换类使用欧拉角。所以我刚刚实现了一个方法来将我的变换类更改为物理引擎变换并返回。它运行良好,但我发现了一个奇怪的错误。
我得到的行为:
当旋转的偏航(欧拉矢量的第二个元素)高于 90 度时,它不会使对象在 y 轴上再旋转并开始与俯仰和滚动搞乱(奇怪的摇晃从 0 跳到 180 和回来很多)。调试工具显示旋转不会超过 91,但会达到最大 90.0003 我确实将度数转换为弧度。 示例: 为了显示这个错误,我有一个带有 python 脚本旋转它的立方体:
from TOEngine import *
class rotate:
direction = vec3(0,10,0)
def Start(self):
pass
def Update(self,deltaTime):
transform.Rotate(self.direction*deltaTime*5)
pass
引擎本身是用 cpp 编写的,但我有一个使用嵌入式 python 的脚本系统。TOEngine 只是我的模块,脚本本身只是每帧旋转立方体。立方体从 0 , 0 , 0 旋转开始,旋转很好,但停止,90 度偏航并开始摇晃。
这仅在启用物理系统时发生,因此我知道错误必须在使用 glm 将旋转从 euler 传输到 quat 并返回每一帧的方法中。
这是实际有问题的代码:
void RigidBody::SetTransform(Transform transform)
{
glm::vec3 axis = transform.rotation;
rigidbody->setGlobalPose(PxTransform(*(PxVec3*)&transform.position,*(PxQuat*)&glm::quat(glm::radians(transform.rotation))));//Attention Over Here
}
Transform RigidBody::GetTransform()
{
auto t = rigidbody->getGlobalPose();
return Transform(*(glm::vec3*)&t.p, glm::degrees(glm::eulerAngles(*(glm::quat*)&t.q)), entity->transform.scale);
}
避免奇怪的类型双关语 PxQuat 与 glm::quat 基本相同,PxVec3 与 glm::vec3 基本相同。我希望此代码通过将旋转从欧拉角角度更改为带有弧度的 quat(困难部分),在物理引擎变换类和我的变换类之间传输。
以及物理系统内部:
void PreUpdate(float deltaTime)override { //Set Physics simulation changes to the scene
mScene->fetchResults(true);
for (auto entity : Events::scene->entities)
for (auto component : entity->components)
if (component->GetName() == "RigidBody")
entity->transform = ((RigidBody*)component)->GetTransform(); //This is running on the cube entity
}
void PostUpdate(float deltaTime)override { //Set Scene changes To Physics simulation
for (auto entity : Events::scene->entities)
for (auto component : entity->components)
if (component->GetName() == "RigidBody")
((RigidBody*)component)->SetTransform(entity->transform);//This is running on the cube entity
mScene->simulate(deltaTime);
}
PreUpdate 在每帧更新之前运行 PostUpdate 在每帧更新之后运行。顾名思义,Update 方法(在上面的脚本中显示)在更新上运行...(在 PreUpdate 和 PostUpdate 之间)。立方体有一个刚体组件。我期望得到:一个旋转的立方体,当它达到 90 度偏航时不会停止旋转。
我知道这个有点复杂。我尽力解释了我认为问题在于将欧拉角更改为四元组的错误。
解决方案
关于从PxQuat
to的转换glm::quat
,请阅读https://en.cppreference.com/w/cpp/language/explicit_cast和https://en.cppreference.com/w/cpp/language/reinterpret_cast上的文档并查找 undefined页面上的行为reinterpret_cast
。据我所知,c 风格的演员表不能保证有效,甚至是不可取的。虽然我在这一点上离题了,但请记住,您有两个选项可用于此转换。
glm::quat glmQuat = GenerateQuat();
physx::PxQuat someQuat = *(physx::PxQuat*)(&glmQuat); //< (1)
physx::PxQuat someOtherQuat = ConvertGlmQuatToPxQuat(glmQuat); //< (2)
(1) 此选项有可能导致未定义的行为,但更重要的是,您没有保存副本。该语句肯定会导致 1 个复制构造函数调用。
(2) 由于返回值优化,此选项也将导致physx::PxQuat
.
因此,实际上,通过采用选项 (1),您不会节省任何成本,而是冒着不确定行为的风险。使用选项 (2),成本相同,但代码现在符合标准。现在回到原点。
我通常会尽我所能避免使用欧拉角,因为它们容易出错并且比四元数更令人困惑。也就是说,这是一个简单的测试,您可以设置它来测试从欧拉角转换的四元数(暂时不要使用 physx)。
您需要生成以下方法。
glm::mat3 CreateRotationMatrix(glm::vec3 rotationDegrees);
glm::mat3 CreateRotationMatrix(glm::quat inputQuat);
glm::quat ConvertEulerAnglesToQuat(glm::vec3 rotationDegrees);
然后您的测试伪代码如下所示。
for (auto angles : allPossibleAngleCombinations) {
auto expectedRotationMatrix = CreateRotationMatrix(angles);
auto convertedQuat = ConvertEulerAnglesToQuat(angles);
auto actualRotationMatrix = CreateRotationMatrix(convertedQuat);
ASSERT(expectedRotationMatrix, actualRotationMatrix);
}
只有当这个测试为你通过时,你才能看下一个将它们转换为PxQuat
. 我猜这个测试对你来说会失败。我要给出的一个建议是输入角度之一(取决于约定)需要限制范围。比如说,如果您将偏航角限制在 -90、90 度之间,那么您的测试可能会成功。这是因为,欧拉角的非唯一组合可能导致相同的旋转矩阵。
推荐阅读
- excel - 修复标签顺序问题
- ios - React Native 存档错误:无法找到 React Native 文件。确保在您的项目依赖项中安装了“react-native”模块
- arduino - 有没有办法为一个输出找到所有可能的输入?
- kotlin - 如何修复错误“java.lang.IllegalStateException:RecyclerView 不能为空”
- django - 在 Django 中命名基于类的视图有错误的方法吗?
- c# - SQL 中的多个 IN 子句
- javascript - 我的车把模板在渲染时自动终止
- javascript - 滚动页面滚动点击非常高的页面 - 无法正常工作
- python - Tensorflow 客户化培训 - ValueError:变量密集/内核/亚当/不存在?
- visual-studio-code - 使用 rust (rls) 扩展在 VS Code 中构建 rust 的最佳实践是什么?