Local geoid height approximation and interpolation using moving least squares approach

In this paper an introduction of the moving least squares approach is presented in the context of data approximation and interpolation problems in Geodesy. An application of this method is presented for geoid height approximation and interpolation using different polynomial basis functions for the approximant and interpolant, respectively, in a regular grid of geoid height data in the region 16.0417 degrees <= phi <= 47.9583 degrees and 36.0417 degrees <= lambda <= 69.9582 degrees, with increment 0.0833 degrees in both latitudal and longitudal directions. The results of approximation and interpolation are then compared with the geoid height data from GPS-Levelling approach. Using the standard deviation of the difference of the results, it is shown that the planar interpolant, with reciprocal of distance as weight function, is the best choice in this local approximation and interpolation problem. (C) 2020 Institute of Seismology, China Earthquake Administration, etc. Production and hosting by Elsevier B.V. on behalf of KeAi Communications Co., Ltd.