import numpy
in_data = numpy.array([[25.89125824, 25.88840675], [25.90930748, 
    25.90640068]], dtype=numpy.float32)
Scipy map_coordinates bilinear interpolation compared to interp and IDL interpolate