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