code: def pdf_interp_coords(sphere, rradius, origin):
      """Precompute coordinates for ODF calculation from the PDF
  
      Parameters
      ----------
      sphere : object,
              Sphere
      rradius : array, shape (N,)
              line interpolation points
      origin : array, shape (3,)
              center of the grid
  
      """
      interp_coords = rradius * sphere.vertices[np.newaxis].T
      origin = np.reshape(origin, [-1, 1, 1])
      interp_coords = origin + interp_coords
      return interp_coords
  
imports: import numpy
  from scipy.ndimage import map_coordinates
  from scipy.fftpack import fftn, fftshift, ifftshift
  from dipy.reconst.odf import OdfModel, OdfFit
  from dipy.reconst.cache import Cache
  from dipy.reconst.multi_voxel import multi_voxel_fit
label_desc: Use image node_id 4 for calling a global function with example usage: pdf_interp_coords(sphere, rradius, origin) and returns: interp_coords