Angles between two n-dimensional vectors in Python

by

Note: all of the other answers here will fail if the two vectors have either the same direction (ex, (1, 0, 0), (1, 0, 0)) or opposite directions (ex, (-1, 0, 0), (1, 0, 0)).

Here is a function which will correctly handle these cases:

import numpy as np

def unit_vector(vector):
    """ Returns the unit vector of the vector.  """
    return vector / np.linalg.norm(vector)

def angle_between(v1, v2):
    """ Returns the angle in radians between vectors 'v1' and 'v2'::

            >>> angle_between((1, 0, 0), (0, 1, 0))
            1.5707963267948966
            >>> angle_between((1, 0, 0), (1, 0, 0))
            0.0
            >>> angle_between((1, 0, 0), (-1, 0, 0))
            3.141592653589793
    """
    v1_u = unit_vector(v1)
    v2_u = unit_vector(v2)
    return np.arccos(np.clip(np.dot(v1_u, v2_u), -1.0, 1.0))

More Related Contents:

  1. numpy matrix vector multiplication [duplicate]
  2. How to know the angle between two vectors?
  3. Subtracting 2 lists in Python
  4. Rotation of 3D vector?
  5. pygame 2 dimensional movement of an enemy towards the player, how to calculate x and y velocity?
  6. Overloading Addition, Subtraction, and Multiplication Operators
  7. Use vector2 in pygame. Collide with the window frame and restrict the ball to the rectangular area
  8. plot decision boundary matplotlib
  9. Distance between point and a line (from two points)
  10. How to create an empty R vector to add new items
  11. Concatenating column vectors using numpy arrays
  12. How to convert a string in the desired manner?
  13. How to serve static files in Flask
  14. Sibling package imports
  15. proper name for python * operator?
  16. Is it possible to break a long line to multiple lines in Python? [duplicate]
  17. Cosine Similarity between 2 Number Lists
  18. How to copy a class?
  19. Converting epoch time with milliseconds to datetime
  20. How to convert an int to a hex string?
  21. Splitting timestamp column into separate date and time columns
  22. How to fill matplotlib bars with a gradient?
  23. How does the python socket.recv() method know that the end of the message has been reached?
  24. TensorFlow saving into/loading a graph from a file
  25. Get coordinates of local maxima in 2D array above certain value
  26. How to get element-wise matrix multiplication (Hadamard product) in numpy?
  27. Python urlparse — extract domain name without subdomain
  28. Combine 3 separate numpy arrays to an RGB image in Python
  29. How to change all the dictionary keys in a for loop with d.items()?
  30. Detecting computer/program shutdown in Python?

Leave a Comment