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 ...