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'))
Билинейная интерполяция Scipy map_coordinates по сравнению с интерполяцией interp и IDL