code: def generate_dataset(shape):
      """Generate a dataset with the described shape"""
      dim = len(shape)
      idx = np.reshape(np.indices(shape), (dim, -1)).T
      return idx
  
imports: import nibabel.gifti
  import numpy
  from nibabel import Nifti1Image
  from numpy.testing import assert_almost_equal, assert_array_equal
  from ..discrete_domain import domain_from_binary_array, domain_from_image, domain_from_mesh, grid_domain_from_binary_array, grid_domain_from_image, grid_domain_from_shape, smatrix_from_3d_array, smatrix_from_nd_array, smatrix_from_nd_idx
label_desc: Use image node_id 1 for calling a global function with example usage: generate_dataset(shape) and returns: idx