I think it is a little more straightforward if you work in Cartesian coordinates.
Take an earth of radius one, for simplicity.
Then we can take the longitudes as pml without loss of generality. Then
the halfway point between the two points, projected to the surface of the earth, will have the highest latitude (equivalently the highest z component).
The two points are (cosphicosl,cosphisinl,sinphi) and
(cosphicosl,−cosphisinl,sinphi), and the mid point is
(cosphicosl,0,sinphi). To project to the surface, we divide by
the norm, to get a z component (on our unit earth) of
sindelta=sinphioversqrt(cosphicosl)2+sin2phi=tanphioversqrtcos2l+tan2phi.
To obtain the arctan, we first need cosdelta, which we get from
cosdelta=sqrt1−sin2delta=cosloversqrtcos2l+tan2phi, and hence
tandelta=tanphiovercosl, which is the desired result.
No comments:
Post a Comment