In Geographic Information Science, polynomial methods such as linear estimation and non-polynomial methods including Inverse Distance Weighting and Kriging have been used for elevation data interpolation. In this paper, 3D data interpolation using linear and non-linear homotopy continuation as well as advanced polynomial interpolation methods are researched. Continuous deformations that reconstruct straight lines or algebraic curves between any pair of 3D data are presented. The implemented topological mathematical algorithm for 3D elevation data interpolation is compared to Inverse Distance Weighting and Triangulated Irregular Network (TIN) methods. The presented linear and non-linear mathematical algorithms show better results compared to Inverse Distance Weighting and TIN in terms of Root Mean Square Error and L-infinity.