code: def make_grid(dims, subsampling=(1, 1, 1), borders=(0, 0, 0)):
      slices = [
          slice(b, d - b, s)
          for d, s, b in zip(dims, subsampling, borders)
      ]
      xyz = np.mgrid[slices]
      xyz = np.rollaxis(xyz, 0, 4)
      xyz = np.reshape(xyz, [np.prod(xyz.shape[0:-1]), 3])
      return xyz
  
imports: import os
  import warnings
  import numpy
  from nibabel import io_orientation
  from nibabel.affines import apply_affine
  from ...core.image.image_spaces import as_xyz_image, make_xyz_image, xyz_affine
  from ...io.nibcompat import get_header
  from ..slicetiming import timefuncs
  from ._registration import _cspline_sample3d, _cspline_sample4d, _cspline_transform
  from .affine import Affine, Rigid
  from .optimizer import configure_optimizer, use_derivatives
  from .type_check import check_type, check_type_and_shape
label_desc: Use image node_id 3 for calling a global function with example usage: make_grid(dims, subsampling, borders) and returns: xyz