In this thesis, we introduce new interpolation methods for two dimensional data via constructing harmonic functions passing through the given data. Two ways to construct the harmonic function are introduced: (1) constructing the harmonic function via the heat equation, and (2) constructing the harmonic function via the boundary element method