The creation of a three-dimensional geoinformation model of aerodrome is considered. An approach is proposed, consisting in supplementing a two-dimensional vector map with the values of the third coordinate, which are calculated using data from a digital terrain model. To implement this approach, the problem of interpolation of a function of two variables with a large number of nodal points is solved. The number of nodes where elevation values are set in a digital terrain model can be several thousand for each dimension. In this case, the use of global interpolation methods, for example, splines, in which function values at all nodal points are used to determine the values of a function at a given point falling into a grid cell, is computationally inefficient, since it requires a lot of machine time. It seems more promising to use a local approach to the interpolation problem. In this approach, the value of the function at an arbitrary point is determined by the values of the function at the nodal points of the grid closest to the specified point. Finite formulas are obtained that do not require solving equations, and allow interpolating functions for cases of linear and smooth approximation. An algorithm and a calculation program have been developed for a regular two-dimensional grid of nodes of a digital terrain model using the publicly available SRTM database. The results of numerical calculations of test cases with a comparison of linear and smooth approximations are presented. Using the proposed approach, a three-dimensional model of the Sochi aerodrome was created to simulate the movement of aircraft on its surface. The proposed approach and the interpolation method can be used to calculate the elevation of given points and to create three-dimensional cartographic models of the earth's surface.