Transformation from Cartesian to Geodetic Coordinates Accelerated by Halley’s Method

By using Halley’s third-order formula to find the root of a non-linear equation, we develop a new iterative procedure to solve an irrational form of the “latitude equation”, the equation to determine the geodetic latitude for given Cartesian coordinates. With a limit to one iteration, starting from zero height, and minimizing the number of divisions by means of the rational form representation of Halley’s formula, we obtain a new non-iterative method to transform Cartesian coordinates to geodetic ones. The new method is sufficiently precise in the sense that the maximum error of the latitude and the relative height is less than 6 micro-arcseconds for the range of height, −10 km ≤ h ≤ 30,000 km. The new method is around 50% faster than our previous method, roughly twice as fast as the well-known Bowring’s method, and much faster than the recently developed methods of Borkowski, Laskowski, Lin and Wang, Jones, Pollard, and Vermeille.