new_result = np.zeros(im.shape)
coords = np.array([yy, xx])
for d in range(im.shape[2]):
    ndimage.map_coordinates(im[:, :, (d)], coords, order=1, prefilter=False,
        output=new_result[:, :, (d)])
toc = time.time()
print('interpolation time:', toc - tic)
Ускорение линейной интерполяции многих пикселей в NumPy