Direct conversion of latitude and height from one ellipsoid to another

Abstract

I suggest a method for converting geodetic height and latitude from one oblate ellipsoid of revolution to another having the same center and symmetry axis. Unlike other approaches, the method here does not obtain height and latitude from Earth-centered, Earth-fixed (ECEF) Cartesian coordinates; this feature allows height conversion with high accuracy even in cases where the data format limits the precision of the latitude data. Height and latitude conversions can be expressed as a Fourier series in even multiples of latitude, with height changes having only cosines and latitude changes having only sines. The absolute difference of the flattenings of the two ellipsoids furnishes a simple upper bound on the maximum absolute latitude change. Conversions between the TOPEX and WGS84 ellipsoids, a practical necessity for the inter-mission comparison of satellite laser and radar altimeter data, illustrate the findings. Because the differences in the flattenings and semi-major axes of these ellipsoids are small, truncating the Fourier series after the term in twice the latitude gives an approximate conversion with an error less than 9 × 10–12 radians of latitude and about 6 × 10–6 m of height.