code: def create_qspace(gtab, origin):
      """create the 3D grid which holds the signal values (q-space)
  
      Parameters
      ----------
      gtab : GradientTable
      origin : (3,) ndarray
          center of qspace
  
      Returns
      -------
      qgrid : ndarray
          qspace coordinates
      """
      # create the q-table from bvecs and bvals
      qtable = create_qtable(gtab, origin)
  
      # center and index in qspace volume
      qgrid = qtable + origin
      return qgrid.astype("i8")
  
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 1 for calling a global function with example usage: create_qspace(gtab, origin) and returns: qgrid