The aim of the paper is the comparison of the least squares prediction presented by Heiskanen and Moritz (1967) in the classical handbook “Physical Geodesy” with the geostatistical method of simple kriging as well as in case of Gaussian random fields their equivalence to conditional expectation. The paper contains also short notes on the extension of simple kriging to ordinary kriging by dropping the assumption of known mean value of a random field as well as some necessary information on random fields, covariance function and semivariogram function. The semivariogram is emphasized in the paper, for two reasons. Firstly, the semivariogram describes broader class of phenomena, and for the second order stationary processes it is equivalent to the covariance function. Secondly, the analysis of different kinds of phenomena in terms of covariance is more common. Thus, it is worth introducing another function describing spatial continuity and variability. For the ease of presentation all the considerations were limited to the Euclidean space (thus, for limited areas) although with some extra effort they can be extended to manifolds like sphere, ellipsoid, etc.

A new method to transform from Cartesian to geodetic coordinates is presented. It is based on the solution of a system of nonlinear equations with respect to the coordinates of the point projected onto the ellipsoid along the normal. Newton’s method and a modification of Newton’s method were applied to give third-order convergence. The method developed was compared to some well known iterative techniques. All methods were tested on three ellipsoidal height ranges: namely, (-10 – 10 km) (terrestrial), (20 – 1000 km), and (1000 – 36000 km) (satellite). One iteration of the presented method, implemented with the third-order convergence modified Newton’s method, is necessary to obtain a satisfactory level of accuracy for the geodetic latitude ( σ φ < 0.0004”) and height ( σ h < 10 − 6 km, i.e. less than a millimetre) for all the heights tested. The method is slightly slower than the method of Fukushima (2006) and Fukushima’s (1999) fast implementation of Bowring’s (1976) method.