code: def box(shape, edges):
      data = np.zeros(shape)
      sl = []
      for i in range(len(shape)):
          sl.append(slice(edges[i][0], edges[i][1], 1))
      data[tuple(sl)] = 1
      return data.astype(np.int_)
  
imports: from itertools import combinations
  import numpy
  import numpy.linalg
  import pytest
  from numpy.testing import assert_almost_equal, assert_array_equal
  from ..None import intvol
label_desc: Use image node_id 3 for calling a global function with example usage: box(shape, edges) and returns: data