code: def pdf_odf(Pr, rradius, interp_coords):
      r"""Calculates the real ODF from the diffusion propagator(PDF) Pr
  
      Parameters
      ----------
      Pr : array, shape (X, X, X)
          probability density function
      rradius : array, shape (N,)
          interpolation range on the radius
      interp_coords : array, shape (3, M, N)
          coordinates in the pdf for interpolating the odf
      """
      PrIs = map_coordinates(Pr, interp_coords, order=1)
      odf = (PrIs * rradius**2).sum(-1)
      return odf
  
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 5 for calling a global function with example usage: pdf_odf(Pr, rradius, interp_coords) and returns: odf