iOS OpenGL ES 2.0四元数旋转Slerp到XYZ位置

时间:2013-06-05 00:32:24

标签: iphone ios opengl-es opengl-es-2.0 glkit

我正在关注四元数教程:http://www.raywenderlich.com/12667/how-to-rotate-a-3d-object-using-touches-with-opengl并尝试将地球仪旋转到某个XYZ位置。我有一个初始四元数,并在地球表面生成一个随机的XYZ位置。我将XYZ位置传递给以下函数。我们的想法是使用GLKMatrix4MakeLookAt生成一个lookAt向量,并从lookAt矩阵中为slerp步骤定义结束四元数。

- (void)rotateToLocationX:(float)x andY:(float)y andZ:(float)z {

    // Turn on the interpolation for smooth rotation
    _slerping = YES; // Begin auto rotating to this location
    _slerpCur = 0;
    _slerpMax = 1.0;
    _slerpStart = _quat;

    // The eye location is defined by the look at location multiplied by this modifier
    float modifier = 1.0;

    // Create a look at vector for which we will create a GLK4Matrix from
    float xEye = x;
    float yEye = y;
    float zEye = z;
    //NSLog(@"%f %f %f %f %f %f",xEye, yEye, zEye, x, y, z);
    _currentSatelliteLocation = GLKMatrix4MakeLookAt(xEye, yEye, zEye, 0, 0, 0, 0, 1, 0);
    _currentSatelliteLocation = GLKMatrix4Multiply(_currentSatelliteLocation,self.effect.transform.modelviewMatrix);

    // Turn our 4x4 matrix into a quat and use it to mark the end point of our interpolation
    //_currentSatelliteLocation = GLKMatrix4Translate(_currentSatelliteLocation, 0.0f, 0.0f, GLOBAL_EARTH_Z_LOCATION);
    _slerpEnd = GLKQuaternionMakeWithMatrix4(_currentSatelliteLocation);

    // Print info on the quat
    GLKVector3 vec = GLKQuaternionAxis(_slerpEnd);
    float angle = GLKQuaternionAngle(_slerpEnd);
    //NSLog(@"%f %f %f %f",vec.x,vec.y,vec.z,angle);

    NSLog(@"Quat end:");
    [self printMatrix:_currentSatelliteLocation];
    //[self printMatrix:self.effect.transform.modelviewMatrix];

}

插值工作,我得到一个平滑的旋转,但结束位置永远不是我输入的XYZ - 我知道这是因为我的地球是一个球体,我正在从Lat Lon计算XYZ。我希望在旋转之后,从地球表面上的纬度/经度位置向地球中心直接向下看“lookAt”向量。我认为它可能与up向量有关但我已经尝试了所有有意义的东西。

我做错了什么 - 如何定义最终的四元数,当我完成旋转时,向下看向地球表面上的XYZ向量?谢谢!

2 个答案:

答案 0 :(得分:1)

以下是您的意思: 你的地球中心是(0,0,0),半径是R,起始位置是(0,0,R),你的最终位置是(0,R,0),所以围绕X-asix旋转地球90度? 如果是这样,只需将lookat功能眼睛位置设置到您的最终位置,查看地球中心的参数。

m_target.x = 0.0f;
m_target.y = 0.0f;
m_target.z = 1.0f;

m_right.x = 1.0f;
m_right.y = 0.0f;
m_right.z = 0.0f;

m_up.x = 0.0f;
m_up.y = 1.0f;
m_up.z = 0.0f;
void CCamera::RotateX( float amount )
{
    Point3D target = m_target;
    Point3D up = m_up;

    amount = amount / 180 * PI;

    m_target.x = (cos(PI / 2 - amount) * up.x) + (cos(amount) * target.x);
    m_target.y = (cos(PI / 2 - amount) * up.y) + (cos(amount) * target.y);
    m_target.z = (cos(PI / 2 - amount) * up.z) + (cos(amount) * target.z);

    m_up.x = (cos(amount) * up.x) + (cos(PI / 2 + amount) * target.x);
    m_up.y = (cos(amount) * up.y) + (cos(PI / 2 + amount) * target.y);
    m_up.z = (cos(amount) * up.z) + (cos(PI / 2 + amount) * target.z);

    Normalize(m_target);
    Normalize(m_up);
}

void CCamera::RotateY( float amount )
{
    Point3D target = m_target;
    Point3D right = m_right;

    amount = amount / 180 * PI;

    m_target.x = (cos(PI / 2 + amount) * right.x) + (cos(amount) * target.x);
    m_target.y = (cos(PI / 2 + amount) * right.y) + (cos(amount) * target.y);
    m_target.z = (cos(PI / 2 + amount) * right.z) + (cos(amount) * target.z);

    m_right.x  = (cos(amount) * right.x) + (cos(PI / 2 - amount) * target.x);
    m_right.y  = (cos(amount) * right.y) + (cos(PI / 2 - amount) * target.y);
    m_right.z  = (cos(amount) * right.z) + (cos(PI / 2 - amount) * target.z);

    Normalize(m_target);
    Normalize(m_right);
}

void CCamera::RotateZ( float amount )
{
    Point3D right = m_right;
    Point3D up = m_up;

    amount = amount / 180 * PI;

    m_up.x = (cos(amount) * up.x) + (cos(PI / 2 - amount) * right.x);
    m_up.y = (cos(amount) * up.y) + (cos(PI / 2 - amount) * right.y);
    m_up.z = (cos(amount) * up.z) + (cos(PI / 2 - amount) * right.z);

    m_right.x = (cos(PI / 2 + amount) * up.x) + (cos(amount) * right.x);
    m_right.y = (cos(PI / 2 + amount) * up.y) + (cos(amount) * right.y);
    m_right.z = (cos(PI / 2 + amount) * up.z) + (cos(amount) * right.z);

    Normalize(m_right);
    Normalize(m_up);
}

void CCamera::Normalize( Point3D &p )
{
    float length = sqrt(p.x * p.x + p.y * p.y + p.z * p.z);
    if (1 == length || 0 == length)
    {
        return;
    }

    float scaleFactor = 1.0 / length;
    p.x *= scaleFactor;
    p.y *= scaleFactor;
    p.z *= scaleFactor;
}

答案 1 :(得分:0)

这个问题的答案是以下rotateTo函数的组合以及对Ray(http://www.raywenderlich.com/12667/how-to-rotate-a-3d-object-using-touches-with-opengl)教程中代码的更改。正如该文章的评论之一所说,在GLKQuaternion中有一个任意因子2.0乘以Q_rot = GLKQuaternionMakeWithAngleAndVector3Axis(angle * 2.0,axis);.删除“2”并使用以下函数创建_slerpEnd - 之后地球将平滑旋转到指定的XYZ。

// Rotate the globe using Slerp interpolation to an XYZ coordinate
- (void)rotateToLocationX:(float)x andY:(float)y andZ:(float)z {

    // Turn on the interpolation for smooth rotation
    _slerping = YES; // Begin auto rotating to this location
    _slerpCur = 0;
    _slerpMax = 1.0;
    _slerpStart = _quat;

    // Create a look at vector for which we will create a GLK4Matrix from
    float xEye = x;
    float yEye = y;
    float zEye = z;
    _currentSatelliteLocation = GLKMatrix4MakeLookAt(xEye, yEye, zEye, 0, 0, 0, 0, 1, 0);

    // Turn our 4x4 matrix into a quat and use it to mark the end point of our interpolation
    _slerpEnd = GLKQuaternionMakeWithMatrix4(_currentSatelliteLocation);

}