An algorithm for interpolation by piece-wise polynomial approximation has been developed and implemented. The hypersurface approximation is carried out by a quadratic surface defined over a two dimensional space of the digital picture in the neighbourhood of the point to be interpolated using orthogonal polynomials as basis functions. The kernel for interpolation has been evaluated with the help of polynomial basis functions; this enables efficient implementation of the algorithm.