code: def test_integrate_2d():
      """test integration in 2d"""
      toto = np.random.rand(*shape)
      ddom = domain_from_binary_array(toto)
      ftoto = np.ravel(toto)
      f2 = np.vstack((ftoto, ftoto)).T
      ddom.set_feature("data", f2)
      ts = np.ones(2) * toto.sum()
      assert_almost_equal(ddom.integrate("data"), ts)
  
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 23 for calling a global function with example usage: test_integrate_2d() without return types