astropy.coordinates
Position Angle (East of North) between two points on a sphere.
Angle
Quantity
float
Longitude and latitude of the two points. Quantities should be in angular units; floats in radians.
The (positive) position angle of the vector pointing from position 1 to position 2. If any of the angles are arrays, this will contain an array following the appropriate numpy broadcasting rules.
numpy