Constrain pitch, yaw, & roll

Question

I've a rotation represented as a quaternion and am trying to constrain the pitch, yaw, & roll axes. I tried doing so thusly:

public struct Orientation
{
    public Vector3 up, forward;

    public Orientation(Vector3 up, Vector3 forward)
    {
        this.up = up;
        this.forward = forward;
    }
}

public static Orientation[] orientations = new Orientation[3]
{
    new Orientation(Vector3.right, Vector3.up),
    new Orientation(Vector3.up, Vector3.forward),
    new Orientation(Vector3.forward, Vector3.right)
};

public enum Axis
{
    Pitch,
    Yaw,
    Roll
};

private Vector3 ConstrainAxis(Vector3 vector, Axis axis, float from, float to)
{
    Orientation orientation = orientations[(int)axis];

    float theta = (to - from) * 0.5F;

    Vector3 cons = Quaternion.AngleAxis(from + theta, orientation.up) * orientation.forward;
    Vector3 proj = Vector3.ProjectOnPlane(vector, orientation.up);

    return ConstrainVector(cons.normalized, proj.normalized, theta);
}

private Vector3 ConstrainVector(Vector3 from, Vector3 to, float angle)
{
    float theta = Mathf.Abs(angle / Vector3.Angle(from, to));

    if(theta < 1.0F)
    {
        return Vector3.Slerp(from, to, theta);
    }

    return to;
}

Which turned out to be nothing more than an over-complicated way of constraining the individual components of an euler angle representation, of which both are subject to a strange jittering issue (gimbal lock related?).

What is the best approach to constraining these axes?

Constraining Euler Angels sounds like incorrect task. It is not fit to psysical motion , and probably is used to constrain only 2 axes. Is it possible that your task represent constraints of swing twist constraint of joint? — minorlogic
Indeed it does. I have an IK system where I need to constrain the joints. — Tannz0rz
Than you can decompose rotation into twist swing ? and apply constraints like in this tip alinenormoyle.com/weblog/?p=726 — minorlogic
That's very nice, but I can't say I understand how I would go about applying a constraint to a quaternion in this fashion. Is it as simple as constraining the real part w? — Tannz0rz
you can constraint quat with angle axis or directly constrain the magnitude of quat "vector" part. (don't forget , it is sin(0.5 * angle), and recalculate W with sqrt(1- vector_part.magnitude())) — minorlogic

minorlogic minorlogic · Accepted Answer · 2015-09-29T14:43:44

For joint constraints it is common practice to use "swing twist" parametrization. To represent current rotation as "swing twist" for quaternions, theare are good decomposition http://www.alinenormoyle.com/weblog/?p=726

And constraint for "swing" and "twist" can be done with quaternions.

if we want to constrain swing to +-30 degrees , pseudocode looks like

Quaternion swing;
const double maxMagnitude = sin(0.5 * toRad(30));
const double maxMagnitudeW = sqrt(1.0 - maxMagnitude * maxMagnitude);
if (swing.vec().normSqr() > maxMagnitude * maxMagnitude)
{
    swing.vec() = swing.vec().normalized() * maxMagnitude;
    swing.w() = maxMagnitudeW;
}

Constrain pitch, yaw, & roll

2 Answers