How do you get the azimuth and elevation from one enu vector to another enu vector?
A link to a formula or piece of code would be helpful. I'm not getting much information at all when searching.
Asked
Active
Viewed 9,519 times
1 Answers
1
You can calculate the azimuth and elevation angles between East-North-Up vectors (x,y,z) and (u,v,w) using the following:
- Subtract the vectors:
(x,y,z) - (u,v,w) = (x-u,y-v,z-w) = (x',y',z') - Compute the azimuth angle:
a = arctan(x'/y') = arctan((x-u)/(y-v)) - Compute the elevation angle:
e = arctan(z'/y') = arctan((z-w)/(y-v))
In Python:
v1 = np.array([3,4,4])
v2 = np.array([1,2,6])
v = v1 - v2
a = np.degrees(np.arctan(v[0]/v[1]))
e = np.degrees(np.arctan(v[2]/v[1]))
print('azimuth = '+str(a)+', elevation = '+str(e))
Output:
azimuth = 45.0, elevation = -45.0
Nathaniel
- 3,230
- 11
- 18
-
2What if `y'` is zero? – Jann Poppinga Nov 23 '20 at 13:03
-
1A bit late perhaps, but use np.arctan2(x, y) instead, it handles y = 0 just fine – alkanen Jan 30 '23 at 08:47
-
1Should z not be `arctan2(z', sqrt((y')^2 + (x')^2))`? Even if you look at the way this picture is drawn seems correct – gkv Feb 27 '23 at 00:33
-
@gkv I think you are correct. In the example the horizontal distance is `sqrt(2^2 + 2^2)`. The elevation angle should be `arctan2(2, sqrt(2^2 + 2^2))) = +/- 35.26` – Michiel van Lohuijzen Aug 14 '23 at 14:29