code: def test_circle():
      sh = (64, 61)
      cr = sh[0] // 2
      cc = sh[1] // 2
      x_0 = np.asarray(range(sh[0]))
      x_1 = np.asarray(range(sh[1]))
      X = np.empty((2,) + sh, dtype=np.float64)
      O = np.ones(sh)
      X[0, ...] = x_0[:, None] * O - cr
      X[1, ...] = x_1[None, :] * O - cc
      nrm = np.sqrt(np.sum(X**2, axis=0))
      for radius in [0, 7, 17, 32]:
          expected = nrm <= radius
          actual = vfu.create_circle(sh[0], sh[1], radius)
          assert_array_almost_equal(actual, expected)
  
imports: import numpy
  from numpy.testing import assert_array_equal, assert_array_almost_equal, assert_almost_equal, assert_equal, assert_raises
  from scipy.ndimage import map_coordinates
  from nibabel.affines import apply_affine, from_matvec
  from dipy.core import geometry
  from dipy.align import floating
  from dipy.align import imwarp
  from dipy.align import vector_fields
  from dipy.align.transforms import regtransforms
  from dipy.align.parzenhist import sample_domain_regular
  from dipy.testing.decorators import set_random_number_generator
label_desc: Use image node_id 5 for calling a global function with example usage: test_circle() without return types