We describe a novel interpolation algorithm to find the optimal image intensity function generating an optimal gray-level estimation of interpolated pixels of digital images. The new approach is based on the proposed image block mapping method and least-square support vector machines (LSSVM) with Gaussian radial basis function (RBF) kernels. With the mapping technique, the interpolation procedure of the LSSVM is actually accomplished in the same input vector space. A number of different scale interpolation experiments are carried out. The experimental results demonstrate that the performance of the proposed algorithm is competitive with many other existing methods, such as cubic, spline, and linear methods. The peak signal-to-noise ratio of the image reconstructed by the proposed algorithm is higher than those obtained by the spline. And the estimated accuracy of the proposed algorithm is similar to that of the cubic algorithm, while the computational requirement is lower than the latter.