code: def test_bundle_shape_profile():
      cingulum_bundles = two_cingulum_bundles()
  
      cb1 = Streamlines(cingulum_bundles[0])
      cb1 = set_number_of_points(cb1, 20)
  
      cb2 = Streamlines(cingulum_bundles[1])
      cb2 = set_number_of_points(cb2, 20)
  
      deformed_bundle, affine_bundle, dists, mp, warp = bundlewarp(
          cb1, cb2
      )
  
      n = 10
      shape_profile, stdv = bundlewarp_shape_analysis(
          affine_bundle, deformed_bundle, no_disks=n
      )
  
      assert_equal(len(shape_profile), n)
      assert_equal(len(stdv), n)
  
      n = 100
      shape_profile, stdv = bundlewarp_shape_analysis(
          affine_bundle, deformed_bundle, no_disks=n
      )
  
      assert_equal(len(shape_profile), n)
      assert_equal(len(stdv), n)
  
imports: from numpy.testing import assert_equal
  from dipy.align.streamwarp import bundlewarp, bundlewarp_shape_analysis, bundlewarp_vector_filed
  from dipy.data import two_cingulum_bundles
  from dipy.tracking.streamline import set_number_of_points, Streamlines
label_desc: Use image node_id 3 for calling a global function with example usage: test_bundle_shape_profile() without return types