code: def Wrotate_element(kt, indi, indj, indk, indl, B):
      r"""Compute the the specified index element of a kurtosis tensor rotated
      to the coordinate system basis B
  
      Parameters
      ----------
      kt : ndarray (x, y, z, 15) or (n, 15)
          Array containing the 15 independent elements of the kurtosis tensor
      indi : int
          Rotated kurtosis tensor element index i (0 for x, 1 for y, 2 for z)
      indj : int
          Rotated kurtosis tensor element index j (0 for x, 1 for y, 2 for z)
      indk : int
          Rotated kurtosis tensor element index k (0 for x, 1 for y, 2 for z)
      indl: int
          Rotated kurtosis tensor element index l (0 for x, 1 for y, 2 for z)
      B: array (x, y, z, 3, 3) or (n, 15)
          Vectors of the basis column-wise oriented
  
      Returns
      -------
      Wre : float
            rotated kurtosis tensor element of index ind_i, ind_j, ind_k, ind_l
  
      Notes
      -----
      It is assumed that initial kurtosis tensor elementes are defined on the
      Cartesian coordinate system.
  
      References
      ----------
      [1] Hui ES, Cheung MM, Qi L, Wu EX, 2008. Towards better MR
      characterization of neural tissues using directional diffusion kurtosis
      analysis. Neuroimage 42(1): 122-34
  
      """
      Wre = 0
  
      xyz = [0, 1, 2]
      for il in xyz:
          for jl in xyz:
              for kl in xyz:
                  for ll in xyz:
                      key = (il + 1) * (jl + 1) * (kl + 1) * (ll + 1)
                      multiplyB = (
                          B[..., il, indi]
                          * B[..., jl, indj]
                          * B[..., kl, indk]
                          * B[..., ll, indl]
                      )
                      Wre = Wre + multiplyB * kt[..., ind_ele[key]]
      return Wre
  
imports: import numpy
  import warnings
  import scipy.optimize
  import dipy.core.sphere
  from dipy.reconst.multi_voxel import multi_voxel_fit
  from dipy.reconst.dti import TensorFit, mean_diffusivity, from_lower_triangular, lower_triangular, decompose_tensor, MIN_POSITIVE_SIGNAL, nlls_fit_tensor, restore_fit_tensor
  from dipy.reconst.utils import dki_design_matrix
  from dipy.reconst.recspeed import local_maxima
  from dipy.reconst.base import ReconstModel
  from dipy.core.ndindex import ndindex
  from dipy.core.geometry import sphere2cart, cart2sphere, perpendicular_directions
  from dipy.core.optimize import PositiveDefiniteLeastSquares
  from dipy.data import get_sphere, get_fnames, load_sdp_constraints
  from dipy.reconst.vec_val_sum import vec_val_vect
  from dipy.core.gradients import check_multi_b
label_desc: Use image node_id 25 for calling a global function with example usage: Wrotate_element(kt, indi, indj, indk, indl, B) and returns: Wre