Abstract Local geoid models presenting higher resolution than global ones are generally derived by a combination of different datasets, integrating individual pure astrogeodetic, gravimetric and GNSS/levelling solutions. To define local geoid, interpolators may be applied starting from dataset height values. It is well known that the accuracy resulting depends not only interpolation method, but...