docs: def remapScipy(im, coords):
      """
      Remap an image using SciPy. See :func:`remap` for parameters.
      """
      height, width = im.shape[0], im.shape[1]
      
      # switch to y,x order
      coords = coords[:,:,::-1]
  
      # make it (h, w, 3, 3)
      coords_channels = np.empty((height, width, 3, 3))
      coords_channel = np.zeros((height, width, 3))
      coords_channel[:,:,:2] = coords
      coords_channels[:,:,0] = coords_channel
      coords_channels[:,:,1] = coords_channel
      coords_channels[:,:,1,2] = 1
      coords_channels[:,:,2] = coords_channel
      coords_channels[:,:,2,2] = 2
      coords = coords_channels
      
      # (3, h, w, 3)
      coords = np.rollaxis(coords, 3)
          
      return map_coordinates(im, coords, order=1)