import numpy
in_data = numpy.array([[25.89125824, 25.88840675], [25.90930748, 
    25.90640068]], dtype=numpy.float32)
print(map_coordinates(in_data.T, [[0.0], [0.125]], order=1, mode='nearest'))
print(interp(in_data, numpy.array([0, 1]), numpy.array([0, 1]), numpy.array
    ([0.0]), numpy.array([0.125]), order=1))
Scipy map_coordinates bilinear interpolation compared to interp and IDL interpolate