34

How would you find the signed angle theta from vector a to b?

And yes, I know that theta = arccos((a.b)/(|a||b|)).

However, this does not contain a sign (i.e. it doesn't distinguish between a clockwise or counterclockwise rotation).

I need something that can tell me the minimum angle to rotate from a to b. A positive sign indicates a rotation from +x-axis towards +y-axis. Conversely, a negative sign indicates a rotation from +x-axis towards -y-axis.

assert angle((1,0),(0,1)) == pi/2.
assert angle((0,1),(1,0)) == -pi/2.
Cerin
  • 50,711
  • 81
  • 269
  • 459

2 Answers2

64

What you want to use is often called the “perp dot product”, that is, find the vector perpendicular to one of the vectors, and then find the dot product with the other vector.

if(a.x*b.y - a.y*b.x < 0)
    angle = -angle;

You can also do this:

angle = atan2( a.x*b.y - a.y*b.x, a.x*b.x + a.y*b.y );
Derek Ledbetter
  • 4,211
  • 3
  • 17
  • 16
37

If you have an atan2() function in your math library of choice:

signed_angle = atan2(b.y,b.x) - atan2(a.y,a.x)
Sparr
  • 7,297
  • 27
  • 46
  • 9
    What about a = (-1,1) and b = (-1,-1), where the answer should be pi/2? You should check if the absolute value is bigger than pi, and then add or subtract 2*pi if it is. – Derek Ledbetter Jan 27 '10 at 21:52
  • @Derek Good catch. I actually discovered this myself while implementing the solution. – Cerin Jan 28 '10 at 13:06
  • 1
    it is not very appropriate for example for computer graphics because it confuse -pi and pi if I have a = {-1, 0} and b = {0, 1}. – user2083364 Oct 17 '13 at 08:26
  • In degree, the result of this answer should be [-180, 180), but some time I discover result like: 358.5. Derek Ledbetter's answer works fine. – Helin Wang Jan 10 '14 at 02:47